Dear Editor,
As a formidable health-care burden, liver disease affects approximately 844 million people worldwide, and around 2 million patients with liver disease die per year (
Marcellin and Kutala, 2018). For patients suffering from end-stage liver disease, liver transplantation is the most ideal and final treatment (
Ferreira-Gonzalez et al., 2022). Unfortunately, up to 35% of the recipients after liver transplantation are subjected to biliary complications (
Ferreira-Gonzalez et al., 2022), such as anastomotic (AS), non-anastomotic biliary strictures (NAS), and bile leakage, which represent the major causes of morbidity and graft failure after liver transplantation (
Brunner et al., 2013). One of the major reasons is that, compared with other liver parenchymal cells, cholangiocytes are more vulnerable to the effects of static cold storage (SCS) (
Ferreira-Gonzalez et al., 2022), the standard clinical method of organ preservation. SCS results in the deterioration of the biliary tract, such as the loss of biliary epithelia, mural necrosis, and the damage of peribiliary vascular plexus (
Brunner et al., 2013). Although more advanced preservation techniques are being developed, such as normo- or hypothermic
ex situ perfusion technologies, the improvement of biliary complications is limited (
Nasralla et al., 2018).
In nature, many mammals can hibernate to survive harsh conditions with food shortage and low ambient temperature (
Mohr et al., 2020). During hibernation, small mammalian hibernators, such as Syrian hamsters, ground squirrels, chipmunks, bats, and so on, can reduce their body temperatures to below 10°C from ~37°C normothermia and repeat many torpor-arousal cycles (
Carey et al., 2003). Thus, their organs are subjected to repetitive cooling-rewarming stresses, yet can be shielded from such detrimental insults (
Dugbartey et al., 2018). Recent studies have shown that hibernator-derived cells are more tolerant to cold stress than cells derived from nonhibernators (e.g., humans and mice) (
Anegawa et al., 2021;
Ou et al., 2018;
Sone et al., 2024). However, these studies mainly depend on primary cells or immortalized cell lines, which cannot fully reflect the cold adaptation mechanism at the level of
in vivo tissues.
Liver-derived cholangiocyte organoids serve as an ideal model for biliary epithelial cells, as they possess the remarkable capacity of self-organization and can maintain the majority of biliary characteristics during
in vitro culture (
Huch et al., 2015). Several studies have proved the superiority of cholangiocyte organoids in modeling biliary cold storage, ischemia and reperfusion injury, and recapitulating cholangiopathy-associated programmed cell death (
Ferreira-Gonzalez et al., 2022). In this study, we attempted to build intrahepatic cholangiocyte organoids (ICOs) from the mammalian hibernator, Syrian hamster (
Mesocricetus auratus), and performed a comparison between Syrian hamster ICOs (shICOs) and mouse ICOs (mICOs) to seek for molecular clues of biliary cold adaptation.
We first determined whether Syrian hamster cholangiocytes could cope with cold stress better than their non-hibernator (i.e., mouse) counterparts could. Using the University of Wisconsin (UW) solution, the standard organ cold preservation solution in clinical practice, we perfused livers of Syrian hamsters and mice and preserved them under SCS at 4°C for different times (i.e., 1, 3 and 5 days) (Fig. 1A). The hematoxylin and eosin (H&E) staining indicated that mouse livers after SCS exhibited severe cold-induced injuries, with focal necrotic areas of hepatocytes, pronounced cellular edema, substantial disruption of the lobular architecture, and particularly the peeling of periportal cells (Figs. 1B and S1A). Moreover, the immunohistochemical (IHC) staining and immunofluorescence (IF) staining for the cholangiocyte marker Keratin19 (KRT19) revealed that the morphology of mouse bile ducts was progressively altered with the extension of SCS (Fig. S1B and S1C). Mouse bile ducts suffered from impaired duct continuity, cell death, and detachment of cholangiocytes into the lumen, akin to those observed in human NAS. The IF staining for inflammatory markers TNF-α and CD68 further demonstrated that the mouse liver (especially cholangiocytes and hepatocytes) presented a pronounced inflammatory response since the first day of SCS (Figs. 1B and S1C). In contrast, the Syrian hamster liver preserved a well-maintained tissue structure and morphology, and no severe bile duct damage and evident inflammatory response were detected even after 5 days of SCS. These results imply that Syrian hamster cholangiocytes exhibit a better cold resistance ability than their mouse counterparts do during SCS, which also suggests Syrian hamster as an ideal model for studying organ cold adaptation.
Next, we attempted to construct shICOs to characterize the cold adaptation mechanisms of Syrian hamster bile ducts. However, the direct use of mouse culture medium failed to support shICOs (Fig. 1C). To identify the molecular differences between Syrian hamster and mouse cholangiocytes, we used single-nucleus RNA sequencing (snRNA-seq) to profile liver transcriptomes of these two species, each with two replicates. In total, 26,384 single nuclei were obtained, and five major cell groups were identified, namely, hepatocytes, cholangiocytes, stellate cells, endothelial cells, and immune cells (Figs. 1D and S2A–D; Table S1). We then performed differential expression analysis and detected a clear difference between Syrian hamster and mouse cholangiocytes (Fig. S2E; Table S2). Gene ontology analysis demonstrated that the upregulated genes of Syrian hamster cholangiocytes were predominantly enriched in terms associated with neuronal development, cAMP, Wnt, and TGF-β signaling pathways (Figs. 1E and S2E). However, the high expression level of
Smurf1 indicated the inhibition of TGF-β signaling in Syrian hamster cholangiocytes, as
Smurf1 encodes a ubiquitin ligase that can ubiquitinate and target TGFBR1 for degradation (
Yan et al., 2011).
Given the above results, we successfully established the culture system of shICOs by adding a composition combination termed as NFRA (N: N2 supplement, for neuronal development; F: forskolin, a cAMP activator; R: R-spondin-1 conditioned medium, for Wnt activation; A: A8301, a TGF-β inhibitor) to the mouse culture medium, and the removal of any one of the four factors would lead to the failure of shICO culture (Figs. 1F and S2F). Notably, shICOs manifested as spherical structures composed of a single layer of cubical epithelium, which could be maintained for a long-term expansion, such as over 20 passages.
Furthermore, the IHC and IF staining for cholangiocyte markers, such as SOX9, KRT19 and EPCAM, suggested their high expression levels in shICOs (Figs. 1G S3A and S3B). Meanwhile, the IF staining for MKI67 also supported the proliferative feature of shICOs. In addition, we conducted snRNA-seq to dissect the cell composition of shICOs and obtained 11,928 single nuclei (Fig. S3C). The integration result of shICO cells and in vivo liver cells revealed that more than 90% of shICO cells were grouped with in vivo cholangiocytes, while other shICO cells were grouped with in vivo hepatocytes (Fig. 1H and 1I). Notably, shICO cells also exhibited consistent expression patterns of marker genes with in vivo hepatocytes (e.g., Hnf4a and Tfr2) and cholangiocytes (e.g., Epcam and Pkhd1) (Fig. S3D and S3E; Table S3). Taken together, these results demonstrate the success of the culture system of shICOs.
Subsequently, we investigated whether shICOs could recapitulate the cold resistance ability of Syrian hamster cholangiocytes. We utilized shICOs to simulate the SCS and reperfusion processes during liver transplantation by cooling (4°C) and rewarming (37°C) experiments, and mICOs were used as controls (Fig. 1J). shICOs and mICOs were initially cultured at 37°C under normal conditions, then were transferred to a cooling condition at 4°C, and finally were rewarmed back to 37°C under normal conditions. Importantly, remarkable differences between shICOs and mICOs were observed (Fig. 1K). After 5 days of cooling, the vast majority of mICO spheroids underwent shrinkage and turned black, and almost no mICOs survived after 2 days of rewarming. In contrast, shICOs were able to preserve their spheroid morphology under cold stress and promptly resume proliferation with cells adhering to the culture vessel wall upon rewarming.
Furthermore, we performed ATP assay to evaluate the cell viability of shICOs and mICOs (Fig. 1L). Upon cooling, ATP levels gradually decreased in mICOs, and only less than 10% was left in mICOs after 5 days of cooling. On the other hand, shICOs only exhibited a sharp decrease of ATP levels after 5 days of cooling, but still maintained a 50% level, and quickly recovered upon rewarming. In addition, propidium iodide (PI) staining of dead cells indicated that 3 days of cooling at 4°C induced widespread cell death in mICOs but not in shICOs (Figs. 1M, 1N and S4A). Therefore, these results indicate that shICOs present superior resistance ability to cooling-rewarming stress compared with mICOs.
Recent investigations have indicated that the cell death triggered by cold exposure manifests the hallmarks of ferroptosis, which is typified by substantial lipid peroxidation facilitated through iron-ion-mediated generation of reactive oxygen species (ROS) (
Sone and Yamaguchi, 2024). Thus, we first assessed the levels of ROS by the live cell fluorescent probe MitoSOX Red that targets mitochondrial superoxide. A significant increase in ROS production was observed in mICOs but not in shICOs after 3 days of cooling at 4°C (Fig. 2A and 2B). We further measured the levels of malondialdehyde (MDA) to assess lipid peroxidation (Fig. 2C). Notably, within the initial day of cold exposure, intense lipid peroxidation was detected in mICOs, with a more than 2-fold elevation in the MDA concentration. Surprisingly, even after 5 days of cooling, the MDA concentration in shICOs was still maintained at a low level. In addition, we also measured the levels of apoptosis with caspase-3/7 activity in both shICOs and mICOs under cold stress (Fig. 2D). Throughout the cooling process, shICOs exhibited no evident apoptosis, whereas mICOs showed significant apoptosis starting from the third day of cold exposure.
To determine whether the results of cold stress-induced cell death obtained from in vitro organoids could truly reflect the results at the in vivo tissue level, we then evaluated the levels of lipid peroxidation and apoptosis in liver tissues of hamsters and mice following 1-, 3-, and 5-day SCS (Fig. 2E). According to the IF staining results for lipid peroxidation marker 4-Hydroxynonenal (4HNE) and apoptosis marker Cleaved Caspase-3 (CC3), mouse livers experienced lipid peroxidation and apoptosis within 1-day and 3-day SCS, respectively, indicating that lipid peroxidation occurs earlier than apoptosis in mouse livers. In contrast, for Syrian hamster livers, lipid peroxidation occurred at a later stage of SCS, and no apoptotic signs were detected even after 5 days of SCS. Furthermore, we examined the impacts of the ferroptosis inhibitor Ferrostatin-1 (Fer-1) and the apoptosis inhibitor Z-VAD(OH)-FMK on the cell viability of mICOs under cooling-rewarming stress (Figs. 2F and S4B). Unlike the weak protective effect of Z-VAD(OH)-FMK, Fer-1 substantially augmented the cell viability of mICOs under cold stress and rapidly assisted mICOs in recovering their vitality after rewarming. To sum up, although mICOs encountered both ferroptosis and apoptosis under cold stress, ferroptosis occurred much earlier and had a greater impact on the survival rate of mICOs. Compared with mICOs, shICOs were able to maintain low levels of ROS and lipid peroxidation, suggesting a better anti-ferroptosis ability.
To explore the molecular mechanisms underlying the disparity in anti-ferroptosis ability between shICOs and mICOs, their transcriptomes were profiled at 6 time points, namely, 37°C normal condition, 6-hour, 3-day, and 5-day 4°C cooling condition, and 6-hour and 2-day 37°C rewarming condition (Fig. S5A). We investigated the expression patterns of the marker genes that are related to the canonical ferroptosis pathway (Fig. S5B). Importantly, aside from the moderate upregulation of the renowned anti-lipid peroxidation gene,
Gpx4), shICOs also exhibited a much higher expression level of
Fth1 that encodes the heavy subunit of ferritin (Fig. 2G). Ferritin stores iron in a form that is not only soluble and non-toxic but also easily retrievable, which is crucial in maintaining iron homeostasis (
Zhang et al., 2021). Moreover, ferritin consists of 24 subunits of the heavy (FTH) and light (FTL) chains, and the variation in the subunit composition may affect its function. Interestingly, the expression level of
Ftl that encodes the light subunit of ferritin was much lower in shICOs than in mICOs, resulting in a remarkably distinct FTH/FTL ratio of ferritin between shICOs and mICOs (Fig. 2H). As ferritin with more heavy chains has a stronger antioxidant activity (
Zhang et al., 2021), the higher FTH/FTL ratio of ferritin in shICOs might suggest their superior antioxidant ability during the cooling-rewarming process.
Furthermore, we also conducted IF staining for FTH1 in mouse and Syrian hamster liver tissues during SCS (Fig. 2I). Notably, FTH1 exhibited much higher expression levels in Syrian hamster livers than mouse livers throughout the SCS process, which was consistent with the organoid results. In addition, shICOs more highly expressed
Slc40a1 (Fig. 2G), which encodes an ion transport protein that can transport Fe
2+ from the intracellular to the extracellular space and is also crucial to the maintenance of iron homeostasis (
Zhang et al., 2018). Based on the higher expression of
Fth1 and
Slc40a1, we inferred a better maintenance of iron homeostasis in shICOs. Thus, we improved the UW solution by adding the iron chelator, deferoxamine (DFO), which might help to chelate the cellular free iron and reduce the generation of toxic ROS. Compared with the standard UW solution, UW solution with DFO indeed improved the cold preservation of mouse bile ducts, with better duct continuity, less detachment of cholangiocytes into the lumen, and the alleviation of ferroptosis and inflammation (Figs. 2J and S5C). Similarly, DFO treatment also improved the cold preservation of rat bile ducts (Fig. S6).
In conclusion, using the mammalian hibernator, Syrian hamster, as an animal model, we demonstrated the superior cold resistance ability in their cholangiocytes and recapitulated such ability
in vitro by building the culture system of shICOs. Recently, several studies have utilized hibernator-derived cells to study cellular cold tolerance (
Anegawa et al., 2021;
Ou et al., 2018;
Sone et al., 2024); however, these cells are much inferior in the self-organization capacity and the maintenance of
in vivo characteristics compared with organoids. We further demonstrated the stronger anti-ferroptosis ability of shICOs than mICOs under cold stress. Although cold exposure-induced cell death exhibits the characteristics of ferroptosis (
Sone and Yamaguchi, 2024), the results obtained by different systems are not exactly the same (
Anegawa et al., 2021;
Sone et al., 2024). The cross-species comparison between shICOs and mICOs revealed that, in addition to
Gpx4,
Fth1 and
Slc40a1 were also upregulated in shICOs during the cooling-rewarming process, implying the importance of iron homeostasis in biliary cold adaptation. Collectively, our findings provide unique insights into improving biliary cold preservation, and our shICO model will be a useful tool for the future study of mammalian hibernation.
The Author(s) 2025. Published by Oxford University Press on behalf of Higher Education Press.