Genome-wide transcriptomics and a machine learning-based random forest model identify novel biomarkers to predict phenotypic variability in Wilson disease

Genotype-phenotype correlation studies fail to explain the phenotypic diversity in Wilson disease (WD) and the missing link may lie in transcriptomics. Therefore, we performed genome-wide transcriptomic profiling and machine-learning-based random forest modeling (RF) to identify the most important genes that could explain phenotypic variability in WD and reveal the molecular pathogenesis of WD. For this purpose, we used bulk RNA sequencing approaches to include 47 genetically confirmed pediatric hepatic WD patients and controls. We observed 321 differentially expressed genes (DEGs) in WD. Among top 100 DEGs, 42 significant DEGs (p 2FC = 4.68) and downregulation of DUSP16 (log2FC = -5.22), and GAPLINC (log2FC = -2.45). Iron homeostasis signaling demonstrated comparable enrichment, involving key genes ATP6V0D2 (log2FC = -4.06), HMOX2 (log2FC = -2.70), HBD (log2FC = 1.96) and HBZ (log2FC = 1.76). These results suggest that children with WD have dysregulation in MAP Kinase and iron homeostasis signaling pathways providing an important insight into the molecular pathogenesis of WD. RF modeling identified a panel of biomarker with ten genes, that could explain the phenotypic variability in WD.