Role of YAP in early ectodermal specification and a Huntington's Disease model of human neurulation
Abstract
The Hippo pathway, a highly conserved signaling cascade that functions as an integrator of molecular signals and biophysical states, ultimately impinges upon the transcription coactivator Yes-associated protein 1 (YAP). Hippo-YAP signaling has been shown to play key roles both at the early embryonic stages of implantation and gastrulation, and later during neurogenesis. To explore YAP’s potential role in neurulation, we used self-organizing neuruloids grown from human embryonic stem cells on micropatterned substrates. We identified YAP activation as a key lineage determinant, first between neuronal ectoderm and nonneuronal ectoderm, and later between epidermis and neural crest, indicating that YAP activity can enhance the effect of BMP4 stimulation and therefore affect ectodermal specification at this developmental stage. Because aberrant Hippo-YAP signaling has been implicated in the pathology of Huntington’s Disease (HD), we used isogenic mutant neuruloids to explore the relationship between signaling and the disease. We found that HD neuruloids demonstrate ectopic activation of gene targets of YAP and that pharmacological reduction of YAP’s transcriptional activity can partially rescue the HD phenotype.
Editor's evaluation
This manuscript harnesses an organoid model of human neurulation to unravel the role of the Hippo signalling pathway in the specification of the three key ectodermal cell types. The authors then investigate how these mechanisms are dysregulated in an organoid model of Huntington's Disease. This work will be of broad interest to readers interested in the process of neurulation and how dysregulation of developmental mechanisms may lead to disease conditions in adulthood.
https://doi.org/10.7554/eLife.73075.sa0Introduction
Although the study of embryogenesis benefits from a broad array of model organisms and the conservation of mechanisms across species, elucidation of the particularities of human development remains challenging. The use of human pluripotent stem cells for the creation of synthetic human organoids has greatly enhanced our ability to investigate human development (Kim et al., 2020). For example, the process of neurulation can be mimicked in vitro through micropattern-based neuruloids (Haremaki et al., 2019). This technique, which allows the generation of hundreds of virtually identical organotypic cultures, offers a powerful opportunity to investigate the molecular and biophysical principles of human neurulation and the associated diseases.
The Hippo pathway is an ancient and highly conserved signaling cascade that operates in numerous cell types and a variety of organisms. In contrast to developmental pathways that require specific ligand-receptor interactions, the Hippo pathway functions as an integrator of molecular signals and biophysical states, including cell polarity, adhesion, GPCR signaling, and mechanical forces sensed through the cytoskeleton. With key roles in development, homeostasis, and regeneration, the pathway interacts with other developmental and regulatory signals including Wnt, Notch, and TGFβ. Despite the apparent complexity, the core cascade is relatively simple, including the parallel mammalian STE20-like kinases 1 and 2 (MST1/2) and mitogen-activated protein kinase kinase kinase kinase (MAP4K). The activity of these enzymes leads to phosphorylation of large tumor suppressor kinases 1 and 2 (LATS1/2), which in turn phosphorylate and limit the nuclear localization of Yes-associated protein 1 (YAP) and the homologous WW domain-containing transcription regulator 1 (TAZ). When they translocate to the nucleus, YAP and TAZ cooperate with TEA-domain transcription factors (TEAD1/2/3/4) to activate the expression of a variety of genes associated with proliferation, cell survival, de-differentiation, and cellular morphogenesis (Davis and Tapon, 2019; Hansen et al., 2015; Juan and Hong, 2016; Moya and Halder, 2019).
Hippo-YAP signaling is essential at various stages of embryonic development, including the zygote, early blastomeres, the inner cell mass (ICM), trophectoderm (TE), and the emergence of individual organs (Wu and Guan, 2021). Because Yap−/− mice die at E8.5 owing to the failure of endothelium formation in the yolk sac, and Yap−/−, Taz−/− embryos fail prior to the morula stage, investigating the potential role of Yap in neurulation in vivo requires complex transgenic strategies (Morin-Kensicki et al., 2006; Nishioka et al., 2009). The use of neuruloids derived from human embryonic stem cells (hESCs) therefore offers a convenient, highly reproducible, and quantitative alternative means of investigating the role of Hippo signaling in neurulation.
Here, we identify YAP activation as a key lineage determinant, first between neuronal ectoderm (NE) and nonneuronal ectoderm (NNE), and later between epidermis and neural crest (NC), in human embryo models. Exploiting the reproducible phenotype of Huntington’s Disease (HD) mutations in neuruloids, we additionally find that HD neuruloids demonstrate stronger YAP’s downstream transcriptional signature and that the pharmacological reduction of YAP’s transcriptional activity partially reverses the HD phenotype.
Results
Differential YAP activity in neural and nonneural ectoderm
In order to study the role of the Hippo pathway during human neurulation in vitro, we applied a neuruloid assay (Haremaki et al., 2019). Pluripotent hESCs were grown on circular, 500-µm diameter circular substrates and converted to ectodermal fates by dual SMAD inhibition for 3 days (SB431542 and LDN193189, SB+LDN, and neural induction) followed by removal of LDN and BMP4 stimulation (SB+BMP4 and neurulation). As previously reported, stimulation until the fourth day after plating, which we designate hereafter as D4, led to the induction of NNE at the edge and NE in the center (Haremaki et al., 2019; Figure 1A).
To assess Hippo-YAP signaling during these early stages of neuruloid formation, we genetically engineered the RUES2 hESC line (National Institutes of Health #0013) to endogenously and biallelically tag the C-terminus of YAP protein with green fluorescent protein (GFP; Figure 1—figure supplement 1 and Figure 1—figure supplement 2; Franklin et al., 2020). During the pluripotency state at the outset of culture, at D0, YAP was highly expressed in the nuclei, and progressively enriched in the cytoplasm upon neural induction. YAP expression was downregulated by D3 (Figure 1B and C).
Live imaging 24 hr after BMP4 stimulation (D4), revealed upregulation of YAP expression at the edge of the colony, in the region where the NNE lineage would emerge. This upregulation at the periphery was BMP4-dependent. At this stage, the upregulation of YAP at the periphery depended on BMP4 stimulation, which might directly regulate YAP expression, for SMAD1, the transcription factor downstream of BMP4 signaling, appeared to be enriched in its promoter and regulatory regions (Figure 1—figure supplement 3A and B). Near the border of the colony, YAP could be observed to undergo sporadic nuclear enrichment (Figure 1C; Figure 1—video 1; Figure 1—video 2, Figure 1—video 3) suggesting dynamic regulation of the protein’s localization in the developing NNE lineage.
YAP activation with TRULI augments differentiation of nonneural ectoderm
Immunoblot analyses of D4 neuruloids confirmed a rise in YAP levels and revealed phosphorylation of YAP at residue S127, consistent with regulation through the Hippo pathway (Figure 1D). Concomitant application of the Lats-kinase inhibitor TRULI (Kastan et al., 2021) suppressed YAP phosphorylation without altering the dynamics of the protein’s concentration during this initial phase of neuruloid self-organization (Figure 1D). The immunoblots also showed that at the pluripotent stage (D0) YAP was expressed as two protein isoforms (Figure 1D; Figure 1—figure supplement 1), both downregulated during neural induction (D1–D3; Figure 1D; Figure 1—figure supplement 4A). The upregulation upon BMP4 induction primarily involved the larger, upper band.
Consistent with the loss of S127 phosphorylation, administration of TRULI in conjunction with BMP4 increased the level of YAP protein and induced its nuclear localization in cells at the edges of neuruloid colonies (Figure 1E; Figure 1—figure supplement 5; Figure 1—video 4, Figure 1—video 5 and Figure 1—video 6). Finally, TRULI-treated neuruloids demonstrated more robust TFAP2A expression, suggesting enhanced NNE induction (Figure 1E). These results together suggest that: YAP expression is upregulated in response to BMP4 at the periphery of the colony and is subsequently regulated in part by the Hippo pathway, contributing directly to induction of the NNE lineage.
YAP activity contributes to differentiation of neural crest and epidermis
After the initial separation of the NE and NNE lineages at D4, maturation to a complete neuruloid occurs by D7 (Haremaki et al., 2019). In the course of this process, the NE organizes into a compact, closed structure at the center of the colony. The NNE differentiates further, forming two principal derivatives: NC, which is found in a radially symmetric arrangement at the periphery of the neuruloid, and epidermal cells (E), which grow atop the culture (Figure 2A). The overall structure models the early stages of human neurulation (Haremaki et al., 2019). On D6, YAP levels rose sharply with a concomitant increase of phospho-S127 YAP (Figure 2B; Figure 1—figure supplement 4B). The effect predominantly involved the larger isoform, which continued to be upregulated following BMP4 induction. Fluorescence-activated cell sorting on D7 demonstrated that enrichment of YAP in NNE (PAX6−) relative to NE (PAX6+) persisted through maturation of the neuruloid (Figure 2C).
We assessed the expression and localization of YAP within the major ectodermal lineages at different times during neuruloid formation (Figure 2D and E; Figure 2—figure supplement 1). This confirmed the exclusion of YAP from the nuclei both of NE (PAX6+) cells and of the emerging neural-crest cells (SOX10+; Figure 2D). YAP was most notably present in the nuclei of cells expressing the early epidermal marker keratin 18 (KRT18; Figure 2E).
Given that YAP activity is enriched first in undifferentiated NNE and subsequently in epidermis, we inquired whether YAP signaling stimulates their differentiation. Taking advantage of a protocol for the direct differentiation of NC (Tchieu et al., 2017), we produced SOX10+ colonies with only sparse KRT18+ cells (Figure 3—figure supplement 1A). The administration of TRULI, which enhanced YAP activity (Figure 1E and Figure 1—figure supplement 2), resulted in the complete loss of NC-like SOX10+ colonies and a significant increase in KRT18 expression (Figure 3—figure supplement 1A). Consistent with these observations, in a converse experiment TRULI treatment boosted KRT18 expression in an epidermal-differentiation protocol (Figure 3—figure supplement 1B; Tchieu et al., 2017). Furthermore, TRULI-treated neuruloids displayed increased differentiation toward epidermis and a reduction in SOX10+ cells (Figure 3). Surprisingly, TRULI treatment also precluded the closure of the central PAX6+ domain, the simulacrum of the neural tube, despite its being devoid of detectable YAP (Figure 3).
These experiments highlight the role of YAP in an in vitro model of human neurulation and suggest that the protein is dynamically regulated through canonical Hippo signaling during the specification of ectodermal lineages. Following an initial downregulation in ectodermal cells, the expression of YAP protein is induced by BMP4 stimulation and its activity supports the induction of the NNE lineage. YAP subsequently contributes directly to the differentiation of NNE by repressing the development of the NC lineage and promoting the emergence of epidermal cells. Excessive YAP activity therefore affects the balance in specification of NC and epidermis and hinders the ability of the central NE domain to fully compact.
Developing HD neuruloids display increased YAP activity
HD is caused by an expansion of CAG repeats in the Huntington locus (HTT). The development of HD can be modeled in RUES2 isogenic lines that are engineered to introduce CAG expansion, generating a unique phenotypic signature in HD-neuruloids (Haremaki et al., 2019). Using this approach, we characterized the phenotype of HD-neuruloids made from RUES2 bearing 56 CAG repeats. Interestingly, HD-neuruloids generate the same signature as TRULI-treated neuruloids: both fail to display closure of the central PAX6+ area (Figure 3). Intrigued by the similarities between the two phenotypes, and in view of the potential roles of YAP in human ectodermal differentiation, we next investigated whether the HD phenotype in neuruloids stems from hyperactivity of YAP.
To directly compare the behavior of YAP in WT and HD neuruloids, we genetically engineered an HD line to express endogenously tagged YAP-GFP (Figure 1—figure supplement 1 and Figure 1—figure supplement 2). The behavior of YAP in WT and HD neuruloids was indistinguishable during the first 3 days of neural induction, when YAP was excluded from nuclei and progressively downregulated (Figure 4—figure supplement 1). Both WT and HD lines likewise upregulated YAP expression at the edges of colonies upon BMP4 stimulation. However, at D4 YAP appeared to be more enriched in nuclei at the most peripheral region of HD-neuruloids compared to WT controls (Figure 4A; Figure 4—video 1, Figure 4—video 2, Figure 4—video 3). A similar result was obtained by immunofluorescence microscopy, which confirmed elevated YAP nuclear signals in HD colonies (Figure 4B; Figure 1—figure supplement 5). Consistent with these observations, immunoblotting analysis of D4 HD neuruloids revealed a reduction in the phosphorylated form of YAP in comparison to WT controls (Figure 4C). These observations suggest increased YAP activity in early HD neuruloids.
As observed after TRULI treatment of WT RUES2 cells (Figure 1E), D4 HD neuruloids displayed more TFAP2A at the edges of colonies (Figure 4D; Figure 4—figure supplement 2). This observation implies that increased YAP activity in early HD neuruloids has early and direct consequences for fate acquisition in the developing NNE. Consistent with elevated YAP activity, by D7, both HD neuruloids and TRULI-treated WT neuruloids displayed a remarkably expanded epidermal lineage and a reduction in the neural-crest population, as well as a failure in NE compaction (Figure 4E).
These results indicate that neuruloids obtained from hESCs carrying an HD mutation display elevated YAP activity. Consistent with our observation that YAP activity regulates NNE and epidermal formation, at D4, the HD neuruloids demonstrate more robust induction of NNE, which results by D7 in significant expansion of the epidermis at the expense of NC and in incomplete closure of the central NE domain.
Neural ectoderm of HD-neuruloids displays ectopic YAP activity
Because we observed elevated YAP activity as early as D4, we characterized this intermediate state by single-cell transcriptomic analysis. The expression profiles of lineage-marker genes in this data set allowed us to annotate the various cell populations (Figure 5—figure supplement 1A). As expected, PAX6 and TFAP2A were mutually exclusive and identified NE and NNE, respectively (Figure 5—figure supplement 1B). The NNE cell cluster could be subdivided into two subpopulations: NC progenitors expressing FOXD3 and early epidermis expressing KRT18. These populations recapitulated the three major human ectodermal lineages: NE, NC, and epidermis (Figure 5—figure supplement 1C and Figure 5—figure supplement 2).
Using the cell type-specific markers described above to identify the progenitors of the major ectodermal lineages, we employed scRNA-seq to assess their relative levels of YAP activity. YAP and key Hippo pathway components were expressed in all three lineages, most strongly in NC and epidermis (Figure 5A and B and Figure 5—figure supplement 3). However, the expression levels of YAP target genes, such as CTGF, TAGLN, and TUBB6, were much higher in epidermis than in the other ectodermal lineages (Figure 5A). To confirm these observations at a whole-genome scale, we analyzed the expression of 364 known direct YAP targets (Zanconato et al., 2015). The majority of these targets were significantly enriched in early epidermal lineage in comparison to other components of D4 neuruloid colonies (Figure 5B; Figure 5—figure supplement 4). These results confirm our initial observation in WT neuruloids that YAP activity is regulated in the ectodermal lineage and predominantly active in the epidermis lineage (Figures 2E and 3, and Figure 3—figure supplement 1).
We next used the scRNA-seq data to compare the molecular profiles of D4 early progenitors of NE, NC, and epidermis in WT and HD-neuruloids. Upon performing an analysis of differential gene expression, we ascertained that the majority of YAP target genes were significantly upregulated in HD compared to WT neuruloids (Figure 5C). This increased expression of YAP target genes was not observed in the epidermal lineage, for which YAP was highly active in both WT and HD-neuruloids. Instead, YAP targets were ectopically activated in the NC and NE lineages of HD neuruloids by comparison to the WT counterparts (Figure 5D). To confirm this observation, we performed gene set-enrichment analyses on a gene set consisting of previously identified YAP target genes (Zanconato et al., 2015). In agreement with the analysis of differential gene expression, this approach demonstrated a significant increase in the expression of YAP targets in NC and especially NE lineages of D4 HD neuruloids (Figure 5E). At higher resolution, the analysis showed that YAP targets were significantly upregulated in the NE population of D4 HD neuruloids, whereas the upregulation in NNE was not significant (Figure 5E). Within the NNE lineage, only the NC cluster displayed a significant upregulation of YAP transcriptional activity. This phenomenon was well illustrated by the behavior of the YAP target gene TUBB6, which was normally expressed specifically in the early epidermal lineage but became active throughout the HD neuruloids (Figure 5F).
These results confirm that HD neuruloids display higher YAP activity during their self-organization. The data also reveal that such HD-associated activation of YAP does not directly affect the epidermal lineage, in which YAP is endogenously active, but primarily perturbs NE and NC, in which YAP becomes ectopically active as a result of the HD mutation.
Effects of HD mutation and YAP activation on cell division rate
Because the Hippo pathway is known to regulate the cell cycle (Totaro et al., 2018; Zanconato et al., 2015), we sought to test whether the hyperactivation of YAP is accompanied by an increase in the rate of cell division in HD neuruloids. We compared the expression of phosphohistone H3 (pH3; J.-Y. Kim et al., 2017) in WT and HD neuruloids, in the absence or presence of TRULI, and determined the number of mitotic cells at different times (D4–D7; Figure 4—figure supplement 3). Enumeration of the mitotic nuclei in the developing neuruloids revealed that YAP activation by TRULI results in an increased proliferation rate, especially at the latest times (D6 and D7). A comparable but less profound effect appeared in HD neuruloids, which proliferated slightly more than WT cells. This result is consistent with a direct biological effect of YAP hyperactivation in HD neuruloids. The finding is also in agreement with reports showing that human cells carrying CAG expansions at the HTT locus display dysregulated polarity in the cell division of neural progenitors (Ruzo et al., 2018; Zhang et al., 2019) and that immortalized fibroblasts derived from HD patients show an increase in proliferation (Hung et al., 2018).
The HD phenotype in neuruloids results from ectopic YAP activation
In view of the increased YAP activity in HD-neuruloids and the associated phenotypic consequences, we inquired whether inhibition of YAP transcriptional activity could rescue the HD phenotype. To block YAP activity, we used verteporfin, which prevents the interaction between YAP and TEADs and thereby inhibits the transcription of target genes (Liu-Chittenden et al., 2012). In agreement with the observation that a minimal level of YAP activity is required for ectodermal induction (Giraldez et al., 2021), administration of verteporfin during the initial 3 days of differentiation resulted in cell death and detachment in both WT and HD samples. In conjunction with BMP4 induction, however, treatment with verteporfin from D3 was not toxic and did not affect the neuruloid phenotype of WT cells (Figure 6A–C; Figure 6—figure supplement 1). This indicates that high level of YAP activity may be dispensable for epidermal specification.
In D7 HD-neuruloids, administration of verteporfin partially rescued aspects of the HD phenotype. Although the extended epidermal population observed in HD samples was maintained, the inhibition of YAP activity partially restored the NC (SOX10+) lineage and the condensation of the central NE (PAX6+) domain (Figure 6A–C; Figure 6—figure supplement 1). This result is in agreement with the observation that in D4 HD-neuruloids, YAP target genes were significantly upregulated in the precursors of NC and NE, but not in the progenitors of the epidermal lineage (Figure 5).
Because YAP activity was elevated during the maturation of HD-neuruloids from D4, we sought to determine the period during which the inhibition of YAP transcriptional activity was most important for rescuing the closure of the central PAX6+ domain. We observed that verteporfin treatment rescued the HD-associated NE expansion if administrated together with BMP4 at D3, whereas it was unable to rescue the phenotype if added after D5 (Figure 6C).
These results indicate that the ectopic activity of YAP in HD neuruloids contributes to the malformation of their central PAX6 domains and that the effect is exerted early in neuruloid development.
Discussion
In numerous organisms and a variety of tissues, the Hippo pathway plays important roles in development, homeostasis, and regeneration (Moya and Halder, 2019). Although exploring the role of this pathway during human development remains a challenge, organoid cultures derived from hESCs offer simplified simulacra for investigation. Here, we have identified YAP as an important determinant in neurulation, during which the protein’s activity is controlled—at least in part—by the Hippo pathway. Our data suggest that YAP plays a role in the specification of early ectodermal lineages, first influencing the differentiation between NNE and NE, and then between epidermis and NC. Because YAP transcriptional activity is already robust at D4 in precursors of the epidermis, the Hippo pathway might play an early role in the specification of NNE during neurulation (Figure 6D).
In addition to the unbalanced differentiation of the NNE, YAP activation also results in the failure of the central NE domain to properly close into a compact structure (Figure 6E). Future investigations will determine whether this phenotype is a cell-intrinsic result of the ectopic YAP activation in the NE lineage, or is instead an indirect consequence of the lack of the NC lineage at its margin, which might be required to provide structural support for the closure of the NE domain. This second scenario would be consistent with a recent study showing that both NE and NNE contribute to the folding morphogenesis during human neurulation (Karzbrun et al., 2021).
In line with previous reports (Yao et al., 2014; Huang et al., 2016), our results illustrate a reciprocal relationship between Hippo-YAP and BMP4-SAMD1/5 signaling during the development of the nervous system. In particular, YAP expression rises upon BMP4 stimulation, and YAP activation in turn enhances the cellular response to BMP4 at various stages: initially by augmenting BMP4-dependent NNE induction, and later by favoring the specification of epidermis, which is known to occur at a high BMP4 concentration (Wilson and Hemmati-Brivanlou, 1995).
Earlier reports showed that YAP is involved in various stages of NC development (Zhao et al., 2021), including the proliferation of NC precursors in animal models such as Xenopus (Gee et al., 2011). Moreover, YAP has been associated with stem cell-derived NC specification in a BMP4-independent manner in vitro (Hindley et al., 2016). These reports contrast with our finding that, in the context of BMP4 induction, excessive YAP activity can skew NNE specification toward an epidermal lineage at the expense of NC. These differences might stem from an epistatic relationship between the Hippo-YAP and BMP pathways during human neurulation, when YAP activity could sensitize cells to BMP stimuli.
In agreement with recent reports showing that YAP is required for early ectodermal specification (Giraldez et al., 2021), the pharmacological inhibition of YAP transcriptional activity with verteporfin during the early phase of neural induction (D0–D3) results in cell death. However, subsequent administration of verteporfin is inconsequential for the generation of WT neuruloids (Figure 6; Figure 6—figure supplement 1). This result indicates that YAP can influence ectodermal specification during neurulation, perhaps by boosting BMP4 signaling, but its activity is not required.
Noting the similarities between TRULI-treated neuruloids and those bearing HD mutations (Haremaki et al., 2019), we explored the potential contribution of YAP dysregulation to the HD phenotype. HD-neuruloids displayed altered YAP regulation that resulted in ectopic activation of its transcriptional program, especially within the NE lineage. This HD-associated YAP hyperactivity results in an unbalanced differentiation of the NNE lineage as well as increased cell proliferation and might contribute directly to the failure of compaction in the central NE, even though this structure is devoid of YAP detectable by immunofluorescence.
Although individuals with HD have not been reported to experience a failure in folding of the neural tube or effects on fertility rate, there is a large body of evidence that connects HTT activity to developmental events. First, homozygote loss of function in the mouse leads to embryonic lethality at the onset of gastrulation, demonstrating that the function of the protein is necessary for early embryogenesis (Nasir et al., 1995). Consistent with that result, no human of the HTT−/− genotype is ever born. Second, in human embryos, the HD mutation impacts development as early as gestational week 13 (Barnat et al., 2020). Moreover, CAG expansion in HD patients causes early fetal anatomical abnormalities including an effect on ventricular volume (Hobbs et al., 2010). The earliest postnatal death due to HD in humans has been reported at 18 months of age (Nicolas et al., 2011). Perhaps the most stringent evidence linking the causality of HD to early embryogenesis is the fact that a pulse of expression of expanded CAG-HTT in a WT mouse leads to the emergence of HD symptoms in the adult (Molero et al., 2016; Wiatr et al., 2018). Third, guided differentiation of hESCs toward mature cortical neurons, self-organizing cerebroids, gastruloids, and neuruloids all point to CAG-length-dependent embryonic phenotypes (Galgoczi et al., 2021; Haremaki et al., 2019; Ruzo et al., 2018). Our study demonstrates that HTT-CAG expansion affects a model for human neurulation and that this phenotype is associated with the dysregulation of the Hippo pathway, with an increase of YAP activity early at D4 and later reduction by D7.
Owing to homeostatic and compensatory mechanisms, development in vivo is likely to be more robust than that in vitro. As a result, the neuruloid offers a sensitive assay for subtle phenomena that might otherwise remain hidden. Because our data suggest that neurulation is perturbed in HD, exploring how the developing human embryo compensates for this challenge could yield insight into the dysfunction manifested postnatally.
Several studies support the role for YAP in the pathophysiology of HD. Postmortem cortical samples from HD patients display diminished YAP activity. Moreover, YAP activators ameliorate symptoms in HD mice and normalize the phenotypes of cultured cells (Mao et al., 2016; Mueller et al., 2018). In our data, YAP activation at D4 was elevated in HD neuruloids compared to WT controls, but at D7, YAP transcription was instead diminished (Figure 5—figure supplement 6). This reversal might reflect compensation for increased YAP activity so that development could proceed normally. The reversal might alternatively represent an end state of failed development. It will be interesting to investigate whether YAP hyperactivation during early neurulation leads to hypoactivity in mature neurons. Lats1/2 knockout mice demonstrate that Hippo signaling later regulates the number and differentiation of neural progenitor (Lavado et al., 2018), and postmortem analysis of HD patients shows significantly increased proliferation of subependymal neural progenitors (Curtis et al., 2005). These results suggest that the effects of HTT mutations on this process through YAP dysregulation bear on the pathophysiology of the disease.
Our study leaves open the question of the connection between Hippo signaling and HD. Huntingtin protein is known to play a role in vesicle trafficking and establishing apical-basal polarity (Caviston and Holzbaur, 2009; Galgoczi et al., 2021; Saudou and Humbert, 2016), a prime regulator of Hippo signaling (Piccolo et al., 2014). The protein is also involved in the establishment and maintenance of intact epithelia, a function that is compromised by HD mutations (Galgoczi et al., 2021). Loss of intercellular adhesion might contribute directly to dysregulation of the Hippo pathway, reflected here as increased YAP activity. In this scenario, it is possible that CAG expansion at the HTT locus effects cell-cell contact, disrupting proper epithelial integrity and in turn dysregulating the Hippo pathway. Furthermore, lack of proper epithelial polarity could also result in the mislocalization of morphogen receptors and thus in increased sensitivity of BMP4/SMAD1/5 signaling, which in turn could affect the Hippo pathway. Future investigations of these mechanisms will reveal the molecular consequences of HD mutations and how mechanical cues interact with morphogenic signaling during embryonic development.
Materials and methods
Cell line
Request a detailed protocolThis work was based on the RUES2 (National Institutes of Health Human Embryonic Stem Cell Registry #0013) human ES cell line, the identity of which was confirmed by STR profiling (performed by Cell Line Genetics). The cells were also tested for mycoplasma and resulted negative.
Cell culture
Request a detailed protocolCells were grown in a HUESM medium that was conditioned with mouse embryonic fibroblasts and supplemented with 20 ng/ml bFGF (MEF-CM; Deglincerti et al., 2016). Cells were grown on tissue culture dishes coated with Geltrex solution (Life Technologies) and tested for mycoplasma at 2-month intervals.
Generation of YAP-GFP hESC lines
Request a detailed protocolCarboxy-terminal tagging of endogenous YAP alleles with GFP was obtained by using a donor plasmid, which was kindly donated by the Liphardt laboratory (Franklin et al., 2020). This contains an upstream homology arm approximately 1 kB in length; a cassette encoding GFP, a P2A site, puromycin resistance, and a stop codon; and a downstream homology arm also about 1 kB long. RUES2 (WT) and isogenic 56CAG (HD) hESC lines were co-transfected with Px330-Cas9-sgYAP (sgRNA sequence targeting YAP locus: TTAGAATTCAGTCTGCCTGA), donor YAP-eGFP-P2A-puromycin resistance, ePB-H2B-mCherry-BSDa, and transposase. Nucleofection was performed with program B-016 on an Amaxa Nucleofector II (Lonza). Cells were then grown for 10 days under selection by 1 μg/ml puromycin and 10 μg/ml blasticidin. Multiple clones from each line were selected for characterization by PCR genotyping with the following oligonucleotides:
GFP integration:
YAP to exon9-1F: CAGGGGTAATTACGGAAGCA
mEGFP_R: CTGAACTTGTGGCCGTTTAC.
Heterozygotic versus homozygotic:
Yap Het/Homo check LHA F: GGTGATACTATCAACCAAAGCACCC.
YAP Het/Homo check R: CATCCATCATCCAAACAGGCTCAC.
Neuruloid culture
Request a detailed protocolTo generate neuruloids (Haremaki et al., 2019), we coated micropatterned glass coverslips (CYTOO Arena A, Arena 500 A, and Arena EMB A) for 3 hr at 37°C with 10 μg/ml recombinant laminin-521 (BioLamina, LN521-05) diluted in PBS+/+ (Gibco). Single-cell suspensions were then incubated for 3 hr at 37°C in HUESM medium supplemented with 20 ng/ml bFGF and 10 μM Rock inhibitor Y27632. The micropattern culture was then washed once with PBS+/+ and incubated with HUESM with 10 μM SB431542 and 0.2 μM LDN 193189. On D3, the medium was replaced with HUESM containing 10 μM SB431542 and 3 ng/ml BMP4. On D5, the medium was replaced with the same fresh medium and then incubated until D7.
NC and epidermis differentiation
Request a detailed protocolDirect differentiation of NC was performed by a published method (Tchieu et al., 2017). RUES2 hESCs were cultured for 3 days in E6 medium (STEMCELL, 05946) containing 1 ng/ml BMP4, 10 μM SB431542, and 600 nM CHIR, then transferred for 4 days to 10 μM SB431542 and 1.5 μM CHIR.
Induction of epidermis was conducted on Transwell (Corning, CLS3413). RUES2 hESCs were cultured for 3 days in E6 medium (STEMCELL, 05946) containing 10 μM SB431542 (top of the transwell) and 10 μM SB431542 and 10 ng/ml BMP4 (bottom of the transwell). After the media had been replaced with 10 μM SB431542 and 0.2 μM LDN193189 (top and bottom), the cells were incubated for an additional 4 days.
Immunofluorescence
View detailed protocolMicropatterned coverslips were fixed with 4% paraformaldehyde (Electron Microscopy Sciences 15713) in warm medium for 30 min, rinsed three times with PBS−/−, and blocked and permeabilized for 30 min with 3% normal donkey serum (Jackson ImmunoResearch 017-000-121) and 0.5% Triton X-100 (Sigma-Aldrich 93443) in PBS−/−. Specimens were incubated with primary antibodies for 1.5 hr at room temperature, washed three times for 5 min each in PBS−/−, incubated with secondary antibodies conjugated with Alexa 488, Alexa 555, Alexa 594, or Alexa 647 (Molecular Probes) at 1/1000 dilution. After 30 min incubation with 100 ng/ml 4′,6-diamidino-2-phenylindole (DAPI; Thermo Fisher Scientific D1306), specimens were washed three times with PBS−/−. Coverslips were mounted on slides using ProLong Gold antifade mounting medium (Molecular Probes P36934).
Antibodies used:
YAP (Santa Cruz Biotechnology 101199);
TFAP2A (DSHB 3B5 concentrate);
PAX6 (BD Biosciences 561462);
SOX10 (R&D Systems AF2864); and
KRT18 (Abcam ab194130).
Imaging and analysis
Request a detailed protocolConfocal images were acquired on a Zeiss Inverted LSM 780 laser-scanning confocal microscope with a 10×, 20×, or 40× water-immersion objective. YAP-GFP Live reporter imaging was performed with a Zeiss AxioObserver z1 or spinning-disk microscope (Cell Voyager CV1000, Yokogawa).
For the acquisition of radial profiles, images were preprocessed to display the maximum-intensity projection (MIP) of at least four z-stacks. MIP images were subsequently imported into Python (czifile 2019.7.2). Individual colonies were cropped, and the radial intensity profile of fluorescence was calculated for each image starting from the center. Results were reported in arbitrary fluorescence units.
For quantitative analysis of single cells (Figure 2D), individual nuclei were segmented in the DAPI channel in each z-slice of a confocal image by a published procedure (Etoc et al., 2016). Twelve images representative of the data were selected to train Ilastik, an interactive learning and segmentation toolkit, to classify each pixel into two categories: nucleus or background. The pixels in the nucleus class were then used to create a binary mask. The original z-slice image was then processed to detect seeds based on a sphere filter with a manually set threshold. Watershedding was then applied from the defined seeds into the mask of the nuclei as defined previously by Ilastik, resulting in a set of segmented nuclei. The median intensity of PAX6, SOX10, or YAP immunofluorescence was then measured for each segmented nucleus.
Single-cell RNA sequencing
Request a detailed protocolMicropatterned glass coverslips with neuruloids 500 μm in diameter were grown in 3 ng/ml BMP4 from RUES2 hESCs and the isogenic 56CAG HD line. On D4, approximately half of the neuruloids on each coverslip were scraped off and treated for 10 min at 37°C with Accutase (Stemcell Technologies). The remaining neuruloids were fixed for immunofluorescence analysis as a quality control. After dissociation, the cells were washed three times in PBS−/− (Gibco) with 0.04% bovine serum albumin and strained through a 40-μm tip strainer (Flowmi 136800040). The number and viability of cells were determined by exclusion of trypan blue on a Countess II automated cell counter. The samples were separately loaded for capture with the Chromium System using Single Cell 3' v2 reagents (10× Genomics). Following cell capture and lysis, cDNA was synthesized and amplified according to the manufacturer’s instructions. The resulting libraries were sequenced on NovaSeq 6000 SP flowcell. The Cell Ranger software pipeline (v2.0.2) was used to create FASTQ files that aligned to the hg19 genome with default parameters. Filtered gene-expression matrices of genes versus detected barcodes (cells) with counts of unique molecular identifiers were generated and used for subsequent analyses. These data are available through NCBI accession number GSE182074.
Single-cell RNA analysis
Request a detailed protocolThe sequencing data were analyzed using Cellranger count (v3.0.2) and aggregated with the Aggr function and default settings. The aggregated datasets were processed with Seurat (v3.2.3) (Stuart et al., 2019). Quality control and filtering were performed with Scater Bioconductor (v1.16.1) to identify and remove poor-quality cells characterized by low library sizes, low numbers of expressed genes, and low or high numbers of mitochondrial reads (McCarthy et al., 2017). Following a cutoff step (counts >6000; genes >2600; 3% <mitochondrial genes <11 %), we obtained about 8000 cells in WT and roughly 6000 cells for HD. Counts of genes in each cell were normalized and log10-transformed. In order to merge the data sets using Seurat’s FindIntegrationAnchors and IntegrateData methods, we identified integration anchors for the WT and HD datasets by means of the top 2000 variable genes and first 20 principal components. The cell cycle phase was estimated by Seurat with default setting and data were rescaled to remove the effects of cell cycle and mitochondrial transcript expression. Clustering of cell was performed with Seurat’s findClusters function by use of the first 20 principal components and with two resolutions—0.5 and 0.1—to yield sets of seven clusters and two clusters, respectively. The seven-cluster set was further aggregated by manual curation of cell-type specific markers into four clusters. Uniform manifold approximation and projection (UMAP) dimension reduction implemented within Seurat was applied to the normalized, integrated data and the cluster sets projected and visualized on the resulting UMAP coordinates. Visualization of normalized expression values on UMAP coordinates and as violin plots were performed using Seurat’s FeaturePlot and VlnPlot, respectively. Differential gene expression and analysis between clusters and samples were performed by MAST (v1.8.2; Finak et al., 2015). YAP targets were retrieved from msigDB C6 gene sets and enrichment with differential expression tests assessed by single-cell GSEA implemented within MAST.
Immunoblotting
Request a detailed protocolProtein samples were resolved in 4%–15% Mini-PROTEAN TGX Stain-Free Gels (for HTT: NuPAGE 3%–8% Tris-Acetate Protein Gels), then transferred to PVDF membranes and incubated with primary antibodies overnight at 4°C followed by secondary antibodies at RT for 1 hr. Pierce ECL Western Blotting Substrate or SuperSignal West Femto Maximum Sensitivity Substrate (Thermo Fisher Scientific) was used for the detection of peroxidase activity.
Antibodies used:
Vinculin (MilliporeSigma 05-386);
GFP (Thermo Fisher Scientific A 11122)
tYAP (Santa Cruz Biotechnology sc 101199)
phYAP (CST 4911 S); and
HTT (MilliporeSigma mab2166).
Data availability
Sequencing data have been deposited in GEO under accession code GSE182074.
-
NCBI Gene Expression OmnibusID GSE182074. scRNAseq profile of Day 4 Neurloids derived from WT and HD isogenic hESC lines.
References
-
Huntington’s disease alters human neurodevelopmentScience (New York, N.Y.) 369:787–793.https://doi.org/10.1126/science.aax3338
-
Huntingtin as an essential integrator of intracellular vesicular traffickingTrends in Cell Biology 19:147–155.https://doi.org/10.1016/j.tcb.2009.01.005
-
Hippo signalling during developmentDevelopment (Cambridge, England) 146:dev167106.https://doi.org/10.1242/dev.167106
-
Concerted localization-resets precede YAP-dependent transcriptionNature Communications 11:4581.https://doi.org/10.1038/s41467-020-18368-x
-
Regulates the Self-organized Fate Patterning of hESCs-Derived GastruloidsDevelopmental Biology 1:34631.https://doi.org/10.1101/2021.03.12.434631
-
YAP and TAZ: a nexus for Hippo signaling and beyondTrends in Cell Biology 25:499–513.https://doi.org/10.1016/j.tcb.2015.05.002
-
Onset and progression of pathologic atrophy in Huntington disease: A longitudinal MR imaging studyAJNR. American Journal of Neuroradiology 31:1036–1041.https://doi.org/10.3174/ajnr.A2018
-
YAP stabilizes SMAD1 and promotes BMP2-induced neocortical astrocytic differentiationDevelopment (Cambridge, England) 143:2398–2409.https://doi.org/10.1242/dev.130658
-
A patient-derived cellular model for Huntington’s disease reveals phenotypes at clinically relevant CAG lengthsMolecular Biology of the Cell 29:2809–2820.https://doi.org/10.1091/mbc.E18-09-0590
-
Human organoids: model systems for human biology and medicineNature Reviews. Molecular Cell Biology 21:571–584.https://doi.org/10.1038/s41580-020-0259-3
-
Scater: Pre-processing, quality control, normalization and visualization of single-cell RNA-seq data in RBioinformatics (Oxford, England) 1:btw777.https://doi.org/10.1093/bioinformatics/btw777
-
Hippo-YAP/TAZ signalling in organ regeneration and regenerative medicineNature Reviews. Molecular Cell Biology 20:211–226.https://doi.org/10.1038/s41580-018-0086-y
-
Juvenile Huntington disease in an 18-month-old boy revealed by global developmental delay and reduced cerebellar volumeAmerican Journal of Medical Genetics. Part A 155A:815–818.https://doi.org/10.1002/ajmg.a.33911
-
The biology of YAP/TAZ: hippo signaling and beyondPhysiological Reviews 94:1287–1312.https://doi.org/10.1152/physrev.00005.2014
-
Chromosomal instability during neurogenesis in Huntington’s diseaseDevelopment (Cambridge, England) 145:dev156844.https://doi.org/10.1242/dev.156844
-
YAP/TAZ upstream signals and downstream responsesNature Cell Biology 20:888–899.https://doi.org/10.1038/s41556-018-0142-z
-
Huntington Disease as a Neurodevelopmental Disorder and Early Signs of the Disease in Stem CellsMolecular Neurobiology 55:3351–3371.https://doi.org/10.1007/s12035-017-0477-7
-
Hippo Signaling in Embryogenesis and DevelopmentTrends in Biochemical Sciences 46:51–63.https://doi.org/10.1016/j.tibs.2020.08.008
-
BMP2-SMAD signaling represses the proliferation of embryonic neural stem cells through YAPThe Journal of Neuroscience 34:12039–12048.https://doi.org/10.1523/JNEUROSCI.0486-14.2014
-
Genome-wide association between YAP/TAZ/TEAD and AP-1 at enhancers drives oncogenic growthNature Cell Biology 17:1218–1227.https://doi.org/10.1038/ncb3216
-
Hippo-Yap Pathway Orchestrates Neural Crest OntogenesisFrontiers in Cell and Developmental Biology 9:706623.https://doi.org/10.3389/fcell.2021.706623
Article and author information
Author details
Funding
CHDI Foundation (A-9423)
- Francesco M Piccolo
- Tomomi Haremaki
- Qingyun Tian
- Tiago L Laundos
- Andrew J Beaudoin
- Ali H Brivanlou
National Institute of General Medical Sciences (T32GM007739)
- Andrew J Beaudoin
European Molecular Biology Organization (EMBO-LTF-254-2019)
- Riccardo De Santis
National Institute on Deafness and Other Communication Disorders (R21DC016984)
- Ksenia Gnedeva
Howard Hughes Medical Institute
- AJ Hudspeth
Fundação para a Ciência e a Tecnologia (PD/BD/127997/2016)
- Tiago L Laundos
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
FMP, TH, QT, TLL, AJB, and AHB were supported by the CHDI foundation (A-9423). NRK was supported by NIGMS Grant T32GM007739. TLL was also supported by FCT PD/BD/127997/2016. RDS was supported by EMBO-LTF-254-2019. KG was supported by NIDCD Grant R21DC016984. AJH is an Investigator of Howard Hughes Medical Institute.
Copyright
© 2022, Piccolo 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,300
- views
-
- 287
- downloads
-
- 12
- 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
-
- Cancer Biology
- Cell Biology
Understanding the cell cycle at the single-cell level is crucial for cellular biology and cancer research. While current methods using fluorescent markers have improved the study of adherent cells, non-adherent cells remain challenging. In this study, we addressed this gap by combining a specialized surface to enhance cell attachment, the FUCCI(CA)2 sensor, an automated image analysis pipeline, and a custom machine learning algorithm. This approach enabled precise measurement of cell cycle phase durations in non-adherent cells. This method was validated in acute myeloid leukemia cell lines NB4 and Kasumi-1, which have unique cell cycle characteristics, and we tested the impact of cell cycle-modulating drugs on NB4 cells. Our cell cycle analysis system, which is also compatible with adherent cells, is fully automated and freely available, providing detailed insights from hundreds of cells under various conditions. This report presents a valuable tool for advancing cancer research and drug development by enabling comprehensive, automated cell cycle analysis in both adherent and non-adherent cells.
-
- Cell Biology
- Developmental Biology
In most murine species, spermatozoa exhibit a falciform apical hook at the head end. The function of the sperm hook is not yet clearly understood. In this study, we investigate the role of the sperm hook in the migration of spermatozoa through the female reproductive tract in Mus musculus (C57BL/6), using a deep tissue imaging custom-built two-photon microscope. Through live reproductive tract imaging, we found evidence indicating that the sperm hook aids in the attachment of spermatozoa to the epithelium and facilitates interactions between spermatozoa and the epithelium during migration in the uterus and oviduct. We also observed synchronized sperm beating, which resulted from the spontaneous unidirectional rearrangement of spermatozoa in the uterus. Based on live imaging of spermatozoa-epithelium interaction dynamics, we propose that the sperm hook plays a crucial role in successful migration through the female reproductive tract by providing anchor-like mechanical support and facilitating interactions between spermatozoa and the female reproductive tract in the house mouse.