Integrative single-cell and spatial transcriptomic analyses identify a pathogenic cholangiocyte niche and TNFRSF12A as therapeutic target for biliary atresia

Abstract

Background & Aims: Biliary Atresia (BA) is a devastating fibro-inflammatory biliary disease that is the leading indication for pediatric liver transplants worldwide. Although cholangiocytes are the primary target cells, the pathogenic mechanisms involving cholangiocytes remain elusive. Here, we aimed to characterize the pathogenic role of cholangiocytes in BA. Approach & Results: Integration of single-cell RNA sequencing of 12 liver tissues (from 9 BA and 3 controls) and the spatial transcriptome of another four liver sections (from 2 BA and 2 controls) provided a comprehensive spatial liver cell atlas of BA. In particular, we identified a cholangiocyte-enriched spatial niche with infiltration of activated hepatic stellate cells, activated portal fibroblasts, macrovascular endothelial cells and TREM2 + macrophages that were elevated in the portal triad of BA. This niche was positively correlated with bile duct profiles, liver fibrosis and poor survival in two independent cohorts of patients with BA. Using integrative bioinformatics analyses to mine the cell-cell communication and regulatory network in BA cholangiocytes, we uncovered the fibro-inflammatory phenotype of cholangiocytes with TNFSF12-TNFRSF12A as a significant signal. Genetic ablation or blockade of TNFRSF12A suppresses liver injury, inflammation, and bile duct profiles in a mouse model of disease. Using human biliary organoids, we revealed that BA organoids expressed higher levels of CCL2 in response to TNFSF12 stimulation and promoted monocyte chemotaxis via the CCL2-CCR2 axis. Conclusions: Pathogenic cholangiocytes enriched niche identifies TNFRSF12A as a potential therapeutic target for BA.