Integrative Analysis Of Multi-omics Data Research Articles

Proteogenomic insights into the biology and treatment of pan-melanoma

Melanoma is one of the most prevalent skin cancers, with high metastatic rates and poor prognosis. Understanding its molecular pathogenesis is crucial for improving its diagnosis and treatment. Integrated analysis of multi-omics data from 207 treatment-naïve melanomas (primary-cutaneous-melanomas (CM, n = 28), primary-acral-melanomas (AM, n = 81), primary-mucosal-melanomas (MM, n = 28), metastatic-melanomas (n = 27), and nevi (n = 43)) provides insights into melanoma biology. Multivariate analysis reveals that PRKDC amplification is a prognostic molecule for melanomas. Further proteogenomic analysis combined with functional experiments reveals that the cis-effect of PRKDC amplification may lead to tumor proliferation through the activation of DNA repair and folate metabolism pathways. Proteome-based stratification of primary melanomas defines three prognosis-related subtypes, namely, the ECM subtype, angiogenesis subtype (with a high metastasis rate), and cell proliferation subtype, which provides an essential framework for the utilization of specific targeted therapies for particular melanoma subtypes. The immune classification identifies three immune subtypes. Further analysis combined with an independent anti-PD-1 treatment cohort reveals that upregulation of the MAPK7-NFKB signaling pathway may facilitate T-cell recruitment and increase the sensitivity of patients to immunotherapy. In contrast, PRKDC may reduce the sensitivity of melanoma patients to immunotherapy by promoting DNA repair in melanoma cells. These results emphasize the clinical value of multi-omics data and have the potential to improve the understanding of melanoma treatment.

Network-based integrative multi-omics approach reveals biosignatures specific to COVID-19 disease phases.

COVID-19 disease is characterized by a spectrum of disease phases (mild, moderate, and severe). Each disease phase is marked by changes in omics profiles with corresponding changes in the expression of features (biosignatures). However, integrative analysis of multiple omics data from different experiments across studies to investigate biosignatures at various disease phases is limited. Exploring an integrative multi-omics profile analysis through a network approach could be used to determine biosignatures associated with specific disease phases and enable the examination of the relationships between the biosignatures. To identify and characterize biosignatures underlying various COVID-19 disease phases in an integrative multi-omics data analysis. We leveraged a multi-omics network-based approach to integrate transcriptomics, metabolomics, proteomics, and lipidomics data. The World Health Organization Ordinal Scale WHO Ordinal Scale was used as a disease severity reference to harmonize COVID-19 patient metadata across two studies with independent data. A unified COVID-19 knowledge graph was constructed by assembling a disease-specific interactome from the literature and databases. Disease-state specific omics-graphs were constructed by integrating multi-omics data with the unified COVID-19 knowledge graph. We expanded on the network layers of multiXrank, a random walk with restart on multilayer network algorithm, to explore disease state omics-specific graphs and perform enrichment analysis. Network analysis revealed the biosignatures involved in inducing chemokines and inflammatory responses as hubs in the severe and moderate disease phases. We observed distinct biosignatures between severe and moderate disease phases as compared to mild-moderate and mild-severe disease phases. Mild COVID-19 cases were characterized by a unique biosignature comprising C-C Motif Chemokine Ligand 4 (CCL4), and Interferon Regulatory Factor 1 (IRF1). Hepatocyte Growth Factor (HGF), Matrix Metallopeptidase 12 (MMP12), Interleukin 10 (IL10), Nuclear Factor Kappa B Subunit 1 (NFKB1), and suberoylcarnitine form hubs in the omics network that characterizes the moderate disease state. The severe cases were marked by biosignatures such as Signal Transducer and Activator of Transcription 1 (STAT1), Superoxide Dismutase 2 (SOD2), HGF, taurine, lysophosphatidylcholine, diacylglycerol, triglycerides, and sphingomyelin that characterize the disease state. This study identified both biosignatures of different omics types enriched in disease-related pathways and their associated interactions (such as protein-protein, protein-transcript, protein-metabolite, transcript-metabolite, and lipid-lipid interactions) that are unique to mild, moderate, and severe COVID-19 disease states. These biosignatures include molecular features that underlie the observed clinical heterogeneity of COVID-19 and emphasize the need for disease-phase-specific treatment strategies. The approach implemented here can be used to find associations between transcripts, proteins, lipids, and metabolites in other diseases.

Predicting regulatory mutations and their target genes by new computational integrative analysis: A study of follicular lymphoma

Mutations in DNA regulatory regions are increasingly being recognized as important drivers of cancer and other complex diseases. These mutations can regulate gene expression by affecting DNA-protein binding and epigenetic profiles, such as DNA methylation in genome regulatory elements. However, identifying mutation hotspots associated with expression regulation and disease progression in non-coding DNA remains a challenge. Unlike most existing approaches that assign a mutation score to individual single nucleotide polymorphisms (SNP), a mutation block (MB)-based approach was introduced in this study to assess the collective impact of a cluster of SNPs on transcription factor-DNA binding affinity, differential gene expression (DEG), and nearby DNA methylation. Moreover, the long-distance target genes of functional MBs were identified using a new permutation-based algorithm that assessed the significance of correlations between DNA methylation at regulatory regions and target gene expression. Two new Python packages were developed. The Differential Methylation Region (DMR-analysis) analysis tool was used to detect DMR and map them to regulatory elements. The second tool, an integrated DMR, DEG, and SNP analysis tool (DDS-analysis), was used to combine the omics data to identify functional MBs and long-distance target genes. Both tools were validated in follicular lymphoma (FL) cohorts, where not only known functional MBs and their target genes (BCL2 and BCL6) were recovered, but also novel genes were found, including CDCA4 and JAG2, which may be associated with FL development. These genes are linked to target gene expression and are significantly correlated with the methylation of nearby DNA sequences in FL. The proposed computational integrative analysis of multiomics data holds promise for identifying regulatory mutations in cancer and other complex diseases.

Mechanistic safety assessment via multi-omic characterisation of systemic pathway perturbations following in vivo MAT2A inhibition.

The tumour suppressor p16/CDKN2A and the metabolic gene, methyl-thio-adenosine phosphorylase (MTAP), are frequently co-deleted in some of the most aggressive and currently untreatable cancers. Cells with MTAP deletion are vulnerable to inhibition of the metabolic enzyme, methionine-adenosyl transferase 2A (MAT2A), and the proteinarginine methyl transferase (PRMT5). This synthetic lethality has paved the way for the rapid development of drugs targeting the MAT2A/PRMT5 axis. MAT2A and its liver- and pancreas-specific isoform, MAT1A, generate the universal methyl donor S-adenosylmethionine (SAM) from ATP and methionine. Given the pleiotropic role SAM plays in methylation of diverse substrates, characterising the extent of SAM depletion and downstream perturbations following MAT2A/MAT1A inhibition (MATi)is critical for safety assessment. We have assessed in vivo target engagement and theresultant systemic phenotype using multi-omic tools to characterise response to a MAT2A inhibitor (AZ'9567). We observed significant SAM depletion and extensive methionine accumulation in the plasma, liver, brain and heart of treated rats, providing the first assessment of both global SAM depletion and evidence of hepatic MAT1A target engagement. An integrative analysis of multi-omic data from liver tissue identified broad perturbations in pathways covering one-carbon metabolism, trans-sulfuration and lipid metabolism. We infer that these pathway-wide perturbations represent adaptive responses to SAM depletion and confer a risk of oxidative stress, hepatic steatosis and an associated disturbance in plasma and cellular lipid homeostasis. The alterations also explain the dramatic increase in plasma and tissue methionine, which could be used as a safety and PD biomarker going forward to the clinic.

Knowledge-guided learning methods for integrative analysis of multi-omics data

Integrative analysis of multi-omics data has the potential to yield valuable and comprehensive insights into the molecular mechanisms underlying complex diseases such as cancer and Alzheimer's disease. However, a number of analytical challenges complicate multi-omics data integration. For instance, -omics data are usually high-dimensional, and sample sizes in multi-omics studies tend to be modest. Furthermore, when genes in an important pathway have relatively weak signal, it can be difficult to detect them individually. There is a growing body of literature on knowledge-guided learning methods that can address these challenges by incorporating biological knowledge such as functional genomics and functional proteomics into multi-omics data analysis. These methods have been shown to outperform their counterparts that do not utilize biological knowledge in tasks including prediction, feature selection, clustering, and dimension reduction. In this review, we survey recently developed methods and applications of knowledge-guided multi-omics data integration methods and discuss future research directions.

Abstract 7355: Integration analysis of lncRNAs and R-loops in 3D genome organization and gene regulation in the cancer genome

Abstract Existing evidence shows that some long non-coding RNAs (lncRNAs) can recruit chromatin regulators and form R-loops to mediate CCTCC-binding factor (CTCF) binding at topologically associated domain (TAD) boundaries to activate transcription of target genes. However, it remains unclear to what extent such lncRNA-dependent chromatin organization events occur genome-wide in normal and cancer cells. By curated public data of lncRNA, CTCF, RAD21 profiles and integrative analysis of multi-omics data using innovative computational approaches, we investigated global associations of lncRNAs and R-loops with activities of transcription factors and chromatin regulators in gene expression regulations across multiple human cell types. We developed a comprehensive computational pipeline for cross-cell-type analysis of lncRNA and protein-coding gene expression profiles and transcription factor binding profiles. We identified a total of 542 candidate genomic regions where lncRNAs play a pivotal role in regulating target genes by facilitating CTCF/cohesin binding. R-loops are enriched at CTCF/cohesin binding sites, suggesting their potential contribution to CTCF binding. Additionally, our findings underscore the significant cell-type specificity of CTCF binding patterns across the genome. This work adds knowledge to cancer genome organization and gene regulation, and builds foundations for further functional and mechanistic studies of lncRNAs and R-loops in the human genome. Citation Format: Shengyuan Wang, Zhenjia Wang, Yuh-Hwa Wang, Hui Li, Peter Todd Stukenburg, Chongzhi Zang. Integration analysis of lncRNAs and R-loops in 3D genome organization and gene regulation in the cancer genome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7355.

Integrative Analysis of Multi-Omic Data for the Characteristics of Endometrial Cancer.

Endometrial cancer (EC) is a frequently diagnosed gynecologic cancer. Identifying reliable prognostic genes for predicting EC onset is crucial for reducing patient morbidity and mortality. Here, a comprehensive strategy with transcriptomic and proteomic data was performed to measure EC's characteristics. Based on the publicly available RNA-seq data, death-associated protein kinase 3, recombination signal-binding protein for the immunoglobulin kappa J region, and myosin light chain 9 were screened out as potential biomarkers that affect the EC patients' prognosis. A linear model was further constructed by multivariate Cox regression for the prediction of the risk of being malignant. From further integrative analysis, exosomes were found to have a highly enriched role that might participate in EC occurrence. The findings were validated by qRT-polymerase chain reaction (PCR) and western blotting. Collectively, we constructed a prognostic-gene-based model for EC prediction and found that exosomes participate in EC incidents, revealing significantly promising support for the diagnosis of EC.

Multi-omics Combined with Machine Learning Facilitating the Diagnosis of Gastric Cancer.

Gastric cancer (GC) is a highly intricate gastrointestinal malignancy. Early detection of gastric cancer forms the cornerstone of precision medicine. Several studies have been conducted to investigate early biomarkers of gastric cancer using genomics, transcriptomics, proteomics, and metabolomics, respectively. However, endogenous substances associated with various omics are concurrently altered during gastric cancer development. Furthermore, environmental exposures and family history can also induce modifications in endogenous substances. Therefore, in this study, we primarily investigated alterations in DNA mutation, DNA methylation, mRNA, lncRNA, miRNA, circRNA, and protein, as well as glucose, amino acid, nucleotide, and lipid metabolism levels in the context of GC development, employing genomics, transcriptomics, proteomics, and metabolomics. Additionally, we elucidate the impact of exposure factors, including HP, EBV, nitrosamines, smoking, alcohol consumption, and family history, on diagnostic biomarkers of gastric cancer. Lastly, we provide a summary of the application of machine learning in integrating multi-omics data. Thus, this review aims to elucidate: i) the biomarkers of gastric cancer related to genomics, transcriptomics, proteomics, and metabolomics; ii) the influence of environmental exposure and family history on multi- -omics data; iii) the integrated analysis of multi-omics data using machine learning techniques.

METTL16 promotes liver cancer stem cell self-renewal via controlling ribosome biogenesis and mRNA translation

BackgroundWhile liver cancer stem cells (CSCs) play a crucial role in hepatocellular carcinoma (HCC) initiation, progression, recurrence, and treatment resistance, the mechanism underlying liver CSC self-renewal remains elusive. We aim to characterize the role of Methyltransferase 16 (METTL16), a recently identified RNA N6-methyladenosine (m6A) methyltransferase, in HCC development/maintenance, CSC stemness, as well as normal hepatogenesis.MethodsLiver-specific Mettl16 conditional KO (cKO) mice were generated to assess its role in HCC pathogenesis and normal hepatogenesis. Hydrodynamic tail-vein injection (HDTVi)-induced de novo hepatocarcinogenesis and xenograft models were utilized to determine the role of METTL16 in HCC initiation and progression. A limiting dilution assay was utilized to evaluate CSC frequency. Functionally essential targets were revealed via integrative analysis of multi-omics data, including RNA-seq, RNA immunoprecipitation (RIP)-seq, and ribosome profiling.ResultsMETTL16 is highly expressed in liver CSCs and its depletion dramatically decreased CSC frequency in vitro and in vivo. Mettl16 KO significantly attenuated HCC initiation and progression, yet only slightly influenced normal hepatogenesis. Mechanistic studies, including high-throughput sequencing, unveiled METTL16 as a key regulator of ribosomal RNA (rRNA) maturation and mRNA translation and identified eukaryotic translation initiation factor 3 subunit a (eIF3a) transcript as a bona-fide target of METTL16 in HCC. In addition, the functionally essential regions of METTL16 were revealed by CRISPR gene tiling scan, which will pave the way for the development of potential inhibitor(s).ConclusionsOur findings highlight the crucial oncogenic role of METTL16 in promoting HCC pathogenesis and enhancing liver CSC self-renewal through augmenting mRNA translation efficiency.

Multi-omics integration with weighted affinity and self-diffusion applied for cancer subtypes identification

BackgroundCharacterizing cancer molecular subtypes is crucial for improving prognosis and individualized treatment. Integrative analysis of multi-omics data has become an important approach for disease subtyping, yielding better understanding of the complex biology. Current multi-omics integration tools and methods for cancer subtyping often suffer challenges of high computational efficiency as well as the problem of weight assignment on data types.ResultsHere, we present an efficient multi-omics integration via weighted affinity and self-diffusion (MOSD) to dissect cancer heterogeneity. MOSD first construct local scaling affinity on each data type and then integrate all affinities by weighted linear combination, followed by the self-diffusion to further improve the patients’ similarities for the downstream clustering analysis. To demonstrate the effectiveness and usefulness for cancer subtyping, we apply MOSD across ten cancer types with three measurements (Gene expression, DNA methylation, miRNA).ConclusionsOur approach exhibits more significant differences in patient survival and computationally efficient benchmarking against several state-of-art integration methods and the identified molecular subtypes reveal strongly biological interpretability. The code as well as its implementation are available in GitHub: https://github.com/DXCODEE/MOSD.

Machine learning combining multi-omics data and network algorithms identifies adrenocortical carcinoma prognostic biomarkers.

Background: Rare endocrine cancers such as Adrenocortical Carcinoma (ACC) present a serious diagnostic and prognostication challenge. The knowledge about ACC pathogenesis is incomplete, and patients have limited therapeutic options. Identification of molecular drivers and effective biomarkers is required for timely diagnosis of the disease and stratify patients to offer the most beneficial treatments. In this study we demonstrate how machine learning methods integrating multi-omics data, in combination with system biology tools, can contribute to the identification of new prognostic biomarkers for ACC. Methods: ACC gene expression and DNA methylation datasets were downloaded from the Xena Browser (GDC TCGA Adrenocortical Carcinoma cohort). A highly correlated multi-omics signature discriminating groups of samples was identified with the data integration analysis for biomarker discovery using latent components (DIABLO) method. Additional regulators of the identified signature were discovered using Clarivate CBDD (Computational Biology for Drug Discovery) network propagation and hidden nodes algorithms on a curated network of molecular interactions (MetaBase™). The discriminative power of the multi-omics signature and their regulators was delineated by training a random forest classifier using 55 samples, by employing a 10-fold cross validation with five iterations. The prognostic value of the identified biomarkers was further assessed on an external ACC dataset obtained from GEO (GSE49280) using the Kaplan-Meier estimator method. An optimal prognostic signature was finally derived using the stepwise Akaike Information Criterion (AIC) that allowed categorization of samples into high and low-risk groups. Results: A multi-omics signature including genes, micro RNA's and methylation sites was generated. Systems biology tools identified additional genes regulating the features included in the multi-omics signature. RNA-seq, miRNA-seq and DNA methylation sets of features revealed a high power to classify patients from stages I-II and stages III-IV, outperforming previously identified prognostic biomarkers. Using an independent dataset, associations of the genes included in the signature with Overall Survival (OS) data demonstrated that patients with differential expression levels of 8 genes and 4 micro RNA's showed a statistically significant decrease in OS. We also found an independent prognostic signature for ACC with potential use in clinical practice, combining 9-gene/micro RNA features, that successfully predicted high-risk ACC cancer patients. Conclusion: Machine learning and integrative analysis of multi-omics data, in combination with Clarivate CBDD systems biology tools, identified a set of biomarkers with high prognostic value for ACC disease. Multi-omics data is a promising resource for the identification of drivers and new prognostic biomarkers in rare diseases that could be used in clinical practice.

Human milk miRNAs associate to maternal dietary nutrients, milk microbiota, infant gut microbiota and growth

Maternal diet influences the milk composition, yet little information is available on the impact of maternal diet on milk miRNAs expression. Further, the association of human milk miRNAs to maternal diet and milk microbiota is not explored. In addition, the role of milk miRNAs on the infant gut microbiota, infant growth and development has not been investigated. Milk samples were collected from 60 healthy lactating women at ≤15d post-partum, HTG transcriptome assay was performed to examine milk miRNA profile. Maternal clinical and dietary clusters information were available and infant anthropometric measures were followed up to one year of age. Milk and infant microbiota were analyzed by 16S rRNA gene sequencing and integrative multi-omics data analysis was performed to identify potential association between microRNA, maternal dietary nutrients and microbiota. Discriminant analysis revealed that the milk miRNAs were clustered into groups according to the maternal protein source. Interestingly, 31 miRNAs were differentially expressed (P adj<0.05) between maternal dietary clusters (Cluster 1: enriched in plant protein and fibers and Cluster 2: enriched in animal protein), with 30 miRNAs downregulated in the plant protein group relative to animal protein group. Pathway analysis revealed that the top enriched pathways (P adj<0.01) were involved in cell growth and proliferation processes. Furthermore, significant features contributing to the clustering were associated with maternal dietary nutrients and milk microbiota (r>0.70). Further, miR-378 and 320 family miRNAs involved in adipogenesis were positively correlated to the infant BMI-z-scores, weight, and weight for length-z-scores at 6 months of age. Maternal dietary source impacts the milk miRNA expression profile. Further, miRNAs were associated with maternal dietary nutrients, milk microbiota and to the infant gut microbiota and infant growth and development. The study is registered in ClinicalTrials.gov. The identification number is NCT03552939.

Multi-omics approach for identifying CNV-associated lncRNA signatures with prognostic value in prostate cancer

BackgroundProstate cancer, the second most prevalent malignancy among men, poses a significant threat to affected patients’ well-being due to its poor prognosis. Novel biomarkers are required to enhance clinical outcomes and tailor personalized treatments. Herein, we describe our research to explore the prognostic value of long non-coding RNAs (lncRNAs) deregulated by copy number variations (CNVs) in prostate cancer. MethodsThe study employed an integrative multi-omics data analysis of the prostate cancer transcriptomic, CNV and methylation datasets to identify prognosis-related subtypes. Subtype-specific expression profiles of protein-coding genes (PCGs) and lncRNAs were determined. We analysed CNV patterns of lncRNAs across the genome to identify subtype-specific lncRNAs with CNV changes. LncRNAs exhibiting significant amplification or deletion and a positive correlation were designated CNV-deregulated lncRNAs. A prognostic risk score model was subsequently developed using these CNV-driven lncRNAs. ResultsSix molecular subtypes of prostate cancer were identified, demonstrating significant differences in prognosis (P = 0.034). The CNV profiles of subtype-specific lncRNAs were examined, revealing their correlation with CNV amplification or deletion. Six lncRNAs (CCAT2, LINC01593, LINC00276, GACAT2, LINC00457, LINC01343) were selected based on significant CNV amplifications or deletions using a rigorous univariate Cox proportional risk regression model. A robust risk score model was developed, stratifying patients into high-risk and low-risk categories. Notably, our prognostic model based on these six lncRNAs exhibited exceptional predictive capabilities for recurrence-free survival (RFS) in prostate cancer patients (P = 0.024). ConclusionsOur study successfully identified a prognostic risk score model comprising six CNV-driven lncRNAs that could potentially be prognostic biomarkers for prostate cancer. These lncRNA signatures are closely associated with RFS, providing promising prospects for improved patient prognostication and personalized therapeutic strategies for novel prostate cancer treatment.

Prioritization of risk genes in colorectal cancer by integrative analysis of multi-omics data and gene networks.

Genome-wide association studies (GWASs) have identified over 140 colorectal cancer (CRC)-associated loci; however, target genes at the majority of loci and underlying molecular mechanisms are poorly understood. Here, we utilized a Bayesian approach, integrative risk gene selector (iRIGS), to prioritize risk genes at CRC GWAS loci by integrating multi-omics data. As a result, a total of 105 high-confidence risk genes (HRGs) were identified, which exhibited strong gene dependencies for CRC and enrichment in the biological processes implicated in CRC. Among the 105 HRGs, CEBPB, located at the 20q13.13 locus, acted as a transcription factor playing critical roles in cancer. Our subsequent assays indicated the tumor promoter function of CEBPB that facilitated CRC cell proliferation by regulating multiple oncogenic pathways such as MAPK, PI3K-Akt, and Ras signaling. Next, by integrating a fine-mapping analysis and three independent case-control studies in Chinese populations consisting of 8,039 cases and 12,775 controls, we elucidated that rs1810503, a putative functional variant regulating CEBPB, was associated with CRC risk (OR=0.90, 95%CI=0.86-0.93, P=1.07×10-7). The association between rs1810503 and CRC risk was further validated in three additional multi-ancestry populations consisting of 24,254 cases and 58,741 controls. Mechanistically, the rs1810503 A to T allele change weakened the enhancer activity in an allele-specific manner to decrease CEBPB expression via long-range promoter-enhancer interactions, mediated by the transcription factor, REST, and thus decreased CRC risk. In summary, our study provides a genetic resource and a generalizable strategy for CRC etiology investigation, and highlights the biological implications of CEBPB in CRC tumorigenesis, shedding new light on the etiology of CRC.

Integrative Analysis of Multi-Omics Data with Deep Learning: Challenges and Opportunities in Bioinformatics.

The advent of high-throughput technologies has ushered in an era of unprecedented data generation in the field of bioinformatics. Omics data, including genomics, transcriptomics, proteomics, and metabolomics, provide comprehensive insights into biological systems, but their integration poses significant challenges. Integrative analysis of multi-omics data holds the promise of unraveling complex biological phenomena and enabling personalized medicine. [1] Deep learning, a subset of machine learning, has gained prominence in bioinformatics due to its ability to automatically extract intricate patterns from large-scale multi-omics datasets. This paper presents an overview of the challenges and opportunities associated with the integrative analysis of multi-omics data using deep learning techniques in bioinformatics.The challenges in multi-omics integration primarily stem from data heterogeneity, dimensionality, and noise. One of the key opportunities presented by deep learning is its ability to capture complex, non-linear relationships in multi-omics data. The paper emphasizes the importance of interpretability and explainability in deep learning models applied to bioinformatics, as they play a crucial role in gaining biological insights and facilitating clinical decision-making. The integration of domain knowledge and biological context is highlighted as a critical aspect of model development. The paper showcases real-world applications of deep learning in multi-omics data integration, such as disease subtype classification, biomarker discovery, and drug response prediction. As the field continues to evolve, addressing these challenges and harnessing the potential of deep learning approaches will pave the way for transformative advancements in our understanding of complex biological systems and the development of precision medicine strategies.

Integrative analysis of multi-omics data reveals a pseudouridine-related lncRNA signature for prediction of glioma prognosis and chemoradiotherapy sensitivity

BackgroundGlioblastoma is the most common type of glioma with a high incidence and poor prognosis, and effective medical treatment remains challenging. Pseudouridine (Ψ) is the first post-transcriptional modification discovered and one of the most abundant modifications to RNA. However, the prognostic value of Ψ-related lncRNAs (ΨrLs) for glioma patients has never been systematically evaluated. This study aims to construct a risk model based on ΨrLs signature and to validate the predictive efficiency of the model. MethodTranscriptomic data, genomic data, and relevant clinical data of glioma patients were extracted from the Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). ΨrLs with significant correlation with Ψ-related genes were identified, and univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression, and multivariate Cox regression were used to further select biomarkers and construct a ΨrLs signature risk model. Then, the expression of lncRNAs of ΨrLs signature in multiple glioma cell lines was detected by qPCR. Further, ROC analysis, stratification analysis, correlation analysis, survival analysis, nomogram, enrichment analysis, immune infiltration analysis, chemoradiotherapy sensitivity analysis, somatic mutation, and recurrent copy number variation (CNV) analysis were used to validate the predictive efficiency of ΨrLs signature in TCGA and CGGA datasets. ResultsA four-lncRNA ΨrLs signature (DNAJC27-AS1, GDNF-AS1, ZBTB20-AS4, and DNMBP-AS1) risk model was constructed. By ROC analysis, stratified analysis, correlation analysis, survival analysis, and nomogram, the signature showed satisfactory predictive efficiency. Functional enrichment analysis revealed the differences in immune-related biological processes between high- and low-risk groups. Immune infiltration analysis showed that the high-risk group had lower tumor purity and higher stromal, immune and ESTIMATE scores. Mitoxantrone was identified as effective drug for low-risk group of glioma patients. Key genes in glioma development, including IDH1, EGFR, PTEN, etc., were differentially mutated between risk groups. The main recurrent CNVs in low-risk groups were 19q13.42 deletion and 7q34 amplification; 10q23.31 deletion and 12q14.1 in the high-risk group. ConclusionsOur study identified a four-lncRNA ΨrLs signature that effectively predicts the prognosis of glioma patients and may serve as a diagnostic tool. Risk scores of glioma patients generated by the signature is associated with immune-related biological processes and chemoradiotherapy sensitivity. These findings may inform the development of more targeted and effective therapies for glioma patients.

Abstract P225: The Reninness Score: Integrative Analysis Of Multi-omic Data To Define Renin Cell Identity.

Renin cell (RC) descendants regain the endocrine renin phenotype to overcome homeostatic threats, a process known as recruitment. The determinants that control the identity, fate, and recruitment of RCs are not well understood. We aim to 1) define the chromatin pattern that determines RC identity; and 2) develop a computational tool to score samples by similarity to that pattern. We compared open chromatin regions (ATAC-seq) between non-renin-expressing (n = 184) and renin-expressing samples (n = 19) from three different sources and stimulation states: 1) primary RCs in basal state (WT); 2) primary RCs chronically stimulated to produce renin (recruited); and 3) constitutively renin-expressing tumoral cells. Differential analyses between the three types of RCs and the non-RCs revealed 1,525 unique open chromatin regions shared in at least two RC groups, including regions associated with RC-specific genes: Ren1 and Akr1b7. Increased accessibility regions in RC groups are enriched in renin-related GO terms ( e.g., cAMP and Notch signaling pathway). The bZip family ( e.g. , Atf3, p &lt; 10 -2147 ) was the most enriched motifs in tumoral cells, and MEF2 family ( e.g. , MEF2a, p &lt; 10 -91 ) in primary cells. Co-enriched motifs (e.g., MEIS1, p &lt; 10 -10 ; NR4A2, p &lt; 10 -9 ) in Ren1 and Akr1b7 regions suggest a group of TFs govern their co-expression. In addition, we identified unique open chromatin regions and enriched motifs ( e.g., E2F4, p &lt; 10 -9 ; SP2, p &lt; 10 -7 ) in recruited RCs, suggesting associated TFs involved in the recruitment.Next, we used a machine-learning approach to calculate a "reninness" score by training a genomic region model with renin and non-renin ATAC-seq data. We identified a cluster of regions (n = 112,992) around the trained renin label, which included 89% (21,820 of 24,612) of the increased accessibility regions from the differential analyses. We calculated the reninness score based on the distance between the trained label embeddings and predicted samples using the trained model. Homogeneous WT RCs had the highest score, followed by renin-expressing tumoral and primary cells. Non-RCs had the lowest score. In conclusion, we identified the chromatin configurations defining the identity and recruitment of RCs and developed a computational score to quantify it.

Multi-omics analysis revealing the interplay between gut microbiome and the host following opioid use

ABSTRACT Opioid crisis is an ongoing epidemic since the past several decades in the United States. Opioid use-associated microbial dysbiosis is emerging as a key regulator of intestinal homeostasis and behavioral responses to opioid. However, the mechanistic insight into the role of microbial community in modulating host response is unavailable. To uncover the role of opioid-induced dysbiosis in disrupting intestinal homeostasis we utilized whole genome sequencing, untargeted metabolomics, and mRNA sequencing to identify changes in microbiome, metabolome, and host transcriptome respectively. Morphine treatment resulted in significant expansion of Parasuterella excrementihominis, Burkholderiales bacterium 1_1_47, Enterococcus faecalis, Enterorhabdus caecimuris and depletion of Lactobacillus johnsonii. These changes correlated with alterations in lipid metabolites and flavonoids. Significant alteration in microbial metabolism (metabolism of lipids, amino acids, vitamins and cofactors) and increased expression of virulence factors and biosynthesis of lipopolysaccharides (LPS) and lipoteichoic acid (LTA) were observed in microbiome of morphine-treated animals. In concurrence with changes in microbiome and metabolome extensive changes in innate and adaptive immune response, lipid metabolism, and gut barrier dysfunction were observed in the host transcriptome. Microbiome depleted mice displayed lower levels of inflammation, immune response and tissue destruction compared to mice harboring a dysbiotic microbiome in response to morphine treatment, thus establishing dysbiotic microbiome as mediator of morphine gut pathophysiology. Integrative analysis of multi-omics data highlighted the associations between Parasutterella excrementihominis, Burkholderiales bacterium 1_1_47, Enterococcus faecalis, Enterorhabdus caecimuris and altered levels of riboflavin, flavonoids, and lipid metabolites including phosphocholines, carnitines, bile acids, and ethanolamines with host gene expression changes involved in inflammation and barrier integrity of intestine. Omic analysis also highlighted the role of probiotic bacteria Lactobacillus johnsonii, metabolites flavonoids and riboflavin that were depleted with morphine as important factors for intestinal homeostasis. This study presents for the first time ever an interactive view of morphine-induced changes in microbial metabolism, strain level gut microbiome analysis and comprehensive view of changes in gut transcriptome. We also identified areas of potential therapeutic interventions to limit microbial dysbiosis and present a unique resource to the opioid research community.

Integrated multiomics analysis of chromosome 19 miRNA cluster in bladder cancer

With 46 microRNAs (miRNAs) embedded tandemly over a distance of ~100 kb, chromosome 19 microRNA cluster (C19MC) is the largest miRNA cluster in the human genome. The C19MC is transcribed from a long noncoding genomic region and is usually expressed simultaneously at a higher level. Hence, we performed an integrative multiomics data analysis to examine C19MC regulation, expression patterns, and their impact on bladder cancer (BCa). We found that 43 members of C19MC were highly expressed in BCa. However, its co-localization with recurrent copy number variation (CNV) gain was not statistically significant to implicate its upregulation. It has been reported that C19MC expression is regulated by a well-established CpG island situated 17.6 kb upstream of the transcription start site, but we found that CpG probes at this island were hypomethylated, which was not statistically significant in the BCa cohort. In addition, the promoter region of C19MC is strongly regulated by a group of seven transcription factors (NR2F6, SREBF1, TBP, GATA3, GABPB1, ETV4, and ZNF444) and five chromatin modifiers (SMC3, KDMA1, EZH2, RAD21, and CHD7). Interestingly, these 12 genes were found to be overexpressed in BCa patients. Further, C19MC targeted 42 tumor suppressor (TS) genes that were downregulated, of which 15 were significantly correlated with patient survival. Our findings suggest that transcription factors and chromatin modifiers at the promoter region may regulate C19MC overexpression. The upregulated C19MC members, transcription regulators, and TS genes can be further exploited as potential diagnostic and prognostic indicators as well as for therapeutic management of BCa.

Identification of USP29 as a key regulator of nucleotide biosynthesis in neuroblastoma through integrative analysis of multi-omics data

ABSTRACT Cancer cells show enhanced nucleotide biosynthesis, which is essential for their unlimited proliferation, but the underlying mechanisms are not entirely clear. Ubiquitin specific peptidase 29 (USP29) was reported to sustain neuroblastoma progression by promoting glycolysis and glutamine catabolism; however, its potential role in regulating nucleotide biosynthesis in tumor cells remains unknown. In this study, we depleted endogenous USP29 in MYCN-amplified neuroblastoma SK-N-BE2 cells by sgRNAs and conducted metabolomic analysis in cells with or without USP29 depletion, we found that USP29 deficiency caused a disorder of intermediates involved in glycolysis and nucleotide biosynthesis. De novo nucleotide biosynthesis was analyzed using 13C6 glucose as a tracer under normoxia and hypoxia. The results indicated that USP29-depleted cells showed inhibition of nucleotide anabolic intermediates derived from glucose, and this inhibition was more significant under hypoxic conditions. Analysis of RNA sequencing data in SK-N-BE2 cells demonstrated that USP29 promoted the gene expression of metabolic enzymes involved in nucleotide anabolism, probably by regulating MYC and E2F downstream pathways. These findings indicated that USP29 is a key regulator of nucleotide biosynthesis in tumor cells.