Identifying and validating molecular subtypes of biliary atresia using multiple high-throughput data integration analysis.

Background: Biliary atresia (BA) is the most common form of severe neonatal obstructive jaundice. The etiology and pathogenesis of BA are multifactorial, and different factors may interact to produce heterogeneous pathological features and clinical outcomes. Despite different pathological features, all patients received the same treatment strategy. This study performed integrative clustering analysis based on multiple high-throughput datasets to identify the molecular subtypes of BA and provide a new treatment strategy for personalized treatment of the different subtypes of BA.

Methods: The RNA sequence dataset GSE122340 in the Gene Expression Omnibus (GEO) database was downloaded; 31 BA and 20 control normal liver tissues were collected at our center for transcriptome sequencing, and clinical and follow-up data of BA patients were available. Molecular subtypes were identified using integrated unsupervised cluster analysis involving gene expression, biliary fibrosis, and immune enrichment scores based on the transcriptome dataset, and the results were validated using independent datasets.

Results: Based on the results of the integrated unsupervised clustering analysis, four molecular subtypes were identified: autoimmune, inflammatory, virus infection-related, and oxidative stress. The autoimmune subtype with a moderate prognosis was dominated by autoimmune responses and morphogenesis, such as the Fc-gamma receptor and Wnt signaling pathway. The biological process of the inflammatory subtype was mainly the inflammatory response, with the best prognosis, youngest age at surgery, and lowest liver stiffness. The virus infection-related subtype had the worst prognosis and was enriched for a variety of biological processes such as viral infection, immunity, anatomical morphogenesis, and epithelial mesenchymal transition. The oxidative stress subtype was characterized by the activation of oxidative stress and various metabolic pathways and had a poor prognosis. The above results were verified independently in the validation sets.

Conclusions: This study identified four molecular subtypes of BA with distinct prognosis and biological processes. According to the pathological characteristics of the different subtypes, individualized perioperative and preoperative treatment may be a new strategy to improve the prognosis of BA.