# A neural crest cell isotropic-to-nematic phase transition in the developing mammalian gut

Nicolas Chevalier, Yanis Ammouche, Anthony Gomis, Lucas Langlois, Thomas Guilbert, Pierre Bourdoncle, Sylvie Dufour

► To cite this version:

Nicolas Chevalier, Yanis Ammouche, Anthony Gomis, Lucas Langlois, Thomas Guilbert, et al.. A neural crest cell isotropic-to-nematic phase transition in the developing mammalian gut. Communications Biology, 2021, 4 (1), pp.770. 10.1038/s42003-021-02333-5 . inserm-03271227

HAL Id: inserm-03271227

https://inserm.hal.science/inserm-03271227

Submitted on 25 Jun 2021

**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>

# A neural crest cell isotropic-to-nematic phase transition in the developing mammalian gut

Nicolas R. Chevalier ID 1✉, Yanis Ammouche<fnref n="1" />, Anthony Gomis<fnref n="1" />, Lucas Langlois<fnref n="1" />, Thomas Guilbert², Pierre Bourdoncle² & Sylvie Dufour ID 3

While the colonization of the embryonic gut by neural crest cells has been the subject of intense scrutiny over the past decades, we are only starting to grasp the morphogenetic transformations of the enteric nervous system happening in the fetal stage. Here, we show that enteric neural crest cell transit during fetal development from an isotropic cell network to a square grid comprised of circumferentially-oriented cell bodies and longitudinally-extending interganglionic fibers. We present ex-vivo dynamic time-lapse imaging of this isotropic-to-nematic phase transition and show that it occurs concomitantly with circular smooth muscle differentiation in all regions of the gastrointestinal tract. Using conditional mutant embryos with enteric neural crest cells depleted of β1-integrins, we show that cell-extracellular matrix anchorage is necessary for ganglia to properly reorient. We demonstrate by whole mount second harmonic generation imaging that fibrous, circularly-spun collagen I fibers are in direct contact with neural crest cells during the orientation transition, providing an ideal orientation template. We conclude that smooth-muscle associated extracellular matrix drives a critical reorientation transition of the enteric nervous system in the mammalian fetus.

<footnote n="1">¹Laboratoire Matière et Systèmes Complexes, Université de Paris/CNRS UMR 7057, Paris, France. ²Institut Cochin, INSERM U1016, CNRS UMR 8104, Université de Paris (UMR-S1016), Paris, France. ³Univ Paris Est Creteil, INSERM, IMRB, Creteil, France. ✉email: nicolas.chevalier@u-paris.fr</footnote>
In vertebrates, most of the ENS derives from vagal enteric neural crest cells (ENCCs). ENCCs invade the gut mesenchyme rostro-caudally early in embryogenesis¹,². Factors affecting the early migration of ENCCs have been the object of considerable scientific scrutiny in the past decades³. Current research efforts have shifted towards post-colonization morphogenetic and physiological developments of the ENS. Recent studies have for instance revealed when different subtypes of ENS neurons differentiate⁴, how the ENS connects to target cells in the gut⁵ and how the submucosal plexus is formed and connects to the extrinsic innervation⁶. Other investigations have unraveled how the ENS starts controlling the motility of the fetal intestine⁷,⁸. The morphogenesis of ENS ganglia has been the object of particular interest. Neuronal cell bodies occupy the central part of ganglia while glial cells are present both inside the ganglion, at its outer periphery, and in the interganglionic fiber tracts⁹. This arrangement as well as the characteristic size of ganglia in embryos emerge as a result of differential cell-cell and cell-extracellular matrix (ECM) adhesion⁹⁻¹². The issue we address in this report is the genesis of the shape (morpho) of the mature murine ENS. The mammalian myenteric plexus (MP) forms a rectangular grid, where ganglia are elongated along the circumference¹³⁻¹⁵, forming “ribs”¹². Here, we reveal for the first time the dynamics of this orientation phase transition and, by using ENCC conditional β1-integrin mutant embryos and *in-toto* second harmonic generation microscopy, reveal that smooth-muscle associated ECM drives this major structural rearrangement of the ENS in the mammalian fetal gut.

## Results

**The fetal mouse enteric nervous system undergoes an orientation phase transition correlated in time with circular smooth muscle (CSM) differentiation.** We used guts obtained from mouse embryos carrying YFP or tdTomato reporter expression in the cytoplasm of all ENCCs (referred as YFP+ and Tomato+, hereafter) and their derivatives (neurons, glial cells)¹⁶,¹⁷. We examined the morphology of the ENS by confocal microscopy in these mice from stages E14.5 through E19.5 in the midgut (Fig. 1a, b), hindgut, and cecum (Fig. S1a, b). We quantitatively assessed cell orientation with the ImageJ “OrientationJ” plugin¹⁸, which yields the average angle of YFP+ structures (0: longitudinal, 90°: circumferential), and the orientation intensity (hereafter called coherence) at the scale of the ganglion size σ (see Materials and Methods). A high coherence indicates a high degree of anisotropy, i.e., that the cell network exhibits a distinct preferential direction. In contrast, a low coherence indicates isotropy, i.e., that the structures are only weakly or not oriented along any direction. Figure 1a shows an example of orientation fields (yellow bars) in E15.5 midgut. The resulting orientation vectors are those of the dominant structure at the scale of the ganglion size, smoothing out smaller structures.

At E14.5, YFP+ cells were distributed homogeneously and isotropically in all parts of the gut (MG *n* = 5, HG *n* = 4, CC *n* = 4). In the following days, YFP+ cells became progressively oriented along the circumference of the GI tract (Fig. 1b, Fig. S1a, b), forming stripes or “ribs”. ENCCs transitioned from a state of low coherence, with a slight longitudinal orientation, to a high coherence, circumferential orientation (Fig. 1c). This transition started at E14.5 in the midgut, and was almost complete by E16.5 (*n* = 4); in the hindgut it started at E15.5 (*n* = 4) and was complete by E17.5 (*n* = 4); in the cecum we found first signs of orientation at E16.5-E17.5 (*n* = 4), and complete circumferential orientation at E19.5 (*n* = 4). We compared the orientation of YFP+ (green, ENCC, Fig. 1a left) and Tuj1+ immunolabeled cells (red, ENCC-derived neurons, Fig. 1a middle) at stage E15.5. Both

labelings overlapped (Fig. 1a right), but YFP delineated cellular bodies while Tuj1 signal revealed the interganglionic fibers. While the YFP+ structure was circumferentially oriented after the nematic transition, interganglionic fibers were predominantly longitudinally oriented at all stages examined (E14.5-E16.5). The resulting structure is a roughly square grid made up of circumferentially oriented ganglia interconnected by longitudinal interganglionic fibers. This was confirmed by high-resolution confocal microscopy (using a Zeiss CSU-W1 microscope) of Tuj1+ immunolabeled whole mounts (Fig. 1d). The higher resolution of this microscope also allowed us to reveal the cellular structure of each ganglion: individual cell bodies are not circumferentially oriented, but the ganglion as the aggregate of these cell bodies is (Fig. 1d). We further monitored the dynamics of the orientation transition by time-lapse imaging E14.5 YFP+ guts in culture for up to 3 days (Fig. 1d, Supplementary Video 1). We found that clusters of YFP+ cells gradually thicken, coalesce and elongate along the circumferential direction while longitudinal connections between adjacent cell clusters are broken down (Supplementary Video 1), resulting in the emergence of rib-like structures. ENCCs transited within 30 h (Fig. 1d, bottom) from an isotropic state to an oriented cell assembly. While orientation in vivo is a robust, systematic phenomenon, we could only observe it in 10-20% of the midgut regions examined (*n* = 10 samples and *n* = 40-50 images along the midgut for each sample) after 2 days in culture (E14.5 + 2). Oriented ganglia could be found in the duodenum, jejunum, or ileum; we did not find a correlation between orientation and the rostro-caudal position along the midgut. Longer culture times resulted in the formation of large aggregates and regions devoid of ENCCs, which indicated a deviation from physiological development.

Because CSM is also circumferentially oriented, we investigated whether the orientation of ENCCs and their derivatives could be driven by CSM. We determined the time of CSM differentiation (CSMD) in the different parts of the GI tract by whole-mount immunohistochemistry for smooth muscle α-actin (Fig. 2a) and by measuring the time of first appearance of contractile waves (motility, Fig. 2b); we have previously demonstrated in the chicken that the emergence is synchronous with CSMD¹⁹. We found that circumferentially oriented α-actin fibers appear at E13.5 in the jejunum, at E14.5 in the hindgut and at E15.5-E16.5 in the cecum (Fig. 2a,c); the onset of contractions was synchronous with the differentiation of CSMD in all regions of the GI tract (Fig. 2b), further confirming this chronology. The time lag between CSM differentiation in the mid- and hindgut and the cecum is consistent with observations in the chicken embryo¹⁹. In every part of the murine gut, orientation of YFP+ cells thus starts 1 day and finishes 3 days after CSMD. This constant time lag suggests that CSM or CSM-associated ECM acts as a cue to guide ENCC circumferential orientation.

**ENCC adhesion through β1-integrins is critical for the orientation transition.** ENCCs express various integrins that control their adhesion to gut ECM²⁰. We previously showed that NCC-specific conditional invalidation of the *Itgb1* gene encoding the β1-integrin subunit (*Itgb1*-cKO) produces aganglionosis of the distal colon²¹ due to defective ENCC-ECM adhesive interactions, altered migration, and exacerbated intercellular adhesion¹¹,¹⁷. At E16.5, the ENS present in the colonized part of the conditional *Itgb1* mutant guts displays a distorted pattern (Fig. 3b) compared to the typical stripe pattern of the ENS in control embryos (Fig. 3a). We quantitatively assessed the orientation of ENCC network in control (*n* = 4) and *Itgb1*-cKO (*n* = 3) E16.5 midguts (Fig. 3c). Control ENCC networks exhibited a circumferential orientation (80.9 ± 2.6°, all results expressed as mean ± SD) with a
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Morphogenesis of the murine fetal ENS. **a** OrientationJ orientation field (yellow vectors) and coherence (α vector length) analysis results for YFP+ cells (left), βIII-tubulin immunolabeling (neurons, interganglionic fibers, middle), and overlap (right). Dashed rectangle: ROI for averaging of orientation & coherence. **b** Orientation of YFP+ cells in the midgut from E14.5 to E16.5. **c** Evolution of the orientation and coherence in the midgut, hindgut, and cecum. Each dot is a different sample (embryo). The dashed blue sigmoid is a guide to the eye, the red arrow indicates the time of circular smooth muscle differentiation (CSMD), see Fig. 2. **d** High-resolution confocal microscopy reveals cell bodies within individual ganglia; βIII-tubulin immunolabeling, max z-projection of confocal stack (**e**) Still shots from Supplementary Video 1 showing progressive circumferential orientation of YFP+ cells in cultured E14.5 jejunum. Bottom graph: orientation and coherence evolution extracted from Supplementary Video 1.</figcaption></figure>

high coherence (0.48 ± 0.05). In contrast, the *Itgb1*-cKO network revealed more longitudinal and scattered orientations (61.7 ± 10.6°), with a lower coherence (0.30 ± 0.08). These differences were even more conspicuous at stage P1 and P14 (see Fig. 2 in Breau et al.²¹). This indicates that β1-integrins mediated ENCC adhesions are required for the proper orientation transition of the ENS.

ENCC are in contact with circular collagen fibers during the orientation transition. We next investigated the nature of the ECM associated with enteric neurons and glia. We used label-free second harmonic generation (SHG) microscopy to reveal the 3D spatial orientation of fibrous collagen relative to that of the Tomato+ ENCC network; this information cannot be gathered from antibody labeling on 2D slices. In the context of the
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Circular smooth muscle differentiation (CSMD) takes place at different times in different parts of the lower GI tract. **a** Whole-mount immunohistochemistry for smooth muscle α-actin. The white arrows indicate the first appearance of circumferential muscle fibers. The scale is indicated at the right and is the same for all images. **b** Average frequency of peristaltic contractions in the mouse midgut, hindgut and cecum at E12.5 (n = 2), E13.5 (n = 20), E14.5 (n = 12), and E15.5 (n = 3), in cycle per minute (cpm). Each data point corresponds to a different embryo (gut). The red arrows indicate the first appearance of contractile waves (motility). Note that in the cecum at E15.5 the contractions were circumferential (i.e., local diameter reduction) only at the base of the cecum, but appeared as "bulk" contractions in the remainder of the cecum. They only became circumferential in the whole cecum as from E16.5. **c** Thin section immunostaining in E16.5 cecum showing full differentiation of the smooth muscle layer at this stage.</figcaption></figure>

intestine, the SHG signal is proportional to the density of collagen I fibers (Materials and Methods).

Collagen I is a permissive α1β1 integrin-mediated adhesion substrate for NCCs²²,²³. We found that circularly oriented collagen fibers were present in the plane and immediately below the myenteric plexus (MP) at the start (E14.5, Fig. 4a, Supplementary Video 2) and after ENCC orientation was complete (E17.5, Fig. 4b, Supplementary Videos 3, 4). This fiber layer was very thin (a few µm) at E14.5 (Fig. 4a), co-located to a large extent with the ENCCs (Fig. 4a, z = 80 µm, merge), but was situated more towards the inside of the intestine, at the boundary between the CSM and the myenteric plexus. SHG microscopy also revealed the collagen-rich serosa (Fig. 4a, z = 80 µm). At E17.5, the collagen fiber layer became much thicker (~50 µm, Fig. 4b). The intensity of the SHG signal just below the myenteric plexus increased 5-fold between E14.5 and E17.5, from 25 ± 5 pixel units (n = 7 stacks from n = 4 guts) to 128 ± 34 (n = 9 stacks from n = 4 guts). At E17.5, collagen fibers not directly in contact with the ENS plexus, in deeper layers of the gut (the submucosa, where most of the collagen is located²⁴), had various orientations; in the duodenum and midgut we found circular, oblique (visible on the left side of Fig. 4b, z = 10 µm), or bidirectional (with fibers at ±45°) orientations. The collagen

distribution in the MP, SMC, and submucosa we report is consistent with immunostaining of collagen I on E18.5 duodenum thin-sections presented by other investigators²⁴. We found similar results in the E17.5 hindgut (Fig. S2). Interestingly, the outer part (closer to the serosa, Fig. S2, z = 0 µm) of the MP in E17.5 hindgut was remarkably less circularly oriented than its inner part (close to the CSM, Fig. S2, z = 10 µm). This indicates that reorientation of the ENS cells is a very local phenomenon, which only occurred in direct contact with the circular ECM fibers. Intense longitudinally oriented collagen fibers were located in the submucosa of the E17.5 hindgut (n = 4/4, Fig. S2, z = 30 µm). ENCCs in the hindgut were not in direct contact with these deeper longitudinal fibers (see gap indicated by arrowhead in Fig. S2, z = 120 µm). We finally compared the fiber pattern revealed by SHG microscopy with the one obtained by collagen I whole-mount antibody immunohistochemistry on E17.5 gut (Fig. S3, Supplementary Video 5). We found that collagen fibers are wound circularly at the level of the myenteric plexus, as found by SHG in Fig. 4b. We could not however reveal collagen deeper in the mucosa, either because the antibody did not penetrate the deeper tissue layers, or because we used a spinning disk microscope that has lower penetration depth than the two-photon SHG method. Antibody staining further labeled
<figure><img src="image_5.png" /><figcaption>**Fig. 3 Conditional β1 integrin knock-out in ENCCs perturbs their circumferential reorientation.** Whole-mount X-Gal staining of (a) control and (b) Itgb1-cKO ENS in E16.5 gut. (c) Orientation and coherence of control and Itgb1-cKO. Each symbol corresponds to a different region analyzed along the midgut (2 to 5 analyzed per embryo); different symbol shapes correspond to different embryos (n = 4 control, n = 3 Itgb1-cKO). \*: statistically significant difference, *p* < 0.05, two-tailed Mann-Whitney test.</figcaption></figure>

individual collagen-rich cells below the myenteric plexus, that we did not observe by SHG microscopy.

## Discussion

We have demonstrated that CSM differentiation and ENS circumferential orientation were temporally tightly correlated, the first preceding the latter by exactly 1 day in all regions of the lower gastro-intestinal tract (Figs. 1, 2). We captured the dynamics of this orientation transition (Fig. 1e, Supplementary Video 1) by ex vivo time-lapse imaging. In analogy to the physics of liquid crystals, the orientation transition can be assimilated to an isotropic-to-nematic transition, where the isotropic phase does not exhibit a preferential direction, whereas the nematic phase is ordered along one direction (the circumference of the gut). We showed that *Itgb1*-cKO ENCCs depleted of β1-integrins display a defective circumferential re-orientation (Fig. 3). While the orientation transition in Itgb1-cKO ENCC is perturbed, it is not completely abolished, suggesting that other types of integrin receptors besides β1 may be involved, or a correct balance between integrins and cadherins¹¹. We revealed circularly oriented collagen I fibers at the MP-CSM interface by second harmonic generation microscopy (SHG, Fig. 4); the density of

these fibers increased during the physiological orientation of the plexus between E14.5 and E17.5. Collagen fibers could provide an ideal template for ENCC orientation, although other ECM molecules might be involved. We conclude that the reorientation of ganglia along the circumference of the intestine is driven by smooth-muscle associated ECM.

Candidate ECM molecules that could serve as a template for ENCC orientation via β1-integrins should fulfill the following criteria: (1) be present in the embryonic gut at the level of the myenteric plexus at the stages investigated (E14.5-E19.5), (2) be permissive to ENCC adhesion, (3) present a circumferential organization, that is of fibrous nature, (4) be associated/produced by the smooth muscle, because we showed that there is a kinetic correlation between the appearance of smooth muscle and the reorientation of ENCCs. Collagen I, III, and V are the main collagen molecules present in the gut, representing respectively 68, 20, and 12% of collagens in the human intestine²⁵. Collagen III and V are generally found associated with collagen I fibers²⁶,²⁷ and for the sake of this discussion we will therefore consider these three fibrous collagen types as one entity. Collagen I have been shown to be permissive to NCC migration²²,²³; in particular, the orientation and translocation of NCC along (parallel to) collagen I fibers has been reported by Davis²⁸. Importantly, collagen I, III, and V have been shown to be produced by fetal human intestinal smooth muscle cells²⁹. We have shown here by second-harmonic generation microscopy that collagen I is present at the right location (myenteric plexus), stage (E14.5-E17.5) and with the proper circumferential orientation. In contrast, collagen IV forms a sheet-like structure which, although permissive to migration²³, is not present at the level of the myenteric plexus in the hindgut when the orientation transition occurs²⁴. Type VI collagen is mostly localized at the basal membrane of the epithelium³⁰, but is also secreted by ENCCs. It however was found to have an inhibitory effect on migration, with overexpression triggering a Hirschsprung type phenotype³¹. Collagen XVIII is only secreted by ENCCs at the colonization wavefront³². We can therefore conclude that, among collagens, collagen I fibers (and associated type III and V) fulfills the criteria of an ideal template ECM for ENCC orientation. We cannot however exclude that other, non-collagenous ECM molecules that support ENCC migration, like fibronectin³¹, could be produced in an oriented way by the CSM and play a role in driving the orientation transition of ENCCs.

Does the orientation transition occur in other species? Whole-mount staining of peripherin and Phox2b in transparized human embryonic midgut at 9 and 13 weeks of development³² (see also http://transparent-human-embryo.com) shows that the stripe pattern of the ENS is present in humans as well: peripherin filaments run parallel to the longitudinal axis whereas Phox2b positive ganglia are oriented circumferentially. The ENS orientation transition in humans probably occurs between weeks 7 and 9 of development, just after CSM differentiation³³,³⁴. In contrast to the mouse colon, the MP of the human colon has a hexagonal honeycomb geometry³⁵⁻³⁷ which resembles more closely that of the chicken hindgut.

When CSMD is inhibited in the chicken midgut, the regular plexus geometry is lost and ENCCs chaotically invade regions closer to the epithelium³⁸. Enteric ganglia in the chicken do not however reorient like in mammals, in spite of the presence of circular collagen I fibers associated with the CSM²². The ENS of the chicken takes on a hexagonal honeycomb geometry shortly after ENCC colonization, and this geometry does not change later in development. Differences in the balance of forces pushing for orientation (ECM-ganglion interactions) and the ones opposing it (cell-cell adhesion and interganglionic fiber tension) are the likely reason for these species-specific variations in the final ENS morphology. Interganglionic fibers in chicken are noticeably
<figure><img src="image_6.png" /><figcaption>**Fig. 4** ENCC and collagen imaging by second-harmonic generation microscopy. **a** E14.5 duodenum, at the start of ENCC reorientation, plane selected from z-stack (Supplementary Video 2) at the level of the myenteric plexus (MP). Collagen fibers are circumferentially oriented. Z = 80 µm is a longitudinal optical section of one border of the gut. The collagen is concentrated at the serosa and at the MP. **b** E17.5 duodenum, after ENCC reorientation. Circumferentially oriented collagen fibers are situated at (z = 0 µm) and just below the MP (z = 10 µm). The strong collagen signal extends further in the submucosa (z = 80 µm, longitudinal optical section of one border of the gut), where it had various orientations (circular, longitudinal or oblique) depending on the sample and sample region examined.</figcaption></figure>

thicker than in mice and longitudinally oriented; the pulling force they exert on ganglia may offset the tendency of ganglia to orient circumferentially with the underlying ECM. An interesting question is the implication of these differences of ENS geometry on the physiology of the intestine. Neurons in the ENS are known to be mechanosensitive³⁹ and we expect that the mechanical stresses (due to the presence of food bolus or spontaneous muscle contractions) experienced by circumferential vs hexagonal ENS lattices will be significantly different.

We shall conclude by advancing the hypothesis that internal mechanical forces are directing the orientation transition via the ECM. It is known that the more important proliferation of epithelial compared to mesenchymal cells induces a circumferential internal mechanical stress⁴⁰ that is present at least as from E12.5 in the murine intestine. This internal stress is responsible for the circular orientation of the CSM in the chicken⁴⁰. Since static

stresses are well-known to orient collagen fibers⁴¹, it is not hard to imagine that circumferential internal stress could orient the CSM-associated ECM, which in turn orients the myenteric plexus by the mechanism outlined here. This points to the potentially long-ranging implications of small internal stresses in organ morphogenesis.

As in central nervous system development⁴², our research stresses the key role played by interactions between enteric neural crest cells and the extra-cellular matrix in driving major structural rearrangements of the enteric nervous system during mammalian fetal development.

**Methods**

**Ethics.** All procedures presented in this manuscript have been approved by the ethics committee ComEth Anses/ENVA/UPEC, APAFiS request number
2019061218208141 (#210725). All experiments were performed in accordance with the ethics guidelines of the INSERM and CNRS.

**Specimens, samples.** Mouse models used in this study are the following. The Cre reporter mice Gt(ROSA)26Sortm1(EYFP)Cos⁴³, and B6.Cg-Gt(ROSA)26Sortm9 (CAG-tdTomato)Hze/J⁴⁴ were refereed as YFPfl/fl and Tomatofl/fl, respectively. A transgenic mouse line in which the transgene is under the control of the 3-kb fragment of the human tissue plasminogen activator (Ht-PA) promoter Tg(PLAT-cre)116Sdu¹⁶, was referred as Ht-PA::Cre. YFPfl/fl or Tomatofl/fl were cross with Ht-PA::Cre to generate embryos carrying the reporter protein in migrating NCCs and their derivatives.

The mice carrying floxed *Itgb1* allele Itgb1tm1⁴⁵ and the mice carrying one *Itgb1* null allele Itgb1tm2⁴⁶ were referred as beta1fl/fl and beta1neo/+, respectively. Ht-PA::Cre were crossed with beta1neo/+ to generate Ht-PA::Cre; beta1neo/+, which were then crossed with beta1fl/fl to generate embryos with Ht-PA::Cre; beta1neo/fl or Ht-PA::Cre; beta1+/fl genotype as previously described¹⁷. Ht-PA::Cre; beta1fl/+ embryos that are heterozygous for *Itgb1* in migrating NCCs were referred to as controls because *Itgb1* heterozygosity has no effect on mice phenotype⁴⁷. Ht-PA::Cre; beta1neo/fl embryos that carry the conditional mutation for *Itgb1* in migrating NCCs were referred to as *Itgb1*-cKO. In both controls and *Itgb1*-cKOs, Cre recombinase-mediated deletion of the beta1fl allele leads to the expression of beta-galactosidase under the control of the *Itgb1* promoter⁴⁶. This allowed us to visualize the NCCs and their derivatives including ENS using Xgal staining. The day of the vaginal plug was considered E0.5.

**Statistics and reproducibility.** All sample numbers indicated in this report correspond to different embryos (guts, biological replicates). Information was collected from 1 or 2 different litters per age. All experiments involving culture were technically replicated at least twice, i.e. they were performed on two different days with fresh samples following the same procedure each time. A minimum of *n* = 4 samples constitutes each group presented in this report. Statistical analysis was performed with the two-tailed Mann-Whitney test and was considered significant at *p* < 0.05.

**Organotypic culture.** E14.5 mice guts were cultured in 1 mL DMEM with 1% penicillin-streptomycin in 35 mm diameter Greiner culture dishes whose bottom was covered with a thin layer of Sylgard. Guts were pinned with 100 µm steel pins (Euronexia) at their extremities (stomach and distal hindgut) and cultured for 72 h in a humidified incubator at 37 °C, in a 5% CO₂ atmosphere. For live time-lapse confocal imaging, the gut was embedded at the surface of a 2% type-VII agarose gel to reduce sample displacement. Live confocal imaging was performed on a DSD2 confocal microscope at x10 magnification, at a frequency of 1 image/10 min, in a humidified 5% CO₂ atmosphere. Z-stacks were collected and the maximum projected. Jitter motion of the gut (Supplementary Video 1) is due to spontaneous CSM contractions that propagate along the sample length. These contractions occur at a smaller time scale (on the order of 1–10 s) than the gradual orientation of the ENS. Fixed-time imaging of guts was performed on an Olympus spinning disk confocal microscope with GFP excitation/emission filters (488/512 nm). Given the small proportion of mouse gut samples (10–20%) that exhibited re-orientation in culture and the limited availability of our confocal microscope for prolonged time-lapse imaging, we were only able to produce one video (Supplementary Video 1) in which we could track the re-orientation of YFP+ cells along the circumference.

**Monitoring of gut motility.** Spontaneous motility of E12.5-E14.5 mice gut was monitored by placing them in 1 mL of DMEM in individual dishes at 37 °C and time-lapse imaging with a binocular (MZ series, Leica) and camera (Stingray) at a 1 Hz frequency for at least 5 min. Space-time diagrams of motility¹⁹ in different regions were extracted with the imageJ “Reslice” function.

**Immunohistochemistry.** For IHC of frozen sections, guts were fixed for 1 h in a 4% PFA in PBS solution, washed in PBS, then let overnight in 30% sucrose in water solutions, and embedded the next day in OCT compound (VWR) on dry ice. Thin (14 µm) slices were cut at -20 °C in a Leica cryotome and deposited on Ther-mofrost glass slides. After rehydration, the slides were blocked for 15 min in a 1% BSA and 0.1% triton in PBS solution, the slides were then incubated overnight in anti α-smooth muscle actin antibody (Abcam, ref5694, dilution 1:2000) and/or anti βIII-tubulin antibody (Tuj-1, Abcam, ref14545, dilution 1:1000) solution composed of 1% BSA in PBS; the following day, after washing, CY3- and A488-conjugated secondary antibodies (ThermoFisher, dilution 1:400 in PBS) were applied for 2 h. The slides were washed in PBS, and immediately imaged in a layer of PBS with a confocal microscope (low-resolution Olympus system or high-resolution Zeiss CSU-W1). For whole-mount Tuj1 staining, after PFA-fixing, the gut was blocked and permeated in 1% BSA and 0.1% triton in PBS overnight, immersed in anti βIII-tubulin antibody, as described above, for one day, washed and immersed in Alexa488 secondary antibody for at least 4 h. For whole-mount collagen I (Sigma SAB1402151) IHC, primary antibody was incubated at 1:100 for 3 days, and

revealed with Alexa647 secondary antibody (1:400, overnight). The detection of beta-galactosidase activity was revealed by whole-mount X-gal staining as described in Breau et al.²¹.

**Harmonic and fluorescence imaging.** Demesenterized guts were embedded at the surface of a 3% agarose gel immersed in PBS. Images were acquired with an upright Leica SP8 DIVE (Leica microsystems Gmbh, Wetzlar, Germany) coupled with a Coherent Discovery (Coherent Inc., Santa Clara, CA, USA) pulsed laser and controlled by LAX imaging software. Both 552 nm continuous laser and 1040 nm pulsed laser were used to acquire tdTomato fluorescence and second harmonic generation, respectively, through a ×25 objective in PBS immersion. Because of the intrinsic high-penetration depth (~100 µm) of two-photon microscopy techniques, SHG can be used to image the orientation of collagen fibers in the relevant optical sections, without the need to actually perform physical sectioning of the sample. In contrast, antibody staining of collagen is most often performed on tissue sections, on which it is very difficult to assess the overall orientation of cellular and extracellular matrix components (e.g. the orientation of ENCCs wound around a cylinder as in Fig. 1 cannot be detected on physical sections). SHG allows for quantitative evaluation of collagen fiber density. SHG is sensitive to collagen type I and type II⁴⁸, but since only type I is found in the gut (where it is the main collagen type expressed²⁵), we image here the collagen I distribution. Comparative images of collagen I distribution by immunostaining and of SHG generation in liver samples have for example shown almost complete overlap⁴⁹. We note that SHG is insensitive to non-fibrillar collagen-like type IV; type III does not generate a measurable SHG signal⁴⁸.

**Image analysis.** Orientation of the ENS in the different segments of the mouse was measured with the ImageJ “Orientation” developed by the Biomedical Imaging Group at EPFL¹⁸. This plugin computes a vector describing the mean orientation of structures of size σ in a region-of-interest (ROI). We used the Gaussian gradient method and a size σ close to that of the ganglia; choosing too small σ traces the contour of individual ganglia and does not give meaningful information about their average orientation. The direction of the vector represents the orientation of the subimage while its length (proportional to the coherence) quantifies the degree of anisotropy. The sides of the gut present an artefactual longitudinal orientation induced by 2D z-stack projection of the 3D cylindrical gut tube. We therefore excluded the sides from the region-of-interest (Fig. 1a dashed white rectangle) for the orientation analysis. We averaged the coherence and direction of all structures present within the ROI.

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

## Data availability

All data generated or analyzed during this study are included in the manuscript and supporting files (Supplementary Material & Supplementary Data 1). Any remaining information can be obtained from the corresponding author upon reasonable request.

Received: 4 August 2020; Accepted: 7 June 2021;
Published online: 23 June 2021

## References

1. Yntema, C. L. & Hammond, W. S. The origin of intrinsic ganglia of trunk viscera from vagal neural crest in the chick embryo. *J. Comp. Neurol.* https://doi.org/10.1002/cne.901010212 (1954).

2. Le Douarin, N. M. & Teillet, M. A. The migration of neural crest cells to the wall of the digestive tract in avian embryo. *J. Embryol. Exp. Morphol.* 30, 31–48 (1973).

3. Nagy, N. & Goldstein, A. M. Enteric nervous system development: a crest cell's journey from neural tube to colon. *Semin. Cell Dev. Biol.* https://doi.org/10.1016/j.semcdb.2017.01.006 (2017).

4. Hao, M. M. & Young, H. M. Development of enteric neuron diversity. *J. Cell. Mol. Med.* 13, 1193–1210 (2009).

5. Hao, M. M. et al. Enteric nervous system assembly: functional integration within the developing gut. *Dev. Biol.* https://doi.org/10.1016/j.ydbio.2016.05.030 (2016).

6. Uesaka, T., Young, H. M., Pachnis, V. & Enomoto, H. Development of the intrinsic and extrinsic innervation of the gut. *Dev. Biol.* 417, 158–167 (2016).

7. Chevalier, N. R. et al. Embryogenesis of the peristaltic reflex. *J. Physiol.* 597, 2785 (2019).

8. Roberts, R. R. et al. The first intestinal motility patterns in fetal mice are not mediated by neurons or interstitial cells of Cajal. *J. Physiol.* 588, 1153–1169 (2010).
9. Rollo, B. N., Zhang, D., Simkin, J. E., Menheniott, T. R. & Newgreen, D. F. Why are enteric ganglia so small? Role of differential adhesion of enteric neurons and enteric neural crest cells. *F1000Research* **4**, 113 (2015).

10. Foty, R. A., Pfleger, C. M., Forgacs, G. & Steinberg, M. S. Surface tensions of embryonic tissues predict their mutual envelopment behavior. *Development* **122**, 1611–1620 (1996).

11. Broders-Bondon, F., Paul-Gilloteaux, P., Carlier, C., Radice, G. L. & Dufour, S. N-cadherin and β1-integrins cooperate during the development of the enteric nervous system. *Dev. Biol.* https://doi.org/10.1016/j.ydbio.2012.02.001 (2012).

12. Newgreen, D. F., Dufour, S., Howard, M. J. & Landman, K. A. Simple rules for a 'simple' nervous system? Molecular and biomathematical approaches to enteric nervous system formation and malformation. *Dev. Biol.* **382**, 305–319 (2013).

13. Lei, J. & Howard, M. J. Targeted deletion of Hand2 in enteric neural precursor cells affects its functions in neurogenesis, neurotransmitter specification and gangliogenesis, causing functional aganglionosis. *Development* **138**, 4789–4800 (2011).

14. Tanano, A. et al. Structural development of PGP9.5-immunopositive myenteric plexus in embryonic rats. *Anat. Embryol. (Berl).* https://doi.org/10.1007/s00429-005-0453-2 (2005).

15. Murata, H., Yagi, M., Kubota, M. & Ushiki, T. Embryonic development of the myenteric nerve plexus in the rat small intestine as revealed by light and scanning electron microscopy. *Arch. Histol. Cytol.* https://doi.org/10.1679/aohc.66.453 (2003).

16. Pietri, T., Eder, O., Blanche, M., Thiery, J. P. & Dufour, S. The human tissue plasminogen activator-Cre mouse: A new tool for targeting specifically neural crest cells and their derivatives in vivo. *Dev. Biol.* https://doi.org/10.1016/S0012-1606(03)00175-1 (2003).

17. Breau, M. A., Dahmani, A., Broders-Bondon, F., Thiery, J.-P. & Dufour, S. Beta1 integrins are required for the invasion of the caecum and proximal hindgut by enteric neural crest cells. *Development* **136**, 2791–2801 (2009).

18. Püspöki, Z., Storath, M., Sage, D. & Unser, M. Transforms and operators for directional bioimage analysis: A survey. *Adv. Anat. Embryol. Cell Biol.* https://doi.org/10.1007/978-3-319-28549-8_3 (2016).

19. Chevalier, N. R., Fleury, V., Dufour, S., Proux-Gillardeaux, V. & Asnacios, A. Emergence and development of gut motility in the chicken embryo. *PLoS ONE* **12**, e0172511 (2017).

20. Dufour, S., Broders-Bondon, F. & Bondurand, N. in *Neural Surface Antigens—fRom Basic Biology Toward Biomedical Applications* (ed. Pruznak, J.) 153–162 (Elsevier Publishers, 2015).

21. Breau, M. A. et al. Lack of beta1 integrins in enteric neural crest cells leads to a Hirschsprung-like phenotype. *Development* **133**, 1725–1734 (2006).

22. Chevalier, N. R. et al. How tissue mechanical properties affect enteric neural crest cell migration. *Sci. Rep.* **6**, 20927 (2016).

23. Perris, R., Syfrig, J., Paulsson, M. & Bronner-Fraser, M. Molecular mechanisms of neural crest cell attachment and migration on types I and IV collagen. *J. Cell Sci.* **106**, 1357–1368 (1993).

24. Ornestad, M. et al. Foxf1 and Foxf2 control murine gut development by limiting mesenchymal Wnt signaling and promoting extracellular matrix production. *Development* **133**, 833–843 (2006).

25. Graham, M. F. et al. Collagen content and types in the intestinal strictures of Crohn's disease. *Gastroenterology* **94**, 257–265 (1988).

26. Wang, C. et al. Type III collagen is a key regulator of the collagen fibrillar structure and biomechanics of articular cartilage and meniscus. *Matrix Biol.* https://doi.org/10.1016/j.matbio.2019.10.001 (2020).

27. Ranjit, S. et al. Imaging fibrosis and separating collagens using second harmonic generation and phasor approach to fluorescence lifetime imaging. *Sci. Rep.* https://doi.org/10.1038/srep13378 (2015).

28. Davis, E. M. Translocation of neural crest cells within a hydrated collagen lattice. *J. Embryol. Exp. Morphol.* **55**, 17–31 (1980).

29. Perr, H. A. et al. Collagen production by human smooth muscle cells isolated during intestinal organogenesis. *Anat. Embryol. (Berl).* https://doi.org/10.1007/BF00185612 (1992).

30. Groulx, J. F. et al. Collagen VI is a basement membrane component that regulates epithelial cell-fibronectin interactions. *Matrix Biol.* https://doi.org/10.1016/j.matbio.2011.03.002 (2011).

31. Soret, R. et al. A collagen VI-dependent pathogenic mechanism for Hirschsprung's disease. *J. Clin. Invest.* **125**, 4483–4496 (2015).

32. Nagy, N. et al. Collagen 18 and agrin are secreted by neural crest cells to remodel their microenvironment and regulate their migration during enteric nervous system development. *Development* **145**, dev160317 (2018).

33. Beaulieu, J., Jutras, S., Durand, J., Vachon, P. H. & Perreault, N. Relationship between tenascin and t-smooth muscle actin expression in the developing human small intestinal mucosa. *Anat. Embryol.* **188**, 149–158 (1993).

34. Romanska, H., Moscoso, G., Polak, H. & Draeger, A. Smooth muscle differentiation during human intestinal development. *Eur. J. Transl. Myol.* **6**, 13–19 (1996).

35. Wedel, T. et al. Organization of the enteric nervous system in the human colon demonstrated by wholemount immunohistochemistry with special

reference to the submucous plexus. *Ann. Anat.* https://doi.org/10.1016/S0940-9602(99)80122-8 (1999).

36. Wester, T., O'Briain, D. S. & Puri, P. Notable postnatal alterations in the myenteric plexus of normal human bowel. *Gut* https://doi.org/10.1136/gut.44.5.666 (1999).

37. Hanani, M., Grossman, S., Nissan, A. & Eid, A. Morphological and quantitative study of the myenteric plexus in the human Tenia coli. *Anat. Rec.* **295**, 1321–1326 (2012).

38. Graham, H. K., Maina, I., Goldstein, A. M. & Nagy, N. Intestinal smooth muscle is required for patterning the enteric nervous system. *J. Anat.* **230**, 567–574 (2017).

39. Mazzuoli-Weber, G. & Schemann, M. Mechanosensitivity in the enteric nervous system. *Front. Cell. Neurosci.* **9**, 408 (2015).

40. Huycke, T. R. et al. Genetic and mechanical regulation of intestinal smooth muscle development. *Cell* **179**, 90–105.e21 (2019).

41. Vader, D., Kabla, A., Weitz, D. & Mahadevan, L. Strain-induced alignment in collagen gels. *PLoS ONE* https://doi.org/10.1371/journal.pone.0005902 (2009).

42. Long, K. R. & Huttner, W. B. How the extracellular matrix shapes neural development. *Open Biol.* **9**, 180216 (2019).

43. Srinivas, S. et al. Cre reporter strains produced by targeted insertion of EYFP and ECFP into the ROSA26 locus. *BMC Dev. Biol.* https://doi.org/10.1186/1471-213X-1-4 (2001).

44. Madisen, L. et al. A robust and high-throughput Cre reporting and characterization system for the whole mouse brain. *Nat. Neurosci.* https://doi.org/10.1038/nn.2467 (2010).

45. Potocnik, A. J., Brakebusch, C. & Fässler, R. Fetal and adult hematopoietic stem cells require β1 integrin function for colonizing fetal liver, spleen, and bone marrow. *Immunity* https://doi.org/10.1016/S1074-7613(00)80216-2 (2000).

46. Fässler, R. & Meyer, M. Consequences of lack of β1 integrin gene expression in mice. *Genes Dev.* https://doi.org/10.1101/gad.9.15.1896 (1995).

47. Stephens, L. E. et al. Deletion of β1 integrins in mice results in inner cell mass failure and peri-implantation lethality. *Genes Dev.* https://doi.org/10.1101/gad.9.15.1883 (1995).

48. Chen, X., Nadiarynkh, O., Plotnikov, S. & Campagnola, P. J. Second harmonic generation microscopy for quantitative analysis of collagen fibrillar structure. *Nat. Protoc.* **7**, 654–669 (2012).

49. Gailhouste, L. et al. Fibrillar collagen scoring by second harmonic microscopy: A new tool in the assessment of liver fibrosis. *J. Hepatol.* https://doi.org/10.1016/j.jhep.2009.12.009 (2010).

## Acknowledgements

The authors thank Benoit Sorre for granting access to the DSD2 confocal microscope, Melina Durande and François Graner for advice on orientation measurements, Don Newgreen for feedback on our manuscript. We thank the Mouse Facility staff of the Institut Mondor de Recherche Biomédicale. We thank the Imagoseine facility for use of the Zeiss CSU-W1 spinning disk microscope. This research was funded by the CNRS Défi Mecanobiologie “MECHENSDEV” grant, by the IdEx Université de Paris ANR-18-IDEX-0001 grant, by the Agence Nationale de la Recherche ANR GASTROMOVE-ANR-19-CE30-0016-01 grant and by the Institut National Pour la Santé de la Recherche Médicale (INSERM).

## Author contributions

N.R.C. performed experiments, conceptualized & supervised the research, produced figures, wrote the paper with support from S.D. and T.G.; Y.A. performed quantitative analysis and produced figures; A.G. and L.L. performed experiments; T.G. and P.B. performed experiments; S.D. genetically engineered the mouse strains, provided samples and critically revised the paper.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42003-021-02333-5.

**Correspondence** and requests for materials should be addressed to N.R.C.

**Peer review information** *Communications Biology* thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: Simona Chera and Anam Akhtar. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_7.png" /></figure>

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

© The Author(s) 2021