Integrated

Communications Biology volume 6, Article number: 310 (2023) Cite this article

2840 Accesses

1 Citations

70 Altmetric

Metrics details

Non-alcoholic fatty liver disease (NAFLD) afflicts a significant percentage of the population; however, no effective treatments have yet been established because of the unsuitability of in vitro assays and animal experimental models. Here, we present an integrated-gut-liver-on-a-chip (iGLC) platform as an in vitro human model of the gut-liver axis (GLA) by co-culturing human gut and liver cell lines interconnected via microfluidics in a closed circulation loop, for the initiation and progression of NAFLD by treatment with free fatty acids (FFAs) for 1 and 7 days, respectively. Co-cultured Caco-2 gut-mimicking cells and HepG2 hepatocyte-like cells demonstrate the protective effects from apoptosis against FFAs treatment, whereas mono-cultured cells exhibit induced apoptosis. Phenotype and gene expression analyses reveal that the FFAs-treated gut and liver cells accumulated intracellular lipid droplets and show an increase in gene expression associated with a cellular response to copper ions and endoplasmic reticulum stress. As an in vitro human GLA model, the iGLC platform may serve as an alternative to animal experiments for investigating the mechanisms of NAFLD.

Non-alcoholic fatty liver disease (NAFLD) is a common, chronic liver disease that leads to hepatic steatosis, cirrhosis, cancer, and cardiovascular diseases1,2,3,4. NAFLD is expected to afflict 33.5% of the United States’ population over 15 years of age by 20305. Currently, liver transplantation is the only method to cure patients with severe liver diseases, and finding donors that match patients is extremely difficult. There is an urgent need for intervention at different stages of fatty liver disease; however, the disease mechanism is largely unknown due to the complicated processes that take place at multiple layers, known as the multiple-hit theory. For example, fat accumulation, oxidative stress, endoplasmic reticulum (ER) stress, and genetic or epigenetic modifications can occur at the cellular level, while insulin resistance and inflammatory responses can occur in multiple organs depending on the individual and environment6. To identify new treatments for NAFLD, a deep understanding of each of these processes is needed, and this accumulated knowledge then needs to be combined.

In this study, we focused on the gut-liver axis (GLA), which is one of the most crucial components for the initiation and progression of NAFLD7,8. The gut is strongly influenced by gut microbiota and dietary carbohydrates, which may accelerate NAFLD9,10,11. Inflammatory products, nutrients, and substances absorbed from food and microbiota via the intestinal barrier are carried by venous blood to the liver. In addition, the products generated by the hepatocytes are transported to the small intestine. Thus, the gut and liver are intricately linked both physiologically and pathologically. GLA dysfunction, including intestinal dysbiosis, bacterial overgrowth, and alteration of mucosa permeability, caused by NAFLD, are potential therapeutic targets;12,13 but no treatments have been made commercially available to date. This is largely because conventional preclinical animal tests do not accurately represent the problems of the multiple-hit theory, lack accessibility to individual organs in living animals, and have species differences. Therefore, establishing a simplified and robust model to study GLA in NAFLD is crucial to obtaining deeper insights into the mechanisms underlying the discovery of new drugs, treatments, and diagnostic tools.

Organs-on-chips (OOCs), which are also known as microphysiological systems (MPSs), hold significant potential for in vitro preclinical tests14,15,16,17 and disease modeling18. Microfluidic technology is the foundation of OOCs because it enables precise control of liquid flow and the three-dimensional architecture of flow channels. These properties give OOCs the ability to control cellular microenvironments spatiotemporally and functionalize tissue cells. The circulation of the cell culture medium can further assist in modeling multi-organ interactions with paracrine and endocrine signaling. OCCs in combination with advanced cellular assays, such as high-content analyses and the omics approach, provide more in-depth insights into biology in a quantitative and multi-parametric manner than animal experiments19. OOCs have been used to recapitulate GLA in vitro and to demonstrate the role of crosstalk via the GLA in pathological situations, including fatty liver disease20,21 and inflammation22, for in vitro pharmacokinetic studies23. However, OOCs need to be further improved to mimic the GLA, and this requires four key features: a closed circulation loop, accessibility to individual chambers, dynamic flow control, and prevention of molecule absorption. The closed circulation loop is necessary for medium circulation to recapitulate inter-tissue interactions in the GLA. Individual accessibility is necessary to introduce tissue cells into the desired chamber and harvest them after treatment without cross-contamination from other cells. The closed circulation loop and individual accessibility may appear to be contradictory features, but both are necessary to investigate the crosstalk between the gut and liver. Dynamic flow control is crucial for obtaining functional tissue cells in vitro, particularly for the gut24. Some OOCs require the use of additional cell culture inserts (e.g., Transwell) to co-culture two or more types of cells separated by a porous membrane. However, due to the macroscale range of the medium, these inserts cannot be utilized for microfabrication and often lack the advantages of microfluidic technology, such as control over the flow dynamics in the cell culture chamber and the cellular microenvironment. These additional inserts often interfere with the microscopic observation of cells as they increase the working distance and light diffraction by the membrane pores. Polydimethylsiloxane (PDMS) is a widely used material for microfluidic cell culture systems due to its biocompatibility, transparency, and elasticity properties. However, preventing PDMS absorption is necessary, as PDMS causes the absorption of hydrophobic molecules, including metabolites, hormones, drug candidates, fatty acids, lipids, and fluorescent indicators, which may influence cellular phenotypes and assay results. Free fatty acids (FFAs) are a critical factor in NAFLD. Although PDMS-based microfluidic cell culture systems25,26, including our previously reported in vitro GLA model27, solve the aforementioned issues of OOCs for GLA, PDMS-based platforms without any treatments for preventing the absorption of hydrophobic molecules are not applicable for recapitulating NAFLD28.

Here, we present an integrated-gut-liver-on-a-chip (iGLC) platform as a simplified in vitro human model of GLA, which can help obtain deeper insights into the underlying mechanisms of NAFLD for the development of new drugs, treatments, and diagnostic tools. The iGLA platform29 has microvalves and a pump to achieve individual accessibility for each cell culture chamber and a closed medium circulation flow to interconnect the gut and liver cells. The integrated micropump controls the perfusion flow to activate cultured gut cells. The iGLC platform does not require additional cell culture inserts and therefore allows high-quality cell monitoring to achieve microscopic single-cell profiling. A simple surface coating with amphipathic molecules prevents the absorption of FFAs into PDMS. We co-cultured gut and liver cells with a closed circulation flow to demonstrate the viability of the iGLC platform as an in vitro human GLA model. We also induced an NAFLD-like cellular state by administering FFAs into the platform for two durations (1 and 7 days) to represent the initial and progressive NAFLD. Finally, we investigated the unique cellular phenotypic changes and associated gene networks for the GLA in the NAFLD-like cellular state by microscopic single-cell profiling in combination with mRNA sequencing (mRNA-seq).

A conceptual illustration of the iGLC platform is shown in Fig. 1a. The design of the microfluidic platform is critical for recapitulating the GLA in vitro. Models based on animal experiments have difficulty investigating the mechanism of NAFLD progression because living organs cannot be connected and disconnected to make simplified multiple-hit-theory models. To elucidate inter-tissue interactions, a single platform is needed that can both mono- and co-culture two or more cell types within the same format. In parallel, cross-contamination needs to be considered to allow the analysis of each cell type and the influence from other cell types. We designed the iGLC platform to meet these requirements (Fig. 1b, c)29. The iGLC platform is made of PDMS with gas permeability, biocompatibility, and transparency and consists of perfusion and control layers. The perfusion layer had three sets of two cell culture chambers (2.1 mm width and 220 µm height) linked by microfluidic channels (200 µm width and 45 µm height). The control layer had a thin and flexible PDMS membrane (200 × 200 µm2 area and 20 µm thickness) to integrate the microvalves and pump, which allowed robust open/close control with our previously reported microfabrication technology (Supplementary Fig. S1)29,30. To investigate the interactions between two or more types of tissue cells, it is necessary to isolate and connect them as desired. Because most OOCs do not have integrated valves and pumps, they require assembling and disassembling tubes among multiple devices, which reduces the accuracy of the control over their interactions, but avoids undesired cross-contamination. The integrated micropump provides more accurate control over the perfused flow than an external pump within the microfluidic device. In addition, it reduces the sample loss caused by extra tubing. The PDMS membranes at the microvalves and pump were actuated by a positive hydraulic pressure (150 kPa) applied by the computer-controlled solenoid valves. The microvalves allow individual cell culture chambers to be accessed without cross-contamination and precise control over sample or reagent introduction. The integrated micropump provides closed-loop medium circulation to recapitulate inter-tissue interactions, such as the GLA, and to regulate the flow rate within 0–20 nL min−1 by tuning its actuating cycle (Supplementary Fig. S2).

a Schematic for NAFLD progression by FFAs via the GLA. b Photograph of an iGLC platform. The perfusion and control layers are colored in pink and blue, respectively. c Illustration of the iGLC platform used for NAFLD. Two cell culture chambers are used for gut cells (Caco-2) and hepatocytes (HepG2) and are linked via a microfluidic channel with a micropump to achieve closed medium circulation of FFAs. Within the microvalves, each cell culture chamber is individually accessible without the risk of cross-contamination. Thus, this setup allows the inter-tissue interaction to be evaluated. A–A’ shows the cross-section of the cell culture chambers for the Caco-2 gut and HepG2 liver cells. B–B’ shows the cross-sectional view of the open and closed integrated microvalves. Because of the elastic PDMS membrane, the normally open valve is closed by applying higher pressure to the microfluidic channel in the control layer.

Prior to the cell culture and FFAs treatments, we coated the PDMS surface of the cell culture chambers with n-dodecyl β-D-maltoside (DDM)31,32 and then Matrigel to reduce the absorption of hydrophobic molecules (e.g., FFAs and AdipoRed fluorescent lipid marker)33,34 and to increase cell adhesion and growth (Supplementary Fig. S3). This is a critical problem when using PDMS-based OOCs for drug discovery and in vitro disease modeling because small-molecule drug candidates and lipids are absorbed by the PDMS before reaching the target cells, and the concentrations differ from those in a microtiter plate. Under these conditions, the results obtained would be unreliable. DDM is a suitable amphipathic molecule that prevents PDMS absorption because its hydrophobic part binds to the PDMS surface and the hydrophilic part prevents the absorption of hydrophobic molecules. However, DDM is not capable of promoting cell adhesion and growth on the PDMS surface. Therefore, we evaluated combining DDM with a subsequent layer of Matrigel. The DDM/Matrigel coating partially prevented the absorption of the AdipoRed fluorescent lipid marker into the PDMS (Supplementary Fig. S4). We also measured the remaining FFAs [a mixture of palmitic acid (PA) and oleic acid (OA) with a molar ratio of 1:2; see Methods] in the cell culture medium after incubation with the coated PDMS (Supplementary Fig. S5). PA was not absorbed by either the non-coated or DDM/Matrigel-coated PDMS after 6 h of incubation at 37 °C, whereas it showed absorptions of 53.7% and 40.9%, respectively, after 24 h of incubation. OA showed absorptions of 22.0% and 17.4%, respectively, after 6 h of incubation, and absorptions of 67.4% and 56.7%, respectively, after 24 h of absorption. These results suggest that although the DDM/Matrigel coating slightly mitigated the absorption of both PA and OA into PDMS after 24 h of incubation, majority of the FFAs were still absorbed. Therefore, the medium was replaced every 6 h during the experiments. However, further improvements are required to prevent PDMS absorption, and other structural materials should be used for OOCs. This has long been a topic of discussion in this research field17,35.

To recapitulate GLA in vitro, Caco-2 and HepG2 cells were introduced into individual cell culture chambers in the iGLC platform without cross-contamination, but with a closed circulation flow (Fig. 2a). To confirm that the iGLC platform allows sustainable cell cultivation, Caco-2 and HepG2 cells were cultured with a medium circulation flow of 15 nL min−1 for 7 days and evaluated using a Calcein AM fluorescent cell viability indicator (Fig. 2b–d and Supplementary Data 1). Although Caco-2 cells generally take approximately 21 days to obtain functionality, perfusion conditions allow Caco-2 cells to express their functions within a shorter time, such as 7 days36,37. Thus, we chose a 7-day cultivation time to recapitulate GLA in vitro. For comparison, we also assessed separately cultured Caco-2 and HepG2 cells in cell culture chambers at a flow rate of 15 nL min−1 under static conditions. Based on computational fluid dynamics (CFD), 3.53 × 10−4 dyne cm−2 of fluid shear stress (FSS) was generated on the cells (Supplementary Fig. S6). The flow condition did not affect the viability of HepG2 cells, which showed minor improvement with a circulating flow. The Caco-2 cells showed increased viability with a circulating flow compared to that with static conditions. The fluorescence intensity of Calcein AM increased approximately 7.4-fold with the circulating flow. Caco-2 cells have previously been reported to show improved functionality and viability when cultured under continuous perfusion on a microfluidic device within 1-week of culture36,37. This proves that the iGLC platform provides better circulation flow conditions for both Caco-2 and HepG2 cells.

a, b Experimental procedure to culture Caco-2 and HepG2 cells in an iGLC platform (see also Supplementary Methods). Briefly, after a microfluidic channel was washed with a fresh cell culture medium, all valves were closed to prevent air contamination into the chip. Two valves next to the cell culture chamber were opened, and cell suspensions of Caco-2 and HepG2 cells were separately introduced into the corresponding chambers. After 1 day of incubation at 37 °C to settle down the cells, the pump was actuated to circulate the medium. The cell culture medium was changed every 6 h. c Phase-contrast (PC) and fluorescent micrographs of Caco-2 and HepG2 cells with closed medium circulation in the iGLC platform and stained with Calcein AM (C_AM) and Hoechst 33258 (Hst). For comparison, we also assessed individually cultured Caco-2 (C:C) and HepG2 (H:H) cells under circulation flow (CF) and static conditions (S). H:C represents co-cultured Caco-2 and HepG2 cells. The scale bar represents 100 µm. d Ridgeline plots showing single-cell profiling of Caco-2 and HepG2 cells stained with Calcein AM. For comparison, we also assessed individually cultured Caco-2 (C:C) and HepG2 (H:H) cells under circulation flow (CF) and static conditions (S). H:C represents co-cultured Caco-2 and HepG2 cells. The p-values were estimated with the Tukey–Kramer test and are presented in Supplementary Table S1 for the Caco-2 and Supplementary Table S2 for HepG2 cells.

To induce NAFLD in the iGLC platform, we used FFAs represented by a mixture of PA and OA, which is typical of Western diets (see Methods)38. As shown in Fig. 2a, prior to the FFAs treatment, Caco-2 and HepG2 cells were cultured in the iGLC platform with Dulbecco's modified Eagle's medium (DMEM) supplemented with 10%(v/v) FBS, and then, replaced to serum-free DMEM for 12 h to ensure cell starvation. A series of FFAs concentrations from 0 to 2 mM were introduced into the iGLC platform, and the cells were incubated for 24 h with medium circulation. To evaluate the FFAs accumulation in the mono-cultured Caco-2 and HepG2 cells, the cells were stained with AdipoRed lipid fluorescent dye to visualize the intracellular lipids (Supplementary Fig. S7 and Supplementary Tables S3 and S4). Both Caco-2 and HepG2 cells showed dose-dependent intracellular FFAs accumulation, and 1 mM FFAs treatment was used for further studies. Lipid droplet accumulation in hepatocytes is a hallmark of NAFLD39 and was observed in our model. During FFAs treatments, we used serum-free medium supplemented with 1% (w/v) BSA to minimize the effects of cell proliferation. Both mono-cultured Caco-2 and HepG2 cells showed a reduction in cell proliferation and did not interfere with cell treatment because of over-cell growth (Supplementary Fig. S8). We then examined two treatment periods, 1 and 7 days, representing NAFLD initiation and progression, respectively. We also confirmed that the FFAs treatment for 7 days increased intracellular lipid accumulation in mono- and co-cultured Caco-2 and HepG2 cells (Fig. 3a, b and Supplementary Data 2). The co-cultured Caco-2 cells showed less lipid accumulation than the mono-cultured cells, but HepG2 cells showed increased lipid accumulation. We also performed Annexin-V staining to evaluate apoptotic status after FFAs treatment. For comparison, the cells were also treated with 1 µM staurosporine (STS), which induces apoptosis. Although PA accumulation has been reported to cause cytotoxicity after 1-day of treatment40,41, the apoptotic staining (Fig. 4a) followed by quantitative single-cell profiling (Fig. 4b and Supplementary Data 3) results suggested that the 7-day FFAs treatment induced apoptosis in mono-cultured Caco-2 and HepG2 cells, whereas 1-day treatment did not induce apoptosis (Supplementary Fig. S9). However, the co-culture of Caco-2 and HepG2 cells showed a reduction in apoptosis. Notably, the expression of albumin (ALB), a functional hepatocyte marker, in mono-cultured HepG2 cells was not changed by FFAs treatment, but co-cultured HepG2 cells demonstrated increased ALB expression (Fig. 4c, d and Supplementary Data 3). This suggested that HepG2 cells co-cultured with Caco-2 cells displayed improved functionality.

a Phase contrast and fluorescent micrographs of Caco-2 and HepG2 cells treated with FFAs (0 and 1 mM) for 7 days stained with AdipoRed lipid fluorescent dye. The scale bars represent 100 µm. b Ridgeline plots to evaluate FFAs accumulation in individual cells [Caco-2 (left) and HepG2 (right)] after FFAs treatment for 7 days. The p-values were estimated with the Tukey–Kramer test and are presented in Supplementary Tables S5, S6 for the Caco-2 and HepG2 cells, respectively. Cells with over 20 of AdipoRed fluorescent intensity were considered as lipid accumulated cells. The data represent the mean ± standard deviation (n = 3).

a Phase contrast and fluorescent micrographs of Caco-2 and HepG2 cells treated with FFAs (0 and 1 mM) stained with the Annexin V apoptotic cell marker. The scale bars represent 100 µm. b Ridgeline plots to evaluate individual apoptotic cells (Caco-2 and HepG2) after FFAs treatment for 7 days. For comparison, cells were treated with 1 µM of staurosporine (STS) for 24 h. The p-values were estimated with the Tukey–Kramer test and are presented in Supplementary Tables S7, S8 for the Caco-2 and HepG2 cells, respectively. Cells with over 0.3 of Annexin V fluorescent intensity were considered as apoptotic cells. c Phase contrast and fluorescent micrographs of HepG2 cells treated with FFAs (0 and 1 mM) stained with the human albumin. The scale bars represent 100 µm. d Ridgeline plots to evaluate individual HepG2 cells with human albumin expression after FFAs treatment for 7 days. The p-values were estimated with the Tukey–Kramer test and are presented in Supplementary Table S9 for the HepG2 cells. Cells with over 100 of ALB fluorescent intensity were considered as ALB-expressing cells. The data represent the mean ± standard deviation (n = 3).

To further evaluate the cellular phenotypes after the FFAs treatment in detail, we stained cells with Calcein AM and Hoechst 33258 fluorescent dyes. This was followed by high-content analysis (HCA) using the CellProfiler computer-guided single-cell analysis software (Fig. 5, Supplementary Data 4, Supplementary Figs. S10, S11)42,43. In general, the Calcein AM and Hoechst 33258 dyes are used for cell viability function and as nuclear markers. Calcein AM staining can also be applied to evaluate cell and nucleus morphologies to obtain cellular parameters for HCA44,45. Moreover, recent advances in single-cell profiling based on high-content images, such as t-distributed stochastic neighbor embedding (t-SNE), have allowed a quantitative understanding of individual cellular status with high-dimensional cellular parameters and visualizing such parameters in a two-dimensional map46. We analyzed microscopic images to quantify 68 cellular parameters (Supplementary Data 4 and Supplementary Table S10) for each cell of Caco-2 and HepG2 cells treated with 1 mM FFAs for 7 days (Fig. 5a, b and Supplementary Data 4). We identified specific parameters that could be distinguished between non-treated and FFAs-treated cells that are associated with the cellular shape (AreaShape_FormFactor, AreaShape_MeanRadius, and AreaShape_EquivalentDiameter), Calcein AM intensity (Intensity_LowerQuartileIntensity, Intensity_MeanIntensity, and Intensity_MaxIntensity), nucleus shape (Nucleus_AreaShape_FormFactor), and Hoechst 33258 nuclei intensity (Intensity_LowerQuaterileIntensity, Intensity_MeanIntensity, Intensity_MaxIntensity, and Intensity_MinIntensityEdge). Notably, for the Caco-2 cells, the most distinguishing features were the cellular shape (AreaShape_EquivalentDiameter), Hoechst intensity (Intensity_MaxIntensity), and Calcein AM intensity (Intensity_MaxIntensity) (Fig. 5a). However, for the HepG2 cells, the most distinguishing feature was the Calcein AM intensity (Intensity_MeanIntensity and Intensity_MaxIntensity), the FFAs-treated HepG2 cells showed a reduction in these features (Fig. 5b). Thus, HCA based on simple cell and nucleus staining can be used to identify minute cellular phenotypic changes upon FFAs treatment that cannot be distinguished using molecular apoptotic markers. Because the iGLC platform does not require the use of cell culture inserts, it allows HCA-based single-cell profiling.

a, b Two-dimensional t-SNE plots of microscopic single-cell profiling of co-cultured Caco-2 (a) and HepG2 (b) treated with 1 mM of FFAs or no treatment and stained with Calcein AM cellular and Hoechst 33258 nuclei markers (n = 3). The most distinguishable cellular parameters (AreaShape_EquivalentDiameter, Intensity_MaxIntensity_Calcein AM and Intensity_MaxIntensity_Hoechst 33258) are shown in the corresponding t-SNE plots, as well as the boxplots in Supplementary Fig. S11. The p-values are presented in Supplementary Tables S11,S12.

To elucidate the effects of FFAs treatments and crosstalk via co-culturing on gene expression, we performed mRNA-seq, followed by PCA on the differentially expressed genes (DEGs), which were defined by an analysis of variance (ANOVA) (Fig. 6a, b, Supplementary Figs. S12, S13, Supplementary Data 5, 6 and Methods). We identified 654 and 1330 genes in Caco-2 and HepG2 cells, respectively, that were differentially expressed among the mono- and co-cultured conditions with or without 1 mM FFAs treatments for 7 days. In Caco-2 cells, the PC1 axis distinguished mono- and co-culture conditions, whereas the PC2 axis distinguished FFAs treatment. Notably, the co-culture conditions showed less of an effect on gene expression altered by FFAs treatments. In contrast, in the case of HepG2 cells, the PC1 axis clearly distinguished HepG2 cells treated with 1 mM FFAs under co-culture conditions from those treated with the other conditions. The PC2 axis in PCA of HepG2 cells also showed changes in gene expression associated with FFAs treatment; however, co-culturing reduced the changes in gene expression following FFAs treatment.

a, b PCA results for DEGs obtained from the Caco-2 (a) and HepG2 (b) experimental sets. The mean PC2 values among the replicates are shown with dotted lines. For these diagrams, the gene name and change direction shown in (a, b) were considered. c, d Heat maps for the DEGs. Z-values of the expression profiles are shown.

The heat map of K-means clustering for DEGs in Caco-2 cells showed multiple clusters with similar upregulated gene expression patterns under co-culture conditions, compared to mono-culture conditions (Clusters 1, 3, 5, and 10 in Fig. 6c). These clusters were related to cell cycle and metabolic processes, as suggested by gene ontology (GO) terms (Supplementary Figs. S14, S15). We also observed downregulated genes specific to the co-culture condition in Caco-2 cells (Clusters 2 and 6); however, the number of genes were too small to obtain sufficient interpretations from GO analysis. In the case of HepG2 cells, Cluster 18 in the heat map for DEGs showed unique upregulated genes in the co-culture condition, whereas only Cluster 16 showed downregulated genes specific to the co-culture condition (Fig. 6d, and Supplementary Figs. S16, S17).

We then compared the effects of FFAs treatment under mono- and coculture conditions on gene expression in Caco-2 and HepG2 cells. The mRNA-seq results revealed that Caco-2 cells showed elevated expression of genes associated with mitosis (Cluster 8 in Figs. 6a, 7a), whereas stress-responsive genes (Cluster 7 in Figs. 6a, 7a) were downregulated. In contrast, HepG2 cells treated with FFAs displayed elevated expression of genes related to the cell cycle (Clusters 12 and 17 in Fig. 7b)47, and suppression of the expression of genes associated with cell-cell adhesion and nucleic acid transport (Cluster 15 in Fig. 7b).

a, b Bar graphs showing gene enrichment related to certain GO terms and pathways for representative K-means clusters of FFAs-treated Caco-2 and HepG2 cells under mono- and co-cultured conditions. The other clusters were shown in Supplementary Figs. S14–S17.

We established a human iGLC platform to investigate GLA response to FFAs in an NAFLD model. Although OOCs have been used to study NAFLD, the iGLC platform has multiple advantages owing to its application of microfabrication technology. Integrating the microvalves and a pump in the platform enables accessibility to individual cell culture chambers without undesired cross-contamination and closed medium circulation to mimic human GLA.

To investigate the effects of FSS on cells caused by closed-circulating flow, we conducted CFD to estimate the FSS generated by closed-circulation flow. The result was 3.53 × 10−4 dyne cm−2 (Supplementary Fig. S6). Under these flow conditions, Caco-2 cells showed increased cell viability (Fig. 2c), while HepG2 cells did not show any change (Fig. 2d). A number of previous reports have suggested that Caco-2 cell functionality can be improved by applying perfusion flow at over 0.0002 dyne cm−2, which is similar to the FSS level we tested48,49,50,51. Additionally, functional activation by FSS in HepG2 cells has been reported, but it requires a higher FSS of over 0.1 dyne cm−2 52,53.

Notably, co-cultured HepG2 and Caco-2 cells significantly changed their gene expression profiles compared to mono-cultured cells (Clusters 1, 3, 5, 10, and 18 in Fig. 6c, d, and Supplementary Figs. S14–S17). The upregulated genes in Caco-2 cells were related to Mitotic G1 phase and G1/S transition, respiratory electron transport, ATP synthesis by chemiosmotic coupling, heat production by uncoupling proteins, positive regulation of spindle checkpoint, and small-molecule catabolic processes. Most of them were associated with the cell cycle and metabolic processes, but the cell proliferation analysis did not show a significant difference between co- and mono-cultured Caco-2 cells (Supplementary Fig. S8). Thus, gene expression in Caco-2 cells co-cultured with HepG2 cells did not affect the cell cycle and proliferation. In contrast, HepG2 cells upregulated genes showed GO terms of ‘plasma membrane-bound cell projection assembly (Cluster 18 in Supplementary Fig. S17), and downregulated genes showed ‘microtubule-based movement’ (Cluster 16 in Supplementary Fig. S17). Similarly, we did not observe any obvious cellular phenotypes associated with changes in gene expression.

To recapitulate NAFLD-like GLA in vitro, FFAs were administered to the cells in the iGLC platform (Fig. 2a). Although lipid accumulation occurred in both Caco-2 and HepG2 cells treated with FFAs for 7 days (Fig. 3a, b), the numbers of apoptotic Caco-2 and HepG2 cells in the co-culture condition were significantly reduced compared to those in mono-culture conditions.

In the case of Caco-2 cells, FFAs treatment for 1 day did not affect the apoptotic cellular status in either condition (Supplementary Fig. S9). Enriched GO terms are presented in Supplementary Fig. S13 is associated with cellular response to copper ions. Copper is an essential trace element for human physiological processes, including enzymatic metabolism54. Dysregulation of intracellular copper ions generates reactive oxygen species (ROS) and metabolic imbalance, which results in DNA damage and apoptosis55,56,57. Our results suggest that Caco-2 cells treated with FFAs for 1 day did not show an apoptotic cellular state with cellular membrane disruption, but initiated DNA damage due to the cellular response to copper ions. In contrast, even though the prolonged 7-day FFAs-treated mono-cultured Caco-2 cells showed an increase in the count of apoptotic cells (22.9 ± 9.9%), the 7-day FFAs-treated Caco-2 cells co-cultured with HepG2 cells had protective effects against lipid-induced apoptosis (8.2 ± 7.2%), and similar levels were observed in non-treated Caco-2 cells (mono-cultured, 7.8 ± 2.4%; co-cultured, 10.3 ± 2.2%). Previously, although intestinal cells treated with FFAs showed apoptosis due to lipid accumulation followed by an increase in the levels of reactive oxygen species58,59, our microscopic single-cell profiling suggests suppression of the apoptotic signaling, which is in agreement with mRNA-seq results shown in Figs. 6a, c and 7a, indicating the down-regulation of ‘Signaling by Rho GTPase,’ the typical apoptosis-associated signaling pathway60.

In addition to examining the apoptotic status of Caco-2 cells, we also further examined the correlation between phenotypic changes in cellular morphology and gene regulation networks for in vitro GLA in the NAFLD-like state by single-cell profiling (Fig. 5 and Supplementary Fig. S10). At the initiation stage of 1-day FFAs treatment, Hoechst 33258 intensity was slightly higher in non-treated Caco-2 cells than in FFAs-treated Caco-2 cells (Supplementary Figs. S10), but at the progression stage by 7-day FFAs treatment, the Hoechst 33258 intensity showed the opposite result, where FFAs-treated Caco-2 cells showed a higher intensity than the non-treated Caco-2 cells (Fig. 5). In general, cellular uptake of Hoechst 33258 can occur either by ABC transporters or by diffusion via a damaged cellular membrane. The mRNA-seq results for Caco-2 cells treated with FFAs for 1 or 7 days (Supplementary Figs. S13–S15) did not indicate any enriched GO terms associated with transporters (i.e., ATP-binding cassette transporter ABCG2)61 that may efflux Hoechst 33258 outside the cells. Thus, the transporters did not affect the Hoechst 33258 intensity. Therefore, FFAs-treatment might damage the cellular membrane, resulting in an increase in Hoechst 33258 intensity, but not induction of apoptosis.

In contrast, FFAs-treated HepG2 cells had the most distinguishable features, such as Calcein AM intensity (Intensity_MeanIntensity and Intensity_MaxIntensity), and FFAs-treated HepG2 cells showed a reduction in these features (Fig. 5b). Considering that the intensity of Calcein AM represents cell viability62, this result suggests that FFAs reduce the viability of co-cultured HepG2 cells. In contrast, we did not observe an increase in apoptotic HepG2 cells under co-culture conditions, as confirmed by Annexin V staining (Fig. 4a, b). Notably, the mRNA-seq results revealed that HepG2 cells co-cultured with Caco-2 cells showed increased gene expression associated with the cell cycle (Clusters 12 and 17 in Fig. 7b). In addition, the expression of genes associated with cell-cell adhesion was reduced in FFAs-treated HepG2 cells under co-culture conditions (Cluster 15 in Fig. 7b). These results suggested that, as FFAs treatment induced mild damage in HepG2 cells co-cultured with Caco-2 cells, HepG2 cells transformed from differentiated liver cells to proliferating cells. Cell cycles and cell division are essential to maintain tissue homeostasis, but in diseases such as NAFLD, they often highlight the imbalance of the cell cycle caused by disrupted metabolic signaling47. Overall, the results with our platform indicate that HepG2 cells co-cultured with Caco-2 cells were activated towards NAFLD-like gene expression networks with FFAs treatment.

To develop an advanced in vitro NAFLD model, the cell types used are critical. In this study, HepG2 human hepatocellular carcinoma cell line and Caco-2 human colorectal adenocarcinoma cell line were used since they have been immortalized and widely used until now. However, such cell lines have genetic mutations and cannot express proper functions as an actual hepatocyte63 and intestinal cells64. In the case of hepatocytes, the HepaRG human hepatocellular carcinoma cell line has been used for liver disease modeling as well as toxicological research due to their relatively higher functionalities than those of HepG2 cells recently, but HepaRG cells also showed differences in lipid metabolism65. Moreover, cell lines generally form a pure cellular population, and cannot represent the cellular diversity. For example, the liver contains hepatocytes and non-parenchymal cells (NPCs) such as Kupffer cells, endothelial cells, and stellate cells. In particular, since the liver microphysiological system (L-MPSs) was reported to recapitulate liver steatosis using NPCs66, NPCs are now recognized for their importance in hepatic disease modeling. In contrast, primary cells obtained from NAFLD patients showed similar cellular diversity and drug metabolic responses67. However, primary cells have limited cell growth capability and often lose their functions in vitro. Moreover, it is difficult to identify healthy donors to obtain a control sample. In addition to primary cells, human pluripotent stem cells (hPSCs), such as embryonic68 and induced pluripotent stem cells69,70,71, can provide multiple cell types from a single-cell source after induction of differentiation by applying growth factors and nutrients. The recent advancement of gene editing technologies, such as CRISPR/Cas systems72, allows the creation of gene-modified hPSCs used for disease modeling, for which the primary cells are not applicable. However, the tissue cells derived from hPSCs could not fully cover the entire cellular diversity nor express their functionality in vivo. Therefore, it remains challenging for hPSCs to fully recapitulate the pathological conditions of NAFLD in vitro. In this study, although hepatic and intestinal cell lines were used to recapitulate GLA and NALFD, they were able to provide new insights into NAFLD-like gene expression networks upon FFAs treatment. Using primary cells or hPSC-derived hepatocytes/intestinal cells with higher functionality and cellular diversity, the iGLC platform will allow the reproduction of the ideal results of in vivo NAFLD initiation and progression.

Moreover, as demonstrated in this study, NAFLD is involved in organ-to-organ interactions. Recent reports have suggested that the gut microbiota significantly impacts NAFLD development and needs to be considered for establishing a new NAFLD model9,10,11. However, each organ often requires optimal culture conditions to express its proper functions. Our current iGLC platform can only use the same cell-culture medium to interconnect two different tissues; therefore, the next generation of iGLC platforms needs to overcome these issues. Recently, we reported the multilayered GLA-MPS27,73 to realize both individual accessibility for each tissue and closed loop medium circulation, for recapitulating inflammatory bowel disease by treating intestinal cells with an inflammatory inducer (e.g., lipopolysaccharide). However, the multilayered GLA-MPS device was made of PDMS and was not suitable for modeling NAFLD without any surface treatment against molecular absorption, such as FFAs28. As shown in Supplementary Fig. S5, because the surface coating for PDMS reduced FFAs absorption, and the multilayered GLA-MPS with the PDMS surface coating will be able to serve as a platform to study NAFLD initiation and progression by treating cells with FFAs in a more quantitative and robust fashion.

Instead, it is also valuable to find alternative materials to PDMS to avoid the well-known concerns of PDMS74,75. Thermoplastic materials (e.g., polystyrene, polymethyl methacrylate, polypropylene, and cyclo-olefin polymers) could be applicable for fabricating OOCs and mass production, but it is difficult to deform membranes for actuating microvalves/pumps or stretching cultured cells76,77, owing to their rigidity. Perfluoropolyether elastomers have recently been used for fabricating OOCs enabling the prevention of the absorption of hydrophobic molecules78,79. Although perfluoropolyether elastomers could be attractive for modeling NAFLD in vitro, such elastic materials still pose a challenge for mass production for industrial applications. Therefore, it is critical to choose the optimal materials for fabricating OOCs to fulfil the requirements of the targeted studies and applications to obtain more accurate results representing human pathophysiological conditions80.

In summary, the iGLC platform represents a new in vitro human model for recapitulating physiological and NAFLD pathological conditions, with a focus on GLA. In the near future, the iGLC platform combined with HCA and omics approach may provide deeper insights into NAFLD development for the establishment of new drugs. The iGLC platform can contribute to the development of drugs for NAFLD, as well as a variety of disorders associated with GLA81, such as inflammatory bowel disease, which do not have any in vitro experimental settings.

The iGLC platform was fabricated from a flexible PDMS (SYLGARD 184, Dow Corning) polymer using a multilayer soft lithography replica molding technique, as previously reported (Supplementary Fig. S1)82,83. Briefly, the control layer consisted of microchannels supplying hydraulic pressure that were cast from a 30 µm thickness of negative photoresist mold (TMMR S2000, Tokyo Ohka Kogyo) and patterned using a standard photolithography tool. The perfusion layer consists of two cell culture chambers (225 µm in height and 2.1 mm in width) interconnected by semi-elliptical microchannels (45 µm in height and 200 µm in width). The mold fabrication process for the perfusion layer followed the multilayer lithography principle, which combines standard photolithography (cell culture chambers) and grayscale lithography (microchannels). First, a negative resist layer (TMMF S2045, Tokyo Ohka Kogyo) with a thickness of 180 µm was patterned on a silicon wafer. Next, a positive resist (PMER P-LA900PM, Tokyo Ohka Kogyo) was spin-coated with a thickness of 45 µm on the wafer. Subsequently, digital micromirror device (DMD)-based grayscale lithography (DL-1000GS/KCH, NanoSystem Solutions) was performed using numerically optimized mask data82 to achieve precise wafer-level mold fabrication. This allowed complete sealing of the microchannels with microvalves and high-efficiency driving of the peristaltic micropumping system. After the mold fabrication, the PDMS base and curing agent were thoroughly mixed at a weight ratio of 10:1. For the control layer, the PDMS mixture was spin-coated with a thickness of 50 µm on the mold to obtain a controlled PDMS thickness of 20 µm at the membrane portion. The PDMS for the control and perfusion layers was cured on a hotplate at 80 °C for 4 min and in a convection oven at 80 °C for 40 min, respectively. The perfusion layer was then peeled off, precisely aligned, and bonded to the control layer using a partial PDMS curing method. The assembled structure was placed in an oven at 80 °C for 2 h and then peeled from the silicon wafer. Finally, the inlet and outlet wells were opened. The assembled device was permanently bonded using O2 plasma (FA-1, SAMCO) onto a microscopic slide glass (25 mm × 75 mm).

The microvalves and micropump were actuated by positive hydraulic pressure via linked control channels. The control channel within the chip was first filled with distilled water using a 1 mL syringe to prevent gas permeation across the PDMS. Metal pins and Teflon tubes (Pilot Corporation) were then used to connect the inlet wells of the control channels and pneumatic system with a compressed nitrogen gas resource (regulated at 0–200 kPa). The pressure actuation and release of the valves were controlled and operated with LabVIEW software (Version 11.0, National Instruments) via solenoid valves (Microfluidic System Works Inc. and LEE Company) using a controller board (VC3 8 controller [ALA Scientific Instruments] and NI USB-6501 [National Instruments]). The micropump consists of three adjacent microvalves with sequential actuation to provide periodic peristaltic motion to generate a medium recirculation flow in the chip.

HepG2 human hepatocellular carcinoma and Caco-2 human colorectal adenocarcinoma cell lines were obtained from the American Type Culture Collection. Cells were maintained in Dulbecco's modified Eagle's medium (DMEM) (Sigma-Aldrich, St. Louis, MO, USA) supplemented with 10% (v/v) fetal bovine serum (FBS, Cell Culture Bioscience), 1% (v/v) nonessential amino acids (Thermo Fisher Scientific), and 1% (v/v) penicillin/streptomycin (Thermo Fisher Scientific) in a humidified incubator at 37 °C with 5% (v/v) CO2. HepG2 and Caco-2 cells were passaged using trypsin/EDTA (0.04% / 0.03%[v/v]) solution every 3 and 5 days and at ratios of 1:5 and 1:10, respectively.

To ensure the absence of mycoplasma contamination in cellular samples, a mycoplasma detection assay was conducted by using the MycoAlert™ mycoplasma detection kit (Lonza, LT07-118) according to the manufacturer's instructions. The positive and negative controls were included in each assay run to ensure the accuracy of the results. By measuring the luminescence before (read A) and after the addition of the MycoAlert™ Substrate (read B), a ratio was used to identify the presence or absence of mycoplasma. The HepG2 and Caco-2 cells have 0.7 and 0.3 B/A ratio, respectively. Then, compared with negative control 0.4 and positive control 18.6, the cellular samples (<0.9) were confirmed no mycoplasma contamination.

Before cell seeding, the platform was sterilized by washing with 70% ethanol and placed under ultraviolet light in a biosafety cabinet for 30 min. Subsequently, the cell culture chambers were coated with 0.1% (w/v) DDM (n-dodecyl β-D-maltoside) in PBS at 4 °C for 24 h and then coated with Matrigel hESC-qualified matrix (Corning) diluted to 1.3% (v/v) in DMEM/F12 (Sigma-Aldrich) at 4 °C for 24 h (Supplementary Fig. S3). After the excess Matrigel was rinsed with DMEM, the chip was placed in an incubator at 37 °C until use.

HepG2 and Caco-2 cells were harvested from the culture flasks with 1 mL of trypsin/EDTA (0.04%/0.03% [v/v]) solution and incubated at 37 °C for 5 min. After centrifugation, the cells were resuspended in fresh cell culture medium at 1.0 × 106 cells mL−1. The microvalves next to the cell culture chambers were closed to prevent cross-contamination during cell seeding. Subsequently, 5 µL of the Caco-2 cell suspension was introduced with a pipette into the well adjacent to the cell culture chamber at 7.0 × 104 cells cm−2. HepG2 cell suspension (5 µL) was introduced into another cell culture chamber. The platform was placed in a humidified incubator at 37 °C and 5% (v/v) CO2. After one day, the cell culture medium was changed to remove floating dead cells. The cell culture medium was changed every 6 h under the control of our custom LabVIEW-based software.

The FFAs treatment solutions were a mixture of PA (Sigma-Aldrich) and OA (Sigma-Aldrich) at a molar ratio of 1:2. To prepare the treatment solutions, PA was dissolved in dimethyl sulfoxide (DMSO; Nacalai Tesque Inc.) at 20 mg mL−1 to serve as the PA stock solution (78 mM). OA was dissolved in DMSO at 100 mg mL−1 to serve as the OA stock solution (354 mM). PA and OA stock solutions were mixed (PA:OA = 1:2) in DMEM (containing 1% BSA fatty acid-free, 1% P/S, and 1 mM nonessential amino acids) to generate a series of FFAs concentrations (0.1, 0.2, 0.5, 1, and 2 mM). Prior to FFAs treatment, both the cell culture chambers were replaced with serum-free DMEM for 12 h of cell starvation. FFAs-containing medium was then introduced into the cell culture chamber. The platform was incubated at 37 °C in a humidified incubator for 7 days, and the medium was exchanged every 6 h.

Calcein AM (Dojindo Molecular Technologies, Inc.) and Annexin V-Alexa Fluor® 647 (BioLegend) dyes were used to stain viable and apoptotic cells, respectively. Hoechst 33258 (Dojindo Molecular Technologies, Inc.) was used to stain the nuclei. The staining solution was prepared by mixing 10 µL Hoechst 33258 (1 mg mL−1 stock concentration), 10 µL Calcein AM (1 mg mL−1 stock concentration), 10 µL Annexin V (50 µg mL−1 stock concentration), 500 µL Annexin V binding buffer (BioLegend), and 500 µL DMEM. The wash buffer comprised 500 µL of Annexin V binding buffer and 500 µL of DMEM. After treatment, the cells were washed twice with fresh DMEM, and 10 µL of staining solution was introduced into a cell culture chamber using a pipette via the adjacent inlet. The cells were then incubated at 37 °C for 30 min. Excess staining solution was removed by washing with 30 µL of wash buffer three times.

Lipid accumulation was visualized using the AdipoRed assay (Lonza), following the manufacturer's protocol. Briefly, 15 µL of AdipoRed assay reagent and 10 µL of Hoechst 33258 were mixed with 1 mL of DMEM to prepare AdipoRed staining solution. The cell culture chambers and channels were washed twice with PBS. AdipoRed staining solution (10 µL) was added and incubated at 37 °C for 15 min. The chambers were washed three times with fresh DMEM solution.

The cells in the cell culture chambers were washed with PBS. The cells were then fixed for 15 min with 4% paraformaldehyde (PFA, FUJIFILM Wako Pure Chemical) in PBS and permeabilized for 30 min with 0.1% (v/v) Triton X-100 in PBS. The cells were then incubated in blocking buffer containing 5% normal goat serum (Vector), 5% normal donkey serum (Wako), 3% BSA (Sigma-Aldrich), and 0.1% Tween-20 (Nacalai Tesque, Inc.) in PBS at 4 °C for 24 h. After blocking, the cells were treated with mouse anti-human albumin IgG (Stock solution 10 µg mL−1, diluted by blocking buffer into 20 ng mL−1, Thermo Fisher Scientific) as the primary antibody in blocking solution for 24 h. The following day, after washing the excess antibodies three times with PBS (containing 0.1% Tween-20), the cells were treated with Alexa Fluor 647-labeled donkey anti-mouse IgG H&L (Stock solution 1 µg mL−1, diluted by 0.1% Tween-20 PBS into 1 ng mL−1, Jackson ImmunoResearch) at 25 °C for 1 h. The nuclei of the cells were stained with 50 µL of a solution of 300 nM 4,6-diamidino-2-phenylindole (DAPI, Dojindo Laboratories) at 25 °C for 30 min and then washed twice with PBS.

The chips were placed on the stage of a Nikon ECLIPSE Ti inverted fluorescence microscope, which was equipped with a CFI plan fluor 10×/0.30 N.A. objective lens (Nikon, Tokyo, Japan), charge-coupled device (CCD) camera (ORCA-R2; Hamamatsu Photonics, Hamamatsu City, Japan), mercury lamp (Intensilight; Nikon), XYZ automated stage (Ti-S-ER motorized stage with encoders; Nikon), and filter cubes for the fluorescence channels (DAPI and GFP HYQ; Nikon). For image acquisition, the exposure times were set to 100 ms for (DAPI) Hoechst 33258, 5 ms for (GFP HYQ) Calcein AM, 2 s for (TRITC) Annexin V, and 100 ms for the (TRITC) AdipoRed assay.

Firstly, 3 µL of the cell culture medium was collected from a chip and mixed with 96.95 µL of the working solution (50% isopropanol, 25% acetonitrile, and 25% water) and 0.05 µL of the internal standard solution (PA-d4:5 mM, OA-d9:5 mM, in isopropanol). The mixture was vigorously mixed for 30 s and centrifuged at 16,000 × g for 10 min at 4 °C. Subsequently, 2 µL of the supernatant was injected and separated on a Nexera UHPLC system (Shimadzu, Kyoto Japan) using a binary gradient with solvent A (50% acetonitrile, 18% isopropanol, and 32% water) and solvent B (90% isopropanol, 10% acetonitrile, 10 mM ammonium formate, and 0.1% formic acid). The gradient program was as follows: 0% B for 14 min, 100% B for 3 min, and 0% B for 3 min. An ACQUITY UPLC BEH C18 column (130 Å, 1.7 µm, 2.1 mm × 100 mm (Waters, Milford, MA)) was used at 40 °C. The UHPLC eluates were infused online to the LC-MS 8030plus (Shimadzu), which was set to negative electrospray ionization (ESI-) mode. The PA and OA responses were observed by pseudo multiple reaction monitoring (pMRM) with transitions at m/z 255.05 > 255.35 and 281.05 > 281.45, respectively. The pMRM transitions for PA-d4 and OA-d9 were 258.95 > 259.45 and 290.10 > 290.40, respectively. The pMRM transitions were optimized, and peak areas were calculated using LabSolutions software (Shimadzu). The PA and OA responses were normalized to those of PA-d4 and OA-d9 for each sample. All measurements were obtained in triplicate, and the average responses were used. Standard curves were generated by measuring blank culture medium supplemented with increasing amounts of PA and OA.

RNA was purified from the cells using the RNeasy Mini Kit (Qiagen, Hilden, Germany). The microvalves were initially closed to prevent cross-contamination between the HepG2 and Caco-2 cells. After the cells were washed with PBS, 10 µL of trypsin/EDTA (0.04% / 0.03%[v/v]) solution was introduced into the cell culture chambers and incubated at 37 °C with 5% CO2 for 10 min. Subsequently, the cells were harvested with a 10-µL pipette and placed in 1.5-mL tubes. Cells were lysed by adding 350 µL of lysis buffer from the kit and 350 µL of 70% (v/v) ethanol to the tubes. Each solution was transferred to an RNeasy Mini spin column placed in a 2-mL collection tube. The column was centrifuged for 15 s at 8000 × g and the flow-through was discarded. Next, 350 µL buffer RW1 was added to the columns and centrifuged. Subsequently, 80 µL DNase digestion buffer was added to the column and incubated at 25 °C for 15 min. Subsequently, 350 µL of buffer RW1 was added to the column tube and centrifuged again. The column was washed twice with 500 µL of buffer RPE, placed in a new 2 mL tube, and centrifuged. The column was then placed in a new 1.5-mL collection tube and 30 µL of RNase-free water was added to the column. This was followed by centrifugation for 1 min at 8000 × g to elute RNA into the collection tube. RNA quality was evaluated using an Agilent 2100 Bioanalyzer (Agilent Technologies Inc., USA).

mRNA-seq was conducted by Takara Bio Inc. or Oxford NANOPORE Technologies. For Takara Bio Inc, briefly, 50 ng of total RNA from each sample was amplified and synthesized to cDNA using SMART-seq (SMART-Seq v4 Ultra Low Input RNA Kit, Takara Bio). The cDNA was used to make a library using the Nextera DNA Flex Library Prep Kit (Illumina), and the cDNA library was sequenced using NovaSeq 6000 (Illumina). For Oxford NANOPORE Technologies, 40 ng of total RNA was diluted with 9 μL of RNase-free water, mixed with VN primer (Oxford NANOPORE Technologies, UK), and 1 μL of 10 mM dNTPs (New England Biolabs Inc. Ipswich, Massachusetts, USA), and incubated at 65 °C for 5 min to prepare the cDNA library. Separately, 4 μL of 5x RT Buffer (Thermo Fisher Scientific), 1 μL of RNaseOUT (Thermo Fisher Scientific), 1 μL of nuclease-free water, and 2 μL of Strand-Switching Primer (Oxford NANOPORE Technologies) were mixed as the strand-switching buffer. The two solutions were mixed at 42 °C for 2 min, and then 1 μL of Maxima H Minus Reverse Transcriptase (Thermo Fisher Scientific) was added. The mixture was incubated at 42 °C for 90 min and 85 °C for 5 min and stored at 4 °C until use as the cDNA library. Exactly 5 μL aliquot of the cDNA library solution was mixed with 25 μL of 2x LongAmp Taq Master Mix (New England Biolabs Inc.), 1.5 μL of Barcode Primers (Oxford NANOPORE Technologies), and 18.5 μL of nuclease-free water. PCR was performed (95 °C for 30 s; 18 cycles of 95 °C for 15 s, 62 °C for 15 s, 65 °C for 50 s, and 65 °C for 6 min) to barcode the cDNA for multiplexing. The PCR products were stored at 4 °C until use. Subsequently, 1 μL of NEB Exonuclease 1 (New England Biolabs Inc.) was added before incubation at 37 °C for 15 min, followed by incubation at 80 °C for 15 min. Amplified DNA was purified and collected in 12 μL of elution buffer (Oxford NANOPORE Technologies) using Agencourt AMPure XP beads (Beckman Coulter Life Sciences, Indianapolis, IN, USA). A BioAnalyzer 2100 with a High Sensitivity DNA Kit (Agilent Technologies) was used to evaluate the amount and quality of barcoded cDNA. Next, 50 fmol of barcoded cDNA was incubated with 1 μl of Rapid Adapter to make up a total volume of 11 μL, and this was incubated for 5 min at 25 °C. For Nanopore sequencing, 12 μL of the prepared DNA library was mixed with 37.5 μL of Sequencing Buffer (Oxford NANOPORE Technologies) and 25.5 μL of Loading Buffer (Oxford NANOPORE Technologies). This solution was added to a Nanopore Flow Cell (v9.4.1) and a sequencing run was performed for 24 h.

Initially, mRNA-seq reads were mapped to the rRNA, tRNA, or mitochondrial genome sequences using BowTie (v.2.1.0)84. The mapped reads were discarded and not used for the following analysis. The remaining reads were mapped to the human genome (GRCh38) with STAR Aligner (2.7.1a)85 using ENCODE options considering gene annotation and Ensembl (ver.98)86. After the mapping to the genome, gene expression values (Transcripts Per Million reads; TPM) were calculated using RSEM (ver. 1.3.0)87. DEGs of the mono-cultured and co-cultured samples were calculated using DEseq2 (ver. 1.8.2)88. If a gene satisfied the following criteria, it was defined as DEG: p < 0.01, abs (log2(Fold Change)) ≥ 0.263, base mean of raw reads ≥ 31, and average TPM in either sample ≥ 1. GO analysis of DEGs was performed using the WEB-based Gene Set Analysis Toolkit (WebGestalt89). ‘Biological Process noRedundant’ was selected for the database, and ‘genome protein-coding’ genes were selected for the reference set. Protein-coding genes among the DEGs were used as the inputs.

To consider FFAs treatment conditions, we employed ANOVA to select DEGs and characterize the samples. A gene was treated as a DEG if p < 0.05 and abs(log2(Fold Change)) ≥ 0.263 for a combination of any two samples. The expression values of the DEG were used for PCA to characterize the samples. To compare DEG sets with FFAs-minus and FFAs-plus under mono-culture conditions, FFAs-minus and FFAs-plus under the co-culture conditions, and mono- and co-cultured samples under the FFAs-minus condition, DEGs with a PC2 loading of ≥2 or ≤0.5 were used. The results were used to assess changes in gene expression related to FFAs treatment. To compare gene expression profiles according to FFAs treatment under the mono- and co-cultured conditions, we employed Gene Set Enrichment Analysis (GSEA90, ver.4.0.4) with DB files (msigdb.v7.1. symbols.gmt). The fold-change values for protein-coding genes were used as inputs. The bar plots show FDR qvalues for the top four terms under both mono- and co-culture conditions. In the analysis, the PCA plots were drawn with the ggplot2 package in R.

Following microscopic image acquisition, CellProfiler software (Broad Institute of Harvard and MIT, Version 3.1.9)43 was used to estimate cellular features (e.g., cell size and fluorescence intensity). After loading the set of stained images (nuclei and cellular phenotype or function), all images were adjusted with "correctilluminationCalcuate" and "correctilluminationApply" modules to reduce the uneven distribution of background signals. The settings provided in the manual were followed. Traditionally, the illumination function is calculated from "Background" with a block size of 20–100 for each image individually. And the Smoothing method was selected as "Median filter," with a filter size of 50. The other options followed the default settings. The mathematical methods of "correctilluminationApply" was "Substract." With the corrected nuclei images, individual cells were identified by using Otsu's method43 in the "IdentifyPrimaryObjects" module as primary objects, followed by "IdentifySecondaryObjects" to evaluate the cellular phenotype/function, automatically. Then, cellular features were calculated using the "MeasureObjectSizeShape" and "MeasureObjectIntensity" modules for both primary and secondary objects. Further analysis of single-cell morphological descriptors was performed using t-SNE in the open-source Orange 3 software (Version 3.23.1; Bioinformatics Laboratory, Faculty of Computer and Information Science, University of Ljubljana, Slovenia)91.

The results expression, statistical tests used, the sample sizes and number of independent experiments are mentioned in the figure legends. All cell-based assays were carried out in triplicate samples with three independent experiments. The Ridgeline plots and t-SNE analysis charts were based on the integration of three independent experiments. The Tukey–Kramer test and Student's t test were performed using R software (ver. 3.5.2; https://www.r-project.org/). For mRNA-seq, all the experiments were conducted with at least two independent experiments.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

The mRNA-seq data has been deposited in the NCBI Gene Expression Omnibus (GEO) under accession number GSE152091 and GSE206417. All the manuscript graphs and charts of Ridgeline plots and t-SNE analysis are available in Supplementary Data 1 to 6.

Chalasani, N. et al. The diagnosis and management of nonalcoholic fatty liver disease: Practice guidance from the American Association for the Study of Liver Diseases. Hepatology 67, 328–357 (2018).

Article PubMed Google Scholar

Younossi, Z. M. Non-alcoholic fatty liver disease—a global public health perspective. J. Hepatol. 70, 531–544 (2019).

Article PubMed Google Scholar

Byrne, C. D. & Targher, G. NAFLD: A multisystem disease. J. Hepatol. 62, S47–S64 (2015).

Article PubMed Google Scholar

Friedman, S. L., Neuschwander-Tetri, B. A., Rinella, M. & Sanyal, A. J. Mechanisms of NAFLD development and therapeutic strategies. Nat. Med. 24, 908–922 (2018).

Article CAS PubMed PubMed Central Google Scholar

Estes, C., Razavi, H., Loomba, R., Younossi, Z. & Sanyal, A. J. Modeling the epidemic of nonalcoholic fatty liver disease demonstrates an exponential increase in burden of disease. Hepatology 67, 123–133 (2018).

Article CAS PubMed Google Scholar

Buzzetti, E., Pinzani, M. & Tsochatzis, E. A. The multiple-hit pathogenesis of non-alcoholic fatty liver disease (NAFLD). Metabolism 65, 1038–1048 (2016).

Article CAS PubMed Google Scholar

Wiest, R., Albillos, A., Trauner, M., Bajaj, J. S. & Jalan, R. Targeting the gut-liver axis in liver disease. J. Hepatol. 67, 1084–1103 (2017).

Article CAS PubMed Google Scholar

Sumida, Y. & Yoneda, M. Current and future pharmacological therapies for NAFLD/NASH. J. Gastroenterol 53, 362–376 (2018).

Article CAS PubMed Google Scholar

Ohtani, N. & Kawada, N. Role of the gut–liver axis in liver inflammation, fibrosis, and cancer: a special focus on the gut microbiota relationship. Hepatol. Commun. 3, 456–470 (2019).

Article PubMed PubMed Central Google Scholar

Leung, C., Rivera, L., Furness, J. B. & Angus, P. W. The role of the gut microbiota in NAFLD. Nat. Rev. Gastroenterol. Hepatol. 13, 412–425 (2016).

Article CAS PubMed Google Scholar

Marchesi, J. R. et al. The gut microbiota and host health: a new clinical frontier. Gut 65, 330–339 (2016).

Article PubMed Google Scholar

Clemente, M. G., Mandato, C., Poeta, M. & Vajro, P. Pediatric non-alcoholic fatty liver disease: Recent solutions, unresolved issues, and future research directions. World J. Gastroenterol. 22, 8078–8093 (2016).

Article PubMed PubMed Central Google Scholar

Rotman, Y. & Sanyal, A. J. Current and upcoming pharmacotherapy for non-alcoholic fatty liver disease. Gut 66, 180–190 (2017).

Article CAS PubMed Google Scholar

Abdalkader, R. & Kamei, K. I. Multi-corneal barrier-on-a-chip to recapitulate eye blinking shear stress forces. Lab Chip 20, 1410–1417 (2020).

Article CAS PubMed Google Scholar

Sung, J. H., Wang, Y. I., Kim, J. H., Lee, J. M. & Shuler, M. L. Application of chemical reaction engineering principles to "body-on-a-chip" systems. AIChE J. 64, 4351–4360 (2018).

Article CAS PubMed PubMed Central Google Scholar

Abaci, H. E. & Shuler, M. L. Human-on-a-chip design strategies and principles for physiologically based pharmacokinetics/pharmacodynamics modeling. Integr. Biol. 7, 383–391 (2015).

Article Google Scholar

Ronaldson-Bouchard, K. & Vunjak-Novakovic, G. Organs-on-a-chip: a fast track for engineered human tissues in drug development. Cell Stem Cell 22, 310–324 (2018).

Article CAS PubMed PubMed Central Google Scholar

Schurdak, M. et al. Applications of the microphysiology systems database for experimental ADME-Tox and disease models. Lab Chip 20, 1472–1492 (2020).

Article CAS PubMed PubMed Central Google Scholar

Grego, S. et al. Systems biology for organotypic cell cultures. ALTEX 34, 301–310 (2017).

Article PubMed Google Scholar

Lee, S. Y. & Sung, J. H. Gut–liver on a chip toward an in vitro model of hepatic steatosis. Biotechnol. Bioeng. 115, 2817–2827 (2018).

Article CAS PubMed Google Scholar

Müller, F. A. & Sturla, S. J. Human in vitro models of nonalcoholic fatty liver disease. Curr. Opin. Toxicol 16, 9–16 (2019).

Article Google Scholar

Chen, W. L. K. et al. Integrated gut/liver microphysiological systems elucidates inflammatory inter-tissue crosstalk. Biotechnol. Bioeng. 114, 2648–2659 (2017).

Article CAS PubMed PubMed Central Google Scholar

Tsamandouras, N. et al. Integrated gut and liver microphysiological systems for quantitative in vitro pharmacokinetic studies. AAPS J. 19, 1499–1512 (2017).

Article CAS PubMed Google Scholar

Costa, J. & Ahluwalia, A. Advances and Current Challenges in Intestinal in vitro Model Engineering: A Digest. Front. Bioeng. Biotechnol 7, 1–14 (2019).

Article Google Scholar

De Gregorio, V. et al. Intestine-liver axis on-chip reveals the intestinal protective role on hepatic damage by emulating ethanol first-pass metabolism. Front. Bioeng. Biotechnol. 8, 163 (2020).

Article PubMed PubMed Central Google Scholar

Jeon, J. W., Lee, S. H., Kim, D. & Sung, J. H. In vitro hepatic steatosis model based on gut–liver-on-a-chip. Biotechnol. Prog. 37, 1–13 (2021).

Article Google Scholar

Yang, J. et al. Gut-liver-axis microphysiological system for studying cellular fluidic shear stress and inter-tissue interaction. Biomicrofluidics 16, 44113 (2022).

Article CAS Google Scholar

Wen, X., Yoshimoto, K., Yamanaka, M., Terada, S. & Kamei, K. In vitro nonalcoholic fatty liver disease model with cyclo-olefin-polymer-based microphysiological systems. Organs-on-a-Chip 3, 100010 (2021).

Article CAS Google Scholar

Kamei, K. I. et al. Integrated heart/cancer on a chip to reproduce the side effects of anti-cancer drugs: In vitro. RSC Adv. 7, 36777–36786 (2017).

Article CAS Google Scholar

Yu, Z. T. F. et al. Integrated microfluidic devices for combinatorial cell-based assays. Biomed. Microdevices 11, 547–555 (2009).

Article CAS PubMed PubMed Central Google Scholar

Kamei, K. I. et al. An integrated microfluidic culture device for quantitative analysis of human embryonic stem cells. Lab Chip 9, 555–563 (2009).

Article CAS PubMed Google Scholar

Huang, B., Wu, H., Kim, S. & Zare, R. N. Coating of poly(dimethylsiloxane) with n-dodecyl-β-D-maltoside to minimize nonspecific protein adsorption. Lab Chip 5, 1005–1007 (2005).

Article CAS PubMed Google Scholar

Kamei, K. I. et al. Phenotypic and transcriptional modulation of human pluripotent stem cells induced by nano/microfabrication materials. Adv. Healthc. Mater. 2, 287–291 (2013).

Article CAS PubMed Google Scholar

Regehr, K. J. et al. Biological implications of polydimethylsiloxane-based microfluidic cell culture. Lab Chip 9, 2132–2139 (2009).

Article CAS PubMed PubMed Central Google Scholar

Wong, I. & Ho, C. M. Surface molecular property modifications for poly(dimethylsiloxane) (PDMS) based microfluidic devices. Microfluid. Nanofluidics 7, 291–306 (2009).

Article CAS PubMed PubMed Central Google Scholar

Kim, H. J. & Ingber, D. E. Gut-on-a-Chip microenvironment induces human intestinal cells to undergo villus differentiation. Integr. Biol. 5, 1130–1140 (2013).

Article CAS Google Scholar

Kulthong, K. et al. Implementation of a dynamic intestinal gut-on-a-chip barrier model for transport studies of lipophilic dioxin congeners. RSC Adv. 8, 32440–32453 (2018).

Article CAS PubMed PubMed Central Google Scholar

Gómez-Lechón, M. J. et al. A human hepatocellular in vitro model to investigate steatosis. Chem. Biol. Interact. 165, 106–116 (2007).

Article PubMed Google Scholar

Marra, F. & Svegliati-Baroni, G. Lipotoxicity and the gut-liver axis in NASH pathogenesis. J. Hepatol. 68, 280–295 (2018).

Article CAS PubMed Google Scholar

Hirsova, P., Ibrabim, S. H., Gores, G. J. & Malhi, H. Thematic review series: Lipotoxicity: Many roads to cell dysfunction and cell death lipotoxic lethal and sublethal stress signaling in hepatocytes: Relevance to NASH pathogenesis. J. Lipid Res. 57, 1758–1770 (2016).

Article CAS PubMed PubMed Central Google Scholar

Ogawa, Y. et al. Palmitate-induced lipotoxicity is crucial for the pathogenesis of nonalcoholic fatty liver disease in cooperation with gut-derived endotoxin. Sci. Rep. 8, 11365 (2018).

Article PubMed PubMed Central Google Scholar

McQuin, C. et al. CellProfiler 3.0: Next-generation image processing for biology. PLoS Biol. 16, e2005970 (2018).

Article PubMed PubMed Central Google Scholar

Carpenter, A. E. et al. CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol. 7, R100 (2006).

Article PubMed PubMed Central Google Scholar

Sirenko, O., Hesley, J., Rusyn, I. & Cromwell, E. F. High-content high-throughput assays for characterizing the viability and morphology of human iPSC-derived neuronal cultures. Assay Drug Dev. Technol. 12, 536–547 (2014).

Article CAS PubMed PubMed Central Google Scholar

Gilbert, D. F. et al. A novel multiplex cell viability assay for high-throughput RNAi screening. PLoS ONE 6, e28338 (2011).

Article CAS PubMed PubMed Central Google Scholar

Van Der Maaten, L. & Hinton, G. Visualizing data using t-SNE. J. Mach. Learn. Res. 9, 2579–2625 (2008).

Google Scholar

Caldez, M. J., Bjorklund, M. & Kaldis, P. Cell cycle regulation in NAFLD: when imbalanced metabolism limits cell division. Hepatol. Int. 14, 463–474 (2020).

Article PubMed Google Scholar

Kulthong, K. et al. Microfluidic chip for culturing intestinal epithelial cell layers: characterization and comparison of drug transport between dynamic and static models. Toxicol. Vitr 65, 104815 (2020).

Article CAS Google Scholar

Villenave, R. et al. Human gut-on-a-chip supports polarized infection of coxsackie B1 virus in vitro. PLoS ONE 12, 1–17 (2017).

Article Google Scholar

Kim, H. J., Huh, D., Hamilton, G. & Ingber, D. E. Human gut-on-a-chip inhabited by microbial flora that experiences intestinal peristalsis-like motions and flow. Lab Chip 12, 2165–2174 (2012).

Article CAS PubMed Google Scholar

Delon, L. C. et al. A systematic investigation of the effect of the fluid shear stress on Caco-2 cells towards the optimization of epithelial organ-on-chip models. Biomaterials 225, 119521 (2019).

Article CAS PubMed Google Scholar

Du, Y. et al. Mimicking liver sinusoidal structures and functions using a 3D-configured microfluidic chip. Lab Chip 17, 782–794 (2017).

Article CAS PubMed Google Scholar

Wang, X. et al. Fluid shear stress promotes autophagy in hepatocellular carcinoma cells. Int. J. Biol. Sci. 14, 1277–1290 (2018).

Article CAS PubMed PubMed Central Google Scholar

Camakaris, J., Voskoboinik, I. & Mercer, J. F. Molecular mechanisms of copper homeostasis. Biochem. Biophys. Res. Commun. 261, 225–232 (1999).

Article CAS PubMed Google Scholar

Khan, S., Zafar, A. & Naseem, I. Redox cycling of copper by coumarin-di(2-picolyl)amine hybrid molecule leads to ROS-mediated modulation of redox scavengers, DNA damage and cell death in diethylnitrosamine induced hepatocellular carcinoma. Bioorg. Chem. 99, 103818 (2020).

Article CAS PubMed Google Scholar

Parsekar, S. U. et al. DNA binding, cleavage and cytotoxicity studies of three mononuclear Cu(II) chloro-complexes containing N–S donor Schiff base ligands. J. Biol. Inorg. Chem. 23, 1331–1349 (2018).

Article CAS PubMed Google Scholar

Foo, J. B. et al. Copper complex derived from S-benzyldithiocarbazate and 3-acetylcoumarin induced apoptosis in breast cancer cell. BioMetals 31, 505–515 (2018).

Article CAS PubMed Google Scholar

Gori, M. et al. Palmitic acid affects intestinal epithelial barrier integrity and permeability in vitro. Antioxidants 9, 417 (2020).

Article CAS PubMed PubMed Central Google Scholar

Storniolo, C. E., Cabral, M., Busquets, M. A., Martín-Venegas, R. & Moreno, J. J. Dual behavior of long-chain fatty acids and their cyclooxygenase/lipoxygenase metabolites on human intestinal Caco-2 cell growth. Front. Pharmacol 11, 529976 (2020).

Article CAS PubMed PubMed Central Google Scholar

Embade, N., Valeron, P. F., Aznar, S., Lopez-Collazo, E. & Lacal, J. C. Apoptosis induced by Rac GTPase correlates with induction of FasL and ceramides production. Mol. Biol. Cell 11, 4347–4358 (2000).

Article CAS PubMed PubMed Central Google Scholar

Simon, S. & Schubert, R. Inhibitory effect of phospholipids on P-glycoprotein: cellular studies in Caco-2, MDCKII mdr1 and MDCKII wildtype cells and P-gp ATPase activity measurements. Biochim. Biophys. Acta—Mol. Cell Biol. Lipids 1821, 1211–1223 (2012).

Article CAS Google Scholar

Berryman, S., Matthews, K., Lee, J. H., Duffy, S. P. & Ma, H. Image-based phenotyping of disaggregated cells using deep learning. Commun. Biol 3, 674 (2020).

Article PubMed PubMed Central Google Scholar

Ramos, M. J., Bandiera, L., Menolascina, F. & Fallowfield, J. A. In vitro models for non-alcoholic fatty liver disease: Emerging platforms and their applications. iScience 25, 103549 (2022).

Article PubMed Google Scholar

Grouls, M. et al. Differential gene expression in iPSC-derived human intestinal epithelial cell layers following exposure to two concentrations of butyrate, propionate and acetate. Sci. Rep. 12, 1–15 (2022).

Article Google Scholar

Green, C. J. et al. Studying non-alcoholic fatty liver disease: the ins and outs of in vivo, ex vivo and in vitro human models. Horm. Mol. Biol. Clin. Investig. 41, 20180038 (2020).

CAS Google Scholar

Kostrzewski, T. et al. Modelling human liver fibrosis in the context of non-alcoholic steatohepatitis using a microphysiological system. Commun. Biol. 4, 1080 (2021).

Article CAS PubMed PubMed Central Google Scholar

Ganji, S. H., Kashyap, M. L. & Kamanna, V. S. Niacin inhibits fat accumulation, oxidative stress, and inflammatory cytokine IL-8 in cultured hepatocytes: Impact on non-alcoholic fatty liver disease. Metabolism 64, 982–990 (2015).

Article CAS PubMed Google Scholar

Thomson, J. A. Embryonic stem cell lines derived from human blastocysts. Science 282, 1145–1147 (1998).

Article CAS PubMed Google Scholar

Takahashi, K. et al. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell 131, 861–872 (2007).

Article CAS PubMed Google Scholar

Yu, J. et al. Induced pluripotent stem cell lines derived from human somatic cells. Science 318, 1917–1920 (2007).

Article CAS PubMed Google Scholar

Lowry, W. E. et al. Generation of human induced pluripotent stem cells from dermal fibroblasts. Proc. Natl. Acad. Sci. USA 105, 2883–2888 (2008).

Article CAS PubMed PubMed Central Google Scholar

Makarova, K. S. et al. Evolutionary classification of CRISPR–Cas systems: a burst of class 2 and derived variants. Nat. Rev. Microbiol 18, 67–83 (2020).

Article CAS PubMed Google Scholar

Yang, J. et al. Multilayered microfluidic device for controllable flow perfusion of gut-liver on a chip. 21st International Conference on Solid-State Sensors, Actuators and Microsystems (Transducers). 176–179 (IEEE, 2021).

van Meer, B. J. et al. Small molecule absorption by PDMS in the context of drug response bioassays. Biochem. Biophys. Res. Commun. 482, 323–328 (2017).

Article PubMed PubMed Central Google Scholar

Berthier, E., Young, E. W. K. & Beebe, D. Engineers are from PDMS-land, Biologists are from Polystyrenia. Lab Chip 12, 1224–1237 (2012).

Article CAS PubMed Google Scholar

Huh, D. et al. Reconstituting organ-level lung functions on a chip. Science 328, 1662–1668 (2010).

Article CAS PubMed PubMed Central Google Scholar

Yoshimoto, K. et al. Recapitulation of human embryonic heartbeat to promote differentiation of hepatic endoderm to hepatoblasts. Front. Bioeng. Biotechnol. 8, 1064 (2020).

Article Google Scholar

Wang, M. et al. Application of perfluoropolyether elastomers in microfluidic drug metabolism assays. Int. J. Pharm 627, 122253 (2022).

Article CAS PubMed Google Scholar

Liao, S., He, Y., Chu, Y., Liao, H. & Wang, Y. Solvent-resistant and fully recyclable perfluoropolyether-based elastomer for microfluidic chip fabrication. J. Mater. Chem. A 7, 16249–16256 (2019).

Article CAS Google Scholar

Campbell, S. B. et al. Beyond polydimethylsiloxane: alternative materials for fabrication of organ-on-a-chip devices and microphysiological systems. ACS Biomater. Sci. Eng. 7, 2880–2899 (2020).

Article PubMed Google Scholar

Albillos, A., de Gottardi, A. & Rescigno, M. The gut-liver axis in liver disease: pathophysiological basis for therapy. J. Hepatol. 72, 558–577 (2020).

Article CAS PubMed Google Scholar

Ma, X. et al. Experimental study of numerical optimization for 3-D microstructuring using DMD-based grayscale lithography. J. Microelectromech. Syst. 24, 1856–1867 (2015).

Article Google Scholar

Kato, Y., Hirai, Y., Kamei, K., Tsuchiya, T. & Tabata, O. Microfluidic device to interconnect multiple organs via fluidic circulation: towards body-on-a-chip. 18th International Conference on Solid-State Sensors, Actuators and Microsystems (Transducers). 1549–1552 (IEEE, 2015).

Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).

Article CAS PubMed PubMed Central Google Scholar

Dobin, A. et al. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).

Article CAS PubMed Google Scholar

Hunt, S. E. et al. Ensembl variation resources. Database 2018, 1–12 (2018).

Article Google Scholar

Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12, 323 (2011).

Article CAS PubMed PubMed Central Google Scholar

Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15, 550 (2014).

Article PubMed PubMed Central Google Scholar

Wang, J., Vasaikar, S., Shi, Z., Greer, M. & Zhang, B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 45, W130–W137 (2017).

Article CAS PubMed PubMed Central Google Scholar

Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 102, 15545–15550 (2005).

Article CAS PubMed PubMed Central Google Scholar

Demšar, J. et al. Orange: Data mining toolbox in python. J. Mach. Learn. Res. 14, 2349–2353 (2013).

Google Scholar

Download references

The authors would like to thank Ms. Miyako Fujita and Dr. Satoshi Imamura for technical assistance. Funding was generously provided by the Japan Society for the Promotion of Science (JSPS) (16K14660, 17H02083, 18KK0306, 19H02572, and 21H01728), Terumo Life Science Foundation, Ebara Hatakeyama Memorial Foundation, Japan Agency for Medical Research and Development (AMED) (17937667), and LiaoNing Revitalization Talents Program (XLYC1902061). MT was supported by the Nakatani Foundation for the Advancement of Measuring Technologies in Biomedical Engineering. Part of this work was supported by the Nanotechnology Platform Project within MEXT, Japan through the Kyoto University Nano Technology Hub. WPI-iCeMS is supported by the World Premier International Research Centre Initiative (WPI), MEXT, Japan. We would also like to thank Editage (www.editage.com) for the English language editing.

Department of Micro Engineering, Kyoto University, Kyotodaigaku-Katsura, Nishikyo-ku, Kyoto, 615-8540, Japan

Jiandong Yang, Yoshikazu Hirai, Marika Trumm, Toshiyuki Tsuchiya & Osamu Tabata

Department of Mechanical Engineering and Science, Kyoto University, Kyotodaigaku-Katsura, Nishikyo-ku, Kyoto, 615-8540, Japan

Yoshikazu Hirai

Medical Research Support Center, Graduate School of Medicine, Kyoto University, Yoshida Konoe-cho, Kyoto, 606-8501, Japan

Kei Iida

Faculty of Science and Engineering, Kindai University, 3-4-1 Kowakae, Higashiosaka, Osaka, 577-8502, Japan

Kei Iida & Shinji Ito

Institute for Integrated Cell-Material Sciences, Kyoto University, Yoshida-Ushinomiya-cho, Sakyo-ku, Kyoto, 606-8501, Japan

Marika Trumm, Shiho Terada, Risako Sakai, Osamu Tabata & Ken-ichiro Kamei

Institute for Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, 69120, Germany

Marika Trumm

Faculty of Engineering/Graduate School of Engineering, Kyoto University of Advanced Science, Gotanda-cho, Yamanouchi, Ukyo-ku, Kyoto, 615-8577, Japan

Osamu Tabata

Wuya College of Innovation, Shenyang Pharmaceutical University, 110016, Liaoning, China

Ken-ichiro Kamei

Department of Pharmaceutics, Shenyang Pharmaceutical University, 110016, Liaoning, China

Ken-ichiro Kamei

Programs of Biology and Bioengineering, Divisions of Science and Engineering, New York University Abu Dhabi, Abu Dhabi, UAE

Ken-ichiro Kamei

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

J.Y., Y.H., T.T., O.T., and K.K. conceptualized the study. J.Y. and Y.H. fabricated the device. J.Y., R.S., and S.T. performed the cell culture experiments and analyses. S.I. conducted the UHPLC-MS/MS experiments. J.Y., R.S., S.T., and K.K. performed the image analysis. J.Y., S.T., and M.T. performed the RNA experiments. I.K. and K.K. performed the RNA-seq analyses. All authors contributed to the data analysis, discussion, and interpretation. J.Y., Y.H., and K.K. wrote and revised the manuscript with inputs from all authors.

Correspondence to Yoshikazu Hirai or Ken-ichiro Kamei.

The authors declare no competing interests.

Communications Biology thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: Anam Akhtar and Christina Karlsson Rosenthal.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Yang, J., Hirai, Y., Iida, K. et al. Integrated-gut-liver-on-a-chip platform as an in vitro human model of non-alcoholic fatty liver disease. Commun Biol 6, 310 (2023). https://doi.org/10.1038/s42003-023-04710-8

Download citation

Received: 25 June 2020

Accepted: 14 March 2023

Published: 23 March 2023

DOI: https://doi.org/10.1038/s42003-023-04710-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.