Differential analysis between somatic mutation and germline variation profiles reveals cancer-related genes

Analysis of somatic mutations in whole blood from 200,618 individuals identifies pervasive positive selection and novel drivers of clonal hematopoiesis

Human aging is marked by the emergence of a tapestry of clonal expansions in dividing tissues, particularly evident in blood as clonal hematopoiesis (CH). CH, linked to cancer risk and aging-related phenotypes, often stems from somatic mutations in a set of established genes. However, the majority of clones lack known drivers. Here we infer gene-level positive selection in whole blood exomes from 200,618 individuals in UK Biobank. We identify 17 additional genes, ZBTB33, ZNF318, ZNF234, SPRED2, SH2B3, SRCAP, SIK3, SRSF1, CHEK2, CCDC115, CCL22, BAX, YLPM1, MYD88, MTA2, MAGEC3 and IGLL5, under positive selection at a population level, and validate this selection pattern in 10,837 whole genomes from single-cell-derived hematopoietic colonies. Clones with mutations in these genes grow in frequency and size with age, comparable to classical CH drivers. They correlate with heightened risk of infection, death and hematological malignancy, highlighting the significance of these additional genes in the aging process.

Familial CCM Genes Might Not Be Main Drivers for Pathogenesis of Sporadic CCMs-Genetic Similarity between Cancers and Vascular Malformations

Cerebral cavernous malformations (CCMs) are abnormally dilated intracranial capillaries that form cerebrovascular lesions with a high risk of hemorrhagic stroke. Recently, several somatic "activating" gain-of-function (GOF) point mutations in PIK3CA (phosphatidylinositol-4, 5-bisphosphate 3-kinase catalytic subunit p110α) were discovered as a dominant mutation in the lesions of sporadic forms of cerebral cavernous malformation (sCCM), raising the possibility that CCMs, like other types of vascular malformations, fall in the PIK3CA-related overgrowth spectrum (PROS). However, this possibility has been challenged with different interpretations. In this review, we will continue our efforts to expound the phenomenon of the coexistence of gain-of-function (GOF) point mutations in the PIK3CA gene and loss-of-function (LOF) mutations in CCM genes in the CCM lesions of sCCM and try to delineate the relationship between mutagenic events with CCM lesions in a temporospatial manner. Since GOF PIK3CA point mutations have been well studied in reproductive cancers, especially breast cancer as a driver oncogene, we will perform a comparative meta-analysis for GOF PIK3CA point mutations in an attempt to demonstrate the genetic similarities shared by both cancers and vascular anomalies.

Increased expression of GPX4 promotes the tumorigenesis of thyroid cancer by inhibiting ferroptosis and predicts poor clinical outcomes

BACKGROUND

Ferroptosis plays a critical role in suppressing cancer progression, and its essential regulator is glutathione peroxidase 4 (GPX4). High GPX4 expression can inhibit accumulation of iron, thus suppressing ferroptosis. However, its function in thyroid cancer has not been fully illuminated. Here, we explore the effect of GPX4 on thyroid cancer tumorigenesis and prognosis.

METHODS

Based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, GPX4 expression was investigated in cancer tissues and adjacent tissues. We determined the biological functions of GPX4-associated differentially expressed genes (DEGs) by using the "clusterProfiler" R package. In addition, the predictive value of GPX4 in thyroid cancer was assessed by using Cox regression analysis and nomograms. Finally, we conducted several in vitro experiments to determine the influence of GPX4 expression on proliferation and ferroptosis in thyroid cancer cells.

RESULTS

GPX4 expression was obviously elevated in thyroid cancer tissues compared with normal tissues. Biological function analysis indicated enrichment in muscle contraction, contractile fiber, metal ion transmembrane transporter activity, and complement and coagulation cascades. GPX4 overexpression was associated with stage T3-T4 and pathologic stage III-IV in thyroid cancer patients. Cox regression analysis indicated that GPX4 may be a risk factor for the overall survival of thyroid cancer patients. In vitro research showed that knockdown of GPX4 suppressed proliferation and induced ferroptosis in thyroid cancer cells.

CONCLUSIONS

GPX4 overexpression in thyroid cancer might play an essential role in tumorigenesis and may have prognostic value for thyroid cancer patients.

Pharmacogenetics of Anticancer Drugs: Clinical Response and Toxicity

Cancer is the most challenging disease for medical professionals to treat. The factors underlying the complicated situation include anticancer drug-associated toxicity, non-specific response, low therapeutic window, variable treatment outcomes, development of drug resistance, treatment complications, and cancer recurrence. The remarkable advancement in biomedical sciences and genetics, over the past few decades, however, is changing the dire situation. The discovery of gene polymorphism, gene expression, biomarkers, particular molecular targets and pathways, and drug-metabolizing enzymes have paved the way for the development and provision of targeted and individualized anticancer treatment. Pharmacogenetics is the study of genetic factors having the potential to affect clinical responses and pharmacokinetic and pharmacodynamic behaviors of drugs. This chapter emphasizes pharmacogenetics of anticancer drugs and its applications in improving treatment outcomes, selectivity, toxicity of the drugs, and discovering and developing personalized anticancer drugs and genetic methods for prediction of drug response and toxicity.

Repulsion and attraction in searching: A hybrid algorithm based on gravitational kernel and vital few for cancer driver gene prediction

By taking a new perspective to combine a machine learning method with an evolutionary algorithm, a new hybrid algorithm is developed to predict cancer driver genes. Firstly, inspired by the search strategy with the capability of global search in evolutionary algorithms, a gravitational kernel is proposed to act on the full range of gene features. Constructed by fusing PPI and mutation features, the gravitational kernel is capable to produce repulsion effects. The candidate genes with greater mutation effects and PPI have higher similarity scores. According to repulsion, the similarity score of these promising genes is larger than ordinary genes, which is beneficial to search for these promising genes. Secondly, inspired by the idea of elite populations related to evolutionary algorithms, the concept of vital few is proposed. Targeted at a local scale, it acts on the candidate genes associated with vital few genes. Under attraction effect, these vital few driver genes attract those with similar mutational effects to them, which leads to greater similarity scores. Lastly, the model and parameters are optimized by using an evolutionary algorithm, so as to obtain the optimal model and parameters for cancer driver gene prediction. Herein, a comparison is performed with six other advanced methods of cancer driver gene prediction. According to the experimental results, the method proposed in this study outperforms these six state-of-the-art algorithms on the pan-oncogene dataset.

KIAA1199 Correlates With Tumor Microenvironment and Immune Infiltration in Lung Adenocarcinoma as a Potential Prognostic Biomarker

Background: KIAA1199 has been considered a key regulator of carcinogenesis. However, the relationship between KIAA1199 and immune infiltrates, as well as its prognostic value in lung adenocarcinoma (LUAD) remains unclear. Methods: The expression of KIAA1199 and its influence on tumor prognosis were analyzed using a series of databases, comprising TIMER, GEPIA, UALCAN, LCE, Prognoscan and Kaplan-Meier Plotter. Further, immunohistochemistry (IHC), western blot (WB) and receiver operating characteristic (ROC) curve analyses were performed to verify our findings. The cBioPortal was used to investigate the genomic alterations of KIAA1199. Prediction of candidate microRNA (miRNAs) and transcription factor (TF) targeting KIAA1199, as well as GO and KEGG analyses, were performed based on LinkedOmics. TIMER and TISIDB databases were used to explore the relationship between KIAA1199 and tumor immune infiltration. Results: High expression of KIAA1199 was identified in LUAD and Lung squamous cell carcinoma (LUSC) patients. High expression of KIAA1199 indicated a worse prognosis in LUAD patients. The results of IHC and WB analyses showed that the expression level of KIAA1199 in tumor tissues was higher than that in adjacent tissues. GO and KEGG analyses indicated KIAA1199 was mainly involved in extracellular matrix (ECM)-receptor interaction and extracellular matrix structure constituent. KIAA1199 was positively correlated with infiltrating levels of CD4+ T cells, macrophages, neutrophil cells, dendritic cells, and showed positive relationship with immune marker subsets expression of a variety of immunosuppressive cells. Conclusion: High expression of KIAA1199 predicts a poor prognosis of LUAD patients. KIAA1199 might exert its carcinogenic role in the tumor microenvironment via participating in the extracellular matrix formation and regulating the infiltration of immune cells in LUAD. The results indicate that KIAA1199 might be a novel biomarker for evaluating prognosis and immune cell infiltration in LUAD.

Genomic data integration and user-defined sample-set extraction for population variant analysis

BACKGROUND

Population variant analysis is of great importance for gathering insights into the links between human genotype and phenotype. The 1000 Genomes Project established a valuable reference for human genetic variation; however, the integrative use of the corresponding data with other datasets within existing repositories and pipelines is not fully supported. Particularly, there is a pressing need for flexible and fast selection of population partitions based on their variant and metadata-related characteristics.

RESULTS

Here, we target general germline or somatic mutation data sources for their seamless inclusion within an interoperable-format repository, supporting integration among them and with other genomic data, as well as their integrated use within bioinformatic workflows. In addition, we provide VarSum, a data summarization service working on sub-populations of interest selected using filters on population metadata and/or variant characteristics. The service is developed as an optimized computational framework with an Application Programming Interface (API) that can be called from within any existing computing pipeline or programming script. Provided example use cases of biological interest show the relevance, power and ease of use of the API functionalities.

CONCLUSIONS

The proposed data integration pipeline and data set extraction and summarization API pave the way for solid computational infrastructures that quickly process cumbersome variation data, and allow biologists and bioinformaticians to easily perform scalable analysis on user-defined partitions of large cohorts from increasingly available genetic variation studies. With the current tendency to large (cross)nation-wide sequencing and variation initiatives, we expect an ever growing need for the kind of computational support hereby proposed.

A Germline Mutation in ATR Is Associated With Lung Adenocarcinoma in Asian Patients

BACKGROUND

Familial lung cancer (FLC) accounts for 8% of lung adenocarcinoma. It is known that a few germline mutations are associated with risk increasing and may provide new screening and treatment option. The goal of this study is to identify an FLC gene among three members of an FLC family.

METHODS

To uncover somatic and embryonic mutations linked with familial lung cancer, whole exome sequencing was done on surgical tissues and peripheral blood from three sisters in a family diagnosed with pulmonary lung adenocarcinoma (LUAD). At the same time, single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing data in public databases were enrolled to identify specific gene expression level.

RESULTS

Ataxia Telangiectasia and Rad3-Related Protein (ATR) gene C.7667C >G (p.T2556S) mutation were found in 3 patients with familial lung cancer. Whole-genome sequencing revealed that the three sisters exhibited similar somatic mutation patterns. Besides ATR mutations, common mutated genes (BRCA1, EGFR, and ROS1) that characterize LUAD were also found in 5 tumor samples. Analysis for the ATR expression in LUAD patients by single-cell sequencing data, we found ATR expression of tumor patients at high level in immune cells when compared with normal patients, but the expression of ATR in stromal cells has the opposite result.

CONCLUSION

We found a germline mutation in the ATR gene in three sisters of a Chinese family affected by familial lung cancer, which may be a genetic factor for lung cancer susceptibility.

Driver gene detection through Bayesian network integration of mutation and expression profiles

MOTIVATION

The identification of mutated driver genes and the corresponding pathways is one of the primary goals in understanding tumorigenesis at the patient level. Integration of multi-dimensional genomic data from existing repositories, e.g., The Cancer Genome Atlas (TCGA), offers an effective way to tackle this issue. In this study, we aimed to leverage the complementary genomic information of individuals and create an integrative framework to identify cancer-related driver genes. Specifically, based on pinpointed differentially expressed genes, variants in somatic mutations and a gene interaction network, we proposed an unsupervised Bayesian network integration (BNI) method to detect driver genes and estimate the disease propagation at the patient and/or cohort levels. This new method first captures inherent structural information to construct a functional gene mutation network and then extracts the driver genes and their controlled downstream modules using the minimum cover subset method.

RESULTS

Using other credible sources (e.g. Cancer Gene Census and Network of Cancer Genes), we validated the driver genes predicted by the BNI method in three TCGA pan-cancer cohorts. The proposed method provides an effective approach to address tumor heterogeneity faced by personalized medicine. The pinpointed drivers warrant further wet laboratory validation.

AVAILABILITY AND IMPLEMENTATION

The supplementary tables and source code can be obtained from https://xavieruniversityoflouisiana.sharefile.com/d-se6df2c8d0ebe4800a3030311efddafe5.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computationally scalable regression modeling for ultrahigh-dimensional omics data with ParProx

Statistical analysis of ultrahigh-dimensional omics scale data has long depended on univariate hypothesis testing. With growing data features and samples, the obvious next step is to establish multivariable association analysis as a routine method to describe genotype-phenotype association. Here we present ParProx, a state-of-the-art implementation to optimize overlapping and non-overlapping group lasso regression models for time-to-event and classification analysis, with selection of variables grouped by biological priors. ParProx enables multivariable model fitting for ultrahigh-dimensional data within an architecture for parallel or distributed computing via latent variable group representation. It thereby aims to produce interpretable regression models consistent with known biological relationships among independent variables, a property often explored post hoc, not during model estimation. Simulation studies clearly demonstrate the scalability of ParProx with graphics processing units in comparison to existing implementations. We illustrate the tool using three different omics data sets featuring moderate to large numbers of variables, where we use genomic regions and biological pathways as variable groups, rendering the selected independent variables directly interpretable with respect to those groups. ParProx is applicable to a wide range of studies using ultrahigh-dimensional omics data, from genome-wide association analysis to multi-omics studies where model estimation is computationally intractable with existing implementation.

A computational and structural analysis of germline and somatic variants affecting the DDR mechanism, and their impact on human diseases

DNA-Damage Response (DDR) proteins are crucial for maintaining the integrity of the genome by identifying and repairing errors in DNA. Variants affecting their function can have severe consequences since failure to repair damaged DNA can result in cells turning cancerous. Here, we compare germline and somatic variants in DDR genes, specifically looking at their locations in the corresponding three-dimensional (3D) structures, Pfam domains, and protein–protein interaction interfaces. We show that somatic variants in metastatic cases are more likely to be found in Pfam domains and protein interaction interfaces than are pathogenic germline variants or variants of unknown significance (VUS). We also show that there are hotspots in the structures of ATM and BRCA2 proteins where pathogenic germline, and recurrent somatic variants from primary and metastatic tumours, cluster together in 3D. Moreover, in the ATM, BRCA1 and BRCA2 genes from prostate cancer patients, the distributions of germline benign, pathogenic, VUS, and recurrent somatic variants differ across Pfam domains. Together, these results provide a better characterisation of the most recurrent affected regions in DDRs and could help in the understanding of individual susceptibility to tumour development.

PertInInt: An Integrative, Analytical Approach to Rapidly Uncover Cancer Driver Genes with Perturbed Interactions and Functionalities

A major challenge in cancer genomics is to identify genes with functional roles in cancer and uncover their mechanisms of action. We introduce an integrative framework that identifies cancer-relevant genes by pinpointing those whose interaction or other functional sites are enriched in somatic mutations across tumors. We derive analytical calculations that enable us to avoid time-prohibitive permutation-based significance tests, making it computationally feasible to simultaneously consider multiple measures of protein site functionality. Our accompanying software, PertInInt, combines knowledge about sites participating in interactions with DNA, RNA, peptides, ions, or small molecules with domain, evolutionary conservation, and gene-level mutation data. When applied to 10,037 tumor samples, PertInInt uncovers both known and newly predicted cancer genes, while additionally revealing what types of interactions or other functionalities are disrupted. PertInInt’s analysis demonstrates that somatic mutations are frequently enriched in interaction sites and domains and implicates interaction perturbation as a pervasive cancer-driving event.

A fast, analytical framework called PertInInt enables efficient integration of multiple measures of protein site functionality—including interaction, domain, and evolutionary conservation—with gene-level mutation data in order to rapidly detect cancer driver genes along with their disrupted functionalities.

uKIN Combines New and Prior Information with Guided Network Propagation to Accurately Identify Disease Genes

Protein interaction networks provide a powerful framework for identifying genes causal for complex genetic diseases. Here, we introduce a general framework, uKIN, that uses prior knowledge of disease-associated genes to guide, within known protein-protein interaction networks, random walks that are initiated from newly identified candidate genes. In large-scale testing across 24 cancer types, we demonstrate that our network propagation approach for integrating both prior and new information not only better identifies cancer driver genes than using either source of information alone but also readily outperforms other state-of-the-art network-based approaches. We also apply our approach to genome-wide association data to identify genes functionally relevant for several complex diseases. Overall, our work suggests that guided network propagation approaches that utilize both prior and new data are a powerful means to identify disease genes. uKIN is freely available for download at: https://github.com/Singh-Lab/uKIN.

Differential Allele-Specific Expression Uncovers Breast Cancer Genes Dysregulated by Cis Noncoding Mutations

Identifying cancer-relevant mutations in noncoding regions is challenging due to the large numbers of such mutations, their low levels of recurrence, and difficulties in interpreting their functional impact. To uncover genes that are dysregulated due to somatic mutations in cis, we build upon the concept of differential allele-specific expression (ASE) and introduce methods to identify genes within an individual's cancer whose ASE differs from what is found in matched normal tissue. When applied to breast cancer tumor samples, our methods detect the known allele-specific effects of copy number variation and nonsense-mediated decay. Further, genes that are found to recurrently exhibit differential ASE across samples are cancer relevant. Genes with cis mutations are enriched for differential ASE, and we find 147 potentially functional noncoding mutations cis to genes that exhibit significant differential ASE. We conclude that differential ASE is a promising means for discovering gene dysregulation due to cis noncoding mutations.

Quantifying gene selection in cancer through protein functional alteration bias

Compiling the catalogue of genes actively involved in cancer is an ongoing endeavor, with profound implications to the understanding and treatment of the disease. An abundance of computational methods have been developed to screening the genome for candidate driver genes based on genomic data of somatic mutations in tumors. Existing methods make many implicit and explicit assumptions about the distribution of random mutations. We present FABRIC, a new framework for quantifying the selection of genes in cancer by assessing the effects of de-novo somatic mutations on protein-coding genes. Using a machine-learning model, we quantified the functional effects of ∼3M somatic mutations extracted from over 10 000 human cancerous samples, and compared them against the effects of all possible single-nucleotide mutations in the coding human genome. We detected 593 protein-coding genes showing statistically significant bias towards harmful mutations. These genes, discovered without any prior knowledge, show an overwhelming overlap with known cancer genes, but also include many overlooked genes. FABRIC is designed to avoid false discoveries by comparing each gene to its own background model using rigorous statistics, making minimal assumptions about the distribution of random somatic mutations. The framework is an open-source project with a simple command-line interface.

LOTUS: A single- and multitask machine learning algorithm for the prediction of cancer driver genes

Cancer driver genes, i.e., oncogenes and tumor suppressor genes, are involved in the acquisition of important functions in tumors, providing a selective growth advantage, allowing uncontrolled proliferation and avoiding apoptosis. It is therefore important to identify these driver genes, both for the fundamental understanding of cancer and to help finding new therapeutic targets or biomarkers. Although the most frequently mutated driver genes have been identified, it is believed that many more remain to be discovered, particularly for driver genes specific to some cancer types. In this paper, we propose a new computational method called LOTUS to predict new driver genes. LOTUS is a machine-learning based approach which allows to integrate various types of data in a versatile manner, including information about gene mutations and protein-protein interactions. In addition, LOTUS can predict cancer driver genes in a pan-cancer setting as well as for specific cancer types, using a multitask learning strategy to share information across cancer types. We empirically show that LOTUS outperforms five other state-of-the-art driver gene prediction methods, both in terms of intrinsic consistency and prediction accuracy, and provide predictions of new cancer genes across many cancer types.

Cancer development is driven by mutations and dysfunction of important, so-called cancer driver genes, that could be targeted by specific therapies. While a number of such cancer genes have already been identified, it is believed that many more remain to be discovered. To help prioritize experimental investigations of candidate genes, several computational methods have been proposed to rank promising candidates based on their mutations in large cohorts of cancer cases, or on their interactions with known driver genes in biological networks. We propose LOTUS, a new computational approach to identify genes with high oncogenic potential. LOTUS implements a machine learning approach to learn an oncogenic potential score from known driver genes, and brings two novelties compared to existing methods. First, it allows to easily combine heterogeneous sources of information into the scoring function, which we illustrate by learning a scoring function from both known mutations in large cancer cohorts and interactions in biological networks. Second, using a multitask learning strategy, it can predict different driver genes for different cancer types, while sharing information between them to improve the prediction for every type. We provide experimental results showing that LOTUS significantly outperforms several state-of-the-art cancer gene prediction software.

From Single Level Analysis to Multi-Omics Integrative Approaches: A Powerful Strategy towards the Precision Oncology

Integration of multi-omics data from different molecular levels with clinical data, as well as epidemiologic risk factors, represents an accurate and promising methodology to understand the complexity of biological systems of human diseases, including cancer. By the extensive use of novel technologic platforms, a large number of multidimensional data can be derived from analysis of health and disease systems. Comprehensive analysis of multi-omics data in an integrated framework, which includes cumulative effects in the context of biological pathways, is therefore eagerly awaited. This strategy could allow the identification of pathway-addiction of cancer cells that may be amenable to therapeutic intervention. However, translation into clinical settings requires an optimized integration of omics data with clinical vision to fully exploit precision cancer medicine. We will discuss the available technical approach and more recent developments in the specific field.

Autologous reference types can confound the detection of somatic mutation in solid cancers

Vast number of somatic mutations has been proved to be affected by the factors of sequencing methods, analysis pipelines and validation methods. We here showed the effect of autologous reference types on the detection of cancer-associated somatic mutations with the somatic single nucleotide variations (SNVs) and clinical data of solid tumors from the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC). The distribution of somatic SNVs was significantly different among groups of autologous references in 6 cancers detected by whole genome sequencing (WGS) and 5 cancers detected by the random sequencing of exonic regions selected from the genome (WXS), especially in protein coding region of 5 cancers with age, gender and TNM adjusted. In addition, only 60.24% (95% CI: 49.65%-70.83%) of the somatic SNVs called from normal blood by WXS were found in those called from normal solid tissue tested by WXS / WGS, while 31.78% (95%CI: 4.14%-59.42%) of the somatic SNVs called from normal tissue adjacent to primary by WXS were found in those from normal blood tested by WXS / WGS. These findings suggested that more representative types of normal tissues should be included in detection of cancer-associated somatic mutations.

Network-Based Coverage of Mutational Profiles Reveals Cancer Genes

A central goal in cancer genomics is to identify the somatic alterations that underpin tumor initiation and progression. While commonly mutated cancer genes are readily identifiable, those that are rarely mutated across samples are difficult to distinguish from the large numbers of other infrequently mutated genes. We introduce a method, nCOP, that considers per-individual mutational profiles within the context of protein-protein interaction networks in order to identify small connected subnetworks of genes that, while not individually frequently mutated, comprise pathways that are altered across (i.e., "cover") a large fraction of individuals. By analyzing 6,038 samples across 24 different cancer types, we demonstrate that nCOP is highly effective in identifying cancer genes, including those with low mutation frequencies. Overall, our work demonstrates that combining per-individual mutational information with interaction networks is a powerful approach for tackling the mutational heterogeneity observed across cancers.