Resolving the full spectrum of human genome variation using Linked-Reads.

Large-scale population analyses coupled with advances in technology have demonstrated that the human genome is more diverse than originally thought. To date, this diversity has largely been uncovered using short-read whole-genome sequencing. However, these short-read approaches fail to give a complete picture of a genome. They struggle to identify structural events, cannot access repetitive regions, and fail to resolve the human genome into haplotypes. Here, we describe an approach that retains long range information while maintaining the advantages of short reads. Starting from ∼1 ng of high molecular weight DNA, we produce barcoded short-read libraries. Novel informatic approaches allow for the barcoded short reads to be associated with their original long molecules producing a novel data type known as "Linked-Reads". This approach allows for simultaneous detection of small and large variants from a single library. In this manuscript, we show the advantages of Linked-Reads over standard short-read approaches for reference-based analysis. Linked-Reads allow mapping to 38 Mb of sequence not accessible to short reads, adding sequence in 423 difficult-to-sequence genes including disease-relevant genes STRC, SMN1, and SMN2 Both Linked-Read whole-genome and whole-exome sequencing identify complex structural variations, including balanced events and single exon deletions and duplications. Further, Linked-Reads extend the region of high-confidence calls by 68.9 Mb. The data presented here show that Linked-Reads provide a scalable approach for comprehensive genome analysis that is not possible using short reads alone.

Haplotyping germline and cancer genomes with high-throughput linked-read sequencing.

Haplotyping of human chromosomes is a prerequisite for cataloguing the full repertoire of genetic variation. We present a microfluidics-based, linked-read sequencing technology that can phase and haplotype germline and cancer genomes using nanograms of input DNA. This high-throughput platform prepares barcoded libraries for short-read sequencing and computationally reconstructs long-range haplotype and structural variant information. We generate haplotype blocks in a nuclear trio that are concordant with expected inheritance patterns and phase a set of structural variants. We also resolve the structure of the EML4-ALK gene fusion in the NCI-H2228 cancer cell line using phased exome sequencing. Finally, we assign genetic aberrations to specific megabase-scale haplotypes generated from whole-genome sequencing of a primary colorectal adenocarcinoma. This approach resolves haplotype information using up to 100 times less genomic DNA than some methods and enables the accurate detection of structural variants.

Real-time DNA sequencing from single polymerase molecules.

Pacific Biosciences has developed a method for real-time sequencing of single DNA molecules (Eid et al., 2009), with intrinsic sequencing rates of several bases per second and read lengths into the kilobase range. Conceptually, this sequencing approach is based on eavesdropping on the activity of DNA polymerase carrying out template-directed DNA polymerization. Performed in a highly parallel operational mode, sequential base additions catalyzed by each polymerase are detected with terminal phosphate-linked, fluorescence-labeled nucleotides. This chapter will first outline the principle of this single-molecule, real-time (SMRT) DNA sequencing method, followed by descriptions of its underlying components and typical sequencing run conditions. Two examples are provided which illustrate that, in addition to the DNA sequence, the dynamics of DNA polymerization from each enzyme molecules is directly accessible: the determination of base-specific kinetic parameters from single-molecule sequencing reads, and the characterization of DNA synthesis rate heterogeneities.

Real-time DNA sequencing from single polymerase molecules.

We present single-molecule, real-time sequencing data obtained from a DNA polymerase performing uninterrupted template-directed synthesis using four distinguishable fluorescently labeled deoxyribonucleoside triphosphates (dNTPs). We detected the temporal order of their enzymatic incorporation into a growing DNA strand with zero-mode waveguide nanostructure arrays, which provide optical observation volume confinement and enable parallel, simultaneous detection of thousands of single-molecule sequencing reactions. Conjugation of fluorophores to the terminal phosphate moiety of the dNTPs allows continuous observation of DNA synthesis over thousands of bases without steric hindrance. The data report directly on polymerase dynamics, revealing distinct polymerization states and pause sites corresponding to DNA secondary structure. Sequence data were aligned with the known reference sequence to assay biophysical parameters of polymerization for each template position. Consensus sequences were generated from the single-molecule reads at 15-fold coverage, showing a median accuracy of 99.3%, with no systematic error beyond fluorophore-dependent error rates.

A rapid fluorescence assay for FtsZ assembly indicates cooperative assembly with a dimer nucleus.

FtsZ is the major cytoskeletal protein operating in bacterial cell division. FtsZ assembles into protofilaments in vitro, and there has been some controversy over whether the assembly is isodesmic or cooperative. Assembly has been assayed previously by sedimentation and light scattering. However, these techniques will under-report small polymers. We have now produced a mutant of Escherichia coli FtsZ, L68W, which gives a 250% increase in tryptophan fluorescence upon polymerization. This provides a real-time assay of polymer that is directly proportional to the concentration of subunit interfaces. FtsZ-L68W is functional for cell division, and should therefore be a valid model for studying the thermodynamics and kinetics of FtsZ assembly. We assayed assembly at pH 7.7 and pH 6.5, in 2.5 mM EDTA. EDTA blocks GTP hydrolysis and should give an assembly reaction that is not complicated by the irreversible hydrolysis step. Assembly kinetics was determined with a stopped-flow device for a range of FtsZ concentrations. When assembly was initiated by adding 0.2 mM GTP, fluorescence increase showed a lag, followed by nucleation, elongation, and a plateau. The assembly curves were fit to a cooperative mechanism that included a monomer activation step, a weak dimer nucleus, and elongation. Fragmentation was absent in the model, another characteristic of cooperative assembly. We are left with an enigma: how can the FtsZ protofilament, which appears to be one-subunit thick, assemble with apparent cooperativity?

Differential and simultaneous adenosine di- and triphosphate binding by MutS.

The roles of ATP binding and hydrolysis in the function of MutS in mismatch repair are poorly understood. As one means of addressing this question, we have determined the affinities and number of adenosine di- and triphosphate binding sites within MutS. Nitrocellulose filter binding assay and equilibrium fluorescence anisotropy measurements have demonstrated that MutS has one high affinity binding site for ADP and one high affinity site for nonhydrolyzable ATP analogues per dimer equivalent. Low concentrations of 5'-adenylylimidodiphosphate (AMPPNP) promote ADP binding and a large excess of AMPPNP is required to displace ADP from the protein. Fluorescence energy transfer and filter binding assays indicate that ADP and nonhydrolyzable ATP analogues can bind simultaneously to adjacent subunits within the MutS oligomer with affinities in the low micromolar range. These findings suggest that the protein exists primarily as the ATP.MutS.ADP ternary complex in solution and that this may be the form of the protein that is involved in DNA encounters in vivo.

Assembly and molecular activities of the MutS tetramer.

Analytical equilibrium ultracentrifugation indicates that Escherichia coli MutS exists as an equilibrating mixture of dimers and tetramers. The association constant for the dimer-to-tetramer transition is 2.1 x 10(7) M-1, indicating that the protein would consist of both dimers and tetramers at physiological concentrations. The carboxyl terminus of MutS is required for tetramer assembly because a previously described 53-amino acid carboxyl-terminal truncation (MutS800) forms a limiting species of a dimer (Obmolova, G., Ban, C., Hsieh, P., and Yang, W. (2000) Nature 407, 703-710; Lamers, M. H., Perrakis, A., Enzlin, J. H., Winterwerp, H. H., de Wind, N., and Sixma, T. K. (2000) Nature 407, 711-717). MutS800 binds a 20-base pair heteroduplex an order of magnitude more weakly than full-length MutS, and at saturating protein concentrations, the heteroduplex-bound mass observed with MutS800 is only half that observed with the full length protein, indicating that the subunit copy number of heteroduplex-bound MutS is twice that of MutS800. Analytical equilibrium ultracentrifugation using a fluorescein-tagged 20-base pair heteroduplex demonstrated that native MutS forms a tetramer on this single site-sized heteroduplex DNA. Equilibrium fluorescence experiments indicated that dimer-to-tetramer assembly promotes mismatch binding by MutS and that the tetramer can bind only a single heteroduplex molecule, implying nonequivalence of the two dimers within the tetramer. Compared with native MutS, the ability of MutS800 to promote MutL-dependent activation of MutH is substantially reduced.