# Biophysical Control of Bile Duct Epithelial Morphogenesis in Natural and Synthetic Scaffolds

Anette Funfak, Latifa Bouzhir, Emilie Gontran, Nicolas Minier, Pascale

Dupuis-Williams, Samy Gobaa

► To cite this version:

Anette Funfak, Latifa Bouzhir, Emilie Gontran, Nicolas Minier, Pascale Dupuis-Williams, et al.. Bio-physical Control of Bile Duct Epithelial Morphogenesis in Natural and Synthetic Scaffolds. Frontiers in Bioengineering and Biotechnology, 2019, 7, pp.417. 10.3389/fbioe.2019.00417. pasteur-02936645

HAL Id: pasteur-02936645

https://pasteur.hal.science/pasteur-02936645

Submitted on 11 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Biophysical Control of Bile Duct Epithelial Morphogenesis in Natural and Synthetic Scaffolds

Anette Funfak¹, Latifa Bouzhir², Emilie Gontran², Nicolas Minier¹,³,
Pascale Dupuis-Williams²,⁴*† and Samy Gobaa¹*†

¹Institut Pasteur, Biomaterials and Microfluidics Core Facility, Paris, France, ²Université Paris-Saclay, UMR-S1174 INSERM, Orsay, France, ³Université de Technologie de Compiègne, Alliance Sorbonne Université, Compiègne, France, ⁴ESPCI, PSL University, Paris, France

**OPEN ACCESS**

**Edited by:**
Elena Martinez,
Institute for Bioengineering of
Catalonia (IBEC), Spain

**Reviewed by:**
Sharon DeMorrow,
University of Texas at Austin,
United States
Romina Mancinelli,
Sapienza University of Rome, Italy

*Correspondence:
Pascale Dupuis-Williams
pascale.dupuis-williams@u-psud.fr
Samy Gobaa
samy.gobaa@pasteur.fr

†These authors have contributed
equally to this work

**Specialty section:**
This article was submitted to
Tissue Engineering and Regenerative
Medicine,
a section of the journal
Frontiers in Bioengineering and
Biotechnology

**Received:** 02 October 2019
**Accepted:** 29 November 2019
**Published:** 13 December 2019

**Citation:**
Funfak A, Bouzhir L, Gontran E, Minier N, Dupuis-Williams P and Gobaa S (2019) Biophysical Control of Bile Duct Epithelial Morphogenesis in Natural and Synthetic Scaffolds. *Front. Bioeng. Biotechnol.* 7:417. doi: 10.3389/fbioe.2019.00417

The integration of bile duct epithelial cells (cholangiocytes) in artificial liver culture systems is important in order to generate more physiologically relevant liver models. Understanding the role of the cellular microenvironment on differentiation, physiology, and organogenesis of cholangiocytes into functional biliary tubes is essential for the development of new liver therapies, notably in the field of cholangiophaties. In this study, we investigated the role of natural or synthetic scaffolds on cholangiocytes cyst growth, lumen formation and polarization. We demonstrated that cholangiocyte cyst formation efficiency can be similar between natural and synthetic matrices provided that the mechanical properties of the hydrogels are matched. When using synthetic matrices, we also tried to understand the impact of elasticity, matrix metalloprotease-mediated degradation and integrin ligand density on cyst morphogenesis. We demonstrated that hydrogel stiffness regulates cyst formation. We found that controlling integrin ligand density was key in the establishment of large polarized cysts of cholangiocytes. The mechanism of lumen formation was found to rely on cell self-organization and proliferation. The formed cholangiocyte organoids showed a good MDR1 (multi drug resistance protein) transport activity. Our study highlights the advantages of fully synthetic scaffold as a tool to develop bile duct models.

**Keywords:** bile duct, microenvironment, organoids, hydrogels, 3D cell culture

## INTRODUCTION

Cellular microenvironment cues, including surface topography, substrate rigidity and biochemical signals determine to a large extent the outcome of many biological processes (Tan et al., 2003; Gobaa et al., 2011; Unadkat et al., 2011). Strikingly, the exposure of stem cells to mechanical cues allows steering the differentiation process (Engler et al., 2006). Gaining control over the mechanical stress landscape of the cell was cornerstone in the development of successful organ-on-chip and strategies (Huh et al., 2013; Benam et al., 2016) and of advanced 3D culture systems (Gjorevski et al., 2016).

Recapitulating a physiological hepatic function *in vitro* is currently a very active line of research. Approaches, based on the use of primary adult hepatic cells (Broutier et al., 2016) or on the controlled differentiation of induced Pluripotent Stem Cells (iPSC) (Takebe et al., 2014) have led to the derivation of new organotypic models. These advanced cell culture systems were developed in order to improve on the physiological relevance of standard 2D culture hepatic cell lines.
Despite great successes (Broutier et al., 2016) in recapitulating the hepatic function *in vitro*, the inclusion of a biliary tree capable of removing the toxic metabolite is yet to be developed. The biliary tree is composed of intra and extra hepatic ducts organized in a complex network of interconnected tubes with luminal diameters ranging from <15 µm to a few mm capable of carrying the bile from the liver to the intestine (Roskams et al., 2004; Boyer, 2013). The bile ducts are lined with the epithelial cholangiocytes which regulate, through secretion and adsorption, the composition and flow of the bile (Bogert and LaRusso, 2007).

The biliary system is itself the target of a vast array of liver pathologies called cholangiopathies (Cheung et al., 2018) that account for a large proportion of liver transplants (Murray et al., 2005). Today's treatments of cholangiopathies consist mostly in either largely ineffective pharmacotherapies or in surgeries that are restricted to the extrahepatic large duct. Living donors of liver grafts often suffer from biliary complications like stricture formation, bile leaks and infection (De Assuncao et al., 2017).

The field of bile duct engineering, previously hindered by the lack of available cells, has recently benefited from progress in the understanding of biliary physiopathology and the availability of cholangiocytes cell lines or methods recapitulating differentiation of hepatic precursor cells, stem cells or iPSCs, using various mixtures of growth factors known to direct bile duct embryogenesis (Dianat et al., 2014; De Assuncao et al., 2015; Ogawa et al., 2015; Sampaziotis et al., 2015; Takayama et al., 2016). What remains to be explored when considering developing innovative liver therapies is to understand how the microenvironment guides the self-organization of cholangiocytes into functional biliary tubes (Shiojiri and Sugiyama, 2004; Raynaud et al., 2011).

Cell culture in 3D has shown that cholangiocytes possess *in vitro* morphogenetic capabilities in natural hydrogels like collagen or Matrigel (Ishida et al., 2001; Tanimizu et al., 2007; Hashimoto et al., 2008; Kido et al., 2015). The scaffolding properties of Matrigel help supporting the 3D reorganization of adult cells. The production of “cholangioids” in Matrigel from healthy cholangiocytes versus issued from patients with primary sclerosing cholangitis (PSC) (Soroka et al., 2019) has recently exemplified their interest for the study of pathogenicity of biliary diseases and their putative use in the identification of therapeutic targets (Loarca et al., 2017). Matrigel embedding is also used as a final step in most cholangiocytes differentiation protocols, to assess their capacity to organize in biliary cysts (Dianat et al., 2014; De Assuncao et al., 2015; Ogawa et al., 2015; Sampaziotis et al., 2015; Takayama et al., 2016). However, relaying on hydrogels based on natural extracellular matrix (ECM) clearly hinders the mechanistic understanding of organogenesis or cysts formation as these scaffolds prohibits the decoupling of the biophysical and biochemical signaling. Similarly, liver stiffness is a key marker of liver pathology (Wells, 2008; Mueller and Sandrin, 2010) and marks the development of chronic fibrotic diseases (Saneyasu et al., 2016). In this context using liver organoids to study disease progression begs for the development of hydrogel systems with tunable mechanical properties. Interestingly, the stiffness-dependent differentiation of cholangiocytes seems to be modulated by the biochemical

nature of the used matrix (Kourouklis et al., 2016). This indicates the high level of intricacy of both the biochemical and biophysical cues for these finely regulated biological processes.

This complex, and sometime undefined, nature of ECM-based hydrogels including Matrigel (MT) hinders the mechanistic understanding of organogenesis in a dish. On the other hand, synthetic hydrogels technology, such as systems based on the polymerization of polyethylene glycol (PEG), underwent many developments in the last decades placing them at the top of the list of the most physiologically relevant yet fully defined substrates for cell culture (Gjorevski et al., 2016). Today, synthetic hydrogel systems gained key functionalities including physiological elastic moduli (Ranga et al., 2014), easy tethering of biochemical ligands (Mosiewicz et al., 2013) and matrix metalloprotease (MMP) mediated proteolytic degradation allowing cell migration (Lutolf and Hubbell, 2003). These engineered functionalities opened the door for the establishment of advanced and physiologically relevant assays in fully synthetic scaffolds including the production of intestinal organoids (Gjorevski et al., 2016), the differentiation of embryoid bodies into neural tubes (Ranga et al., 2016) and the production of epithelial cysts (Enemchukwu et al., 2016). Therefore, the side-by-side comparison of organotypic development in both synthetic (PEG) and natural (MT) hydrogels could help with defining the set of microenvironment features that are essential to the proper self-organization of the seeded cells.

Cholangiocytes, like most other epithelial cells are able, when embedded in natural ECM like Matrigel or collagen, to self-organize into polarized monolayers enclosing a central lumen termed cyst. These structures are good models for studying duct morphogenesis since they share with the tubes the same topology and the same organization of the epithelium, where the polarized cells are connected by adherens and tight junctions, thus ensuring secretion, cohesiveness, and coordinated growth of the monolayer. Here we report on the quantitative impacts of biochemical and biophysical cues on the development of cholangiocyte cysts *in vitro*. We demonstrate that cholangiocytes can organize and produce polarized cysts in both natural and synthetic scaffolds. We also show that the biophysical and biochemical parameters of the extracellular milieu are key for cyst morphogenesis. We also report on the dynamics of lumen formation and how this process is influenced by the extracellular milieu.

# MATERIALS AND METHODS

## Cell Culture

Normal Rat Cholangiocytes (NRC) cell line was obtained from N. LaRusso's laboratory (Vroman and LaRusso, 1996). The cells were cultured on collagen coated T25 flasks and maintained in DMEM F12 medium supplemented with 5% fetal bovine serum (Life Technologies), 1× Antibiotic-Antimycotic (Life Technologies) and growth active factors including Insulin (Life Technologies), Dexomethosane (Sigma), 3,3,5-Trimethoxy-L-thyronine sodium salt (Sigma), Bovine Pituitary Extract (BPE, Life Technologies) were added as described in de Groen et al. (1998).
## Hydrogel Formation

Synthetic hydrogels were formed by Michael-type addition of bicystein-bearing peptides onto star-shaped, vinylsulfone-functionalized, polyethylene glycol macromeres. Lyophilized peptides were ordered from Pepmic (China). Peptide sequences 16A (Ac-GCRD-GPQG↓1AGQ-DRCG-NH2) and 16R (Ac-GCRD-VPMS↓MRGG-DRCG-NH2) were chosen in order to ensure varying MMP sensitivity as established elsewhere (Lutolf and Hubbell, 2003; Patterson and Hubbell, 2010). An RGD peptide (Ac-GCGYGRGDSPG-NH2) was used at various concentrations in order to promote the integrin mediated adhesion of encapsulated cells. Briefly, peptides were resuspended in Milli-Q water to obtain a 12% (w/v) working solution. Before using these peptides to form hydrogel, we controlled their integrity and purity on a Bruker UltrafeXtreme MALDI-TOF/TOF instrumentation (Bruker-Daltonics, Germany). One microliter of peptide at 0.1 mg/ml was deposited on a MTP 384 ground steel target plate with 1 µl of 2,5-Dihydroxybenzoic acid (2,5-DHB) in 50% acetonitrile, 0.1% trifluoroacetic acid as matrix solution. Data were acquired using Flexcontrol software (Bruker-Daltonics, Germany) and shots were recorded in positive ion reflectron mode. Mass spectra were externally calibrated in the m/z range of 700–3,500 Da with a peptide calibration standard (Bruker-Daltonics, Germany) and analyzed with the Flexanalysis software (Bruker). Next, we performed an Ellman assay as per the manufacturer instructions in order to quantify free thiol concentration in each peptide working solution. The difference between theoretical and measured thiol concentration was taken in consideration in order to produce hydrogels with the desired Thiol (Th)/Vinyl Sulfone (VS) stoichiometry. Finally, we measured the peptide net content by hydrolyzing the samples in 1% phenol 6 N HCL for 20 and 48 h at 110°C in presence of known amount of NorLeucine as internal standard. After HCL evaporation the samples are analyzed on a L-8800 Hitachi amino acids analyzer (post-column derivatization with ninhydrin after ion-exchange chromatography separation).

Eight-arm polyethylene glycol-vinyl sulfone or PEG-VS 40 kDa macromeres were ordered as a custom synthesis from NOF (Japan). Lyophilized powder was dissolved in a 0.3 M, pH 7.5, HEPES buffer (Dominique Dutscher, France) at 12% (w/v). Hydrogel sterility was ensured by filtering both peptide and PEG-VS working solutions on a 0.22 µm syringe filter (Millex-GS, Merck). Polyethylene glycol (PEG) hydrogels were produced as described in Lutolf and Hubbell (2003). Rapidly, 4.5% (w/v) containing 100 µM RGD gels were formed by mixing 63 µl PEG-VS macromere (12% w/v) with 103 µl HEPES buffer, 2 µl RGD peptide (10 mM) and 40 µl NRC cells resuspended in HEPES buffer with a cell concentration of $1.5 \times 10^6$ cells/ml. In a second step we added 12 µl of 16R or 16A peptide working solutions (12% w/v). Before gelation, the mixture was pipetted up and down gently. Fifty microliters drops of gel mixture were sandwiched between two mirror-polished PTFE pieces separated with a 1 mm spacer in order to produce 50 µl hydrogel disks. The assemblage was then placed in a Heracell 150i incubator (ThermoFisher Scientific, Waltham, MA USA) at 37°C, 100% relative humidity for 30 min. Upon crosslinking, gel disks were

demolded and transferred to a 12 well plate (Falcon® multiwell plate, Corning) loaded with the adequate cell culture medium.

For Matrigel experiments, NRC cells were first mixed with a 7 mg/ml growth factor reduced Matrigel (Corning) solution and then diluted to 3.5 mg/ml with pre-cooled NRC complete medium on ice with a cell concentration of $3 \times 10^5$ cells/ml. Before the filling of 8-well µ-slides with the Matrigel-cell mixture, the pre-cooled slides were coated with 50 µl 3.5 mg/ml Matrigel solution and polymerized for 15 min at 37°C and 5% CO$_2$. Subsequently 300 µl of Matrigel-cell mixture was pipetted in each well and polymerized for 15 min at 37°C and 5% CO$_2$ before adding cell culture medium.

For both hydrogel systems NRC cells were pre-dissociated by up and down pipetting of the cell suspension followed by a post-dissociation step using a 40 µm cell strainer (Merck). The dissociated cell suspension was then diluted to the aimed cell concentration and mixed with the hydrogel solutions. It is to mention, that the natural and synthetic hydrogel system required a different state of cell dissociation (low for PEG hydrogel and high for Matrigel) due to the system requirements and the inhibition of NRC cell migration in PEG hydrogels. This was achieved by the adjustment of the pre-dissociation step.

## Rheology

The hydrogel shear modulus was obtained by rheometric measurement with a Kinexus Ultra Plus Rheometer (Malvern Panalytical, UK). Briefly, gels of different concentration were casted as cylinders of 4 mm in diameter after swelling in water overnight. The rheometer was fitted with a 4 mm geometry. Each gel disk was compressed to 80% of its nominal height in order to avoid slippage during the measurements. An amplitude sweep was performed to determine the linear viscoelasticity region (LVER) of the measured material. Shear moduli were extracted from a frequency sweep analysis performed with a constant 1% strain. The reported G' values were obtained at 0.5 Hz.

The hydrogel shear modulus for Matrigel was obtained by a direct polymerization of the polymer solution on the working stage of the rheometer. Therefore, samples of 314 µl were added to the center of the pre-cooled (4°C) lower working plate. The geometry of 20 mm was then immediately lowered before the gel started to form, to a working gap of 1 mm. For the polymerization of the polymer suspension, the temperature was increased stepwise with a rate of 5°C/min from 4 to 37°C and kept at 37°C during measurement. A solvent trap cover was used to prevent evaporation effects around the working stage. Rheometric measurements were performed at 0.5 Hz and a constant 0.1% strain.

## Viability/Cytotoxicity Assay

PEG and Matrigel hydrogels were washed with PBS and then incubated for 30 min at 37°C, 5% CO$_2$ with 4 µM Ethidium homodimer-1 and 2 µM Calcein, AM (ThermoFisher scientific) in growth medium. After a second washing step, the gels were imaged as described below. In average 12 images were acquired per condition. Fluorescence intensity was measured after background correction and normalized for the variations in
cyst surface using ImageJ, an image analysis software available on the NCBI website (https://imagej.nih.gov/ij/).

## Functionality Assays

To demonstrate the multidrug resistance protein 1 (MDR1) transporter activity of single or multi-lumen cysts, polyethylene glycol (PEG) and Matrigel hydrogels were treated with 50 µM rhodamine 123 (Rh123, Sigma) in serum free medium for 1 h at 37°C. After a washing step (3 × in serum free medium) the hydrogels were incubated for 2 h at 37°C, 5% CO₂ in fresh complete medium. To inhibit the MDR1 transporter activity the cysts containing hydrogels were incubated with 50 µM verapamil (Sigma-Aldrich) at 37°C for 30 min before adding Rh123. The rhodamine assay was then repeated as described above. For the quantification of the Rh123 assay the number of cysts containing Rh123 in the luminal space and the number of cysts showing Rh123 signal only at the level of the epithelium (blocked cysts) were counted with ImageJ.

## Immunostaining of the Cysts

Immunostaining were performed without extracting the NRC cysts from PEG hydrogel or Matrigel. The PEG hydrogels were washed with PBS and then fixed with 4% paraformaldehyde (PFA) overnight at 4°C. Cysts were permeabilized with 0.5% Triton X-100 in PBS for 1–2 h at 4°C and then blocked with 2% BSA, 0.1% Tween-20 in PBS for 3 h. The primary polyclonal anti-E-cadherin antibody (Thermo Fisher Scientific PA5-32178, 1:400) was diluted in blocking solution and incubated for 2–4 days at 4°C. Unbound primary antibodies were washed by incubating the gels in blocking buffer (3 times, 10 min each). The secondary antibody Goat anti-Rabbit IgG (H+L), Alexa Fluor Plus 647 (Thermo Fisher Scientific A32733, 1:400) was diluted in blocking solution and incubated for 2–4 days at 4°C. Nuclei and f-actin filaments were stained by adding 10 µM Hoechst and 16.5 nM Phalloidin 488 (Thermo Fisher Scientific) in PBS containing 1% BSA. Gels were incubated for 2 days at 4°C. Stained gels were then washed intensively in PBS.

Matrigel hydrogels were washed with PBS and fixed with 4% paraformaldehyde and 5% sucrose in PBS. Cyst were permeabilized with 0.5% Triton X-100 in PBS for 30 min and then blocked with 0.1% BSA and 1% goat serum for 30 min at RT. The primary polyclonal anti-E-cadherin antibody (Thermo Fisher Scientific PA5-32178, 1:400) and 16.5 nM Phalloidin 488 (Thermo Fisher Scientific) were diluted in blocking solution and incubated overnight at 4°C. The Matrigel gels were washed by incubation in PBS containing 0.05% Tween (3 times, 10 min each). The secondary antibody Goat anti-Rabbit IgG (H+L), Alexa Fluor Plus 647 (Thermo Fisher Scientific A32733, 1:400) was diluted in blocking solution and incubated for 1 h at RT. After washing, the gels were mounted with Prolong Gold antifade reagent (Invitrogen Molecular probes) containing DAPI for nuclei staining.

## Imaging

Wide field microscopy was performed on an inverted axio Observer Z1 inverted microscope from Zeiss (Carl Zeiss Microscopy GmbH, Jena, Germany) equipped with an OrcaFlash

4 V2.0 camera (Hamamatsu, Hamamatsu city, Japan) and Plan-Apochromat phase 1 10× (NA, 0.45) or Plan-Apochromat phase 2 20× (NA, 0.5) objectives. In case of live imaging (cyst growth in PEG hydrogel) temperature and CO₂ were ensured with a PeCon incubator equipped with the corresponding CO₂ and temperature modules.

Confocal imaging was performed on an SP8 inverted microscope (Leica Microsystems, Germany) equipped with HC FLUOTAR L 25×/0.95 W 0.17 VISIR objective and a white laser system. Confocal imaging of cyst polarization in Matrigel was acquired with a confocal microscope Nikon Eclipse TE-2000-E equipped with a X20/0.45 Plan Fluor ELWD Ph1 DM objective (∞/0–2 WD 7.4).

## Quantification of Cysts Growth and Morphology

For the quantification of cyst growth, cyst and lumen formation in PEG and Matrigel hydrogels up to 15 z-stacks per condition and time point were acquired. Between 2 and 3 z-stacks were taken from non-overlapping regions within a hydrogel. Z-stacks were performed from the bottom to the top of the hydrogel with fixed height and image sampling number for all condition. Cyst size, number of cysts and number of single- or multi-lumen formed per imaged z-stack was then analyzed using Matlab Image processing toolbox. The cyst formation frequency was determined by relating the cyst number per z-stack to the corresponding initial cell number per z-stack. Lumen formation was also characterized by staining for actin, E-cadherin, and DNA as described above.

## Statistical Analysis

All values are expressed as mean±SEM except for the rheological and viability measurements (mean±SD). Statistical significance was determined using Kruskal-Wallis *t*-test. The dynamics of cyst formation and single lumen production were analyzed in R (3.6.0) or the R project for statistical computing downloadable from https://www.r-project.org/. The fitted curves were generated with the GLM function fed with the models: Cyst frequency ~ scaffold * log(Day) and single lumen ~ scaffold * exp(Day). The obtained *p*-values were corrected for multiple comparisons with the Benjamini and Yekutieli method (BY) (Benjamini and Yekutieli, 2001).

## RESULTS

### Normal Rat Cholangiocytes Form Cysts Upon Culture in 3D Scaffolds

Triggering morphogenetic processes and achieving high levels cell-self organization is best recapitulated when initiated in the presence of a three-dimensional support or scaffold. The initial objective of this study was to identify 3D microenvironments supportive of cholangiocytes self-organization. Based on observations made with other epithelial cell lines we first seeded Normal Rat Cholangiocytes (NRC) in both synthetic, fully defined polyethylene glycol (PEG) hydrogels and in natural extracellular matrix-based Matrigel (MT) (Figure 1A). We quickly confirmed that NRCs can survive and proliferate in both
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Normal Rat Cholangiocyte (NRC) encapsulation and viability assessment. (A) Schematic representation of the used encapsulation procedure showing the differences between PEG (synthetic) and Matrigel (natural) hydrogels in terms of composition. In order to assess the self-organization capabilities of NRCs in synthetic and natural matrices, harvested cells were seeded into PEG hydrogels and Matrigel. Eight-arm PEG macromers bearing a vinyl sulfone function at the end of each arm were crosslinked with different peptides sensitive to matrix metalloprotease degradation in order to produce synthetic hydrogels. These synthetic PEG scaffolds were further functionalized with integrin ligands and proteolytic degradation capabilities in order to ensure NRC survival whereas in Matrigel encapsulated NRCs naturally formed cysts. (B) Brightfield and fluorescent images of live/dead staining of NRCs cultured in PEG and MT for 3 h or 10 days. The NRC cells when maintained for 10 days in culture formed hollow 3D structures, termed cysts. (C) When viability was assessed we found that the obtained structures were viable with PEG hydrogel causing slightly higher cell mortality. Represented PEG hydrogels condition corresponds to 4.5% (w/v) functionalized with 1,200 µM RGD ligand. n = 3. \*Significant at *p* < 0.05 in a Kruskal-Wallis test. Error bars represent SD. Scale bar equal 100 µm.</figcaption></figure>

conditions. Over a culture period of 10 days, cultured NRCs self-organized into hollow spheres or cysts where a very thin layer of epithelial cells marked the boundary between an inner lumen and the extracellular milieu (Video S1). We also observed that unlike in MT, NRCs cannot migrate in PEG and that cyst formation was mostly due to the seeding of small cell aggregates rather than single cells. When assessing cell viability, we found that encapsulated NRC showed globally a good viability in PEG and in MT. At day 10 a modest decrease in viability was observed in PEG, most probably due to the death of the cells that could not form cysts (Figures 1B,C). Overall, the growth kinetics (Video S2) in both materials allowed for the formation of cysts with sizes ranging from 40 to 800 µm. At this stage, we could not observe obvious differences between cysts cultured in either synthetic or natural scaffolds.

## Matrix Stiffness and Integrin Ligand Density Modulate Cyst Formation

After confirming that NRC could self-organize into 3D structures reminiscent of polarized cysts. Unlike unpolarized cell aggregates, the obtained structure clearly showed a large lumen,

cells organized as a monolayer that are staining positive for actin on the apical side and for E-cadherin on the basolateral side. First, we sought to investigate the determinants of cyst formation dynamics across gel composition (natural vs. synthetic matrices), matrix stiffness (or elastic modulus noted *G*') regime and integrin ligand density. Initially, we varied both the PEG concentration and the integrin ligand density (Figure 2A) of the synthetic hydrogel system while documenting the impact of PEG concentration of the elastic modulus (Figure 2B). Unlike ECM-based scaffolds, PEG hydrogels can be produced at varying elastic moduli while maintaining a constant molarity of integrin ligand (Gobaa et al., 2011). In this scenario the impact of increased RGD concentration on elastic modulus remains important and often overlooked (Enemchukwu et al., 2016). In order to quantify the impact of RGD peptides on the mechanical properties of the gel we measured the shear modulus of gels of different concentration while also varying RGD molarity. This showed that the incorporation of 100 µM RGD peptides to the hydrogel network had a very limited impact on the gel stiffness. Increasing the RGD dose up to 1,200 µM led to a 32% decrease of elastic modulus in the stiffest gels (Figure 2B). This RGD-mediated
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Impact of matrix stiffness and integrin ligand density on Normal Rat Cholangiocytes (NRC) cyst formation. (A) Brightfield images of NRC cysts cultured for 10 days. NRCs were encapsulated in PEG hydrogels crosslinked with an MMP sensitive peptide (16R) and varying for concentration (w/v) and integrin ligand density. (B) Quantification of the impact of varying PEG and RGD concentration on the elastic modulus or G' indicating the stiffness of produced PEG hydrogels. There we defined three stiffness domains including soft (2.5% PEG w/v, around 0.5 kPa), intermediate (4.5% PEG w/v, around 3 kPa) and stiff gels (6.5% PEG w/v, >5.5 kPa). The addition of up to 1,200 µM RGD had an overall moderate impact on the elasticity regime. (C,D) Quantification of cyst size and cyst frequency, respectively, at day 10 across soft, intermediate, and stiff polyethylene glycol gels and in Matrigel. (E) Dynamic study showing the evolution of the frequency of cyst production across the studied conditions. *, **, and *** Significant at *p* < 0.05, *p* < 0.01, and *p* < 0.001, respectively. *n* = 4. Error bars represent SEM. Scale bar equal 100 µm.</figcaption></figure>

decrease in stiffness is relatively modest when compared to the impact of varying PEG concentration. Altogether, the measured rheological data showed that varying PEG concentration allowed the production of soft (around 0.5 kPa), intermediate (2–4 kPa) and stiff (around 7 kPa) gels while the measured stiffness of the Matrigel scaffolds was several orders of magnitude lower. Furthermore, we also investigated the impact of gel proteolytic degradation by tuning the matrix metalloproteases sensitivity of PEG hydrogels according to Patterson and Hubbell (2010). There we found that cell survival was independent of the proteolytic degradation of the substrate (Figure S3). We also demonstrated that the most sensitive gels crosslinked with the 16R peptide allowed systematically for a higher efficiency of cyst formation across the entire range of gel stiffness and RGD concentrations (Figure S4).

When we quantified the average cyst size, we found that softer gels favored the formation of larger cysts with the largest cysts being systematically produced in Matrigel (approximatively 20% larger than in the best PEG condition). In addition, this experiment also showed that the presence of RGD, in a dose dependent manner, was critical for the formation of large cysts (Figure 2C). Interestingly, the measurement of cyst sizes across the nine RGD/Stiffness conditions showed that while increasing RGD concentration in soft gels dramatically increased

cyst size, it only had a much lower (but still significant *p* < 0.001) effect in stiff gels. This indicates that stiffness seems to be the parameter that determines the possible maximum cyst size while the addition of RGD allows for the realization of that potential. Conversely, high stiffness clearly overrides the effect of integrin engagement and prohibits the formation of larger NRC cysts. This trend was further confirmed when we looked at the frequency of cyst formation (Figures 2D,E). There we found that increasing the elastic modulus of the hydrogels significantly reduced the frequency of cyst formation. When looking at maximal cyst formation efficiency we found that Matrigel and soft synthetic gels (1,200 µM RGD) had reached comparable levels. We also found that removing RGD from soft PEG drastically limited cyst formation (Figure S5). Finally, when quantifying cyst formation efficiency over 10 days we confirmed that both soft PEG and Matrigel imposed a similar logarithmic increase. Increasing the matrix stiffness was sufficient to significantly (*p* < 0.001) reduce this efficiency while the logarithmic trend was maintained (Figure 2E). Taken together these data clearly show that both, cholangiocytes cyst formation frequency and size are determined to a large extent by variation in the elastic modulus of the substrate. The presence of RGD ligand in high concentration further improved this trend.
Establishment of Cyst Polarity

When cultured in 3D, epithelial cells such as the kidney cell line (MDCK) self-assemble in hollow spheres and recapitulate the morphogenetic program leading to the production of rudimentary epithelial organoids (McAteer et al., 1986). Detailed analysis showed that interactions between key integrins and the extracellular matrix (ECM) regulate the establishment of cell polarity and morphogenesis (Manninen, 2015). In both cases whether NRC cysts can become polarized upon cultured in a synthetic 3D scaffold was investigated first by immunochemistry and rhodamine transport assay. There we found that large NRC cysts stain positive for basolateral E-Cadherin and apical F-actin (Figure 3A; Figure S6). Canonical epithelial polarity was observed in cysts of different dimensions across the tested 3D culture conditions. NRC cysts with inverted polarity could also be observed as described elsewhere (Enemchukwu et al., 2016). However, these observations were very seldom and did appear to correlate with either stiffness or RGD concentration. When looking at the shape of the nuclei of the cholangiocytes lining the cyst wall, we found that in Matrigel these structures seemed to be more elongated and placed tangentially to the formed circle whereas in PEG the nuclei seemed rounder, organized into a mixed mono- and/or double-layer structure (Figure 3A). *In vivo*,

the morphology of the nuclei in the lining of the bile duct seems to vary from round to oval with preferentially a basal localization in the cell. On the functional level we found that, in MT or in PEG hydrogels, formed cysts could accumulate rhodamine 123 in their lumen and that this accumulation was blocked by the addition of verapamil in the culture medium. This clearly indicates that the accumulation of rhodamine is resulting from transmembrane channel protein-multidrug resistance protein-1 (MDR1) activity and not passive diffusion (Figure 3B). Furthermore, we could show that the frequency of rhodamine accumulating cysts was reduced when matrix stiffness was increased (Figure 3C). These data suggest that NRC cysts cultured in either PEG or MT are showing apicobasal polarity and functioning as mature epithelial cysts.

Morphogenesis and Lumen Formation

When imaging the different structures obtained after culturing NRC in PEG and MT hydrogels, we have identified different cell organizations including large cysts including a unique lumen bordered by a thin layer of cells (SL) (Figure 3A; Video S1), spherical cell structures containing multiple gaps (ML) (Figure 4A; Figure S7) and small and compact cell aggregates. We found that the hydrogel recipe was determining to a certain extent the frequency of each structure. Quantifications preformed on different scaffolding conditions showed that decreasing the hydrogel stiffness and increasing the RGD concentration increased the frequency of SL structures (Figures 4B-D). More interestingly the effect of stiffness seemed to plateau for gels between 2.5 and 4.5% (w/v) PEG. We observed that NRCs cultured in 2.5% PEG hydrogel (G' ≈ 0.5 kPa) showed a SL frequency of 73 and 56% for high and low RGD concentrations respectively. The formation of SL was decreased by 20% when using stiffer gels. An even further decrease was observed when using low RGD in the stiffest gels. This indicates that above a certain elasticity threshold (between 0.5 and 2.5 kPa) the luminal organization seems to be rather independent of the hydrogel stiffness whereas RGD density remains an important determinant for the formation of single lumen with normal polarity. Although comparable at day 10, the frequency of SL cyst formation did not appear to follow identical trajectories in PEG vs. MT. Observing the dynamics of SL cyst formation over the course of 10 days revealed that the formation of SL cysts seemed to follow an exponential law in PEG whereas cyst formation dynamics was found more linear in MT. These dynamics reveal that the gap in percent of SL cyst formation is strongly decreasing over time when comparing MT and PEG (Figure 4C). When investigating cyst organization in relation to size we found that ML cysts are significantly smaller than SL cysts. This observation remained true for longer observations (2 weeks plus) (Figure S7). We also found that limiting the dose of RGD clearly restricted the growth of the produced ML cysts (Figure 4D). This experiment also showed that cyst size and organization are not fully correlated and that the two parameters could respond to different microenvironment cues.

<figure><img src="image_5.png" /></figure>

**FIGURE 3** | Phenotypic and functional assessment of Normal Rat Cholangiocytes (NRC) cysts produced in both polyethylene glycol and Matrigel hydrogels. (A) Scanning confocal microphotographs of immunostained NRC cysts after 10 days of culture. Both structures stain positive for actin (apically) and E-cadherin (basolaterally). (B,C) Widefield microphotographs and subsequent quantification showing the accumulation of rhodamine 123 (Rhd. 123) in the lumen of NRC cysts. Adding the channel blocking agent verapamil (VP) reduced this accumulation. \* and \*\*\* Significant at *p* < 0.05 and *p* < 0.001 respectively. *n* = 4. Error bars represent SEM. Scale bar equal 100 µm.
<figure><img src="image_6.png" /></figure>

**FIGURE 4** | Impact of the microenvironment on lumen formation. (A) Immunostaining of f-actin and E-cadherin localization in normal rat cholangiocyte (NRC) cyst with multi-lumen phenotype. The observation of Normal Rat Cholangiocytes (NRC) cysts formation in different matrices showed that some NRC cysts could comprise multiple lumens. (B) Quantification of single lumen formation of NRC cysts during 10 day. The quantifications showed that this phenomenon was prevalent in PEG hydrogel and marginal in Matrigel. We also found that the fraction of single lumen cysts in PEG tended to exponentially increase over time. (C) The analysis of the "multi-lumen" and "single-lumen" phenotype frequency indicated that a soft microenvironment containing 1,200 µM RGD caused the lowest rates of multi-lumen formation. (D) Finally, cyst size analysis showed that variation in RGD concentration seemed to alter the multi-lumen/single-lumen ratio. SL and ML is for single-lumen and multi-lumen, respectively. *, **, and *** Significant at *p* < 0.05, *p* < 0.01, and *p* < 0.001, respectively. *n* = 4. Error bars represent SEM. Scale bar equal 100 µm.

DISCUSSION

The ECM-cholangiocyte interaction is essential for the formation
of the biliary tree. *In vitro*, the ECM is classically mimicked
by natural hydrogels known to promote well-cholangiocyte
differentiation and morphogenesis (Tanimizu et al., 2007).
Understanding the role of the microenvironment cues on bile
duct formation is clearly limited by the type of assays that
can recapitulate this process *in vitro*. On one hand, hydrogels
based on naturel ECM allow to faithfully recapitulate key
morphogenetic processes (Tanimizu et al., 2007). However,
mechanistic studies in these systems are difficult because of the
intertwined biochemical and biophysical signaling. On the other
hand, synthetic hydrogels are more amenable to experiment
design and systematic variation (Ranga et al., 2014) but less
efficient in the support of biological processes. In this work we
clearly benefitted from these developments in order to produce
synthetic, and fully defined scaffolds capable of efficiently
supporting cholangiocyte cyst formation.

Structure wise, epithelial organoids are organized as a
monolayer of polarized epithelial cells surrounding a central
lumen. Two typical structures can emerge from this basic
concept: cysts or tubes. Nevertheless, the precise underlying

morphogenetic mechanism remains elusive and the contribution
of the microenvironment is yet to be quantified (Yu et al.,
2007; Martin-Belmonte et al., 2008). Other studies have shown
that primary biliary mouse epithelial cells (BECs) do form
cysts from a single cell when cultured in Matrigel and that
the resulting cyst formation efficiency is determined by the
initial cell number in the aggregates (Rizki-Safitri et al.,
2018). Similarly, the encapsulation of single MDCK cells was
found sufficient to trigger cyst formation in synthetic and
natural hydrogels (Martin-Belmonte et al., 2008; Enemchukwu
et al., 2016). Interestingly, the cell aggregation route for the
formation of MDCK cysts was also shown pertinent when
these cells were encapsulated at high density. Unlike NRCs
or MDCKs, other cell types showing very limited proliferation
capabilities such as epithelial lung cells (AT II) still produce
cysts with sizes directly linked to initial cell density (Yu
et al., 2007). Our experiments showed that the encapsulation
of aggregates of multiple NRC resulted in cyst formation by
cellular rearrangement whereas single cells systematically failed
to proliferate and organize. Similar findings were made for
NRC cyst formation in MT with the difference that single
cells could migrate and aggregate before proliferating and
forming cysts.
When stiffness was considered, we demonstrated that matrix elasticity variation modulates cyst formation (Figure 2D). Strikingly, cyst formation was found to be very comparable in soft PEG gels and in MT despite the order of magnitude difference in elastic moduli. This is probably due to the fact that both MT and soft PEG conditions are still below the elastic modulus of a healthy liver (~1 kPa) (Yeh et al., 2002). Globally, the cyst formation process was found to tolerate departure from the healthy liver condition. Even recapitulating fibrotic liver conditions (~5 kPa, Stiff PEG) did not fully block the process. Very interestingly this observation is consistent with the fact that an abnormal increase in intrahepatic biliary mass is taking place during the earliest stages of liver fibrogenesis (Schuppan et al., 2018). In comparison, when encapsulating MDCK cells cyst formation is completely blocked at 7 kPa (Enemchukwu et al., 2016). Matrix elasticity was found to play a key role in other systems including intestinal organoids and lung adenocarcinoma spheres (Gill et al., 2012; Raza et al., 2013; Gjorevski et al., 2016; Broguiere et al., 2018). In these systems, relatively high initial stiffness was needed to drive cell compaction. A quick softening over time was also important in order to help with the 3D rearrangement of cells.

When considering the biochemical cues, we first tried to recapitulate physiologically relevant integrin ligand densities by adjusting the RGD peptide concentration in the PEG hydrogels. However, as the *in vivo* density of RGD motifs accessible to intrahepatic cholangiocytes is hard to assess we employed 100 µM as the best approximation based on the fibrin structure (Raeber et al., 2005) and 1,200 µM as the most efficient concentration for cysts derivation in an equivalent MDCK model (Enemchukwu et al., 2016). In these conditions NRC cyst formation frequency was found independent of the level of integrin engagement. This is in good agreement with the work performed on MDCK cells in synthetic hydrogels (Enemchukwu et al., 2016). This way we found that cyst size (growth of cysts) and thus lumen formation relied on integrin engagement (Figures 2C, 4B). The growth in diameter of NRC cysts is known to be fueled by fluid secretion into the cyst lumen and cell proliferation (Neufeld et al., 1992; Doctor et al., 2007). Our results suggest that hydrogel elasticity and cell adhesion can synergistically determine the growth of NRC cysts. When comparing Matrigel and PEG-based systems, Matrigel systematically produced larger cysts. This might be due to the different mechanisms of lumen formation, where NRC cysts cultured in MT achieved cyst polarization at an early time point of cyst formation followed by rapid cyst expansion. The formation of a central lumen is a key feature of epithelial morphogenesis. The mechanism of lumen formation in MDCK cysts grown in natural hydrogels depends on the used ECM composition. MDCK cells cultured in collagen I showed later cell polarization during cyst formation and formed lumen by center cell apoptosis passing through a multi lumen state, whereas MDCK cysts grown in laminin rich Matrigel induced early cell polarization and formed lumen through cell proliferation with almost negligible center cell apoptosis (Martin-Belmonte et al., 2008). This was further confirmed when MDCK cysts were found to form luminal

space even when apoptosis was suppressed, hinting to the involvement of multiple lumen forming mechanisms (O'Brien et al., 2002). In comparison, cysts derived from lung epithelial cells cultured in natural hydrogels formed lumen through cell rearrangement without cell proliferation and cell apoptosis (Yu et al., 2007).

In general, epithelial cysts grown in natural hydrogels clearly show single lumen phenotype predominantly. In synthetic hydrogels different phenotypes were observed for NRC and MDCK cells (Enemchukwu et al., 2016). We demonstrated that NRC cyst grown in PEG hydrogels or in MT both formed lumen with a typical apical-basolateral polarity (Figure 3A). The accumulation of Rhodamine 123 in the produced cysts hints to a proper and functional organization of the cells. Staining for primary cilia and quantifying the response to farnesoid X receptor could help further confirm this aspect. The mechanisms of lumen formation for cysts cultured in PEG hydrogels was found to be based on cell rearrangement, rather than on early cell polarization as shown in MT. Hydrogel stiffness and RGD concentrations were found to be linked to polarization and to the frequency of single lumen formation. A high integrin engagement was shown to enhance the single lumen phenotype over multi-lumen phenotype in NRC cysts for Soft and Stiff PEG hydrogels, respectively. A positive correlation between lumen formation and adhesive ligand density of a synthetic scaffold was also demonstrated in lung adenocarcinoma models (Gill et al., 2012). Not providing polarization cues in synthetic 3D culture systems systematically delays the acquisition of the cyst polarity (Martin-Belmonte et al., 2008; deLeon et al., 2012). For MDCKs a delayed polarization led to the production of multiple lumens that could fuse at later stages. The size of multi-lumen MDCK cysts cultured in natural and synthetic hydrogels was reported to be similar to cysts containing a single lumen (Martin-Belmonte et al., 2008; Enemchukwu et al., 2016). In our hands a majority of NRC cysts with a single lumen derive from multi-lumen structures. In some cases, these multi-lumen cysts persist. In the MDCK system, the multi-lumen cysts phenotype was found associated with the mislocalization of basal components leading to the accumulation of ECM in the luminal space (Torkko et al., 2008). Unlike MDCKs, the NRC multi-lumen cysts were found to be systematically smaller when compared to single-lumen cysts. The suppression of outer layer coat complex II in human intestinal epithelial cells (Caco-2) cultured in 3D resulted in cysts with defective lumen expansion attributed to a faulty ECM signaling when spatially directing cell division. Similar to multi-lumen NRCs, these cysts show a reduced cell number, multi-lumen formation, and smaller sizes (Townley et al., 2012). We speculate that providing NRCs cultured in PEG with a stronger polarization cues by replacing RGD with laminin or laminin derived peptides could significantly stimulate single-lumen formation and expansion of NRC cysts.

In summary, we have established a tunable PEG-based hydrogel system that can be used for cholangiocyte cyst formation and morphogenesis studies. We demonstrated that making these gel cell responsive by adding MMP sensitivity and integrin ligand contributed to meet the cyst
production capacity of NRCs encapsulated in MT. We expect that in the future, this system will be employed to study pathogenicity of biliary diseases and to identify therapeutic targets.

## DATA AVAILABILITY STATEMENT

The datasets generated for this study are available on request to the corresponding author.

## AUTHOR CONTRIBUTIONS

AF carried out the experiments performed in PEG hydrogel including cell culture, hydrogel casting, and cell encapsulation, helped with the main figures, and arranged the supplementary figures. NM Performed the rheology tests on PEG gels. LB and EG carried out the experiments in natural hydrogels, including cell culture, cell encapsulation, microscopy, and rheology tests on Matrigel. SG, PD-W, and AF produced together a first outline of the manuscript. PD-W wrote part of the introduction. SG arranged the four main figures, wrote the abstract, part of the introduction, the results, and the discussion. All the authors reviewed the manuscript.

## REFERENCES

Benam, K. H., Villenave, R., Lucchesi, C., Varone, A., Hubeau, C., Lee, H. H., et al. (2016). Small airway-on-a-chip enables analysis of human lung inflammation and drug responses in vitro. *Nat. Methods* 13, 151–157. doi: 10.1038/nmeth.3697

Benjamini, Y., and Yekutieli, D. (2001). The control of the false discovery rate in multiple testing under dependency. *Ann. Stat.* 29, 1165–1188. doi: 10.1214/ aos/1013699998

Bogert, P. T., and LaRusso, N. F. (2007). Cholangiocyte biology. *Curr. Opin. Gastroenterol.* 23, 299–305. doi: 10.1097/ MOG. 0b013e3280b079fb

Boyer, J. L. (2013). Bile formation and secretion. *Compr. Physiol.* 3, 1035–1078. doi: 10.1002/cphy.c120027

Broguiere, N., Isenmann, L., Hirt, C., Ringel, T., Placzek, S., Cavalli, E., et al. (2018). Growth of epithelial organoids in a defined hydrogel. *Adv. Mater.* 30:e1801621. doi: 10.1002/adma.201801621

Broutier, L., Andersson-Rolf, A., Hindley, C. J., Boj, S. F., Clevers, H., Koo, B. K., et al. (2016). Culture and establishment of self-renewing human and mouse adult liver and pancreas 3D organoids and their genetic manipulation. *Nat. Protoc.* 11, 1724–1743. doi: 10.1038/nprot.2016.097

Cheung, A. C., Lorenzo Pisarello, M. J., and LaRusso, N., F. (2018). Pathobiology of biliary epithelia. *Biochim. Biophys. Acta Mol. Basis Dis.* 1864(4 Pt B), 1220–1231. doi: 10.1016/j.bbadis.2017.06.024

De Assuncao, T. M., Jalan-Sakrikar, N., and Huebert, R. C. (2017). Regenerative medicine and the biliary tree. *Semin. Liver Dis.* 37, 17–27. doi: 10.1055/s-0036-1597818

De Assuncao, T. M., Sun, Y., Jalan-Sakrikar, N., Drinane, M. C., Huang, B. Q., Li, Y., et al. (2015). Development and characterization of human-induced pluripotent stem cell-derived cholangiocytes. *Lab. Invest.* 95, 684–696. doi: 10.1038/labinvest.2015.51

de Groen, P. C., Vroman, B., Laakso, K., and LaRusso, N. F. (1998). Characterization and growth regulation of a rat intrahepatic bile duct epithelial cell line under hormonally defined, serum-free conditions. *In Vitro Cell Dev. Biol. Anim.* 34, 704–710. doi: 10.1007/s11626-998-0066-1

deLeon, O., Puglise, J. M., Liu, F., Smits, J., ter Beest, M. B., and Zegers, M. M. (2012). Pak1 regulates the orientation of apical polarization

## FUNDING

This work received the financial support of both the iLite RHU program (grant ANR ANR-16-RHUS-0005) and the Institut Carnot Pasteur. France-BioImaging infrastructure was supported by the French National Research Agency (ANR 10-INSB-04, Investments for the Future).

## ACKNOWLEDGMENTS

The authors would like thank Sébastien Brûlé from the protein quality control service of the molecular biophysics facility in Institut Pasteur for technical assistance with MALDI experiments. We also would like to extend our thanks to Françoise Baleux for her help with measuring the net peptide. Finally, we acknowledged the support of UTechS PBI (Imagopole) a France-BioImaging infrastructure.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe. 2019.00417/full#supplementary-material

and lumen formation by distinct pathways. *PLoS ONE* 7:e41039. doi: 10.1371/journal.pone.0041039

Dianat, N. H., Dubois-Pot-Schneider, Steichen, C., Desterke, C., Leclerc, P., Raveux, A., et al. (2014). Generation of functional cholangiocyte-like cells from human pluripotent stem cells and HepaRG cells. *Hepatology* 60, 700–714. doi: 10.1002/hep.27165

Doctor, R. B., Johnson, S., Brodsky, K. S., Amura, C. R., Gattone, V., Fitz, J., et al. (2007). Regulated ion transport in mouse liver cyst epithelial cells. *Biochim. Biophys. Acta* 1772, 345–354. doi: 10.1016/j.bbadis.2006.11.006

Enemchukwu, N. O., Cruz-Acuna, R., Bongiorno, T., Johnson, C. T., Garcia, J. R., Sulchek, T., et al. (2016). Synthetic matrices reveal contributions of ECM biophysical and biochemical properties to epithelial morphogenesis. *J. Cell Biol.* 212, 113–124. doi: 10.1083/jcb.201506055

Engler, A. J., Sen, S., Sweeney, H. L., and Discher, D., E. (2006). Matrix elasticity directs stem cell lineage specification. *Cell* 126, 677–689. doi: 10.1016/j.cell.2006.06.044

Gill, B. J., Gibbons, D. L., Roudsari, L. C., Saik, J. E., Rizvi, Z. H., Roybal, J. D., et al. (2012). A synthetic matrix with independently tunable biochemistry and mechanical properties to study epithelial morphogenesis and EMT in a lung adenocarcinoma model. *Cancer Res.* 72, 6013–6023. doi: 10.1158/0008-5472.CAN-12-0895

Gjorevski, N., Sachs, N., Manfrin, A., Giger, S., Bragina, M. E., Ordonez-Moran, P., et al. (2016). Designer matrices for intestinal stem cell and organoid culture. *Nature* 539, 560–564. doi: 10.1038/nature20168

Gobaa, S., Hoehnel, S., Roccio, M., Negro, A., Kobel, S., Lutolf, M., et al. (2011). Artificial niche microarrays for probing single stem cell fate in high throughput. *Nat. Methods* 8, 949–955. doi: 10.1038/nmeth.1732

Hashimoto, W., Sudo, R., Fukasawa, K., Ikeda, M., Mitaka, T., and Tanishita, K. (2008). Ductular network formation by rat biliary epithelial cells in the dynamical culture with collagen gel and dimethylsulfoxide stimulation. *Am. J. Pathol.* 173, 494–506. doi: 10.2353/ajpath.2008.071024

Huh, D., Kim, H. J., Fraser, J. P., Shea, D. E., Khan, M., Bahinski, A., et al. (2013). Microfabrication of human organs-on-chips. *Nat. Protoc.* 8, 2135–2157. doi: 10.1038/nprot.2013.137

Ishida, Y., Smith, S., Wallace, L., Sadamoto, T., Okamoto, M., Auth, M., et al. (2001). Ductular morphogenesis and functional polarization of normal human
biliary epithelial cells in three-dimensional culture. *J. Hepatol.* 35, 2–9. doi: 10.1016/S0168-8278(01)00078-2

Kido, T., Koui, Y., Suzuki, K., Kobayashi, A., Miura, Y., Chern, E. Y., et al. (2015). CPM is a useful cell surface marker to isolate expandable bi-potential liver progenitor cells derived from human iPS cells. *Stem Cell Rep.* 5, 508–515. doi: 10.1016/j.stemcr.2015.08.008

Kourouklis, A. P., Kaylan, K. B., and Underhill, G. H. (2016). Substrate stiffness and matrix composition coordinately control the differentiation of liver progenitor cells. *Biomaterials* 99, 82–94. doi: 10.1016/j.biomaterials.2016.05.016

Loarca, L., De Assuncao, T. M., Jalan-Sakrikar, N., Bronk, S., Krishnan, A., Huang, B., et al. (2017). Development and characterization of cholangioids from normal and diseased human cholangiocytes as an *in vitro* model to study primary sclerosing cholangitis. *Lab. Invest.* 97, 1385–1396. doi: 10.1038/labinvest.2017.63

Lutolf, M. P., and Hubbell, J. A. (2003). Synthesis and physicochemical characterization of end-linked poly(ethylene glycol)-co-peptide hydrogels formed by Michael-type addition. *Biomacromolecules* 4, 713–722. doi: 10.1021/bm025744e

Manninen, A. (2015). Epithelial polarity-generating and integrating signals from the ECM with integrins. *Exp. Cell Res.* 334, 337–349. doi: 10.1016/j.yexcr.2015.01.003

Martin-Belmonte, F., Yu, W., Rodriguez-Fraticelli, A. E., Ewald, A. J., Werb, Z., Alonso, M. A., et al. (2008). Cell-polarity dynamics controls the mechanism of lumen formation in epithelial morphogenesis. *Curr. Biol.* 18, 507–513. doi: 10.1016/j.cub.2008.02.076

McAteer, J. A., Dougherty, G. S., Gardner, K. D. Jr., and Evan, A. P. (1986). Scanning electron microscopy of kidney cells in culture: surface features of polarized epithelia. *Scan. Electron Microsc.* 1986(Pt 3), 1135–1150.

Mosiewicz, K. A., Kolb, L., A. J., van der Vlies, Martino, M. M., Lienemann, P. S., Hubbell, J. A., et al. (2013). *In situ* cell manipulation through enzymatic hydrogel photopatterning. *Nat. Mater.* 12, 1072–1078. doi: 10.1038/nmat3766

Mueller, S., and Sandrin, L. (2010). Liver stiffness: a novel parameter for the diagnosis of liver disease. *Hepat. Med.* 2, 49–67. doi: 10.2147/HMER.S7394

Murray, K. F., Carithers, R. L. Jr., and AASLD (2005). AASLD practice guidelines: evaluation of the patient for liver transplantation. *Hepatology* 41, 1407–1432. doi: 10.1002/hep.20704

Neufeld, T. K., Douglass, D., Grant, M., Ye, M., Silva, F., Nadasdy, T., et al. (1992). *In vitro* formation and expansion of cysts derived from human renal cortex epithelial cells. *Kidney Int.* 41, 1222–1236. doi: 10.1038/ki.1992.184

O'Brien, L. E., Zegers, M. M., and Mostov, K. E. (2002). Opinion: building epithelial architecture: insights from three-dimensional culture models. *Nat. Rev. Mol. Cell. Biol.* 3, 531–537. doi: 10.1038/nrm859

Ogawa, M., Ogawa, S., Bear, C. E., Ahmadi, S., Chin, S., Li, B., et al. (2015). Directed differentiation of cholangiocytes from human pluripotent stem cells. *Nat. Biotechnol.* 33, 853–861. doi: 10.1038/nbt.3294

Patterson, J., and Hubbell, J. A. (2010). Enhanced proteolytic degradation of molecularly engineered PEG hydrogels in response to MMP-1 and MMP-2. *Biomaterials* 31, 7836–7845. doi: 10.1016/j.biomaterials.2010.06.061

Raeber, G. P., Lutolf, M. P., and Hubbell, J. A. (2005). Molecularly engineered PEG hydrogels: a novel model system for proteolytically mediated cell migration. *Biophys. J.* 89, 1374–1388. doi: 10.1529/biophysj.104.050682

Ranga, A., Girgin, M., Meinhardt, A., Eberle, D., Caiazzo, M., Tanaka, E. M., et al. (2016). Neural tube morphogenesis in synthetic 3D microenvironments. *Proc. Natl. Acad. Sci. U.S.A.* 113, E6831–E6839. doi: 10.1073/pnas.1603529113

Ranga, A., Gobaa, S., Okawa, Y., Mosiewicz, K., Negro, A., Lutolf, M., et al. (2014). 3D niche microarrays for systems-level analyses of cell fate. *Nat. Commun.* 5:4324. doi: 10.1038/ncomms5324

Raynaud, P., Carpentier, R., Antoniou, A., and Lemaigre, F. P. (2011). Biliary differentiation and bile duct morphogenesis in development and disease. *Int. J. Biochem. Cell Biol.* 43, 245–256. doi: 10.1016/j.biocel.2009.07.020

Raza, A., Ki, C. S., and Lin, C. C. (2013). The influence of matrix properties on growth and morphogenesis of human pancreatic ductal epithelial cells in 3D. *Biomaterials* 34, 5117–5127. doi: 10.1016/j.biomaterials.2013.03.086

Rizki-Safitri, A., Shinohara, M., Miura, Y., Danoy, M., Tanaka, M., Miyajima, A., et al. (2018). Efficient functional cyst formation of biliary epithelial cells using microwells for potential bile duct organisation *in vitro*. *Sci Rep* 8:11086. doi: 10.1038/s41598-018-29464-w

Roskams, T. A., Theise, N. D., Balabaud, C., Bhagat, G., Bhathal, P. S., Bioulac-Sage, P., et al. (2004). Nomenclature of the finer branches of the biliary tree: canals, ductules, and ductular reactions in human livers. *Hepatology* 39, 1739–1745. doi: 10.1002/hep.20130

Sampaziotis, F., de Brito, M. C., Madrigal, P., Bertero, A., Saeb-Parsy, K., Soares, F. A. C., et al. (2015). Cholangiocytes derived from human induced pluripotent stem cells for disease modeling and drug validation. *Nat. Biotechnol.* 33, 845–852. doi: 10.1038/nbt.3275

Saneyasu, T., Akhtar, R., and Sakai, T. (2016). Molecular cues guiding matrix stiffness in liver fibrosis. *Biomed. Res. Int.* 2016:2646212. doi: 10.1155/2016/2646212

Schuppan, D., Surabattula, R., and Wang, X. Y. (2018). Determinants of fibrosis progression and regression in NASH. *J. Hepatol.* 68, 238–250. doi: 10.1016/j.hep.2017.11.012

Shiojiri, N., and Sugiyama, Y. (2004). Immunolocalization of extracellular matrix components and integrins during mouse liver development. *Hepatology* 40, 346–355. doi: 10.1002/hep.20303

Soroka, C. J., Assis, D. N., and Boyer, J. L. (2019). Patient-derived organoids from human bile: an *in vitro* method to study cholangiopathies. *Methods Mol. Biol.* 1981, 363–372. doi: 10.1007/978-1-4939-9420-5_24

Takayama, K., Mitani, S., Nagamoto, Y., Sakurai, F., Tachibana, M., Taniguchi, Y., et al. (2016). Laminin 411 and 511 promote the cholangiocyte differentiation of human induced pluripotent stem cells. *Biochem. Biophys. Res. Commun.* 474, 91–96. doi: 10.1016/j.bbrc.2016.04.075

Takebe, T., Zhang, R. R., Koike, H., Kimura, M., Yoshizawa, E., Enomura, M., et al. (2014). Generation of a vascularized and functional human liver from an iPSC-derived organ bud transplant. *Nat. Protoc.* 9, 396–409. doi: 10.1038/nprot.2014.020

Tan, J. L., Tien, J., Pирone, D. M., Gray, D. S., Bhadriraju, K., Chen, C., et al. (2003). Cells lying on a bed of microneedles: an approach to isolate mechanical force. *Proc. Natl. Acad. Sci. U.S.A.* 100, 1484–1489. doi: 10.1073/pnas.0235407100

Tanimizu, N., Miyajima, A., and Mostov, K. E. (2007). Liver progenitor cells develop cholangiocyte-type epithelial polarity in three-dimensional culture. *Mol. Biol. Cell* 18, 1472–1479. doi: 10.1091/mbc.e06-09-0848

Torkko, J. M., Manninen, A., Schuck, S., and Simons, K. (2008). Depletion of apical transport proteins perturbs epithelial cyst formation and ciliogenesis. *J. Cell. Sci.* 121(Pt 8), 1193–1203. doi: 10.1242/jcs.015495

Townley, A. K., Schmidt, K., Hodgson, L., and Stephens, D., J. (2012). Epithelial organization and cyst lumen expansion require efficient Sec13-Sec31-driven secretion. *J. Cell. Sci.* 125(Pt 3), 673–684. doi: 10.1242/jcs.091355

Unadkat, H. V., Hulsman, M., Cornelissen, K., Papenburg, B. J., Truckenmuller, R. K., Carpenter, A. E., et al. (2011). An algorithm-based topographical biomaterials library to instruct cell fate. *Proc. Natl. Acad. Sci. U.S.A.* 108, 16565–16570. doi: 10.1073/pnas.1109861108

Vroman, B., and LaRusso, N. F. (1996). Development and characterization of polarized primary cultures of rat intrahepatic bile duct epithelial cells. *Lab. Invest.* 74, 303–313.

Wells, R. G. (2008). The role of matrix stiffness in regulating cell behavior. *Hepatology* 47, 1394–1400. doi: 10.1002/hep.22193

Yeh, W. C., Li, P. C., Jeng, Y. M., Hsu, H. C., Kuo, P. L., Li, M. L., et al. (2002). Elastic modulus measurements of human liver and correlation with pathology. *Ultrasound Med. Biol.* 28, 467–474. doi: 10.1016/S0301-5629(02)00489-1

Yu, W., Fang, X., Ewald, A., Wong, K., Hunt, C. A., Werb, Z., et al. (2007). Formation of cysts by alveolar type II cells in three-dimensional culture reveals a novel mechanism for epithelial morphogenesis. *Mol. Biol. Cell* 18, 1693–1700. doi: 10.1091/mbc.e06-11-1052

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Funfak, Bouzhir, Gontran, Minier, Dupuis-Williams and Gobaa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.