Single-molecule dilution and multiple displacement amplification for molecular haplotyping

"Evaluating the Benefits and Limits of Multiple Displacement Amplification with Whole-Genome Oxford Nanopore Sequencing"

Multiple Displacement Amplification (MDA) outperforms conventional PCR in long fragment and whole genome amplification which makes it attractive to couple with long-read sequencing of samples with limited quantities of DNA to obtain improved genome assemblies. Here, we explore the efficacy and limits of MDA for genome sequence assembly using Oxford Nanopore Technologies (ONT) rapid library preparations and minION sequencing. We successfully generated almost complete genome sequences for all organisms examined, including Cryptosporidium meleagridis, Staphylococcus aureus, Enterococcus faecium, and Escherichia coli, with the ability to generate high-quality data from samples starting with only 0.025 ng of total DNA. Controlled sheared DNA samples exhibited a distinct pattern of size-increase after MDA, which may be associated with the amplification of long, low-abundance fragments present in the assay, as well as generating concatemeric sequences during amplification. To address concatemers, we developed a computational pipeline (CADECT: Concatemer Detection Tool) to identify and remove putative concatemeric sequences. This study highlights the efficacy of MDA in generating high-quality genome assemblies from limited amounts of input DNA. Also, the CADECT pipeline effectively mitigated the impact of concatemeric sequences, enabling the assembly of contiguous sequences even in cases where the input genomic DNA was degraded. These results have significant implications for the study of organisms that are challenging to culture in vitro, such as Cryptosporidium, and for expediting critical results in clinical settings with limited quantities of available genomic DNA.

Phasing DNA Markers Using Digital PCR

Besides quantifying the absolute number of copies of known DNA targets, digital PCR can also be used to assess whether two nonpolymorphic gene sequences or two heterozygous markers reside on the same DNA molecule (i.e., are physically linked). Some useful linkage applications include: phasing variants to define a haplotype; genotyping of inversions; determining the presence of multimarker pathogenic bacteria in a metagenomic sample; and assessing DNA integrity. This chapter describes an efficient and cost-effective method for analyzing linkage of any two genetic sequences up to at least 200 Kb apart, including phasing of heterozygous markers such as that which occur abundantly in the cystic fibrosis transmembrane conductance regulator (CFTR) gene.

Haplotype-Contained PCR Products Analysis by Sequencing with Selective Restriction of Primer Extension

We develop a strategy for haplotype analysis of PCR products that contained two adjacent heterozygous loci using sequencing with specific primers, allele-specific primers, and ddNTP-blocked primers. To validate its feasibility, two sets of PCR products, including two adjacent heterozygous SNPs, UGT1A1⁎6 (rs4148323) and UGT1A1⁎28 (rs8175347), and two adjacent heterozygous SNPs, K1637K (rs11176013) and S1647T (rs11564148), were analyzed. Haplotypes of PCR products, including UGT1A1⁎6 and UGT1A1⁎28, were successfully analyzed by Sanger sequencing with allele-specific primers. Also, haplotypes of PCR products, including K1637K and S1647T, could not be determined by Sanger sequencing with allele-specific primers but were successfully analyzed by pyrosequencing with ddNTP-blocked primers. As a result, this method is able to effectively haplotype two adjacent heterozygous PCR products. It is simple, fast, and irrespective of short read length of pyrosequencing. Overall, we fully hope it will provide a new promising technology to identify haplotypes of conventional PCR products in clinical samples.

Quantitative haplotyping of PCR products by nonsynchronous pyrosequencing with di-base addition

Molecular haplotyping is becoming increasingly important for studying the disease association of a specific allele because of its ability of providing more information than any single nucleotide polymorphism (SNP). Computational analysis and experimental techniques are usually performed for haplotypic determination. However, established methods are not suitable for analyzing haplotypes of massive natural DNA samples. Here we present a simple molecular approach to analyze haplotypes of conventional polymerase chain reaction (PCR) products quantitatively in a single sequencing run. In this approach, specific types and proportions of haplotypes in both individual and pooled samples could be determined by solving equations constructed from nonsynchronous pyrosequencing with di-base addition. Two SNPs (rs11176013 and rs11564148) in the gene for leucine-rich repeat kinase 2 (LRRK2) related to Parkinson's disease were selected as experimental sites. A series of DNA samples, including these two heterozygous loci, were investigated. This approach could accurately identify multiple DNA samples indicating that the approach is likely to be applied for haplotyping of unrestricted conventional PCR products from natural samples, and be especially applicable for analyzing short sequences in clinical diagnosis. Graphical Abstract One DNA sample consisting of 4 different DNA templates with different proportion are sequenced by nonsynchronous pyrosequencing with di-base addition. The number of incorporated nucleotides produced by a single sequencing reaction equals to the total of incorporated nucleotides. Four independent equations are constructed from the pyrograms of nonsynchronous pyrosequencing data. Molecular haplotypes of two adjacent SNPs can be quantitatively identified by solving these equations.

Haplotype-resolved genome sequencing: experimental methods and applications

Human genomes are diploid and, for their complete description and interpretation, it is necessary not only to discover the variation they contain but also to arrange it onto chromosomal haplotypes. Although whole-genome sequencing is becoming increasingly routine, nearly all such individual genomes are mostly unresolved with respect to haplotype, particularly for rare alleles, which remain poorly resolved by inferential methods. Here, we review emerging technologies for experimentally resolving (that is, 'phasing') haplotypes across individual whole-genome sequences. We also discuss computational methods relevant to their implementation, metrics for assessing their accuracy and completeness, and the relevance of haplotype information to applications of genome sequencing in research and clinical medicine.

A rapid molecular approach for chromosomal phasing

Determining the chromosomal phase of pairs of sequence variants - the arrangement of specific alleles as haplotypes - is a routine challenge in molecular genetics. Here we describe Drop-Phase, a molecular method for quickly ascertaining the phase of pairs of DNA sequence variants (separated by 1-200 kb) without cloning or manual single-molecule dilution. In each Drop-Phase reaction, genomic DNA segments are isolated in tens of thousands of nanoliter-sized droplets together with allele-specific fluorescence probes, in a single reaction well. Physically linked alleles partition into the same droplets, revealing their chromosomal phase in the co-distribution of fluorophores across droplets. We demonstrated the accuracy of this method by phasing members of trios (revealing 100% concordance with inheritance information), and demonstrate a common clinical application by phasing CFTR alleles at genomic distances of 11-116 kb in the genomes of cystic fibrosis patients. Drop-Phase is rapid (requiring less than 4 hours), scalable (to hundreds of samples), and effective at long genomic distances (200 kb).

Whole-genome haplotyping by dilution, amplification, and sequencing

Standard whole-genome genotyping technologies are unable to determine haplotypes. Here we describe a method for rapid and cost-effective long-range haplotyping. Genomic DNA is diluted and distributed into multiple aliquots such that each aliquot receives a fraction of a haploid copy. The DNA template in each aliquot is amplified by multiple displacement amplification, converted into barcoded sequencing libraries using Nextera technology, and sequenced in multiplexed pools. To assess the performance of our method, we combined two male genomic DNA samples at equal ratios, resulting in a sample with diploid X chromosomes with known haplotypes. Pools of the multiplexed sequencing libraries were subjected to targeted pull-down of a 1-Mb contiguous region of the X-chromosome Duchenne muscular dystrophy gene. We were able to phase the Duchenne muscular dystrophy region into two contiguous haplotype blocks with a mean length of 494 kb. The haplotypes showed 99% agreement with the consensus base calls made by sequencing the individual DNAs. We subsequently used the strategy to haplotype two human genomes. Standard genomic sequencing to identify all heterozygous SNPs in the sample was combined with dilution-amplification-based sequencing data to resolve the phase of identified heterozygous SNPs. Using this procedure, we were able to phase >95% of the heterozygous SNPs from the diploid sequence data. The N50 for a Yoruba male DNA was 702 kb whereas the N50 for a European female DNA was 358 kb. Therefore, the strategy described here is suitable for haplotyping of a set of targeted regions as well as of the entire genome.

The importance of phase information for human genomics

Contemporary sequencing studies often ignore the diploid nature of the human genome because they do not routinely separate or 'phase' maternally and paternally derived sequence information. However, many findings - both from recent studies and in the more established medical genetics literature - indicate that relationships between human DNA sequence and phenotype, including disease, can be more fully understood with phase information. Thus, the existing technological impediments to obtaining phase information must be overcome if human genomics is to reach its full potential.

Improved multiple displacement amplification with phi29 DNA polymerase for genotyping of single human cells

The ability to genotype multiple loci of single cells would be of significant benefit to investigations of cellular processes such as oncogenesis, meiosis, fertilization, and embryogenesis. We report a simple two-step, single-tube protocol for whole-genome amplification (WGA) from single human cells using components of the GenomiPhi V2 DNA Amplification kit. For the first time, we demonstrate reliable generation of 4-7 microg amplified DNA from a single human cell within 4 h with a minimum amount of artifactual DNA synthesis. DNA amplified from single cells was genotyped for 13 heterozygous short tandem repeats (STRs) and 7 heterozygous single nucleotide polymorphisms (SNPs), and the genotyping results were compared with purified genomic DNA. Accuracy of genotyping (percent of single-cell amplifications genotyped accurately for any particular STR or SNP) varied from 37% to 100% (with an average of 80%) for STRs and from 89% to 100% (averaging 94%) for SNPs. We suggest that the method described in this report is suitable for WGA from single cells, the product of which can be subsequently used for many applications, such as preimplantation genetic analysis (PGD).

Single molecule transcription profiling with AFM

Established techniques for global gene expression profiling, such as microarrays, face fundamental sensitivity constraints. Due to greatly increasing interest in examining minute samples from micro-dissected tissues, including single cells, unorthodox approaches, including molecular nanotechnologies, are being explored in this application. Here, we examine the use of single molecule, ordered restriction mapping, combined with AFM, to measure gene transcription levels from very low abundance samples. We frame the problem mathematically, using coding theory, and present an analysis of the critical error sources that may serve as a guide to designing future studies. We follow with experiments detailing the construction of high density, single molecule, ordered restriction maps from plasmids and from cDNA molecules, using two different enzymes, a result not previously reported. We discuss these results in the context of our calculations.

Whole genome amplification from single cells in preimplantation genetic diagnosis and prenatal diagnosis

The literature on whole genome amplification (WGA) techniques and their application to preimplantation genetic diagnosis (PGD) and prenatal diagnosis is reviewed. General polymerase chain reaction (PCR) fails to provide adequate information from limited cells in PGD and non-invasive prenatal diagnosis. Therefore several WGA techniques, such as primer extension preamplification (PEP) and degenerate oligonucleotide primed PCR (DOP-PCR), have been developed and successfully applied to clinical work during the past decade, especially in PGD and prenatal diagnosis. These techniques can provide ample amplification of genetic sequences from single cells for a series of subsequent PCR analyses such as restriction fragment length polymorphisms (RFLP) and comparative genomic hybridization (CGH), thus opening up a new area for prenatal diagnosis. However, several problems have been reported in the application of these techniques. The ideal WGA technique should have high yield, faithful representation of the original template, complete coverage of the genome, and simply performed procedure. In order to make good use of these techniques in future research and clinical work, it is undoubtedly necessary for an extensive understanding of the merits and pitfalls of these recently developed techniques.

Comparison of two whole genome amplification methods for STR genotyping of LCN and degraded DNA samples

The analysis of LCN or highly degraded DNA samples presents a challenge for forensic science. Improving the quantity and/or quality of samples would greatly increase the profiling success rate from LCN and degraded samples. Whole genome amplification (WGA) is one method that has such potential. Two commercially available WGA kits, GenomePlex and GenomiPhi, were investigated for use on LCN and degraded DNA samples. Both kits amplified genomic DNA, producing microgram quantities from sub-nanogram templates. Profiling success of LCN DNA samples was increased, with improvements of over 700% from 10pg template DNA compared to non-WGA-amplified control samples. The amplification success with degraded DNA was also improved by WGA. Degraded DNA was simulated using restriction enzymes to demonstrate that the application of WGA can result in the typing of STR loci that could not previously be amplified. An increase in artefacts, such as stutter alleles and amplification biases, were observed in many samples. Results show that WGA is capable of increasing both the quality and quantity of DNA, and has the potential to improve profiling success from difficult samples in forensic casework.

Ancient DNA forces reconsideration of evolutionary history of Mediterranean pygmy elephantids

During the Pleistocene pygmy elephantids, some only a quarter of their ancestors' size, were present on Mediterranean islands until about 10,000 years ago (y.a.). Using a new methodology for ancient DNA (aDNA) studies, the whole genomic multiple displacement amplification method, we were able to retrieve cytochrome b (cytb) DNA fragments from 4200 to 800,000 y.a. specimens from island and mainland samples, including pygmy and normal-sized forms. The short DNA sequence (43 bp) retrieved from the 800,000 y.a. sample is one of the oldest DNA fragment ever retrieved. Duplication of the experiments in two laboratories, the occurrence of three diagnostic sites and the results of the phylogenetic analyses strongly support its authenticity. Our results challenge the prevailing view that pygmy elephantids of the eastern Mediterranean originated exclusively from Elephas, suggesting independent histories of dwarfism and the presence of both pygmy mammoths and elephant-like taxa on these islands. Based on our molecular data, the origin of the Tilos and Cyprus elephantids from a lineage within the genus Elephas is confirmed, while the DNA sequence from the Cretan sample falls clearly within the mammoth clade. Thus, the name Mammuthus creticus rather than Elephas creticus, seems to be justified for this form. Our findings also suggest a need to re-evaluate the evolutionary history of the Sicilian/Maltese species, traditionally included in the genus Elephas.

An efficient method for multi-locus molecular haplotyping

Many methods exist for genotyping—revealing which alleles an individual carries at different genetic loci. A harder problem is haplotyping—determining which alleles lie on each of the two homologous chromosomes in a diploid individual. Conventional approaches to haplotyping require the use of several generations to reconstruct haplotypes within a pedigree, or use statistical methods to estimate the prevalence of different haplotypes in a population. Several molecular haplotyping methods have been proposed, but have been limited to small numbers of loci, usually over short distances. Here we demonstrate a method which allows rapid molecular haplotyping of many loci over long distances. The method requires no more genotypings than pedigree methods, but requires no family material. It relies on a procedure to identify and genotype single DNA molecules, and reconstruction of long haplotypes by a ‘tiling’ approach. We demonstrate this by resolving haplotypes in two regions of the human genome, harbouring 20 and 105 single-nucleotide polymorphisms, respectively. The method can be extended to reconstruct haplotypes of arbitrary complexity and length, and can make use of a variety of genotyping platforms. We also argue that this method is applicable in situations which are intractable to conventional approaches.

Multiple displacement amplification to create a long-lasting source of DNA for genetic studies

In many situations there may not be sufficient DNA collected from patient or population cohorts to meet the requirements of genome-wide analysis of SNPs, genomic copy number polymorphisms, or acquired copy number alternations. When the amount of available DNA for genotype analysis is limited, high performance whole-genome amplification (WGA) represents a new development in genetic analysis. It is especially useful for analysis of DNA extracted from stored histology slides, tissue samples, buccal swabs, or blood stains collected on filter paper. The multiple displacement amplification (MDA) method, which relies on isothermal amplification using the DNA polymerase of the bacteriophage phi29, is a recently developed technique for high performance WGA. This review addresses new trends in the technical performance of MDA and its applications to genetic analyses. The main challenge of WGA methods is to obtain balanced and faithful replication of all chromosomal regions without the loss of or preferential amplification of any genomic loci or allele. In multiple comparisons to other WGA methods, MDA appears to be most reliable for genotyping, with the most favorable call rates, best genomic coverage, and lowest amplification bias.

Long-range multilocus haplotype phasing of the MHC

Haplotypes are a powerful tool for identifying the genetic basis of common complex diseases. Disease-association mapping requires molecular methods for haplotyping biallelic SNP variation and highly complex polymorphisms. We developed a method for phasing HLA-A, HLA-B, and HLA-DRB1 alleles on chromosome 6 in unrelated individuals. This method uses the highly polymorphic HLA-B locus to discriminate the two HLA haplotypes in heterozygous individuals and its ideal location 1.4 Mbp telomeric to HLA-DRB1 and 1.2 Mbp centromeric to HLA-A to capture 2-Mbp-long genomic DNA. Genomic DNA representing a single HLA-B-captured haplotype is genotyped for HLA-A and HLA-DRB1 alleles and linkage to HLA-B is established. Proof of principle was established in a large blinded study of phase-known samples. Availability of an efficient method for MHC haplotype phase determination will facilitate the mapping of causative MHC-resident genes in many human diseases and has the potential to be broadened to other polymorphic gene complexes.

[11] Sample Labeling: An Overview

It is not easy to write a critical review of the methods available for labeling RNA and DNA "extracts" for microarray studies. There are a number of reasons for this: Suppliers of the reagents and kits used for this purpose do research and development, quality control, and validation and then they provide a hard-wired, "optimized" product. They often give few details about the compositions of these products, are inclined to put the best face they can on what they sell and gloss over any deficiencies, and have no interest in paying for direct comparisons of their product to those of other companies. These comparisons can be expensive to perform, and there are few good examples in the literature. When comparative experiments have been done, it is not clear that each of the individual methods tested was executed with equal proficiency. Many experiments can be required to determine how best to hybridize any given labeled extract to a particular array and how to block, wash, and postprocess (e.g., stain) the array so that the signal-to-noise ratio is maximized. In addition, authors of comparative studies used different arrays, technical protocols (some of which are out of date), experimental designs, and analyses. Finally, some new techniques, which seem quite promising, have been employed so little that their strengths and shortcomings are difficult to assess.