Enhancing clinical genomic accuracy with panelGC: a novel metric and tool for quantifying and monitoring GC biases in hybridization capture panel sequencing.

Accurate assessment of fragment abundance within a genome is crucial in clinical genomics applications such as the analysis of copy number variation (CNV). However, this task is often hindered by biased coverage in regions with varying guanine-cytosine (GC) content. These biases are particularly exacerbated in hybridization capture sequencing due to GC effects on probe hybridization and polymerase chain reaction (PCR) amplification efficiency. Such GC content-associated variations can exert a negative impact on the fidelity of CNV calling within hybridization capture panels. In this report, we present panelGC, a novel metric, to quantify and monitor GC biases in hybridization capture sequencing data. We establish the efficacy of panelGC, demonstrating its proficiency in identifying and flagging potential procedural anomalies, even in situations where instrument and experimental monitoring data may not be readily accessible. Validation using real-world datasets demonstrates that panelGC enhances the quality control and reliability of hybridization capture panel sequencing.

Data sharing to improve concordance in variant interpretation across laboratories: results from the Canadian Open Genetics Repository.

BACKGROUND: This study aimed to identify and resolve discordant variant interpretations across clinical molecular genetic laboratories through the Canadian Open Genetics Repository (COGR), an online collaborative effort for variant sharing and interpretation. METHODS: Laboratories uploaded variant data to the Franklin Genoox platform. Reports were issued to each laboratory, summarising variants where conflicting classifications with another laboratory were noted. Laboratories could then reassess variants to resolve discordances. Discordance was calculated using a five-tier model (pathogenic (P), likely pathogenic (LP), variant of uncertain significance (VUS), likely benign (LB), benign (B)), a three-tier model (LP/P are positive, VUS are inconclusive, LB/B are negative) and a two-tier model (LP/P are clinically actionable, VUS/LB/B are not). We compared the COGR classifications to automated classifications generated by Franklin. RESULTS: Twelve laboratories submitted classifications for 44 510 unique variants. 2419 variants (5.4%) were classified by two or more laboratories. From baseline to after reassessment, the number of discordant variants decreased from 833 (34.4% of variants reported by two or more laboratories) to 723 (29.9%) based on the five-tier model, 403 (16.7%) to 279 (11.5%) based on the three-tier model and 77 (3.2%) to 37 (1.5%) based on the two-tier model. Compared with the COGR classification, the automated Franklin classifications had 94.5% sensitivity and 96.6% specificity for identifying actionable (P or LP) variants. CONCLUSIONS: The COGR provides a standardised mechanism for laboratories to identify discordant variant interpretations and reduce discordance in genetic test result delivery. Such quality assurance programmes are important as genetic testing is implemented more widely in clinical care.

Data sharing as a national quality improvement program: reporting on BRCA1 and BRCA2 variant-interpretation comparisons through the Canadian Open Genetics Repository (COGR).

PurposeThe purpose of this study was to develop a national program for Canadian diagnostic laboratories to compare DNA-variant interpretations and resolve discordant-variant classifications using the BRCA1 and BRCA2 genes as a case study.MethodsBRCA1 and BRCA2 variant data were uploaded and shared through the Canadian Open Genetics Repository (COGR; http://www.opengenetics.ca). A total of 5,554 variant observations were submitted; classification differences were identified and comparison reports were sent to participating laboratories. Each site had the opportunity to reclassify variants. The data were analyzed before and after the comparison report process to track concordant- or discordant-variant classifications by three different models.ResultsVariant-discordance rates varied by classification model: 38.9% of variants were discordant when using a five-tier model, 26.7% with a three-tier model, and 5.0% with a two-tier model. After the comparison report process, the proportion of discordant variants dropped to 30.7% with the five-tier model, to 14.2% with the three-tier model, and to 0.9% using the two-tier model.ConclusionWe present a Canadian interinstitutional quality improvement program for DNA-variant interpretations. Sharing of variant knowledge by clinical diagnostic laboratories will allow clinicians and patients to make more informed decisions and lead to better patient outcomes.

Real-world cost-effectiveness of panel-based genomic testing to inform therapeutic decisions for metastatic colorectal cancer.

BACKGROUND: Mutations in KRAS and NRAS are associated with a lack of response to cetuximab and panitumumab, two biologics used for third-line therapy of metastatic colorectal cancer (mCRC). In British Columbia, Canada, eligibility for cetuximab or panitumumab was first based on single-gene KRAS testing. OncoPanel, a multi-gene next-generation sequencing panel with both KRAS and NRAS, was introduced in 2016. Our objective was to estimate the real-world cost-effectiveness of OncoPanel versus to single-gene KRAS testing to inform eligibility for cetuximab or panitumumab in mCRC. METHODS: Using population-based administrative health data, we identified a cohort of mCRC patients who had received a KRAS or OncoPanel test, and completed prior chemotherapy in 2010-2019. We matched KRAS- and OncoPanel-tested patients (1:1) using genetic matching to balance baseline covariates. Mean and incremental 3-year costs, survival, and quality-adjusted survival were estimated using inverse-probability-of-censoring weighting and bootstrapping. We conducted scenario-based sensitivity analysis for key costs and assumptions. FINDINGS: All OncoPanel-tested cases (n=371) were matched to a KRAS-tested comparator. In the KRAS and OncoPanel groups, respectively, 55·8 % and 41·2 % of patients were potentially eligible for cetuximab or panitumumab based on mutation status. Incremental cost and effectiveness of OncoPanel were $72 (95 % CI: -6387, 6107), -0·004 life-years (95 % CI: -0·119, 0·113), and -0·011 quality-adjusted life-years (95 % CI: -0·094, 0·075). Reductions in systemic therapy costs were offset by increased costs in other resources. Results were moderately sensitive to time horizon and changes in testing or treatment cost. INTERPRETATION: The use of OncoPanel resulted in more precise targeting of cetuximab and panitumumab, but there was no change in incremental cost or quality-adjusted survival. Understanding the balance of costs achieved in practice can provide insight into the effect of future changes in testing policy, test cost, treatment eligibility, or drug prices in this setting.

Real-World Clinical Outcomes for Patients with EGFR and HER2 Exon 20 Insertion-Mutated Non-Small-Cell Lung Cancer.

(1) Background: Exon 20 insertion mutations (ex20ins) in EGFR and HER2 are uncommon driver mutations in non-small-cell lung cancer (NSCLC), with a poor prognosis and few targeted therapy options, and there are limited real-world data. Here, we report the clinicopathologic features and outcomes for patients with ex20ins NSCLC across British Columbia, Canada. (2) Methods: NSCLC patients with ex20ins in EGFR or HER2 were identified via tumour testing between 1 January 2016 and 31 December 2021 (n = 7233). Data were collected by chart review. Survival analyses were performed using the Kaplan-Meier method using the log-rank test. (3) Results: A total of 131 patients were identified. The median age was 66. Thirty-three percent of patients had brain metastases. For the EGFR cohort, the median OS was 18.6 months for patients who received any systemic therapy (ST) vs. 2.6 months for patients who did not (p < 0.001). Median OS was similar for patients treated with ex20ins-specific tyrosine kinase inhibitors (TKIs) vs. other STs (18.6 vs. 15.9 months; p = 0.463). The median first-line PFS was 4.1 vs. 7.4 months for patients treated with a TKI vs. other ST (p = 0.744). For the HER2 cohort, the median OS was 9.0 months for patients who received any ST vs. 4.9 months for patients who did not (p = 0.015). The median OS was 23.0 months for patients treated with an ex20ins TKI vs. 5.6 months for patients who were not (p = 0.019). The median first-line PFS was 5.4 vs. 2.1 months for patients treated with a TKI vs. other ST (p = 0.343). (4) Conclusions: Overall survival was significantly longer among ex20ins patients who received any systemic therapy vs. those who did not. Overall survival was significantly better among HER2 ex20ins patients who received ex20ins-specific TKIs.

Defining the heterogeneity of unbalanced structural variation underlying breast cancer susceptibility by nanopore genome sequencing.

Germline structural variants (SVs) are challenging to resolve by conventional genetic testing assays. Long-read sequencing has improved the global characterization of SVs, but its sensitivity at cancer susceptibility loci has not been reported. Nanopore long-read genome sequencing was performed for nineteen individuals with pathogenic copy number alterations in BRCA1, BRCA2, CHEK2 and PALB2 identified by prior clinical testing. Fourteen variants, which spanned single exons to whole genes and included a tandem duplication, were accurately represented. Defining the precise breakpoints of SVs in BRCA1 and CHEK2 revealed unforeseen allelic heterogeneity and informed the mechanisms underlying the formation of recurrent deletions. Integrating read-based and statistical phasing further helped define extended haplotypes associated with founder alleles. Long-read sequencing is a sensitive method for characterizing private, recurrent and founder SVs underlying breast cancer susceptibility. Our findings demonstrate the potential for nanopore sequencing as a powerful genetic testing assay in the hereditary cancer setting.

MET exon 14 skipping mutation positive non-small cell lung cancer: Response to systemic therapy.

OBJECTIVES: MET exon 14 skipping is a potentially targetable molecular alteration. The goals of this study were to identify patients treated in British Columbia with MET exon 14 skipping to understand prevalence, biology and response to treatment, and to identify molecular signatures that may predict for response or resistance to targeted MET therapy in the setting of advanced disease. MATERIALS AND METHODS: A retrospective review was completed of patients found to have MET exon 14 skipping alterations between January 2016-September 2019. Information was collected on baseline characteristics, response to systemic treatments, and outcomes. RESULTS: Out of 1934 advanced, non-squamous and never-smoking squamous NSCLC patients tested, 41 patients were found to have MET exon 14 skipping (2.1 %). MET alteration types: 2% CBL binding-domain mutations, 34 % poly-pyrimidine tract deletions, 63 % splice donor mutations or deletions. The most common co-mutation was TP53 (22 %). Thirty-three patients received systemic therapy. Physician-assessed disease control was 68 % among 19 evaluable patients treated with crizotinib, 80 % among 10 evaluable patients treated with platinum-based chemotherapy, and 70 % among 10 evaluable patients treated with immunotherapy. Median time to treatment discontinuation was 3.0, 2.8, and 2.4 months, respectively. Median overall survival for metastatic patients treated with any systemic therapy was 15.4 months. In this small cohort, there were no clear correlations between molecular aberrations and response, time to treatment discontinuation, or survival for crizotinib, chemotherapy, and immunotherapy. CONCLUSION: The prevalence of MET exon 14 skipping in a North American population was 2.1 %. Unlike other targetable mutations, patients were older and more commonly current or former smokers. Patients with MET exon 14 skipping alteration demonstrate disease control with crizotinib, platinum-based chemotherapy and immunotherapy. Co-mutations with TP53 were commonly noted, but correlation between co-mutations and efficacy of therapy were not identified in this cohort.

Assessing Limit of Detection in Clinical Sequencing.

Clinical reporting of solid tumor sequencing requires reliable assessment of the accuracy and reproducibility of each assay. Somatic mutation variant allele fractions may be below 10% in many samples due to sample heterogeneity, tumor clonality, and/or sample degradation in fixatives such as formalin. The toolkits available to the clinical sequencing community for correlating assay design parameters with assay sensitivity remain limited, and large-scale empirical assessments are often relied upon due to the lack of clear theoretical grounding. To address this uncertainty, a theoretical model was developed for predicting the expected variant calling sensitivity for a given library complexity and sequencing depth. Binomial models were found to be appropriate when assay sensitivity was only limited by library complexity or sequencing depth, but functional scaling for library complexity was necessary when both library complexity and sequencing depth were co-limiting. This model was empirically validated with sequencing experiments by using a series of DNA input amounts and sequencing depths. Based on these findings, a workflow is proposed for determining the limiting factors to sensitivity in different assay designs, and the formulas for these scenarios are presented. The approach described here provides designers of clinical assays with the methods to theoretically predict assay design outcomes a priori, potentially reducing burden in clinical tumor assay design and validation efforts.

Use of Treatment-Focused Tumor Sequencing to Screen for Germline Cancer Predisposition.

Next-generation sequencing assays are capable of identifying cancer patients eligible for targeted therapies and can also detect germline variants associated with increased cancer susceptibility. However, these capabilities have yet to be routinely harmonized in a single assay because of challenges with accurately identifying germline variants from tumor-only data. We have developed the Oncology and Hereditary Cancer Program targeted capture panel, which uses tumor tissue to simultaneously screen for both clinically actionable solid tumor variants and germline variants across 45 genes. Validation using 14 tumor specimens, composed of patient samples and cell lines analyzed in triplicate, demonstrated high coverage with sensitive and specific identification of single-nucleotide variants and small insertions and deletions. Average coverage across all targets remained >2000× in 198 additional patient tumor samples. Analysis of 55 formalin-fixed, paraffin-embedded tumor samples for the detection of known germline variants within a subset of cancer-predisposition genes, including one multiexon deletion, yielded a 100% detection rate, demonstrating that germline variants can be reliably detected in tumor samples using a single panel. Combining targetable somatic and actionable germline variants into a single tumor tissue assay represents a streamlined approach that can inform treatment for patients with advanced cancers as well as identify those with potential germline variants who are eligible for confirmatory testing, but would not otherwise have been identified.

Integrating Tumor Sequencing Into Clinical Practice for Patients With Mismatch Repair-Deficient Lynch Syndrome Spectrum Cancers.

INTRODUCTION: Uninformative germline genetic testing presents a challenge to clinical management for patients suspected to have Lynch syndrome, a cancer predisposition syndrome caused by germline variants in the mismatch repair (MMR) genes or EPCAM. METHODS: Among a consecutive series of MMR-deficient Lynch syndrome spectrum cancers identified through immunohistochemistry-based tumor screening, we investigated the clinical utility of tumor sequencing for the molecular diagnosis and management of suspected Lynch syndrome families. MLH1-deficient colorectal cancers were prescreened for BRAF V600E before referral for genetic counseling. Microsatellite instability, MLH1 promoter hypermethylation, and somatic and germline genetic variants in the MMR genes were assessed according to an established clinical protocol. RESULTS: Eighty-four individuals with primarily colorectal (62%) and endometrial (31%) cancers received tumor-normal sequencing as part of routine clinical genetic assessment. Overall, 27% received a molecular diagnosis of Lynch syndrome. Most of the MLH1-deficient tumors were more likely of sporadic origin, mediated by MLH1 promoter hypermethylation in 54% and double somatic genetic alterations in MLH1 (17%). MSH2-deficient, MSH6-deficient, and/or PMS2-deficient tumors could be attributed to pathogenic germline variants in 37% and double somatic events in 28%. Notably, tumor sequencing could explain 49% of cases without causal germline variants, somatic MLH1 promoter hypermethylation, or somatic variants in BRAF. DISCUSSION: Our findings support the integration of tumor sequencing into current Lynch syndrome screening programs to improve clinical management for individuals whose germline testing is uninformative.

Sample Tracking Using Unique Sequence Controls.

Sample tracking and identity are essential when processing multiple samples in parallel. Sequencing applications often involve high sample numbers, and the data are frequently used in a clinical setting. As such, a simple and accurate intrinsic sample tracking process through a sequencing pipeline is essential. Various solutions have been implemented to verify sample identity, including variant detection at the start and end of the pipeline using arrays or genotyping, bioinformatic comparisons, and optical barcoding of samples. None of these approaches are optimal. To establish a more effective approach using genetic barcoding, we developed a panel of unique DNA sequences cloned into a common vector. A unique DNA sequence is added to the sample when it is first received and can be detected by PCR and/or sequencing at any stage of the process. The control sequences are approximately 200 bases long with low identity to any sequence in the National Center for Biotechnology Information nonredundant database (<30 bases) and contain no long homopolymer (>7) stretches. When a spiked next-generation sequencing library is sequenced, sequence reads derived from this control sequence are generated along with the standard sequencing run and are used to confirm sample identity and determine cross-contamination levels. This approach is used in our targeted clinical diagnostic whole-genome and RNA-sequencing pipelines and is an inexpensive, flexible, and platform-agnostic solution.

EGFR circulating tumour DNA testing: identification of predictors of ctDNA detection and implications for survival outcomes.

BACKGROUND: EGFR T790M testing is the standard of care for activating EGFR mutation (EGFRm) non-small cell lung cancer (NSCLC) progressing on 1st/2nd generation TKIs to select patients for osimertinib. Despite sensitive assays, detection of circulating tumour deoxyribonucleic acid (ctDNA) is variable and influenced by clinical factors. The number and location of sites of progressive disease at time of testing were reviewed to explore the effect on EGFR ctDNA detection. The prognostic value of EGFR ctDNA detection on survival outcomes was assessed. METHODS: Following extraction of cell-free DNA from plasma using the QIAamp Circulating Nucleic Acid Kit, custom droplet digital polymerase chair reaction (ddPCR) assays were used to assess EGFR ctDNA using the Bio-Rad QX200 system. The ddPCR assay has a limit of detection of ≤0.15% variant allele fraction. Baseline characteristics and imaging reports at time of EGFR ctDNA testing were reviewed retrospectively for a 1 year period. RESULTS: The study included 177 patients who had an EGFR ctDNA test. Liver (aOR 3.13) or bone (aOR 2.76) progression or 3-5 sites of progression (aOR 2.22) were predictive of EGFR ctDNA detection. The median OS from first ctDNA test after multiple testing iterations was 12.3 m undetectable EGFR ctDNA, 7.6 m for original EGFR mutation only and 24.1 m with T790M (P=0.001). CONCLUSIONS: Patients with liver or bone progression and 3-5 progressing sites are more likely to have informative EGFR ctDNA testing. Detection of EGFR ctDNA is a negative prognostic indicator in the absence of a T790M resistance mutation, potentially reflecting the disease burden in the absence of targeted therapy options.

Establishing a Framework for the Clinical Translation of Germline Findings in Precision Oncology.

Inherited genetic variation has important implications for cancer screening, early diagnosis, and disease prognosis. A role for germline variation has also been described in shaping the molecular landscape, immune response, microenvironment, and treatment response of individual tumors. However, there is a lack of consensus on the handling and analysis of germline information that extends beyond known or suspected cancer susceptibility in large-scale cancer genomics initiatives. As part of the Personalized OncoGenomics program in British Columbia, we performed whole-genome and transcriptome sequencing in paired tumor and normal tissues from advanced cancer patients to characterize the molecular tumor landscape and identify putative targets for therapy. Overall, our experience supports a multidisciplinary and integrative approach to germline data management. This includes a need for broader definitions and standardized recommendations regarding primary and secondary germline findings in precision oncology. Here, we propose a framework for identifying, evaluating, and returning germline variants of potential clinical significance that may have indications for health management beyond cancer risk reduction or prevention in patients and their families.

echinus, required for interommatidial cell sorting and cell death in the Drosophila pupal retina, encodes a protein with homology to ubiquitin-specific proteases.

BACKGROUND: Programmed cell death is used to remove excess cells between ommatidia in the Drosophila pupal retina. This death is required to establish the crystalline, hexagonal packing of ommatidia that characterizes the adult fly eye. In previously described echinus mutants, interommatidial cell sorting, which precedes cell death, occurred relatively normally. Interommatidial cell death was partially suppressed, resulting in adult eyes that contained excess pigment cells, and in which ommatidia were mildly disordered. These results have suggested that echinus functions in the pupal retina primarily to promote interommatidial cell death. RESULTS: We generated a number of new echinus alleles, some likely null mutants. Analysis of these alleles provides evidence that echinus has roles in cell sorting as well as cell death. echinus encodes a protein with homology to ubiquitin-specific proteases. These proteins cleave ubiquitin-conjugated proteins at the ubiquitin C-terminus. The echinus locus encodes multiple splice forms, including two proteins that lack residues thought to be critical for deubiquitination activity. Surprisingly, ubiquitous expression in the eye of versions of Echinus that lack residues critical for ubiquitin specific protease activity, as well as a version predicted to be functional, rescue the echinus loss-of-function phenotype. Finally, genetic interactions were not detected between echinus loss and gain-of-function and a number of known apoptotic regulators. These include Notch, EGFR, the caspases Dronc, Drice, Dcp-1, Dream, the caspase activators, Rpr, Hid, and Grim, the caspase inhibitor DIAP1, and Lozenge or Klumpfuss. CONCLUSION: The echinus locus encodes multiple splice forms of a protein with homology to ubiquitin-specific proteases, but protease activity is unlikely to be required for echinus function, at least when echinus is overexpressed. Characterization of likely echinus null alleles and genetic interactions suggests that echinus acts at a novel point(s) to regulate interommatidial cell sorting and/or cell death in the fly eye.

Homologous Recombination Deficiency and Platinum-Based Therapy Outcomes in Advanced Breast Cancer.

Purpose: Recent studies have identified mutation signatures of homologous recombination deficiency (HRD) in over 20% of breast cancers, as well as pancreatic, ovarian, and gastric cancers. There is an urgent need to understand the clinical implications of HRD signatures. Whereas BRCA1/2 mutations confer sensitivity to platinum-based chemotherapies, it is not yet clear whether mutation signatures can independently predict platinum response.Experimental Design: In this observational study, we sequenced tumor whole genomes (100× depth) and matched normals (60×) of 93 advanced-stage breast cancers (33 platinum-treated). We computed a published metric called HRDetect, independently trained to predict BRCA1/2 status, and assessed its capacity to predict outcomes on platinum-based chemotherapies. Clinical endpoints were overall survival (OS), total duration on platinum-based therapy (TDT), and radiographic evidence of clinical improvement (CI).Results: HRDetect predicted BRCA1/2 status with an area under the curve (AUC) of 0.94 and optimal threshold of 0.7. Elevated HRDetect was also significantly associated with CI on platinum-based therapy (AUC = 0.89; P = 0.006) with the same optimal threshold, even after adjusting for BRCA1/2 mutation status and treatment timing. HRDetect scores over 0.7 were associated with a 3-month extended median TDT (P = 0.0003) and 1.3-year extended median OS (P = 0.04).Conclusions: Our findings not only independently validate HRDetect, but also provide the first evidence of its association with platinum response in advanced breast cancer. We demonstrate that HRD mutation signatures may offer clinically relevant information independently of BRCA1/2 mutation status and hope this work will guide the development of clinical trials. Clin Cancer Res; 23(24); 7521-30. ©2017 AACR.

A set of BAC clones spanning the human genome.

Using the human bacterial artificial chromosome (BAC) fingerprint-based physical map, genome sequence assembly and BAC end sequences, we have generated a fingerprint-validated set of 32 855 BAC clones spanning the human genome. The clone set provides coverage for at least 98% of the human fingerprint map, 99% of the current assembled sequence and has an effective resolving power of 79 kb. We have made the clone set publicly available, anticipating that it will generally facilitate FISH or array-CGH-based identification and characterization of chromosomal alterations relevant to disease.

A population-based review of the feasibility of platinum-based combination chemotherapy after tyrosine kinase inhibition in EGFR mutation positive non-small cell lung cancer patients with advanced disease.

INTRODUCTION: The IPASS trial demonstrated superior progression free survival for Asian, light/never smoking, advanced, pulmonary adenocarcinoma patients treated with first-line gefitinib compared to carboplatin/paclitaxel, of which 59% of those tested were epidermal growth factor receptor (EGFR) mutation positive. In IPASS 39% of gefitinib treated patients went on to receive platin based polychemotherapy. We hypothesized that in a population-based setting fewer patients receive second-line platin based chemotherapy than those enrolled in a clinical trial. METHODS: The Iressa Alliance program provided standardized EGFR mutation testing and appropriate access to gefitinib to all patients in British Columbia with advanced, non squamous non small cell lung cancer (NSCLC). We retrospectively analyzed clinical, pathologic data and outcomes for all patients tested in this program between March 2010 and June 2011. RESULTS: A total of 548 patients were referred for testing and 22% of patients were mutation positive. Baseline characteristics of mutation negative and mutation positive; median age 67/65, male 41%/31%, Asian 15%/51%, never smoker 21%/58%, stage IV 80%/91%. Median overall survival was 12 months in mutation negative patients and not yet reached in mutation positive (p<0.0001). In mutation positive patients 5% of patients had a complete response, 46% partial response, 34% stable disease, 6% progressive disease. Twenty percent of patients continued on gefitinib after radiographic progression and clinical stability. Sixty-one gefitinib treated patients progressed at the time of analysis; 10% of patients received further gefitinib only, 38% platinum based doublet, 8% other chemotherapy and 44% no further treatment. Performance status most strongly predicted for delivery of second line chemotherapy. CONCLUSIONS: This North American population based study shows similar efficacy of gefitinib in mutation positive patients compared to the IPASS trial. Contrary to our hypothesis, delivery of second line chemotherapy was feasible in a significant proportion of gefitinib treated patients.

A clinically validated diagnostic second-generation sequencing assay for detection of hereditary BRCA1 and BRCA2 mutations.

Individuals who inherit mutations in BRCA1 or BRCA2 are predisposed to breast and ovarian cancers. However, identifying mutations in these large genes by conventional dideoxy sequencing in a clinical testing laboratory is both time consuming and costly, and similar challenges exist for other large genes, or sets of genes, with relevance in the clinical setting. Second-generation sequencing technologies have the potential to improve the efficiency and throughput of clinical diagnostic sequencing, once clinically validated methods become available. We have developed a method for detection of variants based on automated small-amplicon PCR followed by sample pooling and sequencing with a second-generation instrument. To demonstrate the suitability of this method for clinical diagnostic sequencing, we analyzed the coding exons and the intron-exon boundaries of BRCA1 and BRCA2 in 91 hereditary breast cancer patient samples. Our method generated high-quality sequence coverage across all targeted regions, with median coverage greater than 4000-fold for each sample in pools of 24. Sensitive and specific automated variant detection, without false-positive or false-negative results, was accomplished with a standard software pipeline using bwa for sequence alignment and samtools for variant detection. We experimentally derived a minimum threshold of 100-fold sequence depth for confident variant detection. The results demonstrate that this method is suitable for sensitive, automatable, high-throughput sequence variant detection in the clinical laboratory.

A BAC clone fingerprinting approach to the detection of human genome rearrangements.

We present a method, called fingerprint profiling (FPP), that uses restriction digest fingerprints of bacterial artificial chromosome clones to detect and classify rearrangements in the human genome. The approach uses alignment of experimental fingerprint patterns to in silico digests of the sequence assembly and is capable of detecting micro-deletions (1-5 kb) and balanced rearrangements. Our method has compelling potential for use as a whole-genome method for the identification and characterization of human genome rearrangements.