ABSTRACT
Between embryonic days 10.5 and 14.5, active proliferation drives rapid elongation of the murine midgut epithelial tube. Within this pseudostratified epithelium, nuclei synthesize DNA near the basal surface and move apically to divide. After mitosis, the majority of daughter cells extend a long, basally oriented filopodial protrusion, building a de novo path along which their nuclei can return to the basal side. WNT5A, which is secreted by surrounding mesenchymal cells, acts as a guidance cue to orchestrate this epithelial pathfinding behavior, but how this signal is received by epithelial cells is unknown. Here, we have investigated two known WNT5A receptors: ROR2 and RYK. We found that epithelial ROR2 is dispensable for midgut elongation. However, loss of Ryk phenocopies the Wnt5a−/− phenotype, perturbing post-mitotic pathfinding and leading to apoptosis. These studies reveal that the ligand-receptor pair WNT5A-RYK acts as a navigation system to instruct filopodial pathfinding, a process that is crucial for continuous cell cycling to fuel rapid midgut elongation.
INTRODUCTION
The adult human small intestine (SI) is a long coiled tube, measuring 5.5 m on average (Weaver et al., 1991). This remarkable length is crucial for providing the large surface area needed for sufficient nutrient absorption. To achieve this great length, active SI elongation begins in utero, which ensures that the SI acquires an adequate length (over 40% of its adult length) at birth to meet postnatal nutritional demands. Infants that are born with significantly shorter SIs, a disorder known as congenital short bowel syndrome (CSBS), fail to thrive (Hasosah et al., 2008), supporting the significance of SI elongation during the fetal stage. However, the causes of CSBS, as well as the detailed mechanisms underlying the rapid elongation of the fetal SI remain largely unexplored.
The SI is derived from the midgut during embryogenesis. In the mouse embryo, the midgut arises from the middle region of the primitive gut tube by E10 (Zorn and Wells, 2009). Shortly afterwards, this tiny midgut primordium (<1.5 mm) begins to rapidly elongate in two phases. In Phase I (E10.5-E14.5), the midgut is a simple tube with a narrow lumen; it rapidly increases its length by 11-fold and slowly gains girth (Wang et al., 2018). In Phase II (E14.5-E18.5), finger-like villi begin to emerge (Freddo et al., 2016; Walton et al., 2016), quickly changing the epithelial structure and expanding the luminal surface area; elongation continues but at a slower rate (<4-fold in four days) (Wang et al., 2019).
During Phase I, the midgut epithelium is pseudostratified and its accelerating elongation is driven by active cell division, which is tied to a process called interkinetic nuclear migration (IKNM) (Grosse et al., 2011; Wang et al., 2018). During IKNM, nuclei migrate in accordance with the cell cycle, entering S phase near the basal surface and traveling apically to divide (Norden, 2017). The constant nuclear movement allows cells to take turns initiating mitotic cell rounding at the very limited apical surface (Miyata et al., 2015). Although efficient, this peculiar cell division behavior nevertheless faces some challenges. After apical mitosis, daughter nuclei need to quickly exit the crowded apical zone and return to the basal side to synthesize DNA for the next round of division. However, only one of the two daughter cells inherits a basal connection from the mother cell; many heirs can directly exploit this basal connection as a ‘conduit’ to directly transport their nuclei to the basal side. However, more than half of the daughter cells use a ‘pathfinding’ strategy: they must first grow a filopodial protrusion to establish a new basal path and then use that path to deliver their nuclei basally. Importantly, WNT5A, which is secreted by the underlying mesenchymal cells, acts as a guidance cue to facilitate efficient pathfinding. Without WNT5A, filopodial pathfinding is impaired; some cells remain trapped apically, leading to increased apoptosis and a shortened midgut (Wang et al., 2018).
How midgut epithelial cells recognize and respond to the mesenchymal WNT5A cue is unknown. WNT5A plays diverse roles by triggering various intracellular signaling cascades upon the engagement of different receptors (Kikuchi et al., 2012). In this study, we specifically focused on the receptors ROR2 and RYK, both of which are single-pass transmembrane receptor tyrosine kinase (RTK) proteins. ROR2 binds WNT5A through its extracellular frizzled like cysteine-rich domain (Oishi et al., 2003) whereas RYK binds WNT5A via its extracellular WNT inhibitory factor (WIF) domain (Keeble et al., 2006; Yoshikawa et al., 2003).
ROR2 is widely appreciated as a critical receptor of WNT5A signaling (Green et al., 2014). Ror2−/− and Wnt5a−/− embryos share similar defects in many outgrowing structures, although the Ror2−/− phenotype is somewhat milder (Oishi et al., 2003; Suzuki et al., 2003; Takeuchi et al., 2000; Yamaguchi et al., 1999). Importantly, Yamada et al. found that the Ror2−/− midgut is shortened and widened; they proposed that ROR2 directs an epithelial convergent extension process that is required for midgut elongation (Yamada et al., 2010). However, more recent evidence suggests that epithelial convergent extension does not occur in Phase I (Wang et al., 2018). Interestingly, ROR2 mediates WNT5A-induced filopodia formation in several types of cultured cells (Nishita et al., 2006). This raises the possibility that ROR2 participates in gut elongation as an epithelial WNT5A receptor, to sense directional cues and instruct the filopodial outgrowth of pathfinding cells. Still, the Ror2−/− gut is not as short as the Wnt5a−/− gut (63% versus 28% of wild type at E13.5) (Wang et al., 2018; Yamada et al., 2010), suggesting that additional receptors may participate.
RYK is a less studied WNT5A receptor. Similar to Wnt5a- and Ror2-null embryos, Ryk−/− embryos have outgrowth defects in multiple organs (Halford et al., 2000; Kugathasan et al., 2018), but no gut phenotype has been described. Importantly, in mouse forebrain, RYK is responsible for establishing an appropriate axon trajectory guided by WNT5A (Keeble et al., 2006). Thus, RYK is another promising candidate for pathfinding cells to sense the WNT5A cue in the pseudostratified midgut epithelium.
Here, we demonstrate that epithelial ROR2 is not required for midgut elongation, contradicting the previous model (Yamada et al., 2010). Instead, using detailed 3D confocal imaging and 2D live imaging, we detect pathfinding aberrations in the Ryk−/− gut epithelium that mirror those seen in the Wnt5a−/− gut, suggesting that the RYK receptor is crucial for the basally oriented extension of filopodial protrusions in response to WNT5A.
RESULTS
Ror2 is expressed in both mouse midgut epithelium and mesenchyme, but epithelial ROR2 is dispensable for gut elongation
To test the hypothesis that ROR2 acts as a WNT5A receptor to facilitate efficient ‘pathfinding’ in Phase I midguts, we first analyzed the ROR2 expression pattern between E10.5 and E14.5 using immunohistochemistry. Consistent with a previous report (Yamada et al., 2010), ROR2 is present in both midgut epithelium and mesenchyme throughout Phase I (Fig. 1A-E). In the epithelium, ROR2 expression is more intense from E10.5 to E12.5, peaking at E12.5 and decreasing thereafter (Fig. 1A-E, Fig. S1A). In the epithelium, ROR2 is absent at the apical surface and restricted to the basolateral membrane of all cells (Fig. 1F,G′). In the mesenchyme, ROR2 is present on the entire cell membrane of nearly all cells, with varying intensity (Fig. 1H,H′). The expression pattern of Ror2 mRNA was consistent with these findings (Fig. 1I-K″).
To better elucidate the role of epithelial ROR2 in midgut elongation, we took advantage of a Ror2 floxed allele (Ho et al., 2012) that can be deleted in specific lineages using appropriate Cre drivers. Shh is expressed in the epithelium of the primitive gut tube as early as E8.5 (Bitgood and Mcmahon, 1995; Mao et al., 2010) and the ShhCre transgene (Harfe et al., 2004) can direct recombination in all gut epithelial cells (Mao et al., 2010). Therefore, we bred Ror2flox/flox; RosamTmG/mTmG mice with Ror2flox/+; ShhCre/+ to achieve early gut epithelium-specific depletion of Ror2. The specificity of Cre activity was confirmed by crossing ShhCre/+ to RosamTmG/mTmG reporter mice (Fig. 1Q,S) and depletion of ROR2 was confirmed by antibody staining (Fig. 1P,R). Surprisingly, after epithelial Ror2 depletion, no gross defects were observed either in embryos or guts at E14.5 (Fig. 1L-O). Measurements demonstrated that, neither embryo size nor average midgut length was changed at E14.5 (Fig. 1T). Even at E18.5, just before birth, no changes were detected (Fig. 1U). Furthermore, we measured the basal circumference of the midgut epithelial tube (Fig. 1V), the epithelial height (Fig. 1W) and epithelial cell number on cross-sections (Fig. 1X), and no morphological changes were observed. Additionally, no significant changes in epithelial mitotic or apoptotic rates (Fig. S1B-D) were detected after epithelial Ror2 depletion. Altogether, these results indicate that epithelial ROR2 is dispensable for midgut elongation.
Mesenchymal ROR2 is required for midgut elongation before Phase I
The above results suggest that loss of Ror2 in other (non-epithelial) compartments may account for the gut lengthening defects previously seen in Ror2−/− mutants (Yamada et al., 2010). We therefore depleted Ror2 in gut mesenchyme, using a Twist2Cre driver, which is active as early as E9.5 and confined to mesodermal tissues (Yu et al., 2003). The specificity of Twist2Cre activity was confirmed using a mTmG reporter (Fig. 2E′,F′) and mesenchymal ROR2 depletion was confirmed by antibody staining (Fig. 2E,F). Although Ror2−/− embryos were reported to be 80-90% of control size (Takeuchi et al., 2000), Ror2flox/flox; Twist2Cre/+ embryos were only slightly smaller (Fig. 2A,C,G); nevertheless, the two models exhibit similar phenotypes, including shortened limbs, tail, maxilla and mandible (Fig. 2C). Importantly, Ror2flox/flox; Twist2Cre/+ midguts measured 73±5% of control length at E14.5 (Fig. 2B,D,G), very close to the documented Ror2−/− midgut length (60-70% of control length). Thus, mesenchymal Ror2 is important for midgut elongation.
Ror2 is expressed in the gut beginning at E9.5 (Yamada et al., 2010) and transcripts remain present throughout the rest of the fetal stage (Fig. S2A-C). To investigate when ROR2 is required for gut elongation, we induced Ror2 depletion in Ror2flox/flox; UbcCreERT2/+; RosamTmG/+ embryos at different time points (Fig. 2N). Strikingly, the administration of tamoxifen from E9 yielded phenotypes closely resembling those of Ror2−/− mutants and a 25% reduction in midgut length (Fig. 2H-M′,U). However, the same dose of tamoxifen starting at E10.5 or later did not decrease midgut length, although a shorter mandible and mild defects in other structures were seen (Fig. 2O-T′,U). These data suggest that mesenchymal ROR2 plays a role in midgut elongation during the initial primordial stage, but not thereafter.
RYK is required for midgut elongation throughout Phase I
After ruling out an essential role for epithelial ROR2 in midgut elongation, we next investigated another potential WNT5A receptor: RYK. Owing to the lack of antibodies that can specifically detect RYK protein by immunostaining, the RNAscope assay was performed to detect Ryk mRNA in the midgut. From E10.5 to E14.5, Ryk expression was relatively homogenous and Ryk mRNAs were present in every epithelial cell and nearly all mesenchymal cells (Fig. 3A-E′). Next, we obtained Ryk−/−, Ryk+/− and wild-type embryos. Consistent with a previous report (Halford et al., 2000), Ryk+/− embryos were indistinguishable from wild type, but Ryk−/− embryos exhibited a flattened midface and a shortened mandible (Fig. 3F,G,J,K,N,O). Importantly, a significant reduction in midgut length was seen in Ryk−/− embryos and this length deficit was progressive throughout Phase I, from 80% of control length at E10.5, to 63.8% at E14.5 (Fig. 3H,I,L,M,P,Q,S). Of note, reduced midgut length was not due to delayed development or to an overall decrease in embryo size, as the Ryk−/− mutant was only slightly smaller than wild type (not significant before E12.5, 98% of control at E14.5) (Fig. 3R). After E14.5 (beginning of Phase II), Ryk mRNAs were still present in the midgut (Fig. S3A-C) but the fraction of Ryk−/− midgut length versus control length remained ∼64-67% and did not decline further (Fig. S3D). Similarly, Wnt5a−/− midguts are short at E10.5, and length deficits increase throughout Phase I (Wang et al., 2018), but further decrements are not seen in Phase II (Fig. S3E). These analyses suggest that RYK, similar to WNT5A, is important for midgut elongation before and throughout Phase I but not during Phase II.
Though the data above show that epithelial ROR2 is dispensable for gut lengthening, it remains possible that RYK and ROR2 play a redundant role in the epithelium, with the presence of RYK compensating when epithelial ROR2 is lost. Therefore, we further depleted Ror2 (after E10.5) in the Ryk−/− background. However, no additional reduction in midgut length was seen in Ryk/Ror2 double knockouts (Fig. S3F), further confirming that RYK plays a major role in gut lengthening, while epithelial ROR2 is not essential for midgut elongation.
As we have previously shown that active proliferation is the major driver for Phase I midgut epithelial tube elongation, we examined the proliferation marker Ki67 (Fig. 3T,U′) and the mitosis marker pHH3 (Fig. 3V,W′) in control and Ryk−/− midgut epithelium. However, no significant differences were found (Fig. 3X), suggesting that the slower elongation of Ryk−/− midguts was not caused by compromised cell division.
RYK is essential for the proper basal connection of epithelial cells during Phase I midgut elongation
To examine how RYK facilitates midgut elongation at the cellular level, we labeled cell clones in Phase I midgut epithelial tube, using the inducible ShhCreERT2 driver and mTmG reporter system. A very low dose of tamoxifen was given to pregnant dams at E12 to activate membrane eGFP (mG) expression in singly dispersed cells in the midgut epithelium, which then went on to divide and form well-separated large clones at E14.5 (Fig. 4A-C). In control (Ryk+/+ and Ryk+/−) clones, all mG cells touched the apical surface and nearly all were connected to the basal surface, except for a few nascent, post-mitotic cells (1.5%, seven out of 460 cells in 45 clones) (Fig. 4D,E,J). In Ryk−/− midguts, 48.5% of clones appeared normal, with all mG cells connected to both apical and basal surfaces (Fig. 4F,K). However, the remaining clones contained one to three cells lacking a basal connection; some of these cells were even displaced into the lumen (Fig. 4G-I,K). In total, 8.9% (27 out of 305) of all mG cells in Ryk−/− clones were unconnected to the basal surface (Fig. 4J), suggesting that RYK is required for efficient basal connection. Of note, these clonal defects were very similar to those seen in large clones in Wnt5a−/− midguts (Wang et al., 2018).
Some RYK-deficient ‘pathfinding’ daughter cells fail to grow a basally directed filopodial protrusion to make a basal connection
To further characterize how some epithelial cells fail to properly connect basally, we captured the 3D morphology of individual epithelial cells at different phases of the cell cycle by generating one- or two-cell clones in control and Ryk−/− midguts (Fig. 5A). Captured cells were arranged in an order based on their nuclear position. Consistent with our previous report (Wang et al., 2018), as the nucleus moves apically for mitosis (Fig. S4A,B-E′), the cell maintains its basal connection, which gradually turns into a thin basal process. This basal process splits in two during mitosis (Fig. S4A,F-G′); later, one of the two processes retracts in most cases (Fig. S4A,H,H′). During cytokinesis, only one of two daughter cells inherits this basal process, while the other is born without a basal connection (Fig. 5B1,C1,F1,G1,H1). No obvious differences in cell behaviors before cytokinesis were seen between the control and Ryk−/− mutants.
After cytokinesis, the two daughter cells begin to return their nuclei basally via one of two modes. In control Mode I (Fig. 5B), one daughter cell (a) directly uses its inherited basal process to transport the nucleus back – the ‘conduit’ strategy – whereas its sister (b) takes a ‘pathfinding’ strategy. The pathfinding daughter (b) first elongates and adopts a pyriform shape (Fig. 5B1,2). Often, before the nucleus of daughter (a) reaches the basal side, daughter (b) has assembled a very fine filopodium (diameter<0.6 µm, Fig. 5B3,3′) at its most basal tip. This nascent filopodium quickly extends basally and slowly increases in girth (0.5-1.6 µm in diameter), becoming a single thick filopodial protrusion (Fig. 5B4,4′). It finally reaches the basal lamina, establishing a path (13-24 µm long) (Fig. 5B5,5′) for the nucleus of daughter (b) to exit the apical zone. In control Mode II (Fig. 5F), one daughter cell possesses a basal process (Fig. 5F1-4) but does not use it as a nuclear return path; rather both daughters (c) and (d) take the ‘pathfinding’ approach and grow a new filopodium (white arrowheads in Fig. 5F2-4) for nuclear return. Interestingly, the unutilized basal process (blue arrowheads in Fig. 5F1-4) persists only in pairs in which neither daughter's filopodial protrusions have yet reached the basal surface, suggesting that the basal process serves as a tether while both daughters are exploring their path, but quickly disappears once one or both daughters successfully establish a de novo basal connection (Fig. 5F5).
In Ryk−/− midgut epithelium, most daughter pairs (>80%) appear normal in both Mode I (Fig. 5C) and Mode II (Fig. 5G). However, abnormal pairs with a range of defects are also seen (Fig. 5D,E,H-J). Abnormal Mode I pairs are divided into two groups. In the first group (Abnormal I_1, Fig. 5D), pathfinding daughter (b) fails to grow a basal filopodial protrusion, even after the nucleus of daughter (a) has reached the basal side (Fig. 5D3). Some daughters (b) project apically rather than basally (Fig. 5D4) or exit the epithelial layer (Fig. 5D5). In the second group (Abnormal I_2, Fig. 5E), daughter (b) extends multiple short protrusions (<3 µm) (Fig. 5E1-3,5,7) or a branched filopodium (Fig. 5E4,6); its nucleus stays at the apical side of the epithelium (Fig. 5E1-5,7) or moves to the lumen (Fig. 5E8). In both groups, the ‘conduit’ daughter (a) seems unaffected.
Abnormal Mode II pairs are divided into three groups. In the first group (Abnormal II_1, Fig. 5H), both daughters retain a relatively round shape with no protrusion (Fig. 5H2) or multiple tiny protrusions (red arrowheads, Fig. 5H1,3). In the second group (Abnormal II_2, Fig. 5I), two daughters make an oddly shaped extension (red arrowhead) that approaches the basal side. In the third group (Abnormal II_3, Fig. 5J), one daughter appears normal or successfully makes a basal connection (white arrowhead), while the other daughter sends a bifurcated (Fig. 5J2), trifurcated (Fig. 5J1), twisted (Fig. 5J3) or truncated protrusion (Fig. 5J4) (red arrowheads), suggesting that two ‘pathfinding’ daughters from the same mother do not always encounter errors simultaneously without RYK. Of note, although some pathfinding cells fail to properly project basally, they still stay apically connected to their sisters.
RYK facilitates efficient basal tethering of ‘pathfinding’ cells, which is essential for a timely nuclear return
All of the above observations were based on a detailed comparison of the shapes of individual static cells within fixed midguts between controls and Ryk−/− mutants. Considering that filopodial pathfinding is highly dynamic, the behaviors and consequences of those abnormally shaped cells were unclear. To address this, we turned to 2D live imaging of cultured gut explants. However, as the active gut peristaltic movements and the culture setting imposed a resolution limit, the time-lapse images were at a lower resolution. Consistent with our previous report (Wang et al., 2018), daughter pairs in control midguts return in either Mode I or II. Mode I pairs (71%, 65 out of 91 pairs observed) return the ‘conduit’ nucleus quickly (average 48 min) and the ‘pathfinding’ nucleus slowly (average 101 min) (Fig. 6A,R). In Mode II, the nuclei of both daughters return at a similar speed (∼100 min) (Fig. 6I,R).
In Ryk−/− midguts, most daughter pairs (83%, 146 out of 176 observed) appear normal in both Mode I (Fig. 6B,C) and Mode II (Fig. 6J,K); Mode I was still more prevalent (100 out of 146). However, some pairs (17%, 30 out of 176) were clearly aberrant in either Mode I (10%) or Mode II (7%). In abnormal Mode I pairs, ‘conduit’ daughters returned on time (40-80 min, average 58 min), but ‘pathfinding’ daughters all failed to return during the recording (Fig. 6D-H,R). The abnormal Mode I pairs can be further characterized into two groups. In group 1 (6%, Fig. 6D-F), ‘pathfinding’ daughters either only grew tiny, short-lived projections (length<2 µm, indicated by red arrowheads in Fig. 6D,E) or had no visible projections (Fig. 6F). In group 2 (4%, Fig. 6G,H), pathfinding daughter cells formed an obvious protrusion (4-15 µm) but this protrusion ceased to grow further to reach the basal lamina (Fig. 6G,H). Some protrusions were bifurcated, the tips of which pointed in different directions, and were later retracted (Fig. 6H).
Three groups of abnormal Mode II pairs (7%, 12 out of 176) were seen (groups 3-5). In group 3 (3%, Fig. 6L-N), both daughter cells actively extended multiple short-lived, tiny spikes (length<2 µm, Fig. 6L,M), corresponding to the 3D cell shape in Fig. 5H3, or no projections (Fig. 6N), corresponding to Fig. 5H2. In group 4 (2%, Fig. 6O), both daughter cells formed a long basal protrusion, which appeared unstable and retracted later, corresponding to cell shapes in Fig. 5I. In group 5 (2%, Fig. 6P,Q), the two sisters in the same pair behaved differently. One daughter grew a long protrusion temporarily whereas the other did not (Fig. 6P). In another pair, one daughter successfully returned but the other projected apically and failed to return (Fig. 6Q). Additionally, with longer recording times, we observed that those pathfinding daughters that fail to return in time eventually fragmented, a sign of apoptosis, in both Mode I and II (yellow arrowhead in Fig. 6F,N), corresponding to cell shapes in Figs Fig. 5E8 and 5I5. Together with quantitative analyses of nuclear return time (Fig. 6R), we conclude that Ryk depletion perturbs the ‘pathfinding’ return, but not the ‘conduit’ return. Although defective pathfinding daughters showed various dynamics, they all failed to extend a stable filopodial protrusion to the basal surface to form a path for nuclear return. These defects closely mimic those aberrant pathfinding cells reported in Wnt5a−/− midguts (Wang et al., 2018).
Increased apoptosis contributes to the length deficit of the Ryk−/− midgut epithelial tube
Cell proliferation is unperturbed in the Ryk−/− midgut epithelium (Fig. 3T-X) and live imaging indicates that failure of the pathfinding process cell leads to apoptosis (Fig. 6F,N), but can the rate of apoptosis alone account for the length deficit seen? To address this, we examined cleaved-caspase 3 staining. Apoptotic fragments were readily observed accumulated at the apical side of the Ryk−/− midgut epithelium (Fig. 7B-C′) but not of the control (Fig. 7A,A′). Indeed, the apoptosis rate was significantly increased to 4.4±4% in the Ryk−/− midgut epithelium, compared with 0.5%±1% in controls at E14.5; a similar increase was also observed at E12.5 (Fig. 7D). Next, we used a previously published mathematical model (Wang et al., 2018) to estimate the net fold change of the epithelial cell population over time (E11.5-E14.5) (Fig. 7E). This proliferation-based growth curve (c) aligns well with the measured wild-type cell population growth curve (a) acquired during this time. As Ryk−/− epithelial proliferation was not affected, we modeled idealized cell population growth, extrapolating from the measured Ryk−/− cell population at E11.5 (d), and applied 5% (e) and 10% (f) apoptosis, respectively to predict growth curves. The application of 5% apoptosis in the model closely mirrored the Ryk−/− experimental growth curve (b). Indeed, ∼4% apoptosis was measured in Ryk−/− midguts (Fig. 7D). Together, these data suggest that the diminished elongation of the Ryk−/− midgut results primarily from augmented apoptosis, secondary to impaired filopodial pathfinding and inability to establish a basal connection.
DISCUSSION
Dynamic spatial and temporal roles for WNT5A, ROR2 and RYK in midgut elongation
Our findings provide a deeper understanding of how WNT5A signaling temporally and spatially regulates midgut lengthening by engaging different receptors (Fig. 7F). Previous studies, based on global knockouts of Ror2 (Yamada et al., 2010) and Wnt5a (Cervantes et al., 2009), suggested that a WNT5A-ROR2 axis is important in regulating convergence-extension cellular movements within the epithelium that drive gut lengthening. However, the epithelial-specific Ror2 deletion data presented here demonstrate that epithelial ROR2 is dispensable for midgut elongation. Rather, mesenchymal ROR2 is required.
Like ROR2, RYK is expressed in both midgut epithelium and mesenchyme. In Ryk−/− midguts, however, clear phenotypes are limited to the epithelial compartment and consist of defective filopodial protrusions by post-mitotic epithelial cells. Though experiments with constitutive Ryk knockout mice alone cannot definitively exclude an additional role for mesenchymal RYK signaling related to midgut elongation, no apparent defects nor increases in apoptosis were noted in the Ryk−/− mesenchymal compartment. Moreover, mathematical modeling establishes that the level of apoptosis measured experimentally in the epithelium can fully account for the observed differences in the elongation rate of Ryk−/− midgut epithelial tubes. Given that RYK is well documented to act in a cell-autonomous manner (Hing et al., 2020; Kim et al., 2019; Roy et al., 2018; Tourette et al., 2014; Zhong et al., 2011), these data together strongly suggest that RYK function is crucial in the epithelium during Phase I midgut elongation.
In addition to acting in distinct sites (epithelial versus mesenchymal), ROR2 and RYK also act in different temporal windows to control gut lengthening (Fig. 7F). ROR2 acts in the mesenchyme prior to E10.5, suggesting that a mesenchymal WNT5A-ROR2 autocrine signaling may mediate the primordial midgut elongation, though mechanistic details await further investigation. After E10.5, the global loss of Ror2 has very little effect on midgut elongation. In contrast, the length deficit of Wnt5a−/− and Ryk−/− midguts increases during Phase I, indicating that both WNT5A and RYK participate throughout the entire Phase I elongation in a WNT5A-RYK paracrine fashion. In Phase II (E14.5 to birth), when the midgut epithelium undergoes expansion initiated by villus formation, WNT5A, RYK and ROR2 do not appear to participate significantly in further elongation. The dynamic spatiotemporal requirement for distinct WNT5A receptors, as illustrated here, indicates that fetal gut elongation occurs through precise coordination of signals in both epithelial and mesenchymal compartments in a phased manner.
RYK mediates filopodial pathfinding in response to an attractant WNT5A cue during Phase I midgut elongation
The loss of Ryk causes epithelial phenotypes (Fig. 7G) that closely mimic those seen in the Wnt5a−/− midgut during Phase I. First, both mutants exhibit identical filopodial navigation errors in three categories: pathfinding daughter cells fail to extend a filopodium, send multiple short protrusions in a random direction or generate an unstable filopodial protrusion (characterized by branching, twisting or adopting other odd shapes) (Fig. 7G). Second, although pathfinding defects are seen, conduit nuclear return is normal. Additionally, cellular behaviors before G1 Phase, including apical nuclear migration, basal process splitting and cytokinesis, are unaffected. Third, in both mutants, failure of basal re-tethering results in the accumulation of apoptotic cells at the apical side of the epithelium. These data together suggest that RYK mediates WNT5A signaling to control directional filopodial growth in the pseudostratified midgut epithelium.
Still, compared with Wnt5a−/− mutants, a smaller population of epithelial cells lack a basal connection in Ryk−/− midguts (13% in Wnt5a−/− versus 8.9% in Ryk−/− large clones), suggesting that RYK is not the only WNT5A receptor to facilitate filopodial pathfinding. We ruled out ROR2 because the loss of epithelial ROR2 causes no defects in gut length. However, another ROR member, ROR1, is also expressed in the gut (Al-Shawi et al., 2001). Ror1−/− mice are viable and appear normal at birth, but Ror1/2 double deletion results in severe defects that phenocopy Wnt5a−/− mutants, suggesting some redundancy between ROR1 and ROR2 in mediating WNT5A signaling (Ho et al., 2012). This raises the possibility that ROR1 and ROR2 work together to transduce WNT5A cue for pathfinding cells. Alternatively, FZD5 is another candidate. It is highly expressed in early mouse midgut (Liu et al., 2008), and it can bind and respond to WNT5A by internalization from the cell surface (Kurayoshi et al., 2007). On the flip side, several other WNTs can also bind and signal through RYK, including canonical WNT1 and WNT3A (Lu et al., 2004), and non-canonical WNT5B (Lin et al., 2010), but none of these WNTs is present in Phase I midgut (Lickert et al., 2001). Thus, WNT5A is probably the primary cue for RYK in this context. Importantly, however, many pathfinding daughters still project their protrusions correctly and directionally without WNT5A or RYK, suggesting the existence of additional unrecognized guidance cues.
Interestingly, WNT5A, which is transduced by RYK, was previously identified as a repulsive cue in both murine corticospinal tract axon pathfinding (Keeble et al., 2006; Liu et al., 2005) and in Drosophila midline axon crossing (Yoshikawa et al., 2003). However, in the midgut epithelium, filopodial protrusions grow towards the source of WNT5A, indicating that WNT5A acts as an attractant. Therefore, WNT5A-RYK can function as either attractant or repellent, depending on the tissue context or responding cell types.
The pathways acting downstream of WNT5A-RYK to promote directional filopodial growth remain unclear. RYK is an atypical RTK receptor that lacks kinase activity (Hovens et al., 1992). It acts more like a scaffold, promoting the binding of other receptors or intracellular proteins (Roy et al., 2018). With this in mind, three possible targets can be entertained. The first candidate is VANGL2. Biochemical data suggest that WNT5A can promote RYK binding to VANGL2, leading to increased VANGL2 stability (Andre et al., 2012). VANGL2 is enriched at tips of filopodia in neurons (Davey et al., 2016; Dos-Santos Carvalho et al., 2020; Shafer et al., 2011) and can regulate filopodial dynamics (Davey et al., 2016; Dos-Santos Carvalho et al., 2020), possibly by modulating the engagement of N-cadherin with the actin cytoskeleton flow (Dos-Santos Carvalho et al., 2020). The second one is c-SRC, a non-receptor tyrosine kinase. The PDZ-binding domain of RYK can promote the binding and activation of intracellular c-SRC (Petrova et al., 2013; Wouda et al., 2008). Local c-SRC kinase activity promotes filopodial initiation and extension at the tips of growth cones (Robles et al., 2005). Additionally, RYK regulates axon outgrowth by activating Ca2+ signaling in a WNT5A-dependent manner (Li et al., 2009). Future exploration of these possibilities in the context of midgut is needed.
Implications for congenital short bowel syndrome (CSBS)
During Phase I, nearly all epithelial cells undergo continuous cell cycling, leading to an exponential growth of the cell population and rapid gut lengthening. Losing a small proportion of cells in Phase I caused by impaired pathfinding can have a robust effect on overall SI length, as seen in shortened Ryk−/− and Wnt5a−/− SIs at birth, and as further confirmed by mathematical modeling. In this vein, it is noteworthy that mutations in FLNA have been identified as the cause of an X-linked CSBS (Gargiulo et al., 2007) but the mechanism remains unclear. FLNA is an actin-binding protein and can modulate filopodial formation (Chiang et al., 2017; Nishita et al., 2006; Ohta et al., 1999). In melanoma cells, WNT5A can remodel the cytoskeleton by promoting cleavage of FLNA through Ca2+-dependent signaling (O'connell et al., 2009). Although it has not been tested to our knowledge, a potential interaction between RYK and FLNA could unite WNT5A/RYK-mediated filopodial dynamics with the pathogenesis of FLNA-related CSBS in humans.
MATERIALS AND METHODS
Experimental model and subject details
Mice and embryo staging
All mouse work was performed in accordance with guidelines from the University of Michigan's Institutional Animal Care and Use Committee and University Laboratory Animal Medicine as well as Duke University's Institutional Animal Care and Use Committee. Mouse strains in this study were: C57BL6/J (Charles River), Ryk+/− from Dr Steven Stacker at Peter MacCallum Cancer Institute (Halford et al., 2000), Ror2fl/fl (JAX 018354) (Ho et al., 2012), ROSAmTmG/+ (JAX 007576) (Muzumdar et al., 2007), ShhCreERT2/+ (JAX 005623) (Harfe et al., 2004), ShhCre/+ (JAX 005622) (Harfe et al., 2004), Twist2Cre/+ (JAX 008712) (Sosic et al., 2003) and UbcCreERT2/+ (JAX 007001) (Ruzankina et al., 2007). Mice (6 weeks to 8 months old) were maintained in facilities with 12 h light/dark cycles. Both males and females were analyzed. Stage-specific embryos were isolated from timed matings based on the observation of a copulatory plug, which represents E0.5. E14.5 and younger embryos were also staged precisely using the staging system of Theiler (1989).
Method details
Genotyping
Genotyping PCRs were performed using GoTaq Flexi PCR kit (Promega M8295) and Taq DNA Polymerase, native kit (Invitrogen 18038). Primers were made from Integrated DNA Technologies. PCR primer sequences and protocols were as listed on the Jackson Laboratory website and in Table S1.
Tamoxifen preparation and administration
Tamoxifen (Sigma, T5648) powder (200 mg) was dissolved into 10 ml fresh corn oil (Sigma, C8267) to make the 20 mg/ml tamoxifen stock. The stock was further diluted into corn oil for a total volume of 200 μl for each gavage. For Ror2 depletion with UbcCreERT2, pregnant females were gavaged with 2 mg tamoxifen at E9 and 2 mg tamoxifen at E10.5, or 2 mg at E10.5 and 2 mg at E11.5. For labeling large clones in midgut epithelial tube, pregnant females were gavaged with 150 μg tamoxifen at E12, and embryos were harvested at E14.5. For labeling 1- or 2-cell clones, pregnant females were gavaged with 150-200 μg tamoxifen 16 h before harvesting at E13.5.
Midgut fixation and immunostaining
Dissected midguts were fixed in 4% paraformaldehyde (PFA) for 2 h at room temperature or overnight at 4°C. For paraffin sectioning, fixed midguts were embedded into histogel (Thermo Scientific HG-4000-012). Histogel blocks were then dehydrated in ethanol, infused with paraffin, embedded in paraffin blocks and sectioned at 5 μm (Wang et al., 2013). Before staining, sections were deparaffinized in xylene, rehydrated, washed with PBS and then boiled in 10 mM sodium citrate (pH 6) for 10 min for antigen retrieval. During staining, sections were blocked with 4% goat serum, 0.1% TritonX-100 in PBS for 30 min, then incubated with primary antibodies overnight at 4°C, rinsed with PBS and incubated with secondary antibodies for 2 h at room temperature (for antibody details, see Table S1). For cryosectioning, fixed midguts were directly embedded into optimal cutting temperature (OCT). OCT blocks were frozen on dry ice and sectioned at 5 μm. The staining protocol was the same as for the paraffin section staining. The stained samples were mounted with Prolong Gold (Life Technologies P36930) and imaged by an Axio Imager, apotome microscope, using 20× NA 0.8 plan-apochromat or 40× NA 1.3 oil EC Plan-Neofluar objective and recommended settings in Zeiss Zen blue software.
3D cell morphology imaging
Dissected midguts were fixed in 4% PFA for 20 min at room temperature and rinsed with PBS. Fixed midguts were aligned on glass slides, immersed in Focus Clear solution (Cedarlane Labs) for 30 min until transparent, and then mounted with Mount Clear solution (Cedarlane Labs). Imaging was performed using a Nikon A-1 confocal microscope or Zeiss 880 inverted confocal microscope using a 40× oil objective (Nikon Plan Fluor 40× NA 1.3 or Zeiss 40× NA 1.3 oil EC Plan-Neofluar) and 488 and 561 nm laser for excitation with pinhole set for 1 Airy unit for red emission range. Z-stacks were acquired and then 3D reconstructed using Imaris software (Bitplane).
2D live imaging of midgut
After removal from embryos, midguts were immediately placed on the top of a clear transwell filter in a glass-bottom Petri dish filled with 1 ml culture medium (10% fetal bovine serum, 1% penicillin-streptomycin-glutamine in DMEM), as previously described (Costantini et al., 2011; Wang et al., 2018). This culture dish was then incubated at 37°C with 5% CO2 for 2 h to allow midguts to attach to the transwell membrane. Time-lapse imaging was performed on a DeltaVision image restoration microscopy system using a 10×/0.40 UPlanSApo dry or 20×/0.75 UPlanSApo dry objective, Photometrics CoolSnap HQ camera, InsightSSI 7 channel solid-state illumination source, and GFP filter set (475/28ex and 525/50em). The imaging data were deconvolved by SoftWoRx software deconvolution algorithm constrained iterative method and analyzed with Fiji software (Schindelin et al., 2012).
RNAscope in situ hybridization
Dissected midguts were fixed in 10% neutral buffered formalin overnight at room temperature and washed with DNAse/RNAse free water (Gibco) for three changes for a total of 2 h. Tissues were then processed for paraffin sectioning as described above. The in situ hybridization assay was performed according to the manufacturer's protocol (ACD: RNAscope 2.5 HD Detection Reagent–BROWN User Manual). Both Ror2 and Ryk probes were generated by ACD and are commercially available (Ror2, 430041; Ryk, 549981). The Ror2 probe targets nucleotides 2204-3125 of mouse Ror2 mRNA (accession number NM_013846.3); the Ryk probe targets nucleotides 301-1433 of mouse Ryk mRNA (accession number NM_013649.3).
Mathematical modeling
Modeling parameters have been described previously (Wang et al., 2018). Briefly, the fold change in the measured epithelial cell population was determined as Fc, Fc=FxyFZ. Fxy, where Fxy is the measured fold change in actual epithelial cell number on the xy plane (nuclei number on cross-sections), while FZ is the measured fold change in cell number along the z-axis, which can be simplified as the measured fold change in length Fl, assuming a negligible change in cell width and volume over time. In the mathematical model, the growth of the epithelial cell population is driven by symmetrical divisions in a uniform population of progenitor cells. As each cell divides every 16 h (Wang et al., 2018), the cell population at the next cycle, Nt+1, will exhibit geometrical growth: Nt+1=kNt, where k is the net cell population fold change. The modeled fold change is fc, fc=Nt/N0. Without apoptosis, 100 cells become 200 cells at the end of one cell cycle, k=2. When 5% apoptosis occurs, five out of 100 cells are lost and only 95 cells divide, resulting in 190 cells after one cell cycle. Mathematically, the model with apoptosis becomes Nt+1=k(1−a)Nt, where a is the apoptosis rate.
Image analysis
Measurements of the midgut and embryo length were performed on GI or embryo images (from at least three independent littermates) using Fiji/ImageJ. Measurements and quantitation of midgut epithelial tube circumferences, epithelial height, nuclei number, mitosis rate and apoptotic rate were also performed on stained sections (from at least three separate midgut specimens) using Fiji/ImageJ. The clonal size, number of cells with or without a basal connection, and number of clones with 0, 1, 2 or 3 basally unattached cells were counted on 3D reconstructed images using Imaris (from at least three independent experiments). Quantitation of returning time for Mode I and II pairs (from at least three independent experiments) was performed on time-lapse images using Fiji/ImageJ.
Quantification and statistical analysis
All statistical analyses and graphs were generated in Excel or GraphPad Prism 6. All quantifications were collected from at least three independent experiments. Statistical comparisons of data between groups (e.g. control versus experimental condition) were made using unpaired nonparametric tests (Mann-Whitney test). P<0.05 was considered significant. Data are represented as mean or mean±s.e.m. On graphs, P-values were reported as follows: * P<0.05, ** P<0.01, *** P<0.001, **** P<0.0001; n.s., not significant.
Image processing
In Fig. 3F-I, a black rectangle was placed behind each image as the background. In Fig. 5B-D, separate images were adjusted to the same scale, cropped, rotated and combined into one panel. One black rectangle is placed behind each panel as the background. In Fig. 6A-Q, images at different time points of live imaging were adjusted to the same scale, cropped, rotated and combined into one panel. One black rectangle is placed behind each panel as the background.
Acknowledgements
We thank Drs Brigid Hogan, Blanche Capel and Kwabena Badu-Nkansah for editing the manuscript; Drs Daniel Lucas and Yuji Mishina for sharing their UbcCreERT2/+ mouse; Dr Kaitlin Basham from Dr Gary Hammer's laboratory and Sha Huang from Dr Jason Spence's laboratory for sharing reagents; Dr Santiago Schnell for mathematical modeling; Ryan Townshend for RT-PCR guidance; Dr Honglai Zhang for RYK western blot; Dr Matt Kofron at CCHMC; Dr Aaron Taylor at the Microscopy and Image Analysis Laboratory at the University of Michigan; and Drs Benjamin Carlson and Yasheng Gao at the Duke Light Microscopy Core Facility for imaging assistance.
Footnotes
Author contributions
Conceptualization: S.W., D.L.G., T.L.; Methodology: S.W., Y.-H.T., L.C., K.W.; Validation: S.W.; Formal analysis: S.W., A.J.T.; Investigation: S.W., J.P.R., A.J.T., E.B.W., J.U.; Resources: D.L.G., T.L., S.A.S., J.R.S.; Data curation: S.W.; Writing - original draft: S.W.; Writing - review & editing: D.L.G., T.L., S.W., J.P.R., S.A.S.; Visualization: S.W.; Supervision: S.W., D.L.G., T.L.; Project administration: S.W., D.L.G., T.L.; Funding acquisition: D.L.G., T.L.
Funding
This work was supported by the National Institutes of Health (R01 DK089933 to D.L.G. and R01 DK117981 to T.L.). Deposited in PMC for release after 12 months.
Peer review history
The peer review history is available online at https://dev.biologists.org/lookup/doi/10.1242/dev.195388.reviewer-comments.pdf
References
Competing interests
The authors declare no competing or financial interests.
Competing interests
The authors declare no competing or financial interests.