An isothermal primer extension method for whole genome amplification of fresh and degraded DNA: applications in comparative genomic hybridization, genotyping and mutation screening

We describe a protocol that uses a bioinformatically optimized primer in an isothermal whole genome amplification (WGA) reaction. Overnight incubation at 37 degrees C efficiently generates several hundred- to several thousand-fold increases in input DNA. The amplified product retains reasonably faithful quantitative representation of unamplified whole genomic DNA (gDNA). We provide protocols for applying this isothermal primer extension WGA protocol in three different techniques of genomic analysis: comparative genomic hybridization (CGH), genotyping at simple tandem repeat (STR) loci and screening for single base mutations in a common monogenic disorder, beta-thalassemia. gDNA extracted from formalin-fixed paraffin-embedded (FFPE) tissues can also be amplified with this protocol.     

An isothermal method for whole genome amplification of fresh and degraded DNA for comparative genomic hybridization, genotyping and mutation detection.

Molecular genotyping has important biomedical and forensic applications. However, limiting amounts of human biological material often yield genomic DNA (gDNA) in insufficient quantity and of poor quality for a reliable analysis. This motivated the development of an efficient whole genome amplification method with quantitatively unbiased representation usable on fresh and degraded gDNA. Amplification of fresh frozen, formalin-fixed paraffin-embedded (FFPE) and DNase-degraded DNA using degenerate oligonucleotide-primed PCR or primer extension amplification using a short primer sequence bioinformatically optimized for coverage of the human genome was compared with amplification using current primers by chromosome-based and BAC-array comparative genomic hybridization (CGH), genotyping at short tandem repeats (STRs) and single base mutation detection. Compared with current primers, genome amplification using the bioinformatically optimized primer was significantly less biased on CGH in self-self hybridizations, and replicated tumour genome copy number aberrations, even from FFPE tissue. STR genotyping could be performed on degraded gDNA amplified using our technique but failed with multiple displacement amplification. Of the 18 different single base mutations 16 (89.5%) were correctly identified by sequencing gDNA amplified from clinical samples using our technique. This simple and efficient isothermal method should be helpful for genetic research and clinical and forensic applications.     

Validation of whole genome amplification for analysis of the p53 tumor suppressor gene in limited amounts of tumor samples

Personalized cancer treatment requires molecular characterization of individual tumor biopsies. These samples are frequently only available in limited quantities hampering genomic analysis. Several whole genome amplification (WGA) protocols have been developed with reported varying representation of genomic regions post amplification. In this study we investigate region dropout using a φ29 polymerase based WGA approach. DNA from 123 lung cancers specimens and corresponding normal tissue were used and evaluated by Sanger sequencing of the p53 exons 5-8. To enable comparative analysis of this scarce material, WGA samples were compared with unamplified material using a pooling strategy of the 123 samples. In addition, a more detailed analysis of exon 7 amplicons were performed followed by extensive cloning and Sanger sequencing. Interestingly, by comparing data from the pooled samples to the individually sequenced exon 7, we demonstrate that mutations are more easily recovered from WGA pools and this was also supported by simulations of different sequencing coverage. Overall this data indicate a limited random loss of genomic regions supporting the use of whole genome amplification for genomic analysis.     

Whole genome amplification of sodium bisulfite-treated DNA allows the accurate estimate of methylated cytosine density in limited DNA resources

Sodium bisulfite modification-based fine mapping of methylated cytosines represents the gold standard technique for DNA methylation studies. A major problem with this approach, however is that it results in considerable DNA degradation, and large quantities of genomic DNA material are needed if numerous genomic regions are to be profiled. In this study, we examined whether whole genome amplification (WGA) techniques can be applied to sodium bisulfite-treated DNA and whether WGA would bias DNA methylation results. Sodium bisulfite-treated DNA was amplified using a standard WGA method: optimized primer-extension preamplification (PEP) with degenerate primers. Following the PCR of bisulfite-treated DNA, the DNA methylation profiles of specific DNA fragments were assessed using three approaches: (i) direct sequencing of the overall product; (ii) the sequencing of cloned PCR products; and (iii) methylation-sensitive single nucleotide primer extension (MS-SNuPE)--and compared with those obtained from bisulfite-treated DNA not subjected to WGA. Our data indicates that the DNA methylation profiles obtained from WGA of sodium bisulfite-treated DNA are consistent with those obtained from non-WGA DNA. The average difference in methylation percentage calculated from the two sets of template using MS-SNuPE was 4%. If our results are replicated on other genomic loci, WGA may become a useful technique in DNA methylation studies.     

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.     

不同单细胞全基因组扩增方法的比较及MALBAC在辅助生殖中的应用

单细胞全基因组扩增(whole genome amplification, WGA)是指在单细胞水平对全基因组进行扩增的新技术,其原理是将分离的单个细胞的微量全基因组DNA进行扩增,获得高覆盖率的完整的基因组后进行高通量测序,用于揭示细胞异质性。目前,WGA方法主要包括引物延伸预扩增(primer extension preamplification PCR, PEP-PCR)、简并寡核苷酸引物PCR (degenerate oligonucleotide primed PCR, DOP-PCR)、多重置换扩增(multiple displacement amplification, MDA)、多次退火环状循环扩增(multiple annealing and looping-based amplification cycles, MALBAC)等。本文对不同的单细胞WGA方法的原理及应用情况分别进行了阐述,并对其扩增效率进行评价和比较,包括基因组覆盖度、均一性、重现性、SNV (single-nucleotide variants)和CNV (copy number variants)检测力等。综合对比不同单细胞WGA方法后发现,MALBAC的扩增均一性最高、等位基因脱扣率最低、重现性最好,且对于CNV和SNV的检测效果最好。本文还阐述了MALBAC技术在人类单精子减数重组、非整倍体分析以及人类卵细胞基因组研究中的应用。     

Single-nucleotide-polymorphism genotyping for whole-genome-amplified samples using automated fluorescence correlation spectroscopy

Whole-genome amplification (WGA) methods were adopted for single-nucleotide-polymorphism (SNP) typing to minimize the amount of genomic DNA that has to be used in typing for thousands of different SNPs in large-scale studies; 5-10 ng of genomic DNA was amplified by a WGA method (improved primer-extension-preamplification-polymerase chain reaction (I-PEP-PCR), degenerated oligonucleotide primer-PCR (DOP-PCR), or multiple displacement amplification (MDA)). Using 1/100 to 1/500 amounts of the whole-genome-amplified products as templates, subsequent analyses were successfully performed. SNPs were genotyped by the sequence-specific primer (SSP)-PCR method followed by fluorescence correlation spectroscopy (FCS). The typing results were evaluated for four different SNPs on tumor necrosis factor receptor 1 and 2 genes (TNFR1 and TNFR2). The genotypes determined by the SSP-FCS method using the WGA products were 100% in concordance with those determined by nucleotide sequencing using genomic DNAs. We have already carried out typing of more than 300 different SNPs and are currently performing 7,500-10,000 typings per day using WGA samples from patients with several common diseases. WGA coupled with FCS allows specific and high-throughput genotyping of thousands of samples for thousands of different SNPs.     

Discovery and experimental analysis of microsatellites in an oil woody plant Camellia chekiangoleosa

Oil camellia trees are important woody plants for the production of high-quality cooking oil. On the contrary to their economic importance, their genetic and genomic resources are very limited, which greatly hamper the genetic studies on oil camellia trees. Microsatellites or simple sequence repeats (SSRs) have great value in many aspects of genetic analyses due to their high polymorphism and codominant inheritance. In this study, we report the large-scale development and characterization of SSR markers derived from genomic sequences of Camellia chekiangoleosa by high-throughput pyrosequencing technology. A total of 1,091,393 genomic shotgun reads were generated using Roche 454 FLX sequencer, the average read length was 319 bp, and the total sequence throughput was 347.9 Mb. These sequences were assembled into 35,315 contigs with total length of 14.8 Mb and the N50 contig size of 770 bp. By analyzing with microsatellite (MISA), a total of 5,844 perfect microsatellites were detected from the assembled sequences. Among them, tetranucleotide repeats were found to be the most frequent microsatellites in the genome of C. chekiangoleosa, and all the dominant repeat motifs for different types of SSRs were detected to be rich in A/T. Experimental analysis with 900 SSR primer pairs revealed that 66 % of them succeeded in PCR amplification. Further investigation with 345 SSR primer pairs showed that a relatively high percentage of primers amplified polymorphic loci (31.9 %). Experimental data also revealed that, overall, long microsatellite repeats (>20 bp) were more variable than the short ones (<20 bp) in the genome of oil camellia tree.     

The Performance of Whole Genome Amplification Methods and Next-Generation Sequencing for Pre-Implantation Genetic Diagnosis of Chromosomal Abnormalities

Reliable and accurate pre-implantation genetic diagnosis(PGD) of patient’s embryos by next-generation sequencing(NGS) is dependent on efficient whole genome amplification(WGA) of a representative biopsy sample. However, the performance of the current state of the art WGA methods has not been evaluated for sequencing. Using low template DNA(15 pg) and single cells, we showed that the two PCR-based WGA systems Sure Plex and MALBAC are superior to the REPLI-g WGA multiple displacement amplification(MDA) system in terms of consistent and reproducible genome coverage and sequence bias across the 24 chromosomes, allowing better normalization of test to reference sequencing data. When copy number variation sequencing(CNV-Seq) was applied to single cell WGA products derived by either Sure Plex or MALBAC amplification, we showed that known disease CNVs in the range of 3e15 Mb could be reliably and accurately detected at the correct genomic positions. These findings indicate that our CNV-Seq pipeline incorporating either Sure Plex or MALBAC as the key initial WGA step is a powerful methodology for clinical PGD to identify euploid embryos in a patient’s cohort for uterine transplantation.     

A Whole Genome Amplification Method to Generate Long Fragments from Low Quantities of Genomic DNA

Several whole genome amplification strategies have been developed to preamplify the entire genome from minimal amounts of DNA for subsequent molecular genetic analysis. However, none of these techniques has proven to amplify long products from very low (nanogram or picogram) quantities of genomic DNA. Here we report a new whole genome amplification protocol using a degenerate primer (DOP-PCR) that generates products up to about 10 kb in length from less than 1 ng genomic template DNA. This new protocol (LL-DOP-PCR) allows in the subsequent PCR the specific amplification, with high fidelity, of DNA fragments that are more than 1 kb in length. LL-DOP-PCR provides significantly better coverage for microsatellites and unique sequences in comparison to a conventional DOP-PCR method.     

The arbuscular mycorrhizal fungus Glomus intraradices is haploid and has a small genome size in the lower limit of eukaryotes

The genome size, complexity, and ploidy of the arbuscular mycorrhizal fungus (AMF) Glomus intraradices was determined using flow cytometry, reassociation kinetics, and genomic reconstruction. Nuclei of G. intraradices from in vitro culture, were analyzed by flow cytometry. The estimated average length of DNA per nucleus was 14.07+/-3.52 Mb. Reassociation kinetics on G. intraradices DNA indicated a haploid genome size of approximately 16.54 Mb, comprising 88.36% single copy DNA, 1.59% repetitive DNA, and 10.05% fold-back DNA. To determine ploidy, the DNA content per nucleus measured by flow cytometry was compared with the genome estimate of reassociation kinetics. G. intraradices was found to have a DNA index (DNA per nucleus per haploid genome size) of approximately 0.9, indicating that it is haploid. Genomic DNA of G. intraradices was also analyzed by genomic reconstruction using four genes (Malate synthase, RecA, Rad32, and Hsp88). Because we used flow cytometry and reassociation kinetics to reveal the genome size of G. intraradices and show that it is haploid, then a similar value for genome size should be found when using genomic reconstruction as long as the genes studied are single copy. The average genome size estimate was 15.74+/-1.69 Mb indicating that these four genes are single copy per haploid genome and per nucleus of G. intraradices. Our results show that the genome size of G. intraradices is much smaller than estimates of other AMF and that the unusually high within-spore genetic variation that is seen in this fungus cannot be due to high ploidy.     

Genome amplification of single sperm using multiple displacement amplification

Sperm typing is an effective way to study recombination rate on a fine scale in regions of interest. There are two strategies for the amplification of single meiotic recombinants: repulsion-phase allele-specific PCR and whole genome amplification (WGA). The former can selectively amplify single recombinant molecules from a batch of sperm but is not scalable for high-throughput operation. Currently, primer extension pre-amplification is the only method used in WGA of single sperm, whereas it has limited capacity to produce high-coverage products enough for the analysis of local recombination rate in multiple large regions. Here, we applied for the first time a recently developed WGA method, multiple displacement amplification (MDA), to amplify single sperm DNA, and demonstrated its great potential for producing high-yield and high-coverage products. In a 50 μl reaction, 76 or 93% of loci can be amplified at least 2500- or 250-fold, respectively, from single sperm DNA, and second-round MDA can further offer >200-fold amplification. The MDA products are usable for a variety of genetic applications, including sequencing and microsatellite marker and single nucleotide polymorphism (SNP) analysis. The use of MDA in single sperm amplification may open a new era for studies on local recombination rates.     

The use of whole genome amplification in the study of human disease

The availability of large amounts of genomic DNA is of critical importance for many of the molecular biology assays used in the analysis of human disease. However, since the amount of patient tissue available is often limited and as particular foci of interest may consist of only a few hundred cells, the yield of DNA is often insufficient for extensive analysis. To address this problem, several whole genome amplification (WGA) methodologies have been developed. Initial WGA approaches were based on the polymerase chain reaction (PCR). However, recent reports have described the use of non-PCR-based linear amplification protocols for WGA. Using these methods, it is possible to generate microgram quantities of DNA starting with as little as 1mg of genomic DNA. This review will provide an overview of WGA approaches and summarize some of the uses for amplified DNA in various high-throughput genetic applications.     

Monodisperse Picoliter Droplets for Low-Bias and Contamination-Free Reactions in Single-Cell Whole Genome Amplification

Whole genome amplification (WGA) is essential for obtaining genome sequences from single bacterial cells because the quantity of template DNA contained in a single cell is very low. Multiple displacement amplification (MDA), using Phi29 DNA polymerase and random primers, is the most widely used method for single-cell WGA. However, single-cell MDA usually results in uneven genome coverage because of amplification bias, background amplification of contaminating DNA, and formation of chimeras by linking of non-contiguous chromosomal regions. Here, we present a novel MDA method, termed droplet MDA, that minimizes amplification bias and amplification of contaminants by using picoliter-sized droplets for compartmentalized WGA reactions. Extracted DNA fragments from a lysed cell in MDA mixture are divided into 10⁵ droplets (67 pL) within minutes via flow through simple microfluidic channels. Compartmentalized genome fragments can be individually amplified in these droplets without the risk of encounter with reagent-borne or environmental contaminants. Following quality assessment of WGA products from single Escherichia coli cells, we showed that droplet MDA minimized unexpected amplification and improved the percentage of genome recovery from 59% to 89%. Our results demonstrate that microfluidic-generated droplets show potential as an efficient tool for effective amplification of low-input DNA for single-cell genomics and greatly reduce the cost and labor investment required for determination of nearly complete genome sequences of uncultured bacteria from environmental samples.     

Improved Multiple Displacement Amplification (iMDA) and Ultraclean Reagents

Background

Next-generation sequencing sample preparation requires nanogram to microgram quantities of DNA; however, many relevant samples are comprised of only a few cells. Genomic analysis of these samples requires a whole genome amplification method that is unbiased and free of exogenous DNA contamination. To address these challenges we have developed protocols for the production of DNA-free consumables including reagents and have improved upon multiple displacement amplification (iMDA).

Results

A specialized ethylene oxide treatment was developed that renders free DNA and DNA present within Gram positive bacterial cells undetectable by qPCR. To reduce DNA contamination in amplification reagents, a combination of ion exchange chromatography, filtration, and lot testing protocols were developed. Our multiple displacement amplification protocol employs a second strand-displacing DNA polymerase, improved buffers, improved reaction conditions and DNA free reagents. The iMDA protocol, when used in combination with DNA-free laboratory consumables and reagents, significantly improved efficiency and accuracy of amplification and sequencing of specimens with moderate to low levels of DNA. The sensitivity and specificity of sequencing of amplified DNA prepared using iMDA was compared to that of DNA obtained with two commercial whole genome amplification kits using 10 fg (~1-2 bacterial cells worth) of bacterial genomic DNA as a template. Analysis showed >99% of the iMDA reads mapped to the template organism whereas only 0.02% of the reads from the commercial kits mapped to the template. To assess the ability of iMDA to achieve balanced genomic coverage, a non-stochastic amount of bacterial genomic DNA (1 pg) was amplified and sequenced, and data obtained were compared to sequencing data obtained directly from genomic DNA. The iMDA DNA and genomic DNA sequencing had comparable coverage 99.98% of the reference genome at ≥1X coverage and 99.9% at ≥5X coverage while maintaining both balance and representation of the genome.

Conclusions

The iMDA protocol in combination with DNA-free laboratory consumables, significantly improved the ability to sequence specimens with low levels of DNA. iMDA has broad utility in metagenomics, diagnostics, ancient DNA analysis, pre-implantation embryo screening, single-cell genomics, whole genome sequencing of unculturable organisms, and forensic applications for both human and microbial targets.     

Suitability of genomic DNA synthesized by strand displacement amplification (SDA) for AFLP analysis: genotyping single spores of arbuscular mycorrhizal (AM) fungi

Limited biological samples of microbial origin often yield insufficient amounts of genomic DNA, making application of standard techniques of genetic analysis, like amplified fragment length polymorphism (AFLP), virtually impossible. The Phi29 DNA polymerase based whole genome amplification (WGA) method has the potential to alleviate this technical bottleneck. In the present work, we have sought to investigate the suitability of genomic DNA synthesized using Phi29 based WGA for AFLP analysis. We first used genomic DNA from Saccharomyces cerevisiae to optimize the protocol for the use of SDA-amplified DNA for AFLP analysis. Based on the optimized protocol we obtained AFLP fingerprints which were indistinguishable from the non-amplified genomic DNA. Finally, AFLP analysis was performed using SDA synthesized genomic DNA from single spores of various species of arbuscular mycorrhizal (AM) fungi. Unique and highly reproducible fingerprints for each species were obtained. The present study introduces the application of WGA-mediated AFLP to AM fungal biology; similarly, our protocol could be useful for other microbial genomes currently not amenable to genetic analysis owing to the paucity of starting template.     

A novel method to perform genomic walks using a combination of single strand DNA circularization and rolling circle amplification

Characterization of regions flanking a known sequence within a genome, known as genome walking, is a cornerstone technique in modern genetic analysis. In the present work we have developed a new PCR-dependent, directional genome walking protocol based on the unique circularization property of a novel DNA ligase, CircLigase. In the first step, PCR based primer extension is performed using a phosphorylated primer, designed to extend from the boundary of the known sequence, into the flanking region. This linear amplification results in the generation of single-stranded (ss) DNA, which is then circularized using CircLigase. Using the hyperbranching activity of Phi29 DNA polymerase, the circular ssDNA is then linearized by rolling circle amplification, resulting in copious amounts of double stranded concatameric DNA. Nested primers are used to amplify the flanking sequence using inverse PCR. The products are resolved on an agarose gel and the bands whose mobility change due to the nested location of the primer combination used are identified, extracted, and cloned into a plasmid vector for sequencing. Empirical proof for this concept was generated on two antimicrobial biosynthetic genes in Pseudomonas sp. LBUM300. Using the hcnB and phlD genes as starting points, ca 1 kb of flanking sequences were successfully isolated. The use of locus specific primers ensured both directionality and specificity of the walks, alleviating the generation of spurious amplicons, typically observed in randomly primed walking protocols. The presented genome walking protocol could be applied to any microbial genome and requires only 100-150 bp of prior sequence information. The proposed methodology does not entail laborious testing of restriction enzymes or adaptor ligation. This is the first report of a successful application of the novel ligase enzyme, CircLigase for genomic walking purposes.     

Comparative analysis of microsatellite loci in chicken and turkey.

Cross-species amplification of 520 chicken microsatellite markers was tested by polymerase chain reaction with genomic DNA of the turkey (Meleagris gallopavo). Each primer pair was tested at six different combinations of annealing temperature and MgCl2 concentration. A total of 280 (54%) of the primer pairs produced amplification products. The majority of these products were similar, if not identical in size to those expected based on the fragment sizes of the corresponding chicken loci. Structure of the dinucleotide repeat and flanking sequences was examined for 13 turkey fragments (amplified with chicken primers) and 5 chicken fragments (amplified with turkey primers). Sequence analysis found a wide array of mutations between species in addition to differences in repeat length. To estimate the usefulness of the amplified loci for genetic mapping in the turkey, allelic polymorphism was determined for 57 of the 280 amplified loci. A total of 20 of 57 markers (35%) were polymorphic with an average of 1.4 alleles per locus. The results of this study suggest that approximately 20% of the chicken microsatellite markers will be useful for mapping the turkey genome.     

Genetic variability of Old Portuguese bread wheat cultivars assayed by IRAP and REMAP markers

Retrotransposons (RTNs) constitute informative molecular markers for plant species as a result of their ability of integrating into a multitude of loci throughout the genome and thereby generating insertional polymorphisms between individuals. Inter-retrotransposon amplified polymorphisms (IRAPs) and the retrotransposon-microsatellite amplified polymorphisms (REMAPs) are marker systems based on long terminal repeats (LTRs) RTNs, developed for plants, that have been widely used for evolution, genetic diversity, DNA fingerprinting of cultivars and varieties, genetic mapping linkage and for detection of genetic rearrangements induced by polyploidisation. In the present study, we aimed to analyse the genetic variability among 48 Old Portuguese bread wheat cultivars using both IRAP and REMAP markers. Five IRAP and six REMAP primer combinations were used. IRAP produced 103 polymorphic fragments in a total of 113 bands. On average, 22.6 bands were amplified per IRAP primer combination. The bands ranged in size from 250 to 5000 bp. The REMAP primer combinations allowed the amplification of 53 bands, 51 of them polymorphic. An average of 8.8 REMAP bands was scored per primer combination. The REMAP bands ranged from 250 to 3000 bp. Both marker systems presented high percentages of polymorphism. However, IRAP markers were suitable for detecting genetic variability at the individual level and did not differentiate higher taxa. The REMAP maker system allowed the clustering by botanical variety and identified most of the homonym bread wheat cultivars.     

Examination of host genome for the presence of integrated fragments of Solenopsis invicta virus 1

A series of oligonucleotide primer pairs covering the entire genome of Solenopsis invicta virus 1 (SINV-1) were used to probe the genome of its host, S. invicta, for integrated fragments of the viral genome. All of the oligonucleotide primer sets yielded amplicons of anticipated size from cDNA created from an RNA template from SINV-1. However, no corresponding amplification was observed when genomic DNA (from 32 colonies of S. invicta) was used as template for the PCR amplifications. Host DNA integrity was verified by amplification of an ant-specific gene, SiGSTS1. The representation of fire ant colonies included both social forms, monogyne and polygyne, and those infected and uninfected with SINV-1. Furthermore, no amplification was observed from genomic DNA from ant samples collected from Argentina or the US. Thus, it appears that SINV-1 genome integration, or a portion therein, has not likely occurred within the S. invicta host genome.