Supervised mutational signatures for obesity and other tissue-specific etiological factors in cancer

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

Identifying somatic fingerprints of cancers defined by germline and environmental risk factors

Numerous studies over the past generation have identified germline variants that increase specific cancer risks. Simultaneously, a revolution in sequencing technology has permitted high-throughput annotations of somatic genomes characterizing individual tumors. However, examining the relationship between germline variants and somatic alteration patterns is hugely challenged by the large numbers of variants in a typical tumor, the rarity of most individual variants, and the heterogeneity of tumor somatic fingerprints. In this article, we propose statistical methodology that frames the investigation of germline-somatic relationships in an interpretable manner. The method uses meta-features embodying biological contexts of individual somatic alterations to implicitly group rare mutations. Our team has used this technique previously through a multilevel regression model to diagnose with high accuracy tumor site of origin. Herein, we further leverage topic models from computational linguistics to achieve interpretable lower-dimensional embeddings of the meta-features. We demonstrate how the method can identify distinctive somatic profiles linked to specific germline variants or environmental risk factors. We illustrate the method using The Cancer Genome Atlas whole-exome sequencing data to characterize somatic tumor fingerprints in breast cancer patients with germline BRCA1/2 mutations and in head and neck cancer patients exposed to human papillomavirus.

Accurate and sensitive mutational signature analysis with MuSiCal

Mutational signature analysis is a recent computational approach for interpreting somatic mutations in the genome. Its application to cancer data has enhanced our understanding of mutational forces driving tumorigenesis and demonstrated its potential to inform prognosis and treatment decisions. However, methodological challenges remain for discovering new signatures and assigning proper weights to existing signatures, thereby hindering broader clinical applications. Here we present Mutational Signature Calculator (MuSiCal), a rigorous analytical framework with algorithms that solve major problems in the standard workflow. Our simulation studies demonstrate that MuSiCal outperforms state-of-the-art algorithms for both signature discovery and assignment. By reanalyzing more than 2,700 cancer genomes, we provide an improved catalog of signatures and their assignments, discover nine indel signatures absent in the current catalog, resolve long-standing issues with the ambiguous 'flat' signatures and give insights into signatures with unknown etiologies. We expect MuSiCal and the improved catalog to be a step towards establishing best practices for mutational signature analysis.

High-throughput sequencing and in-silico analysis confirm pathogenicity of novel MSH3 variants in African American colorectal cancer

The maintenance of DNA sequence integrity is critical to avoid accumulation of cancer-causing mutations. Inactivation of DNA Mismatch Repair (MMR) genes (e.g., MLH1 and MSH2) is common among many cancers, including colorectal cancer (CRC) and is the driver of classic microsatellite instability (MSI) in tumors. Somatic MSH3 alterations have been linked to a specific form of MSI called elevated microsatellite alterations at selected tetranucleotide repeats (EMAST) that is associated with patient poor prognosis and elevated among African American (AA) rectal cancer patients. Genetic variants of MSH3 and their pathogenicity vary among different populations, such as among AA, which are not well-represented in publicly available databases. Targeted exome sequencing of MSH3 among AA CRC samples followed by computational bioinformatic pipeline and molecular dynamic simulation analysis approach confirmed six identified MSH3 variants (c.G1237A, c.C2759T, c.G1397A, c.G2926A, c.C3028T, c.G3241A) that corresponded to MSH3 amino-acid changes (p.E413K; p.S466N; p.S920F; p.E976K; p.H1010Y; p.E1081K). All identified MSH3 variants were non-synonymous, novel, pathogenic, and show loss or gain of hydrogen bonding, ionic bonding, hydrophobic bonding, and disulfide bonding and have a deleterious effect on the structure of MSH3 protein. Some variants were located within the ATPase site of MSH3, affecting ATP hydrolysis that is critical for MSH3's function. Other variants were in the MSH3-MSH2 interacting domain, important for MSH3's binding to MSH2. Overall, our data suggest that these variants among AA CRC patients affect the function of MSH3 making them pathogenic and likely contributing to the development or advancement of CRC among AA. Further clarifying functional studies will be necessary to fully understand the impact of these variants on MSH3 function and CRC development in AA patients.

Obesity-associated changes in molecular biology of primary breast cancer

Obesity is associated with an increased risk of developing breast cancer (BC) and worse prognosis in BC patients, yet its impact on BC biology remains understudied in humans. This study investigates how the biology of untreated primary BC differs according to patients' body mass index (BMI) using data from >2,000 patients. We identify several genomic alterations that are differentially prevalent in overweight or obese patients compared to lean patients. We report evidence supporting an ageing accelerating effect of obesity at the genetic level. We show that BMI-associated differences in bulk transcriptomic profile are subtle, while single cell profiling allows detection of more pronounced changes in different cell compartments. These analyses further reveal an elevated and unresolved inflammation of the BC tumor microenvironment associated with obesity, with distinct characteristics contingent on the estrogen receptor status. Collectively, our analyses imply that obesity is associated with an inflammaging-like phenotype. We conclude that patient adiposity may play a significant role in the heterogeneity of BC and should be considered for BC treatment tailoring.

Evaluating cancer etiology and risk with a mathematical model of tumor evolution

Recent evidence arising from DNA sequencing of healthy human tissues has clearly indicated that our organs accumulate a relevant number of somatic mutations due to normal endogenous mutational processes, in addition to those caused by environmental factors. A deeper understanding of the evolution of this endogenous mutational load is critical for understanding what causes cancer. Here we present a mathematical model of tumor evolution that is able to predict the expected number of endogenous somatic mutations present in various tissue types of a patient at a given age. These predictions are then compared to those observed in patients. We also obtain an improved fitting of the variation in cancer incidence across cancer types, showing that the endogenous mutational processes can explain 4/5 of the variation in cancer risk. Overall, these results offer key insights into cancer etiology, by providing further evidence for the major role these endogenous processes play in cancer.

Theory and Practice of Integrating Machine Learning and Conventional Statistics in Medical Data Analysis

The practice of medical decision making is changing rapidly with the development of innovative computing technologies. The growing interest of data analysis with improvements in big data computer processing methods raises the question of whether machine learning can be integrated with conventional statistics in health research. To help address this knowledge gap, this paper presents a review on the conceptual integration between conventional statistics and machine learning, focusing on the health research. The similarities and differences between the two are compared using mathematical concepts and algorithms. The comparison between conventional statistics and machine learning methods indicates that conventional statistics are the fundamental basis of machine learning, where the black box algorithms are derived from basic mathematics, but are advanced in terms of automated analysis, handling big data and providing interactive visualizations. While the nature of both these methods are different, they are conceptually similar. Based on our review, we conclude that conventional statistics and machine learning are best to be integrated to develop automated data analysis tools. We also strongly believe that machine learning could be explored by health researchers to enhance conventional statistics in decision making for added reliable validation measures.

Genome-wide mutational signatures in low-coverage whole genome sequencing of cell-free DNA

Mutational signatures accumulate in somatic cells as an admixture of endogenous and exogenous processes that occur during an individual’s lifetime. Since dividing cells release cell-free DNA (cfDNA) fragments into the circulation, we hypothesize that plasma cfDNA might reflect mutational signatures. Point mutations in plasma whole genome sequencing (WGS) are challenging to identify through conventional mutation calling due to low sequencing coverage and low mutant allele fractions. In this proof of concept study of plasma WGS at 0.3–1.5x coverage from 215 patients and 227 healthy individuals, we show that both pathological and physiological mutational signatures may be identified in plasma. By applying machine learning to mutation profiles, patients with stage I-IV cancer can be distinguished from healthy individuals with an Area Under the Curve of 0.96. Interrogating mutational processes in plasma may enable earlier cancer detection, and might enable the assessment of cancer risk and etiology.

Detection of mutational signatures in cell-free DNA (cfDNA) is challenging due to low sequence coverage and low mutant allele fractions. Here, the authors identify mutational signatures in plasma whole genome sequencing of cancer patients and use machine learning to distinguish them from healthy individuals.

Mutational signatures and processes in hepatobiliary cancers

The evolutionary history of hepatobiliary cancers is embedded in their genomes. By analysing their catalogue of somatic mutations and the DNA sequence context in which they occur, it is possible to infer the mechanisms underpinning tumorigenesis. These mutational signatures reflect the exogenous and endogenous origins of genetic damage as well as the capacity of hepatobiliary cells to repair and replicate DNA. Genomic analysis of thousands of patients with hepatobiliary cancers has highlighted the diversity of mutagenic processes active in these malignancies, highlighting a prominent source of the inter-cancer-type, inter-patient, intertumour and intratumoural heterogeneity that is observed clinically. However, a substantial proportion of mutational signatures detected in hepatocellular carcinoma and biliary tract cancer remain of unknown cause, emphasizing the important contribution of processes yet to be identified. Exploiting mutational signatures to retrospectively understand hepatobiliary carcinogenesis could advance preventative management of these aggressive tumours as well as potentially predict treatment response and guide the development of therapies targeting tumour evolution.

Premetastatic shifts of endogenous and exogenous mutational processes support consolidative therapy in EGFR-driven lung adenocarcinoma

The progression of cancer is an evolutionary process that is challenging to assess between sampling timepoints. However, investigation of cancer evolution over specific time periods is crucial to the elucidation of key events such as the acquisition of therapeutic resistance and subsequent fatal metastatic spread of therapy-resistant cell populations. Here we apply mutational signature analyses within clinically annotated cancer chronograms to detect and describe the shifting mutational processes caused by both endogenous (e.g. mutator gene mutation) and exogenous (e.g. mutagenic therapeutics) factors between tumor sampling timepoints. In one patient, we find that cisplatin therapy can introduce mutations that confer genetic resistance to subsequent targeted therapy with Erlotinib. In another patient, we trace detection of defective mismatch-repair associated mutational signature SBS3 to the emergence of known driver mutation CTNNB1 S37C. In both of these patients, metastatic lineages emerged from a single ancestral lineage that arose during therapy-a finding that argues for the consideration of local consolidative therapy over other therapeutic approaches in EGFR-positive non-small cell lung cancer. Broadly, these results demonstrate the utility of phylogenetic analysis that incorporates clinical time course and mutational signature deconvolution to inform therapeutic decision making and retrospective assessment of disease etiology.

Bridging Tumorigenesis and Therapy Resistance With a Non-Darwinian and Non-Lamarckian Mechanism of Adaptive Evolution

Although neo-Darwinian (and less often Lamarckian) dynamics are regularly invoked to interpret cancer's multifarious molecular profiles, they shine little light on how tumorigenesis unfolds and often fail to fully capture the frequency and breadth of resistance mechanisms. This uncertainty frames one of the most problematic gaps between science and practice in modern times. Here, we offer a theory of adaptive cancer evolution, which builds on a molecular mechanism that lies outside neo-Darwinian and Lamarckian schemes. This mechanism coherently integrates non-genetic and genetic changes, ecological and evolutionary time scales, and shifts the spotlight away from positive selection towards purifying selection, genetic drift, and the creative-disruptive power of environmental change. The surprisingly simple use-it or lose-it rationale of the proposed theory can help predict molecular dynamics during tumorigenesis. It also provides simple rules of thumb that should help improve therapeutic approaches in cancer.