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

Meredith Yeager and colleagues with the Cancer Genetics Markers of Susceptibility (CGEMS) initiative report a new association to prostate cancer at chromosome 8q24. This defines a new locus, region 4, which shows association to prostate cancer susceptibility independent of previously reported associations at 8q24. We report a genome-wide association study in 10,286 cases and 9,135 controls of European ancestry in the Cancer Genetic Markers of Susceptibility (CGEMS) initiative. We identify a new association with prostate cancer risk on chromosome 8q24 (rs620861, P = 1.3 × 10−10, heterozygote OR = 1.17, 95% CI 1.10–1.24; homozygote OR = 1.33, 95% CI 1.21–1.45). This defines a new locus associated with prostate cancer susceptibility on 8q24.

Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (P_interaction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

Mounting evidence supports the occurrence of accelerating aging among long-term survivors of childhood cancer. We aimed to investigate epigenetic age acceleration (EAA) in survivors and evaluate associations between EAA, treatment exposures, health behaviors, and chronic health conditions (CHCs). Genome-wide methylation data were generated with Infinium EPIC BeadChip on blood-derived DNA from 2139 survivors and 282 frequency matched controls from the St Jude Lifetime Cohort Study. EAAs were estimated as residuals from a linear regression of epigenetic age (Levine’s clock) against chronological age. Adjusted least square mean (ALSM) of EAA was calculated and compared between survivors and controls, across treatment exposures and health behaviors. Associations of EAA with 20 clinically assessed CHCs were evaluated with multivariable piecewise-exponential models. All statistical tests for P values below were 2-sided. EAA was statistically significantly higher in survivors than controls (ALSM = 0.63, 95% confidence interval [CI] = 0.26 to 1.01 vs −3.61, 95% CI = −4.43 to 2.80). In a multivariable model among survivors, statistically significantly higher EAA (P < .05) was observed in those exposed to chest radiotherapy, abdomen or pelvic radiotherapy, alkylating agents, glucocorticoids, or epipodophyllotoxins. Compared with survivors with favorable health behaviors (ALSM = 0.26, 95% CI=−0.36 to 0.87), EAA was statistically significantly higher among survivors with intermediate (ALSM = 1.07, 95% CI = 0.59 to 1.54) or unfavorable health behaviors (ALSM = 1.45, 95% CI = 0.60 to 2.30). In time-to-event analyses, statistically significant associations were identified between EAA tertiles and incidence of 7 CHCs: hypertension (3rd vs 1st tertile, relative rate [RR] = 1.83, 95% CI = 1.17 to 2.83), myocardial infarction (RR = 2.91, 95% CI = 1.27 to 7.21), obesity (RR = 1.39, 95% CI = 1.17 to 1.66), obstructive pulmonary deficit (RR = 1.86, 95% CI = 0.95 to 3.77), peripheral motor neuropathy (RR = 2.89, 95% CI = 1.24 to 6.97), peripheral sensory neuropathy (RR = 2.04, 95% CI = 0.99 to 4.26), and pulmonary diffusion deficits (RR = 2.75, 95% CI = 0.95 to 7.63). EAA is statistically significantly higher in survivors of childhood cancer than in noncancer controls and is associated with specific treatment exposures, unfavorable health behaviors, and presence of specific CHCs.

We present the first comprehensive investigation of clonal hematopoiesis (CH) in 2,860 long-term survivors of pediatric cancer with a median follow-up time of 23.5 years. Deep sequencing over 39 CH-related genes reveals mutations in 15% of the survivors, significantly higher than the 8.5% in 324 community controls. CH in survivors is associated with exposures to alkylating agents, radiation, and bleomycin. Therapy-related CH shows significant enrichment in STAT3, characterized as a CH gene specific to survivors of Hodgkin lymphoma, and TP53. Single-cell profiling of peripheral blood samples revealed STAT3 mutations predominantly present in T cells and contributed by SBS25, a mutational signature associated with procarbazine exposure. Serial sample tracking reveals that larger clone size is a predictor for future expansion of age-related CH clones, whereas therapy-related CH remains stable decades after treatment. These data depict the distinct dynamics of these CH subtypes and support the need for longitudinal monitoring to determine the potential contribution to late effects. Significance: This first comprehensive CH analysis in long-term survivors of pediatric cancer presents the elevated prevalence and therapy exposures/diagnostic spectrum associated with CH. Due to the contrasting dynamics of clonal expansion for age-related versus therapy-related CH, longitudinal monitoring is recommended to ascertain the long-term effects of therapy-induced CH in pediatric cancer survivors. See related commentary by Collord and Behjati, p. 811. This article is highlighted in the In This Issue feature, p. 799

PURPOSE: To investigate cancer treatment plus pathogenic germline mutations (PGMs) in DNA repair genes (DRGs) for identification of childhood cancer survivors at increased risk of subsequent neoplasms (SNs). METHODS: Whole-genome sequencing was performed on blood-derived DNA from survivors in the St Jude Lifetime Cohort. PGMs were evaluated in 127 genes from 6 major DNA repair pathways. Cumulative doses of chemotherapy and body region–specific radiotherapy (RT) were abstracted from medical records. Relative rates (RRs) and 95% CIs of SNs by mutation status were estimated using multivariable piecewise exponential models. RESULTS: Of 4,402 survivors, 495 (11.2%) developed 1,269 SNs. We identified 538 PGMs in 98 DRGs ( POLG, MUTYH, ERCC2, and BRCA2, among others) in 508 (11.5%) survivors. Mutations in homologous recombination (HR) genes were significantly associated with an increased rate of subsequent female breast cancer (RR, 3.7; 95% CI, 1.8 to 7.7), especially among survivors with chest RT ≥ 20 Gy (RR, 4.4; 95% CI, 1.6 to 12.4), or with a cumulative dose of anthracyclines in the second or third tertile (RR, 4.4; 95% CI, 1.7 to 11.4). Mutations in HR genes were also associated with an increased rate of subsequent sarcoma among those who received alkylating agent doses in the third tertile (RR, 14.9; 95% CI, 4.0 to 38.0). Mutations in nucleotide excision repair genes were associated with subsequent thyroid cancer for those treated with neck RT ≥ 30 Gy (RR, 12.9; 95% CI, 1.6 to 46.6) with marginal statistical significance. CONCLUSION: Our study provides novel insights regarding the contribution of genetics, in combination with known treatment-related risks, for the development of SNs. These findings have the potential to facilitate identification of high-risk survivors who may benefit from genetic counseling and/or testing of DRGs, which may further inform personalized cancer surveillance and prevention strategies.