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.     

Comparative genome analysis of Campylobacter jejuni using whole genome DNA microarrays

Whole genome DNA microarrays were constructed and used to investigate genomic diversity in 18 Campylobacter jejuni strains from diverse sources. New algorithms were developed that dynamically determine the boundary between the conserved and variable genes. Seven hypervariable plasticity regions (PR) were identified in the genome (PR1 to PR7) containing 136 genes (50%) of the variable gene pool. When comparisons were made with the sequenced strain NCTC11168, the number of absent or divergent genes ranged from 2.6% (40 genes) to 10.2% (163) and in total 16.3% (269) of the genes were variable. PR1 contains genes important in the utilisation of alternative electron acceptors for respiration and may confer a selective advantage to strains in restricted oxygen environments. PR2, 3 and 7 contain many outer membrane and periplasmic proteins and hypothetical proteins of unknown function that might be linked to phenotypic variation and adaptation to different ecological niches. PR4, 5 and 6 contain genes involved in the production and modification of antigenic surface structures.     

Ordered catenation of sequence-tagged sites and multiplexed SNP genotyping by sequencing

We describe a method for the efficient genotyping of SNPs, involving sequencing of ordered and catenated sequence-tagged sites (OCS). In OCS, short genomic segments, each containing an SNP, are amplified by PCR using primers that carry specially designed extra nucleotides at their 5′-ends. Amplification products are then combined and converted to a concatamer in a defined order by a second round of thermal cycling. The concatenation takes place because the 5′-ends of each amplicon are designed to be complementary to the ends of the presumptive neighboring amplicons. The primer sequences for OCS are chosen using newly developed dedicated software, OCS Optimizer. Using sets of SNPs, we show that at least 10 STSs can be concatenated in a predefined order and all SNPs in the STSs are accurately genotyped by one two-way sequencing reaction.     

High-resolution genotyping via whole genome hybridizations to microarrays containing long oligonucleotide probes

To date, microarray-based genotyping of large, complex plant genomes has been complicated by the need to perform genome complexity reduction to obtain sufficiently strong hybridization signals. Genome complexity reduction techniques are, however, tedious and can introduce unwanted variables into genotyping assays. Here, we report a microarray-based genotyping technology for complex genomes (such as the 2.3 GB maize genome) that does not require genome complexity reduction prior to hybridization. Approximately 200,000 long oligonucleotide probes were identified as being polymorphic between the inbred parents of a mapping population and used to genotype two recombinant inbred lines. While multiple hybridization replicates provided ～97% accuracy, even a single replicate provided ～95% accuracy. Genotyping accuracy was further increased to >99% by utilizing information from adjacent probes. This microarray-based method provides a simple, high-density genotyping approach for large, complex genomes.     

Quantitative PCR assessment of microsatellite and SNP genotyping with variable quality DNA extracts

In this study we developed eight quantitative PCR (qPCR) assays to evaluate the starting copy number of nuclear and mitochondrial DNA fragments ranging from 75 to 350 base-pairs in DNA extracts from Chinook salmon tissues with varying quality. Samples were genotyped with 13 microsatellite and 29 SNP assays and average genotyping success for good, intermediate, and poor quality samples was 96%, 24%, and 24% for microsatellite loci, and 98%, 97%, and 79% for SNPs, respectively. As measured by qPCR, good quality samples had a consistently high number of starting copies across all fragment sizes with little change between the smallest and largest size. In contrast, the intermediate and poor quality samples displayed decreases in starting copy number as fragment size increased, and was most pronounced with poor samples. Logistic regression of genotyping success by starting copy number indicated that in order to achieve at least 90% genotyping success, approximately 1,000 starting copies of nuclear DNA are necessary for microsatellite loci, and as few as 14 starting copies for SNP assays (but we recommend at least 50 copies to reduce genotyping error). While these guidelines apply specifically to Chinook salmon and the genetic markers included in this study, the principles are transferable to other species and markers due to the underlying process associated with template quantity and PCR amplification.     

Quantitative analysis of human DNA sequences by PCR and solid-phase minisequencing

Reliable quantification by PCR requires careful experimental design and conditions, often involving sampling of the PCR reactions at different time points or amplifying multiple dilutions of a standard DNA. We describe here an accurate, quantitative and easily automatizable solid-phase method based on competetive PCR. The PCR products are analyzed by solid-phase minisequencing after capture of biotinylated PCR products in streptavidin-coated microtiter wells and single-nucleotide extension of a specific detection primer by a radioactively labelled nucleotide. The results are expressed as numeric cpm-values, and the incorporated label expresses the relative amount of sequence variants in the original template mixture. We have applied the method to determination of allele frequencies in pooled DNA samples, of mitochondrial heteroplasmy, of gene copy numbers, and to forensic DNA analysis.     

MARA: a novel approach for highly multiplexed locus-specific SNP genotyping using high-density DNA oligonucleotide arrays

We have developed a locus-specific DNA target preparation method for highly multiplexed single nucleotide polymorphism (SNP) genotyping called MARA (Multiplexed Anchored Runoff Amplification). The approach uses a single primer per SNP in conjunction with restriction enzyme digested, adapter-ligated human genomic DNA. Each primer is composed of common sequence at the 5′ end followed by locus-specific sequence at the 3′ end. Following a primary reaction in which locus-specific products are generated, a secondary universal amplification is carried out using a generic primer pair corresponding to the oligonucleotide and genomic DNA adapter sequences. Allele discrimination is achieved by hybridization to high-density DNA oligonucleotide arrays. Initial multiplex reactions containing either 250 primers or 750 primers across nine DNA samples demonstrated an average sample call rate of ~95% for 250- and 750-plex MARA. We have also evaluated >1000- and 4000-primer plex MARA to genotype SNPs from human chromosome 21. We have identified a subset of SNPs corresponding to a primer conversion rate of ~75%, which show an average call rate over 95% and concordance >99% across seven DNA samples. Thus, MARA may potentially improve the throughput of SNP genotyping when coupled with allele discrimination on high-density arrays by allowing levels of multiplexing during target generation that far exceed the capacity of traditional multiplex PCR.     

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.

     

Resolving individuals contributing trace amounts of DNA to highly complex mixtures using high-density SNP genotyping microarrays

We use high-density single nucleotide polymorphism (SNP) genotyping microarrays to demonstrate the ability to accurately and robustly determine whether individuals are in a complex genomic DNA mixture. We first develop a theoretical framework for detecting an individual's presence within a mixture, then show, through simulations, the limits associated with our method, and finally demonstrate experimentally the identification of the presence of genomic DNA of specific individuals within a series of highly complex genomic mixtures, including mixtures where an individual contributes less than 0.1% of the total genomic DNA. These findings shift the perceived utility of SNPs for identifying individual trace contributors within a forensics mixture, and suggest future research efforts into assessing the viability of previously sub-optimal DNA sources due to sample contamination. These findings also suggest that composite statistics across cohorts, such as allele frequency or genotype counts, do not mask identity within genome-wide association studies. The implications of these findings are discussed.     

SNP typing by apyrase-mediated allele-specific primer extension on DNA microarrays

This study reports the development of a microarray-based allele-specific extension method for typing of single nucleotide polymorphisms (SNPs). The use of allele-specific primers has been employed previously to identify single base variations but it is acknowledged that certain mismatches are not refractory to extension. Here we have overcome this limitation by introducing apyrase, a nucleotide-degrading enzyme, to the extension reaction. We have shown previously that DNA polymerases exhibit slower reaction kinetics when extending a mismatched primer compared with a matched primer. This kinetic difference is exploited in the apyrase-mediated allele-specific extension (AMASE) assay, allowing incorporation of nucleotides when the reaction kinetics are fast but degrading the nucleotides before extension when the reaction kinetics are slow. Here we show that five homozygous variants (14% of the total number of variants) that were incorrectly scored in the absence of apyrase were correctly typed when apyrase was included in the extension reaction. AMASE was performed in situ on the oligonucleotide microarrays using fluorescent nucleotides to type 10 SNPs and two indels in 17 individuals generating approximately 200 genotypes. Cluster analysis of these data shows three distinct clusters with clear-cut boundaries. We conclude that SNP typing on oligonucleotide microarrays by AMASE is an efficient, rapid and accurate technique for large-scale genotyping.     

Highly accurate SNP genotyping from historical and low‐quality samples

Historical and other poor‐quality samples are often necessary for population genetics, conservation, and forensics studies. Although there is a long history of using mtDNA from such samples, obtaining and genotyping nuclear loci have been considered difficult and error‐prone at best, and impossible at worst. The primary issues are the amount of nuclear DNA available for genotyping, and the degradation of the DNA into small fragments. Single nucleotide polymorphisms offer potential advantages for assaying nuclear variation in historical and poor‐quality samples, because the amplified fragments can be very small, varying little or not at all in size between alleles, and can be amplified efficiently by polymerase chain reaction (PCR). We present a method for highly multiplexed PCR of SNP loci, followed by dual‐fluorescence genotyping that is very effective for genotyping poor‐quality samples, and can potentially use very little template DNA, regardless of the number of loci to be genotyped. We genotyped 19 SNP loci from DNA extracted from modern and historical bowhead whale tissue, bone and baleen samples. The PCR failure rate was < 1.5%, and the genotyping error rate was 0.1% when DNA samples contained > 10 copies/µL of a 51‐bp nuclear sequence. Among samples with ≤ 10 copies/µL DNA, samples could still be genotyped confidently with appropriate levels of replication from independent multiplex PCRs.     

Quantitative detection of microbial genes by using DNA microarrays

To quantify target genes in biological samples using DNA microarrays, we employed reference DNA to normalize variations in spot size and hybridization. This method was tested using nitrate reductase (nirS), naphthalene dioxygenase (nahA), and Escherichia coli O157 O-antigen biosynthesis genes as model genes and lambda DNA as the reference DNA. We observed a good linearity between the log signal ratio and log DNA concentration ratio at DNA concentrations above the method's detection limit, which was approximately 10 pg. This approach for designing quantitative microarrays and the inferred equation from this study provide a simple and convenient way to estimate the target gene concentration from the hybridization signal ratio.     

Concurrent genotyping of <Emphasis Type="Italic">Helicobacter pylori</Emphasis>virulence genes and human cytokine SNP sites using whole genome amplified DNA derived from minute amounts of gastric biopsy specimen DNA

Background  

Bacterial and cellular genotyping is becoming increasingly important in the diagnosis of infectious diseases. However, difficulties in obtaining sufficient amount of bacterial and cellular DNA extracted from the same human biopsy specimens is often a limiting factor. In this study, total DNA (host and bacterial DNA) was isolated from minute amounts of gastric biopsy specimens and amplified by means of whole genome amplification using the multiple displacement amplification (MDA) technique. Subsequently, MDA-DNA was used for concurrent Helicobacter pylori and human host cellular DNA genotyping analysis using PCR-based methods.     

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.     

Quantitative trait loci for IQ and other complex traits: single-nucleotide polymorphism genotyping using pooled DNA and microarrays

Similar to other complex traits, it is likely that many DNA polymorphisms of small effect size [quantitative trait loci (QTLs)] are responsible for the high heritability of intelligence, in addition to many rare monogenic disorders known to contribute to lowered intelligence. We review the current status of approaches to identify QTLs associations for intelligence employing genome-wide strategies using pooled DNA from many individuals and evaluate the innovative approach of microarray analysis to genotype DNA pools for large numbers of single nucleotide polymorphisms.     

Multiplex SNP genotyping in pooled DNA samples by a four-colour microarray system

We selected 125 candidate single nucleotide polymorphisms (SNPs) in genes belonging to the human type 1 interferon (IFN) gene family and the genes coding for proteins in the main type 1 IFN signalling pathway by screening databases and by in silico comparison of DNA sequences. Using quantitative analysis of pooled DNA samples by solid-phase mini-sequencing, we found that only 20% of the candidate SNPs were polymorphic in the Finnish and Swedish populations. To allow more effective validation of candidate SNPs, we developed a four-colour microarray-based mini-sequencing assay for multiplex, quantitative allele frequency determination in pooled DNA samples. We used cyclic mini-sequencing reactions with primers carrying 5′-tag sequences, followed by capture of the products on microarrays by hybridisation to complementary tag oligonucleotides. Standard curves prepared from mixtures of known amounts of SNP alleles demonstrate the applicability of the system to quantitative analysis, and showed that for about half of the tested SNPs the limit of detection for the minority allele was below 5%. The microarray-based genotyping system established here is universally applicable for genotyping and quantification of any SNP, and the validated system for SNPs in type 1 IFN-related genes should find many applications in genetic studies of this important immunoregulatory pathway.     

SNP genotyping reveals genetic diversity between cultivated landraces and contemporary varieties of tomato

Background

The tomato (Solanum lycopersium L.) is the most widely grown vegetable in the world. It was domesticated in Latin America and Italy and Spain are considered secondary centers of diversification. This food crop has experienced severe genetic bottlenecks and modern breeding activities have been characterized by trait introgression from wild species and divergence in different market classes.

Results

With the aim to examine patterns of polymorphism, characterize population structure and identify putative loci under positive selection, we genotyped 214 tomato accessions (which include cultivated landraces, commercial varieties and wild relatives) using a custom-made Illumina SNP-panel. Most of the 175 successfully scored SNP loci were found to be polymorphic. Population structure analysis and estimates of genetic differentiation indicated that landraces constitute distinct sub-populations. Furthermore, contemporary varieties could be separated in groups (processing, fresh and cherry) that are consistent with the recent breeding aimed at market-class specialization. In addition, at the 95% confidence level, we identified 30, 34 and 37 loci under positive selection between landraces and each of the groups of commercial variety (cherry, processing and fresh market, respectively). Their number and genomic locations imply the presence of some extended regions with high genetic variation between landraces and contemporary varieties.

Conclusions

Our work provides knowledge concerning the level and distribution of genetic variation within cultivated tomato landraces and increases our understanding of the genetic subdivision of contemporary varieties. The data indicate that adaptation and selection have led to a genomic signature in cultivated landraces and that the subpopulation structure of contemporary varieties is shaped by directed breeding and largely of recent origin. The genomic characterization presented here is an essential step towards a future exploitation of the available tomato genetic resources in research and breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-835) contains supplementary material, which is available to authorized users.     

Quantification of trace-level DNA by real-time whole genome amplification

Quantification of trace amounts of DNA is a challenge in analytical applications where the concentration of a target DNA is very low or only limited amounts of samples are available for analysis. PCR-based methods including real-time PCR are highly sensitive and widely used for quantification of low-level DNA samples. However, ordinary PCR methods require at least one copy of a specific gene sequence for amplification and may not work for a sub-genomic amount of DNA. We suggest a real-time whole genome amplification method adopting the degenerate oligonucleotide primed PCR (DOP-PCR) for quantification of sub-genomic amounts of DNA. This approach enabled quantification of sub-picogram amounts of DNA independently of their sequences. When the method was applied to the human placental DNA of which amount was accurately determined by inductively coupled plasma-optical emission spectroscopy (ICP-OES), an accurate and stable quantification capability for DNA samples ranging from 80 fg to 8 ng was obtained. In blind tests of laboratory-prepared DNA samples, measurement accuracies of 7.4%, -2.1%, and -13.9% with analytical precisions around 15% were achieved for 400-pg, 4-pg, and 400-fg DNA samples, respectively. A similar quantification capability was also observed for other DNA species from calf, E. coli, and lambda phage. Therefore, when provided with an appropriate standard DNA, the suggested real-time DOP-PCR method can be used as a universal method for quantification of trace amounts of DNA.