Dataset Information

Influence of low tumor content on tumor mutational burden estimation by whole-exome sequencing and targeted panel sequencing.

ABSTRACT:

Background

Tumor mutational burden (TMB) is a promising biomarker for stratifying patient subpopulation who would benefit from immune checkpoint blockade (ICB) therapies. Although great efforts have been made for standardizing TMB measurement, mutation calling and TMB quantification can be challenging in samples with low tumor content including liquid biopsies. The effect of varying tumor content on TMB estimation by different assay methods has never been systematically investigated.

Method

We established a series of reference standard DNA samples derived from 11 pairs of tumor-normal matched human cell lines across different cancer types. Each tumor cell line was mixed with its matched normal at 0% (control), 1%, 2%, 5%, and 10% mass-to-mass ratio to mimic the clinical samples with low tumor content. TMB of these reference standards was evaluated by both ∼1000× whole-exome sequencing (wesTMB) and targeted panel sequencing (psTMB) at four different vendors. Both regression and classification analyses of TMB were performed for theoretical investigation and clinical practice purposes.

Results

Linear regression model was established that demonstrated in silico psTMB determined by regions of interest (ROI) as a great representative of wesTMB based on TCGA dataset. It was also true in our reference standard samples as the predicted psTMB interval based on the observed wesTMB captured the intended 90% of the in silico psTMB values. Although ∼1000× deep WES was applied, reference standard samples with less than 5% of tumor proportions are below the assay limit of detection (LoD) of wesTMB quantification. However, predicted wesTMB based on observed psTMB accurately classify (>0.97 AUC) for TMB high and low patient stratification even in samples with 2% of tumor content, which is more clinically relevant, as TMB determination should be a qualitative assay for TMB high and low patient classification. One targeted panel sequencing vendor using an optimized blood psTMB pipeline can further classify TMB status accurately (>0.82 AUC) in samples with only 1% of tumor content.

Conclusions

We developed a linear model to establish the quantitative correlation between wesTMB and psTMB. A set of DNA reference standards was produced in aid to standardize TMB measurements in samples with low tumor content across different targeted sequencing panels. This study is a significant contribution aiming to harmonize TMB estimation and extend its future application in clinical samples with low tumor content including liquid biopsy.

SUBMITTER: Zhang W

PROVIDER: S-EPMC8102856 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BACKGROUND:Tumor mutational burden (TMB) measurement is limited by low tumor purity of samples, which can influence prediction of the immunotherapy response, particularly when using whole-exome sequencing-based TMB (wTMB). This issue could be overcome by targeted panel sequencing-based TMB (pTMB) with higher depth of coverage, which remains unexplored. METHODS:We comprehensively reanalyzed four public datasets of immune checkpoint inhibitor (ICI)-treated cohorts (adopting pTMB or wTMB) to test each biomarker's predictive ability for low purity samples (cut-off: 30%). For validation, paired genomic profiling with the same tumor specimens was performed to directly compare wTMB and pTMB in patients with breast cancer (paired-BRCA, n=165) and ICI-treated patients with advanced non-small-cell lung cancer (paired-NSCLC, n=156). RESULTS:Low tumor purity was common (range 30%-45%) in real-world samples from ICI-treated patients. In the survival analyzes of public cohorts, wTMB could not predict the clinical benefit of immunotherapy when tumor purity was low (log-rank p=0.874), whereas pTMB could effectively stratify the survival outcome (log-rank p=0.020). In the paired-BRCA and paired-NSCLC cohorts, pTMB was less affected by tumor purity, with significantly more somatic variants identified at low allele frequency (p<0.001). We found that wTMB was significantly underestimated in low purity samples with a large proportion of clonal variants undetected by whole-exome sequencing. Interestingly, pTMB more accurately predicted progression-free survival (PFS) after immunotherapy than wTMB owing to its superior performance in the low tumor purity subgroup (p=0.054 vs p=0.358). Multivariate analysis revealed pTMB (p=0.016), but not wTMB (p=0.32), as an independent predictor of PFS even in low-purity samples. The net reclassification index using pTMB was 21.7% in the low-purity subgroup (p=0.016). CONCLUSIONS:Our data suggest that TMB characterization with targeted deep sequencing might have potential strength in predicting ICI responsiveness due to its enhanced sensitivity for hard-to-detect variants at low-allele fraction. Therefore, pTMB could act as an invaluable biomarker in the setting of both clinical trials and practice outside of trials based on its reliable performance in mitigating the purity-related bias.

Project description:BACKGROUND:Europeans and American Indians were major genetic ancestry of Hispanics in the U.S. These ancestral groups have markedly different incidence rates and outcomes in many types of cancers. Therefore, the genetic admixture may cause biased genetic association study with cancer susceptibility variants specifically in Hispanics. For example, the incidence rate of liver cancer has been shown with substantial disparity between Hispanic, Asian and non-Hispanic white populations. Currently, ancestry informative marker (AIM) panels have been widely utilized with up to a few hundred ancestry-informative single nucleotide polymorphisms (SNPs) to infer ancestry admixture. Notably, current available AIMs are predominantly located in intron and intergenic regions, while the whole exome sequencing (WES) protocols commonly used in translational research and clinical practice do not cover these markers. Thus, it remains challenging to accurately determine a patient's admixture proportion without additional DNA testing. RESULTS:In this study we designed an unique AIM panel that infers 3-way genetic admixture from three distinct and selective continental populations (African (AFR), European (EUR), and East Asian (EAS)) within evolutionarily conserved exonic regions. Initially, about 1 million exonic SNPs from selective three populations in the 1000 Genomes Project were trimmed by their linkage disequilibrium (LD), restricted to biallelic variants, and finally we optimized to an AIM panel with 250 SNP markers, or the UT-AIM250 panel, using their ancestral informativeness statistics. Comparing to published AIM panels, UT-AIM250 performed better accuracy when we tested with three ancestral populations (accuracy: 0.995?±?0.012 for AFR, 0.997?±?0.007 for EUR, and 0.994?±?0.012 for EAS). We further demonstrated the performance of the UT-AIM250 panel to admixed American (AMR) samples of the 1000 Genomes Project and obtained similar results (AFR, 0.085?±?0.098; EUR, 0.665?±?0.182; and EAS, 0.250?±?0.205) to previously published AIM panels (Phillips-AIM34: AFR, 0.096?±?0.127, EUR, 0.575?±?0.290, and EAS, 0.330?±?0.315; Wei-AIM278: AFR, 0.070?±?0.096, EUR, 0.537?±?0.267, and EAS, 0.393?±?0.300). Subsequently, we applied the UT-AIM250 panel to a clinical dataset of 26 self-reported Hispanic patients in South Texas with hepatocellular carcinoma (HCC). We estimated the admixture proportions using WES data of adjacent non-cancer liver tissues (AFR, 0.065?±?0.043; EUR, 0.594?±?0.150; and EAS, 0.341?±?0.160). Similar admixture proportions were identified from corresponding tumor tissues. In addition, we estimated admixture proportions of The Cancer Genome Atlas (TCGA) collection of hepatocellular carcinoma (TCGA-LIHC) samples (376 patients) using the UT-AIM250 panel. The panel obtained consistent admixture proportions from tumor and matched normal tissues, identified 3 possible incorrectly reported race/ethnicity, and/or provided race/ethnicity determination if necessary. CONCLUSIONS:Here we demonstrated the feasibility of using evolutionarily conserved exonic regions to infer admixture proportions and provided a robust and reliable control for sample collection or patient stratification for genetic analysis. R implementation of UT-AIM250 is available at https://github.com/chenlabgccri/UT-AIM250.