Multiomics data integration is one of the leading frontiers of complex disease research and integrative biology. The advances in single-cell sequencing technologies offer yet another crucial dimension in multiomics research. The single-cell studies enable the study and integration of multiomics data simultaneously in the same cell. We report in this study multiomics data integration in single-cell resolution using Bayesian networks (BNs) in a case study of hepatocellular carcinoma (HCC). A BN encodes the conditional dependencies/independencies of variables using a graphical model with an accompanying joint probability. RNA-seq and Reduced Representation Bisulfite Sequencing data were analyzed separately, and copy number variations were estimated by the hidden Markov model method. Several BN models were constructed to reveal omics' causal and associational relationships. These methods were subjected to a validation study using an independent data set. We show the heterogeneity of the multiple cellular layers of HCC at single-cell omics resolution by identifying best-fitted BN models of 295 genes. We also provide novel insights into the multiomics mechanistic relationships in the human lymphocyte antigen class I genes in HCC. To the best of our knowledge, this is the first study to focus on integrating omics data using a machine learning algorithm, BNs, at the single-cell resolution using a case study of HCC.
BibikovaM, LeJ, BarnesB, et al.Genome-wide DNA methylation profiling using Infinium® assay. Epigenomics, 2009; 1(1):177–200; doi: 10.2217/epi.09.14
3.
BolgerAM, LohseM, UsadelB. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics, 2014; 30(15):2114–2120; doi: 10.1093/bioinformatics/btu170
4.
CadwellCR, ScalaF, LiS, et al.Multimodal profiling of single-cell morphology, electrophysiology, and gene expression using Patch-Seq. Nat Protoc, 2017; 12(12):2531–2553; doi: 10.1038/nprot.2017.120
5.
CeramiE, GaoJ, DogrusozU, et al.The CBio Cancer Genomics Portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov, 2012; 2(5):401–404; doi: 10.1158/2159-8290.CD-12-0095
6.
ChenJ, LiuW, DuJ, et al.Comprehensive genomic and epigenomic analyses on transcriptomic regulation in stomach adenocarcinoma. Front Genet, 2021; 12:778095; doi: 10.3389/fgene.2021.778095
7.
ChenW, ChengP, JiangJ, et al.Epigenomic and genomic analysis of transcriptome modulation in skin cutaneous melanoma. Aging (Albany NY), 2020; 12(13):12703; doi: 10.18632/aging.103115
8.
ColagrandeS, InghilesiAL, AburasS, et al.Challenges of advanced hepatocellular carcinoma. World J Gastroenterol, 2016; 22(34):7645; doi: 10.3748/wjg.v22.i34.7645
9.
Díez-VillanuevaA, JordáM, Carreras-TorresR, et al.Identifying causal models between genetically regulated methylation patterns and gene expression in healthy colon tissue. Clin Epigenetics, 2021; 13(1):162; doi: 10.1186/s13148-021-01148-9
10.
DohmenJ, BaranovskiiA, RonenJ, et al.Identifying tumor cells at the single-cell level using machine learning. Genome Biol, 2022; 23(1):1–23; doi: 10.1186/s13059-022-02683-1
11.
Dongfang WangJG.Integrative clustering methods of multi-omics data for molecule-based cancer classifications. Quant Biol, 2016; 4(1):58–67; doi: 10.1007/s40484-016-0063-4
12.
DrouardG, OllikainenM, MykkänenJ, et al.Multi-omics integration in a twin cohort and predictive modeling of blood pressure values. OMICS, 2022; 26(3):130–141; doi: 10.1089/omi.2021.0201
13.
El-ManzalawyY, HsiehTY, ShivakumarM, et al.Min-redundancy and max-relevance multi-view feature selection for predicting ovarian cancer survival using multi-omics data. BMC Med Genomics, 2018; 11(Suppl 3):71; doi: 10.1186/s12920-018-0388-0
14.
EstellerM.Cancer Epigenomics: DNA methylomes and histone-modification maps. Nat Rev Genet, 2007; 8(4):286–298; doi: 10.1038/nrg2005
15.
FidlerIJ.Commentary on “Tumor Heterogeneity and the Biology of Cancer Invasion and Metastasis.”. Cancer Res, 2016; 76(12):3441–3442; doi: 10.1158/0008-5472.CAN-16-1330
16.
FrauM, BiasiF, FeoF, et al.Prognostic markers and putative therapeutic targets for hepatocellular carcinoma. Mol Aspects Med, 2010; 31(2):179–193; doi: 10.1016/j.mam.2010.02.007
17.
GuanXY, FangY, ShamJS, et al.Recurrent chromosome alterations in hepatocellular carcinoma detected by comparative genomic hybridization. Genes Chromosomes Cancer, 2000; 29(2):110–116; doi: 10.1002/1098-2264(2000)9999:9999<::AID-GCC1022>3.0.CO;2-V
18.
Gutierrez-ArcelusM, LappalainenT, MontgomerySB, et al.Passive and active DNA methylation and the interplay with genetic variation in gene regulation. Elife, 2013; 2:e00523; doi: 10.7554/eLife.00523
19.
HanagiriT, ShigematsuY, KurodaK, et al.Prognostic implications of human leukocyte antigen class I expression in patients who underwent surgical resection for non-small-cell lung cancer. J Surg Res, 2013; 181(2):e57–e63; doi: 10.1016/j.jss.2012.07.029
20.
HanahanD, WeinbergRA. Hallmarks of cancer: The next Generation. Cell, 2011; 144(5):646–674; doi: 10.1016/j.cell.2011.02.013
HeS, TangS. WNT/beta-catenin signaling in the development of liver cancers. Biomed Pharmacother, 2020; 132:110851; doi: 10.1016/j.biopha.2020.110851
23.
HeoH, KimJ-H, LimHJ, et al.DNA methylome and single-cell transcriptome analyses reveal CDA as a potential druggable target for ALK inhibitor–resistant lung cancer therapy. Exp Mol Med, 2022; 54(8):1236–1249; doi: 10.1038/s12276-022-00836-7
24.
HouY, GuoH, CaoC, et al.Single-cell triple omics sequencing reveals genetic, epigenetic, and transcriptomic heterogeneity in hepatocellular carcinomas. Cell Res, 2016; 26(3):304–319; doi: 10.1038/cr.2016.23
25.
HoweyR, ShinSY, ReltonC, et al.Bayesian network analysis incorporating genetic anchors complements conventional Mendelian randomization approaches for exploratory analysis of causal relationships in complex data. PLoS Genet, 2020; 16(3):e1008198; doi: 10.1371/journal.pgen.1008198
26.
HuangP, XuM, HanH, et al.Integrative analysis of epigenome and transcriptome data reveals aberrantly methylated promoters and enhancers in hepatocellular carcinoma. Front Oncol, 2021; 11:769390; doi: 10.3389/fonc.2021.769390
27.
HuangQ, LiJ, ZhengJ, et al.The carcinogenic role of the Notch signaling pathway in the development of hepatocellular carcinoma. J Cancer, 2019; 10(6):1570–1579; doi: 10.7150/jca.26847
28.
HuangS, ChaudharyK, GarmireLX. More is better: recent progress in multi-Omics data integration methods. Front Genet, 2017; 8:84; doi: 10.3389/fgene.2017.00084
29.
KanekoK, IshigamiS, KijimaY, et al.Clinical implication of HLA class I expression in breast cancer. BMC Cancer, 2011; 11:454; doi: 10.1186/1471-2407-11-454
30.
KaufmanDS, SchoonRA, RobertsonMJ, et al.Inhibition of selective signaling events in natural killer cells recognizing major histocompatibility complex class I. Proc Natl Acad Sci U S A, 1995; 92(14):6484–6488; doi: 10.1073/pnas.92.14.6484
31.
KiselevVY, KirschnerK, SchaubMT, et al.SC3: Consensus clustering of single-cell RNA-Seq data. Nat Methods, 2017; 14(5):483–486; doi: 10.1038/nmeth.4236
32.
KongL, LiuP, ZhengM, et al.Multi-omics analysis based on integrated genomics, epigenomics and transcriptomics in pancreatic cancer. Epigenomics, 2020; 12(6):507–524; doi: 10.2217/epi-2019-0374
33.
KruegerF, AndrewsSR. Bismark: A flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics, 2011; 27(11):1571–1572; doi: 10.1093/bioinformatics/btr167
34.
KruegerF, JamesF, EwelsP, et al. FelixKrueger/TrimGalore: V0.6.7-DOI via Zenodo. 2021
35.
LaiD, HaG, ShahS. HMMcopy: Copy Number Prediction with Correction for GC and Mappability Bias for HTS Data. R Package Version 1.38.0; 2022.
36.
LawsonDA, KessenbrockK, DavisRT, et al.Tumour heterogeneity and metastasis at single-cell resolution. Nat Cell Biol, 2018; 20(12):1349–1360; doi: 10.1038/s41556-018-0236-7
37.
LinB, MaY, WuS. Multi-omics and artificial intelligence-guided data integration in chronic liver disease: Prospects and challenges for precision medicine. OMICS, 2022; 26(8):415–421; doi: 10.1089/omi.2022.0079
38.
LinCP, LiuCR, LeeCN, et al.Targeting C-Myc as a novel approach for hepatocellular carcinoma. World J Hepatol, 2010; 2(1):16–20; doi: 10.4254/wjh.v2.i1.16
39.
LiuG, DongC, LiuL. Integrated multiple “-Omics” data reveal subtypes of hepatocellular carcinoma. PLoS One, 2016; 11(11):e0165457; doi: 10.1371/journal.pone.0165457
40.
LiuT, LiuC, YanM, et al.Single cell profiling of primary and paired metastatic lymph node tumors in breast cancer patients. Nat Commun, 2022; 13(1):1–17; doi: 10.1038/s41467-022-34581-2
41.
LorahJ, WomackA. Value of sample size for computation of the Bayesian information criterion (BIC) in multilevel modeling. Behav Res Methods, 2019; 51(1):440–450; doi: 10.3758/s13428-018-1188-3
42.
MadjdZ, SpendloveI, PinderSE, et al.Total loss of MHC class I is an independent indicator of good prognosis in breast cancer. Int J Cancer, 2005; 117(2):248–255; doi: 10.1002/ijc.21163
43.
MarreroJA, KudoM, BronowickiJ-P. The challenge of prognosis and staging for hepatocellular carcinoma. Oncologist, 2010; 15(S4):23–33; doi: 10.1634/theoncologist.2010-S4-23
44.
MizukiN, AndoH, KimuraM, et al.Nucleotide sequence analysis of the HLA class I region spanning the 237-Kb segment around the HLA-B and -C genes. Genomics, 1997; 42(1):55–66; doi: 10.1006/geno.1997.4708
NeubergLG.Causality: Models, reasoning, and inference, by Judea Pearl, Cambridge University Press, 2000. Econ Theory, 2003; 19(4):675–685; doi: 10.1017/S0266466603004109
47.
OnoH, AraiY, FurukawaE, et al.Single-cell DNA and RNA sequencing reveals the dynamics of intra-tumor heterogeneity in a colorectal cancer model. BMC Biol, 2021; 19(1):207; doi: 10.1186/s12915-021-01147-5
48.
PearlJ. Causality: Models, Reasoning, and Inference. Cambridge University Press: Cambridge, United Kingdom; 2000.
49.
PetraliaF, WangP, YangJ, et al.Integrative random forest for gene regulatory network inference. Bioinformatics, 2015; 31(12):i197–i205; doi: 10.1093/bioinformatics/btv268
50.
PollinoCA, HendersonC.Bayesian Networks: A Guide for Their Application in Natural Resource Management and Policy. Landscape Logic, Technical Report; 2010.
51.
PruittKD, TatusovaT, BrownGR, et al.NCBI Reference Sequences (RefSeq): Current status, new features and genome annotation policy. Nucleic Acids Res, 2012; 40(Database Issue):D130–D135; doi: 10.1093/nar/gkr1079
RamnathN, TanD, LiQ, et al.Is downregulation of MHC class I antigen expression in human non-small cell lung cancer associated with prolonged survival?. Cancer Immunol Immunother, 2006; 55(8):891–899; doi: 10.1007/s00262-005-0085-7
54.
RodinAS, BoerwinkleE. Mining genetic epidemiology data with Bayesian networks I: Bayesian networks and example application (plasma ApoE levels). Bioinformatics, 2005; 21(15):3273–3278; doi: 10.1093/bioinformatics/bti813
55.
ScutariM.Learning Bayesian networks with the Bnlearn R package. J Stat Softw, 2010; 35(3):1–22; doi: 10.18637/jss.v035.i03
56.
ShanbhogueAK, PrasadSR, TakahashiN, et al.Recent advances in cytogenetics and molecular biology of adult hepatocellular tumors: Implications for imaging and management. Radiology, 2011; 258(3):673–693; doi: 10.1148/radiol.10100376
57.
ShenY, XiongW, GuQ, et al.Multi-omics integrative analysis uncovers molecular subtypes and mRNAs as therapeutic targets for liver cancer. Front Med (Lausanne), 2021; 8:65635; doi: 10.3389/fmed.2021.654635
58.
SongH, WeinsteinHNW, AllegakoenP, et al.Single-cell analysis of human primary prostate cancer reveals the heterogeneity of tumor-associated epithelial cell states. Nat Commun, 2022; 13(1):1–20; doi: 10.1038/s41467-021-27322-4
59.
SpeetjensFM, de BruinEC, MorreauH, et al.Clinical impact of HLA class I expression in rectal cancer. Cancer Immunol Immunother, 2008; 57(5):601–609; doi: 10.1007/s00262-007-0396-y
60.
StrazzaboscoM, FabrisL. Notch signaling in hepatocellular carcinoma: Guilty in association! Gastroenterology, 2012;143(6):1430–1434; doi: 10.1053/j.gastro.2012.10.025
61.
SunB, LinG, JiD, et al.Dysfunction of sister chromatids separation promotes progression of hepatocellular carcinoma according to analysis of gene expression profiling. Front Physiol, 2018; 9:1019; doi: 10.3389/fphys.2018.01019
62.
TiniG, MarchettiL, PriamiC, et al.Multi-omics integration-a comparison of unsupervised clustering methodologies. Brief Bioinform, 2019; 20(4):1269–1279; doi: 10.1093/bib/bbx167
TownesFW, HicksSC, AryeeMJ, et al.Feature selection and dimension reduction for single-cell RNA-Seq based on a multinomial model. Genome Biol, 2019; 20(1):295; doi: 10.1186/s13059-019-1861-6
65.
TrapnellC, RobertsA, GoffL, et al.Differential gene and transcript expression analysis of RNA-Seq experiments with TopHat and cufflinks. Nat Protoc, 2012; 7(3):562–578; doi: 10.1038/nprot.2012.016
66.
Velazquez-VillarrealEI, MaheshwariS, SorensonJ, et al.Single-cell sequencing of genomic DNA resolves sub-clonal heterogeneity in a melanoma cell line. Commun Biol, 2020; 3(1):318; doi: 10.1038/s42003-020-1044-8
67.
VillanuevaA, AlsinetC, YangerK, et al.Notch signaling is activated in human hepatocellular carcinoma and induces tumor formation in mice. Gastroenterology, 2012; 143(6):1660–1669; doi: 10.1053/j.gastro.2012.09.002
68.
Waly RaphaelS, YangdeZ, YuxiangC. Hepatocellular carcinoma: Focus on different aspects of management. ISRN Oncol, 2012; 2012:421673; doi: 10.5402/2012/421673
69.
WooHG, ChoiJ-H, YoonS, et al.Integrative analysis of genomic and epigenomic regulation of the transcriptome in liver cancer. Nat Commun, 2017; 8(1):839; doi: 10.1038/s41467-017-00991-w
70.
XuC, SuZ. Identification of cell types from single-cell transcriptomes using a novel clustering method. Bioinformatics, 2015; 31(12):1974–1980; doi: 10.1093/bioinformatics/btv088
71.
XuH, ZhuX, XuZ, et al.Non-invasive analysis of genomic copy number variation in patients with hepatocellular carcinoma by next generation DNA sequencing. J Cancer, 2015; 6(3):247–253; doi: 10.7150/jca.10747
72.
YanJ, RisacherSL, ShenL, et al.Network approaches to systems biology analysis of complex disease: Integrative methods for multi-Omics data. Brief Bioinform, 2018; 19(6):1370–1381; doi: 10.1093/bib/bbx066
73.
YangS, LiuG. Targeting the Ras/Raf/MEK/ERK pathway in hepatocellular carcinoma. Oncol Lett, 2017; 13(3):1041–1047; doi: 10.3892/ol.2017.5557
74.
YuG, WangL-G, HanY, et al.ClusterProfiler: An R package for comparing biological themes among gene clusters. OMICS, 2012; 16(5):284–287; doi: 10.1089/omi.2011.0118
75.
ZarayenehN, KoE, OhJH, et al.Integration of multi-Omics data for integrative gene regulatory network inference. Int J Data Min Bioinform, 2017; 18(3):223–239; doi: 10.1504/IJDMB.2017.10008266
76.
ZhangY, LiD, FengF, et al.Progressive and prognosis value of Notch receptors and ligands in hepatocellular carcinoma: A systematic review and meta-analysis. Sci Rep, 2017; 7(1):14809; doi: 10.1038/s41598-017-14897-6
77.
ZhouC, ZhangW, ChenW, et al.Integrated analysis of copy number variations and gene expression profiling in hepatocellular carcinoma. Sci Rep, 2017; 7(1):10570; doi: 10.1038/s41598-017-11029-y
78.
ZhuC, HoY-J, SalomaoMA, et al.Notch activity characterizes a common hepatocellular carcinoma subtype with unique molecular and clinicopathologic features. J Hepatol, 2021; 74(3):613–626; doi: 10.1016/j.jhep.2020.09.032
79.
ZimonjicDB, PopescuNC. Role of DLC1 tumor suppressor gene and MYC oncogene in pathogenesis of human hepatocellular carcinoma: Potential prospects for combined targeted therapeutics (review). Int J Oncol, 2012; 41(2):393–406; doi: 10.3892/ijo.2012.1474
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
