Specific-primer-directed DNA sequencing using automated fluorescence detection.

Automated fluorescence-based DNA sequence analysis offers the possibility to undertake very large scale sequencing projects. Directed strategies, such as the specific-primer-directed sequencing approach ('gene walking'), should prove useful in such projects. Described herein is a study involving the use of this approach in conjunction with automated fluorescence detection on a commercial instrument (ABI 370A DNA sequencer). This includes procedures for the rapid chemical synthesis and purification of labeled primers, the design of primer sequences that are compatible with the commercial analysis software, and automated DNA sequence analysis using such primers. A set of four fluorophore-labeled primers can be reliably synthesized in a twenty-four hour period, and greater than 300 nucleotides of analyzed new sequence obtained using this set in an additional twenty-four hours. Scale-up of these procedures to take advantage of the full capabilities of the sequencer is, at present, too slow and costly to be suitable for routine sequencing, and therefore the use of specific-primers is best suited to the closure of gaps in extended sequence produced using random cloning and sequencing strategies.     

Integration of Wet and Dry Bench Processes Optimizes Targeted Next-generation Sequencing of Low-quality and Low-quantity Tumor Biopsies

All next-generation sequencing (NGS) procedures include assays performed at the laboratory bench ("wet bench") and data analyses conducted using bioinformatics pipelines ("dry bench"). Both elements are essential to produce accurate and reliable results, which are particularly critical for clinical laboratories. Targeted NGS technologies have increasingly found favor in oncology applications to help advance precision medicine objectives, yet the methods often involve disconnected and variable wet and dry bench workflows and uncoordinated reagent sets. In this report, we describe a method for sequencing challenging cancer specimens with a 21-gene panel as an example of a comprehensive targeted NGS system. The system integrates functional DNA quantification and qualification, single-tube multiplexed PCR enrichment, and library purification and normalization using analytically-verified, single-source reagents with a standalone bioinformatics suite. As a result, accurate variant calls from low-quality and low-quantity formalin-fixed, paraffin-embedded (FFPE) and fine-needle aspiration (FNA) tumor biopsies can be achieved. The method can routinely assess cancer-associated variants from an input of 400 amplifiable DNA copies, and is modular in design to accommodate new gene content. Two different types of analytically-defined controls provide quality assurance and help safeguard call accuracy with clinically-relevant samples. A flexible "tag" PCR step embeds platform-specific adaptors and index codes to allow sample barcoding and compatibility with common benchtop NGS instruments. Importantly, the protocol is streamlined and can produce 24 sequence-ready libraries in a single day. Finally, the approach links wet and dry bench processes by incorporating pre-analytical sample quality control results directly into the variant calling algorithms to improve mutation detection accuracy and differentiate false-negative and indeterminate calls. This targeted NGS method uses advances in both wetware and software to achieve high-depth, multiplexed sequencing and sensitive analysis of heterogeneous cancer samples for diagnostic applications.     

Automated sequencing of complete mitochondrial genomes from laser-capture microdissected samples

Mitochondrial DNA mutations have been related to both aging and a variety of diseases such as cancer. Due to the relatively small size of the genome (16 kb) and with the use of automated DNA sequencing, the entire genome can be sequenced from clinical specimens in days. We present a reliable approach to complete mitochondrial genome sequencing from laser-capture microdissected human clinical cancer specimens that overcome the inherent limitations of relatively small tissue samples and partial DNA degradation, which are unavoidable when laser-capture microdissection is used to attain pure populations of cells from heterogeneous tissues obtained from surgical procedures. The acquisition of sufficient template combined with a standard set of 18 pairs of PCR primers allows for the efficient amplification of the genome. Subsequent single-stranded amplification is performed using 36 sequencing primers, and samples are run on an ABI PRISM 3100 Genetic Analyzer. The use of this procedure should allow even investigators with little experience sequencing from clinical specimens success in complete mitochondrial genome sequencing.     

Design considerations for massively parallel sequencing studies of complex human disease

Massively Parallel Sequencing (MPS) allows sequencing of entire exomes and genomes to now be done at reasonable cost, and its utility for identifying genes responsible for rare Mendelian disorders has been demonstrated. However, for a complex disease, study designs need to accommodate substantial degrees of locus, allelic, and phenotypic heterogeneity, as well as complex relationships between genotype and phenotype. Such considerations include careful selection of samples for sequencing and a well-developed strategy for identifying the few "true" disease susceptibility genes from among the many irrelevant genes that will be found to harbor rare variants. To examine these issues we have performed simulation-based analyses in order to compare several strategies for MPS sequencing in complex disease. Factors examined include genetic architecture, sample size, number and relationship of individuals selected for sequencing, and a variety of filters based on variant type, multiple observations of genes and concordance of genetic variants within pedigrees. A two-stage design was assumed where genes from the MPS analysis of high-risk families are evaluated in a secondary screening phase of a larger set of probands with more modest family histories. Designs were evaluated using a cost function that assumes the cost of sequencing the whole exome is 400 times that of sequencing a single candidate gene. Results indicate that while requiring variants to be identified in multiple pedigrees and/or in multiple individuals in the same pedigree are effective strategies for reducing false positives, there is a danger of over-filtering so that most true susceptibility genes are missed. In most cases, sequencing more than two individuals per pedigree results in reduced power without any benefit in terms of reduced overall cost. Further, our results suggest that although no single strategy is optimal, simulations can provide important guidelines for study design.     

Parallel tagged sequencing on the 454 platform

Parallel tagged sequencing (PTS) is a molecular barcoding method designed to adapt the recently developed high-throughput 454 parallel sequencing technology for use with multiple samples. Unlike other barcoding methods, PTS can be applied to any type of double-stranded DNA (dsDNA) sample, including shotgun DNA libraries and pools of PCR products, and requires no amplification or gel purification steps. The method relies on attaching sample-specific barcoding adapters, which include sequence tags and a restriction site, to blunt-end repaired DNA samples by ligation and strand-displacement. After pooling multiple barcoded samples, molecules without sequence tags are effectively excluded from sequencing by dephosphorylation and restriction digestion, and using the tag sequences, the source of each DNA sequence can be traced. This protocol allows for sequencing 300 or more complete mitochondrial genomes on a single 454 GS FLX run, or twenty-five 6-kb plasmid sequences on only one 16th plate region. Most of the reactions can be performed in a multichannel setup on 96-well reaction plates, allowing for processing up to several hundreds of samples in a few days.     

A preprocessor for shotgun assembly of large genomes.

The whole-genome shotgun (WGS) assembly technique has been remarkably successful in efforts to determine the sequence of bases that make up a genome. WGS assembly begins with a large collection of short fragments that have been selected at random from a genome. The sequence of bases at each end of the fragment is determined, albeit imprecisely, resulting in a sequence of letters called a "read." Each letter in a read is assigned a quality value, which estimates the probability that a sequencing error occurred in determining that letter. Reads are typically cut off after about 500 letters, where sequencing errors become endemic. We report on a set of procedures that (1) corrects most of the sequencing errors, (2) changes quality values accordingly, and (3) produces a list of "overlaps," i.e., pairs of reads that plausibly come from overlapping parts of the genome. Our procedures, which we call collectively the "UMD Overlapper," can be run iteratively and as a preprocessor for other assemblers. We tested the UMD Overlapper on Celera's Drosophila reads. When we replaced Celera's overlap procedures in the front end of their assembler, it was able to produce a significantly improved genome.     

Extending assembly of short DNA sequences to handle error

Inexpensive de novo genome sequencing, particularly in organisms with small genomes, is now possible using several new sequencing technologies. Some of these technologies such as that from Illumina's Solexa Sequencing, produce high genomic coverage by generating a very large number of small reads ( approximately 30 bp). While prior work shows that partial assembly can be performed by k-mer extension in error-free reads, this algorithm is unsuccessful with the sequencing error rates found in practice. We present VCAKE (Verified Consensus Assembly by K-mer Extension), a modification of simple k-mer extension that overcomes error by using high depth coverage. Though it is a simple modification of a previous approach, we show significant improvements in assembly results on simulated and experimental datasets that include error. AVAILABILITY: http://152.2.15.114/~labweb/VCAKE     

The construction and analysis of M13 libraries prepared from YAC DNA.

Yeast artificial chromosomes (YACs) provide a powerful way to isolate and map large regions of genomic DNA and their use in genome analysis is now extensive. We modified a series of procedures to produce high quality shotgun libraries from small amounts of YAC DNA. Clones from several different libraries have been sequenced and analyzed for distribution, sequence integrity and degree of contamination from yeast DNA. We describe these procedures and analyses and show that sequencing at about 1-fold coverage, followed by database comparison (survey sequencing) offers a relatively quick method to determine the nature of previously uncharacterized cosmid or YAC clones.     

Challenges in Laboratory Detection of Fungal Pathogens in the Airways of Cystic Fibrosis Patients

Study of the clinical significance of fungal colonization/infection in the airways of cystic fibrosis (CF) patients, especially by filamentous fungi, is challenged by the absence of standardized methodology for the detection and identification of an ever-broadening range of fungal pathogens. Culture-based methods remain the cornerstone diagnostic approaches, but current methods used in many clinical laboratories are insensitive and unstandardized, rendering comparative studies unfeasible. Guidelines for standardized processing of respiratory specimens and for their culture are urgently needed and should include recommendations for specific processing procedures, inoculum density, culture media, incubation temperature and duration of culture. Molecular techniques to detect fungi directly from clinical specimens include panfungal PCR assays, multiplex or pathogen-directed assays, real-time PCR, isothermal methods and probe-based assays. In general, these are used to complement culture. Fungal identification by DNA sequencing methods is often required to identify cultured isolates, but matrix-assisted laser desorption/ionization time-of-flight mass spectrometry is increasingly used as an alternative to DNA sequencing. Genotyping of isolates is undertaken to investigate relatedness between isolates, to pinpoint the infection source and to study the population structure. Methods range from PCR fingerprinting and amplified fragment length polymorphism analysis, to short tandem repeat typing, multilocus sequencing typing (MLST) and whole genome sequencing (WGS). MLST is the current preferred method, whilst WGS offers best case resolution but currently is understudied.     

Rapid and reliable fluorescent cycle sequencing of double-stranded templates.

Automated DNA sequencing is an extremely valuable technique which requires very high quality DNA templates to be carried out successfully. While it has been possible to readily produce large numbers of such templates from M13 or other single-stranded vectors for several years, the sequencing of double-stranded DNA templates using the ABI 373 DNA Sequencer has had a considerably lower success rate. We describe how the combination of a new fluorescent, dideoxy sequencing method, called cycle-sequencing, coupled with modifications to template isolation procedures based on Qiagen columns, makes fluorescent sequencing of double-stranded templates a reliable procedure. From a single five milliliter culture enough DNA can be isolated (up to 20 micrograms) to do 4-8 sequencing reactions, each of which yields 400-500 bases of high quality sequence data. These procedures make the routine use of double-stranded DNA templates a viable strategy in automated DNA sequencing projects.     

The impact of post-alignment processing procedures on whole-exome sequencing data

The use of post-alignment procedures has been suggested to prevent the identification of false-positives in massive DNA sequencing data. Insertions and deletions are most likely to be misinterpreted by variant calling algorithms. Using known genetic variants as references for post-processing pipelines can minimize mismatches. They allow reads to be correctly realigned and recalibrated, resulting in more parsimonious variant calling. In this work, we aim to investigate the impact of using different sets of common variants as references to facilitate variant calling from whole-exome sequencing data. We selected reference variants from common insertions and deletions available within the 1K Genomes project data and from databases from the Latin American Database of Genetic Variation (LatinGen). We used the Genome Analysis Toolkit to perform post-processing procedures like local realignment, quality recalibration procedures, and variant calling in whole exome samples. We identified an increased number of variants from the call set for all groups when no post-processing procedure was performed. We found that there was a higher concordance rate between variants called using 1K Genomes and LatinGen. Therefore, we believe that the increased number of rare variants identified in the analysis without realignment or quality recalibration indicated that they were likely false-positives.     

Next-generation Sequencing of 16S Ribosomal RNA Gene Amplicons

One of the major questions in microbial ecology is “who is there?” This question can be answered using various tools, but one of the long-lasting gold standards is to sequence 16S ribosomal RNA (rRNA) gene amplicons generated by domain-level PCR reactions amplifying from genomic DNA. Traditionally, this was performed by cloning and Sanger (capillary electrophoresis) sequencing of PCR amplicons. The advent of next-generation sequencing has tremendously simplified and increased the sequencing depth for 16S rRNA gene sequencing. The introduction of benchtop sequencers now allows small labs to perform their 16S rRNA sequencing in-house in a matter of days. Here, an approach for 16S rRNA gene amplicon sequencing using a benchtop next-generation sequencer is detailed. The environmental DNA is first amplified by PCR using primers that contain sequencing adapters and barcodes. They are then coupled to spherical particles via emulsion PCR. The particles are loaded on a disposable chip and the chip is inserted in the sequencing machine after which the sequencing is performed. The sequences are retrieved in fastq format, filtered and the barcodes are used to establish the sample membership of the reads. The filtered and binned reads are then further analyzed using publically available tools. An example analysis where the reads were classified with a taxonomy-finding algorithm within the software package Mothur is given. The method outlined here is simple, inexpensive and straightforward and should help smaller labs to take advantage from the ongoing genomic revolution.     

Arthroscopy.

Arthroscopy has reduced the morbidity and period of hospitalisation associated with orthopaedic surgery and has increased the range of procedures that may be performed. From early operations on the knee it has expanded to include procedures for the shoulder, elbow, wrist, hip, ankle, and foot. For some joints the indications for surgery are clear, for others the clinical advantages are still being assessed. This expansion has also led to the recognition of complications, though the incidence is low. Specialist instrumentation has allowed a wide variety of operations previously needing open surgery to be carried out arthroscopically. The repertoire of arthroscopic procedures will undoubtedly continue to expand, and controlled studies are required to validate their efficacy, particularly in the management of degenerative joint diseases.     

A method for determining DNA sequence by labeling the end of the molecule and cleaving at the base. Isolation of DNA fragments, end-labeling, cleavage, electrophoresis in polyacrylamide gel and analysis of results

We elaborate basic chemical principles and current laboratory procedures for sequencing end-labeled DNA by partial cleavage and gel electrophoresis (A. M. Maxam and W. Gilbert, Proc. Natl. Acad. Sci. USA, 1977, v. 74, p. 560-564). We provide step-by-step protocols for 32P-labeling DNA ends, segregating the labeled ends by cutting with a second restriction enzyme or separating strands, partially cleaving the DNA at specific bases with reagents, electrophoresing the labeled products of cleavage on sequencing gels, and interpreting sequencing band patterns. Many of these procedures have been condensed, to make them faster and easier, and some are new. We also discuss sequencing strategies, and suggest a technique which will reduce plasmid or viral DNA to a collection of singly-end-labeled fragments in one day, for efficient sequencing of these chromosomes in 250-nucleotide blocks.     

From micrograms to picograms: quantitative PCR reduces the material demands of high-throughput sequencing

Current efforts to recover the Neandertal and mammoth genomes by 454 DNA sequencing demonstrate the sensitivity of this technology. However, routine 454 sequencing applications still require microgram quantities of initial material. This is due to a lack of effective methods for quantifying 454 sequencing libraries, necessitating expensive and labour-intensive procedures when sequencing ancient DNA and other poor DNA samples. Here we report a 454 sequencing library quantification method based on quantitative PCR that effectively eliminates these limitations. We estimated both the molecule numbers and the fragment size distributions in sequencing libraries derived from Neandertal DNA extracts, SAGE ditags and bonobo genomic DNA, obtaining optimal sequencing yields without performing any titration runs. Using this method, 454 sequencing can routinely be performed from as little as 50 pg of initial material without titration runs, thereby drastically reducing costs while increasing the scope of sample throughput and protocol development on the 454 platform. The method should also apply to Illumina/Solexa and ABI/SOLiD sequencing, and should therefore help to widen the accessibility of all three platforms.     

Fluorescent and radioactive solid phase dideoxy sequencing of pcr products in microtitre plates.

In this paper we describe a rapid method for the direct generation of DNA sequencing templates from phage or bacteria. Sequencing of these PCR products can be performed by radioactive and fluorescent methods. The non-radioactive method has been used to sequence a total of approximately 100 kb of human DNA fragments generated by digestion with HpaII and subsequent cloning. The method depends on direct small scale amplification using a biotinylated primer, and the binding of the product to streptavidin coated magnetic beads. All the procedures are carried out in a microtitre plate thus facilitating the handling of large numbers of clones and has potential for automation.     

Primary structures of Escherichia coli pyruvate formate-lyase and pyruvate-formate-lyase-activating enzyme deduced from the DNA nucleotide sequences

The structural gene of pyruvate formate-lyase (pfl) and that of pyruvate-formate-lyase-activating enzyme were shown to be adjacent on the chromosomal map of Escherichia coli. DNA sequencing was performed along a stretch of 3592 nucleotides to obtain the amino acid sequences of both proteins. The derived primary structures (759 and 245 residues) were confirmed by partial structure analyses on the purified proteins. The open reading frames are separated by a 194-nucleotide stretch, and their flanking regions include signal elements that are compatible with separate control of protein synthesis from the two genes.     

Quantification of mitochondrial DNA mutation load

Mitochondrial DNA (mtDNA) mutations are an important cause of genetic disease and have been proposed to play a role in the ageing process. Quantification of total mtDNA mutation load in ageing tissues is difficult as mutational events are rare in a background of wild-type molecules, and detection of individual mutated molecules is beyond the sensitivity of most sequencing based techniques. The methods currently most commonly used to document the incidence of mtDNA point mutations in ageing include post-PCR cloning, single-molecule PCR and the random mutation capture assay. The mtDNA mutation load obtained by these different techniques varies by orders of magnitude, but direct comparison of the three techniques on the same ageing human tissue has not been performed. We assess the procedures and practicalities involved in each of these three assays and discuss the results obtained by investigation of mutation loads in colonic mucosal biopsies from ten human subjects.     

Solid-phase protein and peptide sequencing using either 4-N,N-dimethylaminoazobenzene 4'-isothiocyanate or phenylisothiocyanate

Automated solid-phase sequencing using 4-N,N-dimethylaminoazobenzene 4'-isothiocyanate (DABITC) double coupling or regular phenylisothiocyanate (PITC) degradation procedures have been investigated. Employing sensitive high-performance liquid chromatography for the identification of amino acid thiohydantoin derivatives (PTH and DABTH), both methods were capable of sequencing immobilized peptides or proteins at the subnanomole levels. In the sequencing program using DABITC, alternate methanol and dichloroethane washes and automated conversion using methanolic HCl containing dithiothreitol were introduced to obtain clean thiazolinones and to ensure high recovery yields of the thiohydantoins. Using regular PITC degradation with a 59-min program, the background peaks of the side products could be reduced to enhance HPLC identification. Peptides or proteins attached to the glass beads or resins via the carboxyl terminii or epsilon-amino groups of lysyl residues could be readily sequenced up to 30 identifiable degradation cycles, where the sequencing is generally terminated due to the increased background components.