Germline Elongator mutations in Sonic Hedgehog medulloblastoma.

Cancer genomics has revealed many genes and core molecular processes that contribute to human malignancies, but the genetic and molecular bases of many rare cancers remains unclear. Genetic predisposition accounts for 5 to 10% of cancer diagnoses in children1,2, and genetic events that cooperate with known somatic driver events are poorly understood. Pathogenic germline variants in established cancer predisposition genes have been recently identified in 5% of patients with the malignant brain tumour medulloblastoma3. Here, by analysing all protein-coding genes, we identify and replicate rare germline loss-of-function variants across ELP1 in 14% of paediatric patients with the medulloblastoma subgroup Sonic Hedgehog (MBSHH). ELP1 was the most common medulloblastoma predisposition gene and increased the prevalence of genetic predisposition to 40% among paediatric patients with MBSHH. Parent-offspring and pedigree analyses identified two families with a history of paediatric medulloblastoma. ELP1-associated medulloblastomas were restricted to the molecular SHHα subtype4 and characterized by universal biallelic inactivation of ELP1 owing to somatic loss of chromosome arm 9q. Most ELP1-associated medulloblastomas also exhibited somatic alterations in PTCH1, which suggests that germline ELP1 loss-of-function variants predispose individuals to tumour development in combination with constitutive activation of SHH signalling. ELP1 is the largest subunit of the evolutionarily conserved Elongator complex, which catalyses translational elongation through tRNA modifications at the wobble (U34) position5,6. Tumours from patients with ELP1-associated MBSHH were characterized by a destabilized Elongator complex, loss of Elongator-dependent tRNA modifications, codon-dependent translational reprogramming, and induction of the unfolded protein response, consistent with loss of protein homeostasis due to Elongator deficiency in model systems7-9. Thus, genetic predisposition to proteome instability may be a determinant in the pathogenesis of paediatric brain cancers. These results support investigation of the role of protein homeostasis in other cancer types and potential for therapeutic interference.

Convergent genetic aberrations in murine and human T lineage acute lymphoblastic leukemias.

The lack of predictive preclinical models is a fundamental barrier to translating knowledge about the molecular pathogenesis of cancer into improved therapies. Insertional mutagenesis (IM) in mice is a robust strategy for generating malignancies that recapitulate the extensive inter- and intra-tumoral genetic heterogeneity found in advanced human cancers. While the central role of "driver" viral insertions in IM models that aberrantly increase the expression of proto-oncogenes or disrupt tumor suppressors has been appreciated for many years, the contributions of cooperating somatic mutations and large chromosomal alterations to tumorigenesis are largely unknown. Integrated genomic studies of T lineage acute lymphoblastic leukemias (T-ALLs) generated by IM in wild-type (WT) and Kras mutant mice reveal frequent point mutations and other recurrent non-insertional genetic alterations that also occur in human T-ALL. These somatic mutations are sensitive and specific markers for defining clonal dynamics and identifying candidate resistance mechanisms in leukemias that relapse after an initial therapeutic response. Primary cancers initiated by IM and resistant clones that emerge during in vivo treatment close key gaps in existing preclinical models, and are robust platforms for investigating the efficacy of new therapies and for elucidating how drug exposure shapes tumor evolution and patterns of resistance.

Lipoblastomas presenting in older children and adults: analysis of 22 cases with identification of novel PLAG1 fusion partners.

Lipoblastomas are benign neoplasms of embryonal white fat that typically present in the first 3 years of life and show a lobular arrangement of maturing adipocytes with variable degrees of myxoid change. We systematically studied the clinicopathologic and genetic features of lipoblastomas arising in older children and adults. Cases with a diagnosis of lipoblastoma or maturing lipoblastoma in patients >3 years of age were retrieved from our archives. Immunostaining for CD34 and desmin and molecular studies (FISH, RNA sequencing) were performed. Twenty-two cases (8F; 14M) were identified in patients ranging from 4 to 44 years of age (median 10 years). Sites included extremity (n = 15), head and neck (n = 4), and trunk (n = 3) with tumor sizes varying from 1.6 to 17.5 cm (median 5). Only three tumors had histologic features of "conventional" lipoblastoma. The majority of tumors (n = 14) were composed of variably sized lobules of mature adipose tissue partitioned by thin fibrous septa ("maturing"). The remaining five cases consisted predominantly of bland spindled to plump ovoid cells embedded in a fibrous stroma, with a vaguely plexiform arrangement of small myxoid and adipocytic nodules ("fibroblastic"). CD34 was diffusely positive in all cases tested (21/21), while desmin immunoreactivity was identified in 12 of 21 cases (diffuse = 7, focal = 5). PLAG1 rearrangements were identified in 13 tumors in the entire cohort (59%), including all 5 fibroblastic tumors. RNA sequencing detected eight PLAG1 fusion partners, of which two were known (CHCHD7 and COL3A1) and six were novel (SRSF3, HNRNPC, PCMTD1, YWHAZ, CTDSP2, and PPP2R2A). Twelve cases had follow-up (1-107 months; median 21 months), and no recurrences were reported. Lipoblastomas may occur in older children and adults and may be difficult to recognize due to their predominantly adipocytic or fibrous appearance. Awareness that lipoblastomas may occur in older patients, careful evaluation for foci showing more typical morphologic features, ancillary immunohistochemistry for CD34 and desmin, and molecular genetic studies to identify PLAG1 rearrangements are the keys to recognizing these tumors.

Loss of oncogenic Notch1 with resistance to a PI3K inhibitor in T-cell leukaemia.

Mutations that deregulate Notch1 and Ras/phosphoinositide 3 kinase (PI3K)/Akt signalling are prevalent in T-cell acute lymphoblastic leukaemia (T-ALL), and often coexist. Here we show that the PI3K inhibitor GDC-0941 is active against primary T-ALLs from wild-type and Kras(G12D) mice, and addition of the MEK inhibitor PD0325901 increases its efficacy. Mice invariably relapsed after treatment with drug-resistant clones, most of which unexpectedly had reduced levels of activated Notch1 protein, downregulated many Notch1 target genes, and exhibited cross-resistance to γ-secretase inhibitors. Multiple resistant primary T-ALLs that emerged in vivo did not contain somatic Notch1 mutations present in the parental leukaemia. Importantly, resistant clones upregulated PI3K signalling. Consistent with these data, inhibiting Notch1 activated the PI3K pathway, providing a likely mechanism for selection against oncogenic Notch1 signalling. These studies validate PI3K as a therapeutic target in T-ALL and raise the unexpected possibility that dual inhibition of PI3K and Notch1 signalling could promote drug resistance in T-ALL.

PTEN action in leukaemia dictated by the tissue microenvironment.

PTEN encodes a lipid phosphatase that is underexpressed in many cancers owing to deletions, mutations or gene silencing. PTEN dephosphorylates phosphatidylinositol (3,4,5)-triphosphate, thereby opposing the activity of class I phosphatidylinositol 3-kinases that mediate growth- and survival-factor signalling through phosphatidylinositol 3-kinase effectors such as AKT and mTOR. To determine whether continued PTEN inactivation is required to maintain malignancy, here we generate an RNA interference-based transgenic mouse model that allows tetracycline-dependent regulation of PTEN in a time- and tissue-specific manner. Postnatal Pten knockdown in the haematopoietic compartment produced highly disseminated T-cell acute lymphoblastic leukaemia. Notably, reactivation of PTEN mainly reduced T-cell leukaemia dissemination but had little effect on tumour load in haematopoietic organs. Leukaemia infiltration into the intestine was dependent on CCR9 G-protein-coupled receptor signalling, which was amplified by PTEN loss. Our results suggest that in the absence of PTEN, G-protein-coupled receptors may have an unanticipated role in driving tumour growth and invasion in an unsupportive environment. They further reveal that the role of PTEN loss in tumour maintenance is not invariant and can be influenced by the tissue microenvironment, thereby producing a form of intratumoral heterogeneity that is independent of cancer genotype.

CONSERTING: integrating copy-number analysis with structural-variation detection.

We developed Copy Number Segmentation by Regression Tree in Next Generation Sequencing (CONSERTING), an algorithm for detecting somatic copy-number alteration (CNA) using whole-genome sequencing (WGS) data. CONSERTING performs iterative analysis of segmentation on the basis of changes in read depth and the detection of localized structural variations, with high accuracy and sensitivity. Analysis of 43 cancer genomes from both pediatric and adult patients revealed novel oncogenic CNAs, complex rearrangements and subclonal CNAs missed by alternative approaches.

Knockdown of Hnrnpa0, a del(5q) gene, alters myeloid cell fate in murine cells through regulation of AU-rich transcripts.

The control of mRNA stability plays a central role in orchestrating gene-regulatory networks in hematopoietic cell growth, differentiation and tumorigenesis. HNRNPA0, which encodes an RNA-binding protein shown to regulate transcript stability via binding to the AU-rich elements of mRNAs, is located within the commonly deleted segment of 5q31.2 in myeloid neoplasms with a del(5q), and is expressed at haploinsufficient levels in these patients. We show that HNRNPA0 is normally highly expressed in hematopoietic stem cells and exhibits dynamic changes in expression during the course of differentiation. To model HNRNPA0 haploinsufficiency, we used RNAi interference in primary murine cells and an experimental cell system, and found that reduced Hnrnpa0 expression leads to a shift from monocytic towards granulocytic differentiation. Microarray-based global expression profiling revealed that Hnrnpa0 knockdown disproportionally impacts AU-rich containing transcripts and alters expression of myeloid specification genes. In therapy-related myeloid neoplasms with a del(5q), AU-rich containing mRNAs are enriched in transcripts that encode proteins associated with increased growth and proliferation. Our findings implicate haploinsufficiency of HNRNPA0 as one of the key initiating mutations in the pathogenesis of myeloid neoplasms with a del(5q), and suggest that therapies that target AU-rich elements warrant consideration in efforts to develop new mechanism-based treatment strategies.

SJPedPanel: A pan-cancer gene panel for childhood malignancies.

Background: Large scale genomics projects have identified driver alterations for most childhood cancers that provide reliable biomarkers for clinical diagnosis and disease monitoring using targeted sequencing. However, there is lack of a comprehensive panel that matches the list of known driver genes. Here we fill this gap by developing SJPedPanel for childhood cancers. Results: SJPedPanel covers 5,275 coding exons of 357 driver genes, 297 introns frequently involved in rearrangements that generate fusion oncoproteins, commonly amplified/deleted regions (e.g., MYCN for neuroblastoma, CDKN2A and PAX5 for B-/T-ALL, and SMARCB1 for AT/RT), and 7,590 polymorphism sites for interrogating tumors with aneuploidy, such as hyperdiploid and hypodiploid B-ALL or 17q gain neuroblastoma. We used driver alterations reported from an established real-time clinical genomics cohort (n=253) to validate this gene panel. Among the 485 pathogenic variants reported, our panel covered 417 variants (86%). For 90 rearrangements responsible for oncogenic fusions, our panel covered 74 events (82%). We re-sequenced 113 previously characterized clinical specimens at an average depth of 2,500X using SJPedPanel and recovered 354 (91%) of the 389 reported pathogenic variants. We then investigated the power of this panel in detecting mutations from specimens with low tumor purity (as low as 0.1%) using cell line-based dilution experiments and discovered that this gene panel enabled us to detect ∼80% variants with allele fraction of 0.2%, while the detection rate decreases to ∼50% when the allele fraction is 0.1%. We finally demonstrate its utility in disease monitoring on clinical specimens collected from AML patients in morphologic remission. Conclusions: SJPedPanel enables the detection of clinically relevant genetic alterations including rearrangements responsible for subtype-defining fusions for childhood cancers by targeted sequencing of ∼0.15% of human genome. It will enhance the analysis of specimens with low tumor burdens for cancer monitoring and early detection.

Malignant Progression of an Ancestral Bone Marrow Clone Harboring a CIC-NUTM2A Fusion in Isolated Myeloid Sarcoma.

Myeloid sarcoma is a rare condition consisting of extramedullary myeloid blasts found in association with acute myeloid leukemia or, in the absence of bone marrow involvement. We identified an infant with isolated myeloid sarcoma whose bone marrow was negative for involvement by flow cytometry. Sequencing revealed the fusion oncogene CIC-NUTM2A and identified the sarcoma to be clonally evolved from the bone marrow, which carried the fusion despite the absence of pathology. Murine modeling confirmed the ability of the fusion to transform hematopoietic cells and identified receptor tyrosine kinase (RTK) signaling activation consistent with disruption of the CIC transcriptional repressor. These findings extend the definition of CIC-rearranged malignancies to include hematologic disease, provide insight into the mechanism of oncogenesis, and demonstrate the importance of molecular analysis and tracking of bone marrow involvement over the course of treatment in myeloid sarcoma, including patients that lack flow cytometric evidence of leukemia at diagnosis. IMPLICATIONS: This study illustrates molecular involvement of phenotypically normal bone marrow in myeloid sarcoma, which has significant implications in clinical care. Further, it extends the definition of CIC-rearrangements to include hematologic malignancies and shows evidence of RTK activation that may be exploited therapeutically in cancer(s) driven by these fusions.

Induction, rapid fixation and retention of mutations in vegetatively propagated banana.

Mutation discovery technologies have enabled the development of reverse genetics for many plant species and allowed sophisticated evaluation of the consequences of mutagenesis. Such methods are relatively straightforward for seed-propagated plants. To develop a platform suitable for vegetatively propagated species, we treated isolated banana shoot apical meristems with the chemical mutagen ethyl methanesulphonate, recovered plantlets and screened for induced mutations. A high density of GC-AT transition mutations were recovered, similar to that reported in seed-propagated polyploids. Through analysis of the inheritance of mutations, we observed that genotypically heterogeneous stem cells resulting from mutagenic treatment are rapidly sorted to fix a single genotype in the meristem. Further, mutant genotypes are stably inherited in subsequent generations. Evaluation of natural nucleotide variation showed the accumulation of potentially deleterious heterozygous alleles, suggesting that mutation induction may uncover recessive traits. This work therefore provides genotypic insights into the fate of totipotent cells after mutagenesis and suggests rapid approaches for mutation-based functional genomics and improvement of vegetatively propagated crops.

Design and Implementation of Improved SARS-CoV-2 Diagnostic Assays To Mitigate the Impact of Genomic Mutations on Target Failure: the Xpert Xpress SARS-CoV-2 Experience.

In 2020, the U.S. Food and Drug Administration (FDA) enabled manufacturers to request emergency use authorization (EUA) to facilitate the rapid authorization of in vitro diagnostic (IVD) platforms for the detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Uncommon SARS-CoV-2 point mutations could cause nucleocapsid (N) gene target failure (NGTF) when using first-generation Xpert Xpress assays, so improvements were designed and implemented. In response to NGTF reports and with consideration of viral genomic information in public databases, the Xpress assays were redesigned to mitigate the impact of SARS-CoV-2 mutations on qualitative assay performance. The second-generation assays include a third gene target (RNA-dependent RNA polymerase [RdRp]) and redundant oligonucleotide probes for the N2 target. First- and second-generation assay performances were evaluated using a challenge set of samples. A second-generation assay with updated oligonucleotide chemistry received FDA EUA in September 2021. A prototype assay with oligonucleotide chemistry similar to that of the second-generation assay with FDA EUA successfully detected all three gene targets (N2, envelope [E], and RdRp) in all challenge samples (100%; 50/50), including variants with N gene mutations (g.29197C>T or g.29200C>T), which caused NGTF in the first-generation assays. Investigation and reporting of IVD target failures, public sharing of viral genomic sequence data, and the FDA EUA pathway were essential components in facilitating a short cycle time from the identification of mutations that impact the performance of an IVD assay to the design and implementation of an improved IVD assay. IMPORTANCE The SARS-CoV-2 genome has mutated during the coronavirus disease 2019 (COVID-19) pandemic. Some of these mutations have impacted the performance of nucleic acid amplification tests like PCR, which are commonly used as diagnostic tools to detect an infection. The U.S. Food and Drug Administration (FDA) emergency use authorization (EUA) process enables the rapid reformulation and regulatory authorization of improved PCRs. In our experience, the identification of SARS-CoV-2 mutations that impact PCR performance, the subsequent development of improved PCR chemistry, and the use of the FDA EUA regulatory pathway led to improved diagnostic performance during the SARS-CoV-2 pandemic that is able to keep pace with the rapidly evolving genome of SARS-CoV-2.

Antigen cross-presentation in young tumor-bearing hosts promotes CD8+ T cell terminal differentiation.

The immune system undergoes a progressive functional remodeling with age. Understanding how age bias shapes antitumor immunity is essential in designing effective immunotherapies, especially for pediatric patients. Here, we explore antitumor CD8+ T cell responses generated in young (prepubescent) and adult (presenescent) mice. Using an MHCI-deficient tumor model, we observed that tumor-reactive CD8+ T cells expanded in young tumor-bearing (TB) mice acquired a terminally differentiated phenotype characterized by overexpression of inhibitory receptors and the transcription factor Tox1. Furthermore, tumor-infiltrating CD8+ T cells from young tumors yielded a poor cytokine response compared with CD8+ T cells infiltrating adult tumors. Young migratory dendritic cells (migDCs) from the draining lymph nodes (dLNs), and mononuclear phagocytic cells (MPCs) infiltrating young tumors, were more competent in capturing and cross-presenting tumor antigen, leading to enhanced priming of CD8+ T cells in dLNs and their subsequent terminal differentiation in the tumors. Single-cell transcriptional profiling of tumor-infiltrating MPCs demonstrated that young MPCs are polarized toward an inflammatory, effector phenotype. Consistent with our observations in young versus adult TB mice, analysis of immune infiltrates from pediatric solid tumors showed a correlation between tumor-infiltrating CD8+ T cells with an exhaustion phenotype and the frequency of PD-L1-expressing monocytes/macrophages. Collectively, these data indicate that a young tissue microenvironment contributes to the generation of an immune response skewed toward a less pliable terminal effector state, thus narrowing the window for immunotherapeutic interventions.

SequencErr: measuring and suppressing sequencer errors in next-generation sequencing data.

BACKGROUND: There is currently no method to precisely measure the errors that occur in the sequencing instrument/sequencer, which is critical for next-generation sequencing applications aimed at discovering the genetic makeup of heterogeneous cellular populations. RESULTS: We propose a novel computational method, SequencErr, to address this challenge by measuring the base correspondence between overlapping regions in forward and reverse reads. An analysis of 3777 public datasets from 75 research institutions in 18 countries revealed the sequencer error rate to be ~ 10 per million (pm) and 1.4% of sequencers and 2.7% of flow cells have error rates > 100 pm. At the flow cell level, error rates are elevated in the bottom surfaces and > 90% of HiSeq and NovaSeq flow cells have at least one outlier error-prone tile. By sequencing a common DNA library on different sequencers, we demonstrate that sequencers with high error rates have reduced overall sequencing accuracy, and removal of outlier error-prone tiles improves sequencing accuracy. We demonstrate that SequencErr can reveal novel insights relative to the popular quality control method FastQC and achieve a 10-fold lower error rate than popular error correction methods including Lighter and Musket. CONCLUSIONS: Our study reveals novel insights into the nature of DNA sequencing errors incurred on DNA sequencers. Our method can be used to assess, calibrate, and monitor sequencer accuracy, and to computationally suppress sequencer errors in existing datasets.

Genomes for Kids: The Scope of Pathogenic Mutations in Pediatric Cancer Revealed by Comprehensive DNA and RNA Sequencing.

Genomic studies of pediatric cancer have primarily focused on specific tumor types or high-risk disease. Here, we used a three-platform sequencing approach, including whole-genome sequencing (WGS), whole-exome sequencing (WES), and RNA sequencing (RNA-seq), to examine tumor and germline genomes from 309 prospectively identified children with newly diagnosed (85%) or relapsed/refractory (15%) cancers, unselected for tumor type. Eighty-six percent of patients harbored diagnostic (53%), prognostic (57%), therapeutically relevant (25%), and/or cancer-predisposing (18%) variants. Inclusion of WGS enabled detection of activating gene fusions and enhancer hijacks (36% and 8% of tumors, respectively), small intragenic deletions (15% of tumors), and mutational signatures revealing of pathogenic variant effects. Evaluation of paired tumor-normal data revealed relevance to tumor development for 55% of pathogenic germline variants. This study demonstrates the power of a three-platform approach that incorporates WGS to interrogate and interpret the full range of genomic variants across newly diagnosed as well as relapsed/refractory pediatric cancers. SIGNIFICANCE: Pediatric cancers are driven by diverse genomic lesions, and sequencing has proven useful in evaluating high-risk and relapsed/refractory cases. We show that combined WGS, WES, and RNA-seq of tumor and paired normal tissues enables identification and characterization of genetic drivers across the full spectrum of pediatric cancers. This article is highlighted in the In This Issue feature, p. 2945.

Primary bone sarcoma with BCOR internal tandem duplication.

BCOR internal tandem duplications (ITDs) and rearrangements are implicated in the oncogenesis of a subset of undifferentiated sarcomas. To date, BCOR ITD sarcomas have been exclusively found in non-appendicular infantile soft tissues, whereas BCOR-rearranged sarcomas occur in both bones and soft tissues affecting a wider patient age range. Little is known about patient outcome in BCOR ITD sarcomas. We present a BCOR-expressing, primary bone, undifferentiated sarcoma case involving an adolescent male's left tibia that, unexpectedly, harbored a BCOR ITD instead of a BCOR rearrangement. Furthermore, the patient achieved a partial histologic response after receiving a Ewing sarcoma chemotherapy regimen. Our case expands the clinical spectrum of BCOR ITD sarcomas and suggests that childhood and adult BCOR-expressing sarcomas with an undifferentiated histology should be considered for both BCOR rearrangement and ITD screening. Accurate BCOR mutation identification in undifferentiated sarcomas is essential to define their clinical spectrum and to develop effective management strategies.

Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations.

Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes. To better understand neuroblastoma pathogenesis, here we analyze whole-genome, whole-exome and/or transcriptome data from 702 neuroblastoma samples. Forty percent of samples harbor at least one recurrent driver gene alteration and most aberrations, including MYCN, ATRX, and TERT alterations, differ in frequency by age. MYCN alterations occur at median 2.3 years of age, TERT at 3.8 years, and ATRX at 5.6 years. COSMIC mutational signature 18, previously associated with reactive oxygen species, is the most common cause of driver point mutations in neuroblastoma, including most ALK and Ras-activating variants. Signature 18 appears early and is continuous throughout disease evolution. Signature 18 is enriched in neuroblastomas with MYCN amplification, 17q gain, and increased expression of mitochondrial ribosome and electron transport-associated genes. Recurrent FGFR1 variants in six patients, and ALK N-terminal structural alterations in five samples, identify additional patients potentially amenable to precision therapy.

Loss of glucocorticoid receptor expression mediates in vivo dexamethasone resistance in T-cell acute lymphoblastic leukemia.

Despite decades of clinical use, mechanisms of glucocorticoid resistance are poorly understood. We treated primary murine T lineage acute lymphoblastic leukemias (T-ALLs) with the glucocorticoid dexamethasone (DEX) alone and in combination with the pan-PI3 kinase inhibitor GDC-0941 and observed a robust response to DEX that was modestly enhanced by GDC-0941. Continuous in vivo treatment invariably resulted in outgrowth of drug-resistant clones, ~30% of which showed markedly reduced glucocorticoid receptor (GR) protein expression. A similar proportion of relapsed human T-ALLs also exhibited low GR protein levels. De novo or preexisting mutations in the gene encoding GR (Nr3c1) occurred in relapsed clones derived from multiple independent parental leukemias. CRISPR/Cas9 gene editing confirmed that loss of GR expression confers DEX resistance. Exposing drug-sensitive T-ALLs to DEX in vivo altered transcript levels of multiple genes, and this response was attenuated in relapsed T-ALLs. These data implicate reduced GR protein expression as a frequent cause of glucocorticoid resistance in T-ALL.

Analysis of error profiles in deep next-generation sequencing data.

BACKGROUND: Sequencing errors are key confounding factors for detecting low-frequency genetic variants that are important for cancer molecular diagnosis, treatment, and surveillance using deep next-generation sequencing (NGS). However, there is a lack of comprehensive understanding of errors introduced at various steps of a conventional NGS workflow, such as sample handling, library preparation, PCR enrichment, and sequencing. In this study, we use current NGS technology to systematically investigate these questions. RESULTS: By evaluating read-specific error distributions, we discover that the substitution error rate can be computationally suppressed to 10-5 to 10-4, which is 10- to 100-fold lower than generally considered achievable (10-3) in the current literature. We then quantify substitution errors attributable to sample handling, library preparation, enrichment PCR, and sequencing by using multiple deep sequencing datasets. We find that error rates differ by nucleotide substitution types, ranging from 10-5 for A>C/T>G, C>A/G>T, and C>G/G>C changes to 10-4 for A>G/T>C changes. Furthermore, C>T/G>A errors exhibit strong sequence context dependency, sample-specific effects dominate elevated C>A/G>T errors, and target-enrichment PCR led to ~ 6-fold increase of overall error rate. We also find that more than 70% of hotspot variants can be detected at 0.1 ~ 0.01% frequency with the current NGS technology by applying in silico error suppression. CONCLUSIONS: We present the first comprehensive analysis of sequencing error sources in conventional NGS workflows. The error profiles revealed by our study highlight new directions for further improving NGS analysis accuracy both experimentally and computationally, ultimately enhancing the precision of deep sequencing.