Multi-omics approaches for drug-response characterization in primary biliary cholangitis and autoimmune hepatitis variant syndrome

Abstract Background Primary biliary cholangitis (PBC) and autoimmune hepatitis (AIH) variant syndrome (VS) exhibit a complex overlap of AIH features with PBC, leading to poorer prognoses than those with PBC or AIH alone. The biomarkers associated with drug response and potential molecular mechanisms in this syndrome have not been fully elucidated. Methods Whole-transcriptome sequencing was employed to discern differentially expressed (DE) RNAs within good responders (GR) and poor responders (PR) among patients with PBC/AIH VS. Subsequent gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted for the identified DE RNAs. Plasma metabolomics was employed to delineate the metabolic profiles distinguishing PR and GR groups. The quantification of immune cell profiles and associated cytokines was achieved through flow cytometry and immunoassay technology. Uni- and multivariable logistic regression analyses were conducted to construct a predictive model for insufficient biochemical response. The performance of the model was assessed by computing the area under the receiver operating characteristic (AUC) curve, sensitivity, and specificity. Findings The analysis identified 224 differentially expressed (DE) mRNAs, 189 DE long non-coding RNAs, 39 DE circular RNAs, and 63 DE microRNAs. Functional pathway analysis revealed enrichment in lipid metabolic pathways and immune response. Metabolomics disclosed dysregulated lipid metabolism and identified PC (18:2/18:2) and PC (16:0/20:3) as predictors. CD4+ T helper (Th) cells, including Th2 cells and regulatory T cells (Tregs), were upregulated in the GR group. Pro-inflammatory cytokines (IFN-γ, TNF-α, IL-9, and IL-17) were downregulated in the GR group, while anti-inflammatory cytokines (IL-10, IL-4, IL-5, and IL-22) were elevated. Regulatory networks were constructed, identifying CACNA1H and ACAA1 as target genes. A predictive model based on these indicators demonstrated an AUC of 0.986 in the primary cohort and an AUC of 0.940 in the validation cohort for predicting complete biochemical response. Conclusion A combined model integrating genomic, metabolic, and cytokinomic features demonstrated high accuracy in predicting insufficient biochemical response in patients with PBC/AIH VS. Early recognition of individuals at elevated risk for insufficient response allows for the prompt initiation of additional treatments.