Allelic imbalance analysis by high-density single-nucleotide polymorphic allele (SNP) array with whole genome amplified DNA

Besides their use in mRNA expression profiling, oligonucleotide microarrays have also been applied to single-nucleotide polymorphism (SNP) and loss of heterozygosity (LOH) or allelic imbalance studies. In this report, we evaluate the reliability of using whole genome amplified DNA for analysis with an oligonucleotide microarray containing 11 560 SNPs to detect allelic imbalance and chromosomal copy number abnormalities. Whole genome SNP analyses were performed with DNA extracted from osteosarcoma tissues and patient-matched blood. SNP calls were then generated by Affymetrix® GeneChip® DNA Analysis Software. In two osteosarcoma cases, using unamplified DNA, we identified 793 and 1070 SNP loci with allelic imbalance, respectively. In a parallel experiment with amplified DNA, 78% and 83% of these SNP loci with allelic imbalance was detected. The average false-positive rate is 13.8%. Furthermore, using the Affymetrix® GeneChip® Chromosome Copy Number Tool to analyze the SNP array data, we were able to detect identical chromosomal regions with gain or loss in both amplified and unamplified DNA at cytoband resolution.     

Quantitative evaluation by minisequencing and microarrays reveals accurate multiplexed SNP genotyping of whole genome amplified DNA

Whole genome amplification (WGA) procedures such as primer extension preamplification (PEP) or multiple displacement amplification (MDA) have the potential to provide an unlimited source of DNA for large-scale genetic studies. We have performed a quantitative evaluation of PEP and MDA for genotyping single nucleotide polymorphisms (SNPs) using multiplex, four-color fluorescent minisequencing in a microarray format. Forty-five SNPs were genotyped and the WGA methods were evaluated with respect to genotyping success, signal-to-noise ratios, power of genotype discrimination, yield and imbalanced amplification of alleles in the MDA product. Both PEP and MDA products provided genotyping results with a high concordance to genomic DNA. For PEP products the power of genotype discrimination was lower than for MDA due to a 2-fold lower signal-to-noise ratio. MDA products were indistinguishable from genomic DNA in all aspects studied. To obtain faithful representation of the SNP alleles at least 0.3 ng DNA should be used per MDA reaction. We conclude that the use of WGA, and MDA in particular, is a highly promising procedure for producing DNA in sufficient amounts even for genome wide SNP mapping studies.     

Assessing the utility of whole genome amplified DNA for next‐generation molecular ecology

DNA quantity can be a hindrance in ecological and evolutionary research programmes due to a range of factors including endangered status of target organisms, available tissue type, and the impact of field conditions on preservation methods. A potential solution to low‐quantity DNA lies in whole genome amplification (WGA) techniques that can substantially increase DNA yield. To date, few studies have rigorously examined sequence bias that might result from WGA and next‐generation sequencing of nonmodel taxa. To address this knowledge deficit, we use multiple displacement amplification (MDA) and double‐digest RAD sequencing on the grey mouse lemur (Microcebus murinus) to quantify bias in genome coverage and SNP calls when compared to raw genomic DNA (gDNA). We focus our efforts in providing baseline estimates of potential bias by following manufacturer's recommendations for starting DNA quantities (>100 ng). Our results are strongly suggestive that MDA enrichment does not introduce systematic bias to genome characterization. SNP calling between samples when genotyping both de‐novo and with a reference genome are highly congruent (>98%) when specifying a minimum threshold of 20X stack depth to call genotypes. Relative genome coverage is also similar between MDA and gDNA, and allelic dropout is not observed. SNP concordance varies based on coverage threshold, with 95% concordance reached at ~12X coverage genotyping de‐novo and ~7X coverage genotyping with the reference genome. These results suggest that MDA may be a suitable solution for next‐generation molecular ecological studies when DNA quantity would otherwise be a limiting factor.     

Effects of DNA mass on multiple displacement whole genome amplification and genotyping performance

Background

Whole genome amplification (WGA) promises to eliminate practical molecular genetic analysis limitations associated with genomic DNA (gDNA) quantity. We evaluated the performance of multiple displacement amplification (MDA) WGA using gDNA extracted from lymphoblastoid cell lines (N = 27) with a range of starting gDNA input of 1–200 ng into the WGA reaction. Yield and composition analysis of whole genome amplified DNA (wgaDNA) was performed using three DNA quantification methods (OD, PicoGreen^® and RT-PCR). Two panels of N = 15 STR (using the AmpFlSTR^® Identifiler^® panel) and N = 49 SNP (TaqMan^®) genotyping assays were performed on each gDNA and wgaDNA sample in duplicate. gDNA and wgaDNA masses of 1, 4 and 20 ng were used in the SNP assays to evaluate the effects of DNA mass on SNP genotyping assay performance. A total of N = 6,880 STR and N = 56,448 SNP genotype attempts provided adequate power to detect differences in STR and SNP genotyping performance between gDNA and wgaDNA, and among wgaDNA produced from a range of gDNA templates inputs.

Results

The proportion of double-stranded wgaDNA and human-specific PCR amplifiable wgaDNA increased with increased gDNA input into the WGA reaction. Increased amounts of gDNA input into the WGA reaction improved wgaDNA genotyping performance. Genotype completion or genotype concordance rates of wgaDNA produced from all gDNA input levels were observed to be reduced compared to gDNA, although the reduction was not always statistically significant. Reduced wgaDNA genotyping performance was primarily due to the increased variance of allelic amplification, resulting in loss of heterozygosity or increased undetermined genotypes. MDA WGA produces wgaDNA from no template control samples; such samples exhibited substantial false-positive genotyping rates.

Conclusion

The amount of gDNA input into the MDA WGA reaction is a critical determinant of genotyping performance of wgaDNA. At least 10 ng of lymphoblastoid gDNA input into MDA WGA is required to obtain wgaDNA TaqMan^® SNP assay genotyping performance equivalent to that of gDNA. Over 100 ng of lymphoblastoid gDNA input into MDA WGA is required to obtain optimal STR genotyping performance using the AmpFlSTR^® Identifiler^® panel from wgaDNA equivalent to that of gDNA.

     

Cubic regression-based degree of correction predicts the performance of whole bisulfitome amplified DNA methylation analysis

Epigenetic mechanisms, including DNA methylation, are important determinants in development and disease. There is a need for technologies capable of detecting small variations in methylation levels in an accurate and reproducible manner, even if only limited amounts of DNA are available (which is the case in many studies in humans). Quantitative methylation analysis of minute DNA amounts after whole bisulfitome amplification (qMAMBA) has been proposed as an alternative, but this technique has not been adequately standardized and no comparative study against conventional methods has been performed, that includes a wide range of methylation percentages and different target assays. We designed an experiment to compare the performance of qMAMBA and bisulfite-treated genomic (non-amplified) DNA pyrosequencing. Reactions were performed in duplicate for each technique in eight different target genes, using nine artificially constructed DNA samples with methylation levels ranging between 0% and 100% with intervals of 12.5%. Cubic polynomial curves were plotted from the experimental results and the real methylation values and the resulting equation was used to estimate new corrected data points. The use of the cubic regression-based correction benefits the accuracy and the power of discrimination in methylation studies. Additionally, dispersion of the new estimated data around a y = x line (R²) served to fix a cutoff that can discriminate, with a single 9-point curve experiment, whether whole bisulfitome amplification and subsequent qMAMBA can produce accurate methylation results. Finally, even with an optimized reagent kit, DNA samples subjected to whole bisulfitome amplification enhance the preferential amplification of unmethylated alleles, and subtle changes in methylation levels cannot be detected confidently.     

Cubic regression-based degree of correction predicts the performance of whole bisulfitome amplified DNA methylation analysis

Epigenetic mechanisms, including DNA methylation, are important determinants in development and disease. There is a need for technologies capable of detecting small variations in methylation levels in an accurate and reproducible manner, even if only limited amounts of DNA are available (which is the case in many studies in humans). Quantitative methylation analysis of minute DNA amounts after whole bisulfitome amplification (qMAMBA) has been proposed as an alternative, but this technique has not been adequately standardized and no comparative study against conventional methods has been performed, that includes a wide range of methylation percentages and different target assays. We designed an experiment to compare the performance of qMAMBA and bisulfite-treated genomic (non-amplified) DNA pyrosequencing. Reactions were performed in duplicate for each technique in eight different target genes, using nine artificially constructed DNA samples with methylation levels ranging between 0% and 100% with intervals of 12.5%. Cubic polynomial curves were plotted from the experimental results and the real methylation values and the resulting equation was used to estimate new corrected data points. The use of the cubic regression-based correction benefits the accuracy and the power of discrimination in methylation studies. Additionally, dispersion of the new estimated data around a y = x line (R²) served to fix a cutoff that can discriminate, with a single 9-point curve experiment, whether whole bisulfitome amplification and subsequent qMAMBA can produce accurate methylation results. Finally, even with an optimized reagent kit, DNA samples subjected to whole bisulfitome amplification enhance the preferential amplification of unmethylated alleles, and subtle changes in methylation levels cannot be detected confidently.     

Optimized design and assessment of whole genome tiling arrays

MOTIVATION: Recent advances in microarray technologies have made it feasible to interrogate whole genomes with tiling arrays and this technique is rapidly becoming one of the most important high-throughput functional genomics assays. For large mammalian genomes, analyzing oligonucleotide tiling array data is complicated by the presence of non-unique sequences on the array, which increases the overall noise in the data and may lead to false positive results due to cross-hybridization. The ability to create custom microarrays using maskless array synthesis has led us to consider ways to optimize array design characteristics for improving data quality and analysis. We have identified a number of design parameters to be optimized including uniqueness of the probe sequences within the whole genome, melting temperature and self-hybridization potential. RESULTS: We introduce the uniqueness score, U, a novel quality measure for oligonucleotide probes and present a method to quickly compute it. We show that U is equivalent to the number of shortest unique substrings in the probe and describe an efficient greedy algorithm to design mammalian whole genome tiling arrays using probes that maximize U. Using the mouse genome, we demonstrate how several optimizations influence the tiling array design characteristics. With a sensible set of parameters, our designs cover 78% of the mouse genome including many regions previously considered 'untilable' due to the presence of repetitive sequence. Finally, we compare our whole genome tiling array designs with commercially available designs. AVAILABILITY: Source code is available under an open source license from http://www.ebi.ac.uk/~graef/arraydesign/.     

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.     

Genotyping of human class I alcohol dehydrogenase. Analysis of enzymatically amplified DNA with allele-specific oligonucleotides

Large inter-individual differences are noted in the susceptibility to alcohol-related problems. Part of this variation may be due to the different isoenzyme patterns of the alcohol-metabolizing enzymes and, consequently, different pharmacokinetics of alcohol degradation. We have used the polymerase chain reaction and oligonucleotide hybridization to amplify and analyze class I alcohol dehydrogenase isoenzyme-specific genomic DNA. The method unambiguously distinguishes between different allelic variants and thus provides a new means of elucidating the alcohol dehydrogenase isoenzyme pattern of humans.     

Simplified method for the collection, storage, and comet assay analysis of DNA damage in whole blood

Single-cell gel electrophoresis (comet assay) is one of the most common methods used to measure oxidatively damaged DNA in peripheral blood mononuclear cells (PBMC), as a biomarker of oxidative stress in vivo. However, storage, extraction, and assay workup of blood samples are associated with a risk of artifactual formation of damage. Previous reports using this approach to study DNA damage in PBMC have, for the most part, required the isolation of PBMC before immediate analysis or freezing in cryopreservative. This is very time-consuming and a significant drain on human resources. Here, we report the successful storage of whole blood in ~ 250 μl volumes, at − 80 °C, without cryopreservative, for up to 1 month without artifactual formation of DNA damage. Furthermore, this blood is amenable for direct use in both the alkaline and the enzyme-modified comet assay, without the need for prior isolation of PBMC. In contrast, storage of larger volumes (e.g., 5 ml) of whole blood leads to an increase in damage with longer term storage even at − 80 °C, unless a cryopreservative is present. Our “small volume” approach may be suitable for archived blood samples, facilitating analysis of biobanks when prior isolation of PBMC has not been performed.     

An information-based sequence distance and its application to whole mitochondrial genome phylogeny

MOTIVATION: Traditional sequence distances require an alignment and therefore are not directly applicable to the problem of whole genome phylogeny where events such as rearrangements make full length alignments impossible. We present a sequence distance that works on unaligned sequences using the information theoretical concept of Kolmogorov complexity and a program to estimate this distance. RESULTS: We establish the mathematical foundations of our distance and illustrate its use by constructing a phylogeny of the Eutherian orders using complete unaligned mitochondrial genomes. This phylogeny is consistent with the commonly accepted one for the Eutherians. A second, larger mammalian dataset is also analyzed, yielding a phylogeny generally consistent with the commonly accepted one for the mammals. AVAILABILITY: The program to estimate our sequence distance, is available at http://www.cs.cityu.edu.hk/~cssamk/gencomp/GenCompress1.htm. The distance matrices used to generate our phylogenies are available at http://www.math.uwaterloo.ca/~mli/distance.html.     

Characterization of miRNAs responsive to exogenous ethylene in grapevine berries at whole genome level

MicroRNAs (miRNAs) are critical regulators of various biological and metabolic processes of plants. Numerous miRNAs and their functions have been identified and analyzed in many plants. However, till now, the involvement of miRNAs in the response of grapevine berries to ethylene has not been reported yet. Here, Solexa technology was employed to deeply sequence small RNA libraries constructed from grapevine berries treated with and without ethylene. A total of 124 known and 78 novel miRNAs were identified. Among these miRNAs, 162 miRNAs were clearly responsive to ethylene, with 55 downregulated, 59 upregulated, and 14 unchanged miRNAs detected only in the control. The other 35 miRNAs responsive to ethylene were induced by ethylene and detected only in the ethylene-treated grapevine materials. Expression analysis of 27 conserved and 26 novel miRNAs revealed that 13 conserved and 18 novel ones were regulated by ethylene during the whole development of grapevine berries. High-throughput sequencing and qRT-PCR assays revealed consistent results on the expression results of ethylene-responsive miRNAs. Moreover, 90 target genes for 34 novel miRNAs were predicted, most of which were involved in responses to various stresses, especially like exogenous ethylene treatment. The identified miRNAs may be mainly involved in grapevine berry development and response to various environmental conditions.     

Identification and genetic mapping of random amplified polymorphic DNA (RAPD) markers to the sheep genome

The random amplified polymorphic DNA (RAPD) assay utilizes the polymerase chain reaction (PCR) and short primers of arbitrary nucleotide sequence to amplify DNA. In this study, the RAPD assay was used to identify and map polymorphic markers in the AgResearch International Mapping Flock (IMF) sheep pedigrees. Sires and dams of eight of the full-sib IMF pedigrees were screened with 131 different 10-mer oligonucleotide primers. An average of 85 RAPD polymorphisms was identified between each parental pair, and 53 markers were contributed to the AgResearch IMF collaboration. Forty-five of the RAPD markers were mapped in the AgResearch IMF genetic linkage map, and at least one marker was located on 17 of the 26 autosomes and both sex chromosomes. Three lines of evidence were used to check for the homology of scored polymorphisms in different pedigrees, pedigree evaluation, segregation analysis, and Southern blot analysis. These results demonstrate that the RAPD assay is a powerful approach for identifying polymorphisms that can be used as markers for constructing a sheep genetic linkage map. Received: 5 October 1995 / Accepted: 16 April 1996     

Extracting DNA from museum bird eggs, and whole genome amplification of archive DNA

We present a comprehensive protocol for extracting DNA from egg membranes and other internal debris recovered from the interior of blown museum bird eggs. A variety of commercially available DNA extraction methods were found to be applicable. DNA sequencing of polymerase chain reaction (PCR) products for a 176‐bp fragment of mitochondrial DNA was successful for most egg samples (> 78%) even though the amount of DNA extracted (mean = 14.71 ± 4.55 ng/µL) was significantly less than that obtained for bird skin samples (mean = 67.88 ± 4.77 ng/µL). For PCR and sequencing of snipe (Gallinago) DNA, we provide eight new primers for the ‘DNA barcode’ region of COI mtDNA. In various combinations, the primers target a range of PCR products sized from 72 bp to the full ‘barcode’ of 751 bp. Not all possible combinations were tested with archive snipe DNA, but we found a significantly better success rate of PCR amplification for a shorter 176‐bp target compared with a larger 288‐bp fragment (67% vs. 39%). Finally, we explored the feasibility of whole genome amplification (WGA) for extending the use of archive DNA in PCR and sequencing applications. Of two WGA approaches, a PCR‐based method was found to be able to amplify whole genomic DNA from archive skins and eggs from museum bird collections. After WGA, significantly more archive egg samples produced visible PCR products on agarose (56.9% before WGA vs. 79.0% after WGA). However, overall sequencing success did not improve significantly (78.8% compared with 83.0%).     

Measuring oxidative damage to DNA and its repair with the comet assay

Background

Single cell gel electrophoresis, or the comet assay, was devised as a sensitive method for detecting DNA strand breaks, at the level of individual cells. A simple modification, incorporating a digestion of DNA with a lesion-specific endonuclease, makes it possible to measure oxidised bases.

Scope of review

With the inclusion of formamidopyrimidine DNA glycosylase to recognise oxidised purines, or Nth (endonuclease III) to detect oxidised pyrimidines, the comet assay has been used extensively in human biomonitoring to monitor oxidative stress, usually in peripheral blood mononuclear cells.

Major conclusions

There is evidence to suggest that the enzymic approach is more accurate than chromatographic methods, when applied to low background levels of base oxidation. However, there are potential problems of over-estimation (because the enzymes are not completely specific) or under-estimation (failure to detect lesions that are close together). Attempts have been made to improve the inter-laboratory reproducibility of the comet assay.

General significance

In addition to measuring DNA damage, the assay can be used to monitor the cellular or in vitro repair of strand breaks or oxidised bases. It also has applications in assessing the antioxidant status of cells. In its various forms, the comet assay is now an invaluable tool in human biomonitoring and genotoxicity testing. This article is part of a Special Issue entitled Current methods to study reactive oxygen species - pros and cons and biophysics of membrane proteins. Guest Editor: Christine Winterbourn.     

DNA metabarcoding multiplexing and validation of data accuracy for diet assessment: application to omnivorous diet

Ecological understanding of the role of consumer–resource interactions in natural food webs is limited by the difficulty of accurately and efficiently determining the complex variety of food types animals have eaten in the field. We developed a method based on DNA metabarcoding multiplexing and next‐generation sequencing to uncover different taxonomic groups of organisms from complex diet samples. We validated this approach on 91 faeces of a large omnivorous mammal, the brown bear, using DNA metabarcoding markers targeting the plant, vertebrate and invertebrate components of the diet. We included internal controls in the experiments and performed PCR replication for accuracy validation in postsequencing data analysis. Using our multiplexing strategy, we significantly simplified the experimental procedure and accurately and concurrently identified different prey DNA corresponding to the targeted taxonomic groups, with ≥60% of taxa of all diet components identified to genus/species level. The systematic application of internal controls and replication was a useful and simple way to evaluate the performance of our experimental procedure, standardize the selection of sequence filtering parameters for each marker data and validate the accuracy of the results. Our general approach can be adapted to the analysis of dietary samples of various predator species in different ecosystems, for a number of conservation and ecological applications entailing large‐scale population level diet assessment through cost‐effective screening of multiple DNA metabarcodes, and the detection of fine dietary variation among samples or individuals and of rare food items.     

Synthesis of universal unmethylated control DNA by nested whole genome amplification with phi29 DNA polymerase

Optimization of highly sensitive methods to detect methylation of CpG islands in gene promoter regions requires adequate methylated and unmethylated control DNA. Whereas universal methylated control DNA is available, universal unmethylated control (UUC) DNA has not been made because demethylase is not available to remove methyl groups from all methylated cytosines. On the basis that DNA synthesized by DNA polymerase does not contain methylated cytosines, we developed a method to create UUC DNA by nested whole genome amplification (WGA) with phi29 DNA polymerase. Contamination of the template genomic DNA in UUC was only 3.1 x 10(-7), below the detection limit of sensitive methods used for methylation studies such as methylation-specific PCR. Assessment of microsatellite markers demonstrated that even nested phi29 WGA achieves highly accurate and homogeneous amplification with very low amounts of genomic DNA as an initial template. The UUC DNA created by nested phi29 WGA is practically very useful for methylation analysis.     

Operation of partial nitrification to nitrite of landfill leachate and its performance with respect to different oxygen conditions

The coupled system of partial nitrification and anaerobic ammonium oxidation (Anammox) is efficient in nitrogen removal from wastewater. In this study, the effect of different oxygen concentrations on partial nitrification performance with a sequencing batch reactor (SBR) was investigated. Results indicate that, partial nitrification of landfill leachate could be successfully achieved under the 1.0–2.0 mg L⁻¹ dissolved oxygen (DO) condition after 118 d long-term operation, and that the effluent is suitable for an Anammox reactor. Further decreasing or increasing the DO concentration, however, would lead to a decay of nitrification performance. Additionally, the MLSS concentration in the reactor increased with increasing DO concentration. Respirometric assays suggest that low DO conditions (<2 mg L⁻¹) favor the ammonia-oxidizing bacteria (AOB) and significantly inhibit nitrite oxidizing bacteria (NOB) and aerobic heterotrophic bacteria (AHB); whereas high DO conditions (>3 mg L⁻¹) allow AHB to dominate and significantly inhibit AOB. Therefore, the optimal condition for partial nitrification of landfill leachate is 1.0–2.0 mg L⁻¹ DO concentration.