Validation of fin cell exposure to Xenopus egg extract
The mesenchymal cell preparation and treatment that were set up in a previous study [31] included plasma membrane permeabilization with digitonin, permeabilized cell exposure to egg extract for 1h, and plasma membrane resealing (Supplementary Fig. S1). Penetration of egg factors per se was not tested here, because this would have required cell fixation. However, all cells displayed the phenotypic characteristics of permeabilized cells and egg extract-treated cells that were described previously [31]: their nuclear membrane was more contrasted after permeabilization than in control cells, their adhesion capacity lessened during egg extract exposure and remained very low during the resealing step and the first 24 h of culture, and they all adopted a round and refracting morphology after resealing. Taken together, our observations indicate that all treated cell batches in the present study did incorporate egg extract.
Treated cells need a suitable culture medium to survive and proliferate.
The culture phase of the treated cells had to be mastered, so that the treated cells could undertake their new cellular program. When the conventional L15 medium was used, we observed that from the second day of culture on, many treated cells displayed a cubic shape (Fig 1, left picture) that contrasted with the much more elongated control cells (Fig 1, central picture). However, after 7 days of culture in L15, the treated cell density decreased, and many cells detached from the culture plate (Fig 1, left pictures), whereas control cells kept proliferating (Fig 1, inset d7). This inability of the treated cells to survive in L15 medium provided a first indication that the egg extract treatment had induced some changes in the treated cell physiology.
Therefore, we sought for a culture medium that would sustain treated cells survival and proliferation, while maintaining their modified state. In mammals, somatic cells treated with egg extract were reported to be cultured in embryonic stem (ES) cells medium containing Leukemia Inhibitory Factor (LIF) and other complements aimed to prevent cell differentiation [29,25,27]. However, fish ES-like cells are known to be independent from LIF (reviewed in [35]). Furthermore, the maintenance of an undifferentiated state in zebrafish and medaka ES-like cultured cells was reported to require a medium enriched with fish serum and species-specific embryo extracts, known as ESM4 medium [35,36]. We therefore tested the ESM4 medium enriched with goldfish embryo extracts (supplementary Table S1). After 2 days in this new culture medium, the cubic shape of the treated cells seen in L15 was maintained in ESM4 (Fig 1, right pictures). The elongated shape of the control cells was not changed either, indicating that ESM4 has no effect of its own on the shape of the cultured cells. Most interestingly, the treated cells cultured in ESM4 were able to proliferate over longer culture time compared to culture in L15 medium. They showed an increased cell density at day 7, and debris and floating cells were no longer observed (Fig 1, right picture). Additionally, they maintained their specific cubic morphology. This demonstrates further the favorable effect of ESM4 medium on the growth of the modified treated cells.
Changes in gene expression eight days after egg extract treatment.
Clustering of the differentially expressed genes (DEGs)
Analysis of the microarray data revealed that 2,286 goldfish genes out of the 52,362 genes on the microarray were differentially expressed between treated and control cells (fold change > 2). Additionally, hierarchical clustering analysis of the differentially expressed genes (DEGs) showed a clear segregation between treated and control samples (Fig 2, upper dendrogram). This demonstrated that the treated cells transcriptome was modified by egg extract treatment and that the consequences were detectable after 8 days of culture. Differentially expressed genes between treated and control cells showed a distribution into two clusters on the heatmap. Cluster I gathers 872 genes (encompassing 38 % of the DEG) that showed upregulation in the treated cultured cells. Cluster II comprises 1 414 genes (62 % of the DEG) that were down regulated in the treated cells. Genes in each cluster are listed in Supplementary Table S2.
Segregation of the treated samples according to egg-extract batches.
Although all treated samples segregated together and showed the same expression profile clustering, it is noteworthy that samples T1 to T3 segregated together, apart from the 4 other samples (T4-T7) (Fig 2, upper dendrogram). One possible explanation lies in the egg-extract batches that were used for the different samples. Indeed, the extracts were all prepared from freshly spawned MII stage eggs, each extract being obtained from the spawn of a different female. We cannot exclude that the individual extracts presented some quality variations one from another, notably because of the instability of MII stage in spawned eggs [37]. In order to validate the egg extract stage, we used two MII markers: Greatwall whose phosphorylated forms prevent mitosis/meiosis exit [38], and Cyclin B whose degradation characterizes mitosis/meiosis exit. Egg extracts arrested at MII stage all displayed a specific western blot profile (supplementary Fig. S2): Greatwall (Gwl) was phosphorylated and stable over incubation time, and cyclin B (CycB) content was high and stable as well. Upon in vitro induction of MII exit by Ca2+, Greatwall was successfully dephosphorylated and cyclin B underwent degradation. Contrarily to these well-defined MII egg extracts, some egg extracts showed Greatwall dephosphorylation and Cyclin B degradation, indicating that they had initiated MII exit (MII late stage). Interestingly, the egg extracts used to treat T1 to T3 samples where in MII stage whereas those of T4 to T7 samples had initiated MII exit to some extent (Supplementary Fig. S2). As a conclusion, the sample segregation in the treated cells was likely related to the extract stages (MII and MII-late). This highlights the importance of a careful characterization of the Xenopus egg extracts. Although the sample number in each category was low, we still performed a fold change analysis between the two groups. We observed that 83% of the DEGs between MII and MII-late extract groups had low fold changes (< 6), and only 52 genes had fold changes above 6, among which only 9 genes were above 20. Besides, no significant or straightforward biological processes were identified via the GO terms analysis, and no marker gene of any specific biological significance emerged from a gene to gene scouting. To conclude, and within the limits of this small sampling, the egg extract stage did not thoroughly affect cellular response, and the two clusters of up- and downregulated genes were observed in all 7 treated cells batches irrespective of the egg extract that was used.
Gene Ontology (GO) analysis of the differentially expressed genes after egg extract treatment.
GO analysis was a perquisite in order to process our DEG list into functions and biological significance. For this purpose, we had first to translate the goldfish gene identifiers into those of the closest species whose genome is well annotated in the GO databases, the zebrafish. This artificially reduced the number of DEG, because the zebrafish did not undergo the genome duplication reported in the Cyprininae sub-family to which goldfish belongs. Only 1533 zebrafish genes (i.e. 67% of total goldfish DEG) were retained for subsequent annotations in GO. Of these, 591 annotated genes were up-regulated (cluster I) and 942 annotated genes were down regulated (cluster II) in the treated cells. GO analysis conducted with the WebGestalt web tool [39] showed that biological regulation and metabolic process were the most represented terms (supplementary Table S3) among biological process GO terms. Surprisingly, no straightforward reprogramming processes such as chromatin remodeling, stem cells, transcription factors, or pluripotency could be emphasized in GO terms after a statistical over-representation analysis. However, several other biological process terms significantly enriched in the GO terms list deserve specific attention
Deregulation of TGFβ and Wnt signaling pathways after egg extract treatment.
Our work is reporting gene expression variation, but GO databases and related publications on gene function report mainly protein functions. It is therefore the protein writing nomenclature that will be used in the following sections dedicated to GO interpretation. The most significant GO term obtained from the cluster of up-regulated genes is the cell surface receptor signaling pathway (Fig 3A). The data mapping showed that this GO term was linked to highly significant child GO terms that are transforming growth factor beta (TGFβ) receptor signaling pathway, and Wnt signaling pathway together with regulation of canonical Wnt signaling pathway. This result was consistent with the KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis performed on the same set of DEG data, that also showed that both TGFβ and Wnt signaling pathways reached a significant level of enrichment among all the database terms (Fig 3B).To add on to the highlighting of these 2 specific pathways, we also observed an enrichment in the GO terms related to the MAPK / ERK cascade (Fig 3A), known to be one of the non-canonical pathways activated by TGFβ [40]. Thus, the GO analysis based on the cluster of up-regulated genes clearly highlighted the TGFβ and Wnt signaling pathways as major ones being affected in the cells exposed to egg extract reprogramming factors.
TGFβ signaling
TGFβ signaling is involved in numerous biological processes related to embryonic development. We then checked the actors of the TGFβ pathway present in our goldfish DEG list (irrespective of their up or down regulation). TGFβ belongs to the superfamily of the growth factors, divided into several subfamilies including TGFβs, and Bone Morphogenetic Proteins (BMPs). For TGFβ signal to be transduced, the TGFβ ligand binds type II receptors. Ligand - type II receptor complex triggers the recruitment of TGFβ type I receptor, and the dimerized receptors subsequently activates specific Smad proteins, able to induce transcription of the TGFβ target genes [41,42]. Beyond signaling pathways involving Smads, known as canonical TGFβ pathways, other pathways independent of Smads are also controlled by TGFβ, including the MAPK Erk1 / ERk2 pathway identified above by GO analysis (Fig. 3A). We therefore analyzed the expression profile of these TGFβ actors and their biological partners.
We found that some TGFβ and BMP ligands together with type I receptors were upregulated in treated cells compared to controls (Table 1, TGFβ Effectors). However, this upregulation is unlikely stimulating the TGFβ signaling pathway, because the key actors binding TGFβ that are the TGFβ type II receptors did not change their expression pattern in treated cells. Besides, most other actors of the TGFβ signaling pathway identified in this study were affected in the direction of a TGFβ signaling inhibition in the treated cells, namely inhibitors upstream of TGFβ signaling that were upregulated in the treated cells (Table 1, TGFβ Inhibitors). Among them, we identified extracellular inhibitors (lft2, nog1, nog2, grem2a, grem2b) and membrane inhibitors (Bambia and Bambib) that are binding to TGFβ and BMP ligands. Such binding prevents TGFβ and BMP to attach to their own receptors, thereby preventing signal transduction activity [43,44]. Beyond these inhibitors, we also found intracellular inhibitors (involved in TGFβ canonical signaling pathway) which included specific smads (smad6a, smad6b, smad7, smad9) and the ubiquitin ligase smurf2 (Table 1). The combined action of Smad7 and Smurf 2 is known to induce TGFβ type I receptor degradation by the proteasome [42,45], leading to inhibition of the TGFβ canonical pathway. Finally, spry1, sry4, and dusp6 genes, inhibiting the MAPK / ERK pathway (non-canonical TGFβ pathway), were also found upregulated in the treated cells (Table 1). In all, our gene to gene analysis of the expressional changes of TGFβ actors, including the non-canonical MAPK / ERK pathway, indicated that fin cells exposure to egg extract induced an overall inhibition of the TGFβ signaling pathway.
Mesenchymal-epithelial transition
In mammals, one consequence of TGFβ signaling inhibition is the induction of a mesenchymal-epithelial transition (MET), considered to be a hallmark of iPSC early phase reprogramming, and described as crucial for reaching pluripotency [46–50]. MET is characterized by the loss of mesenchymal markers and by the activation of genes determining epithelial fate [47]. We therefore investigated whether inhibition of TGFβ signaling in our fish treated cells was also associated with changes in MET marker genes. We found that many mesenchymal marker genes were downregulated, among which several members of the collagen family, matrix-metallo protease (mmp9) and fibronectin (fn1) (Table 2). The fn1 gene was the most strongly affected (- 44 fold change). This was associated with the concomitant upregulation of several epithelial marker genes such as cadherins (pcdh1 cadherin-like 1, pcdh12, cdh18, cdh24b), cytokeratins (krt15, krt18), and cell junction proteins such as pkp3b, cldn5a, tjp1a and cx43 (Table 2). Regarding the gap junction component cx43, it is known to be specifically enriched in epithelial cells and iPSCs, and its ectopic expression and gene upregulation has been associated with an increase in reprogramming efficiency by facilitating MET [51]. Last, the transcription factor zeb1 known to induce EMT (epithelial-mesenchymal transition) [52], ie the reverse of the MET, was downregulated in treated cells (Table 2). In all, the observed expressional changes suggest the initiation of a MET program in the treated cells. This expressional profile is in accordance with the epithelial-like morphology reported above for the treated cells, which were more cubic than the elongated control cells. However, we observed from this DEG analysis and from qPCR analysis that one abundant mesenchymal marker, col1a1a, remained highly expressed in treated cells and was not differentially expressed between the two conditions (relative expression 125.0 ± 52.8 in treated cells, n=7 ; 123.1 ± 25.8 in control cells, n=8). This suggests that MET would be initiated but not terminated in our culture conditions.
Wnt signaling
The second signaling pathway whose terms were enriched in the GO analysis is the Wnt signaling pathway, and particularly the canonical one (Fig. 3). Because β-catenin is a key effector of Wnt signaling, the canonical pathway is referred to as Wnt/β-catenin signaling. The transduction of Wnt signal requires Wnt-induced activation of the receptors complex made of Frizzled (fzd) receptor and low-density lipoprotein co-receptor related 5 or 6 (LRP5/6). In other words, binding of the Wnt ligand to both receptors creates and activates the receptors complex. This initiates a series of molecular events that will protect cytosolic β-catenin from degradation. After nuclear import, β-catenin subsequently triggers the transcription of Wnt target genes by binding to transcription factors belonging to the T-cell factor/Lymphoid enhancer factor (Tcf/Lef) family [53,54].
Our gene to gene analysis of these Wnt-related actors revealed a strong deregulation of the Wnt/β-catenin signaling pathway in egg-extract treated cells (Table 3). Up-regulation of Wnt effectors combined with down- and up-regulation of inhibitors prevented the identification of a straightforward status for the Wnt signaling, be it an activated “on” or inhibited “off” status. In favor of an “on” status is the fact that some secreted Wnt ligands and fzd receptors were up-regulated in treated cells, the expression of the fzd10 receptor being especially strong. Moreover, extracellular Wnt agonists R-spondins (rspo2, rspo3), known to increase fzd receptors availability on the cell surface [53] and to stabilize the LRP5/6 co-receptors [54], were up regulated in treated cells. Additionally, down regulation of the extracellular inhibitors sfrp1a, sfrp2 and dkk1a [53] should be inducing a better availability of the Wnt ligand for fzd receptors.
However, other observations would be more in favor of an “off” status of the Wnt signaling. Firstly, the co-receptors LRP5/6 expression was not changed by the treatment. Despite increase in Wnt ligands and fzd receptors gene expression, LRP5/6 stability would stoichiometrically hamper the formation of the ternary proteic complex Wnt ligand / fzd receptor / LRP5/6 co-receptors that is essential for signal transduction. Secondly, many extracellular inhibitors upstream of the signaling pathway were upregulated in treated cells (Table 3). These included (i) notum1a and frzb, known to prevent Wnt ligand from binding to fzd receptor [54,55], (ii) sclerostin (sost) and dkk1b, which are blocking Wnt-fzd-LRP5/6 complex formation by interacting with LRP5/6 [53] and, (iii) kremen1, a membrane receptor which interacts with dkk1 to increase the removal of the LRP5/6 co-receptors from the cell surface by endocytosis [56]. Finally, the last point concerns the Tcf/Lef transcription factors, known as Tcf1, Tcf3, Tcf4 and Lef1 in mammals, which control the Wnt signaling activity through transcription of the target genes. In this study, Tcf7 (orthologue of Tcf1 in mice) expression was upregulated in treated cells (Table 3). It is known that the action of Tcf1 is triggered by nuclear beta catenin levels (Grainger 2019). Tcf1 acts as a transcriptional activator of Wnt target genes in presence of beta catenin. Conversely, Tcf1 acts as a transcriptional repressor of Wnt signaling in absence of beta catenin (Grainger 2019). In our study, treated cells showed a considerable collapse of fn1 and to a lesser extent a downregulation of pak1, also known as p21 (Table 3). The downregulation of these Wnt target genes associated with Tcf7 upregulation suggest that Tcf7 would act rather as a repressor in treated cells due to low beta catenin levels. All these observations reinforce the hypothesis that the deregulation of the Wnt signaling in treated cells would be rather in “off” configuration.
Altered cell adhesion of the treated cells and link with the deregulation of TGFβ and Wnt signaling pathways
In addition to the change in treated cells morphology, we also reported above a change in their behavior in culture. The cells showed a highly reduced ability to adhere throughout the culture process, and this could be due to changes in some gene expression. And indeed, from the GO analysis of all DEGs between treated and control cells, one biological process GO term highlighted the cell adhesion process (GO: 0007155; P value=3.2560E-08; FDR=1.17E-05). Furthermore, fibronectin (fn1) is a major protein of the extracellular matrix and it provides highly adhesive capacity to the cells by interaction with integrin transmembrane receptors [57]. As shown above, this actor was one of the most highly downregulated gene in our conditions, and this reduced expression is known to be deeply interconnected with the observed inhibition of the TGFβ and Wnt signaling pathways (Table 2, Table 3) [58]. In all, the poor adhesion capacity of our treated cells is another indication of the undergoing reprogramming event triggered by changes in the TGFβ and Wnt signaling pathways expression.
Some pluripotency markers remained silent in the treated cells
The process of somatic reprogramming in iPSCs is generally encompassing two phases [59]: (i) an early or initiation phase during which the somatic cells undergo a MET, lose their mesenchymal characteristics and develop an epithelial phenotype and, (ii) a late maturation phase allowing the reactivation of the pluripotency network. In order to characterize further the reprogramming extent of the treated fin cells, we focused on some marker genes related to pluripotency, previously characterized in goldfish during early development: pou2 (pou5f3 in zebrafish, oct4 in mammals), nanog, sox2 and c-myc [33,34,60]. We observed that none of these genes were identified among the DEGs, and their expression levels remained undetectable on the microarray. These observations were confirmed by qPCR that showed that pou2, nanog, sox2 and c-myc expression was below detection threshold in both treated and control cells.
It was shown previously in goldfish that nanog and pou2 silenced status in fin cells is associated with the hypermethylation of a CpGs locus in their promoter region [33,34]. We also showed recently that after nuclear transfer with non-treated fin cells, these loci underwent a partial and stochastic demethylation in the developing clones [16]. This indicated that embryonic reprogramming relaxed the DNA methylation status of these marker regions to some extent. The methylation profile of nanog and pou2 promoter regions in our treated cells was therefore analyzed, to assess whether some DNA demethylation took place at these marker sites after xenopus egg treatment. This would indeed be a necessary step ahead of any transcription reenabling of these 2 genes. Analysis of the CpG sites in pou2 and nanog promoter regions revealed that they did not underwent any significant demethylation in treated cells (Supplementary Fig. S3). Although the methylation of some CpG sites was lower in treated cells compared to controls, there was no significant differences in the overall DNA methylation rate of pou2 and nanog promoter regions. This indicates that no significant remodeling of the DNA methylation of pou2 and nanog pluripotency genes was triggered following egg-extract treatment, even if some minor variations were observed. In all, the silenced status of pluripotency marker genes associated with the absence of significant DNA methylation remodeling support the idea that treated cells would have been only partially reprogrammed by Xenopus egg-extract treatment. Our cells would not have reached the maturation phase of reprogramming characterized by Oct4 or Nanog and Sox2 re-expression as observed in mammalian somatic cells.
Alteration of de novo lipid biosynthesis in response to egg-extract treatment
Regarding the cluster of downregulated genes, the GO biological processes the most significantly affected by egg-extract are related to lipid metabolism (Fig 4A). Child GO terms targeted biosynthesis of steroid including cholesterol, and biosynthesis of unsaturated fatty acid. This was consistent with KEGG analysis showing the enrichment of the biosynthesis pathways of steroids, unsaturated fatty acids as well as the pathway of fatty acid metabolism (Fig 4B). In this process, acetyl-CoA represents the main precursor for de novo lipid biosynthesis. Produced in the mitochondria after glycolysis, acetyl-coA has to be metabolized into citrate so that it can exit the mitochondria. Once in the cytoplasm, citrate is then converted into lipogenic acetyl-CoA (see the molecular actors of lipogenesis in [61]). A detailed analysis of lipid metabolism genes showed a downregulation of several genes involved in the cytosolic synthesis of acetyl-CoA i.e. slc25a1b, a key mitochondrial transporter of citrate, aclya, which converts cytoplasmic citrate to acetyl-CoA, acss2, which produces acetyl-CoA from acetate, and the acyl transferase acat2 (Table 4). Lipid biosynthesis is also controlled by srebf1/2 transcription factors, whose expression was downregulated in our treated cells. The target genes of these transcription factors were downregulated as well. These included key enzymes for biosynthesis of cholesterol (hmgcs1, hmgcra1, msmo1, fdft1, cyp51, dhcr7) and fatty acid (fasn, sdc, elov1a, elov2, elov5, elov6) (Table 4). Overall, our results clearly indicate that the treated cells have strongly reduced their de novo lipid biosynthesis compared to control cells.