Jacob D. Durtschi

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.

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.

Germline Mutations in NFKB2 Implicate the Noncanonical NF-κB Pathway in the Pathogenesis of Common Variable Immunodeficiency

Common variable immunodeficiency (CVID) is a heterogeneous disorder characterized by antibody deficiency, poor humoral response to antigens, and recurrent infections. To investigate the molecular cause of CVID, we carried out exome sequence analysis of a family diagnosed with CVID and identified a heterozygous frameshift mutation, c.2564delA (p.Lys855Serfs(∗)7), in NFKB2 affecting the C terminus of NF-κB2 (also known as p100/p52 or p100/p49). Subsequent screening of NFKB2 in 33 unrelated CVID-affected individuals uncovered a second heterozygous nonsense mutation, c.2557C>T (p.Arg853(∗)), in one simplex case. Affected individuals in both families presented with an unusual combination of childhood-onset hypogammaglobulinemia with recurrent infections, autoimmune features, and adrenal insufficiency. NF-κB2 is the principal protein involved in the noncanonical NF-κB pathway, is evolutionarily conserved, and functions in peripheral lymphoid organ development, B cell development, and antibody production. In addition, Nfkb2 mouse models demonstrate a CVID-like phenotype with hypogammaglobulinemia and poor humoral response to antigens. Immunoblot analysis and immunofluorescence microscopy of transformed B cells from affected individuals show that the NFKB2 mutations affect phosphorylation and proteasomal processing of p100 and, ultimately, p52 nuclear translocation. These findings describe germline mutations in NFKB2 and establish the noncanonical NF-κB signaling pathway as a genetic etiology for this primary immunodeficiency syndrome.

Phevor Combines Multiple Biomedical Ontologies for Accurate Identification of Disease-Causing Alleles in Single Individuals and Small Nuclear Families

Phevor integrates phenotype, gene function, and disease information with personal genomic data for improved power to identify disease-causing alleles. Phevor works by combining knowledge resident in multiple biomedical ontologies with the outputs of variant-prioritization tools. It does so by using an algorithm that propagates information across and between ontologies. This process enables Phevor to accurately reprioritize potentially damaging alleles identified by variant-prioritization tools in light of gene function, disease, and phenotype knowledge. Phevor is especially useful for single-exome and family-trio-based diagnostic analyses, the most commonly occurring clinical scenarios and ones for which existing personal genome diagnostic tools are most inaccurate and underpowered. Here, we present a series of benchmark analyses illustrating Phevors performance characteristics. Also presented are three recent Utah Genome Project case studies in which Phevor was used to identify disease-causing alleles. Collectively, these results show that Phevor improves diagnostic accuracy not only for individuals presenting with established disease phenotypes but also for those with previously undescribed and atypical disease presentations. Importantly, Phevor is not limited to known diseases or known disease-causing alleles. As we demonstrate, Phevor can also use latent information in ontologies to discover genes and disease-causing alleles not previously associated with disease.

A unified test of linkage analysis and rare-variant association for analysis of pedigree sequence data

High-throughput sequencing of related individuals has become an important tool for studying human disease. However, owing to technical complexity and lack of available tools, most pedigree-based sequencing studies rely on an ad hoc combination of suboptimal analyses. Here we present pedigree-VAAST (pVAAST), a disease-gene identification tool designed for high-throughput sequence data in pedigrees. pVAAST uses a sequence-based model to perform variant and gene-based linkage analysis. Linkage information is then combined with functional prediction and rare variant case-control association information in a unified statistical framework. pVAAST outperformed linkage and rare-variant association tests in simulations and identified disease-causing genes from whole-genome sequence data in three human pedigrees with dominant, recessive and de novo inheritance patterns. The approach is robust to incomplete penetrance and locus heterogeneity and is applicable to a wide variety of genetic traits. pVAAST maintains high power across studies of monogenic, high-penetrance phenotypes in a single pedigree to highly polygenic, common phenotypes involving hundreds of pedigrees.

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).

Identification of Patients with RAG Mutations Previously Diagnosed with Common Variable Immunodeficiency Disorders

PurposeCombined immunodeficiency (CID) presents a unique challenge to clinicians. Two patients presented with the prior clinical diagnosis of common variable immunodeficiency (CVID) disorder marked by an early age of presentation, opportunistic infections, and persistent lymphopenia. Due to the presence of atypical clinical features, next generation sequencing was applied documenting RAG deficiency in both patients.MethodsTwo different genetic analysis techniques were applied in these patients including whole exome sequencing in one patient and the use of a gene panel designed to target genes known to cause primary immunodeficiency disorders (PIDD) in a second patient. Sanger dideoxy sequencing was used to confirm RAG1 mutations in both patients.ResultsTwo young adults with a history of recurrent bacterial sinopulmonary infections, viral infections, and autoimmune disease as well as progressive hypogammaglobulinemia, abnormal antibody responses, lymphopenia and a prior diagnosis of CVID disorder were evaluated. Compound heterozygous mutations in RAG1 (1) c256_257delAA, p86VfsX32 and (2) c1835A>G, pH612R were documented in one patient. Compound heterozygous mutations in RAG1 (1) c.1566G>T, p.W522C and (2) c.2689C>T, p. R897X) were documented in a second patient post-mortem following a fatal opportunistic infection.ConclusionAstute clinical judgment in the evaluation of patients with PIDD is necessary. Atypical clinical findings such as early onset, granulomatous disease, or opportunistic infections should support the consideration of atypical forms of late onset CID secondary to RAG deficiency. Next generation sequencing approaches provide powerful tools in the investigation of these patients and may expedite definitive treatments.

A novel germline PIGA mutation in Ferro-Cerebro-Cutaneous syndrome: A neurodegenerative X-linked epileptic encephalopathy with systemic iron-overload

Three related males presented with a newly recognized x‐linked syndrome associated with neurodegeneration, cutaneous abnormalities, and systemic iron overload. Linkage studies demonstrated that they shared a haplotype on Xp21.3–Xp22.2 and exome sequencing was used to identify candidate variants. Of the segregating variants, only a PIGA mutation segregated with disease in the family. The c.328_330delCCT PIGA variant predicts, p.Leu110del (or c.1030_1032delCTT, p.Leu344del depending on the reference sequence). The unaffected great‐grandfather shared his X allele with the proband but he did not have the PIGA mutation, indicating that the mutation arose de novo in his daughter. A single family with a germline PIGA mutation has been reported; affected males had a phenotype characterized by multiple congenital anomalies and severe neurologic impairment resulting in infantile lethality. In contrast, affected boys in the family described here were born without anomalies and were neurologically normal prior to onset of seizures after 6 months of age, with two surviving to the second decade. PIGA encodes an enzyme in the GPI anchor biosynthesis pathway. An affected individual in the family studied here was deficient in GPI anchor proteins on granulocytes but not erythrocytes. In conclusion, the PIGA mutation in this family likely causes a reduction in GPI anchor protein cell surface expression in various cell types, resulting in the observed pleiotropic phenotype involving central nervous system, skin, and iron metabolism.

A direct comparison of next generation sequencing enrichment methods using an aortopathy gene panel- clinical diagnostics perspective

BackgroundAortopathies are a group of disorders characterized by aneurysms, dilation, and tortuosity of the aorta. Because of the phenotypic overlap and genetic heterogeneity of diseases featuring aortopathy, molecular testing is often required for timely and correct diagnosis of affected individuals. In this setting next generation sequencing (NGS) offers several advantages over traditional molecular techniques.MethodsThe purpose of our study was to compare NGS enrichment methods for a clinical assay targeting the nine genes known to be associated with aortopathy. RainDance emulsion PCR and SureSelect RNA-bait hybridization capture enrichment methods were directly compared by enriching DNA from eight samples. Enriched samples were barcoded, pooled, and sequenced on the Illumina HiSeq2000 platform. Depth of coverage, consistency of coverage across samples, and the overlap of variants identified were assessed. This data was also compared to whole-exome sequencing data from ten individuals.ResultsRead depth was greater and less variable among samples that had been enriched using the RNA-bait hybridization capture enrichment method. In addition, samples enriched by hybridization capture had fewer exons with mean coverage less than 10, reducing the need for followup Sanger sequencing. Variants sets produced were 77% concordant, with both techniques yielding similar numbers of discordant variants.ConclusionsWhen comparing the design flexibility, performance, and cost of the targeted enrichment methods to whole-exome sequencing, the RNA-bait hybridization capture enrichment gene panel offers the better solution for interrogating the aortopathy genes in a clinical laboratory setting.

Mycobacterium chelonae-abscessus Complex Associated with Sinopulmonary Disease, Northeastern USA

Members of the Mycobacterium chelonae-abscessus complex represent Mycobacterium species that cause invasive infections in immunocompetent and immunocompromised hosts. We report the detection of a new pathogen that had been misidentified as M. chelonae with an atypical antimicrobial drug susceptibility profile. The discovery prompted a multicenter investigation of 26 patients. Almost all patients were from the northeastern United States, and most had underlying sinus or pulmonary disease. Infected patients had clinical features similar to those with M. abscessus infections. Taxonomically, the new pathogen shared molecular identity with members of the M. chelonae-abscessus complex. Multilocus DNA target sequencing, DNA-DNA hybridization, and deep multilocus sequencing (43 full-length genes) support a new taxon for these microorganisms. Because most isolates originated in Pennsylvania, we propose the name M. franklinii sp. nov. This investigation underscores the need for accurate identification of Mycobacterium spp. to detect new pathogens implicated in human disease.