Pisces: an accurate and versatile variant caller for somatic and germline next-generation sequencing data

FRα expression in advanced endometrial carcinoma: Clinicopathological, molecular, and prognosis correlates

BACKGROUND

Advanced endometrial carcinoma (EC) is a therapeutic challenge with limited treatment options. Folate receptor alpha (FRα)-targeted antibody-drug conjugates (ADCs) are being explored for EC treatment. However, the molecular profiles of FRα-expressing ECs remain understudied. This study aims to evaluate FRα expression in advanced EC and its associations with clinicopathological features, molecular alterations, and prognosis value.

METHODS

We analyzed a cohort of 111 advanced ECs from our institution (PUMCH). We used immunohistochemistry (IHC) to quantify FRα expression and next-generation sequencing (NGS) to determine genomic mutations. In parallel, a cohort from the TCGA database was included as an external cohort for validation.

RESULTS

In the PUMCH cohort, 21 out of 111 (18.9 %) cases were graded as FRα-high. FRα-high ECs were significantly associated with aggressive histotypes and the p53abn molecular subtype. A significant correlation between FRα and CA125 was demonstrated at the protein level in the PUMCH cohort and the RNA level in the TCGA cohort. Molecularly, FRα-high ECs were more often associated with TP53 mutations, while FRα-low ECs were characterized by PTEN mutations. In survival analysis, high FRα expression was closely correlated with poor outcomes, especially in Stage III ECs. Additionally, FRα could significantly stratify prognosis within a limited cohort of mismatch repair-deficient (MMR-d) ECs.

CONCLUSION

FRα positivity in advanced ECs corresponds with a unique molecular landscape, particularly the p53abn subgroup. The prognostic value of FRα in advanced ECs, especially in MMR-d ECs, warrants clinical attention.

Analysis of DNA from brain tissue on stereo-EEG electrodes reveals mosaic epilepsy-related variants

Somatic mosaic variants contribute to focal epilepsy, with variants often present only in brain tissue and not in blood or other samples typically assayed for genetic testing. Thus, genetic analysis for mosaic variants in focal epilepsy has been limited to patients with drug-resistant epilepsy who undergo surgical resection and have resected brain tissue samples available. Stereo-EEG (sEEG) has become part of the evaluation for many patients with focal drug-resistant epilepsy, and sEEG electrodes provide a potential source of small amounts of brain-derived DNA. We aimed to identify, validate, and assess the distribution of deleterious mosaic variants in epilepsy-associated genes in DNA extracted from trace brain tissue on individual sEEG electrodes. We enrolled a prospective cohort of 10 paediatric patients with drug-resistant epilepsy who had sEEG electrodes implanted for invasive monitoring. We extracted unamplified DNA and in parallel performed whole-genome amplification from trace brain tissue on each sEEG electrode. We also extracted DNA from resected brain tissue and blood/saliva samples where available. We performed deep sequencing (panel and exome) and analysis for candidate germline and mosaic variants. We validated candidate mosaic variants and assessed the variant allele fraction in amplified and unamplified electrode-derived DNA and across electrodes. We extracted unamplified DNA and performed whole-genome amplification from >150 individual electrodes from 10 individuals. Immunohistochemistry confirmed the presence of neurons in the brain tissue on electrodes. Deep sequencing and analysis demonstrated similar depth of coverage between amplified and unamplified DNA samples but significantly more potential mosaic variants in amplified samples. We validated four deleterious mosaic variants in epilepsy-associated genes in electrode-derived DNA in three patients who underwent laser ablation and did not have resected brain tissue samples available. Three of the four variants were detected in both amplified and unamplified electrode-derived DNA, with higher variant allele fraction observed in DNA from electrodes in closest proximity to the electrical seizure focus in one case. We demonstrate that mosaic variants can be identified and validated from DNA extracted from trace brain tissue on individual sEEG electrodes in patients with drug-resistant focal epilepsy, from both unamplified and amplified electrode-derived DNA. Our findings support a relationship between the extent of regional genetic abnormality and electrophysiology and suggest that with further optimization, this minimally invasive diagnostic approach holds promise for advancing precision medicine for patients with drug-resistant epilepsy as part of the surgical evaluation.

Enhanced Error Suppression for Accurate Detection of Low-Frequency Variants

The identification of low-frequency variants remains challenging due to the inevitable high error rates of next-generation sequencing (NGS). Numerous promising strategies employ unique molecular identifiers (UMIs) for error suppression. However, their efficiency depends highly on redundant sequencing and quality control, leading to tremendous read waste and cost inefficiency. Here, we describe a novel approach, enhanced error suppression strategy (EES), that addresses these challenges by (1) optimizing data utilization and reducing read waste by utilizing single-read correction that reserves abundant single reads that complement other single reads or single-strand consensus sequences (SSCSs), and (2) effectively enhancing the accuracy of NGS by employing Bayes' theorem. EES significantly improves variant detection accuracy, achieving a background error rate of less than 4.4 × 10-5 per base pair. Additionally, the data utilization rate is dramatically increased, with a 22.9-fold enhancement in duplex consensus sequence (DCS) recovery compared to traditional methodologies. Furthermore, EES demonstrates superior error suppression performance across various base substitutions. In conclusion, EES represents a significant advancement in detecting low-frequency variants by improving data utilization and reducing sequencing errors. It potentially enhances the sensitivity and accuracy of NGS applications, proving highly valuable in clinical and research contexts where precise variant detection is critical.

Pediatric therapy-related hematologic neoplasms show enrichment for KMT2A rearrangement and lymphoblastic phenotype

In children, therapy-related hematologic neoplasms (t-HN) are uncommon. Many are driven by genetic events independent of clonal hematopoiesis. We sought to understand the clinical and genetic factors of pediatric t-HN in a large independent cohort. Fifty-six t-HN were retrospectively identified. Chromosome microarray, next-generation and/or RNA sequencing were performed. Patients had primary hematologic, solid, or central nervous system tumors. t-HN included myeloid (t-MN) and lymphoblastic (t-ALL) phenotypes. Approximately half of the cases harbored KMTA2A rearrangement (KMT2Ar). Among t-HN without KMT2Ar, genetic drivers were heterogeneous, including diverse fusions or aneuploidy. Approximately 18% harbored 17p deletions and/or TP53 mutations. EFS/OS was not associated with t-HN lineage or KMT2Ar, but HSCT was associated with improved EFS and OS. We detail one of the largest cohorts to date of pediatric t-HN, confirming frequent KMT2Ar and t-ALL.

Development of a Tagmentation-Based Next-Generation Sequencing Clinical Assay as an Alternative to Capillary Electrophoresis-Based Sequencing

BACKGROUND

Next-generation sequencing (NGS) technology enables sample multiplexing for interrogation of multiple regions of interest (ROI). Leveraging this, together with access to affordable NGS platforms, we explored the practicality of moving capillary electrophoresis (CE), noncapillary electrophoresis and single-gene testing to NGS. In this work, we evaluated the iSeq 100's capacity to validate 89 samples at once.

METHODS

Genomic DNA was extracted from 89 archival samples of varying specimen types. Polymerase chain reaction (PCR) was done with in house primers, library preparation with the Nextera XT Library Preparation Kit and cleaning up with paramagnetic beads. The sequencing was performed on one Illumina iSeq 100 flow cell.

RESULTS

With our workflow, 88 out of 89 samples were accurately sequenced with variant alleles identified. One sample of the 88 samples was initially discordant because the primers used were in a heterozygous deletion region. Upon redesigning of primers, the sample proved concordant.

CONCLUSIONS

The iSeq-Nextera workflow proved accurate. However, variant allele frequencys generated by the Nextera are not precise.

Comparison of Results from Two Commercially Available In-House Tissue-Based Comprehensive Genomic Profiling Solutions: Research Use Only AVENIO Tumor Tissue Comprehensive Genomic Profiling Kit and TruSight Oncology 500 Assay

Increased adoption of personalized medicine has brought comprehensive genomic profiling (CGP) to the forefront. However, differences in assay, bioinformatics, and reporting systems and lack of understanding of their complex interplay are a challenge for implementation and achieving uniformity in CGP testing. Two commercially available, tissue-based, in-house CGP assays were compared, in combination with a tertiary analysis solution in a research use only (RUO) context: the AVENIO Tumor Tissue CGP RUO Kit paired with navify Mutation Profiler (RUO) software and the TruSight Oncology 500 RUO assay paired with PierianDx Clinical Genomics Workspace software. Agreements and differences between the assays were assessed for short variants, copy number alterations, rearrangements, tumor mutational burden, and microsatellite instability, including variant categorization and clinical trial-matching (CTM) recommendations. Results showed good overall agreement for short variant, known gene fusion, and microsatellite instability detection. Important differences were obtained in tumor mutational burden scoring, copy number alteration detection, and CTM. Differences in variant and biomarker detection could be explained by bioinformatic approaches to variant calling, filtering, tiering, and normalization; differences in CTM, by underlying reported variants and conceptual differences in system parameters. Thus, distinctions between different approaches may lead to inconsistent results. Complexities in calling, filtering, and interpreting variants illustrate key considerations for implementation of any high-quality CGP in the laboratory and bringing uniformity to genomic insight results.

Association Between Clonal Hematopoiesis and Left Ventricular Reverse Remodeling in Nonischemic Dilated Cardiomyopathy

Although clonal hematopoiesis of indeterminate potential (CHIP) is an adverse prognostic factor for atherosclerotic disease, its impact on nonischemic dilated cardiomyopathy (DCM) is elusive. The authors performed whole-exome sequencing and deep target sequencing among 198 patients with DCM and detected germline mutations in cardiomyopathy-related genes and somatic mutations in CHIP driver genes. Twenty-five CHIP driver mutations were detected in 22 patients with DCM. Ninety-two patients had cardiomyopathy-related pathogenic mutations. Multivariable analysis revealed that CHIP was an independent risk factor of left ventricular reverse remodeling, irrespective of known prognostic factors. CHIP exacerbated cardiac systolic dysfunction and fibrosis in a DCM murine model. The identification of germline and somatic mutations in patients with DCM predicts clinical prognosis.

Analysis of DNA from brain tissue on stereo-EEG electrodes reveals mosaic epilepsy-related variants

Somatic mosaic variants contribute to focal epilepsy, but genetic analysis has been limited to patients with drug-resistant epilepsy (DRE) who undergo surgical resection, as the variants are mainly brain-limited. Stereoelectroencephalography (sEEG) has become part of the evaluation for many patients with focal DRE, and sEEG electrodes provide a potential source of small amounts of brain-derived DNA. We aimed to identify, validate, and assess the distribution of potentially clinically relevant mosaic variants in DNA extracted from trace brain tissue on individual sEEG electrodes. We enrolled a prospective cohort of eleven pediatric patients with DRE who had sEEG electrodes implanted for invasive monitoring, one of whom was previously reported. We extracted unamplified DNA from the trace brain tissue on each sEEG electrode and also performed whole-genome amplification for each sample. We extracted DNA from resected brain tissue and blood/saliva samples where available. We performed deep panel and exome sequencing on a subset of samples from each case and analysis for potentially clinically relevant candidate germline and mosaic variants. We validated candidate mosaic variants using amplicon sequencing and assessed the variant allele fraction (VAF) in amplified and unamplified electrode-derived DNA and across electrodes. We extracted DNA from >150 individual electrodes from 11 individuals and obtained higher concentrations of whole-genome amplified vs unamplified DNA. Immunohistochemistry confirmed the presence of neurons in the brain tissue on electrodes. Deep sequencing and analysis demonstrated similar depth of coverage between amplified and unamplified samples but significantly more called mosaic variants in amplified samples. In addition to the mosaic PIK3CA variant detected in a previously reported case from our group, we identified and validated four potentially clinically relevant mosaic variants in electrode-derived DNA in three patients who underwent laser ablation and did not have resected brain tissue samples available. The variants were detected in both amplified and unamplified electrode-derived DNA, with higher VAFs observed in DNA from electrodes in closest proximity to the electrical seizure focus in some cases. This study demonstrates that mosaic variants can be identified and validated from DNA extracted from trace brain tissue on individual sEEG electrodes in patients with drug-resistant focal epilepsy and in both amplified and unamplified electrode-derived DNA samples. Our findings support a relationship between the extent of regional genetic abnormality and electrophysiology, and suggest that with further optimization, this minimally invasive diagnostic approach holds promise for advancing precision medicine for patients with DRE as part of the surgical evaluation.

Integrative Single Cell Atlas Revealed Intratumoral Heterogeneity Generation from an Adaptive Epigenetic Cell State in Human Bladder Urothelial Carcinoma

Intratumor heterogeneity (ITH) of bladder cancer (BLCA) contributes to therapy resistance and immune evasion affecting clinical prognosis. The molecular and cellular mechanisms contributing to BLCA ITH generation remain elusive. It is found that a TM4SF1-positive cancer subpopulation (TPCS) can generate ITH in BLCA, evidenced by integrative single cell atlas analysis. Extensive profiling of the epigenome and transcriptome of all stages of BLCA revealed their evolutionary trajectories. Distinct ancestor cells gave rise to low-grade noninvasive and high-grade invasive BLCA. Epigenome reprograming led to transcriptional heterogeneity in BLCA. During early oncogenesis, epithelial-to-mesenchymal transition generated TPCS. TPCS has stem-cell-like properties and exhibited transcriptional plasticity, priming the development of transcriptionally heterogeneous descendent cell lineages. Moreover, TPCS prevalence in tumor is associated with advanced stage cancer and poor prognosis. The results of this study suggested that bladder cancer interacts with its environment by acquiring a stem cell-like epigenomic landscape, which might generate ITH without additional genetic diversification.

Lysosomal genes contribute to Parkinson's disease near agriculture with high intensity pesticide use

Parkinson's disease (PD), the second most common neurodegenerative disorder, develops sporadically, likely through a combination of polygenic and environmental factors. Previous studies associate pesticide exposure and genes involved in lysosomal function with PD risk. We evaluated the frequency of variants in lysosomal function genes among patients from the Parkinson's, Environment, and Genes (PEG) study with ambient pesticide exposure from agricultural sources. 757 PD patients, primarily of White European/non-Hispanic ancestry (75%), were screened for variants in 85 genes using a custom amplicon panel. Variant enrichment was calculated against the Genome Aggregation Database (gnomAD). Enriched exonic variants were prioritized by exposure to a cluster of pesticides used on cotton and severity of disease progression in a subset of 386 patients subdivided by race/ethnicity. Gene enrichment analysis identified 36 variants in 26 genes in PEG PD patients. Twelve of the identified genes (12/26, 46%) had multiple enriched variants and/or a single enriched variant present in multiple individuals, representing 61% (22/36) of the observed variation in the cohort. The majority of enriched variants (26/36, 72%) were found in genes contributing to lysosomal function, particularly autophagy, and were bioinformatically deemed functionally deleterious (31/36, 86%). We conclude that, in this study, variants in genes associated with lysosomal function, notably autophagy, were enriched in PD patients exposed to agricultural pesticides suggesting that altered lysosomal function may generate an underlying susceptibility for developing PD with pesticide exposure. Further study of gene-environment interactions targeting lysosomal function may improve understanding of PD risk in individuals exposed to pesticides.

ASCL1 Drives Tolerance to Osimertinib in EGFR Mutant Lung Cancer in Permissive Cellular Contexts

The majority of EGFR mutant lung adenocarcinomas respond well to EGFR tyrosine kinase inhibitors (TKI). However, most of these responses are partial, with drug-tolerant residual disease remaining even at the time of maximal response. This residual disease can ultimately lead to relapses, which eventually develop in most patients. To investigate the cellular and molecular properties of residual tumor cells in vivo, we leveraged patient-derived xenograft (PDX) models of EGFR mutant lung cancer. Subcutaneous EGFR mutant PDXs were treated with the third-generation TKI osimertinib until maximal tumor regression. Residual tissue inevitably harbored tumor cells that were transcriptionally distinct from bulk pretreatment tumor. Single-cell transcriptional profiling provided evidence of cells matching the profiles of drug-tolerant cells present in the pretreatment tumor. In one of the PDXs analyzed, osimertinib treatment caused dramatic transcriptomic changes that featured upregulation of the neuroendocrine lineage transcription factor ASCL1. Mechanistically, ASCL1 conferred drug tolerance by initiating an epithelial-to-mesenchymal gene-expression program in permissive cellular contexts. This study reveals fundamental insights into the biology of drug tolerance, the plasticity of cells through TKI treatment, and why specific phenotypes are observed only in certain tumors.

SIGNIFICANCE

Analysis of residual disease following tyrosine kinase inhibitor treatment identified heterogeneous and context-specific mechanisms of drug tolerance in lung cancer that could lead to the development of strategies to forestall drug resistance. See related commentary by Rumde and Burns, p. 1188.

Clinical, pathologic and molecular findings in 2 Rottweiler littermates with appendicular osteosarcoma

Appendicular osteosarcoma was diagnosed and treated in a pair of littermate Rottweiler dogs, resulting in distinctly different clinical outcomes despite similar therapy within the context of a prospective, randomized clinical trial (NCI-COTC021/022). Histopathology, immunohistochemistry, mRNA sequencing, and targeted DNA hotspot sequencing techniques were applied to both dogs' tumors to define factors that could underpin their differential response to treatment. We describe the comparison of their clinical, histologic and molecular features, as well as those from a companion cohort of Rottweiler dogs, providing new insight into potential prognostic biomarkers for canine osteosarcoma.

Integrative genomic analysis identifies unique immune environments associated with immunotherapy response in diffuse large B cell lymphoma

Most diffuse large B-cell lymphoma (DLBCL) patients treated with bispecific antibodies (BsAb) or chimeric antigen receptor (CAR) T cells fail to achieve durable treatment responses, underscoring the need for a deeper understanding of mechanisms that regulate the immune environment and response to treatment. Here, an integrative, multi-omic approach was employed to characterize DLBCL immune environments, which effectively segregated DLBCLs into four quadrants - termed DLBCL-immune quadrants (IQ) - defined by cell-of-origin and immune-related gene set expression scores. Recurrent genomic alterations were enriched in each IQ, suggesting that lymphoma cell-intrinsic alterations contribute to orchestrating unique DLBCL immune environments. In relapsed/refractory DLBCL patients, DLBCL-IQ assignment correlated significantly with clinical benefit with the CD20 x CD3 BsAb, mosunetuzumab, but not with CD19-directed CAR T cells. DLBCL-IQ provides a new framework to conceptualize the DLBCL immune landscape and uncovers the differential impact of the endogenous immune environment on outcomes to BsAb and CAR T cell treatment.

Comparing malignant monocytosis across the updated WHO and ICC classifications of 2022

ABSTRACT

The World Health Organization (WHO) classification of hematolymphoid tumors and the International Consensus Classification (ICC) of 2022 introduced major changes to the definition of chronic myelomonocytic leukemia (CMML). To assess its qualitative and quantitative implications for patient care, we started with 3311 established CMML cases (according to WHO 2017 criteria) and included 2130 oligomonocytosis cases fulfilling the new CMML diagnostic criteria. Applying both 2022 classification systems, 356 and 241 of oligomonocytosis cases were newly classified as myelodysplastic (MD)-CMML (WHO and ICC 2022, respectively), most of which were diagnosed as myelodysplastic syndrome (MDS) according to the WHO 2017 classification. Importantly, 1.5 times more oligomonocytosis cases were classified as CMML according to WHO 2022 than based on ICC, because of different diagnostic criteria. Genetic analyses of the newly classified CMML cases showed a distinct mutational profile with strong enrichment of MDS-typical alterations, resulting in a transcriptional subgroup separated from established MD and myeloproliferative CMML. Despite a different cytogenetic, molecular, immunophenotypic, and transcriptional landscape, no differences in overall survival were found between newly classified and established MD-CMML cases. To the best of our knowledge, this study represents the most comprehensive analysis of routine CMML cases to date, both in terms of clinical characterization and transcriptomic analysis, placing newly classified CMML cases on a disease continuum between MDS and previously established CMML.

Neuropathological Applications of Microscopy with Ultraviolet Surface Excitation (MUSE): A Concordance Study of Human Primary and Metastatic Brain Tumors

Whereas traditional histology and light microscopy require multiple steps of formalin fixation, paraffin embedding, and sectioning to generate images for pathologic diagnosis, Microscopy using Ultraviolet Surface Excitation (MUSE) operates through UV excitation on the cut surface of tissue, generating images of high resolution without the need to fix or section tissue and allowing for potential use for downstream molecular tests. Here, we present the first study of the use and suitability of MUSE microscopy for neuropathological samples. MUSE images were generated from surgical biopsy samples of primary and metastatic brain tumor biopsy samples (n = 27), and blinded assessments of diagnoses, tumor grades, and cellular features were compared to corresponding hematoxylin and eosin (H&E) images. A set of MUSE-treated samples subsequently underwent exome and targeted sequencing, and quality metrics were compared to those from fresh frozen specimens. Diagnostic accuracy was relatively high, and DNA and RNA integrity appeared to be preserved for this cohort. This suggests that MUSE may be a reliable method of generating high-quality diagnostic-grade histologic images for neuropathology on a rapid and sample-sparing basis and for subsequent molecular analysis of DNA and RNA.

Somatic cancer driver mutations are enriched and associated with inflammatory states in Alzheimer's disease microglia

Alzheimer's disease (AD) is an age-associated neurodegenerative disorder characterized by progressive neuronal loss and pathological accumulation of the misfolded proteins amyloid-β and tau1,2. Neuroinflammation mediated by microglia and brain-resident macrophages plays a crucial role in AD pathogenesis1-5, though the mechanisms by which age, genes, and other risk factors interact remain largely unknown. Somatic mutations accumulate with age and lead to clonal expansion of many cell types, contributing to cancer and many non-cancer diseases6,7. Here we studied somatic mutation in normal aged and AD brains by three orthogonal methods and in three independent AD cohorts. Analysis of bulk RNA sequencing data from 866 samples from different brain regions revealed significantly higher (~two-fold) overall burdens of somatic single-nucleotide variants (sSNVs) in AD brains compared to age-matched controls. Molecular-barcoded deep (>1000X) gene panel sequencing of 311 prefrontal cortex samples showed enrichment of sSNVs and somatic insertions and deletions (sIndels) in cancer driver genes in AD brain compared to control, with recurrent, and often multiple, mutations in genes implicated in clonal hematopoiesis (CH)8,9. Pathogenic sSNVs were enriched in CSF1R+ microglia of AD brains, and the high proportion of microglia (up to 40%) carrying some sSNVs in cancer driver genes suggests mutation-driven microglial clonal expansion (MiCE). Analysis of single-nucleus RNA sequencing (snRNAseq) from temporal neocortex of 62 additional AD cases and controls exhibited nominally increased mosaic chromosomal alterations (mCAs) associated with CH10,11. Microglia carrying mCA showed upregulated pro-inflammatory genes, resembling the transcriptomic features of disease-associated microglia (DAM) in AD. Our results suggest that somatic driver mutations in microglia are common with normal aging but further enriched in AD brain, driving MiCE with inflammatory and DAM signatures. Our findings provide the first insights into microglial clonal dynamics in AD and identify potential new approaches to AD diagnosis and therapy.

The clinical and genomic landscape of patients with DDX41 variants identified during diagnostic sequencing

Deleterious germ line variants in DDX41 are a common cause of genetic predisposition to hematologic malignancies, particularly myelodysplastic neoplasms (MDS) and acute myeloid leukemia (AML). Targeted next-generation sequencing was performed in a large cohort of sequentially recruited patients with myeloid malignancy, covering DDX41 as well as 30 other genes frequently mutated in myeloid malignancy. Whole genome transcriptome sequencing data was analyzed on a separate cohort of patients with a range of hematologic malignancies to investigate the spectrum of cancer predisposition. Altogether, 5737 patients with myeloid malignancies were studied, with 152 different DDX41 variants detected. Multiple novel variants were detected, including synonymous variants affecting splicing as demonstrated by RNA-sequencing. The presence of a somatic DDX41 variant was highly associated with DDX41 germ line variants in patients with MDS and AML, and we developed a statistical approach to incorporate the co-occurrence of a somatic DDX41 variant into germ line variant classification at a very strong level (as per the American College of Medical Genetics and Genomics/Association for Molecular Pathology guidelines). Using this approach, the MDS cohort contained 108 of 2865 (3.8%) patients with germ line likely pathogenic/pathogenic (LP/P) variants, and the AML cohort 106 of 2157 (4.9%). DDX41 LP/P variants were markedly enriched in patients with AML and MDS compared with those in patients with myeloproliferative neoplasms, B-cell neoplasm, and T- or B-cell acute lymphoblastic leukemia. In summary, we have developed a framework to enhance DDX41 variant curation as well as highlighted the importance of assessment of all types of genomic variants (including synonymous and multiexon deletions) to fully detect the landscape of possible clinically relevant DDX41 variants.

Association of ESR1 germline variants with TP53 somatic variants in breast tumors in a genome-wide study

Background

In breast tumors, somatic mutation frequencies in TP53 and PIK3CA vary by tumor subtype and ancestry. HER2 positive and triple negative breast cancers (TNBC) have a higher frequency of TP53 somatic mutations than other subtypes. PIK3CA mutations are more frequently observed in hormone receptor positive tumors. Emerging data suggest tumor mutation status is associated with germline variants and genetic ancestry. We aimed to identify germline variants that are associated with somatic TP53 or PIK3CA mutation status in breast tumors.

Methods

A genome-wide association study was conducted using breast cancer mutation status of TP53 and PIK3CA and functional mutation categories including TP53 gain of function (GOF) and loss of function mutations and PIK3CA activating/hotspot mutations. The discovery analysis consisted of 2850 European ancestry women from three datasets. Germline variants showing evidence of association with somatic mutations were selected for validation analyses based on predicted function, allele frequency, and proximity to known cancer genes or risk loci. Candidate variants were assessed for association with mutation status in a multi-ancestry validation study, a Malaysian study, and a study of African American/Black women with TNBC.

Results

The discovery Germline x Mutation (GxM) association study found five variants associated with one or more TP53 phenotypes with P values <1×10-6, 33 variants associated with one or more TP53 phenotypes with P values <1×10-5, and 44 variants associated with one or more PIK3CA phenotypes with P values <1×10-5. In the multi-ancestry and Malaysian validation studies, germline ESR1 locus variant, rs9383938, was associated with the presence of TP53 mutations overall (P values 6.8×10-5 and 9.8×10-8, respectively) and TP53 GOF mutations (P value 8.4×10-6). Multiple variants showed suggestive evidence of association with PIK3CA mutation status in the validation studies, but none were significant after correction for multiple comparisons.

Conclusions

We found evidence that germline variants were associated with TP53 and PIK3CA mutation status in breast cancers. Variants near the estrogen receptor alpha gene, ESR1, were significantly associated with overall TP53 mutations and GOF mutations. Larger multi-ancestry studies are needed to confirm these findings and determine if these variants contribute to ancestry-specific differences in mutation frequency.

Somatic Mosaicism in Amyotrophic Lateral Sclerosis and Frontotemporal Dementia Reveals Widespread Degeneration from Focal Mutations

Although mutations in dozens of genes have been implicated in familial forms of amyotrophic lateral sclerosis (fALS) and frontotemporal degeneration (fFTD), most cases of these conditions are sporadic (sALS and sFTD), with no family history, and their etiology remains obscure. We tested the hypothesis that somatic mosaic mutations, present in some but not all cells, might contribute in these cases, by performing ultra-deep, targeted sequencing of 88 genes associated with neurodegenerative diseases in postmortem brain and spinal cord samples from 404 individuals with sALS or sFTD and 144 controls. Known pathogenic germline mutations were found in 20.6% of ALS, and 26.5% of FTD cases. Predicted pathogenic somatic mutations in ALS/FTD genes were observed in 2.7% of sALS and sFTD cases that did not carry known pathogenic or novel germline mutations. Somatic mutations showed low variant allele fraction (typically <2%) and were often restricted to the region of initial discovery, preventing detection through genetic screening in peripheral tissues. Damaging somatic mutations were preferentially enriched in primary motor cortex of sALS and prefrontal cortex of sFTD, mirroring regions most severely affected in each disease. Somatic mutation analysis of bulk RNA-seq data from brain and spinal cord from an additional 143 sALS cases and 23 controls confirmed an overall enrichment of somatic mutations in sALS. Two adult sALS cases were identified bearing pathogenic somatic mutations in DYNC1H1 and LMNA, two genes associated with pediatric motor neuron degeneration. Our study suggests that somatic mutations in fALS/fFTD genes, and in genes associated with more severe diseases in the germline state, contribute to sALS and sFTD, and that mosaic mutations in a small fraction of cells in focal regions of the nervous system can ultimately result in widespread degeneration.

Direct prediction of carbapenem resistance in Pseudomonas aeruginosa by whole genome sequencing and metagenomic sequencing

Carbapenem resistance is a major concern in the management of antibiotic-resistant Pseudomonas aeruginosa infections. The direct prediction of carbapenem-resistant phenotype from genotype in P. aeruginosa isolates and clinical samples would promote timely antibiotic therapy. The complex carbapenem resistance mechanism and the high prevalence of variant-driven carbapenem resistance in P. aeruginosa make it challenging to predict the carbapenem-resistant phenotype through the genotype. In this study, using whole genome sequencing (WGS) data of 1,622 P. aeruginosa isolates followed by machine learning, we screened 16 and 31 key gene features associated with imipenem (IPM) and meropenem (MEM) resistance in P. aeruginosa, including oprD(HIGH), and constructed the resistance prediction models. The areas under the curves of the IPM and MEM resistance prediction models were 0.906 and 0.925, respectively. For the direct prediction of carbapenem resistance in P. aeruginosa from clinical samples by the key gene features selected and prediction models constructed, 72 P. aeruginosa-positive sputum samples were collected and sequenced by metagenomic sequencing (MGS) based on next-generation sequencing (NGS) or Oxford Nanopore Technology (ONT). The prediction applicability of MGS based on NGS outperformed that of MGS based on ONT. In 72 P. aeruginosa-positive sputum samples, 65.0% (26/40) of IPM-insensitive and 63.2% (24/38) of MEM-insensitive P. aeruginosa were directly predicted by NGS-based MGS with positive predictive values of 0.897 and 0.889, respectively. By the direct detection of the key gene features associated with carbapenem resistance of P. aeruginosa, the carbapenem resistance of P. aeruginosa could be directly predicted from cultured isolates by WGS or from clinical samples by NGS-based MGS, which could assist the timely treatment and surveillance of carbapenem-resistant P. aeruginosa.

Malignant Histiocytosis Comprises a Phenotypic Spectrum That Parallels the Lineage Differentiation of Monocytes, Macrophages, Dendritic Cells, and Langerhans Cells

Malignant histiocytoses (MHs), or the 'M group' of the Histiocyte Society classification, are characterized by neoplastic histiocytes with large pleomorphic nuclei. MH encompasses the diagnoses of histiocytic sarcoma, interdigitating dendritic cell sarcoma, and Langerhans cell sarcoma. We aimed to define the phenotypic spectrum of MH and examine the genotypic features across this spectrum. Using immunohistochemistry, we arranged the 22 cases into 4 subtypes that correspond to the lines of differentiation from monocytic and dendritic cell precursors as follows: (1) macrophage (n = 5): CD68+, CD163+, CD14+, and Factor 13a+; (2) monocyte-macrophage (n = 5): CD68+, CD163+, CD14+, S100+, and OCT2+; (3) dendritic cell (n = 6): CD68+, CD11c+, S100+, lysozyme+, ZBTB46+, and CD1a/langerin < 5%; and (4) Langerhans cell (n = 6): CD68+, CD11c+, S100+, ZBTB46+, CD1a+, and langerin+. The phenotypic subtypes align with those seen in low-grade histiocytic neoplasms as follows: MH-macrophage type correlates with Erdheim-Chester disease phenotype; MH-monocyte-macrophage type with Rosai-Dorfman disease phenotype, and MH-Langerhans cell type with Langerhans cell histiocytosis. Activating mutations in MAPK-pathway genes were identified in 80% of MH cases; 29% had mutations in the PI3k-AKT-mTOR pathway and 59% had mutations in epigenetic modulating genes. Strong expression of cyclin D1 was present in all cases, whereas p-ERK and p-AKT were not uniformly expressed. Eight of 22 (36%) MH cases were proven to be clonally related to a prior B-cell lymphoma. Defining the phenotypic spectrum of MH provides a guide to diagnosis and allows further exploration into the potential biological and clinical significance.

Identification of New Potential Prognostic and Predictive Markers in High-Grade Osteosarcoma Using Whole Exome Sequencing

Conventional high-grade osteosarcoma (OS) is the most common primary cancer of bone and it typically affects the extremities of adolescents. OS has a complex karyotype, and molecular mechanisms related to carcinogenesis, progression and resistance to therapy are still largely unknown. For this reason, the current standard of care is associated with considerable adverse effects. In this study, our aim was to identify gene alterations in OS patients using whole exome sequencing (WES) to find new potential prognostic biomarkers and therapeutic targets. We performed WES on formalin-fixed paraffin-embedded (FFPE) biopsy materials collected from 19 patients affected by conventional high-grade OS. The clinical and genetic data were analyzed according to response to therapy, presence of metastasis and disease status. By comparing good and poor responders to neoadjuvant therapy, we detected a clear prevalence of mutations in the ARID1A, CREBBP, BRCA2 and RAD50 genes in poor responders that negatively influence the progression-free survival time. Moreover, higher tumor mutational burden values correlated with worse prognosis. The identification of mutations in ARID1A, CREBBP, BRCA2 and RAD50 may support the use of a more specific therapy for tumors harboring these alterations. In particular, BRCA2 and RAD50 are involved in homologous recombination repair, and could thus be used as specific therapy targets of inhibitors of the enzyme Poly ADP Ribose Polymerase (PARP). Finally, tumor mutational burden is found to be a potential prognostic marker for OS.

Genomic profiling informs diagnoses and treatment in vascular anomalies

Vascular anomalies are malformations or tumors of the blood or lymphatic vasculature and can be life-threatening. Although molecularly targeted therapies can be life-saving, identification of the molecular etiology is often impeded by lack of accessibility to affected tissue samples, mosaicism or insufficient sequencing depth. In a cohort of 356 participants with vascular anomalies, including 104 with primary complex lymphatic anomalies (pCLAs), DNA from CD31+ cells isolated from lymphatic fluid or cell-free DNA from lymphatic fluid or plasma underwent ultra-deep sequencing thereby uncovering pathogenic somatic variants down to a variant allele fraction of 0.15%. A molecular diagnosis, including previously undescribed genetic causes, was obtained in 41% of participants with pCLAs and 72% of participants with other vascular malformations, leading to a new medical therapy for 63% (43/69) of participants and resulting in improvement in 63% (35/55) of participants on therapy. Taken together, these data support the development of liquid biopsy-based diagnostic techniques to identify previously undescribed genotype-phenotype associations and guide medical therapy in individuals with vascular anomalies.

BRCA1-methylated triple negative breast cancers previously exposed to neoadjuvant chemotherapy form RAD51 foci and respond poorly to olaparib

Background

About 15% of Triple-Negative-Breast-Cancer (TNBC) present silencing of the BRCA1 promoter methylation and are assumed to be Homologous Recombination Deficient (HRD). BRCA1-methylated (BRCA1-Me) TNBC could, thus, be eligible to treatment based on PARP-inhibitors or Platinum salts. However, their actual HRD status is discussed, as these tumors are suspected to develop resistance after chemotherapy exposure.

Methods

We interrogated the sensitivity to olaparib vs. carboplatin of 8 TNBC Patient-Derived Xenografts (PDX) models. Four PDX corresponded to BRCA1-Me, of which 3 were previously exposed to NeoAdjuvant-Chemotherapy (NACT). The remaining PDX models corresponded to two BRCA1-mutated (BRCA1-Mut) and two BRCA1-wild type PDX that were respectively included as positive and negative controls. The HRD status of our PDX models was assessed using both genomic signatures and the functional BRCA1 and RAD51 nuclear foci formation assay. To assess HR restoration associated with olaparib resistance, we studied pairs of BRCA1 deficient cell lines and their resistant subclones.

Results

The 3 BRCA1-Me PDX that had been exposed to NACT responded poorly to olaparib, likewise BRCA1-WT PDX. Contrastingly, 3 treatment-naïve BRCA1-deficient PDX (1 BRCA1-Me and 2 BRCA1-mutated) responded to olaparib. Noticeably, the three olaparib-responsive PDX scored negative for BRCA1- and RAD51-foci, whereas all non-responsive PDX models, including the 3 NACT-exposed BRCA1-Me PDX, scored positive for RAD51-foci. This suggested HRD in olaparib responsive PDX, while non-responsive models were HR proficient. These results were consistent with observations in cell lines showing a significant increase of RAD51-foci in olaparib-resistant subclones compared with sensitive parental cells, suggesting HR restoration in these models.

Conclusion

Our results thus support the notion that the actual HRD status of BRCA1-Me TNBC, especially if previously exposed to chemotherapy, may be questioned and should be verified using the BRCA1- and RAD51-foci assay.

Molecular Heterogeneity of Pediatric AML with Atypical Promyelocytes Accumulation in Children—A Single Center Experience

Acute promyelocytic leukemia (APL) pathogenesis is based on RARA gene translocations, which are of high importance in the diagnosis of and proper therapy selection for APL. However, in some cases acute myeloid leukemia (AML) demonstrates APL-like morphological features such as atypical promyelocytes accumulation. This type of AML is characterized by the involvement of other RAR family members or completely different genes. In the present study, we used conventional karyotyping, FISH and high-throughput sequencing in a group of 271 de novo AML with atypical promyelocytes accumulation. Of those, 255 cases were shown to carry a typical chromosomal translocation t(15;17)(q24;q21) with PML::RARA chimeric gene formation (94.1%). Other RARA-positive cases exhibited cryptic PML::RARA fusion without t(15;17)(q24;q21) (1.8%, n = 5) and variant t(5;17)(q35;q21) translocation with NPM1::RARA chimeric gene formation (1.5%, n = 4). However, 7 RARA-negative AMLs with atypical promyelocytes accumulation were also discovered. These cases exhibited TBL1XR1::RARB and KMT2A::SEPT6 fusions as well as mutations, e.g., NPM1 insertion and non-recurrent chromosomal aberrations. Our findings demonstrate the genetic diversity of AML with APL-like morphological features, which is of high importance for successful therapy implementation.

Genotypes and phenotypes heterogeneity in PIK3CA-related overgrowth spectrum and overlapping conditions: 150 novel patients and systematic review of 1007 patients with PIK3CA pathogenetic variants

BACKGROUND

Postzygotic activating PIK3CA variants cause several phenotypes within the PIK3CA-related overgrowth spectrum (PROS). Variant strength, mosaicism level, specific tissue involvement and overlapping disorders are responsible for disease heterogeneity. We explored these factors in 150 novel patients and in an expanded cohort of 1007 PIK3CA-mutated patients, analysing our new data with previous literature to give a comprehensive picture.

METHODS

We performed ultradeep targeted next-generation sequencing (NGS) on DNA from skin biopsy, buccal swab or blood using a panel including phosphatidylinositol 3-kinase/AKT/mammalian target of rapamycin pathway genes and GNAQ, GNA11, RASA1 and TEK. Additionally, 914 patients previously reported were systematically reviewed.

RESULTS

93 of our 150 patients had PIK3CA pathogenetic variants. The merged PROS cohort showed that PIK3CA variants span thorough all gene domains, some were exclusively associated with specific PROS phenotypes: weakly activating variants were associated with central nervous system (CNS) involvement, and strongly activating variants with extra-CNS phenotypes. Among the 57 with a wild-type PIK3CA allele, 11 patients with overgrowth and vascular malformations overlapping PROS had variants in GNAQ, GNA11, RASA1 or TEK.

CONCLUSION

We confirm that (1) molecular diagnostic yield increases when multiple tissues are tested and by enriching NGS panels with genes of overlapping 'vascular' phenotypes; (2) strongly activating PIK3CA variants are found in affected tissue, rarely in blood: conversely, weakly activating mutations more common in blood; (3) weakly activating variants correlate with CNS involvement, strong variants are more common in cases without; (4) patients with vascular malformations overlapping those of PROS can harbour variants in genes other than PIK3CA.

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.

Postzygotic mutations and where to find them - Recent advances and future implications in the field of non-neoplastic somatic mosaicism

The technological progress of massively parallel sequencing (MPS) has triggered a remarkable development in the research on postzygotic mutations. Although the overwhelming majority of studies in the field focus on oncogenesis, non-neoplastic diseases are attracting more and more attention. The aim of this review was to summarize some of the most recent findings in the field of somatic mosaicism in diseases other than neoplastic events. We discuss the abundance and role of postzygotic mutations, with a special emphasis on disorders which occur only in a mosaic form (obligatory mosaic diseases; OMDs). Based on the list of OMDs compiled from the published literature and three databases (OMIM, Orphanet and MosaicBase), we demonstrate the prevalence of cancer-related genes across OMDs and suggest other sources to further explore OMDs and OMD-related genes. Additionally, we comment on some practical aspects related to mosaic diseases, such as approaches to tissue sampling, the MPS coverage required to detect variants at a very low frequency, as well as on bioinformatic and molecular tools dedicated to detect somatic mutations in MPS data.

Deep oncopanel sequencing reveals within block position-dependent quality degradation in FFPE processed samples

BACKGROUND

Clinical laboratories routinely use formalin-fixed paraffin-embedded (FFPE) tissue or cell block cytology samples in oncology panel sequencing to identify mutations that can predict patient response to targeted therapy. To understand the technical error due to FFPE processing, a robustly characterized diploid cell line was used to create FFPE samples with four different pre-tissue processing formalin fixation times. A total of 96 FFPE sections were then distributed to different laboratories for targeted sequencing analysis by four oncopanels, and variants resulting from technical error were identified.

RESULTS

Tissue sections that fail more frequently show low cellularity, lower than recommended library preparation DNA input, or target sequencing depth. Importantly, sections from block surfaces are more likely to show FFPE-specific errors, akin to "edge effects" seen in histology, while the inner samples display no quality degradation related to fixation time.

CONCLUSIONS

To assure reliable results, we recommend avoiding the block surface portion and restricting mutation detection to genomic regions of high confidence.

Advanced Molecular Characterisation in Relapsed and Refractory Paediatric Acute Leukaemia, the Key for Personalised Medicine

Relapsed and refractory (R/r) disease in paediatric acute leukaemia remains the first reason for treatment failure. Advances in molecular characterisation can ameliorate the identification of genetic biomarkers treatment strategies for this disease, especially in high-risk patients. The purpose of this study was to analyse a cohort of R/r children diagnosed with acute lymphoblastic (ALL) or myeloid (AML) leukaemia in order to offer them a targeted treatment if available. Advanced molecular characterisation of 26 patients diagnosed with R/r disease was performed using NGS, MLPA, and RT-qPCR. The clinical relevance of the identified alterations was discussed in a multidisciplinary molecular tumour board (MTB). A total of 18 (69.2%) patients were diagnosed with B-ALL, 4 (15.4%) with T-ALL, 3 (11.5%) with AML and 1 patient (3.8%) with a mixed-phenotype acute leukaemia (MPL). Most of the patients had relapsed disease (88%) at the time of sample collection. A total of 17 patients (65.4%) were found to be carriers of a druggable molecular alteration, 8 of whom (47%) received targeted therapy, 7 (87.5%) of them in addition to hematopoietic stem cell transplantation (HSCT). Treatment response and disease control were achieved in 4 patients (50%). In conclusion, advanced molecular characterisation and MTB can improve treatment and outcome in paediatric R/r acute leukaemias.

Combinatorial and Machine Learning Approaches for Improved Somatic Variant Calling From Formalin-Fixed Paraffin-Embedded Genome Sequence Data

Formalin fixation of paraffin-embedded tissue samples is a well-established method for preserving tissue and is routinely used in clinical settings. Although formalin-fixed, paraffin-embedded (FFPE) tissues are deemed crucial for research and clinical applications, the fixation process results in molecular damage to nucleic acids, thus confounding their use in genome sequence analysis. Methods to improve genomic data quality from FFPE tissues have emerged, but there remains significant room for improvement. Here, we use whole-genome sequencing (WGS) data from matched Fresh Frozen (FF) and FFPE tissue samples to optimize a sensitive and precise FFPE single nucleotide variant (SNV) calling approach. We present methods to reduce the prevalence of false-positive SNVs by applying combinatorial techniques to five publicly available variant callers. We also introduce FFPolish, a novel variant classification method that efficiently classifies FFPE-specific false-positive variants. Our combinatorial and statistical techniques improve precision and F1 scores compared to the results of publicly available tools when tested individually.

Pathogenic variant-based preconception carrier screening in the Israeli Jewish population

Preconception carrier screening allows identification of couples at risk to have offspring with autosomal recessive and X-linked disorders. In a current multiethnic world, screening based on self-reported ancestry has limitations. Here we describe the findings of a comprehensive pan-ethnic variant-based carrier screening, using the Israeli Jewish population as a model. The cohort included 1696 individuals (848 couples) tested with the 'MyScreen' multigene panel. The panel covers 1206 variants spanning 385 genes, known in different Jewish ethnicities and local Arab, Druze and Bedouin populations. Out of these, 205 variants in 143 genes are Jewish founder variants. We identified 859 (50.6%), carriers of at least one variant in 151 genes. Importantly, 569 (66.2%) of carriers could be missed by the current Israeli screening program. In total, 1:40 (2.5%) of carrier couples were identified by the 'MyScreen' panel, compared with 1:144 (0.7%) found by the ethnicity-based screening. Surprisingly, 90 individuals (10.5%) were carriers of variants "unexpected" for their reported origin, and 16 variants were previously unreported in Jewish patients. Our results support the advantages of variant-based comprehensive carrier screening for detection of carriers and at-risk couples in a diverse population with many founder disease-causing variants.

Activated interleukin-7 receptor signaling drives B-cell acute lymphoblastic leukemia in mice

Philadelphia chromosome-like acute lymphoblastic leukemia (Ph-like ALL) is a high-risk subtype of B-ALL often associated with genetic variants that alter cytokine receptor signaling, including mutations in the interleukin-7 receptor (IL7R). To investigate whether IL7R variants are leukemia-initiating, we built mouse models expressing activated Il7r (aIL7R). B-cell intrinsic aIL7R mice developed spontaneous B-ALL, demonstrating sufficiency of Il7r activating mutations in leukemogenesis. Concomitant introduction of a knock-out allele in the associated adapter protein Lnk (encoded by Sh2b3) or a dominant-negative variant of the transcription factor Ikaros (Ikzf1) increased disease penetrance. The resulting murine leukemias displayed monoclonality and recurrent somatic Kras mutations and efficiently engrafted into immunocompetent mice. Phosphoproteomic analyses of aIL7R leukemic cells revealed constitutive Stat5 signaling and B cell receptor (BCR)-like signaling despite the absence of surface pre-BCR. Finally, in vitro treatment of aIL7R leukemic B-cells with Jak, mTOR, or Syk inhibitors blocked growth, confirming that each pathway is active in this mouse model of IL7R-driven B-ALL.

Ideafix: a decision tree-based method for the refinement of variants in FFPE DNA sequencing data

Increasingly, treatment decisions for cancer patients are being made from next-generation sequencing results generated from formalin-fixed and paraffin-embedded (FFPE) biopsies. However, this material is prone to sequence artefacts that cannot be easily identified. In order to address this issue, we designed a machine learning-based algorithm to identify these artefacts using data from >1 600 000 variants from 27 paired FFPE and fresh-frozen breast cancer samples. Using these data, we assembled a series of variant features and evaluated the classification performance of five machine learning algorithms. Using leave-one-sample-out cross-validation, we found that XGBoost (extreme gradient boosting) and random forest obtained AUC (area under the receiver operating characteristic curve) values >0.86. Performance was further tested using two independent datasets that resulted in AUC values of 0.96, whereas a comparison with previously published tools resulted in a maximum AUC value of 0.92. The most discriminating features were read pair orientation bias, genomic context and variant allele frequency. In summary, our results show a promising future for the use of these samples in molecular testing. We built the algorithm into an R package called Ideafix (DEAmination FIXing) that is freely available at https://github.com/mmaitenat/ideafix.

Genetic, Clinicopathological, and Radiological Features of Intrahepatic Cholangiocarcinoma with Ductal Plate Malformation Pattern

Background/Aims

Intrahepatic cholangiocarcinoma (iCCA) with a ductal plate malformation (DPM) pattern is a recently recognized rare variant. The genomic profile of iCCA with DPM pattern needs to be elucidated.

Methods

Cases of iCCA with DPM pattern were retrospectively reviewed based on the medical records, pathology slides, and magnetic resonance imaging (MRI) reports collected between 2010 to 2019 at a single center. Massive parallel sequencing was performed for >500 cancer-related genes.

Results

From a total of 175 iCCAs, five (2.9%) cases of iCCA with DPM pattern were identified. All cases were of the small duct type, and background liver revealed chronic B viral or alcoholic hepatitis. Three iCCAs with DPM pattern harbored MRI features favoring the diagnosis of hepatocellular carcinoma, whereas nonspecific imaging features were observed in two cases. All patients were alive without recurrence during an average follow-up period of 57 months. Sequencing data revealed 64 mutated genes in the five cases, among which FGFR2 and PTPRT were most frequently mutated (three cases each) including an FGFR2-TNC fusion in one case. Mutations in ARID1A and CDKN2A were found in two cases, and mutations in TP53, BAP1, ATM, NF1, and STK11 were observed in one case each. No IDH1, KRAS, or PBRM1 mutations were found.

Conclusions

iCCAs with DPM pattern have different clinico-radio-pathologic and genetic characteristics compared to conventional iCCAs. Moreover, FGFR2 and ARID1A variants were identified. Altogether, these findings further suggest that iCCA with DPM pattern represents a specific subtype of small duct type iCCA.

STIC: Predicting Single Nucleotide Variants and Tumor Purity in Cancer Genome

Single nucleotide variant (SNV) plays an important role in cellular proliferation and tumorigenesis in various types of human cancer. Next-generation sequencing (NGS) has provided high-throughput data at an unprecedented resolution to predict SNVs. Currently, there exist many computational methods for either germline or somatic SNV discovery from NGS data, but very few of them are versatile enough to adapt to any situations. In the absence of matched normal samples, the prediction of somatic SNVs from single-tumor samples becomes considerably challenging, especially when the tumor purity is unknown. Here, we propose a new approach, STIC, to predict somatic SNVs and estimate tumor purity from NGS data without matched normal samples. The main features of STIC include: (1) extracting a set of SNV-relevant features on each site and training the BP neural network algorithm on the features to predict SNVs; (2) creating an iterative process to distinguish somatic SNVs from germline ones by disturbing allele frequency; and (3) establishing a reasonable relationship between tumor purity and allele frequencies of somatic SNVs to accurately estimate the purity. We quantitatively evaluate the performance of STIC on both simulation and real sequencing datasets, the results of which indicate that STIC outperforms competing methods.

Proteotranscriptomic classification and characterization of pancreatic neuroendocrine neoplasms

Pancreatic neuroendocrine neoplasms (PNENs) are biologically and clinically heterogeneous. Here, we use a multi-omics approach to uncover the molecular factors underlying this heterogeneity. Transcriptomic analysis of 84 PNEN specimens, drawn from two cohorts, is substantiated with proteomic profiling and identifies four subgroups: Proliferative, PDX1-high, Alpha cell-like and Stromal/Mesenchymal. The Proliferative subgroup, consisting of both well- and poorly differentiated specimens, is associated with inferior overall survival probability. The PDX1-high and Alpha cell-like subgroups partially resemble previously described subtypes, and we further uncover distinctive metabolism-related features in the Alpha cell-like subgroup. The Stromal/Mesenchymal subgroup exhibits molecular characteristics of YAP1/WWTR1(TAZ) activation suggestive of Hippo signaling pathway involvement in PNENs. Whole-exome sequencing reveals subgroup-enriched mutational differences, supported by activity inference analysis, and identifies hypermorphic proto-oncogene variants in 14.3% of sequenced PNENs. Our study reveals differences in cellular signaling axes that provide potential directions for PNEN patient stratification and treatment strategies.

The prevalence of homologous recombination deficiency (HRD) in various solid tumors and the role of HRD as a single biomarker to immune checkpoint inhibitors

Purpose

Homologous recombination deficiency (HRD) is related to tumorigenesis. Currently, the possibility of HRD as a prognostic biomarker to immune checkpoint inhibitors is unknown. We aimed to investigate whether HRD has potential as a biomarker for immunotherapy.

Methods

The status of homologous recombination deficiency (HRD) was assessed with the next-generation sequencing (NGS) TruSight^™ Oncology 500 assay in 501 patients with advanced solid tumor including gastrointestinal (GI), genitourinary (GU), or rare cancer. Results: among the 501 patients, HRD was observed as follows: 74.7% (347/501) patients; GU cancer (92.0%, 23 of 25), colorectal cancer (CRC) (86.1%, 130 of 151), hepatocellular carcinoma (HCC) (83.3%, 10 of 12), pancreatic cancer (PC) (76.2%, 32 of 42), biliary tract cancer (BTC) (75.0%, 36 of 48), sarcoma (65.0%, 39 of 60), melanoma (52.4%, 11 of 21), other GI cancers (50.0%, 11 of 22), and rare cancer (50.0%, 2 of 4). Sixty-five of the 501 patients had received immune checkpoint inhibitors (ICIs) during the course of the disease. Tumor types of 65 patients treated with ICIs are as follows: melanoma (95.2%, 20 of 21), HCC (33.3%, 4 of 12), rare cancer (25.0%, 1 of 4), GC (12.2%, 14 of 116), BTC (10.4%, 5 of 48), and sarcoma (5.0%, 3 of 60). The most frequently reported mutations were BRCA2 (n = 90), ARID1A (n = 77), ATM (n = 71), BARD1 (n = 67). Patients without HRD exhibited an objective response rate (ORR) of 33.3% (4 of 12), and patients with HRD exhibited an ORR of 34.0% (18 of 53). There was no significant difference in ORR between patients with and without HRD (P = 0.967). Progression-free survival (PFS) was 6.5 months (95% CI 0.000–16.175) in patients without HRD and 4.1 months (95% CI 2.062–6.138) in patients with HRD, revealing no statistical significance (P = 0.441).

Conclusion

Herein, we reported the status of HRD using a cancer-panel for various solid tumor patients in routine clinical practice and demonstrated that HRD as a single biomarker was not sufficient to predict efficacy of ICIs in solid tumor patients.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00432-021-03781-6.

Nanopanel2 calls phased low-frequency variants in Nanopore panel sequencing data

MOTIVATION

Clinical decision making is increasingly guided by accurate and recurrent determination of presence and frequency of (somatic) variants and their haplotype through panel sequencing of disease-relevant genomic regions. Haplotype calling (phasing), however, is difficult and error prone unless variants are located on the same read which limits the ability of short-read sequencing to detect, e.g., co-occurrence of drug-resistance variants. Long-read panel sequencing enables direct phasing of amplicon variants besides having multiple other benefits, however, high error rates of current technologies prevented their applicability in the past.

RESULTS

We have developed Nanopanel2, a variant caller for Nanopore panel sequencing data. Nanopanel2 works directly on base-called FAST5 files and uses allele probability distributions and several other filters to robustly separate true from false positive (FP) calls. It effectively calls SNVs and INDELs with variant allele frequencies as low as 1% and 5% respectively and produces only few low-frequency false-positive calls (∼1 FP call with VAF¡5% per kb amplicon). Haplotype compositions are then determined by direct phasing. Nanopanel2 is the first somatic variant caller for Nanopore data, enabling accurate, fast (turnaround <48h) and cheap (sequencing costs ∼10$/sample) diagnostic workflows.

AVAILABILITY

The data for this study have been deposited at zenodo.org under DOIs accession numbers 4110691 and 4110698. Nanopanel2 is open source and available at https://github.com/popitsch/nanopanel2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comparison of the oncogenomic landscape of canine and feline hemangiosarcoma shows novel parallels with human angiosarcoma

Angiosarcoma (AS) is a highly aggressive tumor of blood and lymphatic vessels in humans that shares many similarities with spontaneously occurring hemangiosarcoma (HSA) in dogs and cats. To investigate the genetic suitability of HSA as a model for AS, we sequenced ∼1000 cancer genes in 41 cases of HSA and matched germline tissue: 15 canine visceral HSAs, 13 canine skin HSAs and 13 feline skin HSAs. Analysis of visceral HSAs from dogs presenting with concurrent splenic and cardiac neoplasms showed that the tumors were not independent primaries, consistent with the highly metastatic nature of HSA. Comparison of HSA to AS revealed that several driver genes were recurrently mutated in both species, such as TP53, PIK3CA, ATRX, GRIN2A and LRP1B. Similar to AS, a UV mutational signature was found in a subset of canine cutaneous HSAs and both species show differing mutational profiles between tissue sites. Our characterization of canine and feline HSA demonstrates many important parallels to AS and provides hope that future studies on these cancers will benefit of all three species.

Integrated mutational landscape analysis of uterine leiomyosarcomas

Uterine leiomyosarcomas (uLMS) are aggressive tumors arising from the smooth muscle layer of the uterus. We analyzed 83 uLMS sample genetics, including 56 from Yale and 27 from The Cancer Genome Atlas (TCGA). Among them, a total of 55 Yale samples including two patient-derived xenografts (PDXs) and 27 TCGA samples have whole-exome sequencing (WES) data; 10 Yale and 27 TCGA samples have RNA-sequencing (RNA-Seq) data; and 11 Yale and 10 TCGA samples have whole-genome sequencing (WGS) data. We found recurrent somatic mutations in TP53, MED12, and PTEN genes. Top somatic mutated genes included TP53, ATRX, PTEN, and MEN1 genes. Somatic copy number variation (CNV) analysis identified 8 copy-number gains, including 5p15.33 (TERT), 8q24.21 (C-MYC), and 17p11.2 (MYOCD, MAP2K4) amplifications and 29 copy-number losses. Fusions involving tumor suppressors or oncogenes were deetected, with most fusions disrupting RB1, TP53, and ATRX/DAXX, and one fusion (ACTG2-ALK) being potentially targetable. WGS results demonstrated that 76% (16 of 21) of the samples harbored chromoplexy and/or chromothripsis. Clinically actionable mutational signatures of homologous-recombination DNA-repair deficiency (HRD) and microsatellite instability (MSI) were identified in 25% (12 of 48) and 2% (1 of 48) of fresh frozen uLMS, respectively. Finally, we found olaparib (PARPi; P = 0.002), GS-626510 (C-MYC/BETi; P < 0.000001 and P = 0.0005), and copanlisib (PIK3CAi; P = 0.0001) monotherapy to significantly inhibit uLMS-PDXs harboring derangements in C-MYC and PTEN/PIK3CA/AKT genes (LEY11) and/or HRD signatures (LEY16) compared to vehicle-treated mice. These findings define the genetic landscape of uLMS and suggest that a subset of uLMS may benefit from existing PARP-, PIK3CA-, and C-MYC/BET-targeted drugs.

Phase 2 Clinical Trial of Infusing Haploidentical K562-mb15-41BBL-Activated and Expanded Natural Killer Cells as Consolidation Therapy for Pediatric Acute Myeloblastic Leukemia

BACKGROUND

Acute myeloid leukemia (AML) accounts for approximately 20% of pediatric leukemia cases; 30% of these patients experience relapse. The antileukemia properties of natural killer (NK) cells and their safety profile have been reported in AML therapy. We proposed a phase 2, open, prospective, multicenter, nonrandomized clinical trial for the adoptive infusion of haploidentical K562-mb15-41BBL-activated and expanded NK (NKAE) cells as a consolidation strategy for children with favorable and intermediate risk AML in first complete remission after chemotherapy (NCT02763475).

PATIENTS AND METHODS

Before the NKAE cell infusion, patients underwent a lymphodepleting regimen. After the NKAE cell infusion, patients were administered low doses (1 × 10⁶/IU/m²) of subcutaneous interleukin-2. The primary study endpoint was AML relapse-free survival. We needed to include 35 patients to demonstrate a 50% reduction in relapses.

RESULTS

Seven patients (median age, 7.4 years; range, 0.78-15.98 years) were administered 13 infusions of NKAE cells, with a median of 36.44 × 10⁶ cells/kg (range, 6.92 × 10⁶ to 193.2 × 10⁶ cells/kg). We observed chimerism in 4 patients (median chimerism, 0.065%; range, 0.05-0.27%). After a median follow-up of 33 months, the disease of 6 patients (85.7%) remained in complete remission. The 3-year overall survival was 83.3% (95% confidence interval, 68.1-98.5), and the cumulative 3-year relapse rate was 28.6% (95% confidence interval, 11.5-45.7). The study was terminated early because of low patient recruitment.

CONCLUSION

This study emphasizes the difficulties in recruiting patients for cell therapy trials, though NKAE cell infusion is safe and feasible. However, we cannot draw any conclusions regarding efficacy because of the small number of included patients and insufficient biological markers.

Low frequency mitochondrial DNA heteroplasmy SNPs in blood, retina, and [RPE+choroid] of age-related macular degeneration subjects

Purpose

Mitochondrial (mt) DNA damage is associated with age-related macular degeneration (AMD) and other human aging diseases. This study was designed to quantify and characterize mtDNA low-frequency heteroplasmy single nucleotide polymorphisms (SNPs) of three different tissues isolated from AMD subjects using Next Generation Sequencing (NGS) technology.

Methods

DNA was extracted from neural retina, [RPE+choroid] and blood from three deceased age-related macular degeneration (AMD) subjects. Entire mitochondrial genomes were analyzed for low-frequency heteroplasmy SNPs using NGS technology that independently sequenced both mtDNA strands. This deep sequencing method (average sequencing depth of 30,000; range 1,000–100,000) can accurately differentiate low-frequency heteroplasmy SNPs from DNA modification artifacts. Twenty-three ‘hot-spot’ heteroplasmy mtDNA SNPs were analyzed in 222 additional blood samples.

Results

Germline homoplasmy SNPs that defined mtDNA haplogroups were consistent in the three tissues of each subject. Analyses of SNPs with <40% heteroplasmy revealed the blood had significantly greater numbers of heteroplasmy SNPs than retina alone (p≤0.05) or retina+choroid combined (p = 0.008). Twenty-three ‘hot-spot’ mtDNA heteroplasmy SNPs were present, with three being non-synonymous (amino acid change). Four ‘hot-spot’ heteroplasmy SNPs (m.1120C>T, m.1284T>C, m.1556C>T, m.7256C>T) were found in additional samples (n = 222). Five heteroplasmy SNPs (m.4104A>G, m.5320C>T, m.5471G>A, m.5474A>G, m.5498A>G) declined with age. Two heteroplasmy SNPs (m.13095T>C, m.13105A>G) increased in AMD compared to Normal samples. In the heteroplasmy SNPs, very few transversion mutations (purine to pyrimidine or vice versa, associated with oxidative damage) were found and the majority were transition changes (purine to purine or pyrimidine to pyrimidine, associated with replication errors).

Conclusion

Within an individual, the blood, retina and [RPE+choroid] contained identical homoplasmy SNPs representing inherited germline mtDNA haplogroup. NGS methodology showed significantly more mtDNA heteroplasmy SNPs in blood compared to retina and [RPE+choroid], suggesting the latter tissues have substantial protection. Significantly higher heteroplasmy levels of m.13095T>C and m.13105A>G may represent potential AMD biomarkers. Finally, high levels of transition mutations suggest that accumulation of heteroplasmic SNPs may occur through replication errors rather than oxidative damage.

MSH2 knock-down shows CTG repeat stability and concomitant upstream demethylation at the DMPK locus in myotonic dystrophy type 1 human embryonic stem cells

Myotonic dystrophy type 1 (DM1) is caused by expansion of a CTG repeat in the DMPK gene, where expansion size and somatic mosaicism correlates with disease severity and age of onset. While it is known that the mismatch repair protein MSH2 contributes to the unstable nature of the repeat, its role on other disease-related features, such as CpG methylation upstream of the repeat, is unknown. In this study, we investigated the effect of an MSH2 knock-down (MSH2KD) on both CTG repeat dynamics and CpG methylation pattern in human embryonic stem cells (hESC) carrying the DM1 mutation. Repeat size in MSH2 wild-type (MSH2WT) and MSH2KD DM1 hESC was determined by PacBio sequencing and CpG methylation by bisulfite massive parallel sequencing. We found stabilization of the CTG repeat concurrent with a gradual loss of methylation upstream of the repeat in MSH2KD cells, while the repeat continued to expand and upstream methylation remained unchanged in MSH2WT control lines. Repeat instability was re-established and biased towards expansions upon MSH2 transgenic re-expression in MSH2KD lines while upstream methylation was not consistently re-established. We hypothesize that the hypermethylation at the mutant DM1 locus is promoted by the MMR machinery and sustained by a constant DNA repair response, establishing a potential mechanistic link between CTG repeat instability and upstream CpG methylation. Our work represents a first step towards understanding how epigenetic alterations and repair pathways connect and contribute to the DM1 pathology.

Optimizing performance of GATK workflows using Apache Arrow In-Memory data framework

BACKGROUND

Immense improvements in sequencing technologies enable producing large amounts of high throughput and cost effective next-generation sequencing (NGS) data. This data needs to be processed efficiently for further downstream analyses. Computing systems need this large amounts of data closer to the processor (with low latency) for fast and efficient processing. However, existing workflows depend heavily on disk storage and access, to process this data incurs huge disk I/O overheads. Previously, due to the cost, volatility and other physical constraints of DRAM memory, it was not feasible to place large amounts of working data sets in memory. However, recent developments in storage-class memory and non-volatile memory technologies have enabled computing systems to place huge data in memory to process it directly from memory to avoid disk I/O bottlenecks. To exploit the benefits of such memory systems efficiently, proper formatted data placement in memory and its high throughput access is necessary by avoiding (de)-serialization and copy overheads in between processes. For this purpose, we use the newly developed Apache Arrow, a cross-language development framework that provides language-independent columnar in-memory data format for efficient in-memory big data analytics. This allows genomics applications developed in different programming languages to communicate in-memory without having to access disk storage and avoiding (de)-serialization and copy overheads.

IMPLEMENTATION

We integrate Apache Arrow in-memory based Sequence Alignment/Map (SAM) format and its shared memory objects store library in widely used genomics high throughput data processing applications like BWA-MEM, Picard and GATK to allow in-memory communication between these applications. In addition, this also allows us to exploit the cache locality of tabular data and parallel processing capabilities through shared memory objects.

RESULTS

Our implementation shows that adopting in-memory SAM representation in genomics high throughput data processing applications results in better system resource utilization, low number of memory accesses due to high cache locality exploitation and parallel scalability due to shared memory objects. Our implementation focuses on the GATK best practices recommended workflows for germline analysis on whole genome sequencing (WGS) and whole exome sequencing (WES) data sets. We compare a number of existing in-memory data placing and sharing techniques like ramDisk and Unix pipes to show how columnar in-memory data representation outperforms both. We achieve a speedup of 4.85x and 4.76x for WGS and WES data, respectively, in overall execution time of variant calling workflows. Similarly, a speedup of 1.45x and 1.27x for these data sets, respectively, is achieved, as compared to the second fastest workflow. In some individual tools, particularly in sorting, duplicates removal and base quality score recalibration the speedup is even more promising.

AVAILABILITY

The code and scripts used in our experiments are available in both container and repository form at: https://github.com/abs-tudelft/ArrowSAM .

Accuracy and reproducibility of somatic point mutation calling in clinical-type targeted sequencing data

Background

Treating cancer depends in part on identifying the mutations driving each patient’s disease. Many clinical laboratories are adopting high-throughput sequencing for assaying patients’ tumours, applying targeted panels to formalin-fixed paraffin-embedded tumour tissues to detect clinically-relevant mutations. While there have been some benchmarking and best practices studies of this scenario, much variant calling work focuses on whole-genome or whole-exome studies, with fresh or fresh-frozen tissue. Thus, definitive guidance on best choices for sequencing platforms, sequencing strategies, and variant calling for clinical variant detection is still being developed.

Methods

Because ground truth for clinical specimens is rarely known, we used the well-characterized Coriell cell lines GM12878 and GM12877 to generate data. We prepared samples to mimic as closely as possible clinical biopsies, including formalin fixation and paraffin embedding. We evaluated two well-known targeted sequencing panels, Illumina’s TruSight 170 hybrid-capture panel and the amplification-based Oncomine Focus panel. Sequencing was performed on an Illumina NextSeq500 and an Ion Torrent PGM respectively. We performed multiple replicates of each assay, to test reproducibility. Finally, we applied four different freely-available somatic single-nucleotide variant (SNV) callers to the data, along with the vendor-recommended callers for each sequencing platform.

Results

We did not observe major differences in variant calling success within the regions that each panel covers, but there were substantial differences between callers. All had high sensitivity for true SNVs, but numerous and non-overlapping false positives. Overriding certain default parameters to make them consistent between callers substantially reduced discrepancies, but still resulted in high false positive rates. Intersecting results from multiple replicates or from different variant callers eliminated most false positives, while maintaining sensitivity.

Conclusions

Reproducibility and accuracy of targeted clinical sequencing results depend less on sequencing platform and panel than on variability between replicates and downstream bioinformatics. Differences in variant callers’ default parameters are a greater influence on algorithm disagreement than other differences between the algorithms. Contrary to typical clinical practice, we recommend employing multiple variant calling pipelines and/or analyzing replicate samples, as this greatly decreases false positive calls.

Molecular characterization of precise in vivo targeted gene integration in human cells using AAVHSC15

Targeted gene integration via precise homologous recombination (HR)-based gene editing has the potential to correct genetic diseases. AAV (adeno-associated virus) can mediate nuclease-free gene integration at a disease-causing locus. Therapeutic application of AAV gene integration requires quantitative molecular characterization of the edited sequence that overcome technical obstacles such as excess episomal vector genomes and lengthy homology arms. Here we describe a novel molecular methodology that utilizes quantitative next-generation sequencing to characterize AAV-mediated targeted insertion and detects the presence of unintended mutations. The methods described here quantify targeted insertion and query the entirety of the target locus for the presence of insertions, deletions, single nucleotide variants (SNVs) and integration of viral components such as inverted terminal repeats (ITR). Using a humanized liver murine model, we demonstrate that hematopoietic stem-cell derived AAVHSC15 mediates in vivo targeted gene integration into human chromosome 12 at the PAH (phenylalanine hydroxylase) locus at 6% frequency, with no sign of co-incident random mutations at or above a lower limit of detection of 0.5% and no ITR sequences at the integration sites. Furthermore, analysis of heterozygous variants across the targeted locus using the methods described shows a pattern of strand cross-over, supportive of an HR mechanism of gene integration with similar efficiencies across two different haplotypes. Rapid advances in the application of AAV-mediated nuclease-free target integration, or gene editing, as a new therapeutic modality requires precise understanding of the efficiency and the nature of the changes being introduced to the target genome at the molecular level. This work provides a framework to be applied to homologous recombination gene editing platforms for assessment of introduced and natural sequence variation across a target site.

Uncovering Low-Level Maternal Gonosomal Mosaicism in X-Linked Agammaglobulinemia: Implications for Genetic Counseling

X-linked agammaglobulinemia (XLA) is a clinically and genetically well-defined immunodeficiency and the most common form of agammaglobulinemia. It is characterized by susceptibility to recurrent bacterial infections, profound hypogammaglobulinemia, and few or no circulating B cells. XLA is caused by mutations in the BTK gene, which encodes Bruton's tyrosine kinase (BTK). Because of its X-linked recessive inheritance pattern, XLA virtually only affects males, and the mother is the carrier of the mutation in 80–85% of the males with this condition. In the remaining 15–20% of the cases, the affected male is considered to have a de novo mutation. Here, we present the case of a child with a diagnosis of XLA caused by a missense mutation in the BTK gene (c.494G>A/p.C165Y). Apparently, his mother was wild type for this gene, which implied that the mutation was de novo, but careful analysis of Sanger electropherograms and the use of high-coverage massive parallel sequencing revealed low-level maternal gonosomal mosaicism. The mutation was detected in various samples from the mother (blood, urine, buccal swab, and vaginal swab) at a low frequency of 2–5%, and the status of the patient's mutation changed from de novo to inherited. This study underscores the importance of accurately establishing the parents' status on detection of an apparently de novo mutation in a patient, as inadvertent low-level mosaicism may lead to misinterpretation of the risk of recurrence, vital for genetic counseling.

Development and Analytical Validation of a DNA Dual-Strand Approach for the US Food and Drug Administration-Approved Next-Generation Sequencing-Based Praxis Extended RAS Panel for Metastatic Colorectal Cancer Samples

A next-generation sequencing method was developed that can distinguish single-stranded modifications from low-frequency somatic mutations present on both strands of DNA in formalin-fixed paraffin-embedded colorectal cancer samples. We applied this method for analytical validation of the Praxis Extended RAS Panel, a US Food and Drug Administration-approved companion diagnostic for panitumumab, on the Illumina MiSeqDx platform. With the use of the TruSeq amplicon workflow, both strands of DNA from the starting material were interrogated independently. Mutations were reported only if found on both strands; artifacts usually present on only one strand would not be reported. A total of 56 mutations were targeted within the KRAS and NRAS genes. A minimum read depth of 1800× per amplicon is required per sample but averaged >30,000× at maximum multiplexing levels. Analytical validation studies were performed to determine the simultaneous detection of mutations on both strands, reproducibility, assay detection level, precision of the assay across various factors, and the impact of interfering substances. In conclusion, this assay can clearly distinguish single-stranded artifacts from low-frequency mutations. Furthermore, the assay is accurate, precise, and reproducible, can achieve consistent detection of a mutation at 5% mutation frequency, exhibits minimal impact from tested interfering substances, and can simultaneously detect 56 mutations in a single run using 10 samples plus controls.