Next-generation sequencing: from basic research to diagnostics.

BACKGROUND: For the past 30 years, the Sanger method has been the dominant approach and gold standard for DNA sequencing. The commercial launch of the first massively parallel pyrosequencing platform in 2005 ushered in the new era of high-throughput genomic analysis now referred to as next-generation sequencing (NGS). CONTENT: This review describes fundamental principles of commercially available NGS platforms. Although the platforms differ in their engineering configurations and sequencing chemistries, they share a technical paradigm in that sequencing of spatially separated, clonally amplified DNA templates or single DNA molecules is performed in a flow cell in a massively parallel manner. Through iterative cycles of polymerase-mediated nucleotide extensions or, in one approach, through successive oligonucleotide ligations, sequence outputs in the range of hundreds of megabases to gigabases are now obtained routinely. Highlighted in this review are the impact of NGS on basic research, bioinformatics considerations, and translation of this technology into clinical diagnostics. Also presented is a view into future technologies, including real-time single-molecule DNA sequencing and nanopore-based sequencing. SUMMARY: In the relatively short time frame since 2005, NGS has fundamentally altered genomics research and allowed investigators to conduct experiments that were previously not technically feasible or affordable. The various technologies that constitute this new paradigm continue to evolve, and further improvements in technology robustness and process streamlining will pave the path for translation into clinical diagnostics.

Targeted gene panel sequencing for the rapid diagnosis of acutely ill infants.

BACKGROUND: Exome/genome sequencing (ES/GS) have been recently used in neonatal and pediatric/cardiac intensive care units (NICU and PICU/CICU) to diagnose and care for acutely ill infants, but the effectiveness of targeted gene panels for these purposes remains unknown. METHODS: RapSeq, a newly developed panel targeting 4,503 disease-causing genes, was employed on selected patients in our NICU/PICU/CICU. Twenty trios were sequenced from October 2015 to March 2017. We assessed diagnostic yield, turnaround times, and clinical consequences. RESULTS: A diagnosis was made in 10/20 neonates (50%); eight had de novo variants (ASXL1, CHD, FBN1, KMT2D, FANCB, FLNA, PAX3), one was a compound heterozygote for CHAT, and one had a maternally inherited GNAS variant. Preliminary reports were generated by 9.6 days (mean); final reports after Sanger sequencing at 16.3 days (mean). In all positive infants, the diagnosis changed management. In a case with congenital myasthenia, diagnosis and treatment occurred at 17 days versus 7 months in a historical control. CONCLUSIONS: This study shows that a gene panel that includes the majority of known disease-causing genes can rapidly identify a diagnosis in a large number of tested infants. Due to simpler deployment and interpretation and lower costs, this approach might represent an alternative to ES/GS in the NICU/PICU/CICU.

Specific versus nonspecific isothermal DNA amplification through thermophilic polymerase and nicking enzyme activities.

Rapid isothermal nucleic acid amplification technologies can enable diagnosis of human pathogens and genetic variations in a simple, inexpensive, user-friendly format. The isothermal exponential amplification reaction (EXPAR) efficiently amplifies short oligonucleotides called triggers in less than 10 min by means of thermostable polymerase and nicking endonuclease activities. We recently demonstrated that this reaction can be coupled with upstream generation of trigger oligonucleotides from a genomic target sequence, and with downstream visual detection using DNA-functionalized gold nanospheres. The utility of EXPAR in clinical diagnostics is, however, limited by a nonspecific background amplification phenomenon, which is further investigated in this report. We found that nonspecific background amplification includes an early phase and a late phase. Observations related to late phase background amplification are in general agreement with literature reports of ab initio DNA synthesis. Early phase background amplification, which limits the sensitivity of EXPAR, differs however from previous reports of nonspecific DNA synthesis. It is observable in the presence of single-stranded oligonucleotides following the EXPAR template design rules and generates the trigger sequence expected for the EXPAR template present in the reaction. It appears to require interaction between the DNA polymerase and the single-stranded EXPAR template. Early phase background amplification can be suppressed or eliminated by physically separating the template and polymerase until the final reaction temperature has been reached, thereby enabling detection of attomolar starting trigger concentrations.

Characterization of aberrant melting peaks in unlabeled probe assays.

An unlabeled probe assay relies on a double-stranded DNA-binding dye to detect and verify target based on amplicon and probe melting. During the development and application of unlabeled probe assays, aberrant melting peaks are sometimes observed that may interfere with assay interpretation. In this report, we investigated the origin of aberrant melting profiles observed in an unlabeled probe assay for exon 10 of the RET gene. It was determined that incomplete 3' blocking of the unlabeled probe allowed polymerase-mediated probe extension resulting in extension products that generated the aberrant melting profiles. This report further examined the blocking ability of the 3' modifications C3 spacer, amino-modified C6, phosphate, inverted dT, and single 3' nucleotide mismatches in unlabeled probe experiments. Although no 3' blocking modifications in these experiments were 100% effective, the amino-modified C6, inverted dT, and C3 spacer provided the best blocking efficiencies (1% or less unblocked), phosphate was not as effective of a block (up to 2% unblocked), and single nucleotide mismatches should be avoided as a 3' blocking modification.

Closed-tube SNP genotyping without labeled probes/a comparison between unlabeled probe and amplicon melting.

Two methods for closed-tube single nucleotide polymorphism (SNP) genotyping without labeled probes have become available: unlabeled probe and amplicon melting. Unlabeled probe and amplicon melting assays were compared using 5 SNPs: human platelet antigens 1, 2, 5, and 15 and a C>T variant located 13910 base pairs (bp) upstream of the lactase gene. LCGreen Plus (Idaho Technology, Salt Lake City, UT) was used as the saturating DNA dye. Unlabeled probe data were readily interpretable and accurate for all amplicon lengths tested. Five targets that ranged in size from 42 to 72 bp were well resolved by amplicon melting on the LightScanner (Idaho Technology) or LightTyper (Roche, Indianapolis, IN) with no errors in genotyping. However, when larger amplicons (206 bp) were used and analyzed on lower resolution instruments (LightTyper and I-Cycler, Bio-Rad, Hercules, CA), the accuracy of amplicon genotyping was only 73% to 77%. When 2 temperature standards were used to bracket the amplicon of interest, the accuracy of amplicon genotyping of SNPs was increased to 100% (LightTyper) and 88% (I-Cycler).

A single-tube nucleic acid extraction, amplification, and detection method using aluminum oxide.

A disposable 0.2-ml polymerase chain reaction (PCR) tube modified with an aluminum oxide membrane (AOM) has been developed for the extraction, amplification, and detection of nucleic acids. To assess the dynamic range of AOM tubes for real-time PCR, quantified herpes simplex virus (HSV) DNA was used to compare AOM tubes to standard PCR tubes. AOM PCR tubes used for amplification and detection of quantified HSV-1 displayed a crossing threshold (C(T)) shift 0.1 cycles greater than PCR tube controls. Experiments with HSV-1-positive cerebrospinal fluid (CSF) examined the extraction, amplification, and detection properties of the AOM tubes compared to the Qiagen DNA blood mini kit. The AOM extraction, amplification, and detection of HSV-1 in CSF displayed differences of less than one C(T) when compared to Qiagen-extracted samples. Experiments testing the AOM method using clinical CSF samples displayed 100% concordance with reported results. AOM tubes have no adverse effects on amplification or fluorescence acquisition by real-time PCR and can be effectively used for the extraction, amplification, and detection of HSV from CSF. The AOM single tube method is a fast, reliable, and reproducible technique for the extraction, amplification, and detection of HSV in CSF.

A high-throughput next-generation sequencing assay for the mitochondrial genome.

Next-generation sequencing (NGS) is an effective method for mitochondrial genome (mtDNA) sequencing and heteroplasmy detection. The following protocol describes an mtDNA enrichment method up to library preparation and sequencing on Illumina NGS platforms. A short command line alignment script is available for download via FTP.

The development of next-generation sequencing assays for the mitochondrial genome and 108 nuclear genes associated with mitochondrial disorders.

Sanger sequencing of multigenic disorders can be technically challenging, time consuming, and prohibitively expensive. High-throughput next-generation sequencing (NGS) can provide a cost-effective method for sequencing targeted genes associated with multigenic disorders. We have developed a NGS clinical targeted gene assay for the mitochondrial genome and for 108 selected nuclear genes associated with mitochondrial disorders. Mitochondrial disorders have a reported incidence of 1 in 5000 live births, encompass a broad range of phenotypes, and are attributed to mutations in the mitochondrial and nuclear genomes. Approximately 20% of mitochondrial disorders result from mutations in mtDNA, with the remaining 80% found in nuclear genes that affect mtDNA levels or mitochondrion protein assembly. In our NGS approach, the 16,569-bp mtDNA is enriched by long-range PCR and the 108 nuclear genes (which represent 1301 amplicons and 680 kb) are enriched by RainDance emulsion PCR. Sequencing is performed on Illumina HiSeq 2000 or MiSeq platforms, and bioinformatics analysis is performed using commercial and in-house developed bioinformatics pipelines. A total of 16 validation and 13 clinical samples were examined. All previously reported variants associated with mitochondrial disorders were found in validation samples, and 5 of the 13 clinical samples were found to have mutations associated with mitochondrial disorders in either the mitochondrial genome or the 108 nuclear genes. All variants were confirmed by Sanger sequencing.

Variant identification in multi-sample pools by illumina genome analyzer sequencing.

Multi-sample pooling and Illumina Genome Analyzer (GA) sequencing allows high throughput sequencing of multiple samples to determine population sequence variation. A preliminary experiment, using the RET proto-oncogene as a model, predicted ≤ 30 samples could be pooled to reliably detect singleton variants without requiring additional confirmation testing. This report used 30 and 50 sample pools to test the hypothesized pooling limit and also to test recent protocol improvements, Illumina GAIIx upgrades, and longer read chemistry. The SequalPrep(TM) method was used to normalize amplicons before pooling. For comparison, a single 'control' sample was run in a different flow cell lane. Data was evaluated by variant read percentages and the subtractive correction method which utilizes the control sample. In total, 59 variants were detected within the pooled samples, which included all 47 known true variants. The 15 known singleton variants due to Sanger sequencing had an average of 1.62 ± 0.26% variant reads for the 30 pool (expected 1.67% for a singleton variant [unique variant within the pool]) and 1.01 ± 0.19% for the 50 pool (expected 1%). The 76 base read lengths had higher error rates than shorter read lengths (33 and 50 base reads), which eliminated the distinction of true singleton variants from background error. This report demonstrated pooling limits from 30 up to 50 samples (depending on error rates and coverage), for reliable singleton variant detection. The presented pooling protocols and analysis methods can be used for variant discovery in other genes, facilitating molecular diagnostic test design and interpretation.

Next generation sequencing for clinical diagnostics-principles and application to targeted resequencing for hypertrophic cardiomyopathy: a paper from the 2009 William Beaumont Hospital Symposium on Molecular Pathology.

During the past five years, new high-throughput DNA sequencing technologies have emerged; these technologies are collectively referred to as next generation sequencing (NGS). By virtue of sequencing clonally amplified DNA templates or single DNA molecules in a massively parallel fashion in a flow cell, NGS provides both qualitative and quantitative sequence data. This combination of information has made NGS the technology of choice for complex genetic analyses that were previously either technically infeasible or cost prohibitive. As a result, NGS has had a fundamental and broad impact on many facets of biomedical research. In contrast, the dissemination of NGS into the clinical diagnostic realm is in its early stages. Though NGS is powerful and can be envisioned to have multiple applications in clinical diagnostics, the technology is currently complex. Successful adoption of NGS into the clinical laboratory will require expertise in both molecular biology techniques and bioinformatics. The current report presents principles that underlie NGS including sequencing library preparation, sequencing chemistries, and an introduction to NGS data analysis. These concepts are subsequently further illustrated by showing representative results from a case study using NGS for targeted resequencing of genes implicated in hypertrophic cardiomyopathy.

Comparison of the Illumina Genome Analyzer and Roche 454 GS FLX for resequencing of hypertrophic cardiomyopathy-associated genes.

Next-generation sequencing (NGS) is widely used in biomedical research, but its adoption has been limited in molecular diagnostics. One application of NGS is the targeted resequencing of genes whose mutations lead to an overlapping clinical phenotype. This study evaluated the comparative performance of the Illumina Genome Analyzer and Roche 454 GS FLX for the resequencing of 16 genes associated with hypertrophic cardiomyopathy (HCM). Using a single human genomic DNA sample enriched by long-range PCR (LR-PCR), 40 GS FLX and 31 Genome Analyzer exon variants were identified using >or=30-fold read-coverage and >or=20% read-percentage selection criteria. Twenty-seven platform concordant variants were Sanger-confirmed. The discordant variants segregated into two categories: variants with read coverages >or=30 on one platform but <30-fold on the alternate platform and variants with read percentages >or=20% on one platform but <20% on the alternate platform. All variants with <30-fold coverage were Sanger-confirmed, suggesting that the coverage criterion of >or=30-fold is too stringent for variant discovery. The variants with <20% read percentage were identified as reference sequence based on Sanger sequencing. These variants were found in homopolymer tracts and short-read misalignments, specifically in genes with high identity. The results of the current study demonstrate the feasibility of combining LR-PCR with the Genome Analyzer or GS FLX for targeted resequencing of HCM-associated genes.

Multi-sample pooling and illumina genome analyzer sequencing methods to determine gene sequence variation for database development.

Determination of sequence variation within a genetic locus to develop clinically relevant databases is critical for molecular assay design and clinical test interpretation, so multisample pooling for Illumina genome analyzer (GA) sequencing was investigated using the RET proto-oncogene as a model. Samples were Sanger-sequenced for RET exons 10, 11, and 13-16. Ten samples with 13 known unique variants ("singleton variants" within the pool) and seven common changes were amplified and then equimolar-pooled before sequencing on a single flow cell lane, generating 36 base reads. For comparison, a single "control" sample was run in a different lane. After alignment, a 24-base quality score-screening threshold and 3; read end trimming of three bases yielded low background error rates with a 27% decrease in aligned read coverage. Sequencing data were evaluated using an established variant detection method (percent variant reads), by the presented subtractive correction method, and with SNPSeeker software. In total, 41 variants (of which 23 were singleton variants) were detected in the 10 pool data, which included all Sanger-identified variants. The 23 singleton variants were detected near the expected 5% allele frequency (average 5.17%+/-0.90% variant reads), well above the highest background error (1.25%). Based on background error rates, read coverage, simulated 30, 40, and 50 sample pool data, expected singleton allele frequencies within pools, and variant detection methods; >or=30 samples (which demonstrated a minimum 1% variant reads for singletons) could be pooled to reliably detect singleton variants by GA sequencing.

Unlabeled probes for the detection and typing of herpes simplex virus.

BACKGROUND: Unlabeled probe detection with a double-stranded DNA (dsDNA) binding dye is one method to detect and confirm target amplification after PCR. Unlabeled probes and amplicon melting have been used to detect small deletions and single-nucleotide polymorphisms in assays where template is in abundance. Unlabeled probes have not been applied to low-level target detection, however. METHODS: Herpes simplex virus (HSV) was chosen as a model to compare the unlabeled probe method to an in-house reference assay using dual-labeled, minor groove binding probes. A saturating dsDNA dye (LCGreen Plus) was used for real-time PCR. HSV-1, HSV-2, and an internal control were differentiated by PCR amplicon and unlabeled probe melting analysis after PCR. RESULTS: The unlabeled probe technique displayed 98% concordance with the reference assay for the detection of HSV from a variety of archived clinical samples (n = 182). HSV typing using unlabeled probes was 99% concordant (n = 104) to sequenced clinical samples and allowed for the detection of sequence polymorphisms in the amplicon and under the probe. CONCLUSIONS: Unlabeled probes and amplicon melting can be used to detect and genotype as few as 10 copies of target per reaction, restricted only by stochastic limitations. The use of unlabeled probes provides an attractive alternative to conventional fluorescence-labeled, probe-based assays for genotyping and detection of HSV and might be useful for other low-copy targets where typing is informative.

Gene silencing in Caenorhabditis elegans by transitive RNA interference.

When a cell is exposed to double-stranded RNA (dsRNA), mRNA from the homologous gene is selectively degraded by a process called RNA interference (RNAi). Here, we provide evidence that dsRNA is amplified in Caenorhabditis elegans to ensure a robust RNAi response. Our data suggest a model in which mRNA targeted by RNAi functions as a template for 5' to 3' synthesis of new dsRNA (termed transitive RNAi). Strikingly, the effect is nonautonomous: dsRNA targeted to a gene expressed in one cell type can lead to transitive RNAi-mediated silencing of a second gene expressed in a distinct cell type. These data suggest dsRNA synthesized in vivo can mediate systemic RNAi.