BMP9 and BMP10 coordinate liver cellular crosstalk to maintain liver health

The liver is the largest solid organ in the body and is primarily composed of HCs, ECs, KCs, and HSCs, which spatially interact and cooperate with each other to maintain liver homeostasis. However, the complexity and molecular mechanisms underlying the crosstalk between these different cell types remain to be revealed. Here, we generated mice with conditional deletion of Bmp9/10 in different liver cell types and demonstrated that HSCs were the major source of BMP9 and BMP10 in the liver. Using transgenic ALK1 (receptor for BMP9/10) reporter mice, we found that ALK1 is expressed on KCs and ECs other than HCs and HSCs, and BMP9/10 secreted by HSCs promotes the differentiation of KCs and ECs and maintain their identity. Pdgfb expression was significantly upregulated in KCs and ECs after BMP9 and BMP10 deletion, ultimately leading to HSCs activation and liver fibrosis. ECs express several angiocrine factors, such as BMP2, BMP6, Wnt2 and Rspo3, to regulate hepatocyte iron metabolism and metabolic zonation. We found that these angiocrine factors were significantly decreased in ECs from Bmp9/10HSC-KO mice, which further resulted in liver iron overload and disruption of HC zonation. In summary, we demonstrated that HSCs play a central role in mediating liver cell‒cell crosstalk via the production of BMP9/10, highlighting the important role of intercellular interaction in organ development and homeostasis.


Introduction
Almost all organs in vertebrates consist of different cell types that undertake their individual functions to maintain organ homeostasis (1), such as vascular endothelial cells supplying oxygen and nutrients to the organ (2) and immune cells, especially tissue-resident macrophages, constantly clearing cell debris and monitoring tissue microenvironments to prevent pathogen infection (3).Usually, these different cell types originate from distinct progenitors and differentiate into specialized cell types via the expression of lineage-specific transcription factors during development (4).Meanwhile, crosstalk between these individual cells within one organ occurs and promotes them to continue to specialize (5,6), which is needed for the maintenance of the organotypic phenotypes of these cells.However, the cell-cell communication within organs is largely unknown.
The liver is composed of four major cell types: hepatocytes (HCs), Kupffer cells (KCs), endothelial cells (ECs), and hepatic stellate cells (HSCs).The liver exemplifies the importance of cellular interactions in maintaining organ integrity (7)(8)(9).For example, secretion of Wnt signaling and BMP2/6 by ECs functions on hepatocytes to control hepatic zonation and iron metabolism (10)(11)(12)(13), respectively.Loss of these signals results in liver dysfunction.Recent advances in RNA-sequencing technology and algorithms have enabled a better understanding of the cellular composition and potential cellular interactions of the liver (14).However, the speculative cues still need experimental validation, and more importantly, it is still unknown what consequences happen if specific cell-cell interactions are interrupted.
HSCs reside in the space of Disse and directly interact with HCs, ECs, and KCs (15).HSCs are well known for their pathological role in liver fibrosis, as HSC activation is a critical step in the development of chronic liver disease (16).However, in addition to the storage of retinoids, the physiological role of HSCs in the normal liver remains unclear.ALK1, encoded by the activin A receptor-like type 1 (Acvrl1) gene, is the receptor of BMP9 and BMP10.Its mutations result in hereditary hemorrhagic telangiectasia (HHT), which causes vascular malformations in multiple organs, including the liver (17).It has been reported that liver sinusoidal EC-specific Alk1deficient mice exhibit severe liver vascular malformations (18), which mimic the liver symptoms observed in HHT patients.We and others demonstrated that ALK1 is needed for KC survival and identity, and the ability of KCs from Alk1-deficient mice to capture bacteria from the bloodstream was significantly impaired (19,20).BMP9 is mainly expressed in the liver.Using antibody staining, BMP9 was reported to be expressed in human hepatocytes (21).Recently, it has been reported that BMP9 and BMP10 have been shown to be mainly expressed by murine HSCs (20,22).However, no study has systematically examined which liver cell types secrete BMP9 and BMP10 and regulate ECs and KCs via these two paracrine factors.Here, we generated mice with conditional knockout of both Bmp9 and Bmp10 in different liver cell types and found that HSCs are the functional source of BMP9 and BMP10.Interestingly, conditional deletion of Bmp9/10 from HSCs (Bmp9/10 HSC-KO ) mice not only resulted in loss of KC and EC identity but also affected HC phenotypes in an indirect way, which ultimately resulted in liver dysfunction, suggesting that HSCs play a central role in orchestrating the interactions between different liver cell types under physiological conditions.

HSCs are the functional source of BMP9 and BMP10 in the liver
To generate mice lacking Bmp9 and Bmp10 in different liver cell types, Bmp9 fl/fl and Bmp10 fl/fl mice were prepared and mated with Alb Cre , Clec4f Cre , Tie2 Cre , or Lrat Cre mice, which have been extensively used in the literature to produce HCs, KCs, ECs, and HSC conditional knockout mice (12,19,(23)(24)(25), respectively.Because BMP9 and BMP10 play a redundant role in regulating KCs (19), conditional double knockout mice (Bmp9 fl/fl Bmp10 fl/fl Alb Cre , Bmp9 fl/fl Bmp10 fl/fl Tie2 Cre , Bmp9 fl/fl Bmp10 fl/fl Clec4f Cre and Bmp9 fl/fl Bmp10 fl/fl Lrat Cre ) were generated.To systematically investigate which liver cell type produces BMP9 and BMP10, we acquired liver tissues from these mice.
Quantitative PCR analysis revealed that Bmp9 and Bmp10 mRNA expression was significantly reduced in the liver tissues from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre (Bmp9/10 HSC-KO ) mice compared to other Cre deleter mice (Figure 1a).Using Cre-mediated reporter mice, we demonstrated that Lrat Cre mice can specifically target HSCs (Figure 1b).BMP10 is mostly expressed in the heart, followed by the liver (26) .However, Bmp10 expression in the heart was not affected in Bmp9/10 HSC-KO mice (Figure 1c).ALK1 is the receptor of BMP9/10.To investigate which liver cell types express ALK1 and respond to BMP9/10, we prepared transgenic ALK1 reporter mice, which can faithfully reflect the expression of ALK1.We found that ALK1 was expressed on ECs and KCs but not HCs and HSCs (Figure 1d), suggesting that BMP9/10 can directly function on ECs and KCs via ALK1.Accordingly, we found that surface markers of liver F4/80 + macrophages from Bmp9/10 HSC-KO mice, including CLEC4F, Tim4 and VSIG4, were almost not expressed (Figure 1e and 1f).In contrast, these markers were not affected in the liver F4/80 + macrophages from Bmp9 fl/fl Bmp10 fl/fl Alb Cre , Bmp9 fl/fl Bmp10 fl/fl Tie2 Cre and Bmp9 fl/fl Bmp10 fl/fl Clec4f Cre mice (Supplementary Figure 2).LYVE1 and CLEC2, the surface markers of CD31 + hepatic ECs, were also reduced in the livers of Bmp9/10 HSC- KO mice (Figure 1g).In addition, we also determined the phenotypes of KCs and ECs from Bmp10 fl/fl Lrat Cre mice to exclude the possibility that the altered phenotypes observed in Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice were due to Cre-mediated cytotoxicity (27) and found that these cells were not affected (Supplementary Figure 3).Collectively, these results suggest that HSCs are the functional source of BMP9 and BMP10 in the liver to regulate KCs and ECs.

The differentiation of liver macrophages from Bmp9/10 HSC-KO mice is blocked in a status of monocyte-derived macrophages
To acquire more information about transcriptome changes in liver macrophages in the Bmp9/10-deficient liver microenvironment, we performed bulk RNA-seq analysis on sorted CD45 + Ly6C -F4/80 + CD64 + liver macrophages from Bmp9/10 HSC-KO mice compared with their controls.Principal component analysis (PCA) suggested that liver macrophages from Bmp9/10 HSC-KO mice were different to that of their controls (Figure 2a).We found that 2424 genes were differentially expressed (FC>2, p-adj<0.05),among which 1313 genes were upregulated and 1111 genes were downregulated in liver macrophages from Bmp9/10 HSC-KO mice compared to their control mice (Figure 2b).
From these differentially expressed genes, we found that genes associated with the identity of KCs (28), such as Fabp7, Cd5l, Cdh5 and Clec4f, were significantly downregulated in liver macrophages from Bmp9/10 HSC-KO mice (Figure 2c).Id1 and Id3, target genes of ALK1 signaling and important regulators of KC differentiation, were also significantly reduced in liver macrophages from Bmp9/10 HSC-KO mice (Figure 2d).
We also observed downregulation of embryonic KC marker Timd4 (encoding Tim4) expression and upregulation of monocytic Cx3cr1 and Ccr2 expression (15,29) in liver macrophages from Bmp9/10 HSC-KO mice (Figure 2d), suggesting that embryonic KCs were lost and replaced by monocyte-derived macrophages/KCs, as a reduced number of liver macrophages was observed in Bmp9/10 HSC-KO mice (Figure 2e).
Upon KC loss, blood monocytes are recruited into the liver, immediately differentiate into F4/80 + monocyte-derived macrophages (MoMs) and then gradually acquire the expression of CLEC2, CLEC4F and VSIG4 to become monocyte-derived KCs (MoKCs) with time (15,28).CLEC2, encoded by the Clec1b gene, was thought to be one of the earliest markers gained by monocyte-derived KCs and can be used to discriminate MoKCs from MoMs (30).In addition, TREML4 expression is immediately upregulated upon blood monocyte entry into the liver (25).Interestingly, we found that nearly all liver macrophages from Bmp9/10 HSC-KO mice did not express CLEC2 and TREML4, but liver macrophages from Alk1 fl/fl Vav1 Cre and Smad4 fl/fl Vav1 Cre mice did (Figure 2f), suggesting that the differentiation of liver macrophages from Bmp9/10 HSC-KO mice is completely blocked in a status of monocytederived macrophages, and the differentiation defect in liver macrophages from Bmp9/10 HSC-KO mice is more pronounced than that from Alk1 fl/fl Vav1 Cre and Smad4 fl/fl Vav1 Cre mice.Our previous report has shown that Alk1 and Smad4 were efficiently deleted in KCs from Alk1 fl/fl Vav1 Cre and Smad4 fl/fl Vav1 Cre mice, respectively (19).To understand this discrepancy, we performed bulk RNA-seq analysis on sorted liver macrophages from Smad4 fl/fl Vav1 Cre mice, as Smad4 functions as a common Smad needed for transcriptional regulation in response to ALK signaling.In addition to ALK1, BMP9 and BMP10 also bind ALK2 and ALK6, respectively, with low affinity to transduce signals (31).Principal component analysis suggested that Smad4-deficient KCs were distinguishable from Smad4-sufficient KCs (Figure 2g).Analysis of the DE genes between Smad4-sufficient and Smad4-deficient KCs identified 1340 DE genes, including many KC signature genes (Figure 2h and 2i).Comparison of the difference between the DE genes in liver macrophages from Bmp9/10 HSC-KO and Smad4 fl/fl Vav1 Cre mice revealed that 1659 of 2424 DE genes in liver macrophages from Bmp9/10 HSC-KO mice did not overlap with DE genes in MoKCs from Smad4 fl/fl Vav1 Cre mice (Figure 2j), which is consistent with the phenomenon described above, indicating that in addition to Smad4, alterations in other molecular mechanisms may affect the phenotypes of liver macrophages in Bmp9/10 HSC-KO mice.
Transcription factors play an important role in the differentiation of tissue-resident macrophages.We compared the nonoverlapping transcription factors between the DE genes in liver macrophages from Bmp9/10 HSC-KO and Smad4 fl/fl Vav1 Cre mice and found mice (Figure 2k and Supplementary Table 1).Among these transcription factors, we focused on Rxra because it has been reported to regulate the development of many embryonic tissue-resident macrophages, including KCs (32), but it is unclear whether RXRa regulates the differentiation of blood monocytes to MoKCs.To test this hypothesis, we prepared Rxra fl/fl Csf1r Cre Clec4f DTR mice, in which DT injection can delete KCs and monocytes recruited into the liver to differentiate into liver macrophages, which can be targeted by the Csf1r Cre strain to eliminate the Rxra gene (Supplementary Figure 4a).We found that RXRα deficiency significantly inhibited the differentiation of blood monocytes to MoKCs, as the proportion of CLEC2 + and VSIG4 + macrophages was significantly reduced in Rxra-deficient macrophages compared with their controls (Supplementary Figure 4b and 4c), suggesting that the differentiation defect of liver macrophages from Bmp9/10 HSC-KO mice may be partially caused by Rxra downregulation.

Endothelial cells from Bmp9/10 HSC-KO mice are transdifferentiated into continuous endothelial cells
Next, to understand how ECs were affected in Bmp9/10 HSC-KO mice, we performed bulk RNA-seq analysis on sorted CD45 -CD31 + endothelial cells from Bmp9/10 HSC-KO mice compared with their controls.We found that 2692 genes were differentially expressed (FC>2, p-adj<0.05),among which 1534 genes were upregulated and 1158 genes were downregulated in ECs from Bmp9/10 HSC-KO mice compared to their control mice (Figure 3a and 3b).Both GATA4 and MAF act as master regulators of hepatic sinusoidal differentiation (33,34), and their expression in ECs can be induced by BMP9 in vitro (34,35).In our transcriptomic data, Gata4 and Maf expression in endothelial cells from Bmp9/10 HSC-KO mice were significantly reduced compared with that in their controls (Figure 3c).Consistent with the phenotypes of Gata4/Maf-deficient endothelial cells, the expression of sinusoidal endothelial cell-specific markers (34) such as Stab 2, Fcgr2b and Mrc1 was significantly downregulated (Figure 3d), while the expression of continuous endothelial cell-specific markers (33) such as Cd34, Emcn, Apln and Aplnr was upregulated in ECs from Bmp9/10 HSC-KO mice (Figure 3e).
Immunostaining revealed obviously increased expression of CD34 and EMCN, and decreased expression of LYVE1 (Figure 3f and 3g).Folw cytometry also confirmed the reduced expression of C-Maf (encoded by Maf), CD32b (encoded by Fcgr2b), CD117 (encoded by Kit) and CD206 (encoded by Mrc1) (Figure 3h and 3i).In addition, genes associated with large arteries (34) were significantly upregulated in ECs from Bmp9/10 HSC-KO mice (Figure 3j).These results suggested that hepatic ECs in Bmp9/10 HSC-KO mice lost their identity and were transdifferentiated into continuous endothelial and large arterial cells.
Liver metabolic zonation and iron metabolism are disrupted in Bmp9/10 HSC-KO mice Angiocrine factors Wnts and BMP2/6 signal were shown to control hepatic metabolic zonation and iron metabolism (10)(11)(12)(13), respectively.The expression of these genes was significantly reduced in ECs from Bmp9/10 HSC-KO mice (Figure 4a).Prussian blue staining revealed that 5 out of 7 Bmp9/10 HSC-KO mice showed mild iron overload in the focal liver (Figure 4b).Wnt pathway activation induces the expression of zonated proteins such as glutamine synthetase.Immunostaining confirmed the disappearance of glutamine synthetase (GS) in the livers of Bmp9/10 HSC-KO mice (Figure 4c), suggesting that liver metabolic zonation was disrupted in the Bmp9/10-deficient liver microenvironment.
To acquire more information regarding how the livers were affected under the condition of BMP9/10 deficiency, we performed bulk RNA-seq analysis on whole liver samples of Bmp9/10 HSC-KO mice and their control mice at the age of 12 weeks.We analyzed the expression of Wnt signal target genes enriched in the pericentral zone and found that most pericentral Wnt signal targets were significantly downregulated in the livers of Bmp9/10 HSC-KO mice, including Glul, Oat, Cyp2e1, Cyp1a2, Cldn2, Axin2, Lect2, Lgr5, Tbx3 and Rnf43 (Figure 4d), supporting the fact that the liver microenvironment in Bmp9/10 HSC-KO mice lacks Wnt2, Wnt9b and Rspo3 signaling.In addition, BMP9 and BMP10 play a redundant role in regulating liver metabolic zonation, as glutamine synthetase was still expressed in the livers of Bmp9 fl/fl Lrat Cre and Bmp10 fl/fl Lrat Cre mice (Figure 4e and 4f).

Liver fibrosis occurs in Bmp9/10 HSC-KO mice
Liver weight and the liver/body weight ratio were significantly reduced in Bmp9/10 HSC- KO mice at the age of 12W compared with control mice (Figure 5a), suggesting that liver development is impaired.Picrosirius red staining revealed increased collagen deposition in the livers of Bmp9/10 HSC-KO mice at the age of 12W (Figure 5b), suggesting that these mice presented liver fibrosis.This phenotype was also seen in the livers of double-knockout mice for Bmp9 (constitutive) and Bmp10 (Bmp10 fl/fl Rosa26 CreERT2 , tamoxifen inducible) (36).PDGFB is a profibrotic cytokine that results in HSC activation and liver fibrosis.In the transcriptomic data, we found that Pdgfb was significantly upregulated in the liver tissue, ECs and liver macrophages of Bmp9/10 HSC-KO mice (Figure 5c).Immunostaining revealed that the expression of desmin and type I collagen, two markers of HSC activation, was upregulated in the liver tissues of Bmp9/10 HSC-KO mice (Figure 5d), suggesting that HSC activity was increased.
Endothelial Gata4 deficiency leads to sinusoidal capillarization (37).In line with this report, the expression of collagen IV, a marker of EC capillarization, was significantly increased in the livers of Bmp9/10 HSC-KO mice, as revealed by immunostaining (Figure 5e) and the transcriptomic data of ECs (Figure 5f).These results suggested that liver fibrosis in Bmp9/10 HSC-KO mice was due to EC capillarization and HSC activation.

Discussion
In this study, we demonstrated that HSCs convey BMP9/10 signals to KCs and ECs to maintain their identity.Loss of BMP9/10 in the liver microenvironment resulted in liver macrophage differentiation defects and a shift in hepatic EC differentiation toward continuous ECs and large arterial cells.Under this condition, physiological crosstalk between ECs and HCs was disrupted and pathological interaction of LSECs/KCs with HSCs was enhanced, which ultimately altered the whole liver microenvironment and resulted in liver dysfunction.
Recently, we and others have demonstrated that ALK1 is essential for KC differentiation, self-renewal, identity, and endocytic ability (19,20).Here, we undertook systematic genetic manipulation to scrutinize the functional source of BMP9/10 in the liver.We found that instead of ECs, KCs and HCs, HSCs are the predominant source of BMP9/10 critical for KC homeostasis.BMP9 is mainly expressed in the liver, whereas BMP10 is mostly expressed in the heart and weakly expressed in the liver.We found that BMP10 expression in the heart was not affected in Bmp9/10 HSC-KO mice, suggesting that the BMP10 signal needed for KCs and ECs was provided by the local liver rather than circulating BMP10 secreted by the heart.Indeed, ALK1 is weakly expressed in KCs and ECs, as the fluorescence signal emitted by KCs and LSECs from Alk1 tdTomato reporter mice was weak when detected by flow cytometry, suggesting that the concentration of circulating BMP10 may not be enough to activate ALK1 signaling on hepatic ECs and KCs to maintain their phenotypes.
Recently, Schmid et al. used Stab 2 Cre mice to delete Alk1 in liver sinusoidal ECs to model liver HHT (18).They found that upon deletion of Alk1, liver sinusoidal ECs lost their identity and were transdifferentiated into large arterial cells.In addition, angiocrine Wnt signaling, such as Wnt2, Wnt9b, and Rspondin3, was impaired in Alk1deficient liver sinusoidal ECs, leading to the disruption of liver metabolic zonation (18).These phenomena were also observed in the Bmp9/10 HSC-KO mouse model.However, liver fibrosis observed in HHT patients did not occur in the livers of Alk1 fl/fl Stab 2 Cre mice (18).Indeed, our transcriptomic data showed that Stab 2 expression was significantly reduced in hepatic ECs of Bmp9/10 HSC-KO mice , which was consistent with the data from Bmp9 KO 129/Ola mice (35), suggesting that ALK1 signaling regulated STAB2 expression.Thus, it is possible that heterozygosity for the Alk1 deletion in ECs may result in reduced expression of STAB2, and the expression of Cre recombinase driven by the Stab 2 promoter accordingly decreases, so that the Alk1-floxed allele in ECs cannot be efficiently deleted by Cre-mediated cyclization, which may explain why liver fibrosis did not occur in Alk1 fl/fl Stab 2 Cre mice.In view of liver fibrosis, Bmp9/10 HSC-KO mice may be a better mouse model to mimic the liver pathology observed in HHT patients (35).In addition, the liver macrophage niche is composed of HCs, ECs and HSCs (15,25).Thus, the altered phenotypes of HCs and ECs from Bmp9/10 HSC-KO mice may be another reason for the differentiation defect of liver macrophages observed in Bmp9/10 HSC-KO mice.
Gata4-deficient mice spontaneously develop liver fibrosis (37).Single-cell RNA-seq analysis of healthy and cirrhotic human liver samples revealed that Gata4 expression was also reduced in ECs from cirrhotic human liver samples (34), suggesting that reduced expression of Gata4 may contribute to murine and human liver fibrosis.Given that BMP9/10 regulated the expression of Gata4 in murine ECs, it is possible that the reduced expression of GATA4 may be due to low expression of BMP9/10 in liver cirrhosis.However, the expression of Bmp9/10 was not detected by single-cell RNAseq of either healthy or cirrhotic human livers (38)(39)(40)(41)

Liver cell suspension preparations, flow Cytometry and antibodies
Briefly, the liver was perfused with approximately 20 ml HBSS containing heparin and EDTA, followed by perfusion with digestion buffer containing collagenase type IV (Sigma) and DNase I (Sigma) for 5 min.The digested livers were then disrupted, and pipetted up and down, and the cell suspension was filtered through a 70μm cell strainer.
Parenchymal cells were separated from nonparenchymal cells by centrifugation at 50g for 5 min for two times.ECs and liver macrophages were further enriched by iodixanol gradient (OptiPrep) as previously described (42).

Quantitative Real-Time PCR
Total RNA was isolated with Trizol (Thermo) and cDNA was synthesized with Prime Script RT Reagent Kit (Takara).Quantitative PCR was performed with a SYBR Green PCR kit (Toyobo, Japan) in CFX Connect Real-time PCR detection system (Bio-Rad).
The specific qPCR primers used are listed in Supplemental Table 2.

RNA-seq analysis
Total RNA of ECs and liver macrophages was isolated with RNeasy plus Mini Kit (Qiagen).Total RNA of the liver tissues was extracted with Trizol (Thermo).Library preparation was performed with UltraTM RNA Library Prep Kit (Illumina) followed by RNA sequencing using NovaSeq 6000 (Illumina).All samples passed quality control based on the results of FastQC.Reads were aligned to the mouse genome assembly mm10.GRCm38 using HISAT2.The aligned reads were counted using FeatureCounts.
Differentially expressed genes were determined with DESeq2 program at p-adj (adjusted p value) < 0.05 and logFC (log fold change) > 1 or logFC where indicated.
All RNA sequencing experiments were performed with three independent biological replicates.The heatmap and volcano maps of DEGs were visualized using the "ggplot2" and "pheatmap" R package.

Immunohistochemistry and Immunofluorescence
For immunohistochemistry, the livers were fixed in 4% PFA for 12 hours, and then was dehydrated with gradient ethanol and embedded in paraffin.Sections (4μm) were cut on a Finesse 325 (Thermo Scientific) and adhered to adhesion microscope slides.The paraffin sections were dewaxed and hydrated with xylene and gradient ethanol, and then boiled with antigen retrieval solution to repair antigen.Endogenous peroxidase was blocked with 0.3% hydrogen peroxide for 30 minutes.Paraffin sections were blocked with phosphate-buffered saline (PBS) containing 3% horse serum for 1h, and incubated with Avidin and Biotin for 15 minutes, followed by staining with anti-Collagen IV antibody (ab19808, Abcam) overnight at 4℃.Sections were incubated with secondary antibody, ABC reagent and NovaRed peroxidase substrate of Vector-kit (Vector Laboratories) following with the manufacture.Images were acquired on a C10730-12 NanoZoomer with a 20× /NA 0.75 objective lens.
For immunofluorescence, the fixed livers were dehydrated in PBS solutions containing 15% and 30% sucrose for 12 hours, and embedded in OCT compound (Sakura Finetek).
The samples were cut into 20-micron slices using the Cryotome FSE (Thermo Scientific) and then attached to Adhesion Microscope slides.Frozen liver sections were then blocked for 2 hours with 1% bovine serum albumin, 0.3 M glycine, and 10% donkey or goat serum before being stained overnight at 4°C with primary antibodies diluted in PBS with 0.2% Tween, secondary antibodies, and Hoechst for 2h at room temperature.

Image processing and analysis
Quantitative analyses of all images were measured using the Fiji package for ImageJ software (NIH) and the Imaris (Oxford Instruments).For coverage measurements of KCs and HSCs, surface modules in Imaris were used to create CD64 and F4/80

C-Maf
due to the low abundance of Bmp9/10 mRNA in HSCs.Thus, the correlation between Bmp9/10 expression and liver fibrosis needs further investigation.Collectively, we demonstrated that BMP9/10 sourced from HSCs function on KCs and ECs to maintain their identity, which further directs HCs to wire metabolic function to meet the organ's needs.Upon disruption of this crosstalk, liver dysfunction ensued.This work illustrates the complexity underlying the crosstalk between different liver cell types.Methods Mice Bmp9 fll+ (Stock No: T026617), Bmp10 fll+ (Stock No: T008480), Tie2 Cre (Stock No: T003764), Lrat Cre (Stock No: T006205) and H11 reporter mice (Stock No: T006163) were obtained from GemPharmatech (Nanjing, China).Alk1 tdTomato mice that an expression cassette encoding tdTomato and a self-cleaving 2A peptide was inserted into the start codon of Alk1 gene were prepared by MODEL ORGANISMS (Shanghai, China).Alb Cre mice (Stock No: 003574) and Csf1r Cre (Stock No: 029206) mice were obtained from Jackson Laboratory.Smad4 fl/fl and R26 YFP mice were kindly provided by Dr. Xiao Yang (Beijing Institute of Lifeomics).Alk1 fl/fl mice were kindly provided by Dr. Zhihong Xu (Fudan University, Shanghai).Vav1 Cre mice were kindly provided by Dr. Bing Liu (Fifth Medical Center of Chinese PLA General Hospital).Clec4f Cre/DTR strain mice were established by Nanjing BioMedical Research Institute of Nanjing University (NBRI) using CRISPR/Cas9-mediated genome editing on a C57/BL/6J background.Mice and their littermates were used between eight-to sixteen-week-old unless otherwise specifically indicated.All mice are C57BL/6 background and were maintained at the SPF facilities of the Beijing Institute of Lifeomics.All experimental procedures in mice were approved by the Institutional Animal Care and Use Committee at the Beijing Institute of Lifeomics.
biomarkers.A threshold was determined by the algorithm for background subtraction and manually set to match the cell signals.The number of cells per unit area was accessed by parameters in Imaris.IHC and Picrosirius red staining can be assessed by color deconvolution and the thresholding tool in ImageJ.Statistics Data are presented as mean ± SEM.All statistics analyses were performed in GraphPad Prism (GraphPad Software).Statistical significance was assessed by the Mann Whitney test.*p < 0.05; **p < 0.01; ***p < 0.001; ns, not significant.Each symbol represents an individual mouse.Number of animals is indicated as ''n''.

Figure legends Figure 1 .
Figure legends

Figure 2 .
Figure 2. The differentiation of liver macrophages was inhibited in Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice.(a, b) PCA and volcano plot of the RNA-seq data of sorted CD45 + Ly6C -F4/80 + CD64 + liver macrophages from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 8-10 weeks (n=3/group).Genes upregulated and downregulated are shown in red and blue, respectively (fold change [FC] > 2, adjusted p[p-adj] < 0.05).(c) Heatmap showing signature genes expressed differentially in liver macrophages from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.(d) Expression counts of indicated genes in

Figure 3 .
Figure 3. Endothelial cells from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice are transdifferentiated to continuous endothelial cells.(a, b) PCA and volcano plot of the RNA-seq data of sorted hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 8-10 weeks (n=3/group).Genes upregulated and downregulated are shown in red and blue, respectively (fold change [FC] > 2, adjusted p[p-adj] < 0.05).(c) Expression counts of Gata4 and Maf genes in ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.***p-adj < 0.001.(d, e) Heatmap showing sinusoidal EC-associated genes (d) and continuous EC-associated genes (e) expressed differentially in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.(f) Immunofluorescence images in liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 12 weeks (n=4/group).CD34 + area was measured (right).Scale bars: 100μm.(g) Immunofluorescence images in liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 8 weeks (n=4-7/group).LYVE1 + and EMCN + area was measured (right).Scale bars: 100μm.(h, i) Flow cytometric expression of CD32b, C-Maf, CD206 and CD117 in hepatic ECs from

Figure 4 .
Figure 4. Liver metabolic zonation and iron metabolism are disrupted in Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice.(a) Heatmap showing the indicated genes expressed differentially in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.(b) Prussian blue in livers from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their controls at the age of 12 weeks (n=7/group).(c) Immunofluorescence images in liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their controls at the age of 12 weeks (n=5/group).Scale bars: 200μm.(d) Heatmap showing the indicated genes expressed differentially in liver tissues from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 12 weeks.(e, f) Immunofluorescence images in liver sections from Bmp9 fl/fl Lrat Cre (e) and Bmp10 fl/fl Lrat Cre (f) mice and their controls (n=2/group) at the age of 13-21 weeks (n=2/group).Scale bars: 200μm.Results represent the mean± SEM.**P < 0.01, by Mann-Whitney test (c, e, f).

Figure 5 . 2 ✱✱✱Figure 1 .
Figure 5. Liver fibrosis occurs in Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice.(a) Liver weight and liver to body weight ratio of Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their littermate controls at the age of 12 weeks.(b) PSR staining of liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their controls at the age of 12 weeks (n=9-10/group).Scale bars: 100μm.(c) Expression counts of Pdgfb gene in liver tissues (left), endothelial cells (middle) and liver macrophages (right) from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their control mice.***p-adj < 0.001.(d) Immunofluorescence images of Desmin and collagen I in liver sections in liver tissues from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 12-28 weeks (n=4-6/group).(e) Collagen IV immunohistochemistry staining of liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their controls at the age of 12 weeks (n=4-5/group).Scale bars:100μm.(f) Expression counts of indicated genes in endothelial cells from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.***p-adj < 0.001.Results represent the mean± SEM.*P < 0.05, **P

Figure 2 .
Figure 2. The differentiation of liver macrophages was inhibited in Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice.(a,b) PCA and volcano plot of the RNA-seq data of sorted liver macrophages from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 8-10 weeks (n=3/group).Genes upregulated and downregulated are shown in red and blue, respectively (fold change [FC] > 2, adjusted p[p-adj] < 0.05).(c) Heatmap showing signature genes expressed differentially in liver macrophages from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.(d) Expression counts of indicated genes in liver macrophages from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.***p-adj< 0.001.(e) Immunofluorescence images of F4/80 + and CD64 + liver macrophages in sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their controls at the age of 12 weeks (n=5/group).Liver macrophages number was measured (right).(f) Representative flow cytometric expression of CLEC2 and TREML4 in liver macrophages from the indicated mice at the age of 8-12 weeks (n=3-4/group).(g,h) PCA and volcano plot of the RNA-seq data of sorted liver macrophages from Smad4 fl/fl Vav1 Cre and control mice at the age of 8-10 weeks (n=3/group).Genes upregulated and downregulated are shown in red and blue, respectively (fold change [FC] > 2, adjusted p[padj] < 0.05).(i) Heatmap showing signature genes expressed differentially in liver macrophages from Smad4 fl/fl Vav1 Cre and control mice.(j, k) Venn diagram showing DE genes (j) and transcription factors (k) specific to liver macrophages of Bmp9 fl/fl Bmp10 fl/fl Lrat Cre , Smad4deficient liver macrophages or shared between both mac populations.Results represent the mean ± SEM. *P < 0.05, by Mann-Whitney test (e).

iFigure 3 .
Figure 3. Endothelial cells from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice are transdifferentiated to continuous endothelial cells.(a,b) PCA and Volcano plot of the RNA-seq data of sorted hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 8-10 weeks (n=3/group).(c) Expression counts of Gata4 and Maf genes in ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.***p-adj < 0.001 (d, e) Heatmap showing sinusoidal EC-associated genes (d) and continuous EC-associated genes (e) expressed differentially in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.(f) Immunofluorescence images in liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 12 weeks (n=4/group).CD34 + area was measured (right).Scale bars: 100μm.(g) Immunofluorescence images in liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 8 weeks (n=4-7/group).LYVE1 + and EMCN + cells were measured (right).Scale bars: 100μm.(h,i) Flow cytometric expression of CD32b, C-Maf, CD206 and CD117 in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and littermate controls at the age of 8-10 weeks (n=4-6/group).(j) Heatmap showing the indicated genes expressed differentially in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.Results represent the mean± SEM.*P < 0.05 and **P <0.01, by Mann-Whitney test (f-i).

Figure 4 .Figure 5 .
Figure 4. Liver metabolic zonation and iron metabolism are disrupted in Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice.(a) Heatmap showing the indicated genes expressed differentially in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.(b) Prussian blue in livers from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their controls at the age of 12 weeks (n=7/group).(c) Immunofluorescence images in liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their controls at the age of 12 weeks (n=5/ group).Scale bars: 200μm.(d) Heatmap showing the indicated genes expressed differentially in liver tissues from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 12 weeks.(e, f) Immunofluorescence images in liver sections from Bmp9 fl/fl Lrat Cre (e) and Bmp10 fl/fl Lrat Cre (f) mice and their controls at the age of 13-21 weeks (n=2/group).Scale bars: 200μm.Results represent the mean± SEM.**P < 0.01, by Mann-Whitney test (c, e, f).
Bmp9 fl/fl Bmp10 fl/ fl Bmp9 fl/fl Bmp10 fl/ fl Alb Cre Bmp9 fl/fl Bmp10 fl/ fl Bmp9 fl/fl Bmp10 fl/ fl Tek Cre Bmp9 fl/fl Bmp10 fl/ fl Bmp9 fl/fl Bmp10 fl/ fl Clec4f Cre fl/fl Lrat Cre and control mice at the age of 8-10 weeks (n=3/group).Genes upregulated and downregulated are shown in red and blue, respectively (fold change [FC] > 2, adjusted p[p-adj] < 0.05).(c) Heatmap showing signature genes expressed differentially in liver macrophages from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.(d) Expression counts of indicated genes in liver macrophages from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.***p-adj< 0.001.(e) Immunofluorescence images of F4/80 + and CD64 + liver macrophages in sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre mice and their controls at the age of 12 weeks (n=5/group).Liver macrophages number was measured (right).Heatmap showing signature genes expressed differentially in liver macrophages from Smad4 fl/fl Vav1 Cre and control mice.(j, k) Venn diagram showing DE genes (j) and transcription factors (k) specific to liver macrophages of Bmp9 fl/fl Bmp10 fl/fl Lrat Cre , Smad4deficient liver macrophages or shared between both mac populations.Results represent the mean ± SEM. *P < 0.05, by Mann-Whitney test (e).
fl/fl Lrat Cre and control mice at the age of 8-10 weeks (n=3/group).(c) Expression counts of Gata4 and Maf genes in ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.***p-adj < 0.001 (d, e) Heatmap showing sinusoidal EC-associated genes (d) and continuous EC-associated genes (e) expressed differentially in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.(f) Immunofluorescence images in liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 12 weeks (n=4/group).CD34 + area was measured (right).Scale bars: 100μm.(g) Immunofluorescence images in liver sections from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice at the age of 8 weeks (n=4-7/group).LYVE1 + and EMCN + cells were measured (right).Scale bars: 100μm.(h,i) Flow cytometric expression of CD32b, C-Maf, CD206 and CD117 in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and littermate controls at the age of 8-10 weeks (n=4-6/group).(j) Heatmap showing the indicated genes expressed differentially in hepatic ECs from Bmp9 fl/fl Bmp10 fl/fl Lrat Cre and control mice.Results represent the mean± SEM.*P < 0.05 and **P <0.01, by Mann-Whitney test (f-i).