Automation of Molecular-Based Analyses: A Primer on Massively Parallel Sequencing

Recent advances in genetics have been enabled by new genetic sequencing techniques called massively parallel sequencing (MPS) or next-generation sequencing. Through the ability to sequence in parallel hundreds of thousands to millions of DNA fragments, the cost and time required for sequencing has dramatically decreased. There are a number of different MPS platforms currently available and being used in Australia. Although they differ in the underlying technology involved, their overall processes are very similar: DNA fragmentation, adaptor ligation, immobilisation, amplification, sequencing reaction and data analysis. MPS is being used in research, translational and increasingly now also in clinical settings. Common applications include sequencing of whole genomes, whole exomes or targeted genes for disease-causing gene discovery, genetic diagnosis and targeted cancer therapy. Even though the revolution that is occurring with MPS is exciting due to its increasing use, improving and emerging technologies and new applications, significant challenges still exist. Particularly challenging issues are the bioinformatics required for data analysis, interpretation of results and the ethical dilemma of ‘incidental findings’.     

DNA sequencing by capillary array electrophoresis with an electric field strength gradient

We examined the effect of an electric field strength gradient on DNA sequencing efficiency using capillary array electrophoresis. Several types of gradients were applied to DNA sequencing and tested in terms of read length and accuracy. Our original method improved the accuracy of DNA sequencing for longer fragments at high temperature.     

Tailed pooled suppression subtractive hybridization (PSSH) adaptors do not alter efficiency

Suppression Subtractive Hybridization (SSH) and its derivative, Pooled Suppression Subtractive hybridization (PSSH), are powerful tools used to study variances larger than ~100 bp in prokaryotic genome structure. The initial steps involve ligating an oligonucleotide of known sequence (the “adaptor”) to a fragmented genome to facilitate amplification, subtraction and downstream sequencing. SSH results in the creation of a library of unique DNA fragments which have been traditionally analyzed via Sanger sequencing. Numerous next generation sequencing technologies have entered the market yet SSH is incompatible with these platforms. This is due to the high level of sequence conservation of the oligonucleotide used for SSH. This rigid adherence is partly because it has yet to be determined if alteration of this oligonucleotide will have a deleterious impact on subtraction efficiency. The subtraction occurs when non-unique fragments are inhibited by a secondary self-pairing structure which requires exact nucleotide sequence. We determine if appending custom sequence to the 5′ terminal ends of these oligonucleotides during the nested PCR stages of PSSH will reduce subtraction efficiency. We compare a pool of ten S. aureus clinical isolates with a standard PSSH and custom tailed-PSSH. We detected no statistically significant difference between their subtraction efficiencies. Our observations suggest that the adaptor’s terminal ends may be labeled during the nested PCR step. This produces libraries labeled with custom sequence. This does not lead to loss of subtraction efficiency and would be invaluable for groups wishing to combine SSH or PSSH with their own downstream applications, such as a high throughput sequencing platform.     

Next generation sequencing and bioinformatic bottlenecks: the current state of metagenomic data analysis

The recent technological advances in next generation sequencing have brought the field closer to the goal of reconstructing all genomes within a community by presenting high throughput sequencing at much lower costs. While these next-generation sequencing technologies have allowed a massive increase in available raw sequence data, there are a number of new informatics challenges and difficulties that must be addressed to improve the current state, and fulfill the promise of, metagenomics.     

Sequence analysis of genes and genomes

A major step towards understanding of the genetic basis of an organism is the complete sequence determination of all genes in its genome. The development of powerful techniques for DNA sequencing has enabled sequencing of large amounts of gene fragments and even complete genomes. Important new techniques for physical mapping, DNA sequencing and sequence analysis have been developed. To increase the throughput, automated procedures for sample preparation and new software for sequence analysis have been applied. This review describes the development of new sequencing methods and the optimisation of sequencing strategies for whole genome and cDNA analysis, as well as discusses issues regarding sequence analysis and annotation.     

The nucleotide and partial amino acid sequence of toxic shock syndrome toxin-1

The nucleotide sequence of toxic shock syndrome toxin-1 (TSST-1) has been determined. In addition, one-third of the predicted amino acid sequence was confirmed by amino acid sequence analysis of cyanogen bromide-generated TSST-1 protein fragments. The DNA sequencing results identified a 708-base pair open reading frame starting with an ATG, 7 base pairs downstream from a Shine-Dalgarno sequence, and terminating at a UAA stop codon. Amino acid analysis of the intact protein defined the NH2 terminus of the mature protein and located the cleavage point for the signal peptide (Ala/Ser). The signal peptide contained the first 40 amino acids and had characteristic structural similarities with other bacterial signal peptides. The coding sequence of the mature protein was 585 base pairs (194 amino acids) in length, and the molecular weight of the predicted protein was 22,049. This is in good agreement with the previously reported molecular weight of TSST-1 (22,000), as determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. NH2-terminal amino acid sequence analysis performed on isolated TSST-1 CNBr fragments determined the position of the peptides in the TSST-1 sequence and verified the predicted amino acid sequence in those positions. Computer analyses of the amino acid sequence showed that TSST-1 has little or no sequence homology with biologically related toxins, streptococcal pyrogenic exotoxin A, and staphylococcal enterotoxins B and C.     

Use of competitive DNA hybridization to identify differences in the genomes of bacteria

Although recent technological advances in DNA sequencing and computational biology now allow scientists to compare entire microbial genomes, comparisons of closely related bacterial species and individual isolates by whole-genome sequencing approaches remains prohibitively expensive for most laboratories. Here we report the development and testing of a biochemical approach for targeted sequencing of only those chromosomal regions that differ between two DNA preparations. The method, designated GFE (genome fragment enrichment) uses competitive solution hybridization and positive selection to obtain genomic DNA fragments that are present in one pool of fragments but not another. Repeated comparisons of the genomes of Enterococcus faecalis and E. faecium led to the identification of 225 putative genome-specific DNA fragments. Species and strain variations within these fragments were confirmed by both experimental and bioinformatic analyses. The E. faecalis genome-specific sequences identified included both a preponderance of those predicted to encode surface-exposed proteins, as well as several previously described unique marker regions embedded within highly conserved rrn operons. The GFE strategy we describe efficiently identified genomic differences between two enterococcal genomes, and will be widely applicable for studying genetic variation among closely related bacterial species.     

Analysis of DNA sequences using a single chemical cleavage procedure

A novel approach to sequence analysis of end-labeled, defined DNA fragments, using a single chemical cleavage procedure and electrophoretic separation in a single lane, has been developed. Prolonged treatment with hot aqueous piperidine results in partial cleavage of the DNA at all positions; the relative propensity for this cleavage is different for the various bases in the DNA. The hydrolysate is resolved on a DNA sequencing gel, and the distribution of radioactivity in the electrophoretic lane is analyzed (a) in terms of differential peak heights of the radioactive bands and (b) in terms of the spacings between successive bands. Simultaneous application of these two base-characteristic criteria allows the deduction of the nucleotide sequence with an accuracy approaching that of the established four-lane methods of DNA sequencing.     

DNA sequencing separations in capillary gels on a modified commercial DNA sequencing instrument

DNA sequencing separations of standard DNA fragments of known sequence have been achieved in small diameter capillary gels electrophoresed and analyzed in parallel in a modified commercial DNA sequencer instrument. DNA sequencing in terms of base-calling accuracy is comparable to conventional slab gels; however, the separations in the capillary were performed somewhat faster and required less sample than those in the slab gel. Advantages of this approach vs. separations on conventional slab gels are discussed.     

Improving ancient DNA read mapping against modern reference genomes

ABSTRACT: BACKGROUND: Next-Generation Sequencing has revolutionized our approach to ancient DNA (aDNA) research, by providing complete genomic sequences of ancient individuals and extinct species. However, the recovery of genetic material from long-dead organisms is still complicated by a number of issues, including post-mortem DNA damage and high levels of environmental contamination. Together with error profiles specific to the type of sequencing platforms used, these specificities could limit our ability to map sequencing reads against modern reference genomes and therefore limit our ability to identify endogenous ancient reads, reducing the efficiency of shotgun sequencing aDNA. RESULTS: In this study, we compare different computational methods for improving the accuracy and sensitivity of aDNA sequence identification, based on shotgun sequencing reads recovered from Pleistocene horse extracts using Illumina GAIIx and Helicos Heliscope platforms. We show that the performance of the Burrows Wheeler Aligner (BWA), that has been developed for mapping of undamaged sequencing reads using platforms with low rates of indel-types of sequencing errors, can be employed at acceptable run-times by modifying default parameters in a platform-specific manner. We also examine if trimming likely damaged positions at read ends can increase the recovery of genuine aDNA fragments and if accurate identification of human contamination can be achieved using a strategy previously suggested based on best hit filtering. We show that combining our different mapping and filtering approaches can increase the number of high-quality endogenous hits recovered by up to 33%. CONCLUSIONS: We have shown that Illumina and Helicos sequences recovered from aDNA extracts could not be aligned to modern reference genomes with the same efficiency unless mapping parameters are optimized for the specific types of errors generated by these platforms and by post-mortem DNA damage. Our findings have important implications for future aDNA research, as we define mapping guidelines that improve our ability to identify genuine aDNA sequences, which in turn could improve the genotyping accuracy of ancient specimens. Our framework provides a significant improvement to the standard procedures used for characterizing ancient genomes, which is challenged by contamination and often low amounts of DNA material.     

Resolving Prokaryotic Taxonomy without rRNA: Longer Oligonucleotide Word Lengths Improve Genome and Metagenome Taxonomic Classification

Oligonucleotide signatures, especially tetranucleotide signatures, have been used as method for homology binning by exploiting an organism’s inherent biases towards the use of specific oligonucleotide words. Tetranucleotide signatures have been especially useful in environmental metagenomics samples as many of these samples contain organisms from poorly classified phyla which cannot be easily identified using traditional homology methods, including NCBI BLAST. This study examines oligonucleotide signatures across 1,424 completed genomes from across the tree of life, substantially expanding upon previous work. A comprehensive analysis of mononucleotide through nonanucleotide word lengths suggests that longer word lengths substantially improve the classification of DNA fragments across a range of sizes of relevance to high throughput sequencing. We find that, at present, heptanucleotide signatures represent an optimal balance between prediction accuracy and computational time for resolving taxonomy using both genomic and metagenomic fragments. We directly compare the ability of tetranucleotide and heptanucleotide world lengths (tetranucleotide signatures are the current standard for oligonucleotide word usage analyses) for taxonomic binning of metagenome reads. We present evidence that heptanucleotide word lengths consistently provide more taxonomic resolving power, particularly in distinguishing between closely related organisms that are often present in metagenomic samples. This implies that longer oligonucleotide word lengths should replace tetranucleotide signatures for most analyses. Finally, we show that the application of longer word lengths to metagenomic datasets leads to more accurate taxonomic binning of DNA scaffolds and have the potential to substantially improve taxonomic assignment and assembly of metagenomic data.     

Recent patents of nanopore DNA sequencing technology: progress and challenges

DNA sequencing techniques witnessed fast development in the last decades, primarily driven by the Human Genome Project. Among the proposed new techniques, Nanopore was considered as a suitable candidate for the single DNA sequencing with ultrahigh speed and very low cost. Several fabrication and modification techniques have been developed to produce robust and well-defined nanopore devices. Many efforts have also been done to apply nanopore to analyze the properties of DNA molecules. By comparing with traditional sequencing techniques, nanopore has demonstrated its distinctive superiorities in main practical issues, such as sample preparation, sequencing speed, cost-effective and read-length. Although challenges still remain, recent researches in improving the capabilities of nanopore have shed a light to achieve its ultimate goal: Sequence individual DNA strand at single nucleotide level. This patent review briefly highlights recent developments and technological achievements for DNA analysis and sequencing at single molecule level, focusing on nanopore based methods.     

Analysis and purification of nucleic acids by ion-pair reversed-phase high-performance liquid chromatography

Sizing of DNA fragments is a routine analysis traditionally performed on agarose or polyacrylamide gels. Electrophoretic analysis is labor-intensive with only limited potential for automation. Recovery of DNA fragments from gels is cumbersome. We present data on automated, size-based separation of DNA fragments by ion-pair reversed-phase high performance liquid chromatography (IP RP HPLC) - DNA chromatography - on the WAVE DNA Fragment Analysis System with the DNASep cartridge. This system is suitable for accurate and rapid sizing of double-stranded (ds) DNA fragments from 50 to ca. 2000 base pairs (bp). Fluorescently labeled DNA fragments are compatible with the technology. Length-dependent separation of dsDNA fragments is sequence independent and retention times are highly reproducible. The resolving capabilities of DNA chromatography are illustrated by the analysis of multiple DNA size markers. Resolved dsDNA fragments are easily collected and are suitable for downstream applications such as sequencing and cloning. DNA chromatography under denaturing conditions with fluorescently labeled DNA fragments offers a means for the separation and purification of individual strands of dsDNA. Analysis of DNA fragments on the WAVE System is highly automated and requires minimal manual intervention. DNA chromatography offers a reliable and automated alternative to gel electrophoresis for the analysis of DNA fragments.     

The influence of fluorescent dye structure on the electrophoretic mobility of end-labeled DNA.

Over the past 10 years, fluorescent end-labeling of DNA fragments has evolved into the preferred method of DNA detection for a wide variety of applications, including DNA sequencing and PCR fragment analysis. One of the advantages inherent in fluorescent detection methods is the ability to perform multi-color analyses. Unfortunately, labeling DNA fragments with different fluorescent tags generally induces disparate relative electrophoretic mobilities for the fragments. Mobility-shift corrections must therefore be applied to the electrophoretic data to compensate for these effects. These corrections may lead to increased errors in the estimation of DNA fragment sizes and reduced confidence in DNA sequence information. Here, we present a systematic study of the relationship between dye structure and the resultant electrophoretic mobility of end-labeled DNA fragments. We have used a cyanine dye family as a paradigm and high-resolution capillary array electrophoresis (CAE) as the instrumentation platform. Our goals are to develop a general understanding of the effects of dyes on DNA electrophoretic mobility and to synthesize a family of DNA end-labels that impart identically matched mobility influences on DNA fragments. Such matched sets could be used in DNA sequencing and fragment sizing applications on capillary electrophoresis instrumentation.     

Mass-spectrometry DNA sequencing

Matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) has been explored widely for DNA sequencing. Compared to gel electrophoresis based sequencing systems, mass spectrometry produces very high resolution of sequencing fragments, rapid separation on microsecond time scales, and completely eliminates compressions associated with gel-based systems. While most of the research efforts have focused on using mass spectrometers to analyze the DNA products from Sanger sequencing or enzymatic digestion reactions, the read lengths attainable are currently insufficient for large-scale de novo sequencing. The advantage of mass-spectrometry sequencing is that one can unambiguously identify frameshift mutations and heterozygous mutations making it an ideal choice for resequencing projects. In these applications, DNA sequencing fragments that are the same length but with different base compositions are generated, which are challenging to consistently distinguish in gel-based sequencing systems. In contrast, MALDI-TOF MS produces mass spectra of these DNA sequencing fragments with nearly digital resolution, allowing accurate determination of the mixed bases. For these reasons mass spectrometry based sequencing has mainly been focused on the detection of frameshift mutations and single nucleotide polymorphisms (SNPs). More recently, assays have been developed to indirectly sequence DNA by first converting it into RNA. These assays take advantage of the increased resolution and detection ability of MALDI-TOF MS for RNA.     

Local de novo assembly of RAD paired-end contigs using short sequencing reads

Despite the power of massively parallel sequencing platforms, a drawback is the short length of the sequence reads produced. We demonstrate that short reads can be locally assembled into longer contigs using paired-end sequencing of restriction-site associated DNA (RAD-PE) fragments. We use this RAD-PE contig approach to identify single nucleotide polymorphisms (SNPs) and determine haplotype structure in threespine stickleback and to sequence E. coli and stickleback genomic DNA with overlapping contigs of several hundred nucleotides. We also demonstrate that adding a circularization step allows the local assembly of contigs up to 5 kilobases (kb) in length. The ease of assembly and accuracy of the individual contigs produced from each RAD site sequence suggests RAD-PE sequencing is a useful way to convert genome-wide short reads into individually-assembled sequences hundreds or thousands of nucleotides long.     

Resurrecting ancient animal genomes: the extinct moa and more

Recently two developments have had a major impact on the field of ancient DNA (aDNA). First, new advances in DNA sequencing, in combination with improved capture/enrichment methods, have resulted in the recovery of orders of magnitude more DNA sequence data from ancient animals. Second, there has been an increase in the range of tissue types employed in aDNA. Hair in particular has proven to be very successful as a source of DNA because of its low levels of contamination and high level of ancient endogenous DNA. These developments have resulted in significant advances in our understanding of recently extinct animals: namely their evolutionary relationships, physiology, and even behaviour. Hair has been used to recover the first complete ancient nuclear genome, that of the extinct woolly mammoth, which then facilitated the expression and functional analysis of haemoglobins. Finally, we speculate on the consequences of these developments for the possibility of recreating extinct animals.     

Methodological improvements of pyrosequencing technology

Pyrosequencing technology is a rather novel DNA sequencing method based on the sequencing-by-synthesis principle. This bioluminometric, real-time DNA sequencing technique employs a cascade of four enzymatic reactions producing sequence peak signals. The method has been proven highly suitable for single nucleotide polymorphism analysis and sequencing of short stretches of DNA. Although the pyrosequencing procedure is relatively straightforward, users may face challenges due to varying parameters in PCR and sequencing primer design, sample preparation and nucleotide dispensation; such challenges are labor and cost intensive. In this study, these issues have been addressed to increase signal quality and assure sequence accuracy.