Lola regulates Drosophila adult midgut homeostasis via non-canonical hippo signaling
Abstract
Tissue homeostasis and regeneration in the Drosophila midgut is regulated by a diverse array of signaling pathways including the Hippo pathway. Hippo signaling restricts intestinal stem cell (ISC) proliferation by sequestering the transcription co-factor Yorkie (Yki) in the cytoplasm, a factor required for rapid ISC proliferation under injury-induced regeneration. Nonetheless, the mechanism of Hippo-mediated midgut homeostasis and whether canonical Hippo signaling is involved in ISC basal proliferation are less characterized. Here we identify Lola as a transcription factor acting downstream of Hippo signaling to restrict ISC proliferation in a Yki-independent manner. Not only that Lola interacts with and is stabilized by the Hippo signaling core kinase Warts (Wts), Lola rescues the enhanced ISC proliferation upon Wts depletion via suppressing Dref and SkpA expressions. Our findings reveal that Lola is a non-canonical Hippo signaling component in regulating midgut homeostasis, providing insights on the mechanism of tissue maintenance and intestinal function.
Introduction
Maintenance of tissue homeostasis, as in the intestinal epithelium, is under complex regulation so to achieve a dynamic balance in terms of the rate for cell turnover (Guo et al., 2016). Uncontrolled tissue regeneration leads to intestinal malignancies such as colorectal cancer (Terzić et al., 2010), whereas a steady and homeostatic condition is preferable to maintain the integrity of the intestinal tissue for its normal function such as food digestion and nutrient absorption (Casali and Batlle, 2009; Guo et al., 2016). The Drosophila adult midgut, functionally equivalent to the mammalian small intestine, consists of a single epithelial layer where mature cell types differentiate apical-basally from the intestinal stem cells (ISCs) scattered along the basal side (Jiang et al., 2016). ISCs undergo asymmetric divisions that give rise to a renewable ISC and a non-dividing immature enteroblast (EB), which further differentiates into either an absorptive enterocyte (EC) or a secretory enteroendocrine (ee) cell (Micchelli and Perrimon, 2006; Ohlstein and Spradling, 2006). Previous studies have shown that both ISCs and EBs, commonly referred as midgut precursors, express the Snail/Slug family transcription factor escargot (Micchelli and Perrimon, 2006). Whereas ISCs are marked by the Notch (N) ligand Delta (Dl) (Ohlstein and Spradling, 2007), EBs can be labeled by a reporter of N signaling, Su(H)Gbe-lacZ (Su(H)-Z), due to the activation by neighboring ISCs (Micchelli and Perrimon, 2006; Ohlstein and Spradling, 2006). Terminally differentiated ECs, labeled by the class II POU domain transcription factor Pdm1 or Brush Border Myosin (MyoIA), acquire their large size and polyploid nuclei via endoreplication (Jiang et al., 2009; Lee et al., 2009). On the other hand, ee cells with small nuclei are specifically stained by Prospero (Pros) (Micchelli and Perrimon, 2006; Ohlstein and Spradling, 2006; Singh et al., 2012). These differential patterns in gene expression allow the identification of distinct cell types in the midgut, and provides strategic means to analyze the features and properties of adult stem cells.
Midgut homeostasis is maintained via a basal level of ISC turnover to replace the cells loss during normal gut function (Antonello et al., 2015; Jiang et al., 2011). Upon stress or injury, however, ISCs undergo rapid proliferation in order to replenish enough cells in a limited time for regeneration (Amcheslavsky et al., 2009; Buchon et al., 2010; Buchon et al., 2009b). These processes are regulated by a number of signaling pathways such as N (Ohlstein and Spradling, 2007), Wingless (Lin et al., 2008), JAK-STAT (Beebe et al., 2010; Jiang et al., 2009) and EGFR (Buchon et al., 2010; Jiang et al., 2011). Particularly, Hippo signaling has been shown to play pivotal roles in both Drosophila midgut homeostasis and regeneration via cell-autonomous and non-cell-autonomous mechanisms (Karpowicz et al., 2010; Ren et al., 2010; Shaw et al., 2010; Staley and Irvine, 2010). As an evolutionarily conserved pathway, Hippo signaling controls organ size by balancing cell proliferation and death (Yin and Zhang, 2011). The pathway consists of a core kinase cascade in which Hippo (Hpo) kinase phosphorylates and activates Warts (Wts) kinase via interaction with the scaffold protein Salvador (Sav). Subsequently, Wts interacts with Mob as tumor suppressor (Mats) to trigger phosphorylation of the transcription coactivator Yorkie (Yki), blocking its translocation to form a complex with the transcription factor Scalloped (Sd) in the nucleus, thus inhibiting downstream signal transduction (Goulev et al., 2008; Harvey et al., 2003; Huang et al., 2005; Justice et al., 1995; Oh and Irvine, 2008; Pantalacci et al., 2003; Udan et al., 2003; Wu et al., 2003; Xu et al., 1995). Despite that Hippo signaling mainly transduces via triggering Wts phosphorylation (Udan et al., 2003; Wu et al., 2003), previous studies indicate that some upstream components regulate the Hippo signaling activity by controlling Wts protein levels. The atypical cadherin Fat (Ft) (Cho et al., 2006), the atypical myosin Dachs (D) together with the LIM domain protein Zyxin (Zyx) (Rauskolb et al., 2011), and the tumor suppressor gene Scribble (Scrib) (Verghese et al., 2012) function as Hippo components via regulating Wts protein stability. During midgut homeostasis, Hippo signaling restricts ISC proliferation by sequestering Yki in the cytoplasm, thereby deactivating downstream signaling. Inactivation of Hpo or Wts leads to enhanced ISC proliferation, same as yki overexpression which activates EGFR and JAK-STAT pathways (in ECs, non-cell-autonomously) or promotes expression of target genes such as bantam microRNA (in ISCs, cell-autonomously) (Houtz et al., 2017; Huang et al., 2014; Nolo et al., 2006; Ren et al., 2010; Shaw et al., 2010; Staley and Irvine, 2010; Thompson and Cohen, 2006). In addition, the Yki-Sd complex is considered as the major mediator for injury-induced midgut regeneration, as loss of Yki in either ISCs or ECs blocks DSS- or infection-stimulated ISC proliferation, respectively (Cai et al., 2010; Karpowicz et al., 2010; Ren et al., 2010). Despite abundant evidence underlying the significance of Hippo signaling during midgut homeostasis and regeneration, neither Yki nor Sd depletion causes defects in ISC basal proliferation, raising the controversies whether the Yki-Sd complex truly executes downstream Hippo signaling during midgut homeostasis as it does in regeneration. longitudinals lacking (lola) is a transcription factor implicated in axon growth and guidance (Crowner et al., 2002; Giniger et al., 1994), ovary cell apoptosis (Bass et al., 2007), germline/neuron stem cell differentiation and maintenance (Davies et al., 2013; Silva et al., 2016; Southall et al., 2014) and wing margin development (Bass et al., 2007; Krupp et al., 2005). It also cooperates with Dl to induce the formation of metastatic tumors and regulates cell fate in Drosophila eye by antagonizing N signaling (Ferres-Marco et al., 2006; Zheng and Carthew, 2008). Despite the presence of more than 20 alternatively-spliced isoforms, the protein product of lola contains a N-terminal BTB or POZ domain and one or two C-terminal zinc finger motifs, indicating that Lola is a Cullin-3 (Cul3)-based substrate adaptor that also binds DNA (Furukawa et al., 2003; Geyer et al., 2003; Goeke et al., 2003; Horiuchi et al., 2003; Ohsako et al., 2003; Pintard et al., 2004). In the present study, we identify Lola as a transcription factor acting downstream of Hippo signaling to restrict ISC proliferation and regulate midgut homeostasis. Our results show that the Hippo signaling pathway component Wts physically interacts with Lola and regulates its stability. Whereas reduced yki expression does not rescue the enhanced ISC proliferation induced by Wts depletion, Lola genetically interacts with Wts and exhibits a suppression effect on the overgrowth. Wts-mediated Lola stability provides a means for Hippo signaling to regulate midgut homeostasis in a manner independent of the canonical Yki-Sd complex. Furthermore, Lola restricts ISC proliferation by suppressing downstream Dref and SkpA expression levels. Taken together, our findings reveal that Lola is a novel signaling effector in regulating ISC basal proliferation and midgut homeostasis via its transcriptional activity and interaction with a central Hippo signaling component.
Results
Yki is dispensable for Wts-mediated ISC proliferation in Drosophila midgut
Given the importance of Hippo signaling in regulating Drosophila adult midgut homeostasis, it is rather surprising that lacking the Hippo effector Yki leads to no defects in ISC basal proliferation (Karpowicz et al., 2010; Ren et al., 2010; Shaw et al., 2010; Staley and Irvine, 2010). To first confirm if Hippo signaling regulates midgut homeostasis via the Yki-Sd canonical signaling, escargot-Gal4 (esg-Gal4) that drives gene expression in ISCs and EBs was employed in combination with a temperature-sensitive Gal4 repressor tub-Gal80ts to restrict the time of expression (esgts) (Micchelli and Perrimon, 2006). As shown in Figure 1A–A’, ISCs and EBs were labeled with esgts-driven GFP expression (esg-GFP). Expression of wts RNAi (esgts > wts RNAi) resulted in an increase in the number of GFP+ cells in the posterior part of the midgut, indicating that reduced wts expression leads to an expansion of the precursor cell population (Figure 1B–B’ and Figure 1—figure supplement 1A). Consistent with previous results, reduced yki expression (esgts > yki RNAi) does not lead to obvious difference in the number of GFP+ cells, suggesting that Yki is not required for overall midgut homeostasis (Figure 1C–C’). Immunostaining with antibodies against phospho-Histone 3 (p-H3) that marks mitotic cells derived from ISCs indicated an obvious increase in p-H3+ cell number when downregulating wts expression in ISCs and EBs (esgts > wts RNAi, Figure 1E), suggesting that the expansion in the precursor cell population in Figure 1B is due to enhanced ISC proliferation. Unexpectedly, co-expression of yki and wts RNAi suppressed neither the growth of precursor cell population nor the enhanced proliferation induced by wts RNAi, suggesting that Wts mediates ISC proliferation independently of Yki (Figure 1D–E and Figure 1—figure supplement 1A–F).
Next, GFP+ clones of wts mutant allele wtsx1 were generated using the mosaic analysis with a repressible cell marker (MARCM) system (Lee and Luo, 2001). Interestingly, the number of Dl-labeled ISCs was increased in wtsx1 clones, indicating that ISCs are over proliferative (Figure 1F–G’ and I). Consistent to our wts RNAi results, simultaneous expression of yki RNAi in these clones did not restore the elevated Dl+ cell number (Figure 1H–I), suggesting that the reduced yki expression does not block wtsx1-induced ISC proliferation. These results are inconsistent with previous findings that Yki works downstream of Wts and regulates the output of canonical Hippo signaling, raising concerns whether different mechanisms exist between gut and other tissues, or the yki RNAi line is not efficient enough in downregulating yki expression.
To address the inconsistency, we first sought to determine if Wts and Yki work in a tissue-specific manner. Using Drosophila wing as a model, expression of wts RNAi by MS1096-Gal4 resulted in larger wing size, a phenotype rescued by yki RNAi co-expression. These results suggest that Yki functions together with Wts in controlling the wing size, a mechanism that differs from gut homeostasis (Figure 1—figure supplement 1G–K). These results are consistent with previous findings and indicate that yki RNAi is effective in downregulating yki expression. To further support our hypothesis, a separate MARCM clone approach independent of yki RNAi was employed using yki or sd mutant allele ykiB5 or sd ∆B1. Clones were induced in parallel in adult flies raised at 25 °C for 2 or 4 days after induction. Clone sizes were compared by measuring the number of GFP+ cells per clone. Under this condition, we found that clones expressing wts RNAi were bigger than the control, and the size was not or only slightly affected by the presence of ykiB5 or sd ∆B1, respectively (Figure 1J–N and Figure 1—figure supplement 1L–P). Taken together, these results indicate a non-canonical mechanism regulating midgut homeostasis by Hippo signaling. The regulatory mechanism depends on Wts function, yet is independent of the Yki-Sd transcriptional complex in Drosophila midgut.
Lola is required for ISC proliferation and midgut homeostasis
Based on our results, we propose that Wts-mediated Hippo signaling regulates midgut homeostasis independently of Yki. To identify factors that might act downstream of Wts in mediating midgut homeostasis, a genetic screen using esgts to drive RNAi expression in precursors was conducted. Interestingly, expression of two independent RNAi lines VDRC12574 (Figure 2) and NIG12052 R-1 (Figure 2—figure supplement 1), both target and abolish lola expression (Figure 2—figure supplement 1A), resulted in an increased number of GFP+ and p-H3+ cells in the posterior midgut (Figure 2A–C and Figure 2—figure supplement 1B–D). Moreover, the p-H3+ cells were dually labeled with Dl, suggesting that these cells are proliferative ISCs (Figure 2A–C). These results indicate that reducing lola expression in gut precursor cells promotes ISC proliferation. In addition, MARCM lola5D2 (a hypomorphic allele) mutant clones were bigger in size compared to the control (Figure 2D–E’’’), indicating a growth advantage due to ISC hyperproliferation when lola expression is reduced. As shown in Figure 2E–E’’’, these lola5D2 mutant clones contained Pdm1+ ECs and Pros+ ee cells, suggesting that Lola activity is not required for ISC differentiation.
In addition to the cell-autonomous role of Lola in precursor cells, we assess the possibility that Lola functions non-cell-autonomously by expressing lola RNAi in ECs. MyoIAts is a driver combining tub-Gal80ts and Myo IA-Gal4, an enhancer trap inserted in the gut specific brush border myosin IA gene (Morgan et al., 1994). Myo IAts-driven GFP expression (MyoIA-GFP) labels ECs (Figure 2F–F’’). Compared with the control groups, expressing either lola RNAi in ECs dramatically promoted ISC proliferation, as indicated by the increased Dl+ and p-H3+ cell number (Figure 2F–H and Figure 2—figure supplement 1E–G). Moreover, reducing lola expression in ECs resulted in gut hypertrophy associated with multi-layer of intestinal cells (Figure 2I–J’’), possibly a consequence of hyperproliferation. Based on the findings that ECs interact with ISCs and regulate their proliferation (Gervais and Bardin, 2017), these results suggest a non-cell-autonomous role of Lola in regulating ISC proliferation during midgut homeostasis. Importantly, reducing lola expression in ECs did not induce non-specific apoptosis in the posterior midgut (Figure 2—figure supplement 1H–I’’), indicating that lola RNAi promotes ISC proliferation via a specific signaling pathway. Altogether, Lola is required in both precursor cells and ECs for ISC proliferation and midgut homeostasis.
Loss of Lola activates EGFR, JAK-STAT, and microRNA bantam signaling
Due to functional similarities between Lola and Wts in mediating midgut homeostasis, we next sought to determine whether Lola, like Wts, functions as a component in Hippo signaling. It has been shown that inactivation of Hippo signaling activates EGFR, JAK-STAT, and microRNA bantam pathways to stimulate ISC proliferation (Karpowicz et al., 2010; Ren et al., 2010; Shaw et al., 2010; Staley and Irvine, 2010). Interestingly, Lola depletion in either precursors or ECs caused a significant increase in the levels of dpERK, the dephosphorylated active form of MAPK (Gabay et al., 1997), and the mRNA levels of EGFR ligands Spitz (Spi), Vein (Vn), and Krn, all indicating EGFR activation (Figure 3A–E). Moreover, Lola depletion in either precursors or ECs activates JAK-STAT signaling as indicated by the elevated expressions of JAK-STAT ligand Upd3 and Stat-GFP, a multimerized Stat92E reporter driving the expression of a destabilized GFP to monitor JAK-STAT activity (Bach et al., 2007) (Figure 3F–M). Consistently, mRNA levels of Upds (Upd, Upd2, and Upd3), and the JAK-STAT target Socs36E were elevated dramatically upon Lola depletion (Figure 3N). Furthermore, expression levels of the microRNA bantam, a canonical downstream target of Hippo signaling essential for cell-autonomous ISC proliferation (Huang et al., 2014), were also elevated upon Lola depletion (Figure 3O–P’’’). These results demonstrate that loss of Lola activates EGFR, JAK-STAT, and microRNA bantam signaling, suggesting that Lola regulates ISC proliferation likely via mechanisms involving Hippo signaling.
Lola regulates midgut homeostasis independently of Yki-Sd
We next investigate if Lola function depends on the core Yki-Sd complex. A driver of ubiquitous Actin-Gal4 combined with tub-Gal80ts (Actints) was used to manipulate lola expression ubiquitously in a temporal manner. Real-time PCR analysis of total RNAs from Actints > lola RNAi guts revealed no significant difference in yki mRNA levels between control and RNAi samples, suggesting that Lola does not affect yki transcription (Figure 4A). No Yki protein accumulation was detected in lola5D2 mutant gut clones, suggesting that Lola does not affect Yki protein levels either (Figure 4B–C’). Genetic interaction analysis indicated that co-expression of either yki or sd RNAi with lola RNAi in ISCs did not suppress the increase in p-H3+ cell number induced by lola RNAi alone (Figure 4D and Figure 4—figure supplement 1A–F’). Furthermore, MARCM clones of double mutant lola5D2 and ykiB5 exhibited no significant difference in clone size compared to lola5D2 (Figure 4E–H’). MARCM sd∆B1 mutant clones expressing lola RNAi were similarly larger in size, suggesting that Sd depletion does not suppress lola RNAi-induced ISC proliferation (Figure 4I–L’). Taken together, these results indicate that Lola restricts ISC proliferation independently of the Yki-Sd transcriptional complex.
Lola restricts ISC proliferation by suppressing Dref and SkpA expression levels
To further elucidate the mechanism of Lola regulating ISC proliferation and midgut homeostasis, we turn to its function as a DNA binding protein and transcriptional suppressor (Southall et al., 2014). Chromatin immunoprecipitation assay followed by high-throughput sequencing (ChIP-seq) was performed to identify genes suppressed directly by Lola. Thousands of Lola-associated chromatin binding sites were identified in cultured S2 cells. Analysis of the Lola binding profiles revealed that Lola mainly binds to the regions around transcription start sites (TSS) and promoters. Among the top 500 defined scored peaks (Supplementary file-Table 1), genes related to stem cell mitosis and differentiation were manually selected for further analysis. Results from ChIP assay and real-time PCR showed dramatic enrichment of Lola in promoter regions of Oli, ec, Dref, and 3’UTR regions of pdm3 and SkpA. No significant binding (only a basal level of binding) was detected between Lola and the non-binding region (genebody, gb) (Figure 5A). As a negative control, Lola knockdown by dsRNA significantly reduced the enrichment (Figure 5A). To validate these results in vivo, total RNAs from Actints > lola RNAi guts were collected and analyzed. Similarly, mRNA levels of pdm3, ec, Dref, and SkpA were upregulated when reducing lola expression, suggesting that Lola suppresses the expression of these genes (Figure 5B). To further confirm that Lola represses the transcription of these genes, Lola binding regions (1: long and 2: short) of Dref or SkpA, two target genes among all others consistently affected by Lola, were cloned into upstream or downstream of the luciferase gene and the reporter activity was assayed. As shown in Figure 5C, activities of these reporters were dramatically inhibited by Lola but not Lola∆ZF12, a truncated form of Lola that does not bind to DNA. Activities of 3xSd luc, a dual reporter reflecting Sd-Yki transcriptional activity (Zhang et al., 2008), was not affected by Lola, demonstrating the specificity of Lola repression. In addition, Dref protein levels were significantly elevated in lola5D2 mutant clones (Figure 4D–F). In summary, these results suggest that Lola regulates the expression of Dref and SkpA, and represses their transcription.
To verify whether these genes are responsible for the ISC hyperproliferation induced by lola RNAi, ISC proliferation was monitored in midguts co-expressing RNAi targeting pdm3, ec, Dref or SkpA and lola. As expected, knock down of either Dref or SkpA significantly rescued the increase in GFP+ and p-H3+ cell number induced by lola RNAi (Figure 5G and Figure 5—figure supplement 1A–H’), indicating that Lola restricts ISC proliferation via regulating Dref or SkpA expression levels. On the other hand, knock down of either Pdm3 or ec did not rescue the increased proliferative ISC number induced by lola RNAi (data not shown), indicating that these target genes might be involved in other parts of Lola function. Based on these results, we conclude that Lola functions as a transcription factor to suppress the expression of mitosis-related genes Dref and SkpA, hence controlling ISC proliferation and midgut homeostasis.
Wts interacts with lola and affects lola turnover in vitro
Given that both Lola and Wts regulate ISC proliferation, and Lola deactivates similar downstream pathways as Hippo signaling, biochemical approaches were first taken to study the relationship between Hippo signaling components and Lola in cultured S2 cells. Hippo signaling components such as Hpo, Wts, Yki, Sb, or SdBP engineered with different N-terminal epitope tags for identification (Fg-Hpo, Myc-Wts, Fg-Yki, HA-Sd, or Myc-SbBP) were co-expressed with Lola tagged with a N-terminal Flag epitope (Fg-Lola). Interestingly, Fg-Lola protein levels were significantly increased only in the presence of Myc-Wts (Figure 6A). In addition, Fg-Lola was increasingly stabilized when adding increasing amount of Myc-Wts, indicating that Wts regulates Lola protein levels in a dosage-dependent manner (Figure 6B). Furthermore, Fg-Lola protein half-life in the presence or absence of Myc-Wts was monitored using the nascent protein synthesis inhibitor cycloheximide (CHX). As shown in Figure 6C–D, Fg-Lola proteins exhibited a rapid turnover rate with a half-life of approximately 2 hr, while Myc-Wts prolonged Fg-Lola protein half-life to about 6 hr. To investigate how Lola protein stability is regulated, S2 cells expressing Fg-Lola were treated with specific UPS inhibitors (MG132 or MG101) or lysosomal inhibitors (E64 or Leupeptin) for 3 hr followed by CHX treatment for additional 6 hr. As shown in Figure 6—figure supplement 1A, UPS inhibitors, but not lysosome inhibitors, protected Fg-Lola from degradation, indicating that UPS plays a major role in regulating Lola protein stability. In addition, Lola ubiquitination levels were reduced in the presence of Wts as shown by results from the in vivo ubiquitination assay (Figure 6—figure supplement 1B). These results demonstrate that Wts regulates Lola protein stability by protecting Lola from UPS-mediated degradation, hence the possibility of Hippo signaling regulating Lola function.
To further explore possible interaction between Wts and Lola, co-immunoprecipitation (Co-IP) assays in S2 cells expressing Myc-Wts and Fg-Lola were conducted. As shown in Figure 6E, Lola and Wts reciprocally interacted, validating a physical interaction between two proteins. Co-IP analysis of S2 cells expressing Myc-Wts using the anti-Lola antibody for IP pulled down Myc-Wts, further demonstrating that Wts interacts with endogenous Lola (Figure 6F). To map which Lola region binds to Wts, two truncated Lola protein variants, each engineered with a Flag tag, were constructed: the N-terminal Lola (amino acid 1 to 370) containing the BTB domain (Fg-LolaN) and the C-terminal Lola (amino acid 371 to 748) containing the two zinc finger motifs (ZF1 and ZF2, Fg-LolaC) (Figure 6—figure supplement 2A). Co-IP results indicated that Fg-LolaC, but not Fg-LolaN, co-immunoprecipitated with Myc-Wts (Figure 6—figure supplement 2B), suggesting that Wts interacts with Lola C-terminus. Deletion of either ZF1 or ZF2 or both in Lola C-terminus did not affect binding between Wts and Lola, suggesting that Wts-Lola interaction does not require DNA binding (Figure 6—figure supplement 2A and C). Furthermore, ChIP-qPCR analysis revealed no detectable change of Lola enrichment on Dref or SkpA binding regions when co-expressing Wts in S2 cells (Figure 6G), indicating that Wts does not affect Lola-DNA binding. In toto, these results suggest that Wts interacts with Lola C-terminus independently of Lola binding to DNA.
Given that Wts is a Serine/Threonine protein kinase, we next sought to determine if Lola is a Wts substrate and phosphorylated by Wts. Surprisingly, Lola exhibited no shift in molecular weight when analyzed by the phos-tag gels (Figure 6E), suggesting that Wts interacts with Lola in a manner independent of phosphorylation. To further validate this conclusion, a construct expressing Wts KD (WtsK743R, a kinase dead form of Wts) (Lai et al., 2005) was generated. Not only that Wts KD interacted with Lola by Co-IP (Figure 6—figure supplement 2C), Wts KD stabilized Lola similarly as the wild-type Wts (data not shown). These results suggest that Wts kinase activity is not essential for binding to Lola and does not affect Lola protein stability. On the other hand, previous studies have suggested that Wts interacts and phosphorylates Yki (Goulev et al., 2008; Huang et al., 2005; Oh and Irvine, 2008). In our hands, both total amount and the S168 phosphorylation level of Yki co-immunoprecipitated with Wts remained unaffected in the presence or absence of Fg-Lola, suggesting that Wts-Lola interaction does not affect Wts-Yki interaction and Wts-mediated Yki phosphorylation (Figure 6—figure supplement 2D). Taken together, these results demonstrate that Wts binds to Lola directly to moderate its turnover rate, a process independent of Wts kinase activity, Lola DNA binding, Wts-Yki interaction, and Wts-mediated Yki phosphorylation.
It is noteworthy to mention that both Wts and Lola proteins have been predicted to contain the nuclear export signal (NES) (http://www.cbs.dtu.dk/services/NetNES/) and the nuclear localization signal (NLS) (http://nls-mapper.iab.keio.ac.jp/cgi-bin/NLS_Mapper_form.cgi), indicating that Wts possibly localizes to the nucleus in addition to the cytoplasm, whereas Lola localizes to the cytoplasm in addition to the nucleus. Our results showed that expression of Wts fused with a C-terminal nuclear localization signal (NLS) (WtsNLS) caused a remarkable increase in Lola protein levels, suggesting that WtsNLS stabilizes Lola in the nucleus (Figure 6—figure supplement 1C–D), whereas the expression of Wts fused with a C-terminal nuclear export signal (NES) (WtsNES) only reduced the increased Lola protein levels partially (Figure 6—figure supplement 1C–D), suggesting the possibility of Lola translocating to cytoplasm where it is protected by Wts. These results indicate that Wts interacts with Lola both in the cytoplasm and nucleus.
Wts regulates lola protein levels in vivo
In addition to our evidence that Wts interacts with and stabilizes Lola in S2 cells, endogenous Lola protein levels were examined using antibodies against Lola. Immunostainings showed that Lola is ubiquitously expressed in all cell types in the midgut and absent in the lola5D2 MARCM clones (Figure 6—figure supplement 3A–B’’’). Consistent with the cell culture data, a dramatic decrease in the anti-Lola staining was also detected in both flip-out clones expressing wts RNAi using the driver hsflp; act >CD2>Gal4; UAS-Dicer2, UAS-GFP (AG4; Dcr2GFP > wts RNAi) and wtsx1 MARCM clones in midguts (Figure 6H–J and Figure 6—figure supplement 3C–D’’’). Furthermore, endogenous Lola protein levels were analyzed by Western blot using gut extracts that ubiquitously express lola, lola RNAi, wts, or wts RNAi. Similar to lola overexpression, wts overexpression increased Lola protein levels (Figure 6K). Collectively, these results indicate that Wts regulates Lola protein levels in vivo.
Lola is essential for Wts-mediated ISC proliferation and midgut homeostasis
Based on the results that Wts modulates Lola protein stability and Lola is required for midgut homeostasis, genetic interaction between Wts and Lola in regulating midgut homeostasis was analyzed. Whereas reduced yki expression did not rescue the increase in the GFP+ and p-H3+ cell number induced by esgts > wts RNAi (Figure 1A–E), lola overexpression dramatically reduced this increase (Figure 7A and Figure 7—figure supplement 1A–D’). Similar reduction was detected in midguts expressing lola and wts RNAi in ECs (MyoIAts > lola; wts RNAi, Figure 7—figure supplement 1E). These results indicate that lola overexpression blocks the hyperproliferation induced by wts RNAi both cell-autonomously and non-cell-autonomously. To further validate our observation, expressing lola in GFP+wtsx1 MARCM clones that exhibited the overgrowth phenotype dramatically reduced the clone size (Figure 7B–F). On the other hand, wts overexpression in lola5D2 MARCM clones did not affect the increased clone size (Figure 7—figure supplement 2A–D). Taken together, these results suggest that Lola is essential for Wts function in midgut homeostasis.
Furthermore, Dref and SkpA mRNA levels were also elevated in wtsx1 clonal midguts (Figure 7—figure supplement 1F). Genetic interaction analysis between wts and Dref or SkpA showed that co-expression of SkpA or Dref RNAi in wts RNAi midguts dramatically reduced the increase in the GFP+ and p-H3+ cell number (Figure 7A, Figure 7—figure supplement 1E and G–H’), indicating that SkpA and Dref are both required for Wts function in midgut homeostasis. These results support our hypothesis that Wts stabilizes Lola such that SkpA and Dref expressions are transcriptionally suppressed, thereby restricting ISC proliferation.
Finally, RNA sequencing (RNA-seq) analysis using RNAs isolated from adult guts expressing wts RNAi, lola RNAi or yki using Actints identified a collection of 2116 differentially expressed genes (DEGs), using the R package DESeq2 (padj <0.05) (Figure 7G, Supplementary file 2) (Love et al., 2014). Based on the DEG heatmap and hierarchical cluster analysis, a high degree of similarity was detected between samples from wts and lola RNAi, suggesting strong functional correlation between Wts and Lola, but not Yki and Lola, in regulating midgut homeostasis. On the other hand, lola overexpression did not rescue the elevated GFP+ and p-H3+ cell number induced by yki overexpression (Figure 7—figure supplement 3A–E), reinforcing the notion that Lola restricts ISC proliferation in an Yki-independent manner. Taken together, these results suggest that while Wts regulates Lola and Yki in different means, Lola are functionally similar with Wts. Yki and Lola regulate the expression of different target genes, and direct different transcription programs in the process of Hippo-mediated ISC proliferation and midgut homeostasis (Figure 7H).
Discussion
Mechanisms that regulate Drosophila midgut homeostasis and regeneration are complex and under a series of delicate controls. The Hippo pathway, which output is mainly transduced by the effector Yki-Sd complex, has been shown to be crucial for both midgut homeostasis and regeneration. In the present study, we identify a novel transcription factor Lola acting downstream of Hippo signaling to restrict ISC proliferation. The Hippo component Wts interacts with and stabilizes Lola in a Yki-independent manner. Lola then suppresses downstream Dref and SkpA expression levels to regulate ISC proliferation. Our results suggest that Lola is an effector of a non-canonical Hippo signaling pathway essential for ISC proliferation and midgut homeostasis (Figure 7H).
Lola is a new player in Hippo-mediated ISC proliferation and midgut homeostasis
Our present findings suggest that Lola regulates ISC proliferation and midgut homeostasis. Reduced lola expression in both ISCs and ECs promotes ISC proliferation, indicating that Lola functions both cell-autonomously and non-cell-autonomously. These precursors remain correctly differentiated when Lola is absent, suggesting that Lola only affects ISC proliferation, but not differentiation. Similarly, inactivation of Hippo signaling also causes enhanced ISC proliferation, indicating the possibility that Lola functions downstream of Hippo. Consistent to the characterized phenotype for Hippo signaling inactivation, downregulating lola expression in ECs activates EGFR and JAK-STAT signaling, and Lola regulates Hippo downstream bantam expression levels in ISCs. Taken together, these results suggest that Lola regulates ISC proliferation and homeostasis likely via mechanisms involving Hippo signaling.
Our biochemical results strongly support that Lola is a Hippo component. Not only that Lola interacts with and is stabilized by Wts, our in vivo analysis of ISC proliferation also indicates that lola overexpression rescues the enhanced ISC proliferation induced upon Wts depletion, a phenotype not affected by Yki. These findings indicate that Lola functions downstream of Hippo signaling via interaction with Wts. In addition, four lines of evidence suggest that Lola functions independently of the Yki-Sd complex: first, yki mRNA and protein levels are not affected by Lola, suggesting that yki is not a Lola transcriptional target. Genetic interaction analysis also shows that neither Yki nor Sd suppresses the enhanced ISC proliferation induced by lola RNAi. Next, RNA-seq transcriptional profiles of samples expressing lola RNAi strongly correlate with the ones expressing wts RNAi, but not the ones expressing yki, indicating that Yki and Lola work downstream of Wts to mediate different outputs in terms of regulating ISC proliferation. Furthermore, Lola-Wts interaction does not affect Wts-Yki binding and Wts-mediated Yki phosphorylation, suggesting Wts regulates Yki and Lola via separate means. Finally, co-expression of lola and yki does not rescue the enhanced ISC proliferation induced by yki overexpression, indicating that Lola and Yki work independently from each other. These results demonstrate that Lola is a novel effector of Hippo signaling in regulating ISC proliferation and midgut homeostasis; Lola functions by interacting with Wts and in a manner independent of the Yki-Sd complex.
Mechanisms of Lola-mediated ISC proliferation and midgut homeostasis
As a transcriptional suppressor, it is intriguing to speculate that Lola regulates downstream gene expression by suppressing their transcription levels, hence affecting ISC proliferation. Using ChIP-seq analysis and luciferase reporter assays, two distinct gene targets of Lola were identified: Dref and SkpA. Lola enrichment was detected on the Dref promoter and the 3’UTR region of SkpA, respectively, indicating that the repression modes of Lola on Dref and SkpA are different. Transcription of these two genes are also directly repressed by lola. Downregulation of Dref or SkpA expression in gut precursors expressing either wts or lola RNAi rescues the enhanced ISC proliferation, providing the genetic evidence that Lola regulates ISC proliferation vis suppressing Dref and SkpA expression. Wts-mediated ISC proliferation also requires Lola and its downstream transcription targets Dref and SkpA. Moreover, co-expression of Dref RNAi and yki does not rescue the enhanced ISC proliferation induced by yki overexpression, reconfirming the Yki-independent mechanism of Lola repressing Dref/SkpA signaling in regulating ISC proliferation and midgut homeostasis (data not shown).
Wts regulates lola stability
Due to its important role in ISC proliferation and midgut homeostasis, Lola protein levels are expected to be under precise and delicate controls. Interestingly, Wts physically interacts with Lola and Lola protein levels are regulated by Wts both in vitro and in vivo. It has been proposed that Wts and Mats colocalize in the cytoplasm or at the cell cortex, thereby restricting Yki in the cytoplasm (Dong et al., 2007; Oh and Irvine, 2008; Sun et al., 2015). Exceptions have been reported, however, that activated LATS1/2 (Wts in mammals) is accumulated in the nucleus, where it phosphorylates YAP (Yki in mammals) (Li et al., 2014). Based on these observations, it is feasible to speculate that Wts translocates to the nucleus and interacts with Lola, therefore regulating Lola protein stability. Our cell culture data further indicate a possibility of Lola translocation into the cytoplasm, and Wts interacts with Lola both in cytoplasm and nucleus. Despite the possibility of Wts translocation into the nucleus (or Lola translocation into the cytoplasm), Wts is unlikely to regulate Lola stability via phosphorylation, an event happened for Yki in the cytoplasm based on our results that Wts kinase activity is not required for Wts-Lola interaction. Furthermore, Wts-Lola interaction is independent of Wts-Yki interaction and Wts-mediated Yki phosphorylation, demonstrating that Wts controls Yki and Lola activity via separate and distinct mechanisms.
Our biochemical results provide further insights on possible mechanisms of how Wts controls Lola stability. Wts interacts with the C-terminal Lola, where the DNA binding zinc finger motifs reside. Deletion of zinc finger motifs does not affect Wts-Lola binding, nor Wts affects Lola enrichment on Dref or SkpA binding regions. These results suggest that Wts-Lola interaction and Lola binding to DNA are two independent events. Wts possibly interacts with Lola at the C-terminus in regions outside of the zinc finger motifs, structuring conformational changes which indirectly affect interaction between proteins and other parts of Lola such as the N-terminal BTB domain. It has been shown that BTB domain-containing proteins serve as both a linker and substrate adaptor within the Cul3-based E3 ligases (Furukawa et al., 2003; Geyer et al., 2003; Pintard et al., 2004). Proteins containing either the Leucine-rich repeats (LRR) or WD40 domain, such as the F-box protein in the SCF E3 complex, have been shown to mediate the degradation of BTB domain substrate adaptor (Wimuttisuk et al., 2014). Interestingly, our results indicate that Lola is degraded via a UPS-dependent mechanism and its ubiquitination is affected by the presence of Wts. Based on these findings, it is possible that Wts-Lola interaction modulates the interaction between Lola and other LRR or WD40-containing proteins, hence preventing Lola ubiquitination and degradation by the UPS system.
Hippo-Lola signaling in normal midgut homeostasis
During stress or injury-induced regeneration, ISCs rapidly proliferate in order to replenish the cell loss in a short time (Amcheslavsky et al., 2009; Buchon et al., 2009a; Buchon et al., 2010; Buchon et al., 2009b). Abundant evidence has suggested that Yki is a central Hippo effector in the regeneration process. Not only that Yki protein levels and transcription levels of Hippo downstream target genes are upregulated during regeneration, loss of Yki in either precursors or ECs blocks DSS- or infection-induced ISC proliferation, respectively. yki overexpression also leads to activation of EGFR and JAK-STAT signaling pathways (Karpowicz et al., 2010; Ren et al., 2010; Shaw et al., 2010; Staley and Irvine, 2010). These results suggest that Yki and Hippo signaling are important for midgut regeneration.
Unlike regeneration, normal midgut homeostasis only requires a basal level of cell turnover and ISCs proliferate to a minimum extent to maintain the need (Antonello et al., 2015; Jiang et al., 2011). During normal condition, Hippo signaling restricts ISC proliferation, thus serving an inhibitory role in regulating ISC proliferation and midgut homeostasis. Inactivation of Wts or Hpo, similarly as yki overexpression, leads to inactivation of Hippo signaling, thus enhancing ISC proliferation (Karpowicz et al., 2010; Ren et al., 2010; Shaw et al., 2010; Staley and Irvine, 2010). Nonetheless, Yki inactivation does not affect ISC proliferation, raising the argument that Yki is not as prominently needed in midgut homeostasis as in regeneration. Conserved and similarly in mammals, depletion of Yki homologs YAP1 and/or TAZ protein in the intestine does not affect normal homeostasis, indicating that YAP1 (Barry et al., 2013; Cai et al., 2010; Camargo et al., 2007) and TAZ (Azzolin et al., 2014) might be dispensable under normal conditions. Our results that Lola mediates ISC proliferation via non-canonical Hippo signaling resolve the argument by supporting a complementary yet equally important role for Lola during Hippo-mediated ISC proliferation and midgut homeostasis. Upstream components might trigger Hippo signaling via controlling Wts protein levels instead of Wts phosphorylation, hence regulating Lola stability. Whereas Yki is critically needed during injury-induced regeneration, Lola plays a more fundamental and maintenance role during homeostasis, a process that Yki is dispensable. Collectively, our work uncovers a novel mechanism that Hippo regulates ISC proliferation via Lola-mediated non-canonical downstream signaling; Hippo-Lola signaling controls ISC basal proliferation and midgut homeostasis, whereas Yki in the canonical Hippo signaling controls rapid ISC proliferation during regeneration.
In the present study, we uncover a novel transcription factor Lola that functions downstream of Hippo signaling in regulating Drosophila midgut homeostasis. During normal homeostasis and maintenance, the Wts-Lola-Dref/SkpA signaling axis serves as a critical mediator restricting ISC proliferation, adding another layer of complexity in the regulatory mechanism of ISC proliferation. Considering the importance of Hippo signaling in intestinal diseases and tumorigenesis, our findings provide new insights on developing potential biomarkers or strategies of therapeutic targets for anti-cancer research.
Materials and methods
Plasmids and cloning
Request a detailed protocolThe full length lola DNA fragment (lola-RD: 2247 bp) was amplified from Drosophila cDNA (BDGP DGC clone LD28033) by PCR. LolaN and C are truncated lola variants with 1–1110 bp and 1111–2247 bp, respectively. Lola ∆ZF1, ∆ZF2, and ∆ZF12 are truncated lola variants deleting 1435–1498 bp, 1525–1596 bp, and 1435–1596 bp, respectively. For expression in S2 cells and flies, the DNA fragments mentioned above were cloned in frame with the Flag-tag in the pUAST-Fg vector according to the standard protocols. The NLS and NES sequences were listed as follow:
NLS: 5’-CCTAAGAAGAAGAGGAAGGTT-3’
NES: 5’-CTTCAGCTACCACCGCTTGAGAGACTTACTCTT-3’
Drosophila Stocks and Genetics
Request a detailed protocolThe fly stocks used in this study are listed in the Key Resources Table.
Temperature-controlled expression
Request a detailed protocolThe experiment using esg-Gal4/UAS-GFP; tubGal80ts and MyoIA-Gal4/UAS-GFP; tubGal80ts were crossed and cultured at 18°C to restrict Gal4 activity. The F1 adult flies were shifted to 29°C for 5 days to induce transgene expression. 10 female adult flies at least were dissected for each genotype, followed by immunostaining, microscopy, and statistical analysis.
MARCM clonal analysis
Request a detailed protocolGFP positive mutant clones were generated using the MARCM system. Flies were crossed and raised at 25°C, F1 adult flies with designated genotypes were subjected to heat shock at 37°C for 1 hr and then cultured at 25°C or 18°C for days specified before dissection. Clones of more than 10 midguts were analyzed in each group.
Immunohistochemistry
Request a detailed protocolAll immunostaining experiments were done on midguts of female flies. The guts were dissected and fixed in 4% formaldehyde (Sigma) for 1 hr and washed three times in PBS supplemented with 0.1% Triton X-100 (PBS-T). The guts were incubated in the primary antibody diluted in PBS-T overnight at 4°C, followed by three washes with PBS-T. After incubation with a second antibody at room temperature for 1 hr, midguts were washed for 3 times and mounted in PBS/glycerol medium with DAPI. Primary antibodies used in this study are listed in the Key Resources Table.
Microscopy and statistical analysis
Request a detailed protocolFluorescent microscopy was performed on a Leica LAS SP8 confocal microscope; confocal images were obtained using the Leica AF Lite system. Confocal images from the basal layer of the posterior midguts (sub-region R5a, https://flygut.epfl.ch/overview) where ISCs can be clearly visualized were taken under 40 × objective. Single layer image is shown. The quantification of p-H3+ cell per guts was undertaken by counting the p-H3+ cell numbers over the whole gut of indicated genotypes. The p-H3+ cell number quantification data were statistically present as average with the standard error of mean (SEM) and p-values of significance is calculated by Student’s T-test (tails = 2, Two-sample unequal variance): * is p<0.05, ** is p<0.01, *** is p<0.001, ns is no significance with p>0.05.
To quantify the expression of indicated protein in Figure 5 and Figure 6, the intensity of the indicated proteins and GFP signal in the view region was analyzed using the Leica AF Lite system. For each group, 50 (GFP−) of which were quantified as wild type and 50 (GFP+) were quantified as lola5D2 or wtsx1. The quantification data were statistically present as average with the standard error of mean (SEM) and p-values of significance is calculated by Student’s T-test (tails = 2, Two-sample unequal variance): * is p<0.05, ** is p<0.01, *** is p<0.001, ns is no significance with p>0.05.
Real-time PCR
Request a detailed protocolTotal RNAs were extracted from 10 midguts of 5-day-old female flies or 20 wing discs from third-instar larvae with indicated genotypes using Trizol Reagent (Invitrogen), and the cDNAs were synthesized using ReverTra Ace synthesis kit (Toyobo). Real-time PCR was performed using the SYBR Green Real-time PCR Master Mix (Toyobo) reagent with the ABI7500 System. Results were repeated for three independent biological replicates. RpL32 was used as a normalized control. The RT-PCR data were statistically present as average with the standard error of mean (SEM) and p-values of significance were calculated by Student’s T-test (tails = 2, Two-sample unequal variance) in Excel: * is p<0.05, ** is p<0.01, *** is p<0.001, ns is no significance with p>0.05.
Cell line
Request a detailed protocolS2 Drosophila cells from ATCC Ref: CRL-1963 were used. The identity has been athenticated (STR profiling). No mycoplasm contamination.
Cell culture, Transfection, Immunoprecipitation, Ubiquitination assay, Western blot and immunostaining
Request a detailed protocolS2 cells were cultured in Drosophila Schneider’s Medium (Invitrogen) with 10% fetal bovine serum, 100 U/ml of penicillin, and 100 mg/ml of Streptomycin. Plasmid transfection was carried out using X-tremeGENE HP DNA Transfection Reagent (Sigma) according to manufacturer’s instructions. A ubiquitin-Gal4 construct and the indicated constructs (cloned into pUAST expression vectors) were co-transfected for all transfection experiments. CHX (50 mg/ml, Sigma) was used to inhibit nascent protein synthesis. MG132 (50 mg/ml, Sigma) and MG-101 (50 mg/ml, Sigma) were used to inhibit the UPS activity. E64 (50 mM; Sigma) and Leupeptin (50 mM; Sigma) were used to inhibit lysosome function. Immunoprecipitation, Ubiquitination assay, Western blot and immunostaining were performed according to standard protocols as previously described (Guo et al., 2013; Zhang et al., 2008; Zhang et al., 2006).
RNA interference in S2 Drosophila cells
Request a detailed protocolFor RNAi in S2 cells, primers were designed as follows: Renilla-dsRNA-F (use as control) (5’- taatacgactcaatagggatgacgtvaaaagtttac −3’); Renilla-dsRNA-R (use as control) (5’- taatacgactcaatagggagactacatccggtttacc −3’); lola-dsRNA-F (5’- taatacgactcaatagggatggatgacgatcagcagtt −3’); lola-dsRNA-R (5’- taatacgactcaatagggcggctgccggtccgctggac −3’).
PCR reactions were used as templates for in vitro RNA production (in vitro Transcription T7 kit, Takara), and RNAs were purified using isopropyl-alcohol. dsRNAs were then annealed (68°C for 10 min and 37°C for 30 min). To perform knockdown experiment, S2 cells were diluted into 1 × 106 cells ml−1 with serum free medium for 1 hr starvation with 15 μg dsRNA.
ChIP
Request a detailed protocolChIP assays were performed using S2 cells. Cells were cross-linked for 15 min at 37°C in 1 ml of 1% formaldehyde in PBS buffer. The cross-linking was stopped by adding Glycine to a final concentration of 0.125M. Next, fixed cells were washed for 10 min in 1 ml ice-cold PBS for three times. Cells were sonicated in 1 ml sonication buffer (50 mM Hepes-KOH, pH 7.5, 140 mM NaCl, 1 mM EDTA, pH 8.0, 1% Triton X-100, 1% sodium deoxycholate, 0.1% SDS, and proteinase cocktail) with a Bioruptor sonicator. The sonication yielded genomic DNA fragments with an average size of about 200 bp. After centrifugation, lysates were incubated with 4 μg antibody for 4 hr (or overnight); 40 μl protein A/G PLUS agarose (Santa Cruz) was then added and incubated for another 4 hr (or overnight) on a rotator at 4°C. Beads were washed three times with the ChIP wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, pH 8.0, 150 mM NaCl, and 20 mM Tris-Cl, pH 8.0), and then washed again with ChIP final wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, pH 8.0, 500 mM NaCl, and 20 mM Tris-Cl, pH 8.0). Genomic DNA was eluted with elution buffer (1% SDS and 100 mM NaHCO3) at 65°C for 30 min. 5 M NaCl was added to a final concentration of 200 mM for further incubation at 65°C for 4 hr (or overnight). Then, 0.5 M EDTA and 20 mg/ml proteinase K were added until their final concentrations were 5 mM and 0.25 mg/ml, respectively. The resulted mixture was incubated at 55°C for 2 hr to digest the protein. Genomic DNA was purified with a DNA purification kit (QIAGEN) and subjected to high throughput sequencing or real-time PCR.
The ChIP-seq data were deposited in the Gene Expression Omnibus (accession number GSE136999). The list of top 500 peaks of Lola ChIP-seq profile is shown in Supplementary file 1 in detail.
Luciferase reporter assay
Request a detailed protocolTo generate the Dref-luc-1 or Dref-luc-2 reporter genes, 2 DNA segments covering Lola-binding region of Dref were subcloned between BglII and KpnI sites of the pGL3-Promoter vector (Dref-luc-1: −304 to +43, Dref-luc-2: −211 to −3), respectively. For SkpA-luc-1 or SkpA-luc-2, 2 DNA segments covering Lola-binding region of SkpA were subcloned downstream of the luciferaser gene into the BamHI site of the pGL3-Promoter vector (SkpA-luc-1: +1384 to +1716, SkpA-luc-2: +1458 to +1604), respectively.
For luciferase reporter assays, S2 cells were transfected with indicated reporters and copia-renilla luciferase reporter constructs in 24 well plate together with Fg-Lola or Fg-Lola∆ZF12 expressing constructs. Cells were incubated for 48 hr after transfection. The reporter assay was performed using the Dual-Luciferase reporter assay system (Promega). Dual-Luciferase measurements were performed in triplicate using GloMax-Multi JR Single-Tube Multimode Reader. The DLR data were statistically present as average with the standard error of mean (SEM) and p-values of significance is calculated by Student’s T-test (tails = 2, Two-sample unequal variance) in Excel: * is p<0.05, ** is p<0.01, *** is p<0.001, ns is no significance with p>0.05.
RNA-seq and statistical analysis
Request a detailed protocolThe transcriptomes were generated by RNA sequencing (RNA-seq) analysis using RNAs isolated from adult guts expressing wts RNAi, lola RNAi or yki using Actints. 30 midguts dissected from indicated genotypes incubated for 5 days at 29°C were dissected to extract total RNAs per sample for RNA-Seq experiment.
After assessing RNA quality with Agilent Bioanalyzer, mRNAs were enriched by poly-A pull-down from total RNA samples (3 ug mRNA per sample, RIN >7 required). Then, sequencing libraries constructed with Illumina TruSeq Stranded mRNA Sample Prep Kit were sequenced using Illumina HiSeq X ten. We sequenced about 26M paired-end 150 bp reads for each sample.
After quality check with fastqc (version 0.11.8), we mapped the reads to an index of Drosophila melanogaster reference genome (BDGP6) using Hisat2 (version 2.1.0) (Kim et al., 2015). Then we converted sam files to bam files with samtools (version 1.9; Li et al., 2009). The aligned reads were assigned to genes using annotations from Ensembl (Drosophila_melanogaster.BDGP6.94.gtf) and HTseq-count (version 0.11.2) (Anders et al., 2015). The differentially expressed genes (DEGs) were identified using R package DESeq2 (version 1.25.12) (Love et al., 2014). Compared with the control group, genes with padj (adjusted p-value, corrected p-value after Multiple Comparisons by method 'BH' which DESeq2 provided)<0.05 were defined as DEGs. The union of all the DEGs from 3 groups of samples emerged into a cellection of 2116 DEGs.
The heatmap of the DEGs and hierarchical cluster analysis were employed to examine the correlation of gene expression pattern from indicated samples. We converted gene counts to Z-score to center and scaled the numeric matrix by R generic function ‘scale’. The R package ComplexHeatmap (version 1.20.0) was used to draw the heatmap of the DEGs’ expression from indicated samples (Gu et al., 2016). We used euclidean method to calculate distance between these two vectors in Hierarchical cluster analysis, and ward.D2 is used as the agglomeration method.
The raw RNA-seq data was deposited in the Sequence Read Archive (accession number SRP220236). The expression of 2116 DEGs indicated by reads counts of the RNA-seq profiles are shown in Supplementary file 2.
Data availability
Sequencing data have been deposited in GEO under accession code GSE136999, and SRA under accession code SRP220236. All data generated or analysed during this study are included in the manuscript.
-
Sequence Read ArchiveID SRP220236. Wts/Lola/Yki-induced intestinal stem cell (ISC) overproliferation affects gene expression in fly midgut.
-
NCBI Gene Expression OmnibusID GSE136999. Genome-wide binding of Lola in S2 cells.
References
-
GFP reporters detect the activation of the Drosophila JAK/STAT pathway in vivoGene Expression Patterns 7:323–331.https://doi.org/10.1016/j.modgep.2006.08.003
-
The axon guidance gene lola is required for programmed cell death in the Drosophila ovaryDevelopmental Biology 304:771–785.https://doi.org/10.1016/j.ydbio.2007.01.029
-
YAP1 increases organ size and expands undifferentiated progenitor cellsCurrent Biology 17:2054–2060.https://doi.org/10.1016/j.cub.2007.10.039
-
Intestinal stem cells in mammals and DrosophilaCell Stem Cell 4:124–127.https://doi.org/10.1016/j.stem.2009.01.009
-
Delineation of a fat tumor suppressor pathwayNature Genetics 38:1142–1150.https://doi.org/10.1038/ng1887
-
Lola regulates midline crossing of CNS axons in DrosophilaDevelopment 129:1317–1325.
-
Targeting of protein ubiquitination by BTB-Cullin 3-Roc1 ubiquitin ligasesNature Cell Biology 5:1001–1007.https://doi.org/10.1038/ncb1056
-
Tissue homeostasis and aging: new insight from the fly intestineCurrent Opinion in Cell Biology 48:97–105.https://doi.org/10.1016/j.ceb.2017.06.005
-
Lola encodes a putative transcription factor required for axon growth and guidance in DrosophilaDevelopment 120:1385–1398.
-
Maintenance of the adult Drosophila intestine: all roads lead to homeostasisCurrent Opinion in Genetics & Development 40:81–86.https://doi.org/10.1016/j.gde.2016.06.009
-
Intestinal stem cell response to injury: lessons from DrosophilaCellular and Molecular Life Sciences 73:3337–3349.https://doi.org/10.1007/s00018-016-2235-9
-
HISAT: a fast spliced aligner with low memory requirementsNature Methods 12:357–360.https://doi.org/10.1038/nmeth.3317
-
Mosaic analysis with a repressible cell marker (MARCM) for Drosophila neural developmentTrends in Neurosciences 24:251–254.https://doi.org/10.1016/S0166-2236(00)01791-4
-
The sequence alignment/Map format and SAMtoolsBioinformatics 25:2078–2079.https://doi.org/10.1093/bioinformatics/btp352
-
The molecular cloning and characterization of Drosophila Melanogaster myosin-IA and myosin-IBJournal of Molecular Biology 239:347–356.https://doi.org/10.1006/jmbi.1994.1376
-
The bantam microRNA is a target of the hippo tumor-suppressor pathwayCurrent Biology 16:1895–1904.https://doi.org/10.1016/j.cub.2006.08.057
-
In vivo regulation of yorkie phosphorylation and localizationDevelopment 135:1081–1088.https://doi.org/10.1242/dev.015255
-
The salvador partner hippo promotes apoptosis and cell-cycle exit in DrosophilaNature Cell Biology 5:921–927.https://doi.org/10.1038/ncb1051
-
Cullin-based ubiquitin ligases: cul3-btb complexes join the familyThe EMBO Journal 23:1681–1687.https://doi.org/10.1038/sj.emboj.7600186
-
Zyxin links fat signaling to the hippo pathwayPLOS Biology 9:e1000624.https://doi.org/10.1371/journal.pbio.1000624
-
Generation and staining of intestinal stem cell lineage in adult midgutMethods in Molecular Biology 879:47–69.https://doi.org/10.1007/978-1-61779-815-3_4
-
Dedifferentiation of neurons precedes tumor formation in lola mutantsDevelopmental Cell 28:685–696.https://doi.org/10.1016/j.devcel.2014.01.030
-
Localization of hippo signalling complexes and warts activation in vivoNature Communications 6:8402.https://doi.org/10.1038/ncomms9402
-
Inflammation and colon cancerGastroenterology 138:2101–2114.https://doi.org/10.1053/j.gastro.2010.01.058
-
Hippo promotes proliferation arrest and apoptosis in the salvador/Warts pathwayNature Cell Biology 5:914–920.https://doi.org/10.1038/ncb1050
-
Identifying tumor suppressors in genetic mosaics: the Drosophila lats gene encodes a putative protein kinaseDevelopment 121:1053–1063.
-
Hippo signaling: a hub of growth control, tumor suppression and pluripotency maintenanceJournal of Genetics and Genomics 38:471–481.https://doi.org/10.1016/j.jgg.2011.09.009
-
Lola regulates cell fate by antagonizing notch induction in the Drosophila eyeMechanisms of Development 125:18–29.https://doi.org/10.1016/j.mod.2007.10.007
Article and author information
Author details
Funding
Chinese Academy of Sciences (Strategic Priority Research Program XDB19000000)
- Lei Zhang
National Natural Science Foundation of China (No. 31530043)
- Lei Zhang
National Natural Science Foundation of China (No. 31625017)
- Lei Zhang
National Natural Science Foundation of China (No. 31871039)
- Margaret S Ho
National Natural Science Foundation of China (No. 31900551)
- Jinhui Li
National Natural Science Foundation of China (2017YFA0103601)
- Lei Zhang
National Natural Science Foundation of China (2019YFA0802001)
- Lei Zhang
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Jin Jiang, Rongwen Xi, Xiaohang Yang, Masamitsu Yamaguchi, Vienna Stock Center, Bloomington Stock Center, NIG, Developmental Studies Hybridoma Bank for supplying antibodies and fly stocks. We would like to thank Gang Wei for the generous help on high-throughput sequencing data analysis. This work is supported by National Key Research and Development Program of China (2017YFA0103601, 2019YFA0802001), National Natural Science Foundation of China (No. 31530043, No. 31625017, No. 31871039, and No. 31900551), ‘Strategic Priority Research Program’ of the Chinese Academy of Sciences (XDB19000000), Shanghai Leading Talents Program.
Copyright
© 2020, Hao et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,056
- views
-
- 418
- downloads
-
- 4
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Cells called alveolar myofibroblasts, which have a central role in the development of the lung after birth, receive an orchestrated input from a range of different signaling pathways.
-
- Developmental Biology
Premature infants with bronchopulmonary dysplasia (BPD) have impaired alveolar gas exchange due to alveolar simplification and dysmorphic pulmonary vasculature. Advances in clinical care have improved survival for infants with BPD, but the overall incidence of BPD remains unchanged because we lack specific therapies to prevent this disease. Recent work has suggested a role for increased transforming growth factor-beta (TGFβ) signaling and myofibroblast populations in BPD pathogenesis, but the functional significance of each remains unclear. Here, we utilize multiple murine models of alveolar simplification and comparative single-cell RNA sequencing to identify shared mechanisms that could contribute to BPD pathogenesis. Single-cell RNA sequencing reveals a profound loss of myofibroblasts in two models of BPD and identifies gene expression signatures of increased TGFβ signaling, cell cycle arrest, and impaired proliferation in myofibroblasts. Using pharmacologic and genetic approaches, we find no evidence that increased TGFβ signaling in the lung mesenchyme contributes to alveolar simplification. In contrast, this is likely a failed compensatory response, since none of our approaches to inhibit TGFβ signaling protect mice from alveolar simplification due to hyperoxia while several make simplification worse. In contrast, we find that impaired myofibroblast proliferation is a central feature in several murine models of BPD, and we show that inhibiting myofibroblast proliferation is sufficient to cause pathologic alveolar simplification. Our results underscore the importance of impaired myofibroblast proliferation as a central feature of alveolar simplification and suggest that efforts to reverse this process could have therapeutic value in BPD.