Structural variations (SVs) are wide-spread in human genomes and may have important implications in disease-related and evolutionary studies. High-throughput sequencing (HTS) has become a major platform for SV detection and simulation serves as a powerful and cost-effective approach for benchmarking SV detection algorithms. Accurate performance assessment by simulation requires the simulator capable of generating simulation data with all important features of real data, such GC biases in HTS data and various complexities in tumor data. However, no available package has systematically addressed all issues in data simulation for SV benchmarking.Pysim-sv is a package for simulating HTS data to evaluate performance of SV detection algorithms. Pysim-sv can introduce a wide spectrum of germline and somatic genomic variations. The package contains functionalities to simulate tumor data with aneuploidy and heterogeneous subclones, which is very useful in assessing algorithm performance in tumor studies. Furthermore, Pysim-sv can introduce GC-bias, the most important and prevalent bias in HTS data, in the simulated HTS data.Pysim-sv provides an unbiased toolkit for evaluating HTS-based SV detection algorithms.
Hepatocellular carcinoma (HCC) is a cancer of substantial morphologic, genetic and phenotypic diversity. Yet we do not understand the relationship between intratumor heterogeneity and the associated morphologic/histological characteristics of the tumor. Using single-cell whole-genome sequencing to profile 96 tumor cells (30-36 each) and 15 normal liver cells (5 each), collected from three male patients with HBV-associated HCC, we confirmed that copy number variations occur early in hepatocarcinogenesis but thereafter remain relatively stable throughout tumor progression. Importantly, we showed that specific HCCs can be of monoclonal or polyclonal origins. Tumors with confluent multinodular morphology are the typical polyclonal tumors and display the highest intratumor heterogeneity. In addition to mutational and copy number profiles, we dissected the clonal origins of HCC using HBV-derived foreign genomic markers. In monoclonal HCC, all the tumor single cells exhibit the same HBV integrations, indicating that HBV integration is an early driver event and remains extremely stable during tumor progression. In addition, our results indicated that both models of metastasis, late dissemination and early seeding, have a role in HCC progression. Notably, early intrahepatic spreading of the initiating clone leads to the formation of synchronous multifocal tumors. Meanwhile, we identified a potential driver gene ZNF717 in HCC, which exhibits a high frequency of mutation at both single-cell and population levels, as a tumor suppressor acting through regulating the IL-6/STAT3 pathway. These findings highlight multiple distinct tumor evolutionary mechanisms in HCC, which suggests the need for specific treatment strategies.
Abstract Mutation signature analysis has been used to infer the contributions of various DNA mutagenic-repair events in individual cancer genomes. Here, we build a statistical framework using a multinomial distribution to assign individual mutations to their cognate mutation signatures. We applied it to 47 million somatic mutations in 1925 publicly available cancer genomes to obtain a mutation signature map at the resolution of individual somatic mutations. Based on mutation signature-level genetic-epigenetic correlative analyses, mutations with transcriptional and replicative strand asymmetries show different enrichment patterns across genomes, and “transcribed” chromatin states and gene boundaries are particularly vulnerable to transcription-coupled repair activities. While causative processes of cancer-driving mutations can be diverse, as shown for converging effects of multiple mutational processes on TP53 mutations, the substantial fraction of recurrently mutated amino acids points to specific mutational processes, e.g., age-related C-to-T transition for KRAS p.G12 mutations. Our investigation of evolutionary trajectories with respect to mutation signatures further revealed that candidate pairs of early- vs. late-operative mutation processes in cancer genomes represent evolutionary dynamics of multiple mutational processes in the shaping of cancer genomes. We also observed that the local mutation clusters of kataegis often include mutations arising from multiple mutational processes, suggestive of a locally synchronous impact of multiple mutational processes on cancer genomes. Taken together, our examination of the genome-wide landscape of mutation signatures at the resolution of individual somatic mutations shows the spatially and temporally distinct mutagenesis-repair-replication histories of various mutational processes and their effects on shaping cancer genomes.
The inheritance of recurrent patellar dislocation (RPD) is known, but the susceptible gene remains unidentified. Here, we performed the first whole exome sequencing (WES) cohort study to identify the susceptible genes. The results showed eight genes were associated with this disease. Notably, the carboxypeptidase D (CPD) gene showed the highest relevance based on its gene function and tissue expression. Single-cell sequencing results indicate that the CPD gene is involved in the pathophysiological process of RPD through granulocytes. Implicated pathways include nuclear factor kappa B (NF-κB), mitogen-activated protein kinase (MAPK), and Wnt/β-catenin signaling, potentially influencing CPD's role in RPD pathogenesis. This study identified the susceptible gene and investigates the potential pathogenesis of RPD, which provided a new prospect for the understanding of RPD. Besides, it would offer the theoretical basis for disease prevention and genetic counseling.
We performed an integrated genomic, transcriptomic and proteomic characterization of 373 endometrial carcinomas using array- and sequencing-based technologies. Uterine serous tumours and ∼25% of high-grade endometrioid tumours had extensive copy number alterations, few DNA methylation changes, low oestrogen receptor/progesterone receptor levels, and frequent TP53 mutations. Most endometrioid tumours had few copy number alterations or TP53 mutations, but frequent mutations in PTEN, CTNNB1, PIK3CA, ARID1A and KRAS and novel mutations in the SWI/SNF chromatin remodelling complex gene ARID5B. A subset of endometrioid tumours that we identified had a markedly increased transversion mutation frequency and newly identified hotspot mutations in POLE. Our results classified endometrial cancers into four categories: POLE ultramutated, microsatellite instability hypermutated, copy-number low, and copy-number high. Uterine serous carcinomas share genomic features with ovarian serous and basal-like breast carcinomas. We demonstrated that the genomic features of endometrial carcinomas permit a reclassification that may affect post-surgical adjuvant treatment for women with aggressive tumours. An integrative genomic analysis of several hundred endometrial carcinomas shows that a minority of tumour samples carry copy number alterations or TP53 mutations and many contain key cancer-related gene mutations, such as those involved in canonical pathways and chromatin remodelling; a reclassification of endometrial tumours into four distinct types is proposed, which may have an effect on patient treatment regimes. This paper from The Cancer Genome Atlas Research Network presents an in-depth genome-wide analysis of endometrial (uterine) carcinomas from more than 350 patients. Based on a series of genomic features including newly identified hotspot mutations in the DNA polymerase gene POLE, and novel mutations in the ARID5B DNA-binding protein, the authors propose a reclassification of endometrial tumours into four distinct types. This might have clinical relevance for post-surgical adjuvant treatment of women with aggressive tumours.
Abstract Liver metastasis, the leading cause of colorectal cancer mortality, exhibits a highly heterogeneous and suppressive immune microenvironment. Here, we sequenced 97 matched samples by using single-cell RNA sequencing and spatial transcriptomics. Strikingly, the metastatic microenvironment underwent remarkable spatial reprogramming of immunosuppressive cells such as MRC1+ CCL18+ M2-like macrophages. We further developed scMetabolism, a computational pipeline for quantifying single-cell metabolism, and observed that those macrophages harbored enhanced metabolic activity. Interestingly, neoadjuvant chemotherapy could block this status and restore the antitumor immune balance in responsive patients, whereas the nonresponsive patients deteriorated into a more suppressive one. Our work described the immune evolution of metastasis and uncovered the black box of how tumors respond to neoadjuvant chemotherapy. Significance: We present a single-cell and spatial atlas of colorectal liver metastasis and found the highly metabolically activated MRC1+ CCL18+ M2-like macrophages in metastatic sites. Efficient neoadjuvant chemotherapy can slow down such metabolic activation, raising the possibility to target metabolism pathways in metastasis. This article is highlighted in the In This Issue feature, p. 1
'/%3e%3cuse xlink:href='%23a' transform='rotate(23.036 .093 25.536)'/%3e%3cuse xlink:href='%23a' transform='rotate(45.87 1.273 16.18)'/%3e%3cuse xlink:href='%23a' transform='rotate(69.945 .996 12.078)'/%3e%3cuse xlink:href='%23a' transform='rotate(20.66 -19.689 31.932)'/%3e%3c/svg%3e)
