Establishment and validation of CCL 4 -induced ALI and CLI in mice.
ALI and recovery model, as well as CLI and the recovery model of mice induced by CCL4, were established and validated (Fig. 1).
In acutely injured mice (2D), serum ALT and AST levels increased, hepatocyte necrosis and inflammatory cell infiltration can be observed around the lobular central vein. After five days of recovery (7D), the elevated ALT and AST levels returned to the baseline of control for ALI (AC) group, and histological changes were also recovered (Fig. 1A-B). For CLI, repeated CCL4 treatment induced a slight elevation of serum ALT (8W) and then recovered to baseline of control for CLI (CC) 4 weeks after the cease of CCL4 treatment (12W) (Fig. 1C). Damaged hepatocytes and centrilobular contracture can be found in the liver of CLI (8W), with mild inflammatory cell infiltration. Sirius red staining showed obviously collagen deposition and formation of pseudolobules in 8W CLI mice. These morphological changes were alleviated in recovery mice (12W) (Fig. 1D).
Identification and characterization of serum sEVs from ALI and CLI mice.
The isolated particles were spherical or cup-shaped as observed by transmission electron microscopy (TEM) (Fig. 2A); the exosomal protein markers, including CD63, CD81, and CD9, were all highly expressed (Fig. 2B) as determined by Western blots.
Nanoparticles tracking analysis (NTA) showed that the mean diameter of the particles ranged from 90.2 to 127.8 nm. Compared with control group; the number of particles increased in 2D group, but decreased in 7D group; the particle diameters were smaller in 2D group, but larger in 7D group. We also noticed that the size distribution of sEVs widened with multiple peaks in the 7D group (Fig. 2C). For CLI mice, compared with the control group, the particle concentrations of 8W and 12W groups showed a trend of reduction, although there was no statistical difference. The size distribution of sEVs expanded with multiple peaks in both 8W and 12W groups, and the particle diameters were larger in both 8W and 12W groups (Fig. 2D). The TEM examination revealed that the multiple peaks reflected the aggregation or fusion of mice serum sEVs presented in 7D, 8W, and 12W samples (Fig. 2C-D).
These findings suggested that both ALI and CLI changed the number and morphology of mice serum sEVs, and even when the visible histological changes of liver recovered in 7D and 12W, the changes in particle number and morphology of serum sEVs persisted.
Comparison of the sRNA components in serum sEVs from ALI and CLI mice.
A dramatic increase in both total RNA and sRNA contents was observed in serum sEVs upon ALI (2D (Fig. 3A-B). For CLI, there was no significant difference in either sEVs RNA or small RNA contents among each group. RNA sequencing (RNA-seq) and annotation revealed that each pair of biological repeats had at least 94.97% common sequences in clean reads (Fig S1). In the control groups (AC and CC), tRNA was the dominant sRNA species in serum sEVs, followed by rRNA and miRNA (Fig. 3C). A most remarkable change of serum sEVs sRNAs was the increase of miRNA proportion in both ALI and CLI mice. Compared to the control groups, the proportion of miRNAs increased to more than four folds higher upon ALI, and then returned to baseline in 7D; the proportion of miRNAs increased to almost three folds in CLI mice, and later partially recovered in 12W group. With the increase of miRNA, the proportion of tRNA decreased (Fig. 3D).
The miRNA expression profiles of serum sEVs from ALI and CLI mice.
For ALI and CLI groups, 467 and 488 detectable sEVs (transcript per million reads, TPM ≥ 5.0) miRNAs were obtained, respectively. The biological replications were highly correlated in each group (Fig. 4A). RNA-seq data were further validated by qRT-PCR. Differentially expressed miRNAs of different abundance were selected for validation (Fig. 4B).
Compared to the AC group, 91 miRNAs were up-regulated, and 85 miRNAs were down-regulated (Fold Change ≥ 2.0, P < 0.05) in the 2D group (Fig. 4C, Table S1). Most of these miRNAs recovered to the baseline levels in the 7D group, only eight up-regulated miRNAs and 11 down-regulated miRNAs were detected (Fold Change ≥ 2.0, P < 0.05) (Fig. 4C, Table S2). The cumulative distribution frequency (CDF) plot showed that the top five up-regulated miRNAs (miR-148a-3p, miR-122-5p, miR-192-5p, miR-22-3p, and miR-21a-5p) in 2D group accounted for up to 84.27% of all detectable miRNAs (Fig. 4C). Compared with the CC group, only 13 miRNAs were up-regulated, and six miRNAs were down-regulated (fold change ≥ 2.0, P < 0.05) in the 8W group (Fig. 4D, Table S3); 8 miRNAs were up-regulated, and three miRNAs were down-regulated (fold change ≥ 2.0, P < 0.05) in 12W group (Fig. 4D, Table S4). The CDF analysis showed that the top three up-regulated miRNAs (miR-122-5p, miR-192-5p, and miR-22-3p) in 8W group constituted up to 43.48% of all detectable miRNAs (Fig S2). These findings suggested that both ALI and CLI induced changes in serum sEV miRNAs composition. The changes were caused by the differential expression of a small number of miRNAs with high abundance.
The biological significance of these differentially expressed serum sEV miRNAs in liver injury was explored by Ingenuity Pathway Analysis (IPA, www.qiagen.com/ingenuity). As for the 176 differentially expressed serum sEV miRNAs upon ALI, hepatic steatosis was the most significant hepatotoxicity function, followed by liver inflammation. As for the 19 differentially expressed serum sEV miRNAs upon CLI, hepatocellular carcinoma was the most significant, followed by liver hyperplasia (Fig. 4E).
The liver is the main contributor to the differentially expressed serum sEV miRNAs during AIL and CLI.
We were interested in know the contribution of liver cells to the change of serum sEV miRNAs upon ALI and CLI. Firstly, the liver miRNA expression profile for wild-type male C57BL/6 mice was established using BRB-Array Tool v4.6.0 (https://brb.nci.nih.gov) based on the RNA-seq data from GSE78792 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE78792) [17] (Table S5). Differentially expressed serum sEV miRNAs from ALI and CLI mice were compared to the liver miRNA expression profile. Among the top ten most abundant liver miRNAs, seven were the top ten increased miRNAs in serum sEVs upon ALI (Table 1) and constituted up to 84.70% of increased serum sEV miRNAs; as for CLI, three were raised in serum sEVs. These three miRNAs were exactly the top three miRNAs that increased upon CLI (Table 1) and constituted up to 60.56% of total increased serum sEV miRNAs. These findings suggested the liver as the primary contributor for the up-regulated serum sEV miRNAs during ALI and CLI and convinced that serum sEV miRNA test could be a reliable and sensitive way to monitor either acute or chronic liver injury.
Table 1
Liver enriched miRNAs ranked high in serum sEV miRNAs from AIL and CIL model mice
GSE78792 Liver | Rank | ALI | ALI-Rank | | GSE78792 Liver | Rank | ALI | ALI-Rank |
miR-192 | 1 | miR-192-5p | 3 | | miR-148a | 4 | miR-148a-3p | 1 |
miR-22 | 2 | miR-22-3p | 4 | | miR-122 | 7 | miR-122-5p | 2 |
miR-30a | 3 | miR-30a-5p | 7 | | miR-192 | 1 | miR-192-5p | 3 |
miR-148a | 4 | miR-148a-3p | 1 | | miR-22 | 2 | miR-22-3p | 4 |
miR-21a | 5 | miR-21a-5p | 5 | | miR-21a | 5 | miR-21a-5p | 5 |
miR-26a-2 | 6 | miR-26a-5p | 105 | | miR-30a | 3 | miR-30a-5p | 7 |
miR-122 | 7 | miR-122-5p | 2 | | miR-27b | 10 | miR-27b-3p | 10 |
miR-10a | 8 | n | n | | miR-143 | 9 | miR-143-3p | 92 |
miR-143 | 9 | miR-143-3p | 92 | | miR-26a-2 | 6 | miR-26a-5p | 105 |
miR-27b | 10 | miR-27b-3p | 10 | | miR-10a | 8 | n | n |
GSE78792 Liver | Rank | CLI | CLI-Rank | | GSE78792 Liver | Rank | CLI | CLI-Rank |
miR-192 | 1 | miR-192-5p | 2 | | miR-122 | 7 | miR-122-5p | 1 |
miR-22 | 2 | miR-22-3p | 3 | | miR-192 | 1 | miR-192-5p | 2 |
miR-30a | 3 | n | n | | miR-22 | 2 | miR-22-3p | 3 |
miR-148a | 4 | n | n | | miR-148a | 4 | n | n |
miR-21a | 5 | n | n | | miR-21a | 5 | n | n |
miR-26a-2 | 6 | n | n | | miR-30a | 3 | n | n |
miR-122 | 7 | miR-122-5p | 1 | | miR-26a-2 | 6 | n | n |
miR-10a | 8 | n | n | | miR-10a | 8 | n | n |
miR-143 | 9 | n | n | | miR-143 | 9 | n | n |
miR-27b | 10 | n | n | | miR-27b | 10 | n | n |
ALI, acute liver injury; CLI, chronic liver injury; n, not applicable. |
ALI and CLI signatures based on serum sEV miRNA profiling
To identify serum sEV miRNA signatures for liver injury, we compared the differentially expressed serum sEV miRNAs in various stages of ALI and CLI. Compared to vehicle control, eight miRNAs were up-regulated, and two miRNAs downregulated significantly during the acute injury stage and sustained through the chronic phase (Fold change ≥ 2.0, P < 0.05), which can serve as common liver injury signatures. Besides, 166 miRNAs changed significantly during acute injury, and nine miRNAs changed dramatically during chronic injury. These miRNAs with high abundance (TPM_Mean > 1000), have the potential to be acute or chronic liver injury signatures. The complete lists of these potential ALI and CLI serum sEV miRNA signatures were provided in Table S6. According to their abundance, up to the top 20 miRNAs were listed in Table 2. For the recovery stage, 19 miRNAs changed significantly in the ALI group, 11 miRNAs changed significantly in the CLI group. Some of these miRNAs overlapped with differentially expressed serum sEV miRNAs in corresponding acute or chronic injury stages (Table 2).
Table 2
Common and specific miRNAs for acute and chronic liver injury and recovery in mice model
ALI-Up Specific | CLI-Up Specific | Common-Up | ALI-Down Specific | CLI-Down Specific | Common-Down |
miR-148a-3p | miR-409-3p | miR-122-5p | miR-143-3p | miR-6239 | miR-5112 |
miR-21a-5p | miR-139-3p | miR-192-5p | miR-451a | miR-8109 | miR-9-5p |
miR-92a-3p | miR-690 | miR-22-3p | miR-206-3p | miR-10b-3p | |
miR-30a-5p | miR-211-5p | miR-22-5p | let-7c-5p | miR-8117 | |
miR-140-3p | miR-1941-5p | miR-708-3p | miR-128-3p | | |
miR-29a-3p | | miR-193a-5p | let-7b-5p | | |
miR-27b-3p | | miR-344d-3p | miR-215-5p | | |
miR-101b-3p | | miR-193a-3p | miR-27a-3p | | |
miR-423-3p | | | miR-184-3p | | |
miR-194-5p | | | miR-423-5p | | |
miR-30a-3p | | | miR-126a-3p | | |
miR-1839-5p | | | miR-5106 | | |
miR-152-3p | | | miR-26a-5p | | |
let-7d-3p | | | miR-142a-5p | | |
miR-7a-5p | | | miR-181a-5p | | |
miR-221-3p | | | miR-150-5p | | |
miR-203-3p | | | miR-200a-3p | | |
miR-185-5p | | | miR-100-5p | | |
miR-3535 | | | miR-375-3p | | |
miR-28a-3p | | | miR-139-5p | | |
ALI-R-Up Specific | CLI-R-Up Specific | Common-Up | ALI-R-Down Specific | CLI-R-Down Specific | Common Down |
miR-146b-5p | miR-122-5p | | miR-292a-5p | miR-6239 | |
miR-223-5p | miR-26b-5p | | miR-6538 | miR-8109 | |
miR-27a-5p | miR-142a-3p | | miR-374b-5p | miR-7689-3p | |
miR-147-3p | miR-103-3p | | miR-503-5p | | |
miR-615-3p | miR-374b-5p | | miR-7b-5p | | |
miR-1943-5p | miR-211-5p | | miR-450b-5p | | |
miR-181a-2-3p | miR-337-5p | | miR-431-3p | | |
| miR-1941-5p | | miR-431-5p | | |
| | | miR-292a-3p | | |
| | | miR-450b-3p | | |
| | | miR-181c-3p | | |
ALI, acute liver injury stage, CLI, chronic liver injury stage, ALI-R, acute liver injury-recovery stage, CLI-R, chronic liver injury-recovery stage, Fold change ≥ 2.0, P < 0.05; the bold font, TPM_Mean > 1000; the blue font, overlapping miRNAs in injury and recovery stages. |
Comparison of the potential ALI and CLI serum sEV miRNA signatures to human liver diseases (HLD)-related circulating miRNAs
To explore the biological significance of the ALI and CLI serum sEV miRNA signatures in HLDs, we performed a systemic review of abnormally expressed circulating miRNAs reported in various HLDs. A total of 299 and 257 studies were identified from Pubmed (https://pubmed.ncbi.nlm.nih.gov ) and WOS (http://apps.webofknowledge.com/) database, respectively (Fig. 5A). Data were retrieved from 14 studies, including drug-induced liver injury (DILI) [18], chronic hepatitis B (CHB) [19–23], chronic hepatitis C (CHC) [19, 23–25], non-alcoholic fatty liver disease (NAFLD) [26], non-alcoholic steatohepatitis (NASH) [23], liver cirrhosis (LC) [22, 27, 28] and HCC [19–22, 26–31]. Details of the 14 articles were summarized in Table S7.
Totally 269 abnormally expressed circulating miRNAs related to HLDs were extracted, which appeared ≥ four times were defined as high frequency (Fig. 5A). Of the 30 high-frequency miRNAs (Table S7), 23 were detected in ALI and CLI serum sEVs (Fig. 5B), 14 miRNAs overlapped with ALI and CLI serum sEV signatures, and 12 miRNAs showed the same expression trend. Of the 12 miRNAs, three miRNAs (miR-122-5p, miR-192-5p, and miR-22-3p) were identified as common injury signatures that increased in both ALI and CLI mice and thus have the potential to serve as common signatures for either acute or chronic liver injury. The other 9 miRNAs were identified as ALI signatures (miR-21a-5p, miR-92a-3p, miR-194-5p, miR-17-5p and miR-19b-3p increased, and miR-451a, miR-27a-3p, miR-26a-5p, and miR-223-3p decreased), and may reflect acute or active liver injury (Fig. 5C).
Serum sEVs from liver injury mice induced hepatic macrophages reprogramming
Primary mouse hepatic macrophages were isolated, purified and incubated with SYTO labeled serum sEVs. After 24 h, green fluorescence can be observed in most macrophages in both control and ALI (2D) serum sEVs incubation groups (Fig. 6A). These observations proved that serum sEVs could be uptaken by hepatic macrophages. Serum sEVs from ALI mice (2D) could accelerate the adhesion of hepatic macrophages (Fig. 6B). As determined by qRT-PCR, normal serum sEVs increased the M1-like gene (IL-1B and TNFA) expression, and decrease some of the M2-like gene (IL10RB, CD163, ARG1, and CD206) expression; ALI serum sEVs showed similar effects. However, compared to normal serum sEVs, ALI serum sEVs tended to decrease M1-like IL-1B but increase the M2-like CD163 gene expression (Fig. 6C).
To further dissect the effects of ALI serum sEVs on monocyte-derived and resident hepatic macrophage subgroups, multiple-color flow cytometric analyses were performed to access the expression of M1-like CD86 and the M2-like CD206 in CD11b+F4/80Low monocyte-derived and CD11b+F4/80High resident macrophages ([32] and our unpublished data). It turned out that ALI serum sEVs decreased both CD86 and CD206 expression in CD11b+F4/80Low subgroup, but increased CD206 expression in CD11b+F4/80High subgroup (Fig. 6D-E). These findings indicated that ALI serum sEVs might induce depolarization of CD11b+F4/80Low monocyte-derived macrophages, but M2 differentiation of CD11b+F4/80High resident macrophages.