Algorithmic improvements for discovery of germline copy number variants in next-generation sequencing data.

BACKGROUND: Copy number variants (CNVs) play a significant role in human heredity and disease. However, sensitive and specific characterization of germline CNVs from NGS data has remained challenging, particularly for hybridization-capture data in which read counts are the primary source of copy number information. RESULTS: We describe two algorithmic adaptations that improve CNV detection accuracy in a Hidden Markov Model (HMM) context. First, we present a method for computing target- and copy number-specific emission distributions. Second, we demonstrate that the Pointwise Maximum a posteriori (PMAP) HMM decoding procedure yields improved sensitivity for small CNV calls compared to the more common Viterbi HMM decoder. We develop a prototype implementation, called Cobalt, and compare it to other CNV detection tools using sets of simulated and previously detected CNVs with sizes spanning a single exon to a full chromosome. CONCLUSIONS: In both the simulation and previously detected CNV studies Cobalt shows similar sensitivity but significantly fewer false positive detections compared to other callers. Overall sensitivity is 80-90% for deletion CNVs spanning 1-4 targets and 90-100% for larger deletion events, while sensitivity is somewhat lower for small duplication CNVs.

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 Phevor's 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.

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.

Identification of patients with RAG mutations previously diagnosed with common variable immunodeficiency disorders.

PURPOSE: Combined 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. METHODS: Two 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. RESULTS: Two 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. CONCLUSION: Astute 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.

VarBin, a novel method for classifying true and false positive variants in NGS data.

BACKGROUND: Variant discovery for rare genetic diseases using Illumina genome or exome sequencing involves screening of up to millions of variants to find only the one or few causative variant(s). Sequencing or alignment errors create "false positive" variants, which are often retained in the variant screening process. Methods to remove false positive variants often retain many false positive variants. This report presents VarBin, a method to prioritize variants based on a false positive variant likelihood prediction. METHODS: VarBin uses the Genome Analysis Toolkit variant calling software to calculate the variant-to-wild type genotype likelihood ratio at each variant change and position divided by read depth. The resulting Phred-scaled, likelihood-ratio by depth (PLRD) was used to segregate variants into 4 Bins with Bin 1 variants most likely true and Bin 4 most likely false positive. PLRD values were calculated for a proband of interest and 41 additional Illumina HiSeq, exome and whole genome samples (proband's family or unrelated samples). At variant sites without apparent sequencing or alignment error, wild type/non-variant calls cluster near -3 PLRD and variant calls typically cluster above 10 PLRD. Sites with systematic variant calling problems (evident by variant quality scores and biases as well as displayed on the iGV viewer) tend to have higher and more variable wild type/non-variant PLRD values. Depending on the separation of a proband's variant PLRD value from the cluster of wild type/non-variant PLRD values for background samples at the same variant change and position, the VarBin method's classification is assigned to each proband variant (Bin 1 to Bin 4). RESULTS: To assess VarBin performance, Sanger sequencing was performed on 98 variants in the proband and background samples. True variants were confirmed in 97% of Bin 1 variants, 30% of Bin 2, and 0% of Bin 3/Bin 4. CONCLUSIONS: These data indicate that VarBin correctly classifies the majority of true variants as Bin 1 and Bin 3/4 contained only false positive variants. The "uncertain" Bin 2 contained both true and false positive variants. Future work will further differentiate the variants in Bin 2.

Morphologic Characteristics and Mutational Analysis of Fumarate Hydratase Deficient Kidney and Smooth Muscle Tumors.

OBJECTIVES: Fumarate hydratase (FH)-deficient tumors can occur due to germline or somatic mutations and have distinctive morphologic features. The aims of this study are to refine morphologic criteria and identify mutations in FH-deficient smooth muscle tumors (SMTs). METHODS: The morphology of SMTs and kidney tumors submitted to a national reference laboratory for FH immunohistochemistry (IHC) was reviewed by two gynecologic and two genitourinary pathologists, respectively. Fisher exact test was used for analysis. Fourteen SMTs were sequenced using the Illumina TruSight Oncology 500 Assay. RESULTS: Twenty-two kidney tumors (5 FH deficient) and 51 SMTs (27 FH deficient) were reviewed. FH-deficient kidney tumors exclusively showed cord-like growth, rhabdoid change, and absence of coagulative tumor necrosis and psammoma bodies. FH-deficient SMTs were significantly more likely to have staghorn vessels, eosinophilic cytoplasmic inclusions, schwannoma-like areas, or hereditary leiomyomatosis and renal cell cancer-like nuclei (P < .05 for each). Seven of 14 sequenced SMTs showed mutations of the FH gene and no other driver mutations. CONCLUSIONS: FH-deficient SMTs submitted for FH immunohistochemistry (IHC) showed distinct morphology. Although FH IHC is used for screening of FH-deficient tumors, FH mutations were identified in only 50% of FH-deficient SMTs. This highlights the need for additional exploration of mechanisms of FH protein loss in tumors lacking FH mutations.

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.

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.

Genome-Wide Copy Number Variation Detection Using NGS: Data Analysis and Interpretation.

Copy number variants (CNVs) and copy neutral loss of heterozygosity (CN-LOH) represent important types of genomic abnormalities in cancer. Genomic DNA microarray serves as the current gold standard method for detecting genome-wide CNVs and CN-LOH. However, as next-generation sequencing (NGS) is widely used to detect gene variants in clinical testing, the ability of NGS to detect CNVs and CN-LOH has also been demonstrated. This chapter describes a protocol for detecting genome-wide large somatic CNVs and CN-LOH using a single nucleotide polymorphism (SNP) sequencing backbone. When combined with a targeted gene mutation panel, this strategy allows for simultaneous detection of somatic gene mutations and genome-wide CNVs and CN-LOH.

Novel PLP1 Mutations Identified With Next-Generation Sequencing Expand the Spectrum of PLP1-Associated Leukodystrophy Clinical Phenotypes.

Next-generation sequencing was performed for 2 families with an undiagnosed neurologic disease. Analysis revealed X-linked mutations in the proteolipid protein 1 (PLP1) gene, which is associated with X-linked Pelizaeus-Merzbacher disease and Spastic Paraplegia type 2. In family A, the novel PLP1 missense mutation c.617T>A (p.M206K) was hemizygous in the 2 affected male children and heterozygous in the mother. In family B, the novel de novoPLP1 frameshift mutation c.359_369del (p.G120fs) was hemizygous in the affected male child. Although PLP1 mutations have been reported to cause an increasingly wide range of phenotypes inclusive of the dystonia, spastic paraparesis, motor neuronopathy, and leukodystrophy observed in our patients, atypical features included the cerebrospinal fluid deficiency of neurotransmitter and pterin metabolites and the delayed appearance of myelin abnormalities on neuroimaging studies. Next-generation sequencing studies provided a diagnosis for these families with complex leukodystrophy disease phenotypes, which expanded the spectrum of PLP1-associated leukodystrophy clinical phenotypes.

Unlabeled oligonucleotides as internal temperature controls for genotyping by amplicon melting.

Amplicon melting is a closed-tube method for genotyping that does not require probes, real-time analysis, or allele-specific polymerase chain reaction. However, correct differentiation of homozygous mutant and wild-type samples by melting temperature (Tm) requires high-resolution melting and closely controlled reaction conditions. When three different DNA extraction methods were used to isolate DNA from whole blood, amplicon Tm differences of 0.03 to 0.39 degrees C attributable to the extractions were observed. To correct for solution chemistry differences between samples, complementary unlabeled oligonucleotides were included as internal temperature controls to shift and scale the temperature axis of derivative melting plots. This adjustment was applied to a duplex amplicon melting assay for the methylenetetrahydrofolate reductase variants 1298A>C and 677C>T. High- and low-temperature controls bracketing the amplicon melting region decreased the Tm SD within homozygous genotypes by 47 to 82%. The amplicon melting assay was 100% concordant to an adjacent hybridization probe (HybProbe) melting assay when temperature controls were included, whereas a 3% error rate was observed without temperature correction. In conclusion, internal temperature controls increase the accuracy of genotyping by high-resolution amplicon melting and should also improve results on lower resolution instruments.

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

Increased sensitivity of bacterial detection in cerebrospinal fluid by fluorescent staining on low-fluorescence membrane filters.

A membrane-filter-based, fluorescent Gram stain method for bacterial detection in cerebrospinal fluid samples was developed and evaluated as a rapid, sensitive alternative to standard Gram stain protocols. A recently developed, modified version of the aluminium oxide membrane Anopore with low-fluorescence optical properties showed superior performance in this application. Other aspects of the fluorescent Gram stain system that were evaluated include membrane filter selection, strategies to reduce fluorescence fading and the effect of patient blood cells on bacterial detection in the fluorescently stained cerebrospinal fluid samples. The combination of the membrane filter's bacteria-concentrating ability and absolute retention along with high-contrast, fluorescent Gram discriminating dyes enabled rapid bacterial detection and Gram discrimination, with a 1-1.5 order of magnitude increase in the bacterial concentration limit of detection.

Next-generation sequencing of custom amplicons to improve coverage of HaloPlex multigene panels.

Next-generation sequencing (NGS) of multigene panels performed for genetic clinical diagnostics requires 100% coverage of all targeted genes. In the genetic diagnostics laboratory, coverage gaps are typically filled with Sanger sequencing after NGS data are collected and analyzed. Libraries prepared using the hybridization-based custom capture HaloPlex method are covered at ~98% and include gaps in coverage because of the location of the restriction enzyme sites used for fragmentation and differences in the designed and actual library insert size. We describe a method for improving the coverage of HaloPlex libraries by generating a set of amplicons spanning known low-coverage regions that are pooled, indexed by sample, and sequenced together with the HaloPlex libraries. This approach reduces the number of post-NGS Sanger sequencing reactions required and complements any NGS library preparation method when complete gene coverage is necessary.

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.

Developing genome and exome sequencing for candidate gene identification in inherited disorders: an integrated technical and bioinformatics approach.

CONTEXT: Advances in sequencing technology with the commercialization of next-generation sequencing (NGS) has substantially increased the feasibility of sequencing human genomes and exomes. Next-generation sequencing has been successfully applied to the discovery of disease-causing genes in rare, inherited disorders. By necessity, the advent of NGS has fostered the concurrent development of bioinformatics approaches to expeditiously analyze the large data sets generated. Next-generation sequencing has been used for important discoveries in the research setting and is now being implemented into the clinical diagnostic arena. OBJECTIVE: To review the current literature on technical and bioinformatics approaches for exome and genome sequencing and highlight examples of successful disease gene discovery in inherited disorders. To discuss the challenges for implementing NGS in the clinical research and diagnostic arenas. DATA SOURCES: Literature review and authors' experience. CONCLUSIONS: Next-generation sequencing approaches are powerful and require an investment in infrastructure and personnel expertise for effective use; however, the potential for improvement of patient care through faster and more accurate molecular diagnoses is high.

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

BACKGROUND: Aortopathies 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. METHODS: The 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. RESULTS: Read 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. CONCLUSIONS: When 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.