iGC-an integrated analysis package of gene expression and copy number alteration

From driver genes to gene families: A computational analysis of oncogenic mutations and ubiquitination anomalies in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a widespread primary liver cancer with a high fatality rate. Despite several genes with oncogenic effects in HCC have been identified, many remain undiscovered. In this study, we conducted a comprehensive computational analysis to explore the involvement of genes within the same families as known driver genes in HCC. Specifically, we expanded the concept beyond single-gene mutations to encompass gene families sharing homologous structures, integrating various omics data to comprehensively understand gene abnormalities in cancer. Our analysis identified 74 domains with an enriched mutation burden, 404 domain mutation hotspots, and 233 dysregulated driver genes. We observed that specific low-frequency somatic mutations may contribute to HCC occurrence, potentially overlooked by single-gene algorithms. Furthermore, we systematically analyzed how abnormalities in the ubiquitinated proteasome system (UPS) impact HCC, finding that abnormal genes in E3, E2, DUB families, and Degron genes often result in HCC by affecting the stability of oncogenic or tumor suppressor proteins. In conclusion, expanding the exploration of driver genes to include gene families with homologous structures emerges as a promising strategy for uncovering additional oncogenic alterations in HCC.

Functional Copy-Number Alterations as Diagnostic and Prognostic Biomarkers in Neuroendocrine Tumors

Functional copy-number alterations (fCNAs) are DNA copy-number changes with concordant differential gene expression. These are less likely to be bystander genetic lesions and could serve as robust and reproducible tumor biomarkers. To identify candidate fCNAs in neuroendocrine tumors (NETs), we integrated chromosomal microarray (CMA) and RNA-seq differential gene-expression data from 31 pancreatic (pNETs) and 33 small-bowel neuroendocrine tumors (sbNETs). Tumors were resected from 47 early-disease-progression (<24 months) and 17 late-disease-progression (>24 months) patients. Candidate fCNAs that accurately differentiated these groups in this discovery cohort were then replicated using fluorescence in situ hybridization (FISH) on formalin-fixed, paraffin-embedded (FFPE) tissues in a larger validation cohort of 60 pNETs and 82 sbNETs (52 early- and 65 late-disease-progression samples). Logistic regression analysis revealed the predictive ability of these biomarkers, as well as the assay-performance metrics of sensitivity, specificity, and area under the curve. Our results indicate that copy-number changes at chromosomal loci 4p16.3, 7q31.2, 9p21.3, 17q12, 18q21.2, and 19q12 may be used as diagnostic and prognostic NET biomarkers. This involves a rapid, cost-effective approach to determine the primary tumor site for patients with metastatic liver NETs and to guide risk-stratified therapeutic decisions.

Mutational landscape of HSP family on human breast cancer

Breast cancer (BRCA) is a prevalent malignancy with the highest incidence among females. BRCA can be categorized into five intrinsic molecular subtypes (LumA, LumB, HER2, Basal, and Normal), each characterized by varying molecular and clinical features determined by the expression of intrinsic genes (PAM50). The Heat Shock Protein (HSP) family is composed of 95 genes evolutionary conservated, they have critical roles in proteostasis in both normal and cancerous processes. Many studies have linked HSP to the development and spread of cancer. They modulate the activity of multiple proteins expressed by oncogenes and anti-oncogenes through a range of interactions. In this study, we evaluate the mutational changes that HSP undergoes in BRCA mainly from the TCGA database. We observe that Copy Number Variations (CNV) are the more frequent events analyzed surpassing the occurrence of point mutations, indels, and translation start site mutations. The Basal subtype showcased the highest count of amplified CNV, including subtype-specific changes, whereas the Luminals tumors accumulated the greatest number of deletion CNV. Meanwhile, the HER2 subtype exhibited a comparatively lower frequency of CNV alterations when compared to the other subtypes. This study integrates CNV and expression data, finding associations between these two variables and the influence of CNV on the deregulation of HSP expression. To enhance the role of HSP as a risk predictor in BRCA, we succeeded in identifying CNV profiles as a prognostic marker. We included Artificial Intelligence to improve the clustering of patients, and we achieved a molecular CNV signature as a significant risk factor independent of known classic markers, including molecular subtypes PAM50. This research enhances the comprehension of HSP DNA alterations in BRCA and its relation with predicting the risk of affected individuals providing insights to develop guide personalized treatment strategies.

Machine Learning Methods for Gene Selection in Uveal Melanoma

Uveal melanoma (UM) is the most common primary intraocular malignancy with a limited five-year survival for metastatic patients. Limited therapeutic treatments are currently available for metastatic disease, even if the genomics of this tumor has been deeply studied using next-generation sequencing (NGS) and functional experiments. The profound knowledge of the molecular features that characterize this tumor has not led to the development of efficacious therapies, and the survival of metastatic patients has not changed for decades. Several bioinformatics methods have been applied to mine NGS tumor data in order to unveil tumor biology and detect possible molecular targets for new therapies. Each application can be single domain based while others are more focused on data integration from multiple genomics domains (as gene expression and methylation data). Examples of single domain approaches include differentially expressed gene (DEG) analysis on gene expression data with statistical methods such as SAM (significance analysis of microarray) or gene prioritization with complex algorithms such as deep learning. Data fusion or integration methods merge multiple domains of information to define new clusters of patients or to detect relevant genes, according to multiple NGS data. In this work, we compare different strategies to detect relevant genes for metastatic disease prediction in the TCGA uveal melanoma (UVM) dataset. Detected targets are validated with multi-gene score analysis on a larger UM microarray dataset.

Empirical comparison and analysis of machine learning-based approaches for druggable protein identification

Efficiently and precisely identifying drug targets is crucial for developing and discovering potential medications. While conventional experimental approaches can accurately pinpoint these targets, they suffer from time constraints and are not easily adaptable to high-throughput processes. On the other hand, computational approaches, particularly those utilizing machine learning (ML), offer an efficient means to accelerate the prediction of druggable proteins based solely on their primary sequences. Recently, several state-of-the-art computational methods have been developed for predicting and analyzing druggable proteins. These computational methods showed high diversity in terms of benchmark datasets, feature extraction schemes, ML algorithms, evaluation strategies and webserver/software usability. Thus, our objective is to reexamine these computational approaches and conduct a comprehensive assessment of their strengths and weaknesses across multiple aspects. In this study, we deliver the first comprehensive survey regarding the state-of-the-art computational approaches for in silico prediction of druggable proteins. First, we provided information regarding the existing benchmark datasets and the types of ML methods employed. Second, we investigated the effectiveness of these computational methods in druggable protein identification for each benchmark dataset. Third, we summarized the important features used in this field and the existing webserver/software. Finally, we addressed the present constraints of the existing methods and offer valuable guidance to the scientific community in designing and developing novel prediction models. We anticipate that this comprehensive review will provide crucial information for the development of more accurate and efficient druggable protein predictors.

Multi-Omics Analysis Reveals Clinical Value and Possible Mechanisms of ATAD1 Down-Regulation in Human Prostate Adenocarcinoma

Prostate adenocarcinoma (PRAD) is the most common histological subtype of prostate cancer. Post-treatment biochemical recurrence is a challenging issue. ATAD1 (ATPase Family AAA Domain Containing 1) plays a vital role in mitochondrial proteostasis and apoptosis activity, while its clinical value in PRAD and its impact on the tumor microenvironment (TME) remain unanswered. In this study, we aimed to investigate the clinical value and possible mechanisms of ATAD1 in PRAD via multi-omics analysis. Using cBioPortal, we confirmed that ATAD1 alteration was associated with gene expression and unfavorable DFS. Deep deletion predominantly occurred in PRAD. By integrating DriverDBv3 and GEPIA2, we noted ATAD1 downregulation in PRAD tissues compared to normal tissues, associated with unfavorable DFS in PRAD patients. DNA repair genes ATM, PARP1and BRCA2 had positive associations with ATAD1 expression. We found that the generalization value of ATAD1 could be applied to other cancers such as KIRC and UCEC. In addition, LinkedOmics identified that the functional involvement of ATAD1 participates in mitochondrial structure and cell cycle progression. Using TIMER analysis, we demonstrated that ATAD1 downregulation correlated with an immunosuppressive TME. Furthermore, we accessed a GSE55062 dataset on UALCAN and discovered the involvement of ERG-mediated transcriptional repression on ATAD1 downregulation. Cross-association screening of shATAD1 efficacy vs. altered mRNAs identified 190 perturbed mRNAs. Then, functional enrichment analysis using the Metascape omics tool recognized that shATAD1-perturbed mRNAs are primarily in charge of the activation of Wnt/β-catenin pathway and lipid metabolic processes. In conclusion, multi-omics results reveal that ATAD1 downregulation is a clinical biomarker for pathological diagnosis and prognosis for patients with PRAD. Reduced ATAD1 may be involved in the enhanced activity of mitochondria and cell cycle, as well as possibly shaping an immunosuppressive TME. ERG serves as an upstream transcriptional repressor of ATAD1. Downstream mechanisms of ATAD1 are involved in Wnt/β-catenin pathway and lipid metabolic processes.

Multi-omics data integration approaches for precision oncology

Next-generation sequencing (NGS) has been pivotal to enhance the molecular characterization of human malignancies, allowing multiple omics data types to be available for cancer researchers and practitioners. In this context,...

High prevalence of APOA1/C3/A4/A5 alterations in luminal breast cancers among young women in East Asia

In East Asia, the breast cancer incidence rate among women aged <50 years has rapidly increased. Emerging tumors are distinctly characterized by a high prevalence of estrogen receptor (ER)-positive/human epidermal growth factor receptor (HER2)-negative cancer. In the present study, we identified unique genetic alterations in these emerging tumors. We analyzed gene copy number variations (CNVs) in breast tumors from 120 Taiwanese patients, and obtained public datasets of CNV and gene expression (GE). The data regarding CNV and GE were separately compared between East Asian and Western patients, and the overlapping genes identified in the comparisons were explored to identify the gene-gene interaction networks. In the age <50 years/ER + /HER2- subgroup, tumors of East Asian patients exhibited a higher frequency of copy number loss in APOA1/C3/A4/A5, a lipid-metabolizing gene cluster (33 vs. 10%, P < .001) and lower APOA1/C3/A4/A5 expressions than tumors of Western patients. These copy number loss related- and GE-related results were validated in another Taiwanese cohort and in two GE datasets, respectively. The copy number loss was significantly associated with poor survival among Western patients, but not among East Asian patients. Lower APOA1, APOC3, and APOA5 expressions were associated with higher ESTIMATE immune scores, indicating an abundance of tumor-infiltrating immune cells. In conclusion, APOA1/C3/A4/A5 copy number loss was more prevalent in luminal breast tumors among East Asian women aged <50 years, and its immunomodulatory effect on the tumor microenvironment possibly plays various roles in the tumor biology of East Asian patients.

TPQCI: A topology potential-based method to quantify functional influence of copy number variations

Copy number variation (CNV) is a major type of chromosomal structural variation that play important roles in many diseases including cancers. Due to genome instability, a large number of CNV events can be detected in diseases such as cancer. Therefore, it is important to identify the functionally important CNVs in diseases, which currently still poses a challenge in genomics. One of the critical steps to solve the problem is to define the influence of CNV. In this paper, we provide a topology potential based method, TPQCI, to quantify this kind of influence by integrating statistics, gene regulatory associations, and biological function information. We used this metric to detect functionally enriched genes on genomic segments with CNV in breast cancer and multiple myeloma and discovered biological functions influenced by CNV. Our results demonstrate that, by using our proposed TPQCI metric, we can detect disease-specific genes that are influenced by CNVs. Source codes of TPQCI are provided in Github (https://github.com/usos/TPQCI).

CNA2Subpathway: identification of dysregulated subpathway driven by copy number alterations in cancer

Biological pathways reflect the key cellular mechanisms that dictate disease states, drug response and altered cellular function. The local areas of pathways are defined as subpathways (SPs), whose dysfunction has been reported to be associated with the occurrence and development of cancer. With the development of high-throughput sequencing technology, identifying dysfunctional SPs by using multi-omics data has become possible. Moreover, the SPs are not isolated in the biological system but interact with each other. Here, we propose a network-based calculated method, CNA2Subpathway, to identify dysfunctional SPs is driven by somatic copy number alterations (CNAs) in cancer through integrating pathway topology information, multi-omics data and SP crosstalk. This provides a novel way of SP analysis by using the SP interactions in the system biological level. Using data sets from breast cancer and head and neck cancer, we validate the effectiveness of CNA2Subpathway in identifying cancer-relevant SPs driven by the somatic CNAs, which are also shown to be associated with cancer immune and prognosis of patients. We further compare our results with five pathway or SP analysis methods based on CNA and gene expression data without considering SP crosstalk. With these analyses, we show that CNA2Subpathway could help to uncover dysfunctional SPs underlying cancer via the use of SP crosstalk. CNA2Subpathway is developed as an R-based tool, which is freely available on GitHub (https://github.com/hanjunwei-lab/CNA2Subpathway).

Multi-Omics Approaches to Improve Mitochondrial Disease Diagnosis: Challenges, Advances, and Perspectives

Mitochondrial diseases (MD) are rare disorders caused by deficiency of the mitochondrial respiratory chain, which provides energy in each cell. They are characterized by a high clinical and genetic heterogeneity and in most patients, the responsible gene is unknown. Diagnosis is based on the identification of the causative gene that allows genetic counseling, prenatal diagnosis, understanding of pathological mechanisms, and personalized therapeutic approaches. Despite the emergence of Next Generation Sequencing (NGS), to date, more than one out of two patients has no diagnosis in the absence of identification of the responsible gene. Technologies currently used for detecting causal variants (genetic alterations) is far from complete, leading many variants of unknown significance (VUS) and mainly based on the use of whole exome sequencing thus neglecting the identification of non-coding variants. The complexity of human genome and its regulation at multiple levels has led biologists to develop several assays to interrogate the different aspects of biological processes. While one-dimension single omics investigation offers a peek of this complex system, the combination of different omics data allows the discovery of coherent signatures. The community of computational biologists and bioinformaticians, in order to integrate data from different omics, has developed several approaches and tools. However, it is difficult to understand which suits the best to predict diverse phenotypic outcome. First attempts to use multi-omics approaches showed an improvement of the diagnostic power. However, we are far from a complete understanding of MD and their diagnosis. After reviewing multi-omics algorithms developed in the latest years, we are proposing here a novel data-driven classification and we will discuss how multi-omics will change and improve the diagnosis of MD. Due to the growing use of multi-omics approaches in MD, we foresee that this work will contribute to set up good practices to perform multi-omics data integration to improve the prediction of phenotypic outcomes and the diagnostic power of MD.

DriverDBv3: a multi-omics database for cancer driver gene research

An integrative multi-omics database is needed urgently, because focusing only on analysis of one-dimensional data falls far short of providing an understanding of cancer. Previously, we presented DriverDB, a cancer driver gene database that applies published bioinformatics algorithms to identify driver genes/mutations. The updated DriverDBv3 database (http://ngs.ym.edu.tw/driverdb) is designed to interpret cancer omics’ sophisticated information with concise data visualization. To offer diverse insights into molecular dysregulation/dysfunction events, we incorporated computational tools to define CNV and methylation drivers. Further, four new features, CNV, Methylation, Survival, and miRNA, allow users to explore the relations from two perspectives in the ‘Cancer’ and ‘Gene’ sections. The ‘Survival’ panel offers not only significant survival genes, but gene pairs synergistic effects determine. A fresh function, ‘Survival Analysis’ in ‘Customized-analysis,’ allows users to investigate the co-occurring events in user-defined gene(s) by mutation status or by expression in a specific patient group. Moreover, we redesigned the web interface and provided interactive figures to interpret cancer omics’ sophisticated information, and also constructed a Summary panel in the ‘Cancer’ and ‘Gene’ sections to visualize the features on multi-omics levels concisely. DriverDBv3 seeks to improve the study of integrative cancer omics data by identifying driver genes and contributes to cancer biology.

A comparative study of multi-omics integration tools for cancer driver gene identification and tumour subtyping

Oncogenesis and cancer can arise as a consequence of a wide range of genomic aberrations including mutations, copy number alterations, expression changes and epigenetic modifications encompassing multiple omics layers. Integrating genomic, transcriptomic, proteomic and epigenomic datasets via multi-omics analysis provides the opportunity to derive a deeper and holistic understanding of the development and progression of cancer. There are two primary approaches to integrating multi-omics data: multi-staged (focused on identifying genes driving cancer) and meta-dimensional (focused on establishing clinically relevant tumour or sample classifications). A number of ready-to-use bioinformatics tools are available to perform both multi-staged and meta-dimensional integration of multi-omics data. In this study, we compared nine different integration tools using real and simulated cancer datasets. The performance of the multi-staged integration tools were assessed at the gene, function and pathway levels, while meta-dimensional integration tools were assessed based on the sample classification performance. Additionally, we discuss the influence of factors such as data representation, sample size, signal and noise on multi-omics data integration. Our results provide current and much needed guidance regarding selection and use of the most appropriate and best performing multi-omics integration tools.

Identification of candidate cancer drivers by integrative Epi-DNA and Gene Expression (iEDGE) data analysis

The emergence of large-scale multi-omics data warrants method development for data integration. Genomic studies from cancer patients have identified epigenetic and genetic regulators - such as methylation marks, somatic mutations, and somatic copy number alterations (SCNAs), among others - as predictive features of cancer outcome. However, identification of "driver genes" associated with a given alteration remains a challenge. To this end, we developed a computational tool, iEDGE, to model cis and trans effects of (epi-)DNA alterations and identify potential cis driver genes, where cis and trans genes denote those genes falling within and outside the genomic boundaries of a given (epi-)genetic alteration, respectively. iEDGE first identifies the cis and trans gene expression signatures associated with the presence/absence of a particular epi-DNA alteration across samples. It then applies tests of statistical mediation to determine the cis genes predictive of the trans gene expression. Finally, cis and trans effects are annotated by pathway enrichment analysis to gain insights into the underlying regulatory networks. We used iEDGE to perform integrative analysis of SCNAs and gene expression data from breast cancer and 18 additional cancer types included in The Cancer Genome Atlas (TCGA). Notably, cis gene drivers identified by iEDGE were found to be significantly enriched for known driver genes from multiple compendia of validated oncogenes and tumor suppressors, suggesting that the remainder are of equal importance. Furthermore, predicted drivers were enriched for functionally relevant cancer genes with amplification-driven dependencies, which are of potential prognostic and therapeutic value. All the analyses results are accessible at https://montilab.bu.edu/iEDGE. In summary, integrative analysis of SCNAs and gene expression using iEDGE successfully identified known cancer driver genes and putative cancer therapeutic targets across 19 cancer types in the TCGA. The proposed method can easily be applied to the integration of gene expression profiles with other epi-DNA assays in a variety of disease contexts.

Using R and Bioconductor in Clinical Genomics and Transcriptomics

Bioinformatics pipelines are essential in the analysis of genomic and transcriptomic data generated by next-generation sequencing (NGS). Recent guidelines emphasize the need for rigorous validation and assessment of robustness, reproducibility, and quality of NGS analytic pipelines intended for clinical use. Software tools written in the R statistical language and, in particular, the set of tools available in the Bioconductor repository are widely used in research bioinformatics; and these frameworks offer several advantages for use in clinical bioinformatics, including the breath of available tools, modular nature of software packages, ease of installation, enforcement of interoperability, version control, and short learning curve. This review provides an introduction to R and Bioconductor software, its advantages and limitations for clinical bioinformatics, and illustrative examples of tools that can be used in various steps of NGS analysis.

Integrated landscape of copy number variation and RNA expression associated with nodal metastasis in invasive ductal breast carcinoma

Background

Lymph node metastasis (NM) in breast cancer is a clinical predictor of patient outcomes, but how its genetic underpinnings contribute to aggressive phenotypes is unclear. Our objective was to create the first landscape analysis of CNV-associated NM in ductal breast cancer. To assess the role of copy number variations (CNVs) in NM, we compared CNVs and/or associated mRNA expression in primary tumors of patients with NM to those without metastasis.

Results

We found CNV loss in chromosomes 1, 3, 9, 18, and 19 and gains in chromosomes 5, 8, 12, 14, 16-17, and 20 that were associated with NM and replicated in both databases. In primary tumors, per-gene CNVs associated with NM were ten times more frequent than mRNA expression; however, there were few CNV-driven changes in mRNA expression that differed by nodal status. Overlapping regions of CNV changes and mRNA expression were evident for the CTAGE5 gene. In 8q12, 11q13-14, 20q1, and 17q14-24 regions, there were gene-specific gains in CNV-driven mRNA expression associated with NM.

Methods

Data on CNV and mRNA expression from the TCGA and the METABRIC consortium of breast ductal carcinoma were utilized to identify CNV-based features associated with NM. Within each dataset, associations were compared across omic platforms to identify CNV-driven variations in gene expression. Only replications across both datasets were considered as determinants of NM.

Conclusions

Gains in CTAGE5, NDUFC2, EIF4EBP1, and PSCA genes and their expression may aid in early diagnosis of metastatic breast carcinoma and have potential as therapeutic targets.