The development of PMPs
We aimed to build a high-throughput organoid culture platform (500 microwells in 1cm2), and Fig. 1a summarized the overall strategy. Briefly, we developed PMPs by transcribed silicon mold structures accurately. As the microwell of PMP had a special inverted pyramid structure (the side length of tightly arranged microwells is ~ 450 µm and the depth is ~ 400 µm, as shown in Fig. 1b and Fig. 1c), which was difficult to broadly transcribe in PDMS using conventional molding machine. We achieved highly accuracy by employing three-step transcription method (see Experimental Procedures). Further, PMPs could be punched to fit with any type of multi-well culture plates, flasks, and petri dishes. Hepatic organoids were generated from only HepaRG cells for PMP characterization. The size of organoids could be controlled by adjusting input cell density (as shown in Fig. 1d). 90, 110, 130, and 150 µm diameter organoids were generated by seeding 400, 600, 800, and 1000 cells in each microwell in day4, respectively. Organoids' controllable size allows more flexibility in optimizing printing parameters.
Compared to low attachment plates, organoids generated by PMPs were consistent in size and shape, uniform within and between independent experiments (Fig. 2a). Uniform generation of organoids made it possible to construct reproducible liver tissues. Regulated changes in reactive oxygen species (ROS) activities are important to maintain normal development of hepatoblasts23. As shown in Fig. 2b, intracellular ROS level of hepatic organoids was significantly enhanced in PMP groups compared to Aggrewell™ groups. Intracellular ROS activities are influenced by changes in environmental oxygen concentration24. The high oxygen permeability of PDMS makes it ideal for use in hepatic organoid culture. Increasing evidence suggested that liver homeostasis and function depend on mechanical properties (e.g., stiffness and shear stress) 25,26. We optimized curing temperature and time of PDMS to balance transcription accuracy and mechanical properties. The stiffness of PMPs was 3.3 ± 0.07 MPa. Aggrewell™ plates were stiffer than PMPs by three orders of magnitude (Fig. 2c). Hepatic organoids cultured in PMPs with lower stiffness and high oxygen permeability demonstrated higher cell activity, upregulated ALB and MRP2 gene expression than Aggrewell™ groups, suggesting that PMPs are more suitable for hepatic organoids generation (Fig. 2d, e).
Self-assembly and multicellular crosstalk of HOBBs
Throughout the development of the liver, multicellular crosstalk occurs. Hepatoblasts and endothelial cells interact during liver bud development27. Emergence of liver bud requires inductive signals from endothelial cells. Prior to liver bud expansion, endothelial cells line the nascent sinusoids and physically interact with the hepatoblasts21. Hepatoblasts were also identified as a positive stimulator of sinusoid morphogenesis and maturation28. For the development of liver, crosstalk between hepatoblasts and endothelial cells is essential. Recent advances in developing in vitro liver tissue models show that function of hepatoblasts is altered when co-culturing with endothelial cells13,29,30. Although HUVECs are not native to liver, many have used them due to their ease of access and established protocols31. To generate HOBBs, we co-cultured HepaRG cells and HUVECs with PMPs upon mixture ratio (Supplementary Fig. 1). Comparing cell activity and gene expression of each group, we selected the group that had a cell ratio of HepaRG cells to HUVECs of 2:1 (Supplementary Fig. 1).
By comparing hepatic organoids and HOBBs in day1 (Fig. 3a, Supplementary Fig. 2a), it was found that HUVECs accelerated organoids formation, implying interaction and crosstalk between HepaRG cells and HUVECs. Cell aggregation to form organoids is a multi-step process characterized by successive changes in cell morphology and gene expression, including regulation of actin cytoskeleton, cell-cell assembly, remodeling of extracellular matrix, and cell-matrix adhesion32. Genes regulating actin cytoskeleton signaling pathways such as MYL6, ZYX, ARPC4, TMBS4X, CFL1, PFN1, ACTG1, and ACTGB were upregulated in HOBB group versus 2D group, proving that cells underwent morphological changes and regulation of actin cytoskeleton (Fig. 3b, Supplementary Fig. 2b). Cell–cell adhesion through molecules such as ICAM1, JAM1, JAM2 and Nectin4 genes was required for organoids formation, their upregulation suggesting tight assembly of cells in HOBBs group (Fig. 3b). CLDNs are major structural components of tight junctions (TJs) and regarded as cellular backbone, tighten the paracellular cleft of epithelium and endothelia against unwanted passage of solutes. Epithelial and endothelial cell-cell junction genes CLDN3, CLDN 4, CLDN 5 and CLDN 14 expression of HOBBs were upregulated (Fig. 3b). The above change indicated HUVECs and HepaRG cells within HOBBs were tightly connected, establishing the basis for crosstalk between them. Extracellular matrix (ECM) secretion played a significant role in organoid morphology33. HOBBs exhibited enhanced matrix secretion and higher expression levels of COL9A3, COL9A2, COL6A3, and LAMB3, and the corresponding integrins ITGA2, ITGA5, ITGB5 and other ECM receptors MMP9 were also upregulated (Fig. 3b). According to these findings, HOBBs assembled dynamically, remodeled the extracellular matrix, formed tight cell-cell and cell-matrix connections. HOBBs possessed more bionic niches than 2D environment.
Constant confocal wide-field imaging and fluorescence-based quantitative assessment of collected HOBBs generated from multiple stages (day 0–4) revealed the migration of HUVECs-GFP (Fig. 3c). HOBBs’ fluorescence intensity decreased inside and increased at the edge. In order to elucidate HepaRG cells and HUVECs locations within HOBBs, we performed co-immunostaining with ALB and CD31. ALB+ HepaRG cells were found within HOBBs, and CD31+ HUVECs surrounded HOBBs (Fig. 3d). Gene set enrichment analysis (GSEA) of RNAseq data showed upregulation of VEGF signaling pathway of HOBB groups (Fig. 3e). Due to the VEGF signaling pathway regulates endothelial cell migration, its upregulation explained HUVEC migration-associated gene changes.
As mentioned above, we used HepaRG cells, proliferative hepatoblast-like cells, and HUVECs to generate HOBBs. It is essential for the crosstalk between hepatoblasts and endothelial cells to occur during liver development, as the secretion of cytokines and growth factors is influenced by each other. The upregulated expression of VEGFR2, vWF, CD31, LYVE1, HGF, Stab2, CD32b and VEGFA which participate in hepatoblasts-endothelial cells crosstalk, was observed in HOBB groups (Fig. 3f). Hepatoblasts work as a positive stimulator of endothelial cells morphogenesis and maturation, in which HGF together with VEGF play a decisive role21,34. HGF and VEGF significantly promote expression of sinusoidal endothelial cell markers, such as CD32b, LYVE1, and Stab2, and suppress the expression of angioblast cell markers, such as VEGFR2, VWF, and CD3120,28. We hypothesize that the crosstalk between HepaRG cells and HUVECs was similar to that between hepatoblasts and endothelial cells at liver bud stage, and HOBBs developed dynamically.
Based on qRT-PCR results, HOBB groups showed higher ALB, MRP2, BSEP expression and reduced AFP expression relative to 2D groups (Fig. 3g). Biosynthesis of amino acids, glycolysis-gluconeogenesis signaling pathway enriched by GSEA were upregulated in HOBB groups (Supplementary Fig. 2c). HOBBs had more mature liver function than 2D cells. SOX9, CK19, HEY1, and HES1, landmark genes expressed during bile duct development, were upregulated, suggesting the HOBBs facilitated cholangiocytes differentiation (Fig. 3h). So far, we developed reproducible building blocks with primary hepatobiliary function and multicellular crosstalk for manufacturing large-scale liver tissues.
Reproducible bioprinting of HABs
ECM is complex macromolecular structural network surrounding cells, which is not only a physical scaffold, but also a crucial modulator of biologic processes including cell attachment, migration, differentiation, repair, and development35. Recent studies indicated that ECM modulates hepatic development, regeneration, and even maintenance of the normal architecture and differentiated state36. Collagen Ⅰ is one of the most abundant ECM components supporting hepatocytes in liver cords37. Laminin mainly locates near bile ducts, and is sufficient as a component of basal lamina for the commitment of bipotential liver progenitors to cholangiocytes. Especially, α5-containing laminin is necessary for formation of mature bile duct structures, and laminin 511 improves differentiation of cholangiocytes38. Moreover, collagen I and laminin 511 is coincident with angiogenesis, and commonly used as biomaterials for angiogenesis assay39–41. To promote bile duct morphogenesis and angiogenesis, the liver tissues were built using gelatin-alginate-based bioink (AG) containing collagen I and laminin 511. We added 0.5 mg/ml collagen I and 10 µg/ml laminin 511 to gelatin-alginate-based bioink, developing new quadratic bioink, i.e., AGCL bioink, as shown in Fig. 4a. With AGCL bioink, the 3D-printed constructs had a clear and stable structure with interconnected channels and networks, and the filaments were uniform and smooth (Fig. 4a). SEM image of bioink demonstrated that it was porous, penetrating and scattered with protein fibers, facilitating cell metabolism, migration, and adhesion (Fig. 4b). Integrins constitute a major group of receptors for extracellular matrix components42. Currently collagen-binding I domain-containing integrins are known, namely ITGA2 and ITGB142. Laminin 511 binding integrins are ITGA6 and ITGB141. According to qRT-PCR results, the AGCL groups expressed higher levels of ITGA2, ITGB1, and ITGA6 than AG groups in day 14 (Supplementary Fig. 5), indicating cells were interconnected with collagen Ⅰ and laminin 511 in HABs.
Different from nozzle-free bioprinting technologies (e.g., stereolithography, laser-direct-writing), the microextrusion approach relies on a microneedle to deliver bioinks and print constructions, where the printability of bioink plays a crucial role43. Printability is highly dependent on the rheological properties of bioink, so it is necessary to test its rheological behaviour. Prior to gelation, both AG and AGCL bioink showed a typical fluid-like behaviour (G′′ > G′) (Fig. 4c). Upon cooling, G′ for both solutions increased rapidly and eventually crosses over G′′ showing characteristics of a gel-like structure (Fig. 4c). AGCL bioink's gelation temperature (where G′ and G′′ cross over) was 4°C higher than AG bioink at 27°C, as shown in Fig. 4c. This suggested that printing of AGCL bioink should be conducted at higher temperature to ensure the proper gelation statuses. A print temperature of 27°C had no significant effect on cell activity. Long-term mechanical stability of bioink is required for high-throughput 3D bioprinting. We evaluated the changes in G' and G'' of AGCL bioink through time sweep tests, and it maintained stable viscoelasticity within 30 minutes. In particular, the G' for AGCL was 22.37 ± 0.35 Pa and G'' was 9.15 ± 0.18 Pa with a smoother curve (Fig. 4d). Another key characteristic of bioink for microextrusion printing is shear-thinning. As shown in Fig. 4e, shear force continuously decreased AGCL's viscosity, and it had good shear thinning properties which can prevent clogging and decreasing shear stress from damaging cells. With appropriate printability, AGCL bioink can be used in microextrusion bioprinting and layer-by-layer fabrication.
Creating reproducible liver tissue models requires not only uniform building blocks, but also stable printing process. HOBBs are prone to settling down during long-time printing, so we developed three-level temperature control system which possess precise temperature control (± 0.1℃) to prevent the settlement. In brief, three parts controlled the crosslinking status of bioink: syringe, needle, and substrate. The bioink in syringe was kept in a gel state to prevent the settlement of HOBBs, and microfilaments extruded from nozzle were kept proper crosslinked to print a clear outline and structure. Then, the printer's platform temperature was controlled to maintain a low temperature for supporting bioprinted constructions. HOBBs within constructions were uniformly distributed, allowing for reproducible liver tissue fabrication.
We printed larger HOBBs and found high cell viability after printing (Fig. 4f). This indicates our printing system (including bioink, printing process, and printing parameters) was well adaptable for HOBBs. During culture process, HOBBs within constructions gradually grew larger, and by day 14, their size was averaging 300µm (Supplementary Fig. 6). Due to excessive size and bioink blockage, cells cannot exchange substances resulting in advanced necrosis (Fig. 4f and Supplementary Fig. 7a). The maximum distance between cells and capillaries supplied nutrients and oxygen is 200 µm in vivo. Taking into account growth pattern of HOBBs in scaffolds, we eventually selected HOBBs with 90-µm initial size, and the bioprinting density is 1.06 ± 0.15×107cells/ml (15000 organoids/ml) (Supplementary Fig. 3, 4). The viability of HOBBs remained high until day 14 without necrosis or apoptosis, satisfying in vitro culture requirement of large-scale liver tissues for long period of time (Fig. 4f).
Differentiation and liver function validation of HABs
HepRG cells, exhibiting differentiated cells merging with cholangiocyte and hepatocyte phenotype under DMSO treatment, are frequently described44. The bipotent cell line HepaRG cells constitutes an efficient surrogate of liver function, yet its differentiated status relies on high concentrations of DMSO and is time-consuming, which may compromise the study of drug metabolism and limit the applicability of hepatic model45. Herein, in this work, only a concentration of 0.5% DMSO was required to induce the differentiation of HepaRG cells within 14 days.
HOBBs in scaffolds expanded gradually based on light microscopic examination at different stages, as shown in Supplementary Fig. 6 and Fig. 5a. Cells started sprouting from HOBBs into scaffolds on day 7, and the sprouting cells formed a network by day 14 (Fig. 5a). The cell density had reached 3.33×107/ml at day14 (Fig. 5b). According to SEM image, there is a clear difference between matrix with and without HOBBs, indicating that the HOBBs remodeled the matrix (Fig. 5c). Different morphologies indicated different cell types as shown in Fig. 5a. To determine types of cells were sprouting and aggregating in scaffolds, we performed ALB/CD31 co-immunostaining. Results showed that ALB+/CD31+ cells constituted the aggregated organoids in HABs, but neither ALB nor CD31 were staining cells sprouted in scaffolds (as shown by arrows in Fig. 5e). PAS staining is a tool to mark carbohydrates, and commonly used to evaluate glycogen storage and function of mature hepatocytes46. As expected, HABs exhibited high levels of glycogen storage, acquired a characteristic hepatocyte-like morphology (Fig. 5f). Cryosectioning immunostains further revealed that HABs expressed higher mature liver markers, including ALB, CYP3A4 and MRP2, in comparison with HOBB groups (Fig. 5d).
Based on KEGG pathway analysis of liver function in HAB versus HOBB groups, we found enriched Hippo signaling pathways, NOD-like receptor signaling pathways, FoxO signaling pathways, and thyroid hormone signaling pathways were significantly upregulated, which are important for maintaining hepatocyte function and essential for determining their fate (Fig. 5i). GSEA revealed that the signaling pathway involved in gluconeogenesis, drug metabolism, fatty acid metabolism, biosynthesis of unsaturated fatty acids, N_glycan biosynthesis and histidine metabolism were enhanced in HAB versus HOBB groups, and the expression of related genes was also upregulated (Supplementary Fig. 7b, c). CD31 and VWF levels were reduced in HAB group, while Stab2 and F8 were upregulated, identifying characteristics of convergence to LSECs (Fig. 5g). Gluconeogenesis, urea synthesis, and fatty acid degradation related genes involved in the tricarboxylic acid cycle (TCA) were enhanced in HAB groups compared to HOBB groups. Clearly, the HABs were in constant dynamic change and displayed relatively mature liver function.
HABs with bile duct morphogenesis
HepaRG cells can differentiate into both hepatocytes and cholangiocytes. As mentioned previously, organoids aggregated within HABs were hepatocytes and HUVECs. HOBBs were composed of HepaRG cells and HUVECs, supporting the hypothesis that the cells spreading within HABs were cholangiocytes. We performed CK19/VE-Cad co-immunostaining and identified the spread cells as CK19+/VE-Cad− cholangiocytes (Fig. 6a). Cryosectioning immunostaining revealed the spread cells within HABs expressed cholangiocytes markers as well, including SOX9, Ep-CAM and Cilia (Fig. 6c). Occasionally, cholangiocytes formed tubes with a central lumen as shown in red box (Fig. 6a). Further, fluorescent bile salt was transported by functional cholangiocytes incubated with cholyl-lysylfluorescein (CLF) (Fig. 6b). RNA-seq results revealed upregulation of ABCC4, ABCB1, PRKACB, PRKACB, ADCY3, ATP1B1, ATP1A1, SLC4A4 genes associated with bile acid salt transport and secretion (Fig. 6f). It was suggested that cholangiocytes within HABs had potential to form functional bile ducts. Biliary system is made up of intrahepatic bile duct (IHBD), extrahepatic bile duct (EHBD) and the gallbladder. HAB groups had noticeably higher expression level in intrahepatic bile duct (IHBD) marker genes than extrahepatic bile duct (EHBD) and gallbladder (Fig. 6d). The cholangiocytes interwoven with the hepatocytes were closer to the IHBD cells. FGF signal pathway, TGF-β signal pathway, Notch signal pathway and Wnt signal pathway were main pathways to regulate cholangiocytes differentiation and bile duct morphogenesis in vivo. In HAB groups, their signature genes CK19, SOX9, HES1, HES4, Notch1, BMP4, Wnt5, HNF1β, β-catenin, Jagged2, Smad4, TGFβ2, TGFβ3 and TGFβR2 were altered (Fig. 6h). SOX9, CK19, Hes1, Hes4 and Notch1 were also upregulated in the HAB groups based on qRT-PCR results (Fig. 6e). These data suggested that appropriate signaling pathways were activated to be involved in biliary differentiation. Altogether, part of HepaRG cells differentiated into functional cholangiocytes successfully.
Bile duct development can be divided into two steps, i.e., cell fate decision and morphogenesis27. As shown in Fig. 6a and c, part of HepaRG cells differentiated into cholangiocytes successfully. Then, cholangiocytes form ductal plates, and reorganize into bile duct tubules, which depend matrix remodeling and cellular self-assembly47. Gene Ontology (GO) enrichment analysis displayed that HAB groups upregulated in cell assembly and matrix organization related signaling pathway, including occluding junction, apical junction assembly, adherens junction assembly, focal adhesion assembly, extracellular matrix organization, and tube morphogenesis and branching related signaling pathways, including epithelial cell migration, regulation of cell morphogenesis, morphogenesis of a branching epithelium, morphogenesis of an epithelial sheet, epithelial tube morphogenesis, branching morphogenesis of an epithelial tube compared to HOBB groups (Fig. 6g). Based on GO enrichment analysis, cholangiocytes in HABs were continually assembling and remodeling extracellular matrix, indicating HABs’ potential for bile duct morphogenesis.
As IHBD develop from bipotent hepatoblasts around the portal vein, angiogenic signaling and crosstalk between them are vital to its development48. Gene Ontology (GO) enrichment analysis showed that HAB groups were positively associated with angiogenesis (Fig. 7b). VEGF/VEGFR-2 signaling pathway links bile duct morphogenesis to angiogenesis49. VEGFA, and other angiogenic signals in VEGF/VEGFR-2 signaling pathway, such as FAK, ERK, are enhanced expressed in HAB groups (Fig. 7c). Ang-1 represents a group of tyrosine kinase receptor ligands of VEGFA, which primarily play a role in developmental vascular remodeling and angiogenesis, also upregulated (Fig. 7c). The potential of HABs to form blood vessels was demonstrated. Enriched PPI network, which consisted of significantly enriched angiogenesis and biliary development signaling pathways, further evidenced crosstalk between bile duct morphogenesis and angiogenesis in HABs (Fig. 7d).