Scalable detection of technically challenging variants through modified next-generation sequencing.

BACKGROUND: Some clinically important genetic variants are not easily evaluated with next-generation sequencing (NGS) methods due to technical challenges arising from high- similarity copies (e.g., PMS2, SMN1/SMN2, GBA1, HBA1/HBA2, CYP21A2), repetitive short sequences (e.g., ARX polyalanine repeats, FMR1 AGG interruptions in CGG repeats, CFTR poly-T/TG repeats), and other complexities (e.g., MSH2 Boland inversions). METHODS: We customized our NGS processes to detect the technically challenging variants mentioned above with adaptations including target enrichment and bioinformatic masking of similar sequences. Adaptations were validated with samples of known genotypes. RESULTS: Our adaptations provided high-sensitivity and high-specificity detection for most of the variants and provided a high-sensitivity primary assay to be followed with orthogonal disambiguation for the others. The sensitivity of the NGS adaptations was 100% for all of the technically challenging variants. Specificity was 100% for those in PMS2, GBA1, SMN1/SMN2, and HBA1/HBA2, and for the MSH2 Boland inversion; 97.8%-100% for CYP21A2 variants; and 85.7% for ARX polyalanine repeats. CONCLUSIONS: NGS assays can detect technically challenging variants when chemistries and bioinformatics are jointly refined. The adaptations described support a scalable, cost-effective path to identifying all clinically relevant variants within a single sample.

One in seven pathogenic variants can be challenging to detect by NGS: an analysis of 450,000 patients with implications for clinical sensitivity and genetic test implementation.

PURPOSE: To evaluate the impact of technically challenging variants on the implementation, validation, and diagnostic yield of commonly used clinical genetic tests. Such variants include large indels, small copy-number variants (CNVs), complex alterations, and variants in low-complexity or segmentally duplicated regions. METHODS: An interlaboratory pilot study used synthetic specimens to assess detection of challenging variant types by various next-generation sequencing (NGS)-based workflows. One well-performing workflow was further validated and used in clinician-ordered testing of more than 450,000 patients. RESULTS: In the interlaboratory study, only 2 of 13 challenging variants were detected by all 10 workflows, and just 3 workflows detected all 13. Limitations were also observed among 11 less-challenging indels. In clinical testing, 21.6% of patients carried one or more pathogenic variants, of which 13.8% (17,561) were classified as technically challenging. These variants were of diverse types, affecting 556 of 1,217 genes across hereditary cancer, cardiovascular, neurological, pediatric, reproductive carrier screening, and other indicated tests. CONCLUSION: The analytic and clinical sensitivity of NGS workflows can vary considerably, particularly for prevalent, technically challenging variants. This can have important implications for the design and validation of tests (by laboratories) and the selection of tests (by clinicians) for a wide range of clinical indications.

Misattributed parentage as an unanticipated finding during exome/genome sequencing: current clinical laboratory practices and an opportunity for standardization.

PURPOSE: Clinical laboratories performing exome or genome sequencing (ES/GS) are familiar with the challenges associated with proper consenting for and reporting of medically actionable secondary findings based on recommendations from the American College of Medical Genetics and Genomics (ACMG). Misattributed parentage is another type of unanticipated finding a laboratory may encounter during family-based ES/GS; however, there are currently no professional recommendations related to the proper consenting for and reporting of misattributed parentage encountered during ES/GS. METHODS: We surveyed 10 clinical laboratories offering family-based ES/GS regarding their consent language, discovery, and reporting of misattributed parentage. RESULTS: Many laboratories have already developed their own practices/policies for these issues, which do not necessarily agree with those from other labs. CONCLUSION: There are several other possibilities besides true misattributed parentage that could result in similar laboratory findings, and laboratories often feel they lack sufficient information to make formal conclusions on a report regarding the true genetic relatedness of the submitted samples. However, understanding the genetic relatedness (or lack thereof) of the samples submitted for family-based ES/GS has medical relevance. Therefore, professional recommendations for the appropriate handling of suspected misattributed parentage encountered during ES/GS are needed to help standardize current clinical laboratory practices.

Copy-number variants in clinical genome sequencing: deployment and interpretation for rare and undiagnosed disease.

PURPOSE: Current diagnostic testing for genetic disorders involves serial use of specialized assays spanning multiple technologies. In principle, genome sequencing (GS) can detect all genomic pathogenic variant types on a single platform. Here we evaluate copy-number variant (CNV) calling as part of a clinically accredited GS test. METHODS: We performed analytical validation of CNV calling on 17 reference samples, compared the sensitivity of GS-based variants with those from a clinical microarray, and set a bound on precision using orthogonal technologies. We developed a protocol for family-based analysis of GS-based CNV calls, and deployed this across a clinical cohort of 79 rare and undiagnosed cases. RESULTS: We found that CNV calls from GS are at least as sensitive as those from microarrays, while only creating a modest increase in the number of variants interpreted (~10 CNVs per case). We identified clinically significant CNVs in 15% of the first 79 cases analyzed, all of which were confirmed by an orthogonal approach. The pipeline also enabled discovery of a uniparental disomy (UPD) and a 50% mosaic trisomy 14. Directed analysis of select CNVs enabled breakpoint level resolution of genomic rearrangements and phasing of de novo CNVs. CONCLUSION: Robust identification of CNVs by GS is possible within a clinical testing environment.

Successful Application of Whole Genome Sequencing in a Medical Genetics Clinic.

A pilot program was initiated using whole genome sequencing (WGS) to diagnose suspected genetic disorders in the Genetics Clinic at Children's Hospital of Wisconsin. Twenty-two patients underwent WGS between 2010 and 2013. Initially, we obtained a 14% (3/22) diagnosis rate over 2 years; with subsequent reanalysis, this increased to 36% (8/22). Disease causing variants were identified in SKIV2L, CECR1, DGKE, PYCR2, RYR1, PDGFRB, EFTUD2, and BCS1L. In 75% (6/8) of diagnosed cases, the diagnosis affected treatment and/or medical surveillance. Additionally, one case demonstrated a homozygous A18V variant in VLDLR that appears to be associated with a previously undescribed phenotype.

A survey of current practices for genomic sequencing test interpretation and reporting processes in US laboratories.

PURPOSE: While the diagnostic success of genomic sequencing expands, the complexity of this testing should not be overlooked. Numerous laboratory processes are required to support the identification, interpretation, and reporting of clinically significant variants. This study aimed to examine the workflow and reporting procedures among US laboratories to highlight shared practices and identify areas in need of standardization. METHODS: Surveys and follow-up interviews were conducted with laboratories offering exome and/or genome sequencing to support a research program or for routine clinical services. The 73-item survey elicited multiple choice and free-text responses that were later clarified with phone interviews. RESULTS: Twenty-one laboratories participated. Practices highly concordant across all groups included consent documentation, multiperson case review, and enabling patient opt-out of incidental or secondary findings analysis. Noted divergence included use of phenotypic data to inform case analysis and interpretation and reporting of case-specific quality metrics and methods. Few laboratory policies detailed procedures for data reanalysis, data sharing, or patient access to data. CONCLUSION: This study provides an overview of practices and policies of experienced exome and genome sequencing laboratories. The results enable broader consideration of which practices are becoming standard approaches, where divergence remains, and areas of development in best practice guidelines that may be helpful.Genet Med advance online publication 03 Novemeber 2016.

Whole exome sequencing in patients with white matter abnormalities.

Here we report whole exome sequencing (WES) on a cohort of 71 patients with persistently unresolved white matter abnormalities with a suspected diagnosis of leukodystrophy or genetic leukoencephalopathy. WES analyses were performed on trio, or greater, family groups. Diagnostic pathogenic variants were identified in 35% (25 of 71) of patients. Potentially pathogenic variants were identified in clinically relevant genes in a further 7% (5 of 71) of cases, giving a total yield of clinical diagnoses in 42% of individuals. These findings provide evidence that WES can substantially decrease the number of unresolved white matter cases. Ann Neurol 2016;79:1031-1037.

Reducing false-positive incidental findings with ensemble genotyping and logistic regression based variant filtering methods.

As whole genome sequencing (WGS) uncovers variants associated with rare and common diseases, an immediate challenge is to minimize false-positive findings due to sequencing and variant calling errors. False positives can be reduced by combining results from orthogonal sequencing methods, but costly. Here, we present variant filtering approaches using logistic regression (LR) and ensemble genotyping to minimize false positives without sacrificing sensitivity. We evaluated the methods using paired WGS datasets of an extended family prepared using two sequencing platforms and a validated set of variants in NA12878. Using LR or ensemble genotyping based filtering, false-negative rates were significantly reduced by 1.1- to 17.8-fold at the same levels of false discovery rates (5.4% for heterozygous and 4.5% for homozygous single nucleotide variants (SNVs); 30.0% for heterozygous and 18.7% for homozygous insertions; 25.2% for heterozygous and 16.6% for homozygous deletions) compared to the filtering based on genotype quality scores. Moreover, ensemble genotyping excluded > 98% (105,080 of 107,167) of false positives while retaining > 95% (897 of 937) of true positives in de novo mutation (DNM) discovery in NA12878, and performed better than a consensus method using two sequencing platforms. Our proposed methods were effective in prioritizing phenotype-associated variants, and an ensemble genotyping would be essential to minimize false-positive DNM candidates.

Opportunities and challenges associated with clinical diagnostic genome sequencing: a report of the Association for Molecular Pathology.

This report of the Whole Genome Analysis group of the Association for Molecular Pathology illuminates the opportunities and challenges associated with clinical diagnostic genome sequencing. With the reality of clinical application of next-generation sequencing, technical aspects of molecular testing can be accomplished at greater speed and with higher volume, while much information is obtained. Although this testing is a next logical step for molecular pathology laboratories, the potential impact on the diagnostic process and clinical correlations is extraordinary and clinical interpretation will be challenging. We review the rapidly evolving technologies; provide application examples; discuss aspects of clinical utility, ethics, and consent; and address the analytic, postanalytic, and professional implications.

BAP1 loss defines a new class of renal cell carcinoma.

The molecular pathogenesis of renal cell carcinoma (RCC) is poorly understood. Whole-genome and exome sequencing followed by innovative tumorgraft analyses (to accurately determine mutant allele ratios) identified several putative two-hit tumor suppressor genes, including BAP1. The BAP1 protein, a nuclear deubiquitinase, is inactivated in 15% of clear cell RCCs. BAP1 cofractionates with and binds to HCF-1 in tumorgrafts. Mutations disrupting the HCF-1 binding motif impair BAP1-mediated suppression of cell proliferation but not deubiquitination of monoubiquitinated histone 2A lysine 119 (H2AK119ub1). BAP1 loss sensitizes RCC cells in vitro to genotoxic stress. Notably, mutations in BAP1 and PBRM1 anticorrelate in tumors (P = 3 × 10(-5)), [corrected] and combined loss of BAP1 and PBRM1 in a few RCCs was associated with rhabdoid features (q = 0.0007). BAP1 and PBRM1 regulate seemingly different gene expression programs, and BAP1 loss was associated with high tumor grade (q = 0.0005). Our results establish the foundation for an integrated pathological and molecular genetic classification of RCC, paving the way for subtype-specific treatments exploiting genetic vulnerabilities.

Platelet hexosaminidase a enzyme assay effectively detects carriers missed by targeted DNA mutation analysis.

Biochemical testing of hexosaminidase A (HexA) enzyme activity has been available for decades and has the ability to detect almost all Tay-Sachs disease (TSD) carriers, irrespective of ethnic background. This is increasingly important, as the gene pool of those who identify as Ashkenazi Jewish is diversifying. Here we describe the analysis of a cohort of 4,325 individuals arising from large carrier screening programs and tested by the serum and/or platelet HexA enzyme assays and by targeted DNA mutation analysis. Our results continue to support the platelet assay as a highly effective method for TSD carrier screening, with a low inconclusive rate and the ability to detect possible disease-causing mutation carriers that would have been missed by targeted DNA mutation analysis. Sequence analysis performed on one such platelet assay carrier, who had one non-Ashkenazi Jewish parent, identified the amino acid change Thr259Ala (A775G). Based on crystallographic modeling, this change is predicted to be deleterious, as threonine 259 is positioned proximal to the HexA alpha subunit active site and helps to stabilize key residues therein. Accordingly, if individuals are screened for TSD in broad-based programs by targeted molecular testing alone, they must be made aware that there is a more sensitive and inexpensive test available that can identify additional carriers. Alternatively, the enzyme assays can be offered as a first tier test, especially when screening individuals of mixed or non-Jewish ancestry.

The role of the F508C mutation in congenital bilateral absence of the vas deferens.

PURPOSE: Congenital bilateral absence of the vas deferens is a pathologic condition associated with normal spermatogenesis, azoospermia, and lack of both vasa deferentia. A significant association between mutations in the cystic fibrosis transmembrane conductance regulator gene among men with congenital bilateral absence of the vas deferens has been established. The objective of this study was to determine whether the F508C variant in the cystic fibrosis transmembrane conductance regulator gene has a significant effect on congenital bilateral absence of the vas deferens prevalence, when present in conjunction with a second cystic fibrosis transmembrane conductance regulator disease causing mutation. METHODS AND RESULTS: We compared the frequency of F508C in male subjects submitted for diagnostic testing on suspicion of cystic fibrosis or during cystic fibrosis carrier screening, to men with a clinical diagnosis of congenital bilateral absence of the vas deferens. Although frequencies of F508C did not vary significantly between 850 individuals undergoing cystic fibrosis carrier screening and those submitted for diagnostic testing on suspicion of cystic fibrosis, the frequency of F508C in the congenital bilateral absence of the vas deferens population was significantly higher than expected (chi2 = 6.95, corrected P = 0.0486). CONCLUSION: We conclude that the F508C variant in cystic fibrosis transmembrane conductance regulator may represent a pathogenic defect and lead to congenital bilateral absence of the vas deferens when combined with a second cystic fibrosis transmembrane conductance regulator mutation.

Patterns of synonymous codon usage in Drosophila melanogaster genes with sex-biased expression.

The nonrandom use of synonymous codons (codon bias) is a well-established phenomenon in Drosophila. Recent reports suggest that levels of codon bias differ among genes that are differentially expressed between the sexes, with male-expressed genes showing less codon bias than female-expressed genes. To examine the relationship between sex-biased gene expression and level of codon bias on a genomic scale, we surveyed synonymous codon usage in 7276 D. melanogaster genes that were classified as male-, female-, or non-sex-biased in their expression in microarray experiments. We found that male-biased genes have significantly less codon bias than both female- and non-sex-biased genes. This pattern holds for both germline and somatically expressed genes. Furthermore, we find a significantly negative correlation between level of codon bias and degree of sex-biased expression for male-biased genes. In contrast, female-biased genes do not differ from non-sex-biased genes in their level of codon bias and show a significantly positive correlation between codon bias and degree of sex-biased expression. These observations cannot be explained by differences in chromosomal distribution, mutational processes, recombinational environment, gene length, or absolute expression level among genes of the different expression classes. We propose that the observed codon bias differences result from differences in selection at synonymous and/or linked nonsynonymous sites between genes with male- and female-biased expression.

Population genetic analysis of Bartonella bacilliformis isolates from areas of peru where Carrion's disease is endemic and epidemic.

Carrion's disease is caused by infection with the alpha-proteobacterium Bartonella bacilliformis. Distribution of the disease is considered coincident with the distribution of its known vector, the sand fly Lutzomyia verrucarum. Recent epidemics of B. bacilliformis infections associated with atypical symptomatology in nonendemic regions have raised questions regarding the historic and present distribution of this bacterium and the scope of disease that infection causes. Phylogenetic relationships and genomic diversity of 18 B. bacilliformis isolates (10 isolates from a region where Carrion's disease is epidemic, Cuzco, Peru, and 8 isolates from a region where Carrion's disease is endemic, Caraz, Peru) were assessed using genomic data generated by infrequent restriction site PCR and gene sequence analysis of the flagellin gltA and ialB genes. A population genetic analysis of the genomic diversity suggests that what was once considered an epidemic region of Peru did not result from the recent introduction of B. bacilliformis.

Molecular evolution of sex-biased genes in Drosophila.

Studies of morphology, interspecific hybridization, protein/DNA sequences, and levels of gene expression have suggested that sex-related characters (particularly those involved in male reproduction) evolve rapidly relative to non-sex-related characters. Here we report a general comparison of evolutionary rates of sex-biased genes using data from cDNA microarray experiments and comparative genomic studies of Drosophila. Comparisons of nonsynonymous/synonymous substitution rates (d(N)/d(S)) between species of the D. melanogaster subgroup revealed that genes with male-biased expression had significantly faster rates of evolution than genes with female-biased or unbiased expression. The difference was caused primarily by a higher d(N) in the male-biased genes. The same pattern was observed for comparisons among more distantly related species. In comparisons between D. melanogaster and D. pseudoobscura, genes with highly biased male expression were significantly more divergent than genes with highly biased female expression. In many cases, orthologs of D. melanogaster male-biased genes could not be identified in D. pseudoobscura through a Blast search. In contrast to the male-biased genes, there was no clear evidence for accelerated rates of evolution in female-biased genes, and most comparisons indicated a reduced rate of evolution in female-biased genes relative to unbiased genes. Male-biased genes did not show an increased ratio of nonsynonymous/synonymous polymorphism within D. melanogaster, and comparisons of polymorphism/divergence ratios suggest that the rapid evolution of male-biased genes is caused by positive selection.

Enhanced selection for MHC diversity in social tuco-tucos.

To explore the effects of behavior and demography on balancing selection at major histocompatibility complex (MHC) loci, we examined allelic diversity at exon 2 of the MHC class II DQbeta locus in a social and a solitary species of tuco-tuco (Rodentia: Ctenomyidae: Ctenomys), both of which occur in the same valley in southwestern Argentina. By comparing patterns of diversity at this MHC gene to the diversity evident at fifteen microsatellite loci, we demonstrate that balancing selection at the DQbeta locus is enhanced in the social species compared to its solitary congener. These findings have intriguing implications for the role of behavioral and demographic parameters in maintaining diversity at MHC loci.