ABSTRACT
Although Wnt signaling is clearly important for the intestinal epithelial homeostasis, the relevance of various sources of Wnt ligands themselves remains incompletely understood. Blocking the release of Wnt in distinct stromal cell types suggests obligatory functions of several stromal cell sources and yields different observations. The physiological contribution of epithelial Wnt to tissue homeostasis remains unclear. We show here that blocking epithelial Wnts affects colonic Reg4+ epithelial cell differentiation and impairs colonic epithelial regeneration after injury in mice. Single-cell RNA analysis of intestinal stroma showed that the majority of Wnt-producing cells were contained in transgelin (Tagln+) and smooth muscle actin α2 (Acta2+) expressing populations. We genetically attenuated Wnt production from these stromal cells using Tagln-Cre and Acta2-CreER drivers, and found that blockage of Wnt release from either epithelium or Tagln+ and Acta2+ stromal cells impaired colonic epithelial healing after chemical-induced injury. Aggregated blockage of Wnt release from both epithelium and Tagln+ or Acta2+ stromal cells drastically diminished epithelial repair, increasing morbidity and mortality. These results from two uncharacterized stromal populations suggested that colonic recovery from colitis-like injury depends on multiple Wnt-producing sources.
INTRODUCTION
Wingless-type MMTV integration site family members (Wnts) are secreted glycolipoproteins engaging in short-range signaling (Willert et al., 2003; Willert and Nusse, 2012). Newly synthesized Wnt proteins are lipidated in endoplasmic reticulum by an acyltransferase, Porcupine (Porcn) (Takada et al., 2006), then transported by Wntless (Wls, or Gpr177 in mice) transmembrane protein for secretion (Bänziger et al., 2006; Bartscherer et al., 2006). Binding of extracellular Wnts to their Frizzled receptors in ligand-receiving cells initiate a signal cascade (Schulte, 2010; Schulte and Bryja, 2007; Wu and Nusse, 2002), inactivating the glycogen synthase kinase 3 (GSK3)-containing cytoplasmic destruction machinery and stabilizing β-Catenin (Cadigan and Peifer, 2009; Huang and He, 2008; MacDonald et al., 2009). Binding of β-Catenin to Tcf factors in nucleus drives a gene program related to cellular growth and proliferation (He et al., 2004; MacDonald and He, 2012; Tamai et al., 2000; Wehrli et al., 2000). Abnormal activation of canonical Wnt pathway in most human colon cancers highlights its profound influence on cellular behavior (Angers and Moon, 2009; Clevers and Nusse, 2012; de Lau et al., 2007; MacDonald et al., 2009; Nusse et al., 2008; Polakis, 2007; Reya and Clevers, 2005). Noncanonical Wnt molecules also activate pathways such as Rho family small GTPases to control cellular migration (Boutros and Mlodzik, 1999; Eaton et al., 1996; Fanto et al., 2000; Habas et al., 2003, 2001; Sakamori et al., 2014; Strutt et al., 1997; Wallingford and Habas, 2005).
The development and maintenance of intestinal epithelium depend on intact Wnt signaling. Removing β-catenin or Tcf4 from mouse intestinal epithelia impaired crypt formation causing lethality (Fevr et al., 2007; Korinek et al., 1998; van de Wetering et al., 2002; van Es et al., 2012). Systemic overexpression of a Wnt antagonist dickkopf homolog 1 induced severe damages to adult intestines (Kuhnert et al., 2004; Pinto et al., 2003). Although Wnt signaling is clearly important, the importance of Wnts themselves is not as clear. Paneth cells were first shown as the primary source of Wnts in small intestine (Gregorieff et al., 2005); ex vivo cultured Lgr5+ intestinal stem cells (ISCs) (Barker et al., 2007), when associated with Paneth cells, showed enhanced enteroid-forming capability (Sato et al., 2011). Enteroids lacking Wnt3 could not be properly passaged (Farin et al., 2012), while enteroids lacking Wls or Porcn failed to grow ex vivo (Kabiri et al., 2014; Valenta et al., 2016). The growth of above epithelial Wnt deficient enteroids could nevertheless be rescued via supplementing stromal cells or stromal Wnt proteins to the culture (Kabiri et al., 2014; Valenta et al., 2016; Farin et al., 2012). These observations, combined with the overall lack of in vivo phenotype in the absence of epithelial Wnt3 (Farin et al., 2012), Vps35 (de Groot et al., 2013), Porcn (Kabiri et al., 2014; San Roman et al., 2014) or Wls (Valenta et al., 2016), collectively suggested that although epithelial Wnts were required to support ISCs ex vivo, their physiological relevance remained unclear in the presence of a potentially redundant stromal source of Wnt. The question regarding the importance of Wnts themselves was further highlighted by the observation that R-spondin was more powerful than Wnts in activating Wnt signaling (Yan et al., 2017).
Existing literature suggested a complexity and functional divergence of various stromal sources of Wnt. The stroma of the intestinal lamina propria consist myofibroblasts, fibroblasts, vascular and visceral smooth muscles, as well as innate and professional immune cells (Powell et al., 2011). Some of these cell lineages express Wnt2b, Wnt4 and Wnt5a (Farin et al., 2012; Gregorieff et al., 2005; Miyoshi et al., 2012). Myh11-CreER-mediated ablation of Porcn in intestinal myofibroblasts alone or in combination with epithelial Porcn deletion did not perturb intestinal epithelial homeostasis (San Roman et al., 2014), marking myofibroblasts a non-essential stromal Wnt source. However, deleting Porcn from myofibroblasts expressing Pdgf receptor α (Pdgfrα) blocked intestinal crypt formation (Greicius et al., 2018). Foxl1-expressing cells were subsequently identified as Pdgfrα+ telocytes that immediately juxtaposed epithelium, produced Wnts, R-spondin 3 and gremlin proteins, but lacked myofibroblast markers. Genetic ablation of Foxl1+ cells or inhibition of Wnt production by Foxl1-expressing cells, disrupted ISC homeostasis in both small intestine and colon (Aoki et al., 2016; Shoshkes-Carmel et al., 2018). McCarthy and co-workers recently identified a gremlin 1-producing CD81+ Pdgfrαlow population, referred to as trophocytes, present immediately below crypts and contributing to ISC maintenance in vivo (McCarthy et al., 2020). In the colon, CD34+; α-smooth muscle actin-negative pericryptal stromal cells expressing Wnt2b, R-spondin 1 and gremlin 1, were shown to support Lgr5+ ISCs ex vivo (Stzepourginski et al., 2016). Preventing Wnt secretion from a heterogenous Gli1+ stromal population disrupted colonic crypt homeostasis (Degirmenci et al., 2018). In small intestines, however, these Wnt2b-enriched Gli1+ cells were proposed to be a reserve Wnt source for epithelial Wnt deprivation. In addition to above homeostasis conditions, Wnt5a+ mesenchymal cells were also shown to promote colonic crypt regeneration after biopsy injury (Miyoshi et al., 2012), whereas Wnt deprivation from intestinal mononuclear phagocytes exacerbated irradiation-induced mucosal damage (Saha et al., 2016).
The above studies, on the one hand, illustrated the ISC-supporting roles of Wnt-producing epithelial and stromal cell types, but, on the other hand, led to perplexing questions. First, if Wnts from specific stromal cell type are obligatory for ISC homeostasis in vivo, then what is the physiological purpose of making epithelial Wnts, especially those from Paneth cells have the prime access to ISCs in the small intestine? Second, transcriptomes of Gli1+ or Pdgfrα+ cells appeared to intersect those of myofibroblasts, yet different Cre drivers yielded distinct, sometimes opposite, outcomes regarding whether, when and where an intestinal phenotype was observed (Degirmenci et al., 2018; San Roman et al., 2014; Greicius et al., 2018; Shoshkes-Carmel et al., 2018). It is also plausible that: (1) instead of the existence of a specific obligatory Wnt-supplying cell type, the quantity of Wnt molecules collected from various adjacent sources in the niche might be crucial for intestinal homeostasis; and (2) epithelial Wnt production may be essential during mucosal injury that disrupts epithelial-mesenchymal interactions. Based on these hypotheses, we describe the under-appreciated role of epithelial Wnts in homeostasis and during various injuries. A single-cell (sc) RNA survey of intestinal stroma revealed that mesenchymal Wnt producers were constrained to cell clusters expressing transgelin (Tagln, also known as SM22α) (Dong et al., 2012) or Acta2. By using Tagln-Cre and Acta2-CreER drivers to abrogate Wnt productions in these stromal subsets, we demonstrate that the abundance of Wnts collected from multiple tissue sources at the regenerating niche may be crucial for colonic healing after injury.
RESULTS
Blocking epithelial Wnts impaired colonic Reg4 cell differentiation
Global or tissue-specific ablation of Wls (Gpr177 in mouse) in various animals caused phenotypes resembling loss of Wnt signaling (Bnziger et al., 2006; Bartscherer et al., 2006; Fu et al., 2009) due to its indispensable function for Wnt secretion (Ching and Nusse, 2006; Das et al., 2012; Port and Basler, 2010). Using constitutive (Villin-Cre) or inducible (Villin-CreER) intestinal epithelial cell (IEC) drivers, we performed Gpr177 ablation to abrogate epithelial Wnt production. Mice lacking IEC-specific Gpr177 did not show histological abnormality in small intestines (Fig. 1A) or alterations in Paneth or mitotic cell numbers in the crypts (Fig. 1B,C, quantified from three mice for each genotype). The deletion of Gpr177 specifically in IECs was further assessed by immunogold labeling of endogenous Gpr177 in wild-type and Gpr177L/−;Vil-Cre intestines. Gpr177+ gold particles were diminished by 90% in endoplasmic reticulum (ER), 91% in non-ER cytoplasmic region and 93% in plasma membrane of Gpr177-deficient Paneth cells (Fig. S1A,B, quantified in Fig. S1C). No obvious change was detected in stromal cells of the same tissues (Fig. S1A,B).
Despite the lack of histological abnormality, Gpr177-deficient colonic epithelium had a significant reduction of Reg4+ cells, previously reported as deep crypt secretory (DCS) cells, in both Vil-Cre- and Vil-CreER-driven homozygotes (Fig. 1D,E). This change has not been previously reported and was specific to Reg4+ cells. No change was found in the number of Muc2+ goblet cells in the same tissues (Fig. 1D,E, Fig. S1D). The number of nuclear β-catenin+ cells in both small intestine and colonic epithelia were equivalent between wild-type and Gpr177-deficient epithelium (Fig. 1F,G, Fig. S1E). Quantitative RT-PCR supported the immunohistochemistry finding that there was no transcriptional alteration in canonical Wnt targets or in Paneth cell genes (Fig. S1F).
Gpr177L/L;Vil-Cre mouse crypts failed to develop into enteroids in regular ENR medium ex vivo (Fig. S1G). The enteroid growth could be partially restored if recombinant Wnt3a or CHIR (the GSK3 inhibitor) was supplemented at seeding (Fig. S1G). We then cultured enteroids from non-induced adult Gpr177L/−;Villin-CreER mice in ENR medium for 4 days, followed by 4-OHT administration for 12 h to induce Gpr177 deletion. Addition of 4-OHT, but not vehicle, caused growth arrest of nearly 90% Gpr177L/−;Vil-CreER enteroids 5 days after 4-OHT treatment (Fig. 1H,I). Recombinant Wnt3a proteins, when added into ENR medium upon 4-OHT withdrawal, significantly improved the survival and growth of Gpr177L/−;Vil-CreER enteroids, judged by enteroid budding morphologies and by EdU incorporation (Fig. 1H,I). These results documented a previously undescribed inhibition of Reg4 cell differentiation in colon upon epithelial Wnt inhibition.
Tagln+ stromal domain encompasses subepithelial Wnt producers
Previous studies investigating sources of stromal Wnts using Myh11, Gli1, Pdgfrα, Cd81 or Foxl1 as markers suggested partial overlap among these cell populations, but the extent was not examined in detail. To gain a systematic understanding of the overall distribution of Wnt producers within intestinal stroma, we conducted a single-cell RNA sequencing (scRNA-seq) analysis of mouse total ileal stromal tissue. Uniform manifold approximation and projection (UMAP) analyses revealed 22 cell clusters encompassing major lamina propria cell types (Fig. 2A). By examining all mouse Wnt genes, we were able to map six major Wnt-expressing cell clusters (Fig. 2A). Cluster 9 was identified as the largest Wnt4- and Wnt5a-expressing cluster co-expressing Foxl1, Pdgfra and Gli1 (Fig. 2B). This cluster also robustly expressed Tagln, Acta2 and Myh11, genes that are expressed by fibroblast, myofibroblast and smooth muscle cells. Other Wnt4- and Wnt5a-expressing clusters are 5, 11, 14 and 21, where a small number of Foxl1+ cells were found. Clusters 5, 11 and 21 co-expressed high levels of Pdgfra while cluster 14 did not. Gli was highly expressed in 5 and 14, but not as highly in clusters 11 and 21. However, all these clusters were marked by Tagln and Acta2 (Fig. 2B). In contrast to the Wnt4- and Wnt5a- expressing clusters, cluster 15 was defined by Wnt2b-expressing Pdgfralow trophocytes, as observed by the presence of Cd81 and Grem1 (McCarthy et al., 2020), and the absence of Foxl1. In addition, this cluster expressed lower abundance of Tagln, Acta2 and Myh11 (Fig. 2B). Examination of Tagln-expressing Wnt-producing clusters 5, 9, 11, 14 and 21, in light of Pdgfra, Gli1 and Foxl1 suggested that they were substantially overlapped by Tagln-, Acta2- and Myh11-expressing cell populations (Fig. 2C).
Compared with the population that do not express Wnt2b, Wnt4 or Wnt5a, cells expressing a high abundance of these ligands also expressed abundant Tagln and Acta2 compared with Gli1 and Foxl1 (Fig. 2D). We then used a single-cell deep constrained clustering (scDCC) approach, a newly developed deep learning model that guides the clustering by giving weight to marker genes (here, we used Wnt5a and Wnt4, respectively). These Wnt-based clusterings demonstrated a high abundance of Tagln, Acta2, Pdgfra, Myh11, Gli1 and Foxl1, and a low abundance of Cd81 in Wnt5a-expressing and Wnt4-expressing stromal clusters (circled) (Fig. S2A,B).
Compared with the heterogeneous nature of Gli1+ cells (Degirmenci et al., 2018), Foxl1+ telocytes represented a specific cell population shown to support ISC homeostasis (Shoshkes-Carmel et al., 2018). The majority of Foxl1-expressing cells also expressed Tagln (Fig. 2E), and Tagln was not detected in 16 Foxl1+ cells. Examining the transcriptome of Foxl1+/Tagln− and Foxl1+/Tagln+ cells revealed their high expression of Wnt4, Wnt5a, Pdgfra and Gli1 (Fig. 2F). Gene-specific analysis (GSA) identified only six genes differentially expressed between Foxl1+/Tagln− and Foxl1+/Tagln+ populations, while 614 genes were differentially expressed between Foxl1−/Tagln+ and Foxl1+/Tagln+ populations (Fig. S2C-G), suggesting that Foxl1 defines a strict subset of Tagln+ cell population. Overall, these analyses revealed that Tagln+ domains encompassed the major stromal Wnt producers, and there was a high degree of overlap among reported stromal cell lineages, except for Pdgfralow trophocytes (McCarthy et al., 2020).
Tagln-Cre and Acta2-CreER targeted intestinal stromal domains
Examination of endogenous Tagln expression revealed a mesenchymal domain throughout the intestinal tract (Fig. 3A), consistent with its reported expression in fibroblasts, myofibroblasts and smooth muscles (Katajisto et al., 2008; Kuhbandner et al., 2000; Dong et al., 2012). In duodenum and ileum, we found that subset of Tagln-expressing cells was in immediate proximity of crypt cells, wrapping the epithelial glands and extending into villus mesenchymal region (Fig. 3A). A similar pattern was found in colon (Fig. 3B).
We crossed Tagln-Cre mouse (Boucher et al., 2003; Holtwick et al., 2002) to Rosa26R reporter mice (Srinivas et al., 2001), and found that Tagln-Cre activated reporter expression in a pattern similar to endogenous Tagln (Fig. 3C, compare with A). We next stained for endogenous Wnt proteins in Tagln-Cre;Rosa26REYFP mouse intestinal tissues. In both small intestine and colon, a subset of Tagln-Cre expressing cells (EYFP+) co-expressed Wnt2b (white arrows, Fig. 3D,E). Some Wnt2b-expressing Tagln+ cells showed a large cell body with long processes and were located underneath the epithelium (Fig. 3D), resembling what was described for telocytes. Close examination revealed Wnt2b+ vesicular puncta in the vicinity of cells with strong cytoplasmic Wnt2b staining (arrowheads in Fig. 3E, Fig. S3A). A subset of Tagln-Cre-expressing cells expressed endogenous Wnt5 (Fig. 3F,G), and tended to show a different localization and cell shapes. Close examination revealed both perinuclear and vesicular localization of Wnt5 (Fig. S3B). Examination of GFP positivity for Wnt2b+ or Wnt5+ cells suggested that the majority of them were labeled by Tagln-Cre (Fig. 3H), consistent with scRNA-seq results (Fig. 2C,D).
In addition to Tagln-Cre, we also examined the Acta2-CreER driver (a knock-in allele referred to as Acta2-CreER hereafter) (Manieri et al., 2015) and observed a slightly narrow Cre-expressing stromal domain, with scattered EYFP+ cells detected in the lamina propria of villi in Acta2-CreER;Rosa26REYFP mice (Fig. 3I). The relatively constrained Acta2-CreER activity reflected a potential mosaic recombination by this tamoxifen-inducible line, as opposed to the constitutive expression of Tagln-Cre. We then used these two stromal drivers to assess the contribution of Wnts produced by Tagln+ and Acta2+ cells to intestinal homeostasis, which have not been studied before.
Preventing Wnt production in epithelium and in Tagln+ or Acta2+ cells diminished colonic Reg4 cells
We crossed Gpr177 flox mice to Tagln-Cre or Acta2-CreER mice. Gpr177L/L; Tagln-Cre mice did not display global intestinal defects. Similarly, tamoxifen-treated Gpr177L/L; Acta2-CreER mice did not exhibit any defects, observations reminiscent of a previous report using Myh11-Cre to delete Porcn in myofibroblasts (San Roman et al., 2014). We then developed Gpr177L/−;Vil-Cre;Tagln-Cre as well as Gpr177L/−;Vil-Cre;Acta2-CreER mice to remove Gpr177 from both epithelium and Tagln+ or Acta2+ stromal cells. In situ hybridization for Gpr177 showed effective gene ablation in both epithelial and stromal compartments (Fig. S4A,B).
Gpr177L/−;Vil-Cre;Tagln-Cre mouse intestines appeared histologically normal (Fig. 4A) but with reduced crypt base columnar cells (CBCs) marked by Olfm4 by in situ (Fig. 4C, quantified in 4B), Paneth cells (Fig. 4D, quantified in 4B) and proliferative cells marked by 30 min EdU labeling (Fig. 4E, quantified in 4B) or by PCNA staining (Fig. S4C). No change was detected for nuclear β-catenin+ cells (Fig. 4F, quantified in 4B). Examination of Gpr177L/−;Vil-Cre;Tagln-Cre mouse colons showed a notable alteration of glandular morphology (Fig. 4G), accompanied by a diminished number of Reg4+ cells (Fig. 4H, quantified in 4I). Goblet cells remained unchanged in the same tissue (Fig. 4H, quantified in 4I, also see Fig. S4D).
A similar reduction of colonic Reg4+ cells was observed in Gpr177L/−;Vil-Cre;Acta2-CreER mice (Fig. 4H), but to a less extent compared with Gpr177L/−;Vil-Cre;Tagln-Cre mice (Fig. 4I). However, Gpr177L/−;Vil-Cre;Acta2-CreER mice did not show other changes in CBCs, Paneth cells or proliferative cells (Fig. 4A-F). Quantitative RT-PCR analysis showed a significant reduction in Gpr177 transcripts in Gpr177L/−;Vil-Cre;Acta2-CreER intestine (Fig. S4E), suggesting the lack of phenotypic difference was not due to an inefficient deletion.
Quantitative RT-PCR analysis detected a significantly reduced Mmp7 in Gpr177L/−;Vil-Cre;Tagln-Cre intestines (Fig. 4J), and reduced Defa5 in Gpr177L/−;Vil-Cre;Acta2-CreER intestines (Fig. 4K). Mmp7 and Defa5 represented signature markers of mature Paneth cells. These results suggested an additive impact of Wnts from epithelium and Tagln+ or Acta2+ stromal cells on Reg4+ cell differentiation.
Blocking Wnts released from the epithelium and from Tagln+ cells modestly affected regeneration following irradiation
Whole-body irradiation rapidly eliminates proliferating intestinal stem cells (ISCs) (Potten, 1977) and elicits epithelial regeneration (Yu et al., 2018; Zhang et al., 2020). We asked whether Wnts produced in epithelium and Tagln stromal cells participated in irradiation-induced regeneration. Adult control, Gpr177L/L;Vil-Cre and Gpr177L/L;Vil-Cre;Tagln-Cre mice (all littermates) received 12 Gy whole-body irradiation (n=6 for each genotype, from at least three different litters). All mice showed a similar trend of body weight loss starting from day 3 (Fig. 5A). On day 5, wild-type and Gpr177L/L; Vil-Cre mice showed a weight gain whereas Gpr177L/L;Vil-Cre;Tagln-Cre mice displayed the smallest weight gain (P<0.05, Fig. 5A). Upon sacrifice, Gpr177L/L;Vil-Cre;Tagln-Cre intestines showed a reduced number of regenerative crypts (Fig. 5B, quantified in 5D), proliferative cells (Fig. 5C,D), CBCs and Paneth cells (Fig. 5D,E). Gpr177L/L;Vil-Cre mice exhibited only a moderate reduction of Paneth cells (P<0.05, Fig. 5D,E). Nuclear β-catenin+ cell numbers were significantly reduced in crypts of regenerating Gpr177L/L;Vil-Cre;Tagln-Cre intestines (two morphologically different regions shown in Fig. 5F,G). The observed modest reduction of body weight and epithelial regeneration in Gpr177L/L;Vil-Cre;Tagln-Cre mice collectively indicated a resilience to abrogated Wnt production from the epithelium and Tagln+ stromal cells under the context of irradiation.
Colonic healing from DSS-induced mucosal injury requires Wnts from epithelial cells, and Tagln+ and Acta2+ stromal cells
Compared with irradiation that primarily eliminated cycling ISCs, dextran sulfate sodium (DSS) induces mucosal demolishing damage, characterized by epithelial death, ulceration, impairment of glandular architecture and inflammation. Upon DSS withdrawal, colonic epithelia heal through regeneration (Yu et al., 2020). As Wnt signaling was shown to play a role in DSS-induced epithelial recovery (Koch et al., 2011), we further asked whether colonic epithelial restoration required Wnts from epithelial or Tagln+ and Acta2+ cells. We treated Gpr177L/L;Vil-Cre (n=10) and Gpr177L/L;Vil-Cre;Tagln-Cre mice (n=6) with 2.5% DSS in drinking water for 7 days and returned them to regular water for another 7 days to allow recovery (Fig. 6A). We used littermates that were Cre negative, as well as mice expressing Cre only, and tamoxifen-injected wild-type and Vil-CreER mice (non-littermates) as controls (n=13). Strikingly, compared with control mice, Gpr177L/L;Vil-Cre and Gpr177L/L;Vil-Cre;Tagln-Cre mice showed significantly more body weight losses (Fig. 6B). All control mice and most Gpr177L/L;Vil-Cre mice survived the entire 14-day procedure, whereas Gpr177L/L;Vil-Cre;Tagln-Cre mice exhibited mortality from day 9 (Fig. 6C), an early point of recovery. Gpr177L/L;Vil-Cre;Tagln-Cre survivors did not exhibit body weight gain, which was observed for wild-type and Gpr177L/L;Vil-Cre mice during recovery phase (Fig. 6B). Gpr177L/L;Vil-Cre;Tagln-Cre survivors showed shortened colons upon sacrifice on day 14 (Fig. 6A).
Pathological evaluation revealed greater epithelial damages in Gpr177L/L;Vil-Cre mice than wild type, and the worst colitis in Gpr177L/L;Vil-Cre;Tagln-Cre mice on day 14 (Fig. 6D,E). Colitis scoring based on ulceration, crypt death, immune cell infiltration, thickness of colonic wall and loss of goblet cell, reached the same conclusion for day 11 (Fig. S5B).
Control mouse colonic mucosa collected on day 11 already exhibited prominent wound-healing activities, including the formation of nascent crypts flanking the wound, epithelial proliferation and migration towards the surfaces of the wound (red brackets designate ulceration, Fig. 6D). Control mouse mucosa showed near-complete healing on day 14, exemplified as a restoration of goblet cell-filled glandular structures (Fig. 6F). In contrast, Gpr177L/L;Vil-Cre and Gpr177L/L;Vil-Cre;Tagln-Cre mice had significantly larger and more ulcerations at the beginning and end of recovery (Fig. 6D,G). On day 11, Gpr177L/L;Vil-Cre;Tagln-Cre mouse colons showed continuous ulceration, massive immune cell infiltration and a severely thickened colonic wall (Fig. 6D). On day 14, Gpr177L/L;Vil-Cre;Tagln-Cre survivors had scattered colonic crypts with diminished goblet cell differentiation (Fig. 6D,F).
The impact of Gpr177 removal from intestinal epithelium or Tagln+ stromal cells on colonic epithelial resolution demonstrated here was never described. We therefore also performed DSS experiments on Gpr177L/L;Vil-CreER mice to validate this finding. Tamoxifen-induced Gpr177L/L;Vil-CreER mice showed a similarly exacerbated colitis morbidity to that seen in Gpr177L/L;Vil-Cre mice (Fig. 6B,H,I). The above changes were not due to toxicity of Cre or CreER, as mice expressing Cre or CreER only did not exhibit difference from wild-type mice (Fig. S5A). In addition, we used a single administration of tamoxifen (75 mg/kg) in all our CreER experiments, and this was well below the reported toxicity-inducing dose (Bohin et al., 2018; Huh et al., 2010). Indeed, wild-type mice or CreER mice treated with this tamoxifen dose did not show difference in terms of body weight recovery in response to DSS treatment (Fig. S5A).
As Acta2-CreER targeted a relatively restricted stromal domain (Fig. 3I), we asked whether Wnts produced by Acta2+ stromal cells also contributed to colonic wound healing. Surprisingly, both Gpr177L/L;Acta2-CreER (n=6) and Gpr177L/L;Vil-Cre;Acta2-CreER (n=6) mice showed increased morbidity (Fig. 6B,C). We noted that Gpr177L/L;Vil-Cre;Acta2-CreER mice showed a worse body weight recovery than Gpr177L/L;Acta2-CreER mice, starting from day 9 (Fig. 6B). The slower recovery was consistent with their greater colitis shown at pathological levels (Fig. 6H,I). These analyses suggested that Wnts produced by epithelial cells, and by Tagln+ and Acta2+ cells were all important for colonic wound healing and the disease severities were proportionally correlated with additive Wnt inhibition from different tissue compartments.
Wnts from wound-associated Tagln+ cells are required for β-catenin nuclear localization in regenerative epithelium
DSS treatment disrupted the colonic mucosal microenvironment by severely demolishing glandular structure and epithelial-stromal interaction. We examined the localization of Tagln+ cells in DSS-treated colonic mucosa 4 days after withdrawing DSS treatment, a time point when epithelial recovery was prominent (Fig. 7A). Compared with colons at steady states (Fig. 3A), DSS-treated colons contained increased number of Tagln+ cells surrounding the hyperplastic crypts (white arrows, Fig. 7B). Immediately flanking the ulcered region (asterisk in Fig. 7A), Tagln-expressing cells were detected at the base of nascent crypt (white arrows, Fig. 7C) and underneath the single layer of wound-covering epithelium (white arrows, Fig. 7D).
Within these regenerating epithelia associated with Tagln+ cells, epithelial cells with nuclear β-catenin were found (red arrows, Fig. 7E-G). We found a significant reduction in such nuclear β-catenin+ regenerative epithelial cells in Gpr177L/L;Vil-Cre (Fig. 7H-J) and Gpr177L/L;Vil-Cre;Tagln-Cre colons (Fig. 7K-M). Compared with Gpr177L/L;Vil-Cre colons, Gpr177L/L;Vil-Cre;Tagln-Cre colons showed a further reduction of nuclear β-catenin+ epithelial cells (Fig. 7N). As Tagln+ cells were prominently observed near the injured epithelia (Fig. 7O,P), these changes were apparently caused by lack of Wnt production in the niche, rather than lack of Tagln+ cells. We further examined three canonical Wnt targets, CD44, Myc and Axin2, by immunostaining (Fig. S6A-C). Gpr177L/L;Vil-Cre;Tagln-Cre colons showed the lowest CD44 and Axin2 in regenerative epithelia (Fig. S6B and C). These results collectively suggested that, in absence of epithelial Wnts, inhibition of Wnt production from Tagln+ or Acta2+ cells diminished nuclear β-catenin accumulation and canonical Wnt signaling in regenerative epithelial cells.
DISCUSSION
By ablating Wnt secretion from two distinct stromal populations via Tagln-Cre and Acta2-CreER drivers, or from epithelium using Villin-Cre or Villin-CreER, we demonstrated that diminishing Wnt production from either epithelial or these particular stromal subsets compromised the regenerative capability of colonic epithelium. We further showed that aggregated Wnt deprivation from epithelium as well as from Tagln+ or Acta2+ subsets exacerbated the mucosal damage. In addition to documenting a previously undescribed role of epithelial Wnts in colonic Reg4 cell differentiation and injury-induced regeneration, our findings supported a model that the accumulated Wnt molecules from multiple tissue sources collected at ISC niche may be crucial for injury-induced regeneration.
Literature documenting the growth-driving effects of various Wnts on intestinal stem cell renewal and enteroid growth. Farin et al. and Valenta et al. showed that the addition of Wnt2b was sufficient to restore growth of Wnt3Δ/Δ enteroids and Wls-deficient enteroids (Farin et al., 2012; Valenta et al., 2016). Wnt5a treatment also resulted in the growth of Wls-deficient crypts into small spheroids (Valenta et al., 2016). Kabiri et al. showed that stromal cells promoted survival and growth of Porcn-deficient enteroids and supported epithelial regeneration in IEC-specific Porcn-deficient mice post-radiation injury (Kabiri et al., 2014). Deleting Porcn in telocytes reduced epithelial proliferation in the small and large intestine (Shoshkes-Carmel et al., 2018). Furthermore, mice with Wnt5a deleted in telocytes showed shortened intestines, an observation that was attributed to increased apoptosis in the neonatal mice (Kondo and Kaestner, 2021). Upon infection with rotavirus (RV), an enteric pathogen that affects differentiated cells, the number of CBCs expanded and the proliferating cells migrated faster. Mice lacking epithelial Wnt did not exhibit a similar response post-infection, suggesting the role of epithelial Wnts in inducing a regenerative response to infection (Zou et al., 2018). In et al. (2020) demonstrated a mosaic expression of Wnt2b in specific epithelial cells of human colonic crypts and undifferentiated colonoids. In response to E. coli-secreted cytotoxin EspP-induced injury to colonoids, there was an increase in Wnt2b-expressing cells (In et al., 2020). The infected colonoids failed to survive in the absence of epithelial Wnts when treated with a Porcn inhibitor, suggesting the importance of epithelial Wnts for colonic regeneration in response to injury. In the context of these studies, we delineated the functional contribution of Wnts from two uncharacterized stromal populations (Tagln+ and SMA+) to the homeostasis of small and large intestines.
In homeostasis, we found that either constitutive (by Villin-Cre) or inducible (by Villin-CreER) deletion of epithelial Wls/Gpr177 led to reduction of Reg4+ deep crypt secretory (DCS) cells in the colon, an observation that was not reported before. This change was specific to DCS cells as goblet cells and enterocyte numbers were not affected. DCS cells were originally identified as Paneth cell equivalents that support colonic stem cells (Sasaki et al., 2016). Compared with ablation of epithelial Wnt alone, removing Wnt productions from both epithelium and Tagln+ cells resulted in a near-complete inhibition of DCS differentiation, suggesting that the differentiation or maturation of colonic DCS cells were proportionally sensitive to a cumulative loss of Wnt availability in the niche. Previous DCS ablation studies documented their indispensable role in maintaining colonic stem cell homeostasis in vivo (Sasaki et al., 2016). Interestingly, the different extent of DCS cell loss in Gpr177L/L;Vil-Cre and in Gpr177L/L;Vil-Cre;Tagln-Cre colons positively correlated with different degrees of colonic damage manifested by these DSS-treated homozygotes. Wnt-dependent DCS differentiation identified here offered a mechanism that might contribute to the impaired colonic epithelial repair after Wnt deprivation from epithelial and Tagln+/Acta2+ cell populations.
Bulk RNA-seq and scRNA-seq studies used sorted Foxl1+ or Gli1+ cells for transcriptomic analysis (Degirmenci et al., 2018; Shoshkes-Carmel et al., 2018). By contrast, our scRNA study surveyed the intact stroma without selecting or filtering out a specific subepithelial population. The analysis revealed substantial overlaps among Pdgfrα+, Gli1+ and Foxl1+ cell clusters, all of which robustly co-expressed Tagln, Acta2 and Myh11. Interestingly, our scRNA analysis also confirmed the Wnt2b-expressing Cd81+ Pdgrfalo cluster (McCarthy et al., 2020), which does not express abundant Tagln, Acta2 and Myh11. Ninety three percent of Foxl1+ cells expressed Tagln, a gene encoding an actin-binding protein. Tagln-expressing cells define the widest spectrum of mesenchymal domain, including vascular and visceral smooth muscles, myofibroblasts and fibroblasts (Dong et al., 2012). This mesenchymal cell population has previously been shown to regulate intestinal epithelial proliferation and polyposis (Katajisto et al., 2008). We demonstrated that endogenous Tagln+ cells wrapped around crypts in steady-state small intestine and colon. In DSS-treated colon, Tagln+ cells were immediately associated with the regenerating epithelium flanking the wound. Based on our Cre reporter analysis, Tagln-Cre conceivably targeted more Wnt-producing cells than Acta2-CreER, possibly owing to their constitutive versus inducible expression features. In line with this, combined Wnt inhibition by Villin-Cre and Tagln-Cre caused the strongest colonic damage, reflected by the highest morbidity and lethality, and a continuous loss of epithelium even in late recovery stage. In addition to chemical injury, we also observed moderately impaired regeneration following irradiation. The stronger impact elicited by DSS could be due to its robust and comprehensive epithelium-damaging effect that was not limited only to cycling ISCs in the way that the irradiation was.
Although Acta2-CreER marked a slightly constrained stromal subset, compared with controls, mice lacking Wnt production from these cells also exhibited a pronounced colonic damage, suggesting that even these smaller stromal subsets actively contributed to epithelial regeneration during injury. Interestingly, this same Acta2+ cell population was shown to migrate to colonic wound sites and, when mucosally transplanted, improved colonic epithelial healing after biopsy injury (Manieri et al., 2015). Our data suggested that the presence of Acta2+ cells near epithelial wound might promote colonic healing through Wnt-dependent mechanisms.
Blocking Wnt production from Foxl1+ cells caused severe damage to epithelium in both small intestine and colon in unchallenged mice (Shoshkes-Carmel et al., 2018), whereas blocking Wnts from Gli1+ cells resulted in collapse of colonic crypts and death (Degirmenci et al., 2018). Our scRNA analysis suggested that Foxl1+ cells represented a fraction of Tagln+ domain, and that the majority of Foxl1+ cells expressed Tagln. In homeostasis, we did not observe pronounced epithelial damage in small intestines or colons when using various single or combinatory Cre strategies. Our study aligns with the accumulating evidence suggesting that Foxl1+ cells may not function as canonical Wnt producers, in spite of expressing Wnt4 and Wnt5a (McCarthy et al., 2020). It is plausible that the Cd81+ Pdgfrαlow trophocytes, which have been shown to maintain intestinal stem cells in vivo (McCarthy et al., 2020), may not be effectively deleted by Tagln-Cre and Acta2-CreER drivers reported in this study. Nevertheless, we revealed that blocking Wnt from either the Tagln+ or Acta2+ compartment attenuated colonic epithelial regeneration upon DSS injury. Thus, the non-primary Wnt-producing stromal cells also played a role in supporting injury-induced colonic epithelial regeneration. Taken together, we conclude that the differentiation and regeneration of colonic epithelium are sensitive to Wnt abundance in ISC niche collected from multiple tissue sources.
MATERIALS AND METHODS
Mice and treatment
Gpr177LoxP/LoxP (Gpr177L/L), Gpr177+/−, Acta2-CreER, Vil-Cre and Villin-CreER mice have been described previously (el Marjou et al., 2004; Fu et al., 2011, 2009; Madison et al., 2002; Manieri et al., 2015). Tagln-Cre mice (Boucher et al., 2003; Holtwick et al., 2002) were obtained from the Jackson Laboratory (stock 004746). For tamoxifen treatment, adult mice were intraperitoneally injected with a single dose of tamoxifen (75 mg/kg body weight). For total body irradiation, mice were subjected to 12 Gy γ-irradiation, followed by monitoring of body weight daily. For 5-ethynyl-2′-deoxyuridine (EdU) labeling, EdU (2.5 mg/kg) (Thermo Fisher Scientific, C10338) was injected 30 min before mice were sacrificed at the indicated days (Yu et al., 2018). Mice were maintained on a C57BL/6 and 129 mixed background. All experiments were performed on littermates unless explicitly stated in results. Data were quantified from a minimum of three different litters. Numbers of animals used were described in figure legend. All experimental procedures with mice were approved by Rutgers Institutional Animal Care and Use Committee.
Immunohistochemistry and immunofluorescence staining
Procedures and antibodies for immunohistochemistry and immunofluorescence have been described previously (Das et al., 2015; Feng et al., 2017). In brief, the intestinal tissues were fixed in 4% PFA at 4°C overnight, embedded in paraffin wax and cut into 5 μm sections. For immunofluorescence staining, the slides were rehydrated, subjected to antigen retrieval (citrate acid buffer, pH 6) and incubated in the blocking buffer (PBS containing 0.1% Triton X-100 with 2% normal serum and 2% BSA) at room temperature for 1 h. Then the slides were probed with the primary antibodies at 4°C overnight. Lysozyme (Biogenex, AR024-5R, 1:1), GFP (Invitrogen, A11122 and Abcam, ab6673, 1:500), E-cadherin (BD Transduction Laboratories, 610182, 1:500), pHH3 (Millipore, 06-570, 1:500), Reg4 (R&D Systems, AF1379, 1:100), Muc2 (Santa Cruz Biotechnology, sc-15334, 1:200), EdU (Thermo Fisher Scientific, C10338), Wnt5a/b (Cell Signaling Technology, 2530s), Wnt2b (Abcam, ab50575) and Sm22 (Abcam, ab14106, 1:250) staining was performed according to the manufacturer's instructions. The next day, slides were washed three times with PBS and then incubated with Alexa-labeled secondary antibody (Thermo Fisher Scientific, 1:1000) in the dark for 1-2 h. The samples were counterstained with Topro-3 (Thermo Fisher Scientific, T3605, 1:300) and mounted with the ProLong Gold anti-fade mounting media (Thermo Fisher Scientific, P36930). Images were taken using a Zeiss laser scanning microscope. For immunohistochemistry staining, the slides were hydrated, subjected to antigen retrieval (citrate acid buffer, pH 6) and quenched for endogenous peroxidase. The samples were incubated in the blocking buffer at room temperature for 1 h and then incubated with Lysozyme (Biogenex, AR024-5R, 1:8), β-catenin (BD Transduction Laboratories, 610153, 1:1000), PCNA (Santa Cruz Biotechnology, sc-56, 1:1000), Myc (Millipore Sigma, CBL430, 1:100), CD44 (Cell Signaling Technology, 3570, 1:100) and Axin2 (Thermo Fisher Scientific, MA5-32646, 1:100) antibodies at 4°C overnight. The next day, slides were washed three times with PBS and then incubated with biotin-labeled secondary antibody (Vector Laboratories, 1:500) for 1-2 h. Staining was developed using an ABC Kit (Vector Laboratories, PK-6100) and followed by the DAB kit (Vector Laboratories, SK-4100). The samples were counterstained with Hematoxylin (Vector Laboratories, H-3404) and mounted with cytoseal mounting media (Thermo Fisher Scientific, 8310-4). Images were taken under the Nikon (TE2000-U) light microscope.
Quantitative real-time PCR
Quantitative RT-PCR has been described previously (Das et al., 2015; Gao and Kaestner, 2010; Gao et al., 2009; Sakamori et al., 2012, 2014), with a list of primers provided in Table S1. Threshold cycle (Ct) values obtained for each gene were normalized to Ct values obtained for either β-actin or hypoxanthine-guanine phosphoribosyl transferase (Hprt). Data were obtained from three independent biological samples with three technical replicates.
RNA in situ hybridization
The procedure was followed as previously described (Gregorieff and Clevers, 2010) with the following modifications. In brief, cryo-sections were cut at 8 μm, thawed at room temperature for 15 min and post-fixed in 4% PFA. After incubating the tissue sections in 0.25% acetic anhydride solution, the sections were incubated in the hybridization solution with the probes, as mentioned below, for 24-48 h. The concentration of Olfm4 probe is 1 μg/ml and Gpr177 probe is 2 μg/ml. The sections were incubated with antibody overnight at 4°C and then in BM-Purple (Roche, 11 442 074 001) for up to 24 h in the dark. Slides were counterstained with Nuclear Fast Red (Vector Laboratories, H-3403) and mounted with Cytoseal mounting media (Thermo Fisher Scientific, 8310-4). Images were taken with a Nikon (TE2000-U) light microscope.
Organoid culture
Procedures for crypt isolation and organoid culture were as previously described (Das et al., 2015; Sato and Clevers, 2013; Sato et al., 2009). Organoids isolated from Gpr177L/L;Vil-Cre mice were cultured in ENR media for 7 days and counted daily. Wnt3a (100 ng/ml and 200 ng/ml) and CHIR (3 μM) were added to the epidermal growth factor/Noggin/R-spondin1 (ENR) media at day 1. Organoids isolated from Gpr177L/−;Vil-CreER mice were cultured in ENR media for 5 days. Vehicle or 4-OHT (0.5 μM) was added to the media for 12 h on day 6. The organoids were then cultured in ENR media without or with Wnt3a (100 ng/ml). The organoids were counted daily and imaged on day 6. EdU (10 μm) was added into the culture media for 30 min on day 6. The organoids were fixed in 4% PFA and the EdU staining was performed according to the manufacturer's instructions.
DSS-induced experimental colitis
The mice were administered with 2.5% dextran sulfate sodium (DSS, Colitis grade, 36-50 kDa, MP Biologics, SKU 0216011080) in tap water for 7 days and recovered with tap water for another 7 days. The body weight of individual mice was monitored daily and the graph was plotted using GraphPad Prism version 8.3. Colon tissues were harvested for pathological analysis. Colonic damage was scored blindly by a GI pathologist as described previously (Chassaing et al., 2014) in Fig. 6E,I. Briefly, scores (0-4) were assigned based on the severity of epithelial injury and leukocyte infiltration into the mucosa, submucosa and muscularis. These three scores were multiplied by an extended factor to assess the extent of the change: 1 for focal, 2 for patchy and 3 for diffuse, and summed to achieve the final score out of the maximal score of 36. In Fig. S5, colitis scoring was performed by an independent experimentalist using a different method described previously (Chinen et al., 2011).
Single-cell dissociation and sorting
The distal ileum of a wild-type mouse was harvested and briefly rinsed in ice-cold PBS (Fisher Scientific, SH30256LS). The tissue was opened longitudinally, rinsed in ice-cold PBS and further sliced to 1-2 mm pieces. The pieces were rinsed in 30 ml of ice-cold PBS by inverting the Falcon tube 10-15 times. The pieces were then transferred into fresh 30 ml of ice-cold PBS. This step was repeated until the solution remained clear after inverting. Each piece was then transferred to 30 ml crypt isolation buffer, containing 5 mM EDTA (Invitrogen, AM9260G), 2% BSA (Sigma, A3294) and HBSS Ca/Mg-free (Sigma, H9394), and shaken at 37°C for 15 min. The tube was subsequently vigorously shaken to release the epithelial layer from the pieces, after which the solution was discarded. This was repeated in another 30 ml of crypt isolation buffer. The pieces were then incubated in DMEM/F12 medium (Thermo Fisher, 12634-010) for 10 min on a Petri dish at room temperature. The solution was discarded, and the pieces were thoroughly minced. 10 ml of digestion buffer [100 U/ml Collagenase II (Thermo Fisher, 17101015), 500 U/ml DNase I (QIAGEN, 79254) and HBSS Ca/Mg free] was added to the minced tissue and shaken at 37°C for 30 min. The tissue was completely digested, and the solution was passed through a 70 µM cell strainer into a 50 ml Falcon tube. The tube was centrifuged at 200 g for 10 min. The pellet was resuspended in 1 ml of FACS solution (2% BSA in PBS). The single-cell suspensions were stained with DAPI prior to sorting. Single-cell suspensions were subjected to sorting by BD Biosciences Aria II Flow Cytometer (BD FACSAria II). Single viable lamina propria lymphocytes and mesenchymal cells were gated by forward scatter, side scatter and negative staining for DAPI.
Droplet-based scRNA-seq
Cell number and viability were determined by a propidium iodide-based fluorescence assay using Moxi GO II System, Orflo Prod#MXG102 (ORFLO Technologies). Droplet-based single-cell partitioning and single-cell RNA-seq libraries were generated using the Chromium Single-Cell 3′ Reagent v3 Kit (10X Genomics) on the 10X Chromium Controller as per the manufacturer's protocol. Live cells in single-cell suspension with 95% viability were mixed with gel beads, reverse transcription reagents and partitioning oil into a single-cell 3′ Chip and loaded onto 10X Chromium Controller for GEM generation. Briefly, the protocol includes reverse transcription, cleanup, cDNA amplification, fragmentation, end repair and A-tail prep, adapter ligation, and incorporation of sample indices into finished libraries, which are compatible with Illumina next-generation sequencing platforms. Sample quantification and quality control were determined using Qubit Fluorometer (Invitrogen, Life Technologies) and TapeStation (Agilent Technologies), respectively. cDNA libraries were sequenced on Illumina NovaSeq 6000 sequencer (Illumina) with a configuration of 28/8/0/91 bp for cell barcode, sample barcode, i5 index, and mRNA reads, respectively, as recommended by 10X Genomics. The Chromium Single Cell Software is used to analyze and visualize single-cell 3′ RNA-seq data produced by the 10X Chromium Platform. The 10X chromium software package includes Cell Ranger Pipelines and Loupe Cell Browser. Cell Ranger pipelines use raw 10X single-cell sequencing data from an Illumina sequencer and perform demultiplexing, unique barcode processing and single-cell 3′ gene counting.
For Fig. 2A-D, the wild-type count matrix files generated by Cell Ranger were transformed into the count matrix by R (version 3.6.1). Count matrices for different groups were combined and used as the input for Seurat package (version 3.1.0) (Butler et al., 2018; Stuart et al., 2019). Data normalization was performed by Seurat using the ‘LogNormalize’ method. Briefly, the gene counts for each cell were divided by the total counts for that cell and multiplied by the scale factor (10,000). The data were then natural log transformed. Cells with fewer than 200 expressed genes and genes expressed in fewer than three cells were removed from the downstream analyses. The top 2000 variable genes were used to perform PCA (primary component analysis). Ten out of 20 PCs were selected for doing clustering and Uniform Manifold Approximation and Projection (UMAP). DE (differential expression) was performed by Seurat using the Wilcoxon test. DEGs with P≤0.05, logFC>0.2 and percentage of expression in each group >0.25 were visualized by generating a heat map. Based on the expression of Wnts indicated in Fig. 2B, clusters 4, 9, 11, 14 and 21 were further selected to visualize the expression of key Wnts and stromal markers mentioned in Fig. 2C.
For Fig. S2A,B, single-cell deep constrained clustering (scDCC) was performed for Wnt4 and Wnt5a. The top 2000 most highly variable genes were identified using Seurat with the selection method of vst that calculates feature variance by fitting the variance-mean relationship. The expression data of these genes was then fed to scDCC, a deep learning model that requires two inputs: expression data and marker genes. In this study, scDCC was based on marker genes Wnt4 and Wnt5a. The restriction of binding genes was based on their distances from each other. If the distance between two genes was among the top 10% farthest distances, they were marked as not being connected, whereas if their distance was smaller than 90% of distances of other gene pairs, they were marked as must be connected. After scDCC reduced the dimensions of the scRNA-seq data, t-SNE was applied for visualization.
For Fig. 2E,F and Fig. S2C-G, the wild-type count matrix file generated by Cell Ranger was imported to Partek flow software for secondary data analysis. The single-cell analysis was generated using Partek Flow, version 9.0, build 9.0.20.0202. For each cell, quality control metrics such as number of read counts, number of detected genes and percentage of mitochondrial reads were calculated and used to filter out low-quality cells. The data were then log2+1 normalized based on counts per million and the genes that were not expressed in at least 99% of the cells were excluded from the dataset. Principal components (PCs) were calculated for the samples and the scree plot was used to determine the number of PCs to be used. Based on the PC value, cells were divided into clusters in an unsupervised fashion. The ANOVA test was used to identify differentially expressed genes (DEGs) for each cluster by comparing genes in all the clusters, filtering genes that were upregulated by 1.5-fold and then sorting them by ascending P-value. The clusters were then visualized using t-Distributed Stochastic Neighbor Embedding (t-SNE). Based on the expression of Tagln and Foxl1, the cells were selected and classified. DEGs were identified for the three populations, Foxl1−/Tagln+, Foxl1+/Tagln+ and Foxl1+/Tagln−, and the signature DEG genes (FDR set up ≤1e-8, fold change>2) were visualized by generating a heat map and by generating violin plots for genes of interest.
Transmission EM analysis and Gpr177 immunogold labeling
TEM procedures have been described previously (Gao and Kaestner, 2010; Sakamori et al., 2012). For Gpr177 immunogold labeling and EM analysis, duodenal and jejunal tissues were dissected from wild-type and Gpr177L/; Vil-Cre mice and immediately fixed as ∼1 mm fragments in 2.5% paraformaldehyde in cacodylate buffer (pH 7.4) overnight. The tissue was then sliced to 100-200 μm thickness on a vibratome and frozen between two brass ‘top-hats’ in a HPM010 (Abra Fluid) at 5000 PSI at −180°C. Next, the frozen tissue was transferred to frozen glass-distilled 100% acetone and dehydrated at −90°C for 48 h. The tissue was then infiltrated with HM-20 lowicryl and polymerized with 360 nm light at −50°C in a dry nitrogen environment. Tissue sections cut 60 nm thick containing Paneth cells were immunolabeled with Gpr177 antibody (Fu et al., 2009) at 1:250 to 1:50 in 5% BSA, 0.1% ice-cold water fish gelatin in PBS (pH 7.4). No primary antibody control and Gpr177-deficient cells were used in initial tests to optimize labeling conditions. Images shown were immunolabeled with 1:100 dilution of primary antibody. Stable antigen-antibody complexes were detected with protein A conjugated to either 15 nm gold colloids or 20 nm gold colloids. Imaging was performed with a FEI Tecnai-12 microscope at 80 keV using a nominal magnification of 6500×. Montage images were collected using Serial EM and stitched together with the IMOD subroutine Blendmont (Kremer et al., 1996; Mastronarde, 2005). Each image has a pixel dimension of ∼3 nm such that each spherical gold particle should fill five pixels and resizing of the images should provide information to reveal the approximate volume the gold would occupy in the images. Immunogold particles were counted manually, excluding particles within nuclear area. Area and perimeter of individual Paneth cells were measured using Photoshop CS. Number of immunogold particles in each subcellular compartment per unit area (for ER and Golgi) or per unit length (for plasma membrane) were calculated and compared by t-test.
Quantifications and statistical analyses
Quantification was made from at least three animals per genotype for each experiment, as previously described (Das et al., 2015; Feng et al., 2017). All comparisons were made among littermates. 50-100 crypts of small intestine and about 50 colonic crypts were counted per marker per genotype. For the number of regenerative crypts, all crypts were counted within 50 villi regardless of their sectioned shape. Length of wound was measured along multiple wound beds by ImageJ and the average length was calculated. Statistical significance was analyzed using Student's t-test, two-tailed, unpaired or one-way ANOVA. *P<0.05; **P<0.01; ***P<0.001. In Fig. 2, the statistical significance for body weights was analyzed using two-way ANOVA and the survival curve was analyzed using Gehan-Breslow-Wilcoxon test. Graphs were made using GraphPad Prism version 7.02, 8.3 and 9.
Acknowledgements
We thank Juan Flores, Sheila Bandyopadhya and Ivor Joseph for their helpful discussion of the work.
Footnotes
Author contributions
Conceptualization: S.D., Q.F., N.G.; Methodology: S.D., Q.F., I.B., X.L., H.L., O.P.-C.; Software: I.B., X.L., H.L.; Validation: S.D., Q.F., I.B., X.L., H.L., S.Y., X.Z., Y.L.; Formal analysis: S.D., Q.F., I.B., O.P.-C., S.Y., X.Z., Y.L., Z.W., L.Z., T.C.W.; Investigation: S.D., Q.F., I.B., S.Y., X.Z., Y.L., Z.W., E.M.B., M.P.V., W.H., L.Z., T.C.W.; Resources: O.P.-C., M.P.V., W.H., L.Z., T.C.W., N.G.; Data curation: S.D., Q.F., I.B., X.L., H.L.; Writing - original draft: S.D., Q.F., I.B., E.M.B., N.G.; Writing - review & editing: I.B., E.M.B., N.G.; Visualization: S.D., Q.F., I.B., X.L., H.L.; Supervision: E.M.B., N.G.; Project administration: N.G.; Funding acquisition: N.G.
Funding
This work was supported by National Institutes of Health grants (R01DK102934, R01AT010243 and R01DK119198; R01DE015654 and R01DE026936 to W.H.), an American Cancer Society Scholar Award (RSG-15-060-01-TBE), National Science Foundation/Biology/Instrument Development for Biological Research grants (1353890 and 195282), a Rutgers Initiative for Multidisciplinary Research Teams (IMRT) award to N.G. and with a New Jersey Commission on Cancer Research Postdoctoral Fellowship (DFHS16PPC045 to Q.F.). Deposited in PMC for release after 12 months.
Data availability
The single-cell RNA-seq data have been deposited in GEO under accession number GSE151152.
Peer review history
The peer review history is available online at https://journals.biologists.com/dev/article-lookup/doi/10.1242/dev.199587.
References
Competing interests
The authors declare no competing or financial interests.