Tao Wu

Copy number alteration (CNA) is a major type of cancer genome alteration that drives cancer progression. CNA signature analysis can reveal underlying etiology and provide biomarkers for cancer treatment, and existing CNA signature analyses are all performed using bulk tissue samples. However, CNA usually affects a large proportion of genome, and the CNA profile of bulk sample does not reflect the actual CNA profiles of the individual cancer cells of the sample, especially in tumors with high heterogeneity, such as hepatocellular carcinoma (HCC). Furthermore, the evolutionary trajectory of CNA mutational processes still remains elusive. Here we build a method to comprehensively analyze the CNA signatures of HCC from single-cell and bulk sample perspective, revealing patterns and potential noise signals from the usually performed bulk tissue CNA signature analysis. Single-cell signature analysis delineates the evolutionary trajectory of HCC CNA signatures, and different CNA signatures consistently emerge in different HCC evolution stages. Single-cell CNA signatures show robust performance in patient prognosis and drug sensitivity prediction. This work not only reveals specific considerations in analyzing CNA signature derived from bulk tissue but also depicts CNA evolution process and provides potential biomarkers for the prognosis and treatment of HCC patients.

Low-coverage whole-genome sequencing (lcWGS) provides a cost-effective method for genome-wide copy number variation (CNV) profiling, yet its technical limitations and analytical variability require systematic evaluation. We benchmarked five CNV detection tools using simulated and real-world datasets, focusing on sequencing depth, formalin-fixed paraffin-embedded (FFPE) artifacts, tumor purity, multi-center reproducibility, and signature-level stability. Our results demonstrate that ichorCNA outperformed other tools in precision and runtime at high purity (≥50%), making it the optimal choice for lcWGS-based workflows. Prolonged FFPE fixation induced artifactual short-segment CNVs due to formalin-driven DNA fragmentation, a bias none of the tools could computationally correct, necessitating strict fixation time control or prioritization of fresh-frozen samples. Multi-center analysis revealed high reproducibility for the same tool across sequencing facilities, but comparisons between different tools showed low concordance. Copy number features extracted by the Wang et al. method exhibited superior stability across conditions compared with the Steele et al. method and the Tao et al. method. This study establishes actionable guidelines for lcWGS: prioritize ichorCNA (ensuring ≥50% tumor purity), optimize FFPE protocol, and use Wang et al. features to ensure robust copy number profiling in precision oncology.

Copy number alteration (CNA) is a major type of cancer genome alteration that drives cancer progression. CNA signature analysis can reveal underlying etiology and provide biomarkers for cancer treatment, and existing CNA signature analyzes are all performed using bulk tissue samples. However CNA usually affect large proportion of genome, and the CNA profile of bulk sample does not reflect the actual CNA profiles of the individual cancer cells of the sample, especially in tumors with high heterogeneity, such as hepatocellular carcinoma (HCC). Furthermore, the evolutionary trajectory of CNA mutational processes still remain elusive. Here we build a method to comprehensively analyze the CNA signatures of HCC from single-cell and bulk sample perspective, revealing patterns and potential noise signals from the usually performed bulk tissue CNA signature analysis. Single-cell signature analysis delineated the evolutionary trajectory of HCC CNA signatures, and different CNA signatures consistently emerge in different HCC evolution stages. Single-cell CNA signatures show robust performance in patient prognosis and drug sensitivity prediction. This work not only reveals specific considerations in analyzing CNA signature derived from bulk tissue but also depicts CNA evolution process and provides potential biomarkers for the prognosis and treatment of HCC patients.

Homologous recombination deficiency (HRD) renders cancer cells vulnerable to unrepaired double-strand breaks and is an important therapeutic target as exemplified by the clinical efficacy of poly ADP-ribose polymerase (PARP) inhibitors as well as the platinum chemotherapy drugs applied to HRD patients. However, it remains a challenge to predict HRD status precisely and economically. Copy number alteration (CNA), as a pervasive trait of human cancers, can be extracted from a variety of data sources, including whole genome sequencing (WGS), SNP array, and panel sequencing, and thus can be easily applied clinically. Here we systematically evaluate the predictive performance of various CNA features and signatures in HRD prediction and build a gradient boosting machine model (HRD_CNA) for pan-cancer HRD prediction based on these CNA features. CNA features BP10MB[1] (The number of breakpoints per 10MB of DNA is 1) and SS[ > 7 & <=8] (The log10-based size of segments is greater than 7 and less than or equal to 8) are identified as the most important features in HRD prediction. HRD_CNAsuggests the biallelic inactivation ofBRCA1,BRCA2,PALB2,RAD51C,RAD51D, andBARD1as the major genetic basis for human HRD, and may also be applied to effectively validate the pathogenicity of BRCA1/2 variants of uncertain significance (VUS). Together, this study provides a robust tool for cost-effective HRD prediction and also demonstrates the applicability of CNA features and signatures in cancer precision medicine.

DNA methylation analysis has been applied to determine the primary site of cancer; however, robust and accurate prediction of cancer types with a minimum number of sites is still a significant scientific challenge. To build an accurate and robust cancer type prediction tool with a minimum number of DNA methylation sites, we internally benchmarked different DNA methylation site selection and ranking procedures, as well as different classification models. We used The Cancer Genome Atlas dataset (26 cancer types with 8296 samples) to train and test models and used an independent dataset (17 cancer types with 2738 samples) for model validation. A deep neural network model using a combined feature selection procedure (named MethyDeep) can predict 26 cancer types using 30 methylation sites with superior performance compared with the known methods for both primary and metastatic cancers in independent validation datasets. In conclusion, MethyDeep is an accurate and robust cancer type predictor with the minimum number of DNA methylation sites; it could help the cost-effective clarification of cancer of unknown primary patients and the liquid biopsy-based early screening of cancers.

Immunoediting includes three temporally distinct stages, termed elimination, equilibrium, and escape, and has been proposed to explain the interactions between cancer cells and the immune system during the evolution of cancer. However, the status of immunoediting in cancer remains unclear, and the existence of neoantigen depletion in untreated cancer has been debated. Here we developed a distribution pattern-based method for quantifying neoantigen-mediated negative selection in cancer evolution. The method can provide a robust and reliable quantification for immunoediting signal in individual patients with cancer. Moreover, this method demonstrated the prevalence of immunoediting in the immunotherapy-naive cancer genome. The elimination and escape stages of immunoediting can be quantified separately, where tumor types with strong immunoediting-elimination exhibit a weak immunoediting-escape signal, and vice versa. The quantified immunoediting-elimination signal was predictive of clinical response to cancer immunotherapy. Collectively, immunoediting quantification provides an evolutionary perspective for evaluating the antigenicity of neoantigens and reveals a potential biomarker for precision immunotherapy in cancer.Significance:Quantification of neoantigen-mediated negative selection in cancer progression reveals distinct features of cancer immunoediting and can serve as a potential biomarker to predict immunotherapy response.

Supplementary Figure from Quantification of Neoantigen-Mediated Immunoediting in Cancer Evolution

Supplementary Figure from Quantification of Neoantigen-Mediated Immunoediting in Cancer Evolution

Immunoediting includes three temporally distinct stages, termed elimination, equilibrium, and escape, and has been proposed to explain the interactions between cancer cells and the immune system during the evolution of cancer. However, the status of immunoediting in cancer remains unclear, and the existence of neoantigen depletion in untreated cancer has been debated. Here we developed a distribution pattern-based method for quantifying neoantigen-mediated negative selection in cancer evolution. The method can provide a robust and reliable quantification for immunoediting signal in individual patients with cancer. Moreover, this method demonstrated the prevalence of immunoediting in the immunotherapy-naive cancer genome. The elimination and escape stages of immunoediting can be quantified separately, where tumor types with strong immunoediting-elimination exhibit a weak immunoediting-escape signal, and vice versa. The quantified immunoediting-elimination signal was predictive of clinical response to cancer immunotherapy. Collectively, immunoediting quantification provides an evolutionary perspective for evaluating the antigenicity of neoantigens and reveals a potential biomarker for precision immunotherapy in cancer.Significance:Quantification of neoantigen-mediated negative selection in cancer progression reveals distinct features of cancer immunoediting and can serve as a potential biomarker to predict immunotherapy response.