Self-Organizing Human Induced Pluripotent Stem Cell Hepatocyte 3D Organoids Inform the Biology of the Pleiotropic TRIB1 Gene.

Establishment of a physiologically relevant human hepatocyte-like cell system for in vitro translational research has been hampered by the limited availability of cell models that accurately reflect human biology and the pathophysiology of human disease. Here we report a robust, reproducible, and scalable protocol for the generation of hepatic organoids from human induced pluripotent stem cells (hiPSCs) using short exposure to nonengineered matrices. These hepatic organoids follow defined stages of hepatic development and express higher levels of early (hepatocyte nuclear factor 4A [HNF4A], prospero-related homeobox 1 [PROX1]) and mature hepatic and metabolic markers (albumin, asialoglycoprotein receptor 1 [ASGR1], CCAAT/enhancer binding protein α [C/EBPα]) than two-dimensional (2D) hepatocyte-like cells (HLCs) at day 20 of differentiation. We used this model to explore the biology of the pleiotropic TRIB1 (Tribbles-1) gene associated with a number of metabolic traits, including nonalcoholic fatty liver disease and plasma lipids. We used genome editing to delete the TRIB1 gene in hiPSCs and compared TRIB1-deleted iPSC-HLCs to isogenic iPSC-HLCs under both 2D culture and three-dimensional (3D) organoid conditions. Under conventional 2D culture conditions, TRIB1-deficient HLCs showed maturation defects, with decreased expression of late-stage hepatic and lipogenesis markers. In contrast, when cultured as 3D hepatic organoids, the differentiation defects were rescued, and a clear lipid-related phenotype was noted in the TRIB1-deficient induced pluripotent stem cell HLCs.
Conclusion: This work supports the potential of genome-edited hiPSC-derived hepatic 3D organoids in exploring human hepatocyte biology, including the functional interrogation of genes identified through human genetic investigation.

two‐dimensional

three‐dimensional

alpha‐fetoprotein

acetyl‐coenzyme A carboxylase

albumin

apolipoprotein B

apolipoprotein CIII

asialoglycoprotein receptor 1

bone morphogenetic protein 4

complementary DNA

CCAAT/enhancer binding protein α

4′,6‐diamidino‐2‐phenylindole

definitive endoderm

diacylglycerol acyltransferase

delta‐like 1 homolog

Dulbecco’s phosphate buffered saline

fluorescence‐activated cell sorting

fatty acid synthase

fibroblast growth factor 2

forkhead box A2

glyceraldehyde 3‐phosphate dehydrogenase

green fluorescent protein

genome‐wide association studies

human CCAAT/enhancer binding protein α

human growth factor

hematopoietically expressed homeobox protein

human‐induced pluripotent stem cell

hepatocyte nuclear factor 4 alpha

induced pluripotent stem cell

knockout

low‐density lipoprotein

liver‐specific KO mice

multidrug resistance protein 1

messenger RNA

periodic acid–Schiff

polymerase chain reaction

paraformaldehyde

prospero‐related homeobox 1

room temperature

stearoyl‐coenzyme A desaturase

single‐stranded oligo nucleotide

triglyceride

Tribbles‐1

wild‐type

tight junction protein 1

Large‐scale genome‐wide association studies (GWAS) have identified thousands of genomic loci associated with human traits and diseases, most of which lack a clear causal gene or mechanism.( , , ) Model systems that reflect human biology are needed to systematically examine these findings and translate them into meaningful biology. Human induced pluripotent stem cells (hiPSCs) offer the opportunity to provide an inexhaustible supply of human cells that reflect naturally occurring genetic variation or can be genome‐edited for experimental purposes.( ) Genetic variation in liver‐expressed genes and the effect on phenotype and disease has been of major interest, given the major impact of the many physiological processes performed by hepatocytes. Differentiation of hiPSCs into hepatocyte‐like cells has proved to be an important approach to understanding human biology, and there are several established protocols available to date.( , , ) These hiPSC‐derived hepatocytes in a standard two‐dimensional (2D) monolayer culture system efficiently express hepatic markers and mimic some hepatocyte functions, enabling their use in drug toxicity screening( ) and understanding biology.( ) However, this approach suffers from specification defects. Most of the derived monolayer hepatocytes either lack expression of some key transcription factors or express them at suppressed levels, in turn making them metabolically immature.( , , ) Another limitation is that many gene mutations in this system are studied in isolation, which does not reflect the true underlying biology.( , )

Recently, an induced pluripotent stem cell (iPSC)–derived hepatocyte three‐dimensional (3D) organoid culture system has been reported to establish heterotypic interactions and mimic the in vivo environment more closely in terms of distinct cellular behavior.( , , , , , , ) However, most of the 3D approaches involve complex bioengineered matrices and synthetic scaffolds that need skilled handling expertise. Moreover, these different 3D models have not been rigorously evaluated in terms of their ability to functionally assess genes at GWAS loci in a way that reflects human biology more faithfully than conventional 2D iPSC hepatocytes.

To address these limitations, we developed a robust, growth factor–defined, scalable hepatic 3D organoid generation method from hiPSCs with short exposure to a laminin‐enriched, nonengineered matrix for the purpose of interrogating lipid GWAS loci. We then applied it to understanding the role of the pleiotropic lipid/metabolic GWAS gene TRIB1. Variants at the TRIB1 locus are associated with liver transaminases and nonalcoholic fatty liver disease (NAFLD), plasma triglycerides, total cholesterol, low‐density lipoprotein (LDL) cholesterol, high‐density lipoprotein cholesterol, adiponectin, and coronary artery disease.( , , ) Studies in mice have indicated that the protein product of TRIB1, namely Tribbles‐1, posttranslationally regulates the degradation of CCAAT/enhancer binding protein α (C/EBPα) in the liver, which in turn contributes to the dysregulation of lipid‐related gene expression and thus lipid phenotypes.( ) However, the role of TRIB1 in liver lipid biology remains uncertain, and an in vitro model is needed for more detailed functional investigation. We used CRISPR/Cas9 to delete TRIB1 in a human iPSC line, and show that iPSC‐derived 3D hepatic organoids lacking TRIB1 more faithfully expressed an expected lipid phenotype than 2D hepatocytes derived from the same iPSCs. This study emphasizes the potential advantages of hiPSC‐derived 3D hepatic organoids to interrogate the effect of genetic variation on human liver lipid metabolism.

Results

Generation of Self‐Organizing Functional Hepatic Organoids From hiPSCs

To generate hepatic organoids from hiPSCs, we developed a simple approach by making modest modifications to a conventional growth factor–induced protocol and providing them the 3D environment to overcome specification and maturation defects. Detailed methodology is described in the Methods section, but in brief, hiPSCs were passaged and aggregated in 96‐well U‐bottom low‐attachment plates. Later, laminin‐rich matrix was added only in the definitive endoderm (DE) stage for 4 days. To achieve homogeneous and robust DE induction, the DE Kit (STEMCELL Technologies, Vancouver, Canada) was used, which resulted in robust DE differentiation and greater than 90% forkhead box A2+ (FOXA2+) and GATA binding protein 4+ (GATA4+) cells (Supporting Fig. S1). On day 5, these organoids were moved to a medium containing bone morphogenetic protein 4 (BMP4), fibroblast growth factor 2 (FGF2), and keratinocyte growth factor (KGF) to promote foregut differentiation. On day 10, the developing foregut organoids were shifted to a medium containing human growth factor (HGF) and 0.5% DMSO to promote hepatoblast formation. For maturation on day 15, aggregates were moved to commercial hepatocyte culture medium along with Oncostatin M and Dexamethasone (R&D Systems, Minneapolis, MN) (Fig. 1A). To compare the functional capacity of the organoids with 2D monolayer cultures, a previously reported timeline was followed,( , , , ) and organoids were harvested at day 20 for comparison with 2D HLCs (Supporting Fig. S2). Immunostaining showed robust expression of early and mature hepatic markers such as hepatocyte nuclear factor 4 alpha (HNF4A), alpha‐fetoprotein (AFP), and albumin (ALB), in day 20 organoids (Fig. 1Bi‐iii). These organoids also possess characteristic hepatocyte features and express mature markers including arginase, a canalicular marker multidrug resistance protein 1 (MDR1), and tubular marker HNF1β (Fig. 1Biv‐vi). Functional assays were performed, which indicated efficient BODIPY‐LDL (Fisher Scientific, Waltham, MA) and CDFDA (Invitrogen, Carlsbad, CA) uptake (Fig. 1Ci‐ii). Glycogen and lipid droplet accumulation was also observed as analyzed by periodic acid–Schiff (PAS), Oil Red O (ORO) staining with corresponding electron microscopy images (Fig. 1Ciii‐vi), establishing the generation of functional hiPSC‐derived hepatic organoids.

FIG. 1

Hepatic organoids generated from hiPSCs exhibit liver‐like features more effectively than 2D HLCs. (A) Schematic showing directed differentiation of 3D hepatic organoids from hiPSCs. (B) Representative images of immunostaining showing expression of early (i‐ii; HNF4A, green; AFP, red) and mature (iii‐v; ALB, green; Arginase, green; MDR1, green) hepatic and tubular (vi; HNF1β, pink) markers at day 20 of organoid differentiation. (Ci‐ii) LDL‐uptake and bile‐uptake assays were performed with hepatic organoids at day 20, showing uptake of fluorescent BODIPY‐LDL (red) and fluorescent substrate CDFDA (green). Hepatic organoids exhibit functional liver‐like properties such as glycogen accumulation as shown by PAS staining (pink) (iii‐iv) and lipid droplet formation as shown by ORO (red) (v‐vi) with corresponding electron microscopy images showing glycogen and lipid droplets. (D) Hepatic organoids efficiently express early, late, and metabolic markers at day 20 organoids versus 2D HLCs calculated relative to GAPDH (reference gene) using the 2−ddCt method, normalized to a mean expression level in 2D HLCs. Data are representative of three independent experiments. (E) Increased HNF4A+ cells in day 20 organoids versus 2D HLCs when measured by FACS. (F‐H) Hepatic organoids secrete higher amounts of ALB, APOB, and APOCIII versus 2D HLCs when compared at day 20 of differentiation by enzyme‐linked immunosorbent assay. Values represent the mean ± SEM; P values were calculated with Welch’s t tests. DAPI (4′,6‐diamidino‐2‐phenylindole) is shown in blue in all of the immunofluorescence images. Scale bar: 125 µm (Biii,v; Ci‐iii,v) and 275 µm (Bi,ii,iv,vi); and electron microscopy: 500 nm (Civ,vi).

Hepatic Organoids Are More Effectively Differentiated Than Monolayer 2D HLC Cultures

To determine the differentiation efficiency of hiPSCs to hepatocytes, organoids were analyzed for gene expression of early hepatic markers such as HNF4A, AFP, and prospero‐related homeobox 1 (PROX1); maturation markers such as ALB, asialoglycoprotein receptor 1 (ASGR1), and C/EBPα; and metabolic markers including apolipoprotein B (APOB), apolipoprotein CIII (APOC3), angiopoietin‐like 3 (ANGPTL3), microsomal triglyceride transfer protein (MTTP), cytochrome P450 3A4 (CYP3A4), and cholesterol 7 alpha‐hydroxylase (CYP7A1) at day 20 of differentiation (Fig. 1D). Hepatic organoids express these markers at a higher level than 2D HLCs. Similar observations were made with a K3 hiPSC line derived from a different genetic background (Supporting Fig. S3). Interestingly, hepatic 3D organoids showed enhanced expression of canonical hepatocyte and metabolic genes than 2D HLCs and even primary human hepatocytes, but not when compared with human liver samples (Supporting Fig. S3A‐C). Fluorescence‐activated cell sorting (FACS) quantification on day 20 hepatic 3D organoids demonstrated an approximate 1.5‐fold higher HNF4A compared with day 20 2D HLCs (Fig. 1E). Following continuation to day 30, we noted an approximate 2‐5‐fold increase in expression of hepatic markers over day 20 3D organoids (Supporting Fig. S1B); in marked contrast, the 2D HLCs started to dedifferentiate after day 20 and began to detach around day 22‐25.( ) Secretory proteins such as albumin, APOC3, and APOB were measured in the media, and were observed to be significantly higher in the media from hepatic organoids versus 2D HLCs (Fig. 1F‐H).

Hepatic Organoids for Modeling the Role of TRIB1

To determine whether hiPSC hepatocytes can recapitulate TRIB1 lipid biology, we generated TRIB1 knockout (KO) hiPSCs using genome engineering. We took our in‐house‐derived and fully characterized wild‐type (WT) hiPSC line, Sv20. Using CRISPR/Cas9, a stop codon was introduced in the first exon by frameshift mutation, eight amino acids after the start codon. Homology‐directed repair was used by designing a single‐stranded oligo nucleotide (ssODN), which introduced a stop codon and HindIII restriction site to screen the clones (Fig. 2A). Three clones each for WT and KO were selected for further characterization. KO clones showed restriction by HindIII, and later Sanger sequencing confirmed the insertion of ssODN (Fig. 2B,C). Due to the lack of availability of any reliable human TRIB1 antibody, we performed Sanger sequencing on the generated complementary DNA (cDNA) to confirm the genotype of KO clones. Fully characterized WT and KO clones were progressed for differentiation studies (Fig. 2D and Supporting Fig. S4).

FIG. 2

TRIB1 KO 2D HLCs failed to mimic TRIB1 lipid biology versus 3D hepatic organoids. (A) Strategy for generation of TRIB1 KO hiPSCs by CRISPR‐Cas9 using ssODN knock‐in approach. A stop codon was introduced after eight amino acids from the start codon in exon 1. Hind III site was introduced to screen for ssODN‐edited clones. (B) Representative image shows WT and HindIII positive TRIB1 heterozygous and KO clones in undifferentiated hiPSCs. (C) Sequence confirmation shows inclusion of knocked‐in ssODN in genomic DNA and cDNA versus WT clone. (D) Schematic shows parallel differentiation of WT and TRIB1 KO hiPSCs to 2D and 3D hepatocytes, which were used for characterization. (E,F) TRIB1 KO 2D HLCs exhibit lower levels of hepatic markers than WT, whereas 3D hepatic organoids exhibit no defect in hepatic differentiation in WT and KO as shown by quantitative PCR (E) and immunostaining (F). (G) TRIB1 KO 2D HLCs showed decreased HNF4A+ (i) and ASGR1+ (ii) cells as quantified by immunofluorescence image quantification and FACS sorting, respectively. (H) TRIB1 KO 2D HLCs express low DNL (DGAT1, DGAT2, FASN, ACACA, and SCD) markers versus WT, whereas TRIB1 KO 3D hepatic organoids express higher DNL markers and TG (I), mimicking TRIB1 lipid biology. Values represent the mean ± SEM shown for one representative experiment, calculated relative to GAPDH (reference gene) using the 2−ddCt method, normalized to the mean expression level in WT cells (n = 4 WT and four KO wells; two wells per clone with 24 organoids per well. For the 2D HLCs, n = 3 WT and seven KO wells; and n = 1 WT and three KO clones). P values were calculated with Welch’s t tests. DAPI is shown in blue in all of the immunofluorescence images. Scale bar: 200 µm (Fi,ii), 75 µm (Fiii,iv), 125 µm (Fv,vi), and 275 µm (Fvii‐xii). Abbreviation: Het, heterozygous.

Both WT and KO hiPSCs were differentiated to hepatocytes in parallel from the same culture dish, following 2D and 3D protocols. To assess the differentiation efficiency, at day 20 cells were harvested for gene expression and lipid phenotyping. Significant decreases in the levels of hepatic markers such as HNF4A, AFP, ALB, and ASGR1 were observed in TRIB1 KO 2D HLCs when compared with WT 2D HLCs, whereas TRIB1 KO 3D hepatic organoids exhibited no such defect in differentiation compared with WT 3D organoids (Fig. 2E). Similar observations were made using immunofluorescence for hepatic markers (Fig. 2Fi‐vi,vii‐xii), which were further quantified by image quantification and FACS sorting for HNF4A+ and ASGR1+ cells (Fig. 2Gi‐ii).

Trib1 liver‐specific KO mice have increased expression of genes involved in lipid synthesis and increased liver triglyceride (TG) mass.( ) We characterized the expression of lipid synthesis genes such as diacylglycerol acyltransferase (DGAT)1, DGAT2, fatty acid synthase (FASN), acetyl‐coenzyme A carboxylase (ACACA), and stearoyl‐coenzyme A desaturase (SCD)1 in TRIB1 KO hiPSC‐derived hepatocytes. At day 20, expression of lipid synthesis genes was significantly reduced in TRIB1 KO 2D HLCs compared with WT cultures, consistent with the reduction in other canonical hepatocyte genes (Fig. 2H). Furthermore, the amount of TG in the TRIB1 KO 2D HLCs was markedly lower than in WT 2D cultures (Fig. 2I). In marked contrast, TRIB1 KO 3D organoids displayed significantly increased expression of lipid synthesis genes including DGAT1, DGAT2, and ACACA and a trend toward an increase in FASN and SCD compared with WT 3D organoids (Fig. 2H). The accumulated TGs were dramatically higher in the TRIB1 KO 3D organoids than in TRIB1 KO 2D HLCs and trended to being slightly higher compared with WT 3D organoids (Fig. 2I). Secreted TGs in the media of TRIB1 KO 3D organoids also trended higher than in the media of WT 3D organoids (Supporting Fig. S5).

TRIB1 KO 2D HLCs, but Not 3D Organoids, Exhibit Developmental Arrest Associated With Low Expression of C/EBPα

To determine how early TRIB1 affects hepatocyte differentiation in the 2D differentiation protocol, we measured gene expression at different stages of differentiation. At day 5, no change in the expression of DE stage marked by FOXA2, SRY (sex determining region Y)‐box (SOX17), and chemokine (C‐X‐C motif) receptor 4 (CXCR4) was observed (Fig. 3A). At day 12, a trend toward an increased expression of hepatoblast markers, such as hematopoietically expressed homeobox protein (hHEX) and delta‐like 1 homolog (DLK1), was observed in TRIB1 KO 2D HLCs over WT day 12 HLCs (Fig. 3B). However, at day 20, a significant reduction in hepatic markers like HNF4A to AFP, ALB, and ASGR1 was observed, and increased expression of hepatoblast markers (hHEX, DLK1, and PROX1) was observed in day 20 HLCs (Figs. 2E and 3B). No such increase in expression of hepatoblast markers was observed in 3D TRIB1 KO organoids at day 20 of differentiation (Fig. 3C).

FIG. 3

TRIB1 KO 2D HLCs exhibit developmental arrest due to the inherent low expression of C/EBPα. (A) TRIB1 KO 2D HLCs have equivalent DE markers (FOXA2, SOX17, and CXCR4) versus WT HLCs at day 5 of differentiation. (B) Increased expression of progenitor markers (hepatoblasts: hHEX, DLK1, and PROX1) were observed in TRIB1 KO HLCs at day 12 but predominantly at day 20 of differentiation versus WT HLCs. (C) Equivalent hepatoblast markers were observed in WT and TRIB1 KO 3D organoids at day 20 of differentiation. (D) Two‐dimensional HLCs have low expression of C/EBPα and C/EBPβ, which are further reduced in TRIB1 KO HLCs, whereas 3D organoids express adequate levels of C/EBPα and C/EBPβ, which increase in TRIB1 KO organoids as shown by quantitative PCR (D) and immunostaining for C/EBPα (E) (green) at day 20 of differentiation. Values represent the mean ± SEM. Data are representative of independent experiments, calculated relative to GAPDH (reference gene) using the 2−ddCt method, normalized to the mean expression level in WT cells (n = 4 WT and four KO wells; two wells per clone with 24 organoids per well. For the 2D HLCs, n = 8 WT and eight KO wells; and n = 2 WT and two KO clones). For C/EBPα and C/EBPβ quantitative PCR, normalized to the mean expression level in WT 2D HLCs. P values were calculated with Welch’s t tests. DAPI is shown in blue in all of the immunofluorescence images. Scale bar: 125 µm (Ei,ii) and 275 µm (Eiii,iv).

TRIB1 regulates levels of C/EBPα posttranslationally, and we reported increased levels of C/EBPα protein in Trib1 liver‐specific KO mice (LsKO).( ) TRIB1 KO 2D HLCs exhibited decreased C/EBPα messenger RNA (mRNA) (Fig. 3D) and protein (3.9%) (Fig. 3E and Supporting Fig. S6) compared with WT 2D HLCs (6.7%). In contrast, TRIB1 KO 3D organoids displayed dramatically increased C/EBPα protein (45.7%) (Fig. 3E and Supporting Fig. S6) compared with WT 3D organoids (20.4%), despite C/EBPα mRNA. Similar observations were made for C/EBPß mRNA (Fig. 3D). Together, these data suggest the developmental arrest of hepatoblasts to mature hepatocytes in TRIB1 KO 2D HLCs, but that this does not exist in TRIB1 KO 3D organoids, which recapitulate the expected phenotype of increased C/EBPα protein without an increase in C/EBPß mRNA.

Overexpression of C/EBPα Rescues Hepatic and TG Phenotypes of TRIB1 KO hiPSC Hepatocytes

To determine whether increased C/EBPα protein was sufficient to rescue the maturation of hepatocytes and in turn drive increased lipid synthesis, we overexpressed human C/EBPα in both WT and TRIB1 KO hiPSC‐derived hepatocytes using an inducible lentivirus system (Fig. 4A and Supporting Fig. S7). The hiPSCs were differentiated to 2D HLCs and 3D hepatic organoids in parallel. hC/EBPα expression was induced starting at day 13 of differentiation and continued for 7 days. At day 20, cells were harvested for gene expression and lipid phenotyping (Fig. 4A). As expected, hC/EBPα‐induced WT and TRIB1 KO 2D HLCs had increased expression of C/EBPα (Fig. 4B) and increased C/EBPα protein (Fig. 4C); notably, while the TRIB1 KO 2D HLCs without overexpression of C/EBPα had reduced C/EBPα protein compared with WT 2D HLCs, those with overexpression had markedly more C/EBPα protein than the overexpressing WT 2D HLCs (Fig. 4C). C/EBPα overexpression “rescued” expression of hepatic markers HNF4A, AFP, ALB and ASGR1, and decreased the expression of hepatoblast markers at day 20 of differentiation (Fig. 4D). The lipid synthesis genes showed significantly increased expression in the C/EBPα‐overexpressing WT and TRIB1 KO 2D HLCs (Fig. 4E). Consistent with this, TG synthesis and TG mass were both significantly increased in the C/EBPα‐overexpressing WT and TRIB1 KO 2D HLCs compared with control cells (Fig. 4F). Interestingly, both TG synthesis and TG mass were significantly higher in the C/EBPα‐overexpressing TRIB1 KO 2D HLCs compared with the C/EBPα‐overexpressing WT 2D HLCs, despite being lower in the control cells (Fig. 4F). These data suggest that the low inherent expression of C/EBPα in 2D HLCs restricts the formation of mature hepatocytes and can be overcome by forced C/EBPα overexpression.

FIG. 4

Overexpression of C/EBPα in WT and TRIB1 KO cells rescues the TRIB1 lipid phenotype in 2D HLCs. (A) Schematic showing experimental strategy in which WT and TRIB1 KO hiPSCs were exposed to empty vector and C/EBPα lentivirus. Differentiation was performed in parallel to 2D and 3D. At day 13, cells were induced with Dox for 7 days and at day 20, cells were harvested for downstream analysis. Shown here are the data for 2D HLCs; data for 3D are shown in Supporting Fig. S8. (B,C) C/EBPα was decreased in empty lentivirus‐transduced TRIB1 KO 2D HLCs versus WT HLCs, which were rescued both at the RNA and protein level on overexpression of C/EBPα. (D) Rescue of hepatic markers (HNF4A, AFP, ALB, and ASGR1) and decrease in progenitor markers (PROX1 and hHEX) were observed in day 20 WT and TRIB1 KO HLCs on overexpression of C/EBPα. Rescue of DNL gene expression (E) and TG phenotype (F) was observed in 2D TRIB1 KO HLCs on overexpressing C/EBPα. Values represent the mean ± SEM. Fold change calculated relative to GAPDH (reference gene) using the 2−ddCt method, normalized to the mean expression level in WT cells transduced with empty lenti vector as control (n = 6 WT + empty vector, four KO + empty vector; and n = 3 WT + C/EBPα and six KO + C/EBPα wells). P values were calculated with Welch’s t tests.

A similar observation was made with the 3D hepatic organoids, in which overexpression of hC/EBPα further increased the expression of hepatic (ALB) and lipid synthesis (DGAT1, DGAT2, FASN, ACACA, and SCD) genes as well as TG synthesis, and suppressed the expression of hepatoblast progenitor markers (hHEX and PROX1) at day 20 of differentiation (Supporting Fig. S8).

Discussion

Here we describe a 3D hiPSC‐derived hepatic organoid model that temporally recapitulates the stages of hepatocyte development more effectively than 2D hiPSC‐derived hepatocytes. We show that hiPSCs deleted in TRIB1 and differentiated to monolayer 2D HLCs are immature and have little of the expected lipid phenotype, but when differentiated to 3D, hiPSC‐derived hepatic organoids express more mature hepatocyte biology and recapitulate the expected biology of hepatocyte TRIB1 deficiency based on mouse studies, providing even further insights into the role of hepatocyte TRIB1 in lipid metabolism.

The early development and dynamic metabolic functions of a complex organ like the liver cannot be fully understood by analysis of a single type of cell in isolation, nor can the pathogenesis of genetic diseases affecting liver development and its function be fully recapitulated using individual cells.( , , ) Recently, liver organoids have been generated from a myriad of cells including progenitor cells obtained from human liver tissue,( ) hiPSCs,( , , ) a pool of hiPSC‐derived progenitors( ) with endothelial cells, mesenchymal stromal cells, human umbilical vein endothelial cells involving complex engineered matrices( ) or perfused systems,( ) or aggregation only in later stages of differentiation.( , ) We made the same observation that organoids have enhanced expression of hepatic markers and are more mature than 2D cultured cells, as reported previously.( , , ) Here we describe an efficient protocol using hiPSCs for hepatic 3D organoid differentiation by using short exposure to a nonengineered matrix. In contrast to the engineered matrices such as polyethylene glycol, polyglycolic/poly‐l‐lactic acid, glycosaminoglycans, pullulan, synthetic scaffolds with collagen, fibrinogen and laminin commonly used for 3D organoid studies,( , ) we used a short exposure to nonengineered laminin‐enriched matrigel and CultureX (CultureX, Cambridge, MA) as a matrix for hepatic organoid differentiation. Studies are ongoing to further characterize these 3D organoids and improve the strategies to better obtain defined populations of nonparenchymal cells. One needs further detailed analysis about possible central cell death, its relationship to organoid size, and the potential use of high oxygen conditions for application of this 3D approach in understanding the liver injury disease modeling.

We used this approach to explore the pleiotropic cardiometabolic gene TRIB1. GWAS has identified a locus at chromosome 8q24, near the TRIB1 gene, at which variants are associated with NAFLD, plasma lipids, adiponectin, and coronary artery disease.( ) Previously, we showed that targeted disruption of TRIB1 in mice results in an increased TG phenotype.( ) However, considering the appreciable differences in mice versus human lipid biology and the absence, thus far, of any individual with a homozygous KO mutation in the Tribbles locus, led us to generate TRIB1 KO hiPSC using genome editing to understand its role in human hepatic and lipid biology. This approach is far better and offers several advantages over using teratogenically transformed hepatoma cell lines and primary cells, which dedifferentiate in a culture dish, or tissues from diseased biopsies, which do not reflect the true metabolic state.( ) Moreover, of the other model systems discussed previously, only hiPSCs provide the opportunity to understand the role of a gene or rare mutation in early development.

To model TRIB1 lipid biology in a dish, WT and KO hiPSCs were differentiated to 2D HLCs, and surprisingly, TRIB1 KO HLCs showed reduced expression of hepatic markers and de novo lipogenesis (DNL) markers over WT HLCs in contrast to our earlier observations made with TRIB1 LsKO mice.( ) Interestingly, no defect in the induction of DE from hiPSCs was observed. However, expression of hepatoblast markers was observed to be maintained at a higher level in TRIB1 KO HLCs, even in later stages of differentiation when compared with WT HLCs. We also observed a decrease in the expression of C/EBPα in TRIB1 KO 2D HLCs over WT HLCs in contrast to TRIB1 LsKO mice.

Organoid systems have previously been reported to offer several advantages regarding maturation and function of iPSC‐derived cell types for several organ systems including hepatocytes.( , ) Taking advantage of our 3D hepatic organoid method, we differentiated WT and TRIB1 KO hiPSCs in parallel and observed no decrease in the hepatic markers in contrast to the 2D HLCs. However, a trend toward an increase in gene expression of hepatic markers as early as HNF4A was observed in TRIB1 KO organoids over WT organoids. Following this increase in C/EBPα, its downstream target lipid synthesis gene expression (DGAT1, DGAT2, FASN, ACACA, and SCD) and TGs were observed in TRIB1 KO organoids, mimicking the phenotype observed in TRIB1 LsKO mice.( ) Considering the shortcomings of conventional lipid measurement assays, stringent fluorescent‐based quantitative approaches are required to measure changes in miniscule lipid accumulation in small cultures like organoids. As reported previously in mice, an increase in the expression of genes responsible for fatty acid synthesis (FASN, ACACA, and SCD) and incorporation of di acyl to tri acyl glycerol (DGAT1 and DGAT2) were observed. This further reflects that TRIB1 is not responsible for hepatic specification, but is involved with hepatocyte maturation by affecting C/EBPα. This establishes a role of a metabolic gene, TRIB1, in hepatocyte development.

Monolayer 2D HLCs have been explored extensively for their robustness in understanding development and as an essential tool in drug discovery, but they lack the expression of critical late‐stage hepatic markers and are therefore considered to be metabolically immature in nature, either lacking or expressing suppressed levels of vital specification markers such as PROX1 and C/EBPα.( ) Leveraging our previously generated 2D HLC database,( ) consisting of approximately 90 HLCs of diverse background, we analyzed the gene expression of PROX1 and C/EBPα and found them to be expressed at significantly reduced levels when compared with other liver hepatic and metabolic genes such as ALB and APOC3 (Supporting Fig. S9). Interestingly, these HLCs also express low levels of peroxisome proliferator‐activated receptor (PPARα and PPAR delta) (Supporting Fig. S9), expressed robustly in metabolically mature hepatocytes. This highlights that although 2D cultures exhibit robust expression of hepatocyte markers and are a good model for some studies, they still lack the expression of key hepatocyte markers. Therefore, the need for 3D organoid cultures prevails to overcome this limitation and to produce metabolically mature hepatocyte‐like cells.

TRIB1 is reported to posttranscriptionally regulate levels of C/EBPα by ubiquitination.( , ) It has been observed that 2D HLCs express only low levels of C/EBPα, and in the absence of TRIB1, there is a further decrease in C/EBPα RNA. This affects the overall protein turnover of C/EBPα, which is reflected as reduced levels of the C/EBPα protein. This decrease in C/EBPα resulted in maturation defects of TRIB1 KO 2D HLCs, which fail to up‐regulate late‐stage hepatic markers and maintain the hepatoblast stage markers, as previously shown in primary hepatocytes( ) and C/EBPα KO mice.( ) Along similar lines, overexpression of C/EBPα rescued the maturation defects in TRIB1 KO 2D HLCs, further emphasizing that the TRIB1‐C/EBPα relationship holds true in hiPSC hepatocytes, which was responsible for influencing TG levels.( ) However, overexpression of C/EBPα in WT HLCs also increased hepatic and DNL gene expression, indicating the attainment of a more mature state, emphasizing the definite need of C/EBPα for maturation of hepatocytes.( , , )

Previous reports have established that C/EBPα regulates HNF4A, and both C/EBPα and HNF4A co‐operatively act to regulate downstream genes.( , ) Low levels of C/EBPα in 2D HLCs, in turn, resulted in decreased expression of HNF4A, which further explains the effect on hepatocyte differentiation. Recently, primary hepatocytes, in which TRIB1 was knocked down using small interfering RNA,( ) showed a decrease in expression of HNF4A, which further mechanistically explains the influence on hepatocyte differentiation in 2D HLCs at the level of HNF4A in TRIB1 KO 2D HLCs.

Supporting our observation in hiPSC hepatic organoids, recent reports in HepG2 and rat primary hepatocytes have shown that hepatic genes such as PROX1, HNF4α, and C/EBPα are induced in 3D hepatocyte cultures.( , ) TRIB1 liver‐specific KO mice has increased C/EBPα protein (due to reduced protein degradation) and decreased C/EBPα transcript (due to negative feedback from the increased protein).( ) We also know that TRIB1 is a transcriptional target of C/EBPα, leading to a negative feedback loop in which excess C/EBPα induces TRIB1, which then promotes the degradation of C/EBPα. In 2D cultures, we suggest that poor expression of C/EBPα fails to effectively induce TRIB1 expression, whereas consistent with the recent reports, we suggest that C/EBPα is more efficiently induced in 3D organoids, which more efficiently activates TRIB1 expression. This led us to speculate that in 3D organoids, enhanced expression of hepatic markers—especially C/EBPα—is maintained, which allowed the Trib1‐ C/EBPα feedback loop to establish and recapitulate the phenotype.

In summary, we report a simplified and defined 3D iPSC‐hepatocyte organoid system. We show, using TRIB1‐deleted iPSCs as proof of principle, that this 3D system more faithfully recapitulates the expected biology of TRIB1 deletion in hepatocytes than the 2D system. In this model, the effects of TRIB1 KO are more clearly seen in the 3D organoids than in the 2D cultures. This study establishes the importance of having a robust 3D model system to understand the role of genes, which are poorly reflected by conventional monolayer cultures. Furthermore, this organoid approach has the potential to permit expansion of our understanding of genetic diseases affecting the liver and has substantial implications for drug discovery.

Materials and Methods

hiPSC Culture

Established hiPSC lines were maintained in feeder‐free culture conditions in either mTeSR1 (STEMCELL Technologies) or iPS‐Brew (Miltenyi Biotec, Bergisch Gladbach, Germany) on Geltrex coated dishes (Thermo Fisher Scientific, Waltham, MA) as described previously.( ) In‐house‐derived hiPSC line Sv20 derived from peripheral blood mononuclear cells and the K3 hiPSC line produced from foreskin fibroblasts were used for the studies, which were characterized previously in detail.( , )

Hepatocyte Differentiation

2D Monolayer

hiPSC lines were passaged with Accutase (Innovative Cell Technologies, San Diego, CA) for 5 minutes at room temperature (RT) and plated for differentiation on Geltrex coated dishes following Stephen Duncan’s protocol with modifications.( ) For efficient DE differentiation, the DE kit (STEMCELL Technologies) was used for the first 4 days in ambient O2/5% CO2. Later, cells were cultured in RPMI‐B27 (with insulin) supplemented with BMP4 (20 ng/mL; PeproTech, Rocky Hill, NJ) and FGF2 (10 ng/mL; R&D Systems) for 5 days in 5% O2/5% CO2, followed by RPMI‐B‐27 (with insulin) supplemented with recombinant human HGF (20 ng/mL; PeproTech) for 5 days in 5% O2/5% CO2. The last stage required culturing cells in HCM Hepatocyte Culture Medium (Lonza Group AG, Basel, Switzerland) that was supplemented with recombinant human Oncostatin M (20 ng/mL; R&D Systems) for 5 days in ambient O2/5% CO2.( )

3D Organoid

hiPSC lines were passaged with Accutase for 5 minutes at RT and plated in 96‐well U‐bottom low adherent plates using the DE kit with 10 µM Rock inhibitor and continued in ambient O2/5% CO2 throughout the differentiation. Approximately 3,000‐5,000 cells per 100 µL per 96‐well plate gave efficient aggregation to organoids with no to minimal evidence of central darkening suggestive of cell death. However, standardization from line to line is necessary to make sure the organoid size does not become limiting, in which case the potential for central cell death is higher. After 24 hours, half the medium was removed, and fresh medium with 4% matrigel and CultureX was added and cultured for 2 more days with daily medium addition. At day 5, aggregates were harvested and cultured in low‐adherent 6‐well plates in RPMI‐B27 (with insulin) supplemented with BMP4 (20 ng/mL), FGF2 (10 ng/mL), and KGF (20 ng/mL; PeproTech) for 5 days, followed by RPMI‐B‐27 (with insulin) supplemented with recombinant human HGF (20 ng/mL; PeproTech) and 0.5% DMSO for 5 days. Later, aggregates were cultured in HCM Hepatocyte Culture Medium and supplemented with recombinant human Oncostatin M (20 ng/mL) and Dexamethasone for 5 days.

CRISPR‐Mediated Genome Engineering of hiPSCs

Guide RNAs were designed and evaluated for potential off‐target activity using the online CRISPR Design tool at http://crispr.mit.edu. Genome editing was performed using SpCas9_GFP (Addgene plasmid #48138), which co‐expresses a human codon–optimized Cas9 nuclease and green fluorescent protein (GFP) using a viral 2A sequence. To generate the knockout using a knock‐in approach, a 200‐bp ssODN was designed in Exon 1 with a stop codon introduced after eight amino acids and a HindIII site, to confirm the targeted clones.

For nucleofection, cells in a 70%‐80% confluent 10‐cm plate were dissociated into single cells with Accutase, resuspended in phosphate‐buffered saline, and combined with 5 µg SpCas9‐GFP and 1 µg of ssODN C*G*C*A*GCGTCTTCCCGCGCGGATCCCGGGACTTAAAAAGCCGGGGCCACCCCGGCCCAGGACGGGATGCGGGTCGGTCCGGTGCGCTCTTAGTGATAAGCTTCTCGCAGCCCCGCGGCCCGGCCCTGCTCTTCCCAGCCACCCGAGGCGTCCCGGCCAAACGCCTGCTGGACGCCGACGACGCGGC*G*G*C*T. A single pulse was delivered at 250 V/500 mF (Lonza), and the cells were recovered and plated in mTeSR1 with 10 μM ROCK inhibitor (Y‐27632; Cayman Chemical Co., Ann Arbor, MI). Cells were dissociated with Accutase 48 hours after electroporation, and GFP‐positive cells were isolated by FACS (FACS AriaII; BD Biosciences, San Jose, CA) and replated onto 10‐cm Geltrex‐coated plates (~10,000 cells/plate) with conditioned medium and 10 µM ROCK inhibitor to facilitate recovery. Single cell–derived clones were picked after 2 weeks, expanded, and evaluated for presence of HindIII site by restriction digestion after polymerase chain reaction (PCR) amplification of a small region surrounding the targeting site, using primers F‐5′AAGACCAGTCTGCAAACTCCAT3′ and R‐ 5′ CTGAGGAAAGGAACAGATAGAGAAAG 3′. Clones were further confirmed with Sanger sequencing of the PCR products, before selecting them for differentiation.

Lentiviral Mediated Overexpression of hC/EBPα

To generate the lentivirus construct, hC/EBPα cDNA was purchased from the DNASU Plasmid Repository (Tempe, AZ) and cloned into pEn_TTmcs vector with Acc65I and XhoI, and gel purified the 3.8‐kb vector band. C/EBPα was ligated as insert and transformed into DH5α followed by digestion with XbaI and NcoI. To generate the empty pEn_TTmcs control vector, the XbaI‐digested vector was recircularized. Vectors were confirmed with Sanger sequencing and transferred to a pSLIK‐puro lenti backbone using Gateway LR clonase (Thermo Fisher Scientific).

Human embryonic kidney cells were seeded in 10‐cm lysine‐coated dishes and transfected using Lipofectamine 3000 with 12 µg lenti backbone, 4.5 µg PMD2, and 7.7 µg pSPAX2 in the presence of Polybrene (Sigma Aldrich, St. Louis, MO) according to the manufacturer’s protocol. Media were collected every 24 hours following transfection for 3 consecutive days, filtered with 0.45‐µ filters and concentrated 10 times, then added to hiPSCs. After 24 hours the regular iPS‐Brew medium was changed to hiPSCs and selected with Puromycin. Puromycin‐resistant cells were passaged with Accutase for 5 minutes to single cells and replated on 10‐cm Geltrex dishes for expansion. Single cell–derived clones were picked manually and cultured in the presence and absence of Doxycycline (Dox). hiPSCs were then analyzed for induction of C/EBPα by gene‐expression analysis.

Quantitative Real‐Time PCR

RNA isolation was performed using Qiagen miRNeasy and Zymo kits (Hilden, Germany) according to the manufacturer’s instructions, and cDNA was prepared using a High Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific) with an equimolar mixture of random hexamers and oligo‐dT. Gene expression was measured using the TaqMan Gene Expression Assays (Thermo Fisher Scientific) in 10 µL quantitative real‐time PCR reaction with 20 ng cDNA performed in technical duplicate. Reactions were carried out on a 7900HT QuantStudio Fast Real‐Time PCR System (Thermo Fisher Scientific). Relative expression levels were quantified by the 2–DDCt method using glyceraldehyde 3‐phosphate dehydrogenase (GAPDH) as reference gene.

Immunofluorescence Staining

Cells were fixed for 15 minutes at RT with 4% paraformaldehyde (PFA) and incubated overnight with antibodies anti‐FOXA2 (sc‐6554), anti‐GATA4 (sc‐1237), anti‐HNF4A (sc‐6556), anti‐AFP (A‐8452), anti‐albumin (CL‐2513A), anti‐C/EBPα (8178S), anti‐MDR1, anti‐arginase, anti‐HNF1ß (12533‐1‐AP), and anti‐ZO‐1 (40‐2200) as described previously.( )

Flow Cytometry

Day 20 hepatocytes were dissociated using enzyme mixture (Accutase:collagenase:trypsin in 1: 1:1), washed with Permeabilization Wash Buffer (BioLegend, San Diego, CA) and stained for antibodies against asialoglycoprotein receptor 1 (563655) for 1 hour in the dark. Flow cytometry was performed using an Accuri Flow Cytometer (BD Biosciences) with unstained cells and only secondary stained as controls.

Bile Acid and LDL Uptake

The organoids were washed with Dulbecco’s phosphate buffered saline (DPBS), then incubated with 5 μg/mL CDFDA and 5 μg/mL Dil‐Ac‐LDL (Invitrogen) following the manufacturer’s instructions, and washed again before examination with fluorescence microscopy. Later, organoids were fixed with 4% PFA and mounted in optimal cutting temperature for sectioning.

PAS Staining

The organoids were fixed with 4% PFA in DPBS and the sections were stained with a PAS staining kit from Sigma‐Aldrich according to the manufacturer’s instructions.

ORO Staining

For detecting accumulated lipids, organoids were fixed with 4% PFA in DPBS and the sections were treated with 60% isopropanol for 5 minutes. Then the isopropanol was removed and ORO (Sigma‐Aldrich) was added following the manufacturer’s instructions and incubated for 15 minutes at RT. ORO was then removed and sections were rinsed with water until clear before imaging.

Enzyme‐Linked Immunosorbent Assay for Secretory Proteins

The human albumin, APOB, and APOC3 content in the cultured media was measured using enzyme‐linked immunosorbent assay kits from Bethyl Laboratories (Montgomery, TX), MABtech AB (Nacka Strand, Sweden), and Thermo Fisher Scientific, respectively, following the manufacturers’ instructions. The secretion was normalized to total cell RNA isolated from Zymo kits.

TG Measurement

Total TG was measured by colorimetric enzymatic assay. Briefly, lipids were extracted from day 20 HLCs and organoids using (3:2) hexane:isopropanol, and pellet was resuspended in water after 15% Triton‐chloroform treatment. Absorbance was measured at 500 nm after incubating 10 µL of extracted lipid or 10 µL of media with 190 µL TG Infinity enzymatic solution for 15 minutes at 37°C (Thermo Fisher Scientific). To measure newly synthesized TGs, day 20 cells were pre‐incubated with [3H]‐1,2,3‐Glycerol (10 µCi/mL) in the presence of Oleic Acid (0.4 mM) for 4 hours. To measure radiolabeled TG, lipids were extracted from cells, fractionated by thin‐layer chromatography, and quantified by beta counter (Beckman Coulter, Brea, CA). TG counts in cells was normalized to the total cellular protein for 2D HLCs, as measured by BCA. For organoids, equal number and sized organoids per treatment were used for the assay.( )

Statistics

All of the quantitative data represent mean + SEM. The average gene‐expression levels or changes in lipid levels were compared among groups using the Welch’s t test, due to the groups having unequal sizes. Graphs were plotted using GraphPad Prism 8 for Mac OS X (GraphPad Software, San Diego, CA).

Fig S1

Click here for additional data file.

Fig S2

Click here for additional data file.

Fig S3

Click here for additional data file.

Fig S4

Click here for additional data file.

Fig S5

Click here for additional data file.

Fig S6

Click here for additional data file.

Fig S7

Click here for additional data file.

Fig S8

Click here for additional data file.

Fig S9

Click here for additional data file.