Society: AASLD

Background: Chronic cholestatic liver diseases, such as primary sclerosing cholangitis (PSC), often lead to end-stage liver disease requiring liver transplantation. PSC also significantly increases the risk of cholangiocarcinoma (CCA). There is currently no effective therapeutic agent for PSC due to a limited understanding of its pathogenesis. Exosomes carry a variety of bioactive molecules, including miRNAs, lipids and proteins. The dysregulation of miRNAs is closely related to various pathological aspects of PSC. In this study, we aimed to investigate the expression profiles and function of exosomal miRNAs under cholestatic conditions.
Methods: One-year-old Mdr2^-/- mice (FVB background) and age-/gender- matched wild-type (WT) mice were used in this study. The serum exosomes were isolated using qEV isolation system (from IZON). Total RNA was isolated from the serum exosomes and liver tissues using Trizol. Total RNA transcriptome was determined by RNAseq. The miRNA profiles in the exosomes and livers were analyzed using the NanoString nCounter® miRNA Expression panel. TargetScan and mirnet2.0 were used to identify the potential target genes. Protein-protein interaction (PPI) was constructed based on the STRING database. Hub genes under cholestatic conditions were visualized by Cytoscape software. The human CCA patient data were downloaded from the Cancer Genome Atlas (TCGA) database.
Results: Bioinformatic analysis identified 74 differentially expressed miRNAs (DE-miRs) and 137 DE-miRs in the male and female Mdr2^-/- mice compared to corresponding WT controls, respectively. The identified DE-miRs were used to predict the target genes. GO and KEGG analysis further showed that these DE-miRs were linked to multiple signaling pathways related to cholestatic liver injury. We also identified 431 up-, and 1106 down- DE genes (DEGs) in the Mdr2^-/- mice compared to WT mice. From the PPI network analysis of the DEGs, the top 100 hug genes were determined. We identified 11 upregulated DE-miRs in Mdr2^-/- mice, which overlapped with the predicted miRs based on the the top 100 hug genes. Interestingly, among the overlapping DE-miRs, two miRs, miR-16 and let-7e, were also upregulated in the serum exosomes of Mdr2^-/- mice. Among the 100 hug genes, we identified nine downregulated genes that were potential targets of miR-16 and let-7e. Further analysis of TCGA CCA data showed upregulation of miR-16 and let-7e, and 7 out of the 9 identified targeted genes were also significantly downregulated in CCA compared to healthy controls.
Conclusion: This study identified potential DE-miRs in the liver and serum exosomes which may contribute to cholestatic liver injury. This study supports the possible use of exosomal miRs as diagnostic and prognostic markers and therapeutic targets for cholestatic liver diseases.

<b>Fig 1. Liver DE-miRs and DEGs identification in Mdr2<sup>-/-</sup> mice compared to corresponding WT controls, target prediction and analysis of these DE-miRs and DEGs, and the overlapping miRNAs in cholestatic liver injury.</b><br /> (A)<b> </b>Liver DE-miRs expression profiles between male Mdr2<sup>-/-</sup> and WT controls. (B) Liver DE-miRs expression profiles between female Mdr2<sup>-/- </sup>and WT controls. (C)<b> </b>Liver DEGs expression profiles between Mdr2<sup>-/-</sup> and WT controls. (D) The top 50 hug genes of the PPI network of upregulated DEGs and the top 50 hug genes of the PPI network of downregulated DEGs. (E) Predicted network of miRNAs – 100 identified hug genes of DEGs between Mdr2<sup>-/- </sup>and WT controls. (F) Identification of 11 overlapping upregulated miRNAs in cholestatic liver injury. FC, fold change; miR, micro RNA. DE, differentially expressed, Mdr2<sup>-/-</sup> , multidrug resistance 2 gene knockout; WT, wild type.

Fig 1. Liver DE-miRs and DEGs identification in Mdr2^-/- mice compared to corresponding WT controls, target prediction and analysis of these DE-miRs and DEGs, and the overlapping miRNAs in cholestatic liver injury.
(A) Liver DE-miRs expression profiles between male Mdr2^-/- and WT controls. (B) Liver DE-miRs expression profiles between female Mdr2^-/-and WT controls. (C) Liver DEGs expression profiles between Mdr2^-/- and WT controls. (D) The top 50 hug genes of the PPI network of upregulated DEGs and the top 50 hug genes of the PPI network of downregulated DEGs. (E) Predicted network of miRNAs – 100 identified hug genes of DEGs between Mdr2^-/-and WT controls. (F) Identification of 11 overlapping upregulated miRNAs in cholestatic liver injury. FC, fold change; miR, micro RNA. DE, differentially expressed, Mdr2^-/- , multidrug resistance 2 gene knockout; WT, wild type.

<b>Fig 2. Exosome DE-miRs identification in Mdr2<sup>-/- </sup>mice compared to corresponding WT controls, the overlapping DE-miRs in the liver and serum exosomes, and comparison of the expression levels of these DE-miRs and target genes between normal and CHOL samples in TCGA CCA database.</b><br /> (A) Exosome DE-miRs expression profiles between female Mdr2<sup>-/-</sup> and WT controls. (B) Identification of 2 overlapping upregulated DE-miRs in the serum exosomes which may contribute to cholestatic liver injury. (C) The expression levels of these 2 overlapping upregulated DE-miRs between normal and CHOL samples. (D) Heat map, the expression levels of these target hug genes of 2 overlapping upregulated DE-miRs between normal and CHOL samples.

Fig 2. Exosome DE-miRs identification in Mdr2^-/-mice compared to corresponding WT controls, the overlapping DE-miRs in the liver and serum exosomes, and comparison of the expression levels of these DE-miRs and target genes between normal and CHOL samples in TCGA CCA database.
(A) Exosome DE-miRs expression profiles between female Mdr2^-/- and WT controls. (B) Identification of 2 overlapping upregulated DE-miRs in the serum exosomes which may contribute to cholestatic liver injury. (C) The expression levels of these 2 overlapping upregulated DE-miRs between normal and CHOL samples. (D) Heat map, the expression levels of these target hug genes of 2 overlapping upregulated DE-miRs between normal and CHOL samples.

Sayed Obaidullah Aseem

Appearances