Single-cell and spatial transcriptomics reveal the fibrosis-related immune landscape of biliary atresia

Spatial genomics: mapping human steatotic liver disease

Metabolic dysfunction-associated steatotic liver disease (MASLD, formerly known as non-alcoholic fatty liver disease) is a leading cause of chronic liver disease worldwide. MASLD can progress to metabolic dysfunction-associated steatohepatitis (MASH, formerly known as non-alcoholic steatohepatitis) with subsequent liver cirrhosis and hepatocellular carcinoma formation. The advent of current technologies such as single-cell and single-nuclei RNA sequencing have transformed our understanding of the liver in homeostasis and disease. The next frontier is contextualizing this single-cell information in its native spatial orientation. This understanding will markedly accelerate discovery science in hepatology, resulting in a further step-change in our knowledge of liver biology and pathobiology. In this Review, we discuss up-to-date knowledge of MASLD development and progression and how the burgeoning field of spatial genomics is driving exciting new developments in our understanding of human liver disease pathogenesis and therapeutic target identification.

Single-cell sequencing combined with spatial transcriptomics reveals that the IRF7 gene in M1 macrophages inhibits the occurrence of pancreatic cancer by regulating lipid metabolism-related mechanisms

AIM

The main focus of this study is to explore the molecular mechanism of IRF7 regulation on RPS18 transcription in M1-type macrophages in pancreatic adenocarcinoma (PAAD) tissue, as well as the transfer of RPS18 by IRF7 via exosomes to PAAD cells and the regulation of ILF3 expression.

METHODS

By utilising single-cell RNA sequencing (scRNA-seq) data and spatial transcriptomics (ST) data from the Gene Expression Omnibus database, we identified distinct cell types with significant expression differences in PAAD tissue. Among these cell types, we identified those closely associated with lipid metabolism. The differentially expressed genes within these cell types were analysed, and target genes relevant to prognosis were identified. Flow cytometry was employed to assess the expression levels of target genes in M1 and M2 macrophages. Cell lines with target gene knockout were constructed using CRISPR/Cas9 editing technology, and cell lines with target gene knockdown and overexpression were established using lentiviral vectors. Additionally, a co-culture model of exosomes derived from M1 macrophages with PAAD cells was developed. The impact of M1 macrophage-derived exosomes on the lipid metabolism of PAAD cells in the model was evaluated through metabolomics analysis. The effects of M1 macrophage-derived exosomes on the viability, proliferation, division, migration and apoptosis of PAAD cells were assessed using MTT assay, flow cytometry, EdU assay, wound healing assay, Transwell assay and TUNEL staining. Furthermore, a mouse PAAD orthotopic implantation model was established, and bioluminescence imaging was utilised to assess the influence of M1 macrophage-derived exosomes on the intratumoural formation capacity of PAAD cells, as well as measuring tumour weight and volume. The expression of proliferation-associated proteins in tumour tissues was examined using immunohistochemistry.

RESULTS

Through combined analysis of scRNA-seq and ST technologies, we discovered a close association between M1 macrophages in PAAD samples and lipid metabolism signals, as well as a negative correlation between M1 macrophages and cancer cells. The construction of a prognostic risk score model identified RPS18 and IRF7 as two prognostically relevant genes in M1 macrophages, exhibiting negative and positive correlations, respectively. Mechanistically, it was found that IRF7 in M1 macrophages can inhibit the transcription of RPS18, reducing the transfer of RPS18 to PAAD cells via exosomes, consequently affecting the expression of ILF3 in PAAD cells. IRF7/RPS18 in M1 macrophages can also suppress lipid metabolism, cell viability, proliferation, migration, invasion and intratumoural formation capacity of PAAD cells, while promoting cell apoptosis.

CONCLUSION

Overexpression of IRF7 in M1 macrophages may inhibit RPS18 transcription, reduce the transfer of RPS18 from M1 macrophage-derived exosomes to PAAD cells, thereby suppressing ILF3 expression in PAAD cells, inhibiting the lipid metabolism pathway, and curtailing the viability, proliferation, migration, invasion of PAAD cells, as well as enhancing cell apoptosis, ultimately inhibiting tumour formation in PAAD cells in vivo. Targeting IRF7/RPS18 in M1 macrophages could represent a promising immunotherapeutic approach for PAAD in the future.

Integrated multi-omics profiling landscape of organising pneumonia

BACKGROUND

Organising pneumonia (OP) is one of the most common and lethal diseases in the category of interstitial pneumonia, along with lung cancer. Reprogramming of lipid metabolism is a newly recognized hallmark of many diseases including cancer, cardiovascular disorders, as well as liver fibrosis and sclerosis. Increased levels of ceramides composed of sphingosine and fatty acid, are implicated in the development of both acute and chronic lung diseases. However, their pathophysiological significance in OP is unclear. The aim of this study was to investigate the role of lipid metabolism reprogramming in OP, focusing on inflammation and fibrosis.

METHODS

Comprehensive multi-omics profiling approaches, including single-cell RNA sequencing, Visium CytAssist spatial transcriptomics, proteomics, metabolomics and mass spectrometry, were employed to analyze the tissues. OP mice model was utilized and molecular mechanisms were investigated in macrophages.

RESULTS

The results revealed a significant association between OP and lipid metabolism reprogramming, characterized by an abnormal expression of several genes related to lipid metabolism, including CD36, SCD1, and CES1 mainly in macrophages. CD36 deficiency in alveolar macrophages, led to an increased expression of C16/24 ceramides that accumulated in mitochondria, resulting in mitophagy or mitochondrial dysfunction. The number of alveolar macrophages in OP was significantly reduced, which was probably due to the ferroptosis signaling pathway involving GSH/SLC3A2/GPX4 through CD36 downregulation in OP. Furthermore, macrophage secretion of DPP7 and FABP4 influenced epithelial cell fibrosis.

CONCLUSIONS

CD36 inhibited the ferroptosis pathway involving SLC3A2/GPX4 in alveolar macrophages of OP tissue by regulating lipid metabolism, thus representing a new anti-ferroptosis and anti-fibrosis effect of CD36 mediated, at least in part, by ceramides.

HIGHLIGHTS

Our findings reveal a significant association between organising pneumonia and lipid metabolism reprogramming and will make a substantial contribution to the understanding of the mechanism of organising pneumonia in patients.

Biliary atresia

Biliary atresia (BA) is a progressive inflammatory fibrosclerosing disease of the biliary system and a major cause of neonatal cholestasis. It affects 1:5,000-20,000 live births, with the highest incidence in Asia. The pathogenesis is still unknown, but emerging research suggests a role for ciliary dysfunction, redox stress and hypoxia. The study of the underlying mechanisms can be conceptualized along the likely prenatal timing of an initial insult and the distinction between the injury and prenatal and postnatal responses to injury. Although still speculative, these emerging concepts, new diagnostic tools and early diagnosis might enable neoadjuvant therapy (possibly aimed at oxidative stress) before a Kasai portoenterostomy (KPE). This is particularly important, as timely KPE restores bile flow in only 50-75% of patients of whom many subsequently develop cholangitis, portal hypertension and progressive fibrosis; 60-75% of patients require liver transplantation by the age of 18 years. Early diagnosis, multidisciplinary management, centralization of surgery and optimized interventions for complications after KPE lead to better survival. Postoperative corticosteroid use has shown benefits, whereas the role of other adjuvant therapies remains to be evaluated. Continued research to better understand disease mechanisms is necessary to develop innovative treatments, including adjuvant therapies targeting the immune response, regenerative medicine approaches and new clinical tests to improve patient outcomes.

Plac8-ERK pathway modulation of monocyte function in sepsis

Sepsis, a life-threatening condition caused by infection, is characterized by the dysregulation of immune responses and activation of monocytes. Plac8, a protein, has been implicated in various inflammatory conditions. This study aimed to investigate the effect of Plac8 upregulation on monocyte proliferation and activation in sepsis patients. Peripheral blood samples were collected from healthy individuals and sepsis patients. Monocytes were stimulated with lipopolysaccharide (LPS) to create an in vitro sepsis model, while a murine sepsis model was established using cecal ligation and puncture (CLP). The levels of monocyte markers, proliferation index (PI), and pro-inflammatory cytokines were assessed using flow cytometry and qPCR, respectively. Plac8 and phosphorylated ERK protein levels were determined by western blot, and TNF-α, IL-6, and IL-10 levels were quantified using ELISA. The CCK-8 assay was used to evaluate PBMC proliferation and activation. The results showed that Plac8 was highly expressed in sepsis models, promoting the survival, proliferation, and activation of monocytes. Plac8 upregulation activated the ERK pathway, leading to increased phosphorylation of ERK protein and elevated levels of CD14, CD16, TNF-α, IL-6, Plac8, and IL-10. In sepsis mice, Plac8 overexpression similarly activated the ERK pathway and promoted the survival, proliferation, and activation of monocytes. In conclusion, the upregulation of Plac8 enhances the activation of the ERK pathway and promotes monocyte proliferation and activation in sepsis patients.

Pkd1l1-deficiency drives biliary atresia through ciliary dysfunction in biliary epithelial cells

BACKGROUND & AIMS

Syndromic biliary atresia is a cholangiopathy characterized by fibro-obliterative changes in the extrahepatic bile duct (EHBD) and congenital malformations including laterality defects. The etiology remains elusive and faithful animal models are lacking. Genetic syndromes provide important clues regarding the pathogenic mechanisms underlying the disease. We investigated the role of the gene Pkd1l1 in the pathophysiology of syndromic biliary atresia.

METHODS

Constitutive and conditional Pkd1l1 knockout mice were generated to explore genetic pathology as a cause of syndromic biliary atresia. We investigated congenital malformations, EHBD and liver pathology, EHBD gene expression, and biliary epithelial cell turnover. Biliary drainage was functionally assessed with cholangiography. Histology and serum chemistries were assessed after DDC (3,5-diethoxycarbony l-1,4-dihydrocollidine) diet treatment and inhibition of the ciliary signaling effector GLI1.

RESULTS

Pkd1l1-deficient mice exhibited congenital anomalies including malrotation and heterotaxy. Pkd1l1-deficient EHBDs were hypertrophic and fibrotic. Pkd1l1-deficient EHBDs were patent but displayed delayed biliary drainage. Pkd1l1-deficient livers exhibited ductular reaction and periportal fibrosis. After DDC treatment, Pkd1l1-deficient mice exhibited EHBD obstruction and advanced liver fibrosis. Pkd1l1-deficient mice had increased expression of fibrosis and extracellular matrix remodeling genes (Tgfα, Cdkn1a, Hb-egf, Fgfr3, Pdgfc, Mmp12, and Mmp15) and decreased expression of genes mediating ciliary signaling (Gli1, Gli2, Ptch1, and Ptch2). Primary cilia were reduced on biliary epithelial cells and altered expression of ciliogenesis genes occurred in Pkd1l1-deficient mice. Small molecule inhibition of the ciliary signaling effector GLI1 with Gant61 recapitulated Pkd1l1-deficiency.

CONCLUSIONS

Pkd1l1 loss causes both laterality defects and fibro-proliferative EHBD transformation through disrupted ciliary signaling, phenocopying syndromic biliary atresia. Pkd1l1-deficient mice function as an authentic genetic model for study of the pathogenesis of biliary atresia.

IMPACT AND IMPLICATIONS

The syndromic form of biliary atresia is characterized by fibro-obliteration of extrahepatic bile ducts and is often accompanied by laterality defects. The etiology is unknown, but Pkd1l1 was identified as a potential genetic candidate for syndromic biliary atresia. We found that loss of the ciliary gene Pkd1l1 contributes to hepatobiliary pathology in biliary atresia, exhibited by bile duct hypertrophy, reduced biliary drainage, and liver fibrosis in Pkd1l1-deficient mice. Pkd1l1-deficient mice serve as a genetic model of biliary atresia and reveal ciliopathy as an etiology of biliary atresia. This model will help scientists uncover new therapeutic approaches for patients with biliary atresia, while pediatric hepatologists should validate the diagnostic utility of PKD1L1 variants.

Foundation models in molecular biology

Determining correlations between molecules at various levels is an important topic in molecular biology. Large language models have demonstrated a remarkable ability to capture correlations from large amounts of data in the field of natural language processing as well as image generation, and correlations captured from data using large language models can also be applicable to solving a wide range of specific tasks, hence large language models are also referred to as foundation models. The massive amount of data that exists in the field of molecular biology provides an excellent basis for the development of foundation models, and the recent emergence of foundation models in the field of molecular biology has really pushed the entire field forward. We summarize the foundation models developed based on RNA sequence data, DNA sequence data, protein sequence data, single-cell transcriptome data, and spatial transcriptome data respectively, and further discuss the research directions for the development of foundation models in molecular biology.

MIF as a potential diagnostic and prognostic biomarker for triple-negative breast cancer that correlates with the polarization of M2 macrophages

Macrophage migration inhibitory factor (MIF) is a proinflammatory cytokine that plays a crucial role in antitumor immunity. However, the role of MIF in influencing the tumor microenvironment (TME) and prognosis of triple-negative breast cancer (TNBC) remains to be elucidated. Using R, we analyzed single-cell RNA sequencing (scRNA-seq) data of 41 567 cells from 10 TNBC tumor samples and spatial transcriptomic data from two patients. Relationships between MIF expression and immune cell infiltration, clinicopathological stage, and survival prognosis were determined using samples from The Cancer Genome Atlas (TCGA) and validated in a clinical cohort using immunohistochemistry. Analysis of scRNA-seq data revealed that MIF secreted by epithelial cells in TNBC patients could regulate the polarization of macrophages into the M2 phenotype, which plays a key role in modulating the TME. Spatial transcriptomic data also showed that epithelial cells (tumor cells) and MIF were proximally located. Analysis of TCGA samples confirmed that tumor tissues of patients with high MIF expression were enriched with M2 macrophages and showed a higher T stage. High MIF expression was significantly associated with poor patient prognosis. Immunohistochemical staining showed high MIF expression was associated with younger patients and worse clinicopathological staging. MIF secreted by epithelial cells may represent a potential biomarker for the diagnosis and prognosis of TNBC and may promote TNBC invasion by remodeling the tumor immune microenvironment.

Identification of ligand and receptor interactions in CKD and MASH through the integration of single cell and spatial transcriptomics

BACKGROUND

Chronic Kidney Disease (CKD) and Metabolic dysfunction-associated steatohepatitis (MASH) are metabolic fibroinflammatory diseases. Combining single-cell (scRNAseq) and spatial transcriptomics (ST) could give unprecedented molecular disease understanding at single-cell resolution. A more comprehensive analysis of the cell-specific ligand-receptor (L-R) interactions could provide pivotal information about signaling pathways in CKD and MASH. To achieve this, we created an integrative analysis framework in CKD and MASH from two available human cohorts.

RESULTS

The analytical framework identified L-R pairs involved in cellular crosstalk in CKD and MASH. Interactions between cell types identified using scRNAseq data were validated by checking the spatial co-presence using the ST data and the co-expression of the communicating targets. Multiple L-R protein pairs identified are known key players in CKD and MASH, while others are novel potential targets previously observed only in animal models.

CONCLUSION

Our study highlights the importance of integrating different modalities of transcriptomic data for a better understanding of the molecular mechanisms. The combination of single-cell resolution from scRNAseq data, combined with tissue slide investigations and visualization of cell-cell interactions obtained through ST, paves the way for the identification of future potential therapeutic targets and developing effective therapies.

Emerging Roles of Spatial Transcriptomics in Liver Research

Spatial transcriptomics, leveraging sequencing- and imaging-based techniques, has emerged as a groundbreaking technology for mapping gene expression within the complex architectures of tissues. This approach provides an in-depth understanding of cellular and molecular dynamics across various states of healthy and diseased livers. Through the integration of sophisticated bioinformatics strategies, it enables detailed exploration of cellular heterogeneity, transitions in cell states, and intricate cell-cell interactions with remarkable precision. In liver research, spatial transcriptomics has been particularly revelatory, identifying distinct zonated functions of hepatocytes that are crucial for understanding the metabolic and detoxification processes of the liver. Moreover, this technology has unveiled new insights into the pathogenesis of liver diseases, such as the role of lipid-associated macrophages in steatosis and endothelial cell signals in liver regeneration and repair. In the domain of liver cancer, spatial transcriptomics has proven instrumental in delineating intratumor heterogeneity, identifying supportive microenvironmental niches and revealing the complex interplay between tumor cells and the immune system as well as susceptibility to immune checkpoint inhibitors. In conclusion, spatial transcriptomics represents a significant advance in hepatology, promising to enhance our understanding and treatment of liver diseases.

A spatiotemporal atlas of cholestatic injury and repair in mice

Cholestatic liver injuries, characterized by regional damage around the bile ductular region, lack curative therapies and cause considerable mortality. Here we generated a high-definition spatiotemporal atlas of gene expression during cholestatic injury and repair in mice by integrating spatial enhanced resolution omics sequencing and single-cell transcriptomics. Spatiotemporal analyses revealed a key role of cholangiocyte-driven signaling correlating with the periportal damage-repair response. Cholangiocytes express genes related to recruitment and differentiation of lipid-associated macrophages, which generate feedback signals enhancing ductular reaction. Moreover, cholangiocytes express high TGFβ in association with the conversion of liver progenitor-like cells into cholangiocytes during injury and the dampened proliferation of periportal hepatocytes during recovery. Notably, Atoh8 restricts hepatocyte proliferation during 3,5-diethoxycarbonyl-1,4-dihydro-collidin damage and is quickly downregulated after injury withdrawal, allowing hepatocytes to respond to growth signals. Our findings lay a keystone for in-depth studies of cellular dynamics and molecular mechanisms of cholestatic injuries, which may further develop into therapies for cholangiopathies.

Technical optimization of spatially resolved single-cell transcriptomic datasets to study clinical liver disease

Single cell and spatially resolved 'omic' techniques have enabled deep characterization of clinical pathologies that remain poorly understood, providing unprecedented insights into molecular mechanisms of disease. However, transcriptomic platforms are costly, limiting sample size, which increases the possibility of pre-analytical variables such as tissue processing and storage procedures impacting RNA quality and downstream analyses. Furthermore, spatial transcriptomics have not yet reached single cell resolution, leading to the development of multiple deconvolution methods to predict individual cell types within each transcriptome 'spot' on tissue sections. In this study, we performed spatial transcriptomics and single nucleus RNA sequencing (snRNAseq) on matched specimens from patients with either histologically normal or advanced fibrosis to establish important aspects of tissue handling, data processing, and downstream analyses of biobanked liver samples. We observed that tissue preservation technique impacts transcriptomic data, especially in fibrotic liver. Single cell mapping of the spatial transcriptome using paired snRNAseq data generated a spatially resolved, single cell dataset with 24 unique liver cell phenotypes. We determined that cell-cell interactions predicted using ligand-receptor analysis of snRNAseq data poorly correlated with cellular relationships identified using spatial transcriptomics. Our study provides a framework for generating spatially resolved, single cell datasets to study gene expression and cell-cell interactions in biobanked clinical samples with advanced liver disease.

Long non-coding RNA growth arrest-specific 5 inhibits liver fibrogenesis in biliary atresia by interacting with microRNA-222 and repressing IGF1/AKT signaling

Background

Long non-coding RNA growth arrest-specific 5 (lncRNA GAS5) has been shown to inhibit liver fibrosis through serving as a competing endogenous RNA for microRNA-222 (miR-222). Progressive liver fibrosis is a typical characteristic of biliary atresia (BA). However, the role of GAS5/miR-222 and its underlying mechanisms remain largely unknown in BA.

Methods

The expression of GAS5 was determined in the liver and primary hepatic stellate cells (HSCs) of BA patients. Then, the effects of GAS5 on the activation and proliferation of HSCs were evaluated. Furthermore, the interaction between GAS5 and miR-222 was investigated by a luciferase gene report assay. Next, the effects of IGF1/AKT signaling were determined to clarify the downstream mechanism of GAS5. Finally, GAS5 administration was performed to explore its role in an experimental BA mouse model.

Results

GAS5 expression was decreased in liver tissues and HSCs of BA patients, and was inversely correlated with liver fibrosis in BA. Up-regulation of GAS5 in LX-2 cells significantly reduced smooth muscle α-actin (α-SMA) and collagen 1a1 (COL1A1) expression, inhibited cell proliferation and clone formation ability, induced S phase increase, and promoted cell apoptosis. Moreover, GAS5 was negatively regulated by miR-222, which promoted HSCs activation and proliferation, and was positively correlated with liver fibrosis in BA. Additionally, the expressions of IGF1, p-PI3K, and p-AKT were decreased when LX-2 cells over-expressed GAS5, whereas knockdown of IGF1 or AKT significantly decreased α-SMA and COL1A1 expression, suppressed cell proliferation, and enhanced cell apoptosis in LX-2 cells. Furthermore, GAS5 administration significantly increased apoptosis and reduced liver fibrosis, α-SMA and COL1A1 expressions in liver tissues of BA mice.

Conclusions

GAS5 inhibited liver fibrosis in BA by interacting with miR-222 and regulating IGF1/AKT signaling, which may be a therapeutic target to alleviate liver fibrosis in BA.

Single cell RNA-sequencing analysis reveals that N-acetylcysteine partially reverses hepatic immune dysfunction in biliary atresia

Background & Aims

Our previous study indicated that CD177+ neutrophil activation has a vital role in the pathogenesis of biliary atresia (BA), which is partially ameliorated by N-acetylcysteine (NAC) treatment. Here, we evaluated the clinical efficacy of NAC treatment and profiled liver-resident immune cells via single cell RNA-sequencing (scRNA-seq) analysis to provide a comprehensive immune landscape of NAC-derived immune regulation.

Methods

A pilot clinical study was conducted to evaluate the potential effects of intravenous NAC treatment on infants with BA, and a 3-month follow-up was carried out to assess treatment efficacy. scRNA-seq analysis of liver CD45+ immune cells in the control (n = 4), BA (n = 6), and BA + NAC (n = 6) groups was performed and the effects on innate cells, including neutrophil and monocyte-macrophage subsets, and lymphoid cells were evaluated.

Results

Intravenous NAC treatment demonstrated beneficial efficacy for infants with BA by improving bilirubin metabolism and bile acid flow. Two hepatic neutrophil subsets of innate cells were identified by scRNA-seq analysis. NAC treatment suppressed oxidative phosphorylation and reactive oxygen species production in immature neutrophils, which were transcriptionally and functionally similar to CD177+ neutrophils. We also observed the suppression of hepatic monocyte-mediated inflammation, decreased levels of oxidative phosphorylation, and M1 polarisation in Kupffer-like macrophages by NAC. In lymphoid cells, enhancement of humoral immune responses and attenuation of cellular immune responses were observed after NAC treatment. Moreover, cell-cell interaction analysis showed that innate/adaptive proinflammatory responses were downregulated by NAC.

Conclusions

Our clinical and scRNA-seq data demonstrated that intravenous NAC treatment partially reversed liver immune dysfunction, alleviated the proinflammatory responses in BA by targeting innate cells, and exhibited beneficial clinical efficacy.

Impact and implications

BA is a serious liver disease that affects newborns and has no effective drug treatment. In this study, scRNA-seq showed that NAC treatment can partially reverse the immune dysfunction of neutrophil extracellular trap-releasing CD177+ neutrophils and Kupffer cells, and lower the inflammatory responses of other innate immune cells in BA. In consequence, intravenous NAC treatment improved the clinical outcomes of patients with BA in term of bilirubin metabolism.

Crosstalk between alternative splicing and inflammatory bowel disease: Basic mechanisms, biotechnological progresses and future perspectives

BACKGROUND

Alternative splicing (AS) is an omnipresent regulatory mechanism of gene expression that enables the generation of diverse splice isoforms from a single gene. Recently, AS events have gained considerable momentum in the pathogenesis of inflammatory bowel disease (IBD).

METHODS

Our review has summarized the complex process of RNA splicing, and firstly highlighted the potential involved molecules that target aberrant splicing events in IBD. The quantitative transcriptome analyses such as microarrays, next-generation sequencing (NGS) for AS events in IBD have been also discussed.

RESULTS

Available evidence suggests that some abnormal splicing RNAs can lead to multiple intestinal disorders during the onset of IBD as well as the progression to colitis-associated cancer (CAC), including gut microbiota perturbations, intestinal barrier dysfunctions, innate/adaptive immune dysregulations, pro-fibrosis activation and some other risk factors. Moreover, current data show that the advanced technologies, including microarrays and NGS, have been pioneeringly employed to screen the AS candidates and elucidate the potential regulatory mechanisms of IBD. Besides, other biotechnological progresses such as the applications of third-generation sequencing (TGS), single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST), will be desired with great expectations.

CONCLUSIONS

To our knowledge, the current review is the first one to evaluate the potential regulatory mechanisms of AS events in IBD. The expanding list of aberrantly spliced genes in IBD along with the developed technologies provide us new clues to how IBD develops, and how these important AS events can be explored for future treatment.

Single-cell analysis reveals the immune heterogeneity and interactions in lungs undergoing hepatic ischemia-reperfusion

Hepatic ischemia-reperfusion IR (HIR) is an unavoidable pathophysiological process during liver transplantation, resulting in systematic sterile inflammation and remote organ injury. Acute lung injury (ALI) is a serious complication after liver transplantation with high postoperative morbidity and mortality. However, the underlying mechanism is still unclear. To assess the phenotype and plasticity of various cell types in the lung tissue microenvironment after HIR at the single-cell level, single-cell RNA sequencing (scRNA-seq) was performed using the lungs from HIR-induced mice. In our results, we identified 23 cell types in the lungs after HIR and found that this highly complex ecosystem was formed by subpopulations of bone marrow-derived cells that signaled each other and mediated inflammatory responses in different states and different intervals. We described the unique transcriptional profiles of lung cell clusters and discovered two novel cell subtypes (Tspo+Endothelial cells and Vcan+ monocytes), as well as the endothelial cell-immune cell and immune cell-T cell clusters interactome. In addition, we found that S100 calcium binding protein (S100a8/a9), specifically and highly expressed in immune cell clusters of lung tissues and exhibited detrimental effects. Finally, the cellular landscape of the lung tissues after HIR was established, highlighting the heterogeneity and cellular interactions between major immune cells in HIR-induced lungs. Our findings provided new insights into the mechanisms of HIR-induced ALI and offered potential therapeutic target to prevent ALI after liver transplantation.

Technical optimization of spatially resolved single-cell transcriptomic datasets to study clinical liver disease

Single cell and spatially resolved 'omic' techniques have enabled deep characterization of clinical pathologies that remain poorly understood, providing unprecedented insights into molecular mechanisms of disease. However, transcriptomic platforms are costly, limiting sample size, which increases the possibility of pre-analytical variables such as tissue processing and storage procedures impacting RNA quality and downstream analyses. Furthermore, spatial transcriptomics have not yet reached single cell resolution, leading to the development of multiple deconvolution methods to predict individual cell types within each transcriptome 'spot' on tissue sections. In this study, we performed spatial transcriptomics and single nucleus RNA sequencing (snRNASeq) on matched specimens from patients with either histologically normal or advanced fibrosis to establish important aspects of tissue handling, data processing, and downstream analyses of biobanked liver samples. We observed that tissue preservation technique impacts transcriptomic data, especially in fibrotic liver. Deconvolution of the spatial transcriptome using paired snRNASeq data generated a spatially resolved, single cell dataset with 24 unique liver cell phenotypes. We determined that cell-cell interactions predicted using ligand-receptor analysis of snRNASeq data poorly correlated with celullar relationships identified using spatial transcriptomics. Our study provides a framework for generating spatially resolved, single cell datasets to study gene expression and cell-cell interactions in biobanked clinical samples with advanced liver disease.

Novel cell subtypes of SPP1 + S100P+, MS4A1-SPP1 + S100P+ were key subpopulations in intrahepatic cholangiocarcinoma

BACKGROUND

In this study, we integrated single-cell RNA sequencing (scRNA-seq) data to investigate cell heterogeneity and utilized MSigDB and CIBERSORTx to explore the pathways of major cell types and the relationships between different cell subtypes. Subsequently, we explored the correlation of cell subtypes with survival and used Gene Set Enrichment Analysis (GSEA) analyses to assess the pathways associated with the infiltration of specific cell subtypes. Finally, multiplex immunohistochemistry in tissue microarray cohort were performed to validate differences in protein level and their correlation with survival.

RESULTS

iCCA presented a unique immune ecosystem, with increased proportions of Epi (epithelial)-SPP1-2, Epi-S100P-1, Epi-DN (double negative for SPP1 and S100P expression)-1, Epi-DN-2, Epi-DP (double positive for SPP1 and S100P expression)-1, Plasma B-3, Plasma B-2, B-HSPA1A-1, B-HSPA1A-2 cells, and decreased proportions of B-MS4A1. High level of Epi-DN-2, Epi-SPP1-1, Epi-SPP1-2, B-MS4A1, and low level of Epi-DB-1, Epi-S100P-1, and Epi-S100P-2 was significantly associated with longer overall survival (OS), and high level of B-MS4A1_Low_Epi-DN-2_Low was associated with the shortest OS. Moreover, the results of MsigDB and GSEA suggest that bile acid metabolism is a crucial process in iCCA. Finally, we found that S100P+, SPP1+, SPP1 + S100P+, and MS4A1-SPP1 + S100P+ were highly expressed, whereas MS4A1 was lowly expressed in iCCA, and patients with high level of S100P+, SPP1 + S100P+, and MS4A1-SPP1 + S100P+ exhibited shorter survival.

CONCLUSIONS

We identified the cell heterogeneity of iCCA, found that iCCA is a unique immune ecosystem with many cell subtypes, and showed that the novel cell subtypes of SPP1 + S100P+ and MS4A1-SPP1 + S100P+ were key subpopulations in iCCA.

Single-cell sequencing of a novel model of neonatal bile duct ligation in mice identifies macrophage heterogeneity in obstructive cholestasis

Macrophages (MΦ) play a role in neonatal etiologies of obstructive cholestasis, however, the role for precise MΦ subsets remains poorly defined. We developed a neonatal murine model of bile duct ligation (BDL) to characterize etiology-specific differences in neonatal cholestatic MΦ polarization. Neonatal BDL surgery was performed on female BALB/c mice at 10 days of life (DOL) with sham laparotomy as controls. Comparison was made to the Rhesus Rotavirus (RRV)-induced murine model of biliary atresia (BA). Evaluation of changes at day 7 after surgery (BDL and sham groups) and murine BA (DOL14) included laboratory data, histology (H&E, anti-CD45 and anti-CK19 staining), flow cytometry of MΦ subsets by MHCII and Ly6c expression, and single cell RNA-sequencing (scRNA-seq) analysis. Neonatal BDL achieved a 90% survival rate; mice had elevated bile acids, bilirubin, and alanine aminotransferase (ALT) versus controls (p < 0.05 for all). Histology demonstrated hepatocellular injury, CD45+ portal infiltrate, and CK19+ bile duct proliferation in neonatal BDL. Comparison to murine BA showed increased ALT in neonatal BDL despite no difference in histology Ishak score. Neonatal BDL had significantly lower MHCII-Ly6c+ MΦ versus murine BA, however, scRNA-seq identified greater etiology-specific MΦ heterogeneity with increased endocytosis in neonatal BDL MΦ versus cellular killing in murine BA MΦ. We generated an innovative murine model of neonatal obstructive cholestasis with low mortality. This model enabled comparison to murine BA to define etiology-specific cholestatic MΦ function. Further comparisons to human data may enable development of immune modulatory therapies to improve patient outcomes.

Analysis of intercellular communication in the osteosarcoma microenvironment based on single cell sequencing data

Osteosarcoma (OS) is the most common primary bone cancer in children and young adults, patient survival rates have not improved in recent decades. To further understand the interrelationship between different cell types in the tumor microenvironment of osteosarcoma, we comprehensively analyzed single-cell sequencing data from six patients with untreated osteosarcoma. Nine major cell types were identified from a total of 46,046 cells based on unbiased clustering of gene expression profiles and canonical markers. Osteosarcoma from different patients display heterogeneity in cellular composition. Myeloid cells were the most commonly represented cell type, followed by osteoblastic and TILs. Copy number variation (CNV) results identified amplifications and deletions in malignant osteoblastic cells and fibroblasts. Trajectory analysis based on RNA velocity showed that osteoclasts in the OS microenvironment could be differentiated from myeloid cells. Furthermore, we explored the intercellular communications in OS microenvironment and identified multiple ligand-receptor pairs between myeloid cells, osteoblastic cells and their cells, including 21 ligand-receptor pair genes that significantly associated with survival outcomes. Importantly, we found chemotherapy may have an effect on cellular communication in the OS microenvironment by analyzing single-cell sequencing data from seven primary osteosarcoma patients who received chemotherapy. We believe these observations will improve our understanding of potential mechanisms of microenvironment contributions to OS progression and help identify potential targets for new treatment development in the future.

Macrophage in liver Fibrosis: Identities and mechanisms

Liver fibrosis is a chronic disease characterized by the deposition of extracellular matrix and continuous loss of tissues that perform liver functions. Macrophages are crucial modulators of innate immunity and play important roles in liver fibrogenesis. Macrophages comprise heterogeneous subpopulations that exhibit different cellular functions. Understanding the identity and function of these cells is essential for understanding the mechanisms of liver fibrogenesis. According to different definitions, liver macrophages are divided into M1/M2 macrophages or monocyte-derived macrophages/Kupffer cells. Classic M1/M2 phenotyping corresponds to pro- or anti-inflammatory effects, and, therefore, influences the degree of fibrosis in later phases. In contrast, the origin of the macrophages is closely associated with their replenishment and activation during liver fibrosis. These two classifications of macrophages depict the function and dynamics of liver-infiltrating macrophages. However, neither description properly elucidates the positive or negative role of macrophages in liver fibrosis. Critical tissue cells mediating liver fibrosis include hepatic stellate cells and hepatic fibroblasts, with hepatic stellate cells being of particular interest because of their close association with macrophages in liver fibrosis. However, the molecular biological descriptions of macrophages are inconsistent between mice and humans, warranting further investigations. In liver fibrosis, macrophages can secrete various pro-fibrotic cytokines, such as TGF-β, Galectin-3 and interleukins (ILs), and fibrosis-inhibiting cytokines, such as IL10. These different secretions may be associated with the specific identity and spatiotemporal characteristics of macrophages. Furthermore, during fibrosis dissipation, macrophages may degrade extracellular matrix by secreting matrix metalloproteinases (MMPs). Notably, using macrophages as therapeutic targets in liver fibrosis has been explored. The current therapeutic approaches for liver fibrosis can by categorized as follows: treatment with macrophage-related molecules and macrophage infusion therapy. Although there have been limited studies, macrophages have shown reliable potential for liver fibrosis treatment. In this review, we focu on the identity and function of macrophages and their relationship to the progression and regression of liver fibrosis.