Intestinal hormone-producing cells represent the largest endocrine system in the body, but remarkably little is known about enteroendocrine cell type specification in the embryo and adult. We analyzed stage- and cell type-specific deletions of Nkx2.2 and its functional domains in order to characterize its role in the development and maintenance of enteroendocrine cell lineages in the mouse duodenum and colon. Although Nkx2.2 regulates enteroendocrine cell specification in the duodenum at all stages examined, it controls the differentiation of progressively fewer enteroendocrine cell populations when deleted from Ngn3+ progenitor cells or in the adult duodenum. During embryonic development Nkx2.2 regulates all enteroendocrine cell types, except gastrin and preproglucagon. In developing Ngn3+ enteroendocrine progenitor cells, Nkx2.2 is not required for the specification of neuropeptide Y and vasoactive intestinal polypeptide, indicating that a subset of these cell populations derive from an Nkx2.2-independent lineage. In adult duodenum, Nkx2.2 becomes dispensable for cholecystokinin and secretin production. In all stages and Nkx2.2 mutant conditions, serotonin-producing enterochromaffin cells were the most severely reduced enteroendocrine lineage in the duodenum and colon. We determined that the transcription factor Lmx1a is expressed in enterochromaffin cells and functions downstream of Nkx2.2. Lmx1a-deficient mice have reduced expression of Tph1, the rate-limiting enzyme for serotonin biosynthesis. These data clarify the function of Nkx2.2 in the specification and homeostatic maintenance of enteroendocrine populations, and identify Lmx1a as a novel enterochromaffin cell marker that is also essential for the production of the serotonin biosynthetic enzyme Tph1.
The intestinal epithelium comprises five terminally differentiated cell types: the absorptive enterocytes and the secretory Paneth cells, goblet cells, tuft cells and enteroendocrine cells. Enterocytes are the major cell population in the intestine and are important for nutrient absorption. Paneth cells produce antimicrobial peptides and lysozyme, and possibly provide the stem cell niche (Porter et al., 2002; Sato et al., 2011). Goblet cells secrete mucins and thereby establish and maintain the protective mucus layer (Kim and Ho, 2010). Tuft cells comprise a rare cell population marked by doublecortin-like kinase 1 (Dclk1) expression (Gerbe et al., 2011) and are implicated in chemoreception (Gerbe et al., 2012; Sato, 2007). Enteroendocrine cells are the hormone-producing cells in the intestine. Although they represent only 1% of the cells in the intestinal epithelium, they secrete at least fifteen different types of hormones (May and Kaestner, 2010; Rindi et al., 2004) and represent the largest endocrine system in the body. Enteroendocrine cells are found in the small and large intestine and are classified by their location and principal hormone and peptide product. However, most enteroendocrine cells express more than one hormone (Arnes et al., 2012a; Egerod et al., 2012; Habib et al., 2012; Sykaras et al., 2014) and can be identified by chromogranin A (Chga) expression. Since enteroendocrine cells secrete many different hormones, they control a variety of physiological functions in the intestine and body, including gut motility, glucose homeostasis, appetite and food intake.
The serotonin [5-hydroxytryptamine (5-HT)]-producing enterochromaffin cells are the largest enteroendocrine cell population in the intestine. Approximately 90% of the 5-HT in the body is synthesized in the gut, but it is also produced in the CNS. Biosynthesis of 5-HT is a two-step process. The first step involves the conversion of the essential amino acid tryptophan to 5-hydroxytryptophan by the rate-limiting enzyme tryptophan hydroxylase (Tph). Two Tph enzymes have been found to mediate this conversion; Tph1 is expressed in the enterochromaffin cells in the intestine, whereas Tph2 is only found in the brain (Walther et al., 2003). Subsequently, 5-hydroxytryptophan becomes decarboxylated by the enzyme 5-hydroxytryptophan decarboxylase to 5-HT (Manocha and Khan, 2012).
Several transcription factors are known to regulate the enteroendocrine cell lineages. The basic helix-loop-helix (bHLH) protein neurogenin 3 (Ngn3, or Neurog3) is expressed in enteroendocrine progenitor cells and is required for induction of the enteroendocrine cell lineage (Jenny et al., 2002; Lopez-Diaz et al., 2007; Schonhoff et al., 2004). Ngn3−/− mice do not develop enteroendocrine cells in the intestinal epithelium (Jenny et al., 2002). In addition, a number of transcription factors specify subpopulations of enteroendocrine cells downstream of Ngn3, including Arx (Beucher et al., 2012; Du et al., 2012), Foxa1/2 (Ye and Kaestner, 2009), Isl1 (Terry et al., 2014), Insm1 (Gierl et al., 2006), Neurod1 (Mutoh et al., 1997; Naya et al., 1997), Pax4 (Beucher et al., 2012; Larsson et al., 1998) and Pax6 (Larsson et al., 1998). The NK2 homeobox 2 (Nkx2.2) transcription factor also regulates cell fate decisions within the enteroendocrine cell lineage in the embryo (Desai et al., 2008; Wang et al., 2009); however, postnatal lethality of Nkx2.2−/− mice (Briscoe et al., 1999; Sussel et al., 1998) precludes functional analysis of Nkx2.2 in the adult intestine. Since the intestinal epithelium undergoes constant turnover in the adult, we sought to investigate whether Nkx2.2 is required for enteroendocrine cell subtype specification in the adult as well.
In this study, we demonstrate that deletion of Nkx2.2 specifically in the intestinal epithelium in the embryo and the adult, and deletion of Nkx2.2 in Ngn3+ enteroendocrine progenitor cells, results in loss of most enteroendocrine cell types and an increase in the ghrelin (Ghrl) + cell population within the duodenum. Deletion of Nkx2.2 from the large intestine affects only a small number of enteroendocrine cell populations. Interestingly, Ghrl- and 5HT-producing cells are the most affected populations in the duodenum and colon. Overall, the intestine-specific Nkx2.2 deletion displays a developmental phenotype that is similar to that of global Nkx2.2 null mice (Desai et al., 2008; Wang et al., 2009), indicating that the misspecification of enteroendocrine cells is due to intestinal cell-intrinsic functions of Nkx2.2. Deletion of Nkx2.2 from the adult intestinal epithelium did not affect the duodenal expression of cholecystokinin (Cck), gastrin (Gast), neuropeptide Y (Npy), secretin (Sct) and vasoactive intestinal polypeptide (Vip), indicating that Nkx2.2 is dispensable for the homeostatic maintenance of these enteroendocrine lineages in the adult. Additional analysis of the intestinal epithelium in Nkx2.2 mutant mouse models carrying deletions of either the tinman (TN) domain or the NK2-specific domain (SD) revealed discrete functions of these Nkx2.2 regulatory domains in enteroendocrine cell specification. By determining gene changes that were common to the small and large intestine of all Nkx2.2 mutant mice evaluated, we identified Tph1 and the LIM homeobox transcription factor 1 alpha (Lmx1a) to be coordinately downregulated. Gene expression and deletion analyses also revealed Lmx1a to be a novel marker of 5-HT+ enterochromaffin cells, which is essential for 5-HT biosynthesis in the intestine.
Characterization of Nkx2.2Δint mice
Expression of the homeodomain transcription factor Nkx2.2 in the murine intestine begins at embryonic day (E) 15.5 and persists into adulthood (Desai et al., 2008; Wang et al., 2009). To analyze the function of Nkx2.2 in the adult intestine, we specifically deleted Nkx2.2 in the intestinal epithelium using a conditional Nkx2.2 allele (Mastracci et al., 2013) and the VillinCre/+ transgene (Madison et al., 2002). Intestine-specific deletion of Nkx2.2 circumvents the early postnatal lethality of Nkx2.2−/− mice caused by the pancreatic defect (Sussel et al., 1998). Nkx2.2flox/flox;VillinCre/+ or Nkx2.2flox/lacZ;VillinCre/+ mice are referred to hereafter as Nkx2.2Δint mice. To verify that deletion of Nkx2.2 is restricted to the intestine and does not occur in other organs, we performed PCR for the recombined Nkx2.2 allele in several representative tissues. As expected, a recombined product was only detected in intestinal tissues (Fig. S1A). Furthermore, qPCR analysis of the duodenum and colon of 6-week-old adult Nkx2.2Δint mice showed significant ablation of Nkx2.2 in the intestine (Fig. S1B).
Nkx2.2Δint mice at all ages were viable and indistinguishable from their littermate controls, with no significant change in body weight (Fig. S2A). Interestingly, we observed a small but significant increase in the length of the small but not large intestine of 6-week-old Nkx2.2Δint mice (Fig. S2B,C,J). Increases in both the villus and crypt lengths in the small intestine appeared to contribute to the overall change in length (Fig. S2D-I). Interestingly, the change in intestinal length occurred gradually and was transient: there were no length differences in neonatal animals and intestine length had normalized by 19-20 weeks (Fig. S2K,L). The transient manifestation of this phenotype at the post-weaning stage suggests that the change in length might be due to an adaptive effect.
Loss of most enteroendocrine cell populations in the duodenum of Nkx2.2Δint mice
Since the enteroendocrine lineages within the duodenum are well characterized and Nkx2.2 is expressed at the highest levels within this region of the intestine, we chose to focus on the duodenum to analyze the precise molecular changes in the intestinal epithelium of Nkx2.2Δint mice. RNA-Seq analysis of the duodenum of 6-week-old Nkx2.2Δint mice demonstrated that expression of 395 genes was significantly altered compared with controls, with a slightly larger proportion of genes upregulated than downregulated (Fig. 1A). Similar to E18.5 mice carrying a null mutation of Nkx2.2 (Nkx2.2−/−) (Desai et al., 2008; Wang et al., 2009), Nkx2.2Δint mice displayed altered enteroendocrine cell lineages. In the Nkx2.2Δint duodenum, the enteroendocrine cell marker Chga, as well as the hormones Cck, gastric inhibitory polypeptide (Gip), Npy, neurotensin (Nts), Sct, somatostatin (Sst), tachykinin 1 (Tac1), Tph1 and Vip showed significantly decreased expression. The only hormone that demonstrated a higher expression level in the duodenum of Nkx2.2Δint mice was Ghrl (Table 1). Interestingly, the expression of Gast and ppGcg is altered in the small intestine of Nkx2.2−/− mice (Desai et al., 2008), but did not appear changed in the duodenum of Nkx2.2Δint mice (Gast, fold change 0.22, P=0.38; ppGcg, fold change 1.15, P=0.68), suggesting that these changes could be secondary to loss of Nkx2.2 in the CNS or pancreas. In the intestinal epithelium, the ppGcg product is processed to glucagon-like peptide 1 (GLP-1), an incretin hormone that is important for glucose homeostasis. Consistent with unchanged ppGcg expression in the duodenum of Nkx2.2Δint mice, we could not detect a change in blood glucose levels in fed mice or after glucose challenge in an intraperitoneal glucose tolerance test (ipGTT) (Fig. S3).
We confirmed the transcriptome results by qPCR analysis of gene expression in the duodenum and immunofluorescence analysis of the small intestine. The R26RTomato reporter (Madisen et al., 2010) was used to identify recombined areas of the small intestine in Nkx2.2Δint mice (data not shown). Representative images of immunofluorescent stainings of the small intestine from Nkx2.2Δint and control mice were consistent with the changes in transcript levels. In particular, both analyses demonstrated the absence and/or decrease of Chga+ enteroendocrine cells (Fig. 1B,C), including the 5-HT+ and Sst+ subpopulations (Fig. 1D,E,H,I), whereas the Ghrl+ cell number was increased (Fig. 1F,G). With the exception of Sst, the hormonal gene expression changes observed in the duodenum of 6-week-old Nkx2.2Δint mice were also observed in the colon (Fig. 1C,E,G,I). Consistent with the change in Sst+ cells in the duodenum of 6-week-old Nkx2.2Δint mice (Fig. 1H,I), urocortin 3 (Ucn3) was also decreased (Table 2). A recent study demonstrated a reduction in the Sst+ cell number in the pancreas of Ucn3-deficient mice (van der Meulen et al., 2015), suggesting that there is a similar relationship between Ucn3+ and Sst+ cells in both the pancreas and intestine.
Although Nkx2.2 is expressed at lower levels in the colon, we sought to determine whether Nkx2.2 also regulates the enteroendocrine lineages in the large intestine. Comparative transcriptome analysis of the Nkx2.2Δint colon revealed a smaller number of gene changes but, similar to the duodenum, comprising both upregulated and downregulated genes (Fig. 1A). There was a severe reduction in expression of the pan-endocrine genes Chga and Chgb, indicating that the absence of Nkx2.2 also affects enteroendocrine lineages of the colon (Table 1, Fig. 1C). The expression of Tph1 and Ghrl was significantly altered (Table 1, Fig. 1E,G). However, the function of Nkx2.2 appears to be more limited in the colon as many of the hormones that were regulated by Nkx2.2 in the duodenum were not affected by the deletion of Nkx2.2 in the colon (Table 1, Fig. 1I). This suggests that the regulation of these lineages in the colon is independent of Nkx2.2 or, as in the case of Gip and Cck, these hormones are not highly expressed in the colon.
Since the expression of the enteroendocrine progenitor marker Ngn3 was not changed in the duodenum or colon of Nkx2.2Δint mice (Fig. 2A), it is unlikely that the change in enteroendocrine hormone expression is due to a loss of enteroendocrine progenitor cells or of the upstream progenitor populations that contribute to enteroendocrine cell lineages. For example, there is no change in atonal homolog 1 (Atoh1), which is essential for the production of all secretory cells, including enteroendocrine cells (Shroyer et al., 2007; Yang et al., 2001). Furthermore, hairy and enhancer of split 1 (Hes1), which functions upstream of Atoh1 and Ngn3 (Jensen et al., 2000; Kopinke et al., 2011), is also unchanged (Fig. 2A). Expression of even earlier genes, such as Kruppel-like factor 5 (Klf5), which marks early proliferative populations that contribute to proper cellular differentiation (Bell and Shroyer, 2015), and forkhead box A2 (Foxa2) and Klf4, which are expressed in most cells of the upper crypt and villus (Katz et al., 2002; Ye and Kaestner, 2009), is also unchanged (Fig. 2A). The lack of expression changes in these early markers of intestinal progenitor populations suggests that the defect in enteroendocrine cell specification in Nkx2.2Δint mice is likely to be specific to the enteroendocrine lineage and downstream of Ngn3+ progenitor formation.
In addition to the changes in enteroendocrine hormone expression in the duodenum, we detected differences in several transcription factors necessary for the development of specific enteroendocrine cell subtypes. Neurod1 was significantly downregulated in the duodenum of 6-week-old Nkx2.2Δint mice, consistent with the decrease in Cck and Sct expression, the two cell populations regulated by Neurod1 (Naya et al., 1997), whereas aristaless related homeobox (Arx) and Isl1 were highly upregulated (Table 2). Consistent with the increase in Ghrl+ cells, there was elevated expression of membrane bound O-acyltransferase domain containing 4 (Mboat4), the enzyme that is co-expressed with Ghrl and converts the Ghrl peptide to its biologically active, acylated form (Gutierrez et al., 2008,, 2012; Kang et al., 2012) (Table 2).
Intriguingly, there was a striking upregulation of several duodenal genes that have antimicrobial and antiviral activity, and that have been implicated in the host immune response and/or inflammation in the intestine. Among these were resistin like alpha and beta (Retnla and Retnlb) (Artis et al., 2004; Munitz et al., 2009; Wang et al., 2005), defensin beta 1 (Defb1) (Morrison et al., 2002), mast cell protease 9 (Mcpt9) (Friend et al., 1998), interferon-induced protein 44 (Ifi44) (Hallen et al., 2007), vav 1 oncogene (Vav1) (Spurrell et al., 2009), TBC1 domain family, member 23 (Tbc1d23) (De Arras et al., 2012) and several genes encoding 2′-5′ oligoadenylate synthetases (Oas1a, Oas1b, Oas1g, Oas2, Oas3) (Mashimo et al., 2003). Given that expression of Nkx2.2 is restricted to the intestinal epithelium, these gene expression changes are likely to be secondary to the dysregulation of enteroendocrine populations, such as cells expressing 5-HT, Cck and Ghrl, and support the emerging concept that enteroendocrine hormones can play immunomodulatory roles in the gut (Worthington, 2015).
It is also interesting to note that the entire gasdermin C cluster (Gsdmc, Gsdmc2, Gsdmc3, Gsdmc4) was significantly upregulated in the duodenum and colon of 6-week-old Nkx2.2Δint mice (Table 2). Although it has been demonstrated that these genes are expressed in the intestinal epithelium, the function of this subfamily is relatively uncharacterized (Tamura et al., 2007).
Deletion of Nkx2.2 in Ngn3-expressing enteroendocrine progenitor cells and in the adult intestine
To determine whether Nkx2.2 functions in enteroendocrine progenitor cells to regulate subsequent lineage decisions in the duodenum, we deleted Nkx2.2 from Ngn3+ cells using the Ngn3Cre/+ allele (Schonhoff et al., 2004). Since Nkx2.2flox/flox;Ngn3Cre/+ mice (referred to hereafter as Nkx2.2Δprogenitor) die shortly after birth with severe hyperglycemia due to the absence of insulin-producing cells in the pancreas (Sussel, et al., 1998; A. J. Churchill and L.S., unpublished), we examined the duodenum at postnatal day (P) 0 by qPCR. Although expression of Nkx2.2 was only decreased by 50% in the duodenum of Nkx2.2Δprogenitor mice, we observed reduced expression of Chga, Cck, Gast, Gip, Nts, peptide YY (Pyy), Sct, Sst, Tac1 and Tph1. The expression of ppGcg, Npy and Vip was unchanged (Fig. 2B). Similar to Nkx2.2Δint mice, Ghrl expression was significantly upregulated, although to a lesser extent (Fig. 1G, Fig. 2B). Compared with Nkx2.2Δint mice, the significantly downregulated hormones were less drastically affected in Nkx2.2Δprogenitor mice. However, this could be a consequence of lower Ngn3-Cre recombination efficiency; alternatively, it could be due to the fact that Nkx2.2 is expressed in only ∼80% of Ngn3+ cells (Wang et al., 2009).
To determine whether Nkx2.2 is also required for maintenance of the enteroendocrine cell lineages during normal cellular turnover in the adult, we deleted Nkx2.2 from the duodenum of 9- to 10-week-old tamoxifen-injected Nkx2.2flox/flox;VillinCreERT2/+ and Nkx2.2flox/lacZ;VillinCreERT2/+ mice (referred to hereafter as Nkx2.2Δadult), and compared them with 7-week-old Nkx2.2Δadult mice that were not tamoxifen injected. Analysis of the intestine from Nkx2.2floxflox;VillinCreERT2/+ mice showed that Nkx2.2 was efficiently deleted after tamoxifen injection (Fig. S4A). In addition, PCR analysis exclusively for the recombined Nkx2.2flox/+ allele, but not the Nkx2.2lacZ/+ allele, confirmed successful tamoxifen-induced deletion of exon 2 of Nkx2.2 specifically in the duodenum (Fig. S4B). Chga, Ghrl, Gip, Nts, Pyy, Sst, Tac1 and Tph1 displayed a similar expression change as in the developing duodenum of Nkx2.2Δint mice (Table 1, Fig. 2C). Ghrl was similarly upregulated, regardless of the timing of Nkx2.2 deletion from the intestine (Fig. 1G, Fig. 2C; Fig. S5A,B). Interestingly, the expression of Cck, Gast, Npy, Sct and Vip was unchanged in tamoxifen-injected Nkx2.2Δadult mice, suggesting that Nkx2.2 might not be required for the continued production of these enteroendocrine subtypes in the adult. Furthermore, ppGcg was expressed at significantly lower levels in tamoxifen-injected Nkx2.2Δadult mice (Fig. 2C). We conclude that Nkx2.2 is necessary for the maintenance of only a subset of enteroendocrine cell populations in the adult.
Changes in hormone expression after mutating the SD or TN domain of Nkx2.2
Previous studies in the ventral neural tube and pancreas have shown that the TN domain of Nkx2.2 is important for interaction with the transducin-like enhancer of split (Tle) proteins (Muhr et al., 2001; Papizan et al., 2011) to regulate gene expression. The function of the SD domain of Nkx2.2 is unknown but appears to be important for endocrine cell differentiation in the pancreas (J. Levine and L.S., unpublished). To determine whether the TN or SD domains contribute to the distinct functional activities of Nkx2.2 in regulating the various enteroendocrine lineages, we generated VillinCre/+;Nkx2.2flox/TN or Nkx2.2flox/SD mice. These mice, hereafter referred to as TNint or SDint mice, express the respective mutant allele of Nkx2.2 in the intestine and are wild type for Nkx2.2 in other tissues.
Transcriptome analysis of the duodenum of 6-week-old TNint or SDint mice revealed significant gene expression changes: 299 genes were significantly altered in expression in the TNint duodenum versus 206 genes in the duodenum of SDint mice. By comparing the datasets from Nkx2.2Δint, TNint and SDint mice, we found that only eight genes were significantly changed in the duodenum of all three mutant mouse strains (Fig. 3A). Six genes were significantly downregulated (Fig. 3B), including the enteroendocrine cell marker Chga, the enteroendocrine hormone Nts, Tph1, alpha fetoprotein (Afp) and transient receptor potential cation channel, subfamily A, member 1 (Trpa1). All of these genes have been shown to be expressed in enteroendocrine cells (Cho et al., 2014; Rindi et al., 2004; Tyner et al., 1990), suggesting that the TN and SD domains in Nkx2.2 are important for enteroendocrine cell specification, especially for the Nts+ and 5-HT+ cell subtypes (Fig. 3B; Fig. S6D-G). In addition, Ghrl+ cell numbers were also increased in both mutants (Fig. S6A-C,G). However, Ghrl expression was significantly increased in the SDint mice (Fig. 3B) but only trended up in the TNint mutant (fold change 1.78, P=0.08). Interestingly, metallothionein 1 (Mt1) is the only gene that was highly upregulated in all three mutants. Since Mt1 functions as an antioxidant, its upregulation might be due to a secondary response to altered gut hormone ratios. Comparison of Nkx2.2Δint with the TNint or SDint mice revealed that the TN domain alone is important for Pyy, Sct and Tac1 expression, whereas the SD domain alone is necessary for Npy and Sst expression (Fig. 3B).
Lmx1a is expressed downstream of Nkx2.2 and regulates 5-HT production
A notably downregulated gene in the duodenum of Nkx2.2Δint, SDint and TNint mice encoded the transcription factor Lmx1a. Since expression of Tph1, the rate-limiting enzyme for 5-HT biosynthesis, was decreased to a similar degree in mice carrying each of the three mutant Nkx2.2 alleles (Fig. 3B), we hypothesized that Lmx1a might be expressed in 5-HT+ cells in the intestine. Furthermore, Lmx1a expression was significantly downregulated in the duodenum of Nkx2.2Δprogenitor and Nkx2.2Δadult mice, corresponding to the observed decrease in Tph1 expression (Fig. 2B,C, Fig. 4A). Interestingly, Lmx1b, a paralog of Lmx1a, regulates serotonergic neuron development in the brain downstream of Nkx2.2 (Cheng et al., 2003; Ding et al., 2003), suggesting that Lmx1a could be important for the specification of 5-HT+ cells. We performed immunofluorescence staining on the duodenum of 6-week-old Lmx1aeGFP/+ mice and confirmed that Lmx1a is expressed in Chga+ enteroendocrine cells in the epithelium of both the small and large intestine, and is specifically expressed in 5-HT+ cells (Fig. 4B-E). To investigate whether Lmx1a is important for the differentiation of 5-HT+ cells, we assessed the intestinal phenotype of homozygous Lmx1aeGFP/eGFP null mice (Deng et al., 2011). Since these mice die shortly after birth, we analyzed the small intestine of newborn Lmx1aeGFP/eGFP mice. Chga and Tph1 expression was significantly reduced (Fig. 4F). It has been demonstrated that Fev (Pet1) + cells can lineage label 5-HT+ cells; however, deletion of Fev does not affect the formation of the enterochromaffin population (Wang et al., 2010). Interestingly, Fev expression, as well as that of Nkx2.2, was unchanged in the small intestine of newborn Lmx1aeGFP/eGFP mice (Fig. 4F), suggesting that Lmx1a functions independently of Fev and downstream of Nkx2.2 to regulate Tph1 and 5-HT biosynthesis in the gut.
To determine whether Nkx2.2 directly activates Lmx1a, we identified an active Lmx1a enhancer element (mm9; Chr1:169730978-169733204) in the ENCODE dataset that contained two Nkx2.2 consensus binding sites. Since an enterochromaffin cell line is not available, we tested the ability of Nkx2.2 to activate the Lmx1a enhancer in the closely related pancreatic MIN6 cell line (Ishihara et al., 1993). In this cellular context, Nkx2.2 was able to activate the Lmx1a enhancer in a luciferase assay (Fig. 4G). Future studies in a more relevant cellular context will be necessary to confirm whether endogenous Lmx1a is a direct target of Nkx2.2 in the intestinal enterochromaffin cells.
Previous studies have analyzed the function of the homeodomain transcription factor Nkx2.2 in the brain and pancreas (Briscoe et al., 1999; Sussel et al., 1998). Its role in the intestine has only been analyzed in a global knockout during embryonic development (Desai et al., 2008; Wang et al., 2009). In this study, we analyzed the intrinsic function of Nkx2.2 in the duodenum by deleting Nkx2.2 specifically in the intestinal epithelium using VillinCre/+ mice (Madison et al., 2002). Nkx2.2Δint mice are viable, but display a transiently elongated small intestine and a reduction in most enteroendocrine hormones. Despite the dramatic changes in gut hormones, glucose homeostasis remains normal. Interestingly, mice with an intestine-specific deletion of the enteroendocrine progenitor marker Ngn3, which display a loss of all enteroendocrine cells, have a different metabolic phenotype. Ngn3Δint mice are severely growth retarded, frequently die in the first week after birth, have a smaller intestine and have improved glucose clearance (Mellitzer et al., 2010). We postulate that the differences in metabolic phenotypes between these two mouse models might be attributed to the hormone Ghrl, which is highly upregulated in the Nkx2.2Δint mice (this study), but downregulated in the Ngn3Δint mutant (Mellitzer et al., 2010).
The transient increase in length of the small intestine of 6-week-old Nkx2.2Δint mice, as opposed to the decrease seen in the Ngn3Δint mice (Mellitzer et al., 2010), is also of interest. The increased length could be due to compensatory intestinal growth as an adaptive response to the alteration of gut hormone expression. For example, the altered ratios of hormones that either stimulate or inhibit proliferation, such as Gast and Sst (Thomas et al., 2003), could favor excess growth. Alternatively, we observed increased expression of insulin-like growth factor 1 (Igf1) in the Nkx2.2Δint mice at 6 weeks of age (Table 2). Igf1 is positively regulated by luminal nutrients and is able to promote growth of the epithelium of the small intestine. Furthermore, transgenic mice expressing human IGF1 exhibit a longer small intestine as well as increased villus length and crypt depth similar to Nkx2.2Δint mice (Ohneda et al., 1997). A study in rats also showed that Ghrl administration increases serum levels of IGF1, thereby stimulating duodenal growth (Warzecha et al., 2006), suggesting that the increase in Ghrl contributes to the increase in intestinal length in Nkx2.2Δint mice. It is possible that these compensatory responses are not triggered in the Ngn3Δint mice because they have a more uniform lack of all hormones.
Similar to the Nkx2.2 null mice, we demonstrated that most enteroendocrine hormones are significantly reduced in 6-week-old Nkx2.2Δint mice, with the exception of Ghrl, which is highly increased, and Gast and ppGcg, which are unchanged (Fig. S7A). However, in contrast to deletion of Nkx2.2 throughout the intestinal epithelium, Npy and Vip expression was unchanged in the Nkx2.2Δprogenitor mice, suggesting that the differentiation of these two subtypes is independent of Nkx2.2 function in the Ngn3+ progenitor cells (Fig. S7B). Deletion of Nkx2.2 in the adult intestine (Nkx2.2Δadult) also showed that Nkx2.2 plays an essential role postnatally in maintaining enteroendocrine specification during the normal turnover of enteroendocrine cells. However, during normal turnover, Nkx2.2 does not appear to be required for maintaining the Cck+, Gast+, Npy+, Sct+ and Vip+ enteroendocrine cell populations, suggesting that although these cell lineages are specified by Nkx2.2 in the embryo they are maintained in the adult by alternative mechanisms. Alternatively, ppGcg was specifically downregulated in Nkx2.2Δadult mice (Fig. S7C), suggesting there might be distinct regulatory programs for GLP-1+ cells in the adult versus the developing duodenum. Currently, it is not well understood why some enteroendocrine cell populations should be differentially regulated.
Our studies also begin to clarify the position of Nkx2.2 within the known enteroendocrine regulatory pathways. For example, the expression of Neurod1, a transcription factor that is essential for Cck+ and Sct+ cell development (Naya et al., 1997), was severely reduced in Nkx2.2Δint mice, suggesting that Nkx2.2 regulates the Cck+ and Sct+ cell lineages through the regulation of Neurod1. Furthermore, there is a correlation between Nkx2.2 and Arx regulation of Sst, in that Sst is upregulated in Arx-deficient intestine (Beucher et al., 2012; Du et al., 2012) and downregulated when Arx expression is increased in Nkx2.2Δint mice, suggesting that Sst might be regulated by Arx downstream of Nkx2.2.
In addition to elucidating the relationship between known intestinal regulatory proteins, our studies have identified Lmx1a as a novel regulator of the 5-HT signaling pathway in the gut. Enterochromaffin cells are the major source of 5-HT in the body, regulating a variety of processes, including gut motility (Gershon, 2013). Although Lmx1a RNA expression has been reported in the intestine, the function of Lmx1a was not investigated (Makarev and Gorivodsky, 2014). We have demonstrated that Lmx1a is co-expressed with Chga and 5-HT in enterochromaffin cells and is essential for the expression of Tph1, the gut-specific 5HT-synthesizing gene. Interestingly, the Lmx1a paralog Lmx1b regulates Tph2 (Song et al., 2011) and serotonergic neuron development downstream of Nkx2.2 in the brain (Cheng et al., 2003; Song et al., 2011). Since Lmx1b is only expressed at extremely low levels in the intestine (Makarev and Gorivodsky, 2014), it is likely that Lmx1a performs analogous functions downstream of Nkx2.2 for intestinal 5-HT+ cell development. Interestingly, other regulators of 5-HT+ cell development do not appear to be conserved between the intestine and CNS. For example, Fev functions downstream of Nkx2.2 in the brain to specify 5-HT neurons (Cheng et al., 2003; Hendricks et al., 2003). In the intestine, however, expression of Fev was unchanged in Nkx2.2Δint mice. Furthermore, although Fev+ cells lineage-label 5-HT+ cells, Fev-deficient mice do not show a change in Tph1 or Chga expression, or in 5-HT+ cell numbers (Wang et al., 2010). In addition, Ascl1 – another key regulator of hindbrain 5-HT+ cells (Tsarovina et al., 2004) – is not expressed in the intestine. Furthermore, in the brain Ngn3 represses the serotonergic neuron fate (Carcagno et al., 2014), whereas in the intestine Ngn3 is expressed in enteroendocrine progenitor cells and is required for all enteroendocrine lineages, including 5-HT+ enterochromaffin cells. These findings suggest that, although Nkx2.2 and Lmx1a/b may represent conserved essential components of the transcriptional pathway regulating 5-HT+ lineages, other constituents of these pathways in the intestine and CNS have diverged to provide important tissue-specific 5-HT+ cell identities (Fig. 5).
In conclusion, our data show that Nkx2.2 is required for the specification of enteroendocrine cells during development and that it is necessary for the maintenance of most enteroendocrine lineages in the adult (Fig. S7). In addition, we have identified Lmx1a as a novel marker for 5-HT+ cells that is expressed downstream of Nkx2.2 to regulate Tph1 expression, which is analogous to the role of Lmx1b in CNS 5-HT+ cells (Fig. 5). Further investigation of the shared and distinct regulatory pathways of 5-HT+ cells in the intestine and CNS will help elucidate the important regulatory mechanisms that regulate superficially similar cell types in two different tissues.
MATERIALS AND METHODS
Mice were housed and treated in accordance with the animal care protocol (AAAG3206) approved by Columbia University's Institutional Animal Care and Use Committee (IACUC). Mice were maintained on a C57BL/6J background (The Jackson Laboratory). VillinCre/+ [B6.SJL-Tg(Vil-cre)997Gum/J] (Madison et al., 2002) and R26RTomato [B6.Cg-Gt(ROSA)26-Sortm14(CAG-tdTomato)Hze/J] (Madisen et al., 2010) mice were obtained from The Jackson Laboratory. Nkx2.2flox/+, Nkx2.2lacZ/+, Nkx2.2TN/+, Ngn3Cre/+, VillinCreERT2/+ and Lmx1aeGFP/+ mice were described previously (Arnes et al., 2012b; Deng et al., 2011; el Marjou et al., 2004; Mastracci et al., 2013; Papizan et al., 2011; Schonhoff et al., 2004). Nkx2.2SD/+ mice (J. Levine and L.S., unpublished) were genotyped with the following PCR primers: 5′-GCGGCAGCACCGGCAGCCGCA-3′ and 5′-GACAACGTTAACGTTGGGATG-3′. To analyze recombination of the Nkx2.2flox/+ allele in different tissues, the following primers were used, resulting in a 464 bp PCR product: 5′-TCCTTTTAAAAATCTGCCCACGTCT-3′ and 5′-GAGGTCAACTAGGCCTCAACTTGGT-3′. Unless otherwise indicated, adult mice were analyzed at 6 weeks of age. In all experiments, Nkx2.2flox/+, Nkx2.2lacZ/+, Nkx2.2flox/flox or wild-type mice were used as controls.
To delete Nkx2.2 in adult mice, 7-week-old Nkx2.2flox/flox;VillinCreERT2/+ and Nkx2.2flox/lacZ;VillinCreERT2/+ mice were injected intraperitoneally with 100 μl tamoxifen (100 mg/ml; Sigma, T5648) for 5 consecutive days (days 1-5) and analyzed 2 weeks after the last injection (day 19; Fig. 2C). These mice are referred to as Nkx2.2Δadult mice. Seven-week-old Nkx2.2Δadult mice that were not injected with tamoxifen served as controls for the qPCR analysis and were dissected at day 1 (Fig. 2C). Tamoxifen was prepared in corn oil (Sigma, C8267).
To analyze blood glucose levels in the fed state, measurements were obtained at the same time of day while mice were kept on a regular chow diet. Intraperitoneal glucose tolerance tests (ipGTTs) were performed after a 16 h overnight fast, followed by an intraperitoneal glucose injection (2 g/kg body weight). Blood glucose was measured at 0, 15, 30, 45, 60, 90, 120 and 150 min after the glucose injection. Blood glucose measurements were taken with an Accu-Check Compact Plus glucose monitor (Model GT; Roche).
Histology and immunofluorescence
Intestines were cut longitudinally, washed with cold PBS and rolled into ‘swiss rolls' (Moolenbeek and Ruitenberg, 1981). After overnight fixation in 4% paraformaldehyde at 4°C, samples were cryopreserved with 30% sucrose and cryo-embedded in Tissue-Tek O.C.T. (Fisher Scientific, 14-373-65). Sections were cut to 5 μm thickness.
For immunofluorescence staining, sections were incubated for 15 min in 0.3% H2O2, washed in PBS and blocked for 30 min at room temperature with 10% donkey serum (Fisher Scientific, NC9624464) in PBT (PBS with 0.3% Triton X-100). Primary antibodies were diluted in 5% donkey serum in PBT and incubated on the sections overnight at 4°C. The following primary antibodies were used: rabbit anti-chromogranin A (1:500-1000; ImmunoStar, 20085), goat anti-ghrelin (1:200; Santa Cruz Biotechnology, sc-10368), rabbit anti-5-HT (1:200; ImmunoStar, 20079), rabbit anti-Sst (1:200; Phoenix Pharmaceuticals, H-060-03) and rat anti-Sst (1:500; Abcam, ab30788). The GFP signal was detected by direct fluorescence of the protein. After washing with PBT, sections were incubated with appropriate secondary antibodies diluted in 5% donkey serum in PBT for 2 h at room temperature. Secondary antibodies were conjugated with Alexa 488 or Alexa 647 (1:200; Jackson ImmunoResearch). Nuclei were stained with DAPI (1:1000; Invitrogen) for 15 min at room temperature. Sections were mounted with fluorescence mounting medium (Dako, S3023). Images were acquired with either a Zeiss LSM710 confocal microscope (Zen 2012 software) or a Leica DM5500B upright microscope (LAS AF version 126.96.36.19966 software).
To analyze tissue morphology, sections were stained with Alcian Blue (pH 2.5; Sigma, A-3157) to visualize goblet cells and counterstained with Nuclear Fast Red (Vector Laboratories, NC9483816).
Gene expression analysis
The duodenum or colon of adult mice (2 cm, measured from the stomach or caecum), 1 cm of the duodenum or the whole small intestine of newborn mice was dissected and stored at −20°C in RNAlater (Ambion, AM7021) until total RNA was extracted using the RNeasy Mini or Midi Kit (Qiagen, 74106 or 75144). cDNA was prepared with random hexamer primers and the SuperScript III First-Strand Synthesis System (Invitrogen, 18080-051). Quantitative real-time PCR (qPCR) was performed with TaqMan assays (Applied Biosystems; Table S1) and qPCR MasterMix (AnaSpec, RT-QP2X-03-15+) on a CFX96 Touch Real-Time PCR Detection System (Bio-Rad). A standard two-step real-time PCR program was used with an annealing temperature of 61°C and 40 cycles of amplification. All gene expression values were normalized to cyclophilin B (Ppib) and the fold change between wild-type and mutant samples calculated. All samples were analyzed in triplicate.
RNA sequencing was performed by the Columbia Genome Center (Columbia University) on the duodenum (n=3) and colon (n=2) of 6-week-old mice. Libraries were prepared from total RNA [RNA integrity number (RIN)>8] with the TruSeq RNA Preparation Kit (Illumina). Libraries were then sequenced using the HiSeq2000 instrument (Illumina). More than 20 million reads were mapped to the mouse genome (UCSC/mm9) using TopHat (Trapnell et al., 2009) (v2.0.4) with four mismatches and ten maximum multiple hits. Significantly differentially expressed genes were calculated using DEseq (Anders and Huber, 2010). RNA-Seq data have been deposited at GEO under series accession numbers GSE72761, GSE72762, GSE72764 and GSE78902.
A 2.226 kb fragment containing an active Lmx1a enhancer element (mm9; Chr1:169730978-169733204) was cloned into the pGL4.27 luciferase vector (Promega). 1 µg of the experimental vector pGL4.27:Lmx1a enhancer region (Lmx1a-Enh) was co-transfected with 0.1 µg Renilla luciferase vector pRL into MIN6 cells (Ishihara et al., 1993) in triplicate. The MIN6 cells were recently validated in our laboratory by RNA-Seq and tested for contamination. Luciferase activity was measured after 48 h using the Dual Luciferase Assay System (Promega). pcDNA3:myc-Nkx2.2 (pcDNA3:Nkx2.2) has been described previously (Anderson et al., 2009; Raum et al., 2006). Luciferase values were normalized to Renilla activity to account for transfection efficiencies and expressed as fold increase over the empty vector.
Results are expressed as mean±s.e.m. Statistical analysis on qPCR data and measurements was performed using a two-tailed unpaired Student's t-test. P<0.05 was considered significant. Fig. 3A was made with the help of Venn diagram generator (http://www.bioinformatics.lu/venn.php).
We thank S. Robine (Institut Curie, Paris, France) for the VillinCreERT2 mouse line; Giselle Dominguez and Shouhong Xuan for providing technical assistance; the Columbia University Genome Center for performing RNA-Seq; and members of the L.S. lab for helpful discussions.
S.G. designed and performed experiments, analyzed data and wrote the paper. D.C.G. assisted with the immunostaining and cell counting experiments. D.A.B. performed the colon RNA-Seq studies and the luciferase assays. J.M.D. provided Lmx1a mutant mice and edited the manuscript. T.L.M. assisted with generation of the Nkx2.2 floxed allele and edited the manuscript. T.P. and J.E. made the Lmx1aeGFP mice and edited the manuscript. L.S. oversaw the entire project, designed experiments, analyzed data and wrote the paper.
Funding was provided by the National Institutes of Health/National Institute of Diabetes and Digestive and Kidney Diseases [R01 DK082590] and American Diabetes Association [7-11-MN-61]. Deposited in PMC for release after 12 months.
The RNA-Seq data discussed in this publication have been deposited in NCBI Gene Expression Omnibus and are accessible through GEO series accession numbers GSE72761 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE72761), GSE72762 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE72762), GSE72764 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE72764) and GSE78902 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE78902).
The authors declare no competing or financial interests.