- 1Department of Genetics, Microbiology and Statistics, Faculty of Biology, University of Barcelona, Barcelona, Spain
- 2Institute of Biomedicine of the University of Barcelona (IBUB), Barcelona, Spain
Forkhead box (Fox) genes belong to the “winged helix” transcription factor superfamily. The function of some Fox genes is well known, such as the role of foxO in controlling metabolism and longevity and foxA in controlling differentiation of endodermal tissues. However, the role of some Fox factors is not yet well characterized. Such is the case of FoxK genes, which are mainly studied in mammals and have been implicated in diverse processes including cell proliferation, tissue differentiation and carcinogenesis. Planarians are free-living flatworms, whose importance in biomedical research lies in their regeneration capacity. Planarians possess a wide population of pluripotent adult stem cells, called neoblasts, which allow them to regenerate any body part after injury. In a recent study, we identified three foxK paralogs in the genome of Schmidtea mediterranea. In this study, we demonstrate that foxK1 inhibition prevents regeneration of the ectodermal tissues, including the nervous system and the epidermis. These results correlate with foxK1 expression in neoblasts and in neural progenitors. Although the triggering of wound genes expression, polarity reestablishment and proliferation was not affected after foxK1 silencing, the apoptotic response was decreased. Altogether, these results suggest that foxK1 would be required for differentiation and maintenance of ectodermal tissues.
Introduction
Forkhead box (Fox) genes belong to the “winged helix” superfamily of transcription factors, showing a Forkhead DNA-binding domain. Over 2,000 Fox proteins, phylogenetically classified into 25 families (A to S), have been identified in a number of species of fungi and metazoans (Benayoun, Caburet, and Veitia 2011; Pascual-Carreras et al., 2021). Fox genes control essential processes such as cell death, cell cycle and cell differentiation during all stages of development and in adult tissues, although many members show a tissue/stage-specific expression and function (Benayoun, Caburet, and Veitia 2011). The function and the implication of some Fox genes in human diseases is well known, such as the role in the differentiation of endodermal tissues of FoxA factors, the control of the cell cycle by FoxM, the role of FoxO in regulating metabolism and longevity and the contribution to speech acquisition of FoxP (Golson and Kaestner 2016). Some Fox transcription factors, especially FoxA, FoxO and FoxI members, are shown to act as pioneer factors, which open compacted chromatin to facilitate the binding of other transcription factors at enhancer sites (Lalmansingh, 2012). However, the function of some Fox families has not been properly addressed. An example are the FoxK genes, which are characterized by containing a forkhead-associated domain (FHA) besides the Forkhead DNA-binding domain, which mediates its interaction with multiple proteins (Ji et al., 2021). Few studies have been reported on the function of foxK genes, and they are mainly performed in mammals, which show two members (FoxK1 and FoxK2). From these studies, it is known that FoxK proteins are not tissue-specific but are ubiquitously expressed, and they regulate diverse biological processes, including cell proliferation, myogenic differentiation, cell cycle, autophagy, DNA damage and carcinogenesis (Liu et al., 2019).
Planarians are Lophotrochozoans well known for possessing pluripotent adult stem cells (called neoblasts) that confer them an extreme plasticity. They can regenerate any body part and change their size according to food availability (Wagner, Wang, and Reddien 2011; Baguñà 2012; Cebrià, Adell, and Saló 2018; Molina and Cebrià 2021). Previous studies of Fox genes function in planarians have demonstrated their evolutionary conserved role, and provided new data about their cellular function (Scimone, Lapan, and Reddien 2014; Vogg et al., 2014; He et al., 2017; Scimone et al., 2018). For instance, foxA was found to be essential for the maintenance of the pharynx and endodermal tissues (Koinuma et al., 2000; Adler et al., 2014) and foxO showed a conserved role in regulating metabolism during planarian regeneration and homeostasis (Pascual-Carreras et al., 2021). In a recent study, we identified 3 foxK paralogs in the genome of the planarian species S. mediterranea (Smed foxK-1-2-3). In this study, we aimed to analyze the function of Smed-Fox genes to provide new data about the role of this family in a model that, in contrast to mammals, displays stem cell-based tissue maintenance and regeneration.
We found that foxK1, expressed in neoblasts and their progenitors, is required for regeneration of neural and epidermal tissues. We hypothesize that foxK1 may play a role in activating and maintaining lineage specific enhancers.
Results
foxK1 RNAi Planarians Can Properly Regenerate Neither the Eyes nor the Nervous System
In a previous study, we identified 3 paralogs within the FoxK family in S. mediterranea (foxK1-2-3) (Pascual-Carreras et al., 2021). In this study, we have inhibited the three foxK genes by RNAi, showing that RNAi of foxk1 generates the strongest phenotype (Figure 1A) (Supplementary Figure S1A). foxk1 (RNAi) animals cannot regenerate proper anterior or posterior structures and 25% of them die after developing epidermal lesions after a single round of injection and amputation (Figure 1A). Prolonged silencing of foxK1 eventually impaired the viability of all the treated animals.
FIGURE 1. foxK1 silencing impairs planarian regeneration and eye differentiation. (A) foxK1 RNAi animals present smaller anterior and posterior blastemas in comparison to controls. Red arrowhead points to the epidermal injuries observed in foxK1 RNAi animals. Scale bar: 500 µm. (B–D) Aberrant eye regeneration after silencing foxK1. (B) foxK1 RNAi animals present smaller eyes that differentiate closer to the pre-existing tissue in comparison to controls (yellow arrows). Visualization of the eyes with an anti-VC1 immunostaining reveals defects in eye regeneration and an incorrect guidance of visual axons (red arrows). Scale bar: 100 µm. (C) Whole-mount in situ hybridization analysis of the expression of ovo in foxK1 RNAi animals. Scale bar: 200 µm. Quantification of ovo + cells/mm2 shows a significant decrease in the number of eye progenitor cells in foxK1 RNAi animals (n = 3) compared to gfp RNAi animals (n = 6). Values represent the mean ± standard deviation [s.d]. **p-value < 0.01. (D) Quantification of the eye distance/head tip distance ratio shows a significant increase in foxK1 RNAi animals (n = 14) compared to control gfp (RNAi) animals (n = 17). **p-value < 0.01. (A–D): All animals shown and analyzed are at 7 days of regeneration after anterior-posterior amputation. In the schematic planarians, red lines indicate the amputation planes and black squares the region analysed.
FoxK1 RNAi planarians regenerated smaller blastemas and eyes (Figure 1A). In some animals, eyes were merged in the midline (Figure 1B). Immunostaining of photoreceptors demonstrated that visual axons were aberrant and not properly connected to form a correct optic chiasm (Figure 1B). The expression of ovo, a pan-eye marker (Lapan and Reddien 2011), corroborated the decrease in the number of both differentiated and progenitor eye cells (Figure 1C). It also showed that eyes were more separated and closer to the anterior tip in foxK1-silenced planarians (Figure 1D). A phototaxis assay showed that foxK1 RNAi animals do not respond to the light stimulus, suggesting that although they could regenerate some photoreceptor cells (Figure 1B), these were not fully functional (Supplementary Figure S1B).
Analysis of the nervous system of foxK1 RNAi animals through anti-synapsin and anti-α-tubulin immunostaining showed that cephalic ganglia did not regenerate, and only a few axons at the level of the neural nerve cords (VNCs) crossed the midline (Figure 2A). The analysis of the expression of markers of differentiated neural cells such as gpas (Cebrià, Nakazawa, and Mineta 2002; Iglesias et al., 2011), pc2 (Collins et al., 2010) and th (Nishimura et al., 2007; Fraguas et al., 2011) corroborated that cephalic ganglia could not be formed (Figure 2A). The expression of sim, a marker of differentiated and progenitor neural cells (Cowles et al., 2013), was decreased in foxK1 RNAi animals (Supplementary Figure S2A). Quantification of differentiated (sim+) and progenitor (sim+/PIWI1+) neural cells indicated that both populations were decreased in foxK1 RNAi planarians (Figure 2B). The problem in regenerating the nervous system was not anterior specific, but it also affected the regeneration of proper posterior VNCs in the tail (Figure 2C).
FIGURE 2. foxK1 silenced animals show a reduced number of neural progenitors and almost absent differentiated neural structures. (A) Whole-mount immunostaining and in situ hybridization analyses (WISH) of the anterior nervous system regeneration in foxK1 RNAi animals 7 days post amputation. From top to bottom: immunostaining of the brain ganglia visualized with anti-synapsin which labels the brain synapsis; Immunostaining of the brain ganglia visualized with anti-α-tubulin which labels axons; WISH of gpas, pc2 and th which label differentiated neural cells. Scale bar: 100 µm. (B) Whole-mount fluorescence in situ hybridization analysis (WFISH) combined with immunostaining of the anterior blastemas in foxK1 RNAi animals 7 days post amputation. From left to right: WFISH with sim, which labels both neural progenitors and differentiated neurons; immunostaining with anti-PIWI1 which labels stem cells; double sim+/PIWI1+ positive cells, which indicate neural progenitors; Scale bar: 100 µm. Quantification of sim + PIWI + cells/mm2 reveals a decreased number of neural progenitors in foxK1 RNAi animals. Values represent the mean ± standard deviation [s.d] of a mean of 5 samples per condition. ***p-value < 0.001. (C) Whole-mount immunostaining analysis of the posterior nervous system regeneration in foxK1 RNAi animals at 7 days post amputation. From top to bottom: anti-synapsin labelling the synapsis and anti-α-tubulin labeling the axons, showing defects in the posterior VNCs. Scale bar: 100 µm. In the schematic planarians, red lines indicate the amputation planes and black squares the region analysed.
According to the RNAi phenotype, foxK1 was expressed in the nervous system of intact animals and during regeneration (Supplementary Figure S1C). Furthermore, the analysis of the expression pattern of foxK1 in irradiated animals showed that it was expressed in differentiated cells and also in progenitors (Supplementary Figure S1C), which agrees with published single-cell RNAseq databases (Plass et al., 2018; Zeng et al., 2018) (Supplementary Figures S1D–F).
foxK1 RNAi Animals Cannot Regenerate a Proper Epidermis
foxK1 appears expressed in the neural and in the epidermal lineages according to the published single-cell RNAseq data (Supplementary Figure S1C). Furthermore, its inhibition produced epidermal lesions (Figure 1A). Thus, we analyzed deeper the integrity of the epidermal tissue in RNAi animals. Quantification of the nucleus demonstrated a decrease of cell density in the epidermis (Figure 3A). Interestingly, the nucleus appeared bigger in foxK1 RNAi animals, as demonstrated by the quantification of their area (Figure 3A). Furthermore, the cilia of the epidermal cells in the dorsal midline were not properly regenerated (Figure 3B) and, in the ventral part, appeared disorganized (Figure 3B). Immunolabeling of adherens junctions with anti-β-catenin2 antibody indicated that these structures were not properly assembled in foxK1 silenced animals, since the signal was dimmer and the pattern was irregular (Figure 3B).
FIGURE 3. foxK1 silenced planarians show defects in the epidermis. (A) TO-PRO®-3 staining of the epidermis in anterior wounds and quantification showing that the epidermal nuclear density decreases (number of nuclei/mm2 of the area of the confocal section) while the size for each nucleus increases (average size in mm2 of each nucleus in the confocal section of every individual) after foxK1 silencing. Animals at 7 days of regeneration are shown. All images are dorsal views. Values represent the mean ± standard deviation [s.d] of a mean of at least 11 planarians per condition. Scale bar: 50 µm ****p-value < 0.0001. (B) Immunostaining of the epidermis in anterior wounds. Anti-α-tubulin shows the aberrant regeneration of cilia along the dorsal midline and the ventral epidermis. Merge of the immunostaining of anti-β-catenin2 and TO-PRO®-3 in the dorsal region between the eyes showing a misorganization of the epithelial adherens junctions. Animals at 7 days of regeneration are shown. Scale bar: 50 µm. In the schematic planarians, red lines indicate the amputation planes and black squares the region analysed.
On the other hand, other tissues such as the digestive or the excretory systems seemed not to be affected after foxK1 inhibition (Supplementary Figures S2B–D).
Overall, these results suggest that foxK1 is required for maintenance and regeneration of the epidermis in planarians.
Apoptotic Response is Affected in foxK1 RNAi Animals
Regeneration of the missing tissues in planarians is preceded by the triggering of early signals that activate the expression of the early genes as well as the proliferative and apoptotic response (Saló and Baguñà 1984; Wenemoser and Reddien 2010; Wenemoser et al., 2012; Pirotte et al., 2015; Owlarn et al., 2017; Jaenen et al., 2021). Analysis of the expression of the early response gene runt1, as well as the polarity genes notum and wnt1, whose early expression is independent of polarity and occurs after any injury (Wenemoser et al., 2012; Sandmann et al., 2011; Petersen and Reddien 2009, 2011), indicated that the early gene response was not affected after foxK1 inhibition (Figure 4A).
FIGURE 4. Apoptotic response is affected in foxK1 RNAi animals, but the polarity establishment and proliferative response occurs normally. (A) Whole-mount in situ hybridization analysis of the expression of the early wound response genes in controls and foxK1 RNAi animals 6 h post amputation (runt1 and notum) and 12 h post amputation (wnt1). From top to bottom: runt1 in anterior wounds, notum in anterior wounds and wnt1 in posterior wounds. Scale bar: 100 µm. (B) Quantification of mitotic PH3+ immunolabeled cells after silencing foxK1 at several time points after anterior-posterior amputation. h, hours; d, days (*p-value < 0.05; ** p-value < 0.01 ****p-value < 0.0001. Student’s t-test). Values represent the mean ± standard error of the mean [s.e.m] of a mean of at least 4 samples per time point and condition. (C) Quantification of caspase-3 activity after silencing foxK1 at several time points after anterior-posterior amputation. h, hours; d, days (*p-value < 0.05; Student’s t-test). Values represent the mean ± standard error of the mean [s.e.m] of a mean of 3 biological samples per time point and condition. (D) Whole-mount in situ hybridization analysis of the expression of the genes responsible for establishing the anterior and posterior identity in controls and foxK1 RNAi animals 3 days post amputation (notum, wnt1, frz4, post2d and hox4b), 6 days post amputation (wnt11_1 and wnt11_2) and 7 days post amputation (sfrp1). Anterior polarity genes (from top to bottom): notum and srfp1. Posterior polarity genes (from top to bottom): wnt1, frz4, post-2d, hox4b, wnt11-1 and wnt11-2. Scale bar: 100 µm. In the schematic planarians, red lines indicate the amputation planes and black squares the region analysed.
In addition, quantification of the proliferative response near the anterior and the posterior wounds showed that the dynamics of proliferation was not perturbed in foxK1 RNAi animals, and the two mitotic peaks at 6 and 48 h described during planarian regeneration occurred normally (Figure 4B). It must be noted, however, that proliferation appeared to be increased in regenerating animals before amputation, but this did not interfere with a proper proliferative response after injury. On the other side, quantification of the apoptotic response showed a decrease in foxK1 RNAi animals, both at early (4 h) and late phases of regeneration (7 days) (Figure 4C).
These results indicate the only event found deregulated in foxK1 RNAi animals during the early wound response is the activation of apoptosis.
foxK1 RNAi Animals Show a Correct Polarity Establishment
During planarian regeneration, the general early wound response is followed by a pole-specific response, which includes the expression of genes specific to anterior or posterior wounds that are required to specify the identity and to differentiate the proper missing tissue (Wurtzel et al., 2015; Cebrià, Adell, and Saló 2018). To clarify whether the problems of foxK1 RNAi planarian regeneration could be associated with an incorrect pole-specific response, we analyzed several genes specific to anterior and posterior regeneration. The results showed that both anterior and posterior poles were properly formed in foxK1 silenced planarians, since notum and wnt1 were restricted to anterior and posterior wounds, respectively (Figure 4D) (Petersen and Reddien 2011; Adell et al., 2009). post2-d and wnt11-2, two positional control genes (PCG) (Witchley et al., 2013) that are specifically expressed in posterior facing wounds after wnt1 (Adell et al., 2009; Currie et al., 2016; Tewari et al., 2019), also showed a proper expression in foxK1 RNAi planarians (Figure 4D). However, the expression of some anterior and posterior PCG appeared decreased in foxK1 RNAi animals, such as sfrp1 in anterior; and fz4, hox4b and wnt11-1 in posterior (Gurley, Rink, and Alvarado 2008; Petersen and Reddien 2008; Currie et al., 2016).
These results suggest that polarity was not affected in foxK1 silenced planarians, even though the inability to regenerate a correct head and tail was associated with a decreased expression of some markers of anterior and posterior polarity.
Discussion
In this study, we have found that inhibition of foxK1 impedes correct regeneration of missing structures, while inhibition of foxK2 or foxK3 only produces mild defects. The strong phenotype observed after foxK1, but not foxK2 or foxK3 inhibition, could be due to the fact that the last two are mainly expressed in differentiated cells (Pascual-Carreras et al., 2021), while foxK1 is also expressed in progenitors and in neoblasts. foxk1 RNAi animals cannot properly regenerate the nervous system or the epidermis, which correlates with a decrease in the number of differentiated neural and epidermal cells. At least in the case of the nervous system, foxK1 RNAi animals also show a decrease in the number of progenitors (sim+/PIWI1+) cells. However, the number of cycling cells is not decreased. This result could indicate that foxK1 is required for maintenance and differentiation of the neural precursors and may be for specification of neoblasts towards the neural fate. This result agrees with a recent study in human ESCs, in which the FoxK2 transcription factor is found to be bound to thousands of regulatory regions and to act as a pioneer factor, playing a role in enhancer priming during differentiation (Ji et al., 2021). Thus, foxK1 could also act as a pioneer factor during neural determination in planarians. This pioneer role has been already demonstrated for other Fox families, such as FoxA, FoxD or FoxH (Krishnakumar et al., 2016; Charney et al., 2017). In human ESCs, FoxK2 binding to enhancers is maintained as cells differentiate to neural precursor cell types (Ji et al., 2021). This could also be the case of foxK1 in planarians. In contrast to other Fox factors, foxK1 could be recruited not only in enhancer regions prior to their activation, but its binding could be retained as enhancers become activated during differentiation, since both progenitor and differentiated cells are decreased in RNAi animals.
Although the cephalic ganglia are more affected than VNCs in foxK1 RNAi planarians, foxK1 seems to be required for nervous system differentiation in general, since posterior fascicles of the VNCs appear thinner and the decrease in the expression of sim affects to all the body regions.
The defect in the eyes of foxK1 RNAi animals is one of the first traits that can be observed, since they are the only structures clearly identified when planarians regenerate the head. With the analysis of specific markers, we demonstrate that eyes are not properly specified nor differentiated, since both eye precursor and eye differentiated cells are decreased. However, we cannot clarify whether this defect is due to an autonomous role of foxK1 in the eye cells or it is linked to the inability of these animals to regenerate the cephalic ganglia, since the eyes connect to the brain while they both are regenerating. We have not observed the expression of foxK1 in eye cells, but we cannot discard it.
FoxK1 RNAi animals die after developing epidermal lesions, thus it can be due to its role in the maintenance of the epidermal layer. We observed a decrease in epidermal cell density and an increase in the size of the nucleus, which agrees with the role of foxK1 for correct terminal differentiation of epidermal cells. However, we cannot conclude whether foxK1 is also required for maintenance and differentiation of epidermal progenitors, and whether it also acts as a pioneer factor in the epidermal lineage, since the number of epidermal precursors has not been analyzed. It is tempting to speculate that the increase in the size of the nucleus could be related with a decompaction of the chromatin, which would agree with the role of foxK1 as a pioneer factor, although it could also be due to other mechanisms.
FoxK1 expression is not neural or epidermal specific, thus it could be also required for specification or differentiation of other tissues. However, its RNAi inhibition apparently only affects the ectodermal lineages, as the digestive or the excretory system appear normal.
Despite in other models, as Drosophila, foxK has been involved in the control of proliferation of myogenic stem cells (Garry et al., 2000), our results indicate that in a stem cell-based model as planarians the FoxK family is not required for neoblast proliferation. foxK2 and foxK3 are not expressed in neoblasts, and foxK1 is expressed in neoblasts but we have found that the bimodal proliferative response triggered by injury is correct in foxK1 RNAi planarians, and in animals that have already regenerated there are more proliferative cells. This result agrees with the hypothesis that foxK1 is activating neoblast and progenitor enhancers required for determination and differentiation, but it is not required for their proliferation. If foxK1 is required to prime neoblasts, they would maintain their naive and cycling properties when foxK1 is not functional, which would explain the increase in proliferation seen in regenerated animals that must maintain their tissues. However, the triggering of stem cell proliferation after injury is not dependent on foxK1.
foxK1 RNAi animals present a correct early mitotic response and early wound gene expression. However, we have found that inhibition of foxK1 decreases the early apoptotic response at 4 h of regeneration, suggesting that this gene could have a key role in controlling apoptosis. In this study we have not been able to clarify the impact of the early apoptotic decrease in the phenotype observed. However, the fact that the apoptotic activity seems also reduced at later stages of regeneration in foxK1 RNAi planarians suggests that the observed decrease in the number of neural progenitors and epidermal cells is not due to an increase in their cell death.
Polarity establishment is neither affected in foxK1 RNAi planarians. Some PCGs are not correctly expressed, such as sfrp1 in anterior and fz4, hox4b and wnt11-1 in posterior. However, this could be associated with the inability to form a proper head or tail, but not to the incorrect pole formation, since notum and wnt1 are correctly expressed at all regeneration stages. A possibility is that the expression of those PCGs genes could require the differentiation of a nervous system.
While several Fox proteins have been implicated in enhancer activation in a tissue-specific manner, the role of the FoxK has been associated to a range of biological processes apparently unconnected. This could be mainly due to its ubiquitous expression and to having an additional conserved domain, the FHA domain, which binds to multiple proteins. Analysing foxK1 function in planarians we have been able to characterize a specific role of foxK1 in the differentiation of the neural tissue, possibly by acting as a pioneer factor, as other members of the Fox family, which allows and maintains the activation of enhancers during differentiation of neural progenitors. Furthermore, the role of foxK1 in controlling apoptosis deserves further attention, since a functional relationship between FoxK and apoptosis has only been suggested in cancer cells (Liu et al., 2019).
Materials and Methods
Animal Maintenance
S. mediterranea from the asexual clonal line BCN-10 were used for all experiments. Animals were maintained at 20°C in 1X PAM water as previously described (Cebrià et al., 2007). Animals were fed twice per week with organic veal liver and were starved for at least 1 week before experiments. For irradiation experiments, planarians were exposed to 86 Gray (Gy) of γ-irradiation.
RNA Interference
Double-stranded RNA (dsRNA) for Smed-foxK genes were synthesised as previously described (Sánchez Alvarado and Newmark 1999). The region amplified by the foxK1 specific forward primer 5′ AATCTCATTCTTTTATTCCT 3′ and reverse primers 5′ AGTATCGTTCATCAGTCCAT 3′ or 5′ TTGAGCGTATGAATATGGAGG 3’ was used to synthesize the dsRNA. The injection protocol consisted of 1 round of 4 consecutive days of injections or 2 rounds of 3 consecutive days of injections separated by a 4-days interval. A Nanoject II (Drummond Scientific, Broomall, PA, United States) was used to administer 3 injections of 32 nl of dsRNA (1 μg/μl) per day. Control animals were injected with gfp dsRNA. In each round, 1 day after the last injection, planarians were amputated to induce anterior and posterior regeneration.
Single-Cell Sequencing Data
FoxK1,2,3 genes expression profile (dd_Smed_v6_4500, dd_Smed_v6_5767 and dd_Smed_v6_7583, respectively) were obtained from the planaria single-cell database hosted by the Rajewsky lab at the Berlin Institute for Medical Systems Biology of the Max Delbrück Center, Berlin (Plass et al., 2018) and the planaria single-cell database hosted by the Sánchez Alvarado lab at Planosphere website (Zeng et al., 2018).
Whole-Mount In Situ Hybridisation and Whole-Mount Fluorescent In Situ Hybridisation
Whole-mount in situ hybridisation (WISH) (Currie et al., 2016) whole-mount fluorescent in situ hybridisation (WFISH) (King and Newmark 2013) were performed as previously described. Riboprobes for in situ hybridisation were synthesised using the DIG RNA labe ling kit (Sp6/T7, Roche) following the manufacturer’s instructions. Samples were mounted in 70% glycerol/PBS solution.
Immunohistochemistry
Whole-mount immunohistochemistry was performed as previously described (Ross et al., 2015). Treated animals were euthanised by immersion in cold 2% HCl in ultrapure H2O for 5 min, washed with PBS-Tx (PBS +0.3% Triton X-100) at room temperature (RT), and placed in a fixative solution (4% formaldehyde in PBS-Tx) for 15 min at RT with agitation. Subsequently, samples were washed with PBS-Tx and bleached in 6% H2O2 (in PBS-Tx) at RT for 16 h under direct light. Bleached animals were then washed with PBS-Tx, incubated for 2 h in 1% blocking solution (1% BSA in PBS-Tx), and overnight at 4°C in the primary antibody (diluted in blocking solution). The following primary antibodies were used: anti-phospho-histone3 (PH3, Cell Signaling Technology) to detect mitotic cells, diluted 1:300; anti-SYNAPSIN, a pan-neural marker (anti-SYNORF1, Developmental Studies Hybridoma Bank), diluted 1:50; anti-VC-1, specific for planarian photosensitive cells (Sakai et al., 2000), diluted 1:15,000; anti-SMEDWI-1, specific for neoblasts, diluted 1:1,500 (Rouhana et al., 2010) (März, Seebeck, and Bartscherer 2013); polyclonal antibody against Smed-β-catenin-2 used to visualize adherens junctions (Chai et al., 2010), and used at a 1:2,000 dilution; the AA4.3 antibody (Developmental Studies Hybridoma Bank), specific for α-tubulin to visualize the epithelial cilia, diluted 1:20. The following secondary antibodies were used: Alexa-488-conjugated goat anti-mouse (Molecular Probes), diluted 1:400; and Alexa 568-conjugated goat anti-rabbit (Molecular Probes), diluted 1:1,000. Samples were mounted in 70% glycerol/PBS solution. Nuclei were stained with DAPI (1:5000; Sigma-Aldrich) and TO-PRO®-3 (1:3,000, Thermo Fisher Scientific, Waltham, MA, United States).
Phototactic Assay
Phototactic assay was carried out using a modified version of the method described by (Inoue et al., 2004). Planarian behavior was recorded during 120 s using an overhead digital video camera (Canon EOS550D). The behavior analysis software SMART v.2.5.21 (Panlab, Spain) was used to quantify the time spent by the animals in each of the three virtual subdivisions of a transparent container of 60 × 30 × 10 mm, filled with planarian water. To obtain a light gradient, the container was protected by a black screen with a hole that allowed the entrance of 500 lux of white light from one side of the container.
Caspase-3 Activity Assay
The protein concentration of the cell lysates of three biological samples of 4–7 planarians each was measured using BioRad protein reagent. Fluorometric analysis of caspase-3 activity was performed as described previously (González-Estévez et al., 2007) using 16 μg of protein extract after incubation for 2 h at 37°C with the caspase-3 substrate Ac-DEVD-AMC. A Fluostar Optima microplate fluorescence reader (BMG Labtech) and a luminescence spectrophotometer (Perkin- Elmer LS- 50) were used to measure luminescence. The following settings were applied: excitation, 380 nm; emission, 440 nm. Three technical replicates were analysed per condition.
Microscopy, Image Acquisition, and Image Analysis
Live animals were photographed with an sCM EX-3 high end digital microscope camera (DC.3000s, Visual Inspection Technology). WISH, WFISH, and immunostained animals were observed with a Leica MZ16F stereomicroscope. Images were captured with the ProGres C3 camera (Jenoptik) and then processed in Photoshop CS6 for publication. Representative images of WFISH and immunostained animals were captured with confocal laser scanning microscopy (Leica TCS-SPE microscope) and processed in ImageJ2.3.0 and Photoshop CS6 for publication.
Statistical Analysis
All comparisons were performed using the unpaired Student’s t-test, after first confirming data normality and homogeneity using the Shapiro-Wilk test.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
PC-C, MDM, EP-C, FC, and TA contributed to the conception and design of the study. PC-C, MDM, GC, SF, and EP-C performed the experiments and statistical analysis. PC-C, MDM and TA wrote the first draft of the manuscript. ES, FC, and TA obtained funding. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
MDM has received funding from the postdoctoral fellowship programme Beatriu de Pinós, funded by the Secretary of Universities and Research (Government of Catalonia) and by the European Union Horizon 2020 research and innovation programme under Marie Sklodowska-Curie grant agreement No. 801370. EP-C was the recipient of an FPI (Formación del Profesorado Investigador) scholarship from the Spanish Ministerio de Educación y Ciencia (MEC). FC was supported by grants BFU 2015-65704-P and PGC 2018-100747-B-100 from the Ministerio de Ciencia, Innovación y Universidades, Spain. ES and TA received funding from the Ministerio de Educación y Ciencia (grant number BFU 2017-83755-P and PID2020-116372GB-I00). ES, FC and TA benefits from 2017SGR 1455 from AGAUR (Generalitat de Catalunya).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank all members of the ESE, FC and TA laboratories for discussions. We thank H. Orii and Prof. K. Watanabe for providing anti-VC-1. We thank Dra. K. Bartscherer for providing anti-PWI1. Monoclonal anti-SYNORF1 and AA4.3 antibodies were obtained from the Developmental Studies Hybridoma Bank, developed under the auspices of the National Institute of Child Health and Human Development and maintained by the Department of Biological Sciences, University of Iowa, Iowa City, IA, United States.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.808045/full#supplementary-material
Supplementary Figure S1 | RNAi silencing and expression pattern of the three foxK S. mediterranea paralogs. (A) foxK2 or foxK3 silenced planarians regenerate the missing head and tail. Scale bar: 500 µm. (B) Phototactic assay. Schematic representation of the position of treated planarians in each of the three subdivisions of 20 cm each after 120 s of light exposure. Red dots represent the position of foxK1 silenced animals and blue dots the position of gfp RNAi control planarians. Graphical representation of the quantification of the velocity (meter/seconds) of gfp and foxK1 silenced planarians showing the reduced mobility of foxK1 silenced planarians. (C) Whole-mount in situ hybridization of foxK1 in intact, γ-irradiated and planarians at 2 days of regeneration after anterior-posterior amputation. Scale bar: 200 µm. (D) Schematic representation of single-cell transcriptomic expression of foxK genes [according to (Plass et al., 2018)]. It must be noted that in a previous publication (Pascual-Carreras et al., 2021), these plots were not properly assigned, but this is corrected in the present publication. (E) Schematic representation of single-cell transcriptomic expression of foxK1 in all planarian neoblasts (Nb) subtypes [according to (Zeng et al., 2018)]. (F) Relative gene expression of foxK1 and piwi1 for the individual stages during differentiation of the epidermal and neuronal planarian cell lineages (in pseudotime) according to single cell analysis [https://shiny.mdc-berlin.de/psca/, (Plass et al., 2018)]. a.u, arbitrary units. In the schematic planarians, red lines indicate the amputation planes and black squares the region analysed.
Supplementary Figure S2 | Expression pattern of tissue-specific markers in foxK1 RNAi planarians. (A) Whole-mount in situ hybridization of the neural marker sim (Cowles et al., 2013) reveals a decreased number of sim-expressing cells all along the planarian body of foxK1 silenced planarians. Animals at 7 days of regeneration are shown. Scale bar: 100 µm. (B) Whole-mount in situ hybridization of innexin10 (Nogi and Levin 2005) labeling the planarian excretory system reveals correct regeneration of excretory cells in foxK1 silenced planarians. Animals at 7 days of regeneration are shown. Scale bar: 100 µm. (C) Whole-mount in situ hybridization of the gut marker pk (Fraguas et al., 2014) showing that treated planarians regenerate correctly the gut. Whole-mount immunostaining against α-tubulin showing a detail of the newly regenerated posterior gut branches. Animals at 7 days of regeneration are shown. Scale bars: 100 µm.
References
Adell, T., Salò, E., Boutros, M., and Bartscherer, K. (2009). Smed-Evi/Wntless Is Required for β-catenin-dependent And-independent Processes during Planarian Regeneration. Development 136 (6), 905–910. doi:10.1242/dev.033761
Adler, C. E., Seidel, C. W., McKinney, S. A., and Sánchez Alvarado, A. (2014). Selective Amputation of the Pharynx Identifies a FoxA-dependent Regeneration Program in Planaria. ELife 3, 1–22. doi:10.7554/elife.02238
Baguñà, J. (2012). The Planarian Neoblast: the Rambling History of its Origin and Some Current Black Boxes. Int. J. Dev. Biol. 56, 19–37. doi:10.1387/ijdb.113463jb
Benayoun, B. A., Caburet, S., and Veitia, R. A. (2011). Forkhead Transcription Factors: Key Players in Health and Disease. Trends Genet. 27 (6), 224–232. doi:10.1016/j.tig.2011.03.003
Cebrià, F., Adell, T., and Saló, E. (2018). Rebuilding a Planarian: From Early Signaling to Final Shape. Int. J. Dev. Biol. 62 (6–8), 537–550. doi:10.1387/ijdb.180042es
Cebrià, F., Guo, T., Jopek, J., and Newmark, P. A. (2007). Regeneration and Maintenance of the Planarian Midline Is Regulated by a Slit Orthologue. Develop. Biol. 307 (2), 394–406. doi:10.1016/j.ydbio.2007.05.006
Cebrià, F., Nakazawa, M., Mineta, K., Ikeo, K., Gojobori, T., and Agata, K. (2002). Dissecting Planarian Central Nervous System Regeneration by the Expression of Neural-specific Genes. Dev. Growth Differ. 44 (2), 135–146. doi:10.1046/j.1440-169x.2002.00629.x
Chai, G., Ma, C., Bao, K., Zheng, L., Wang, X., Sun, Z., et al. (2010). Complete Functional Segregation of Planarian β-Catenin-1 and -2 in Mediating Wnt Signaling and Cell Adhesion. J. Biol. Chem. 285 (31), 24120–24130. doi:10.1074/jbc.M110.113662
Charney, R. M., Forouzmand, E., Cho, J. S., Cheung, J., Paraiso, K. D., Yasuoka, Y., et al. (2017). Foxh1 Occupies Cis -Regulatory Modules Prior to Dynamic Transcription Factor Interactions Controlling the Mesendoderm Gene Program. Develop. Cel 40 (6), 595–607. e4. doi:10.1016/j.devcel.2017.02.017
Collins, J. J., Hou, X., Romanova, E. V., Lambrus, B. G., Miller, C. M., Saberi, A., et al. (2010). Genome-Wide Analyses Reveal a Role for Peptide Hormones in Planarian Germline Development. Plos Biol. 8 (10), e1000509. doi:10.1371/journal.pbio.1000509
Cowles, M. W., Brown, D. D. R., Nisperos, S. V., Stanley, B. N., Pearson, B. J., and Zayas, R. M. (2013). Genome-Wide Analysis of the BHLH Gene Family in Planarians Identifies Factors Required for Adult Neurogenesis and Neuronal Regeneration. Development (Cambridge) 140 (23), 4691–4702. doi:10.1242/dev.098616
Currie, K. W., Brown, D. D. R., Zhu, S., Xu, C., Voisin, V., Bader, G. D., et al. (2016). HOX Gene Complement and Expression in the Planarian Schmidtea Mediterranea. EvoDevo 7, 7. doi:10.1186/s13227-016-0044-8
Fraguas, S., Barberán, S., and Cebrià, F. (2011). EGFR Signaling Regulates Cell Proliferation, Differentiation and Morphogenesis during Planarian Regeneration and Homeostasis. Develop. Biol. 354 (1), 87–101. doi:10.1016/j.ydbio.2011.03.023
Fraguas, S., Barberán, S., Iglesias, M., Rodríguez-Esteban, G., and Cebrià, F. (2014). egr-4, a Target of EGFR Signaling, Is Required for the Formation of the Brain Primordia and Head Regeneration in planarians. Development (Cambridge, England) 141 (9), 1835–1847. doi:10.1242/dev.101345
Garry, D. J., Meeson, A., Elterman, J., Zhao, Y., Yang, P., Bassel-Duby, R., et al. (2000). Myogenic Stem Cell Function Is Impaired in Mice Lacking the Forkhead/Winged Helix Protein MNF. Proc. Natl. Acad. Sci. 97 (10), 5416–5421. doi:10.1073/pnas.100501197
Golson, M. L., and Kaestner, K. H. (2016). Fox Transcription Factors: From Development to Disease. Development 143 (24), 4558–4570. doi:10.1242/dev.112672
González-Estévez, C., Felix, D. A., Aboobaker, A. A., and Saló, E. (2007). Gtdap-1 Promotes Autophagy and Is Required for Planarian Remodeling during Regeneration and Starvation. Proc. Natl. Acad. Sci. 104 (33), 13373–13378. doi:10.1073/pnas.0703588104
Gurley, K. A., Rink, J. C., and Alvarado, A. S. (2008). β-Catenin Defines Head versus Tail Identity during Planarian Regeneration and Homeostasis. Science 319 (5861), 323–327. doi:10.1126/science.1150029
He, X., Lindsay-Mosher, N., Li, Y., Molinaro, A. M., Pellettieri, J., and Pearson, B. J. (2017). FOX and ETS Family Transcription Factors Regulate the Pigment Cell Lineage in Planarians. Development 144 (24), 4540–4551. doi:10.1242/dev.156349
Iglesias, M., Almuedo-Castillo, M., Aboobaker, A. A., and Saló, E. (2011). Early Planarian Brain Regeneration Is Independent of Blastema Polarity Mediated by the Wnt/β-Catenin Pathway. Develop. Biol. 358 (1), 68–78. doi:10.1016/j.ydbio.2011.07.013
Inoue, T., Kumamoto, H., Okamoto, K., Umesono, Y., Sakai, M., Alvarado, A. S., et al. (2004). Morphological and Functional Recovery of the Planarian Photosensing System during Head Regeneration. Zoolog. Sci. 21 (3), 275–283. doi:10.2108/zsj.21.275
Jaenen, V., Fraguas, S., Bijnens, K., Heleven, M., Artois, T., Romero, R., et al. (2021). Reactive Oxygen Species rescue Regeneration after Silencing the MAPK-ERK Signaling Pathway in Schmidtea Mediterranea. Sci. Rep. 11 (1), 1–16. doi:10.1038/s41598-020-79588-1
Ji, Z., Li, Y., Liu, S. X., and Sharrocks, A. D. (2021). The Forkhead Transcription Factor FOXK2 Premarks Lineage-specific Genes in Human Embryonic Stem Cells for Activation during Differentiation. Nucleic Acids Res. 49 (3), 1345–1363. doi:10.1093/nar/gkaa1281
King, R. S., and Newmark, P. A. (2013). In Situ Hybridization Protocol for Enhanced Detection of Gene Expression in the Planarian Schmidtea Mediterranea. BMC Develop. Biol. 13 (March), 8. doi:10.1186/1471-213X-13-8
Koinuma, S., Umesono, Y., Watanabe, K., and Agata, K. (2000). Planaria FoxA (HNF3) Homologue Is Specifically Expressed in the Pharynx-Forming Cells. Gene 259 (1–2), 171–176. doi:10.1016/S0378-1119(00)00426-1
Krishnakumar, R., Chen, A. F., Pantovich, M. G., Danial, M., Parchem, R. J., Labosky, P. A., et al. (2016). FOXD3 Regulates Pluripotent Stem Cell Potential by Simultaneously Initiating and Repressing Enhancer Activity. Cell Stem Cell 18 (1), 104–117. doi:10.1016/j.stem.2015.10.003
Lalmansingh, A. S., Karmakar, S., Jin, Y., and Nagaich, A. K. (2012). Multiple Modes of Chromatin Remodeling by Forkhead Box Proteins. Biochim. Biophys. Acta (Bba) - Gene Regul. Mech. 1819 (7), 707–715. doi:10.1016/j.bbagrm.2012.02.018
Lapan, S. W., and Reddien, P. W. (2011). Dlx and Sp6-9 Control Optic Cup Regeneration in a Prototypic Eye. Plos Genet. 7 (8), e1002226. doi:10.1371/journal.pgen.1002226
Liu, Y., Ding, W., Ge, H., Ponnusamy, M., Wang, Q., Hao, X., et al. (2019). FOXK Transcription Factors: Regulation and Critical Role in Cancer. Cancer Lett. 458 (38), 1–12. doi:10.1016/j.canlet.2019.05.030
März, M., Seebeck, F., and Bartscherer, K. (2013). A Pitx Transcription Factor Controls the Establishment and Maintenance of the Serotonergic Lineage in Planarians. Development (Cambridge, England) 140 (22), 4499–4509. doi:10.1242/dev.100081
Molina, M. D., and Cebrià, F. (2021). Decoding Stem Cells: An Overview on Planarian Stem Cell Heterogeneity and Lineage Progression. Biomolecules 11 (10), 1532–1623. doi:10.3390/biom11101532
Nishimura, K., Kitamura, Y., Inoue, T., Umesono, Y., Sano, S., Yoshimoto, K., et al. (2007). Reconstruction of Dopaminergic Neural Network and Locomotion Function in Planarian Regenerates. Devel Neurobio 67 (8), 1059–1078. doi:10.1002/dneu.20377
Nogi, T., and Levin, M. (2005). Characterization of Innexin Gene Expression and Functional Roles of Gap-Junctional Communication in Planarian Regeneration. Develop. Biol. 287 (2), 314–335. doi:10.1016/j.ydbio.2005.09.002
Owlarn, S., Klenner, F., Schmidt, D., Rabert, F., Tomasso, A., Reuter, H., et al. (2017). Generic Wound Signals Initiate Regeneration in Missing-Tissue Contexts. Nat. Commun. 8 (1), 1–13. doi:10.1038/s41467-017-02338-x
Pascual-Carreras, E., Herrera-Úbeda, C., Rosselló, M., Coronel-Córdoba, P., Garcia-Fernàndez, J., Saló, E., et al. (2021). Analysis of Fox Genes in Schmidtea Mediterranea Reveals New Families and a Conserved Role of Smed-FoxO in Controlling Cell Death. Sci. Rep. 11 (1), 1–18. doi:10.1038/s41598-020-80627-0
Petersen, C. P., and Reddien, P. W. (2009). A Wound-Induced Wnt Expression Program Controls Planarian Regeneration Polarity. Proc. Natl. Acad. Sci. 106 (40), 17061–17066. doi:10.1073/pnas.0906823106
Petersen, C. P., and Reddien, P. W. (2011). Polarized Notum Activation at Wounds Inhibits Wnt Function to Promote Planarian Head Regeneration. Science 332 (6031), 852–855. doi:10.1126/science.1202143
Petersen, C. P., and Reddien, P. W. (2008). Smed- β Catenin-1 Is Required for Anteroposterior Blastema Polarity in Planarian Regeneration. Science 319 (5861), 327–330. doi:10.1126/science.1149943
Pirotte, N., Stevens, A.-S., Fraguas, S., Plusquin, M., Van Roten, A., Van Belleghem, F., et al. (2015). Reactive Oxygen Species in Planarian Regeneration: An Upstream Necessity for Correct Patterning and Brain Formation. Oxidative Med. Cell Longevity 2015, 1–19. doi:10.1155/2015/392476
Plass, M., Solana, J., Wolf, F. A., Ayoub, S., Misios, A., Glažar, P., et al. (2018). Cell Type Atlas and Lineage Tree of a Whole Complex Animal by Single-Cell Transcriptomics. Science 360 (6391). doi:10.1126/science.aaq1723
Ross, K. G., Omuro, K. C., Taylor, M. R., Munday, R. K., Hubert, A., King, R. S., et al. (2015). Novel Monoclonal Antibodies to Study Tissue Regeneration in Planarians. BMC Dev. Biol. 15 (1). doi:10.1186/s12861-014-0050-9
Rouhana, L., Shibata, N., Nishimura, O., and Agata, K. (2010). Different Requirements for Conserved Post-Transcriptional Regulators in Planarian Regeneration and Stem Cell Maintenance. Develop. Biol. 341 (2), 429–443. doi:10.1016/j.ydbio.2010.02.037
Sakai, F., Agata, K., Orii, H., and Watanabe, K. (2000). Organization and Regeneration Ability of Spontaneous Supernumerary Eyes in Planarians -Eye Regeneration Field and Pathway Selection by Optic Nerves-. Zoolog. Sci. 17 (3), 375–381. doi:10.2108/jsz.17.375
Saló, E., and Baguñà, J. (1984). Regeneration and Pattern Formation in Planarians. I. The Pattern of Mitosis in Anterior and Posterior Regeneration in Dugesia (G) Tigrina, and a New Proposal for Blastema Formation. J. Embryol. Exp. Morphol. 83 (October), 63–80. doi:10.1242/dev.83.1.63
Sánchez Alvarado, A., and Newmark, P. A. (1999). Double-Stranded RNA Specifically Disrupts Gene Expression during Planarian Regeneration. Proc. Natl. Acad. Sci. 96 (9), 5049–5054. doi:10.1073/pnas.96.9.5049
Sandmann, T., Vogg, M. C., Owlarn, S., Boutros, M., and Bartscherer, K. (2011). The Head-Regeneration Transcriptome of the Planarian Schmidtea Mediterranea. Genome Biol. 12 (8), R76. doi:10.1186/gb-2011-12-8-r76
Scimone, M. L., Lapan, S. W., and Reddien, P. W. (2014). A Forkhead Transcription Factor Is Wound-Induced at the Planarian Midline and Required for Anterior Pole Regeneration. Plos Genet. 10 (1), e1003999. doi:10.1371/journal.pgen.1003999
Scimone, M. L., Wurtzel, O., Malecek, K., Fincher, C. T., Oderberg, I. M., Kravarik, K. M., et al. (2018). FoxF-1 Controls Specification of Non-body Wall Muscle and Phagocytic Cells in Planarians. Curr. Biol. 28 (23), 3787–3801. e6. doi:10.1016/j.cub.2018.10.030
Tewari, A. G., Owen, J. H., Petersen, C. P., Wagner, D. E., and Reddien, P. W. (2019). A Small Set of Conserved Genes, Including Sp5 and Hox, Are Activated by Wnt Signaling in the Posterior of Planarians and Acoels. Plos Genet. 15, e1008401. doi:10.1371/journal.pgen.1008401
Vogg, M. C., Owlarn, S., Pérez Rico, Y. A., Xie, J., Suzuki, Y., Gentile, L., et al. (2014). Stem Cell-dependent Formation of a Functional Anterior Regeneration Pole in Planarians Requires Zic and Forkhead Transcription Factors. Develop. Biol. 390 (2), 136–148. doi:10.1016/j.ydbio.2014.03.016
Wagner, D. E., Wang, I. E., and Reddien, P. W. (2011). Clonogenic Neoblasts Are Pluripotent Adult Stem Cells that Underlie Planarian Regeneration. Science 332 (6031), 811–816. doi:10.1126/science.1203983
Wenemoser, D., Lapan, S. W., Wilkinson, A. W., Bell, G. W., and Reddien, P. W. (2012). A Molecular Wound Response Program Associated with Regeneration Initiation in Planarians. Genes Dev. 26 (9), 988–1002. doi:10.1101/gad.187377.112
Wenemoser, D., and Reddien, P. W. (2010). Planarian Regeneration Involves Distinct Stem Cell Responses to Wounds and Tissue Absence. Develop. Biol. 344 (2), 979–991. doi:10.1016/j.ydbio.2010.06.017
Witchley, J. N., Mayer, M., Wagner, D. E., Owen, J. H., and Reddien, P. W. (2013). Muscle Cells Provide Instructions for Planarian Regeneration. Cel Rep. 4 (4), 633–641. doi:10.1016/j.celrep.2013.07.022
Wurtzel, O., Cote, L. E., Poirier, A., Satija, R., Regev, A., and Reddien, P. W. (2015). A Generic and Cell-type-specific Wound Response Precedes Regeneration in Planarians. Develop. Cel 35 (5), 632–645. doi:10.1016/j.devcel.2015.11.004
Keywords: FOXK, planarian, nervous system, regeneration, stem cell, differentiation
Citation: Coronel-Córdoba P, Molina MD, Cardona G, Fraguas S, Pascual-Carreras E, Saló E, Cebrià F and Adell T (2022) FoxK1 is Required for Ectodermal Cell Differentiation During Planarian Regeneration. Front. Cell Dev. Biol. 10:808045. doi: 10.3389/fcell.2022.808045
Received: 02 November 2021; Accepted: 25 January 2022;
Published: 22 February 2022.
Edited by:
Rosa Barrio, CIC bioGUNE, SpainReviewed by:
Ricardo M. Zayas, San Diego State University, United StatesJose Maria Martin-Duran, Queen Mary University of London, United Kingdom
Copyright © 2022 Coronel-Córdoba, Molina, Cardona, Fraguas, Pascual-Carreras, Saló, Cebrià and Adell. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Teresa Adell, dGFkZWxsY0B1Yi5lZHU=; Francesc Cebrià, ZmNlYnJpYUB1Yi5lZHU=
†These authors have contributed equally to this work and share first authorship
‡These authors share last authorship