Selective activation of FZD7 promotes mesendodermal differentiation of human pluripotent stem cells
Abstract
WNT proteins are secreted symmetry breaking signals that interact with cell surface receptors of the FZD family to regulate a multitude of developmental processes. Studying selectivity between WNTs and FZDs has been hampered by the paucity of purified WNT proteins and by their apparent non-selective interactions with the FZD receptors. Here, we describe an engineered protein, called F7L6, comprised of antibody-derived single-chain variable fragments, that selectively binds to human FZD7 and the co-receptor LRP6. F7L6 potently activates WNT/β-catenin signaling in a manner similar to Wnt3a. In contrast to Wnt3a, F7L6 engages only FZD7 and none of the other FZD proteins. Treatment of human pluripotent stem (hPS) cells with F7L6 initiates transcriptional programs similar to those observed during primitive streak formation and subsequent gastrulation in the mammalian embryo. This demonstrates that selective engagement and activation of FZD7 signaling is sufficient to promote mesendodermal differentiation of hPS cells.
Introduction
WNTs are highly conserved, lipid-modified secreted proteins with a broad range of activities throughout development and during adult tissue homeostasis (reviewed in Clevers et al., 2014). Deregulated WNT activity is associated with many pathologies, including degenerative and age-related diseases and cancer (reviewed in Nusse and Clevers, 2017). With the human genome encoding 19 WNTs and an equally large number of WNT receptors (Frizzled [FZD]1–10, LRP5 and 6, ROR1 and 2, RYK, PTK7 and more; reviewed in Niehrs, 2012), relatively little is known about signaling specificities between WNT ligands and their receptors. In experimental settings, activation of WNT signaling with either purified recombinant WNT proteins, such as Wnt3a, or with small molecule agonists, such as GSK3 inhibitors, frequently elicits the desired downstream transcriptional effect. Such observations have led to the prevailing and oversimplified view that specific WNT-receptor interactions are not as critical as the ensuing downstream signaling event. Because of the lack of selectivity between WNTs and their receptors in vitro and the paucity of purified and biologically active WNT proteins, it has not been possible to address whether engagement of a single WNT receptor is sufficient to elicit the complex biological processes affected by WNT.
Human pluripotent stem (hPS) cells provide a powerful in vitro system to study early processes of human development. The role of WNT signaling in the regulation of pluripotency is dependent on both the developmental stage of the cells and on the level of signaling. Experiments in mouse embryonic stem cells, which reside in the naive state (reviewed in Smith, 2017), have indicated a role for WNT/β-catenin signaling in self renewal and maintenance of pluripotency (Hao et al., 2006; Ogawa et al., 2006; Sato et al., 2004; ten Berge et al., 2011). In contrast, in hPS cells, which reside in the primed state and resemble epiblast-derived stem cells (Tesar et al., 2007), WNT/β-catenin signaling drives mesendodermal lineage specification (Davidson et al., 2012; Sumi et al., 2008; Tsakiridis et al., 2014), and an analysis of transcriptome-wide gene expression changes over a time-series treatment of hPS cells with Wnt3a reveals a transcriptional response reminiscent of primitive streak formation and ensuing gastrulation (Huggins et al., 2017). This is consistent with observations that short-term activation of Wnt/β-catenin signaling, via addition of Wnt3a, promotes differentiation toward definitive endoderm (DE) (D'Amour et al., 2006). Furthermore, a low level of endogenous WNT activity is detectable in undifferentiated hPS cells (Blauwkamp et al., 2012; Jiang et al., 2013) and is required for reprogramming of human fibroblasts to an induced pluripotent stem (iPS) cell state (Ross et al., 2014).
The WNT receptor FZD7 plays a prominent role in the regulation of pluripotency of hPS cells (Fernandez et al., 2014; Melchior et al., 2008), and its downregulation accompanies differentiation and exit from the pluripotent stem cell state. However, since hPS cells express several WNT receptors, including FZD2, 3, and 5 (Fernandez et al., 2014; Huggins et al., 2017), it is unclear whether Wnt3a promotes mesendodermal differentiation through FZD7 or whether other FZDs are required.
Here, we describe the development and characterization of a highly specific antibody and derivative single-chain variable fragment (scFv) to FZD7. Fusing this FZD7-specific scFv to an LRP6-directed scFv creates a bispecific binder, called F7L6, with potent WNT signaling activities that require FZD7 expression. Treatment of hPS cells with F7L6 elicits a transcriptional response similar to that observed for Wnt3a treatment, establishing that signaling through FZD7 is sufficient to promote mesendodermal differentiation.
Results
Design, specificity, and expression of F7L6
WNT/β-catenin signaling is initiated by WNT binding to and heterodimerizing the cell surface receptors FZD and either LRP5 or LRP6. To study the role of one specific FZD protein, FZD7, which is highly expressed in hPS cells (Fernandez et al., 2014; Melchior et al., 2008), we developed a bispecific binder to human FZD7 and LRP6, called F7L6. We generated a single-chain variable fragment (scFv) to FZD7 by repurposing the variable region of our FZD7-targeting antibody (F7-Ab), a chimeric human/mouse IgG1. This antibody was derived from a bacterially produced antigen binding fragment (Fab) described in a previous study (Fernandez et al., 2014). In parallel, we engineered an scFv to the third β-propeller domain of LRP6 based on previously published LRP6 antibodies (Ettenberg et al., 2010). We then engineered F7L6 by conjugating the FZD7-scFv and LRP6-scFv to a human IgG1 Fc to create the FZD7-specific Wnt mimetic F7L6 (complete amino acid sequence is provided in Figure 1—figure supplement 1). Additionally, we generated a FZD7-scFv-Fc (F7) and an LRP6-scFv-Fc (L6) (Figure 1A). We confirmed secretion of F7-Ab, F7, L6, and F7L6 from CHO cells by SDS-PAGE and Coomassie Blue staining and immunoblotting (Figure 1B). Using a dot blot with lysates of HEK293T cells carrying loss-of-function mutations in FZD1, 2, and 7 (F127-KO) (Voloshanenko et al., 2017) or in LRP6 (LRP6-KO) (Grainger et al., 2019) confirmed binding specificity of each of these molecules. FZD7-specific binders (F7 and F7L6) only reacted with lysates upon FZD7 overexpression (note that endogenous FZD7 expression in LRP6 KO is below the level of detection), whereas LRP6-specific binders (L6 and F7L6) only reacted with lysates of cells expressing LRP6 (Figure 1C).
To further characterize the specificity of F7-Ab, we mapped its binding site to the ‘neck’ region between the cysteine-rich domain (CRD) and the first transmembrane domain of FZD7 (Figure 1D). Despite the high degree of homology between human and mouse Fzd7 (hFZD7 and mFzd7, respectively), F7-Ab only reacts with hFZD7 (Figure 1E). Protein alignment (Figure 1D) indicated that this neck region was the only extracellular portion harboring differences between hFZD7 and mFzd7 (disregarding the amino-terminal signal sequence). A single amino acid change at position 188 from leucine to proline (P188) renders hFZD7 non-reactive to F7-Ab, whereas the corresponding amino acid change in mFzd7 (L188) restores F7-Ab reactivity (Figure 1E). Using fusion proteins between glutathione S-transferase (GST) and the FZD7 neck region followed by sequential shortening of the neck region, we mapped the F7-Ab epitope to an eight amino acid stretch containing L188 (Figure 1F). Since FZD proteins are highly conserved and several available FZD antibodies react with multiple FZD proteins (for example, OMP-18R5/Vantictumab reacts with FZD1, 2, 5, 7, and 8), we confirmed that F7-Ab does not cross-react with the other nine human FZD receptors (FZD1-6, 8–10) (Figure 1G). To rule out the possibility that some FZDs fail to react with F7-Ab because they do not reach the cell surface, we performed confocal microscopy. This analysis confirmed that all FZD proteins were expressed on the cell surface (Figure 1—figure supplement 2). Taken together, F7-Ab, and hence F7 and F7L6, are specific to hFZD7 and do not cross-react with any of the other nine human FZD proteins or with mFzd7.
F7L6 activates Wnt/β-catenin signaling
The bispecific binder, F7L6, activates WNT/β-catenin signaling by heterodimerizing FZD7 and LRP6 (Figure 2A). We confirmed F7L6 signaling activity using HEK293T cells stably transfected with the WNT reporter Super TOP-Flash (STF) (Veeman et al., 2003; Figure 2B). As expected, single binders to FZD7 (F7) or LRP6 (L6) did not activate signaling. As is the case with native WNT proteins, activity of F7L6 is potently augmented by addition of R-Spondin1 (RSPO1) (Kim et al., 2008; Figure 2C) with potency in the single-digit nanomolar range. Interestingly, addition of RSPO1 increased F7L6 activity approximately 10-fold, but only augmented the activities of Wnt3a or FLAg FP+P-L1+3, a previously published WNT mimetic (Tao et al., 2019), by twofold (Figure 2—figure supplement 1). A possible reason for this distinction is that F7L6 is selective for FZD7, whereas Wnt3a and FLAg FP+P-L1+3 are capable of interacting with multiple FZDs. Furthermore, both Wnt mimetics retain significant activity at sub-nanomolar concentrations, whereas Wnt3a’s activity is undetectable at such concentrations.
In contrast to WNT, which binds to the CRD of FZD (Janda et al., 2012), F7L6 binds to the neck region between the CRD and first transmembrane domain, allowing us to assess the requirement of the CRD in signaling. Interestingly, FZD7 lacking the CRD (CRD-less FZD7) activates signaling when heterodimerized to LRP6 with F7L6 while Wnt3a does not (Figure 2D), indicating that the CRD is dispensable for signaling. Furthermore, appending the F7-Ab epitope of eight amino acids (see Figure 1F) onto another FZD, FZD2, is sufficient for F7L6 to activate WNT/β-catenin signaling through FZD2 (Figure 2E). These data establish that heterodimerization of FZD7 or FZD2 with LRP6 is sufficient for pathway activation.
Treatment of mouse L-cells with Wnt3a activates downstream signaling, as assayed by cytoplasmic β-catenin stabilization (Shibamoto et al., 1998), whereas treatment with F7L6 does not (Figure 2F), consistent with the lack of reactivity of F7-Ab with mouse Fzd7 (see Figure 1E). However, overexpression of human FZD7 in these cells is sufficient to activate signaling by F7L6 (Figure 2F), thus confirming F7L6 signaling specificity through FZD7. The kinetics of β-catenin stabilization in L-cells is similar for Wnt3a and F7L6, being first detectable between 15 and 30 min after stimulation (Figure 2G). Lastly, using reverse transcription quantitative PCR (RT-qPCR), we demonstrated that F7L6, like Wnt3a, activates expression of the WNT target gene, SP5, in hPS cells, which express FZD7 (Fernandez et al., 2014; Melchior et al., 2008), in a dose-dependent manner (Figure 2H, Figure 2—figure supplement 2). In sum, F7L6 potently induces WNT/β-catenin signaling through FZD7.
F7L6 is a bispecific bivalent (or tetravalent) molecule, potentially capable of simultaneously engaging two FZD7 and two LRP6 molecules, thus leading to receptor oligomerization, which has been proposed to promote intracellular signalosome formation (DeBruine et al., 2017). To address whether FZD7-LRP6 heterodimerization alone is sufficient for signaling, we deleted the Fc portion to generate a single-chain F7L6 (F7L6-sc) that is predicted to engage one FZD7 and one LRP6 receptor. F7L6-sc activated Wnt signaling, albeit to a lesser extent than the tetravalent F7L6 protein (Figure 2I). Interestingly, the maximal response in the STF assay for F7L6 was approximately twice that for F7L6-sc. Size exclusion chromatography confirmed that F7L6-sc is monomeric in solution and does not form larger oligomers (Figure 2—figure supplement 3). These data demonstrate that a 1:1 association of FZD7 and LRP6 is sufficient to activate downstream signaling.
Differential effects of F7L6, Wnt3a, and CHIR on gene expression in hPS cells
Previous studies have shown that modulating WNT/β-catenin signaling is crucial for differentiation of hPS cells into mesendoderm and DE (D'Amour et al., 2006; Huggins et al., 2017; Jiang et al., 2013; Martyn et al., 2018; Yoney et al., 2018). Wnt3a and GSK3-inhibitors (e.g. CHIR99021, CHIR98014, BIO, Li) (Klein and Melton, 1996; Ring et al., 2003; Sato et al., 2004; Stambolic et al., 1996) are frequently used interchangeably as WNT pathway activators in these differentiation protocols (Gertow et al., 2013; Huang et al., 2017; Kumar et al., 2015; Loh et al., 2014; Naujok et al., 2014; Teo et al., 2014). Although Wnt3a is an established activator of WNT/β-catenin signaling, relatively little is known about the receptors that promote WNT-driven differentiation of hPS cells. We hypothesized that FZD7 mediates this process, because it transduces the Wnt3a signal (Dijksterhuis et al., 2015; Voloshanenko et al., 2017) and regulates the pluripotent state (Fernandez et al., 2014; Melchior et al., 2008). Here, we used F7L6 to selectively activate FZD7 signaling in hPS cells and analyzed global changes in gene expression by RNA sequencing (RNA-seq). We also examined the temporal kinetics of gene expression changes in response to a single, continuous dose of F7L6, Wnt3a, or CHIR98014 (CHIR) treatment.
Clustering of significantly differentially expressed genes (1814 genes, Figure 3—figure supplement 1, Supplementary file 1A) according to change in percent maximum reads per kilobase per million mapped reads (RPKM) revealed distinct waves and classes of gene expression, including activation of many known Wnt pathway genes (Figure 3A), activation of mesendodermal and endodermal genes, and downregulation of pluripotency- and ectoderm-associated genes (Figure 3B).
Gene expression changes in response to Wnt3a and F7L6 were overall highly similar, with robust activation of classic markers of mesendoderm (i.e. GSC, MIXL1, SP5 and T) and DE (i.e. LEFTY1, EOMES, CER1, FOXA2, CXCR4, and SOX17) and downregulation of pluripotency markers, such as POU3F1, FZD7, and PODXL (TRA-1–81). However, there were some notable differences between Wnt3a and F7L6 treatments. For example, induction of established WNT target genes, such as DKK1, AXIN2, WLS, WNT3, and NKD1, was more potent with Wnt3a than with F7L6. Also, induction of genes, such as T, SP5, GSC, and FOXA2, was slightly delayed for F7L6 relative to Wnt3a. These differences are potentially due to F7L6 solely engaging FZD7 and LRP6 to promote downstream signaling, while Wnt3a, in contrast, promotes signaling through multiple FZD-LRP5/6 heterodimers.
In contrast to Wnt3a and F7L6, CHIR treatment produced significantly different gene expression profiles, both in kinetics and amplitude (Figure 3A,B). For example, while all three treatments activated known WNT target genes (e.g. AXIN2, DKK1, LEF1, and SP5) and mesendodermal genes (e.g. T, MIXL1, GATA4, and EOMES), CHIR did so more robustly than Wnt3a or F7L6. Such differences are likely attributable to the fact that CHIR acts downstream of the receptor complexes and hence is not restricted by receptor availability. Furthermore, while all three treatments activated expression of Wnt pathway and mesendodermal genes, CHIR failed to activate expression of certain genes, in particular the DE markers FOXA2 and SOX 17 (and to a lesser extent CER and GSC), an effect we validated using RT-qPCR (Figure 3—figure supplement 2). Interestingly, expression of other DE markers, such as CXCR4, which encodes a cell surface receptor commonly used to enrich DE populations using flow cytometry (D'Amour et al., 2005; McGrath et al., 1999), was activated by CHIR (Figure 3B). These results indicate that sustained activation of WNT/β-catenin signaling by GSK3 inhibition with CHIR produces significant differences compared to selective pathway activation with either Wnt3a or F7L6.
F7L6, Wnt3, and CHIR treatment significantly altered expression of 525, 708, and 1646 genes, respectively (Figure 4, Supplementary file 1B–D). A majority of genes (428 of 805 genes) (Supplementary file 1E and F) in the F7L6 and Wnt3a lists overlapped, whereas expression changes of a large number of genes was unique to CHIR (1009 genes) (Supplementary file 1G), highlighting the substantial differences when the pathway is activated at the level of receptor rather than at the level of GSK3 inhibition. Using Gene Set Enrichment Analysis (GSEA) (Mootha et al., 2003; Subramanian et al., 2005) of the 391 genes (Supplementary file 1H) that changed in response to all three treatments indicated enrichment of diverse developmental processes, including Animal Organ Morphogenesis (Gene Ontology [GO]:0009887), Anatomical Structure Formation Involved In Morphogenesis (GO:0048646), Tube Development (GO:0035295), Circulatory System Development (GO:0072359), and Neurogenesis (GO:0022008) (Figure 4, Supplementary file 2A).
Importantly, GSEA of all genes altered by CHIR indicated enrichment in Neurogenesis (GO:0022008) and Neuron Differentiation (GO:0030182), both categories that were not among the top five GO terms for F7L6 and Wnt3a (Figure 4, Supplementary file 2B–D), indicating that CHIR treatment additionally promotes differentiation into ectodermal lineages. The 1009 differentially expressed genes unique to CHIR were additionally enriched in processes such as Biological Adhesion (GO:0022610), Locomotion (GO:0040011), Cell Junction (GO:0030054), and Positive Regulation Of Multicellular Organismal Process (GO:0051240) (Supplementary file 2E). In summary, our RNA-seq analyses confirmed that activation of FZD7 signaling by F7L6 was sufficient to modulate WNT target genes and induce expression of mesendoderm and DE differentiation in hPS cells. Overall, F7L6 and Wnt3a treatments induced similar (although not identical) patterns and temporal kinetics of gene expression, while CHIR caused much wider changes in global gene expression.
Early WNT target gene activation by F7L6 and Wnt3a in hPS cells
As revealed by our RNA-seq analysis, selective FZD7 pathway activation elicits a complex program of gene expression. Recent studies have provided evidence that certain Wnts recruit additional receptors to the Fzd-Lrp5/6 complex to increase signaling specificity. For example, the cell surface proteins Reck and Gpr124 promote Wnt7 signaling through Fzd (Cho et al., 2017; Eubelen et al., 2018). Furthermore, Egfr promotes a specific interaction of Wnt9a with Fzd9 and Lrp5/6 (Grainger et al., 2019). Likewise, Wnt3a may recruit additional, and currently unknown, co-receptors to activate signaling in hPS cells. In contrast, owing to its design, our Wnt mimetic F7L6 only recruits FZD7 and LRP6. To confirm that Wnt3a and F7L6 promote similar downstream signaling events, we more closely examined early WNT target gene activation and identified groups of genes with maximal activation at 3 hr and at 6 hr upon Wnt3a treatment (Figure 5A and B, Supplementary file 3A and B). We used RT-qPCR to validate changes in expression of several genes within these clusters in response to Wnt3a and F7L6 over a 6-hr time course in hPS cells. Both Wnt3a and F7L6 activated these early target genes with similar kinetics, with expression of AXIN2, SOX9, and ARL4A at 3 hr (Figure 5C–E) and NKX1-2, NKD1 and FGF18 at 6 hr (Figure 5F–H). These findings support that heterodimerization of FZD7 and LRP6 by F7L6 is sufficient to elicit a transcriptional response in hPS cells similar to that elicited by Wnt3a and likely does not involve the recruitment of additional receptors.
FZD7 activation promotes differentiation of hPS cells
Our RNA-Seq analysis indicated that F7L6, like Wnt3a, promotes mesendodermal differentiation. To monitor differentiation in real time, we employed human embryonic stem (hES) cells (H9/WA09) carrying reporter genes that mark meso- and endo-dermal differentiation: (1) for mesodermal differentiation, TBXT(T)-GFP, which harbors the gene encoding enhanced green fluorescent protein (eGFP) under control of the TBXT (T) promoter (Kita-Matsuo et al., 2009) and (2) for endodermal differentiation, SOX17-eGFP, which carries the eGFP gene in the SOX17 locus (Wang et al., 2011). We treated these reporter cell lines with a single, continuous dose of F7L6, Wnt3a or CHIR and monitored GFP expression for 5 days (Figure 6A). F7L6 potently activated expression of both reporters (Figure 6B,C, Figure 6—videos 1 and 2), with T-GFP expression first detectable at 24 hr (Figure 6B,D) and SOX17-eGFP at 36 hr (Figure 6C,E). Wnt3a likewise activated expression of both transgenes with similar temporal kinetics, however, their induction was significantly lower than with F7L6 (Figure 6B–E). This difference in potency between F7L6 and Wnt3a is likely due to the fact that Wnt3a is more unstable than F7L6 in these culture conditions. Treatment with CHIR yielded reporter gene expression distinct from both F7L6 and Wnt3a, with an early peak in T-GFP expression that subsequently declines (Figure 6B,D) and no induction of the SOX17-eGFP reporter (Figure 6C,E). This is consistent with our RNA-Seq (Figure 3B) and RT-qPCR results (Figure 3—figure supplement 2), again highlighting the distinct effects on differentiation of GSK3 inhibition versus FZD7 activation with either Wnt3a or F7L6.
FZD7 activation directs endodermal differentiation
Protocols to differentiate hPS cells specifically toward DE have been established and are widely used to generate mature endodermally derived cell populations, such as pancreas (D'Amour et al., 2006), liver (Touboul et al., 2010), and intestine (Workman et al., 2017). DE differentiation of hPS cells is induced by Activin/Nodal signaling, and Wnt3a addition on day 1 of differentiation increases the efficiency of mesendoderm specification and subsequent DE formation (D'Amour et al., 2006; Figure 7A). Using this protocol, we replaced Wnt3a treatment with F7L6 and monitored gene expression by RT-qPCR. As expected, expression of the pluripotency marker OCT4 (POU5F1) declined over the 3 days of differentiation with a more potent effect by F7L6 compared to Wnt3a (Figure 7B). Concurrently, expression of endodermal markers, SOX17, CXCR4 and FOXA2, increased upon treatment with F7L6 to a similar extent as Wnt3a (Figure 7B), indicating that activation of FZD7 alone is sufficient to promote DE formation.
FZD7 activation prevents differentiation toward the hematopoietic lineage
Hematopoietic stem and progenitor cells can be derived from hPS cells using defined culture conditions (Ng et al., 2008; Ng et al., 2005) and can be identified by dual expression of the cell surface markers CD34 and CD45 (Figure 8A). We previously showed that a specific WNT signal involving WNT9A and FZD9 increased the efficiency of this differentiation protocol, as monitored by an increase in the yield of CD34/CD45 double positive cells (Grainger et al., 2019; Richter et al., 2018). Given this highly selective requirement for WNT9A/FZD9 signaling, we reasoned that activation of FZD7 signaling with F7L6 would likely interfere with the differentiation of hPS cells toward the hematopoietic lineage. We found that F7L6 treatment at days 2–4 of differentiation, the treatment window we previously established to be critical for optimal differentiation, significantly reduced the number of CD34/CD45 double positive cells at day 14 relative to untreated or Wnt3a treated (Figure 8B). The percentages of CD34 and CD45 single positive cells were likewise adversely affected by F7L6 (Figure 8C and D). RT-qPCR demonstrated that F7L6 failed to induce expression of hematopoietic markers, including MESP1 and MIXL1 at day 4 and KDR at day 11 of differentiation (Figure 8E). Therefore, selective activation of FZD7 with F7L6 promotes differentiation into certain lineages, such as DE, but hinders differentiation into other lineages, such as the blood lineage.
Discussion
The WNT/β-catenin signaling pathway is well known for its critical roles in stem cell biology, including maintenance of pluri- and multi-potency and differentiation. However, the study of individual ligand-receptor pairs in regulating stem cell behavior has been complicated by two main factors: first, the apparent promiscuity in WNT-receptor interactions and, second, the large number of WNT and FZD gene family members. Here, we used hPS cells to analyze the role of a single WNT receptor complex, comprised of FZD7 and LRP6, in regulating the transition from undifferentiated and pluripotent stem cell to a mesendodermally restricted cell lineage. Specific activation of signaling was achieved using a WNT mimetic, F7L6, which only engaged and thereby heterodimerized FZD7 and LRP6. In contrast to the widely used Wnt3a protein, which engages multiple cell surface receptors, F7L6 exclusively binds to FZD7 and LRP6.
We demonstrated that selective engagement and heterodimerization of FZD7 and LRP6 with F7L6 is sufficient to drive the differentiation program promoted by Wnt3a in hPS cells. FZD7 is the most abundantly expressed FZD receptor in hPS cells, and its knockdown by RNA interference disrupts expression of pluripotency-associated genes, such as POU5F1/OCT4 and NANOG (Fernandez et al., 2014; Melchior et al., 2008). Likewise, inhibition of endogenous WNT processing in hPS cells using PORCN inhibitors also interferes with their pluripotency (Fernandez et al., 2014), suggesting that an endogenous WNT signal mediated by FZD7 is required to maintain hPS cells in an undifferentiated and pluripotent state. On the other hand, ectopic activation of WNT/β-catenin signaling, either with Wnt3a or with F7L6, leads to hPS cell differentiation along the mesendodermal lineage.
Whole transcriptome analysis revealed that Wnt3a and F7L6 elicit nearly identical transcriptional responses associated with early embryonic development, including formation of the primitive streak and gastrulation movements. Interestingly, we found that treatment with a GSK3 inhibitor (CHIR98014), which is commonly used interchangeably with Wnt3a, elicits significantly different effects. Most notably, we observed that over a 2-day treatment window, Wnt3a and F7L6 led to robust expression of the endodermal markers FOXA2 and SOX17, whereas the GSK3 inhibitor did not. Consistent with these observations, Kafri et al. found that GSK3 inhibition elicited significantly different dynamics and kinetics of β-catenin accumulation and localization compared to Wnt3a (Kafri et al., 2016). These findings indicate that specific activation of WNT/FZD signaling is important to recapitulate developmental processes.
The WNT mimetic F7L6 has several advantages over recombinant WNT proteins. First, native WNT proteins interact promiscuously with multiple receptors. For example, Wnt3a interacts with FZD1, 2, 4, 5, 7, 8, 9, and 10 (Dijksterhuis et al., 2015; Voloshanenko et al., 2017). In addition, the co-receptors LRP5 and LRP6 are functionally redundant (Kelly et al., 2004), indicating that WNT proteins can interact equivalently with either one. In contrast, WNT mimetics can be engineered to interact with greater selectivity for cell surface receptors. In vivo, WNT-FZD signaling specificity is regulated in part by co-receptors. For example, Reck, a glycosylphosphatidylinositol-anchored cell surface protein, acts together with Gpr124, a 7-transmembrane protein, to promote Wnt7a-Fzd-Lrp5/6 signaling to promote angiogenesis in the developing central nervous system (Cho et al., 2017; Cho et al., 2019; Eubelen et al., 2018; Vallon et al., 2018). Furthermore, we recently identified a novel function for Egfr in mediating Wnt9a-Fzd9-Lrp5/6 signaling during hematopoietic development in zebrafish (Grainger et al., 2019). The use of engineered proteins capable of engaging specific WNT receptors, such as F7L6, may overcome the need for the recruitment of specificity-conferring co-receptors, such as Reck, Gpr124 and Egfr, and thus replace the need for purification of biologically active WNT proteins.
A second advantage of WNT mimetics is that they can be assembled using recombinant proteins with favorable biochemical properties. Purification of native WNT proteins has proven extremely difficult, requiring 3 to 4 chromatography steps with yields of approximately 0.1 mg per liter of WNT conditioned medium (Willert, 2008). Furthermore, WNT proteins are lipid modified (Takada et al., 2006; Willert et al., 2003), rendering them highly hydrophobic and necessitating the addition of detergents to maintain their solubility in aqueous conditions. Consistent with its poor physicochemical properties, at low concentrations the activity of Wnt3a is rapidly extinguished and nearly undetectable. In contrast, F7L6, which is entirely comprised of engineered immunoglobulin sequences and hence is soluble and stable in standard biologically compatible buffers, retains signaling activity at significantly lower concentrations. In addition, by appending tags, such as Fc or 6x-His, F7L6 can be purified using a single affinity binding step, with an approximate yield of 2.5 mg per liter of conditioned media.
Several other groups have described recombinant proteins similar in design to F7L6 and capable of heterodimerizing FZDs and LRP5/6 (Chen et al., 2020; Dang et al., 2019; Janda et al., 2017; Miao et al., 2020; Tao et al., 2019). In contrast to these other WNT agonists that engage the CRD of FZD, the FZD7-binding arm of F7L6 binds the linker, or neck region, between the CRD and the first transmembrane domain of FZD7. Since this neck region is poorly conserved among the ten FZD receptors, we were able to exclusively engage FZD7 and no other FZD, including mouse Fzd7. This feature of F7L6 allowed us to demonstrate that the CRD is dispensable for signaling, since F7L6, but not Wnt3a, could induce signaling by a CRD-less FZD7. Therefore, aside from docking a WNT protein, the CRD has no additional role in activating the downstream signaling events.
We designed F7L6 as a bispecific bivalent (=tetravalent) molecule capable of simultaneously binding two FZD7 and two LRP6 cell surface receptors. This design possibly permits WNT receptor oligomerization, which in turn promotes intracellular signalosome formation, as previously proposed (DeBruine et al., 2017). Structural studies have revealed that WNTs can dimerize FZD CRDs either in a 1:2 or a 2:2 stoichiometry (Hirai et al., 2019; Nile et al., 2017). These higher order oligomers are hypothesized to promote optimal downstream signaling. Interestingly, a bivalent version, F7L6-sc, capable of forming 1:1 FZD7-LRP6 heterodimers, promotes downstream signaling, indicating that higher order oligomerization is not strictly required for signaling, consistent with previous findings by others (Miao et al., 2020). Of note, signal saturation is approximately twofold greater for the tetravalent F7L6 version compared to the bivalent F7L6-sc version, indicating that oligomerization may control signal amplitude, possibly by polymerization of the intracellular scaffolding molecules DVL and AXIN to drive signalosome formation. Therefore, tetrameric ligands, such as F7L6 and those developed by others (Chen et al., 2020; Tao et al., 2019), may augment and potentiate signaling by promoting receptor clustering. Additional experiments are needed to resolve the contribution of WNT versus DVL/AXIN to receptor oligomerization and signalosome formation.
Despite their potent stem cell activities in vivo and their potential as therapeutics in regenerative medicine, WNT proteins have not yet been approved for any clinical applications, owing in part to their poor physicochemical properties. Certain formulations, such as liposomal packaging (Morrell et al., 2008), have produced more stable and bioactive WNT3A, which is in early-stage clinical trials to treat patients undergoing posterolateral spinal fusion. However, such formulations will require purification of individual WNTs, which, to date, has only been achieved for WNT3A and WNT5A (Mikels and Nusse, 2006; Willert et al., 2003). WNT mimetics, such as the one described here, can be designed to target all WNT receptors and co-receptors and thus offers new opportunities for their development as therapeutics in regenerative medicine.
Materials and methods
Design of F7, L6 and F7L6
Request a detailed protocolThe sequence coding for the heavy (VH) and light (VL) chain of the variably fragment (Fv) portion of the FZD7 antibody (F7-Ab) were identified and used to design a single-chain variable fragment (F7-Ab scFv) that consists of the heavy chain Fv fused to the light chain Fv with a linker peptide (GGGGS)3. The F7-Ab scFv sequence was cloned into the pFuse-hIgG1-Fc2 mammalian expression vector (Invivogen), containing the IL-2 signal sequence and the IgG1 crystallizable fragment (Fc), to form the F7-Ab scFv-Fc (F7). The domains of the F7 construct are in the following order: IL-2 ss, F7-VH, (GGGGS)3, F7-VL, Fc. The coding sequence for the LRP6 scFv (United States Patent No.: US8,883,735, SEQ ID NO: 81 VL and SEQ ID NO: 82 VH) was similarly cloned into the pFuse mammalian expression vector to form the LRP6-Ab scFv-Fc (L6). The domains of the L6 construct are in the following order: IL-2 ss, L6-VL, (GGGGS)4, L6-VH, Fc. The coding sequence for the two scFv were fused with an additional flexible linker and cloned into the pFuse mammalian expression vector to create a bispecific scFv-Fc (F7L6). The domains of the F7L6 construct are in the following order: IL-2 ss, F7-VH, (GGGGS)3, F7-VL, (GGGGS)3, L6-VL, (GGGGS)4, L6-VH, Fc. Additional constructs containing a 6xHis tag between the IL-2 ss and the F7-VH domain with and without the Fc region were generated to express and produce F7L6 and F7L6-sc, respectively (Supplementary file 4).
Expression, purification, and characterization of recombinant proteins
Request a detailed protocolCHO cells (RRID:CVCL_0213) were transfected with plasmids encoding F7, L6, F7L6 and F7L6-sc, drug selected (1000 µg/mL zeocin or 4 µg/mL Puromycin, ThermoFisher Scientific) and expanded as clonal stable lines. Conditioned media (CM) were collected from confluent cultures every 3–4 days. Medium was replenished and CM was collected until the cells no longer adhered to the plate. CM were passed through a 0.22 µm filter (Genesee Scientific) and stored at 4°C until time of purification. CHO cells expressing the recombinant Wnt3a protein were similarly cultured and conditioned medium was collected, filtered, and stored. Wnt3a protein used in these studies was purified by four-step column chromatography as previously described (Willert, 2008). F7, L6, and F7L6 were purified from CM as follows: CM were applied to a protein G Sepharose (Biovision) column. The column was then extensively washed with phosphate buffered saline (PBS), and bound proteins were eluted with 0.1M glycine pH 2.5 and collected as 1 mL fractions into tubes containing 0.1 mL 2M Tris-Cl, pH8. Fractions containing the protein of interest were combined and dialyzed against PBS using dialysis tubing with a MWCO of 10 kDa (Thermo Scientific). After dialysis, the recombinant proteins were sterile filtered and aliquoted and stored at −80°C. The 6xHis-tagged version of F7L6 and F7L6-sc were purified as follows: CM from CHO cells expressing either F7L6 or F7L6-sc were adjusted to contain 30 mM imidazole and applied to a 1 mL immobilized metal affinity column charged with NiCl2 (HiTrap IMAC HP, Cytiva Lifesciences). The column was washed with 30 mM imidazole in PBS. Bound proteins were eluted with a linear gradient from 30 mM to 300 mM imidazole in PBS into 1 mL fraction. Fractions containing F7L6 were combined and dialyzed against PBS. Fractions containing F7L6-sc were further fractionated by size exclusion chromatography (SEC, Superdex 200 10/300 GL, Cytiva Lifesciences) to remove contaminant proteins. After dialysis or gel filtration, the recombinant proteins were aliquoted and stored at −80°C. SEC of F7L6 and F7L6-sc (as shown in Figure 2—figure supplement 3) was performed on a Superdex 200 10/300 GL column (Cytiva Lifesciences) in PBS, and 1 mL fractions were assayed in the presence of RSPO (100 ng/mL) for signaling activity using the STF assay. The previously published Wnt mimetic FLAg FP+P-L61+3 (Tao et al., 2019) was generously provided by Stephane Angers (University of Toronto, Canada).
Mapping of the F7-Ab epitope
Request a detailed protocolThe sequence encoding 42 amino acids (aa) of FZD7 containing the epitope for the F7-Ab (Glycine 168 to Serine 209) was cloned into the bacterial GST expression vector pGEX-4T3 (Cytiva Lifesciences). The epitope was shortened from either end through sequential cloning to produce 10 different length epitopes ranging from the original 42 aa to eight aa. The GST fusion proteins were expressed in BL21 cells and cell pellets were resuspended in protein sample loading buffer (2% SDS, 5% 2-mercaptoethanol, 10% glycerol, 62.5 mM Tris-Cl pH6.8) and boiled for 5 min (min) at 95°C. Proteins were visualized by Coomassie or by immunoblotting using the F7-Ab to determine binding to each epitope length.
Cell lines and culture conditions
Request a detailed protocolAll cell lines used in these studies were directly obtained from either ATCC (CHO, HEK293, L) or WiCell Research Institute (H1, H9). Identity of H1 and H9 was authenticated by karyotyping and functional assays, including EB formation and directed differentiation. All HEK293/HEK293T (RRID:CVCL_0045/RRID:CVCL_0063) lines and mouse L-cells (ATCC-CRL-2648) (RRID:CVCL_4536) were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and penicillin/streptomycin. HEK293/HEK293T cells stably transduced with the Wnt reporter Super TOP-Flash (STF, Addgene Plasmid #12456, RRID:Addgene_12456) were previously described (Bauer et al., 2013). HEK293T cells carrying knockout mutations in LRP6 were previously described (Grainger et al., 2019). HEK293T harboring mutations in multiple FZD genes (FZD1,2,7 [F127-KO] and FZD1,2,4,5,7,8) were kindly provided by Professor M. Boutros, Heidelberg University, Germany (Voloshanenko et al., 2017). CHO cells overexpressing Wnt3a were cultured in DMEM, 10% FBS and penicillin/streptomycin, and Doxycycline (250 ng/mL) was added to induce Wnt3a expression.
All experiments using human pluripotent stem cell lines (hPS cells) were approved under IRB/ESCRO protocol number 100210 (Principal investigator: K.W.). HPS cell lines H1 (WA01, NIH Registration Number 0043, RRID:CVCL_9771) and H9 (WA09, NIH Registration Number 0062, RRID:CVCL_9773) were obtained from WiCell and cultured in E8 culture medium (Chen et al., 2011) on Matrigel (BD Biosciences). The H9 reporter lines, H9 SOX17:GFP (Wang et al., 2011) and H9 T-GFP (Kita-Matsuo et al., 2009) were cultured in E8 culture medium, passaged with TrypLE Express (Gibco), and seeded with Rock inhibitor Y-27631, 5 μM (Tocris). The H9 SOX17:GFP cell line was kindly provided by material transfer agreement by Dr. Seung Kim (Stanford School of Medicine). iPS cell lines were kindly provided by Dr. Dan Kaufman (UC San Diego) (Li et al., 2018). The iPS cells were cultured in mTeSR1 (StemCell Technologies), passaged with TrypLE Express, and seeded with Rock inhibitor. All hPS and iPS cells were fed fresh media daily.
Differentiation of hPS cells
Non-directed differentiation
Request a detailed protocolReporter hPS cell lines, H9 TBXT(T)-eGFP and H9 Sox17-eGFP, were passaged into a Matrigel-coated 12 well plate at 10,000 cells per well. When cells reached 40–60% confluence, media in each plate were replenished and the cells were treated with F7L6 (10 nM), Wnt3a (5 nM), CHIR98014 (250 nM, Sigma-Aldrich), or an equivalent volume of Wnt storage buffer (PBS, 1% CHAPS, 1 M NaCl). Plates were placed in the IncuCyte Life Cell Analysis System (Sartorius) and cultured without media changes for 5 days at 37°C. Phase and GFP images were recorded every 3 hr for a total of 120 hr after treatment. The IncuCyte software was used to process the images and generate the integrated GFP intensity (GCU x μm2/image).
Endoderm differentiation
Request a detailed protocolH1 cells were differentiated to endoderm as previously described (D'Amour et al., 2006). Initiated on days 4–6 after passage (depending on culture density), sequential, daily media changes were made for the entire differentiation protocol. After a brief wash in PBS, cells were cultured in RPMI, Activin A (100 ng/mL, R and D Systems) and a treatment for the first day. Treatments were F7L6 (10 nM), Wnt3a (5 nM), CHIR (250 nM) or an equivalent volume of Wnt storage buffer. The next day the medium was changed to RPMI with 0.2% vol/vol FBS and Activin A (100 ng/mL), and the cells were cultured for 2 additional days. Cells were collected at day three and total RNA was isolated and analyzed by RT-qPCR.
Hematopoietic differentiation
Request a detailed protocolHPS cells were differentiated to the hematopoietic lineage as previously described (Ng et al., 2008; Ng et al., 2005). The day prior to differentiation, iPS cells were passaged with TrypLE Express at a high density in mTeSR with Rock inhibitor. On Day 0 of the differentiation, iPS cells were dissociated into a single cell suspension and plated in a 96-well U-bottom plate (Genesee Scientific) at 3000 cells per well in APEL2 (StemCell Technologies) with 40 ng/mL BMP4 (R and D Systems), 40 ng/mL SCF (R and D Systems), 20 ng/mL VEGF (R and D System), 10 ng/mL FGF-2 (StemCell Technology), and 5 nM Rock inhibitor. On Day 2 of the differentiation, the embryoid bodies (EB) were treated with F7L6 (5 nM), Wnt3a (5 nM), CHIR (250 nM) or an equivalent volume of Wnt storage buffer. On Day 4 of the differentiation, the media was removed from the EBs and replaced with fresh APEL2 with BMP4, SCF, VEGF, and FGF-2. On Day 7 of the differentiation, EBs were transferred to gelatin coated plates. On Day 8 of the differentiation one volume of APEL with BMP4, SCF, and VEGF was added. EBs were dissociated on Day 14 of the differentiation for analysis by flow cytometry.
Immunoblotting by dot blot
Request a detailed protocolCells were lysed in TNT buffer (1% Triton X-100, 150 mM NaCl, 50 mM Tris HCl, pH 8) with protease inhibitors (Roche). Five µL of whole cell lysate (approximately 40 µg) was dotted directly onto nitrocellulose membrane and allowed to dry completely. The nitrocellulose membrane was incubated in blocking buffer (TBST [20 mM Tris-Cl pH8, 150 mM NaCl, 0.2% Tween-20], 1% BSA, 3% non-fat dry milk) for 30–60 min at room temperature (RT), and then incubated in primary antibody for 1 hr at RT or overnight at 4°C. The conditioned media collected from F7, L6, and F7L6 expressing CHO cells was mixed 1:1 with blocking buffer and used as a primary antibody solution. The nitrocellulose membrane was washed three times in TBST, and then incubated in goat anti-human HRP secondary antibody (ThermoFisher Scientific) at 1:10,000 dilution in blocking buffer for 45–60 min at RT. The nitrocellulose membrane was washed three times in TBST before enhanced chemiluminescent (ECL) detection by SuperSignal West Dura Western Blot Substrate (ThermoFisher Scientific) and exposure to autoradiography film.
Immunoblotting by western
Request a detailed protocolTo obtain whole cell lysate, cells were lysed in TNT buffer with protease inhibitors. Protein concentrations were determined using Pierce Coomassie (Bradford) Protein Assay Kit (ThermoFisher Scientific). Protein sample loading buffer (62.5 mM Tris-Cl pH6.8, 2% SDS, 10% glycerol, 5% 2-mercaptoethanol, bromophenol blue) was added to 20 μg total protein. Samples were denatured at 95°C for 5 min unless samples were intended to be blotted for FZD proteins. Cell lysates were resolved by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), transferred to nitrocellulose membrane, and incubated for 30–60 min in blocking buffer. Primary antibody incubations were done overnight at 4°C. Primary antibodies and dilutions used: anti-FZD7 (F7-Ab) (1:1000 of 1 mg/mL stock), V5 tag antibody (GeneTex/GTX628529/1:4,000), monoclonal anti-β-catenin antibody (Sigma-Aldrich/C7207/1:2,000) (RRID_AB_476865), anti-β-actin antibody (Sigma-Aldrich/A2228/1:5,000) (RRID:AB_476697). All western blots were washed three times in TBST prior to incubation in secondary antibody for 45–60 min at RT. Secondary antibodies and dilutions used: goat anti-human IgG HRP-conjugated (ThermoFisher Scientific/62–8420/1:20,000) (RRID:AB_88136), goat anti-mouse IgG HRP-conjugated (Southern Biotech/1030–05/1:20,000) (RRID:AB_2619742). All western blots were washed three times in TBST before protein detection by Millipore Sigma Luminata Forte Western HRP substrate and exposure to autoradiography film.
Super TOP-Flash (STF) luciferase assays
Request a detailed protocolCells were lysed in luciferase assay buffer (100 mM K-PO4 buffer pH 7.8, 0.2% Triton X-100) and transferred to a black-walled 96-well plate. One hundred µL of luciferase assay cocktail was added to each well of lysate 25 mM Tris-Cl pH 7.8, 15 mM MgSO4, 10 mM ATP (Sigma-Aldrich), 65 µM BD D-luciferin-Potassium Salt (Fisher Scientific). Luciferase assay readouts were performed on a Promega GloMax Discover Microplate Reader.
Real-time quantitative polymerase chain reaction (RT-qPCR)
Request a detailed protocolRNA expression was measured by RT-qPCR. RNA was extracted using TRIzol Reagent (ThermoFisher Scientific) and Direct-zol RNA MiniPrep Kit (Zymo Research). cDNA was generated using 50 ng RNA and iScript Reverse Transcription Supermix (Bio-Rad), then diluted 1:10 in UltraPure DNase/RNase-Free Distilled Water (ThermoFisher). RT-qPCR was performed using iTaq Universal SYBR Green Supermix (Bio-Rad) according to the manufacturer’s recommendations, and a two-step amplification CFX_2stepAmp protocol on a Bio-Rad CFX384 Touch Real-Time PCR Detection System. Data was analyzed using the 2−ΔΔCt method (Schefe et al., 2006). All gene expressions were normalized to the expression of RPL13A. The following RT-qPCR primers were used (Gene ID, forward primer, reverse primer):
ARL4A, GGCGATTTAGTCAAGAGGAT, GCTCTTCTCAACACACTACA
AXIN2, TATCCAGTGATGCGCTGACG, CGGTGGGTTCTCGGGAAATG
CXCR4, ACTACACCGAGGAAATGGGCT, CCCACAATGCCAGTTAAGAAGA
FGF18, GAGGAGAACGTGGACTTCCG, ACCTGGATGTGTTTCCCACT
FOXA2, GGAGCAGCTACTATGCAGAGC, CGTGTTCATGCCGTTCATCC
KDR, CCTGTATGGAGGAGGAGGAAGT, CAAATGTTTTTACACTCACAGGCCG
MESP1, CTCTGTTGGAGACCTGGATG, CCTGCTTGCCTCAAAGTG
MIXL1, TCCAGGATCCAGCTTTTATTTTCT, GAGGATAATCTCCGGCCTAGC
NKD1, TCCCAACCTAGAAACCTTAG, AGAAGAAGGAGAAGGAAGAG
NKX1-2, GTAGAAGAGAGGGAATAGGGAGAG, AGCAGCAGAAGTCCAAAGTC
POU5F1, CTTGAATCCCGAATGGAAAGGG, GTGTATATCCCAGGGTGATCCTC
RPL13A, CCTGGAGGAGAAGAGGAAAGAGA, TTGAGGACCTCTGTGTATTTGTCAA
SOX9, GACACAAACATGACCTATCC, GATTCTCCATCATCCTCCAC
SOX17, GTGGACCGCACGGAATTTG, GGAGATTCACACCGGAGTCA
SP5, TCGGACATAGGGACCCAGTT, CTGACGGTGGGAACGGTTTA
Flow cytometry
Request a detailed protocolHPS cell-derived EBs were dissociated with TrypLE Express at 37°C for 5–10 min, periodically triturated using P1000 pipette. Dissociated EBs were resuspended in FACS buffer (PBS, 1 mM EDTA, and 0.5% FBS) and passed through a cell strainer (Corning). Cell suspensions were pelleted at 200xg for 3 min and resuspended in 100 µL of FACS buffer (approximately 1 × 105 cells). Cells were incubated on ice for 30 min in the following primary fluorophore-conjugated antibodies (Antibody/Vendor/Catalog#/Concentration): APC anti-human CD34/Biolegend/343608/1:100 (RRID:AB_2228972); PE anti-human CD45/Biolegend/304008/1:100 (RRID:AB_314396); APC mouse IgG2a k isotype control/Biolegend/400222/1:100; PE mouse IgG1, k isotype control/Biolegend/400112/1:100 (RRID:AB_2847828). Cells were washed with 2 mL of FACS buffer and spun down at 200xg for 3 min. Cells were resuspended in 500 mL of FACS buffer with 0.5 µg/mL DAPI (Cell Signaling Technology). The cell suspensions were analyzed using the FACS Fortessa (BD Biosciences) and the resulting FSC files were processed with the FlowJo software (BD Biosciences).
Immunofluorescence of overexpressed FZD
Request a detailed protocolHEK293T F127-KO cells were plated in a 24 well-plate with coverslips. Cells were transfected with V5-tagged FZD receptors and 2 days post-transfection, cells were fixed with 4% formaldehyde for 15 min at RT. Cells were washed twice with 1X PBS and permeabilized with 1X PBS with 0.5% Triton X-100 for one hour at RT. Following permeabilization, cells were blocked with 0.5% BSA in 1X PBS with 0.1% Triton X-100 (1X PBST) for 1 hr at RT. Cells were stained with V5 antibody (Genetex/GTX628529/1:500) in 1X PBST at 4°C overnight. Cells were washed three times with 1X PBST for 10 min each. Alexa-Fluor 488 Goat anti-Mouse IgG (Invitrogen/A11001/1:1000) (RRID:AB_2534069) in 1X PBST was used as the secondary antibody. Secondary antibody staining was done for 2 hr at RT. Following the secondary staining, cells were washed three times with 1X PBST for 10 min each. Coverslips were mounted on glass slides and imaged on Nikon Eclipse Ti2-E microscope with SR HP Apo TIRF 100x NA 1.49 objective.
Transcriptome analysis (RNA-seq)
Request a detailed protocolTotal RNA from cells was extracted using TRIzol Reagent and Direct-zol RNA MiniPrep Kit, according to manufacturer recommendations. cDNA library preparation and sequencing were done by Novogene Co., Ltd. A 250–300 bp insert cDNA library was generated by using the NEBNext Ultra RNA Library Prep Kit for Illumina (New England Biolabs, Inc). Transcriptome sequencing was performed on an Illumina NovaSeq 6000. TopHat (RRID:SCR_013035) and Cufflinks (RRID:SCR_014597) (Trapnell et al., 2013; Trapnell et al., 2012) were used to perform differential gene expression analysis of RNA-seq experiments. Briefly, sequencing reads were quality filtered, mapped, and aligned to the reference human genome (hg19) with TopHat and Cuffdiff (RRID:SCR_001647) was used to calculate gene expression levels as reads per thousand transcript bases per million reads mapped (RPKM). Statistically significant changes in gene expression were obtained from RPKM values. Genes were clustered by expression pattern and principal component analysis was performed in Genesis (Sturn et al., 2002) (RRID:SCR_015775). Gene ontology was performed using GSEA (Mootha et al., 2003; Subramanian et al., 2005) (RRID:SCR_003199).
Data availability
The RNA-seq and ChIP-seq data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE158121.
-
NCBI Gene Expression OmnibusID GSE158121. Selective activation of FZD7 in human pluripotent stem cells.
References
-
Efficient differentiation of human embryonic stem cells to definitive endodermNature Biotechnology 23:1534–1541.https://doi.org/10.1038/nbt1163
-
Production of pancreatic hormone-expressing endocrine cells from human embryonic stem cellsNature Biotechnology 24:1392–1401.https://doi.org/10.1038/nbt1259
-
Receptor subtype discrimination using extensive shape complementary designed interfacesNature Structural & Molecular Biology 26:407–414.https://doi.org/10.1038/s41594-019-0224-z
-
Assembly and architecture of the Wnt/β-catenin signalosome at the membraneBritish Journal of Pharmacology 174:4564–4574.https://doi.org/10.1111/bph.14048
-
Systematic mapping of WNT-FZD protein interactions reveals functional selectivity by distinct WNT-FZD pairsJournal of Biological Chemistry 290:6789–6798.https://doi.org/10.1074/jbc.M114.612648
-
Crystal structure of a mammalian Wnt-frizzled complexNature Structural & Molecular Biology 26:372–379.https://doi.org/10.1038/s41594-019-0216-z
-
R-Spondin family members regulate the wnt pathway by a common mechanismMolecular Biology of the Cell 19:2588–2596.https://doi.org/10.1091/mbc.e08-02-0187
-
Embryonic expression and function of the chemokine SDF-1 and its receptor, CXCR4Developmental Biology 213:442–456.https://doi.org/10.1006/dbio.1999.9405
-
The WNT receptor FZD7 contributes to self-renewal signaling of human embryonic stem cellsBiological Chemistry 389:897–903.https://doi.org/10.1515/BC.2008.108
-
The complex world of WNT receptor signallingNature Reviews Molecular Cell Biology 13:767–779.https://doi.org/10.1038/nrm3470
-
Synergistic action of Wnt and LIF in maintaining pluripotency of mouse ES cellsBiochemical and Biophysical Research Communications 343:159–166.https://doi.org/10.1016/j.bbrc.2006.02.127
-
Quantitative real-time RT-PCR data analysis: current concepts and the novel "gene expression's CT difference" formulaJournal of Molecular Medicine 84:901–910.https://doi.org/10.1007/s00109-006-0097-6
-
Genesis: cluster analysis of microarray dataBioinformatics 18:207–208.https://doi.org/10.1093/bioinformatics/18.1.207
-
Embryonic stem cells require wnt proteins to prevent differentiation to epiblast stem cellsNature Cell Biology 13:1070–1075.https://doi.org/10.1038/ncb2314
-
Differential analysis of gene regulation at transcript resolution with RNA-seqNature Biotechnology 31:46–53.https://doi.org/10.1038/nbt.2450
-
Isolation and application of bioactive wnt proteinsMethods in Molecular Biology 468:17–29.https://doi.org/10.1007/978-1-59745-249-6_2
Article and author information
Author details
Funding
National Institutes of Health (R35GM134961)
- Karl Willert
National Institutes of Health (S10OD026929)
- Karl Willert
National Cancer Institute (T32 CA067754 graduate student fellowship to Myan Do)
- Myan Do
NIH (S10OD025060)
- Karl Willert
National Institutes of Health (R01HL135205)
- Karl Willert
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
Research reported in this publication was supported by the National Institute of General Medical Sciences of the National Institutes of Health under Award Number R35GM134961 (awarded to KW). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. MD was supported by in part by a Cancer Biology, Informatics, Omics Training Program from the National Cancer Institute (T32 CA067754). This work was additionally made possible by the CIRM Major Facilities grant (FA1-00607) to the Sanford Consortium for Regenerative Medicine. This publication includes data generated at the UCSD Human Embryonic Stem Cell Core Facility utilizing the IncuCyte Zoom Live Cell Analysis System that was purchased with funding from a NIH Shared Instrumentation Grant (S10OD025060), and at the UC San Diego IGM Genomics Center utilizing an Illumina NovaSeq 6000 that was purchased with funding from a NIH Shared Instrumentation Grant (S10OD026929).
Copyright
© 2020, Gumber 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
-
- 1,913
- views
-
- 255
- downloads
-
- 19
- 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
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.
-
- Developmental Biology
We have generated single cell transcriptomic atlases of vomeronasal organs (VNO) from juvenile and adult mice. Combined with spatial molecular imaging, we uncover a distinct, previously unidentified class of cells that express the vomeronasal receptors (VRs) and a population of canonical olfactory sensory neurons in the VNO. High-resolution trajectory and cluster analyses reveal the lineage relationship, spatial distribution of cell types, and a putative cascade of molecular events that specify the V1r, V2r, and OR lineages from a common stem cell population. The expression of vomeronasal and olfactory receptors follow power law distributions, but there is high variability in average expression levels between individual receptor and cell types. Substantial co-expression is found between receptors across clades, from different classes, and between olfactory and VRs, with nearly half from pairs located on the same chromosome. Interestingly, the expression of V2r, but not V1r, genes is associated with various transcription factors, suggesting distinct mechanisms of receptor choice associated with the two cell types. We identify association between transcription factors, surface axon guidance molecules, and individual VRs, thereby uncovering a molecular code that guides the specification of the vomeronasal circuitry. Our study provides a wealth of data on the development and organization of the accessory olfactory system at both cellular and molecular levels to enable a deeper understanding of vomeronasal system function.