Skip to main content

Comparative transcriptomic analysis reveals gene regulation mediated by caspase activity in a chordate organism



Apoptosis is a caspase regulated cell death present in all metazoans defined by a conserved set of morphological features. A well-described function of apoptosis is the removal of excessive cells during development and homeostasis. Recent studies have shown an unexpected signalling property of apoptotic cells, affecting cell fate and/or behaviour of neighbouring cells. In contrast to the apoptotic function of cell elimination, this new role of apoptosis is not well understood but seems caspase-dependent. To deepen our understanding of apoptotic functions, it is necessary to work on a biological model with a predictable apoptosis pattern affecting cell fate and/or behaviour. The tunicate Ciona intestinalis has a bi-phasic life cycle with swimming larvae which undergo metamorphosis after settlement. Previously, we have shown that the tail regression step during metamorphosis, characterized by a predictable polarized apoptotic wave, ensures elimination of most tail cells and controls primordial germ cells survival and migration.


We performed differential transcriptomic analysis between control metamorphosing larvae and larvae treated with the pan-caspase inhibitor Z-VAD-fmk in order to explore the transcriptional control of apoptotic cells on neighbouring cells that survive and migrate. When caspase activity was impaired, genes known to be involved in metamorphosis were downregulated along with other implicated in cell migration and survival molecular pathways.


We propose these results as a confirmation that apoptotic cells can control surrounding cells fate and as a reference database to explore novel apoptotic functions in animals, including those related to migration and differentiation.

Peer Review reports


Apoptosis is a regulated cell death defined by morphological features and depending on caspases [1,2,3]. Being ubiquitous in metazoans [4,5,6,7], the main function of apoptosis is the elimination of unwanted cells [8], such as the interdigital tissues during embryogenesis of tetrapods [9] or larval tissues during metamorphosis of many animals [10]. Apoptosis is also implicated in homoeostasis by fine-tuning the balance between cell death and proliferation like in gonad size control in rodents [11].

While the function of cell elimination by apoptosis is highly documented, apoptosis seems to have a broader morphogenetic function that goes beyond cells elimination. Accumulating evidence in many metazoans suggest that apoptotic cells can emit caspase-dependent signals to their neighbours, modulating their fate (survival/death, differentiation or proliferation) or behaviour (migration) [12, 13]. Apoptotic-induced proliferation was reported in several organisms, such as in the cnidarian Hydra, where apoptosis promotes proliferation of adjacent cells through emission of Wnt3 signalling in a caspase-dependant manner during head regeneration [14]. Similar observations have been made during rodent liver regeneration, where cleavage of iPLA2 by caspases in apoptotic cells led to proliferation of adjacent cells [15]. Apoptotic cells can also induce differentiation, as was reported during the metamorphosis of the cnidarian Hydractinia echinata [16]. Finally, apoptotic cells modulate death or survival balance, as was observed in the Drosophila imaginal disc, by producing the Tie ligand Pvf1, that induces death resistance of neighbouring cells [17]. Apoptotic-dependent migration is well described in mammals, in which apoptotic cells secrete directional signals such as Lysophophatidylcholine or Sphingosine-1-Phosphate, promoting migration of leukocytes [18, 19]. Furthermore, in Xenopus laevis tadpole, an apoptotic-dependent axon guidance was suggested during tail regeneration after injury [20].

While caspase-dependent signals emitted by apoptotic cells start to be characterised, almost no target genes are identified in receptor cells. Evaluating the capacity of caspases to affect gene expression and influence fate and/or behaviour of cells adjacent to apoptosis is the first fundamental step to better characterise this new function in animals.

In this report, we took advantage of the predictable spatiotemporal apoptotic profile of ascidians [21, 22], in particular during the Ciona intestinalis metamorphosis, which is concomitant with cell survival and cell migration [21, 23, 24]. Ciona belongs to the Urochordata, the sister group of vertebrates [25], and has a biphasic life cycle composed of a pelagic larva and a benthic adult. Embryogenesis generates a tadpole-like swimming larva consisting of an anterior trunk and a posterior tail which, after few hours of swimming, settles and metamorphoses into a sessile juvenile. A hallmark of metamorphosis is the apoptotic-dependent tail regression [21] during which apoptosis starts at the tip of the tail and propagates toward the trunk, making it easily anticipated; a unique feature among chordates. Using TUNEL labelling and electron microscopy, we previously described in detail this apoptotic profile. Apoptosis starts in epidermal cells at the most posterior tail tip, then, it affects the notochord cells and the striated muscle cells in a postero-anterior wave. The tunic cells are also randomly affected. On the opposite, two tissues survive and migrate during, or prior to, tail regression; the ventral endodermal strand and the eight primordial germ cells (PGC) [23, 24, 26]. PGC, located at the postero-ventral side of the tail tip, escape from cell death by moving toward the trunk, and are consequently surrounded by apoptotic cells during the tail regression. Interestingly, we previously demonstrated that the pan-caspase inhibitor Z-VAD-Fmk blocks both apoptotic dependent-tail regression and PGC migration, suggesting that PGC migration is regulated by apoptotic cells [23].

To identify genes controlled by caspases activity during tail regression, we performed a comparative transcriptomic analysis at the beginning and at mid-tail regression between control metamorphosing larvae or larvae exposed to the pan-caspase inhibitor Z-VAD-fmk, already successfully used in Ciona [21] (Fig. S1). Based on a de novo transcriptome, we identified 61 genes differentially expressed between the beginning and mid-tail regression in control larvae. Next, we found 65 genes differentially expressed between control and treated larvae (21 and 44 at the beginning and at mid-tail regression, respectively). Interestingly, misregulation of some genes upon caspase inhibition was reported during ascidian metamorphosis, but also reported in cell survival and cell migration in several animals [27,28,29,30]. Our results offer a first database of genes transcriptionally modulated by apoptosis during a dynamic morphogenetic process in a chordate, allowing exploration of novel apoptotic functions in animals.


Genes differentially express between the onset and the mid-tail regression

Globally, the robustness of our analysis is shown by a heat map allowed visualization of result congruence between replicates (Fig. 1). To confirm relevance of our data, we randomly chose 6 genes differentially expressed in our transcriptomic analysis to performed real-time PCR (Fig. 2A). Five of them exhibited an expression profile that was similar for both experiments.

Fig. 1
figure 1

A, Counting of genes differentially expressed between conditions. B, Heatmaps showing genes differentially expressed in control between the beginning of tail regression (TR) and the mid-regression. C, Heatmaps showing genes diffenrentially expressed at the beginning of tail regression between control and Z-VAD-fmk treated larvae. D, Heatmaps showing genes diffenrentially expressed at mid-regression between control and Z-VAD-fmk treated larvae

Fig. 2
figure 2

A, Control of gene expression by real-time PCR. Real-time PCR results are consistent with our differential transcriptomic analysis, with however expression level of KH.C7.314 between beginning and mid-tail regression of control larvae not statistically significant. BC: Tail regression beginning, control. BZ: Tail regression beginning, Z-VAD. MC: Mid-tail regression, control. MZ: Mid-tail regression, Z-VAD. * indicates statistically significant differences. B, In situ hybridization of KH.L132.4 (matrilin-2) gene (in green) combined with immunnodection of Ci-VASA to detect PGC (in red) in premetamorphic larva. Matrilin-2 is expressed in epidermal cells at the tip of the tail, the most anterior expressing epidermal cells are in contact with PGC (higher magnification, white square). Matrilin-2 expression is also detected in the trunk epidermis. White square higher magnification of the PGC region

One of them, KH.C7.314, is less expressed at the mid-tail regression in comparison to the beginning of the process in control larvae in our transcriptomic analysis. Similar profile was observed with real-time PCR, but the difference is less pronounced and not statistically significant. Despite this, control by real-time PCR globally validated our transcriptomic data.

First, we have focused on genes differentially expressed in control metamorphosing larvae between the beginning and the mid-tail regression. Twenty-five genes (41%) were up-regulated from the beginning to the mid-tail regression with the strongest increase recorded for KH.C9.385 (smyd3) (Table 1). Conversely, 36 genes (59%) show decreased expression, with a subset of genes whose expression was no longer detectable such as KH.14.261 (lysozyme-like) (Table 2). Among these identified genes, several were already reported during a previous cDNA chip analysis of Ciona metamorphosis, such as KH.C7.314 (MED30), KH.C11.274 (p-selectin, known as Ci-sushi) which is expressed in almost all the tail, or KH.L132.4 (matrilin-2) expressed in epidermal cells at the tip of the tail in premetamorphic larva (Fig. 2B) [29]. In addition, some of the genes identified by our screen had already been detected by substantial RNA sequencing data coupled with in situ hybridization data available on ANISEED database ( For example, KH.L132.16 (fibrillin-3) starts to be expressed at mid-neurulae with a maximum expression in larva, in which it is detected in all the epidermis including the tail. Furhermore, KH.C10.250 (fox-1, RNA-binding) is expressed in tail muscle at mid-tailbud stage to become ubiquitous in larva [31].

Table 1 Genes upregulated at mid-tail regression in comparison to the beginning
Table 2 Genes downregulated at mid-tail regression in comparison to the beginning

We also observed that some among these were also impacted by Z-VAD-fmk treatment (Tables 3, 4, 5 and 6) such as KH.C6.197 (meta2), KH.C1.10 (LIM), KH.C1.216 (myosin light chain 3), and KH.C11.274 (Ci-sushi) [27,28,29,30]. The detection of genes already known to be expressed in the tail and involved in tail regression (i.e. Ci-sushi, Ci-meta2) is consistent with previous studies on ascidians metamorphosis, validates our experimental procedures and supports the relevance of our database.

Table 3 Genes upregulated by the Z-VAD-fmk at the tail regression beginning
Table 4 Genes downregulated by the Z-VAD-fmk at the tail regression beginning
Table 5 Genes upregulated by the Z-VAD-fmk at mid-tail regression
Table 6 Genes down regulated by the Z-VAD-fmk at mid-tail regression

Caspase inhibition affects gene expression during the tail regression

In total, 65 genes were affected by the inhibition of caspases during Ciona tail regression (Fig. 1). At the beginning of tail regression 12 genes were upregulated by the Z-VAD-fmk treatment (so negatively regulated by caspase activity in physiological conditions) with the greatest expression increase observed for KH.C11.506 (girdin) and KH.S1380.1 (mucin-5 AC) (Table 3), and 9 genes were downregulated (so positively regulated by caspase activity) with extreme effects on KH.C1.66, which is no longer expressed (Table 4).

At mid-tail regression, we identified 38 upregulated genes in the Z-VAD-fmk condition with the highest increase noted for KH.C12.516 (Table 5), and 6 downregulated genes (Table 6) with strongest effect concerning KH.C1.910.

Transcription is necessary for the PGCs migration

We previously demonstrated that Z-VAD-fmk blocked both tail regression and PGC migration, and that apoptotic wave propagation and PGC migration speed were correlated [23]. Here, we show that transcription is necessary for PGC survival and migration (Fig. 3). We exposed metamorphosing larvae to the transcription inhibitor actinomycin-D; in the control (devoid of actinomycin-D), a classical apoptotic profile is observed with PGC movement, as already described [23], and in actinomycin-D-treated larvae PGC migration is blocked and PGC nuclei become TUNEL-positive, indicating destruction as the apoptotic wave progresses (Fig. 3). This result, and our previous studies, argue that survival and migration of PGC depend on transcription that could be partially controlled by caspase activity from apoptotic cells (Fig. 3). The comparative transcriptomic analysis between control and Z-VAD-fmk treated larvae is fundamental to understand this apoptotic dependent migration.

Fig. 3
figure 3

A, Larvae labelled by TUNEL (green) with VASA immunostaining (red) and counterstained with DAPI (blue). PGC of larvae treated with the transcription inhibitor actinomycin-D (1 μg/ml) show impaired migration and undergo apoptosis. B, Representation of apoptotic functions during tail regression. Apoptotic cells emit signals modifying transcription of adjacent cells leading to PGC survival and migration to the trunk


Caspase activity modulates expression of genes implicated in cell death/survival

Interestingly, the 2 genes most inhibited by caspase activity at the beginning of tail regression, girdin and mucin-5 AC, were both reported to promote survival in mammalian cells [32, 33]. Among the other genes negatively impacted by caspase activity, KH.C1.1084, the thioredoxin homologue in mammals, is able to control the activation of the executioner caspase 3, the initiator caspase 9, and the release of the cytochrome c, all characterising the intrinsic apoptosis [34,35,36,37]. In addition, KH.C3.45 (kinase C alpha) homologue is known to promote survival and increases anti-apoptotic Bcl-2 efficiency in cultured cells [38]. Taken together, our database highlights that cells expressed a variety of genes regulating cell death/survival fate in a caspase dependent manner, suggesting that apoptotic cells actively participate in apoptosis propagation during tail regression.

Caspase activity modulates expression of genes implicated in cell migration

We identified numerous genes modulated by caspase activity, with a majority of them upregulated, previously observed in cell migration regulation. It is the case of KH.C1.216 (myosin light chain 3) and KH.C9.494 (myosin light chain kinase) for which homologues are key players in mesenchymal migration [32, 39] and KH.C13.153 (uromodulin-like) which provides guidance cues for neutrophile migration [40]. Furthermore, integrin subunits and components of the Extra Cellular Matrix, crucial for mesenchymal migration, are well represented in our transcriptomic approach. In addition, KH.C1.10 (lasp1) mammalian homolog is necessary for migration of NIH 3 T3 cells [41, 42].

Furthermore, KH.C12.516 (ceramidase) could be involved in tissue migration regulation, as it has been described as participating in the sphingolipid pathways. Interestingly, sphingosine 1-phosphate signalling leads to tissue migration [43], particularly in germline cells (due to to a sphingosine 1-phosphate receptor) in the tunicate Botryllus [44], making this pathway a strong candidate for PGC migration in Ciona. Given that these genes and PGC migration are both controlled by caspases [23], PGC migration can be regulated by apoptotic cells.

Ascidians as crucial biological models to studies novel apoptotic function

This study, coupled with our previous work suggested that apoptosis regulates the migration of primordial germ cells [23], confirm Ciona as a perfect biological model to understand apoptotic-induced migration. In addition, a regulation of cell proliferation by apoptosis was reported in this species [45]. Ascidians are well known to have an impressive regenerative capacity for a chordate, in particular in the context of the reconstruction of the distal body parts after amputation, thanks to stem cell requisition from the branchial sac [46,47,48]. Apoptosis was reported to occur as a transient event at the site of injury, and precedes the regeneration of the neural complex, the oral siphon, or the atrial siphon [49]. In the case of siphon amputation, it has been shown by caspases inhibition experiment that an apoptotic-dependant Wnt signalling allows the cell proliferation, leading to the successful regeneration [45]. In addition, apoptosis appears crucial to maintain the growth of branchial sac and its inhibition lead to a hazardous organization of this tissue, suggesting a fundamental function in homeostasis maintaining. Importantly, stem cells requisition and migration from the branchial sac toward the wounding site seem also dependant on caspases, making it an other potential apoptotic-dependant migration [45]. Taken together, all of these studies make Ciona one of the rare biological models where several cell fates and behaviours, controlled by apoptotic cells in a caspases dependent manner, has been put in evidence or suggested in adults as well as in larvae.

In the colonial ascidian Botryllus schlosseri, apoptosis is known to play on budlets elimination [50,51,52]. In this species, blastogenetic cycle (asexual reproduction of adults) give rise to a bud, itself given next four budlets (reviewed in [53, 54]). Usually, two of them growth despite the two other undergo apoptosis and are reabsorbed. In addition to this function in the blastogenetic cycle, apoptosis was also reported in tunic, epidermis, branchial leaflet, vessel epithelium, circulatory system or gonads in adults. With so many tissues presenting apoptosis in physiological condition, it seems clear that apoptosis should have crucial functions in Botryllus and likely could participate to regulation of other cell fates and behaviours as it was reported in Ciona. All of these tissues are as many opportunities to deepen or understanding of apoptosis functions in a same species.


Here, we offer a potential roadmap in order to design future functional studies on Ciona and other animals, including mammals, where apoptosis takes place together with other cell fates and behaviours. Identification of molecular actors needing cleavage for their activation, cell receptors receiving signals from apoptotic cells, or transcription factors could be promising for further research.

The last decade has brought evidence that the cleavage capacity of caspases is a prerequisite to activate molecular actors promoting migration or differentiation in a non-apoptotic context [55,56,57,58,59]. Consequently, Z-VAD-fmk blocks all processes led by caspases, making it crucial to distinguish the level of implication of apoptotic and non-dying cells for future research.

Finally, our study provides an additional and novel chapter in the research of apoptotic functions beyond their role as a simple cell elimination mechanism, opening opportunities in the field of medical science [60].


Sample collection

Adult Ciona intestinalis were collected in the field in Brittany (France) by the Centre National Marine Biological Resources, EMBRC-France, Roscoff (Station Biologique de Roscoff, Sorbonne Université, Brittany, France). Individuals were maintained in 35‰ salinity artificial seawater (ASW) at 18 °C in the UMR7138 Evolution Paris Seine laboratory (Sorbonne University, Paris, France). Oocytes and sperm were obtained by dissection and cross-fertilizations were performed in plastic Petri dishes. Embryos were cultured at 18 °C in 0.2 μm filtered ASW with 0.1 M Hepes (#H4034, Sigma-Aldrich; Merck KGaA). Embryos of each fertilisation were divided into two subsets before the first cell division to obtain one control and one treatment condition. Six fertilisations were performed, divided into twelve petri dishes (6 control with their respective 6 treatment conditions). After the treatment, larvae were collected at two time points (see below), and compacted by manual centrifugation before RNA extraction.

Pharmacological treatments and TUNEL labelling

The pan-caspase inhibitor Z-VAD-Fmk (#V116, Sigma-Aldrich; Merck KGaA) was stored in DMSO (#D8418, Sigma-Aldrich; Merck KGaA) at − 20 °C and used at final concentration of 10 μM. After the settlement of more than 50% of larvae on the dish, and when some of them stopped moving (indicating that metamorphosis will begin soon), supernatant was discarded and replaced with filtered ASW containing DMSO as control, or Z-VAD-Fmk diluted in DSMO as treatment. Z-VAD-Fmk was renewed every two hours. The experiment was divided into two time points (with three biological replicates for each one – itself divided between control and treated larva) following control tail regression stages: first when the tail regression just began in the control (beginning), and the second at mid-progression (mid-tail regression) given three controls with three Z-VAD-Fmk equivalents for each point.

In a separate experiment, 1 μg/ml of actinomycin-D (#A1410, Sigma-Aldrich; Merck KGaA) was added to swimming larvae. After settlement and tail regression initiation, paraformaldehyde fixation, TUNEL labelling and Vasa immunostaining were performed as previously described [23].

RNA preparation and sequencing

RNA was extracted using the RNAqueous™-Micro Total RNA Isolation Kit (#AM1931 – ThermoFisher Scientific) according to the manufacturer’s protocol. DNAse treatment was performed using the TURBO DNA-free™ Kit (#AM1907 - ThermoFisher Scientific). RNA was purified and concentrated using the RNeasy MinElute Cleanup Kit (#74204 – Qiagen) according to the manufacturer’s protocol. RNA quantity was measured using a NanoDrop and quality was assessed with the Experion RNA HighSens (#700–7155 – BioRad) from 1 μL (0.1 μL of RNA completed with 0.9 μL of Milli-Q water) following standard instructions. Transcriptomes were sequenced on Illumina Hi-Seq 2000 (150 bp reads) at the Sistemas Genomicos plate-form of the Ascires Biomedical Group (Spain).

De novo transcriptome assembly and differential expression

Raw reads of the twelve Ciona intestinalis samples were checked with FastQC and merged to assemble a de novo transcriptome with Oases [61] and Trinity [62]. Best assemblies from each one (best N50) were merged with Cap3 to obtain the de novo transcriptome. De novo transcriptome contains 85,506 contigs (N50 = 1314, minimum contig length = 102, maximum contig length = 13,998). Then, raw reads of each sample were mapped against the de novo transcriptome with Bowtie2 [63] with standard parameters. Low quality mapping reads (MQ > 20) were removed using Samtools [64] and Picard Tools ( Expression quantification was performed using HTSeq [65] and differential expression evaluated with the DESeq2 package from R. We kept the genes with an expression ratio - (control 1 + control 2 + control 3) / (treatment 1 + treatment 2 + treatment 3) – above 2 (upregulated) and lower than 0.5 (downregulated). Then, Gene Ontology (GO) enrichment was performed with Blast2GO 5.2 and sequences were blasted on the Ciona genome to identify their KyotoGrail KH gene model.

Real-time PCR

We synthetised the cDNA from RNA coming from same sample which was sending for sequencing using the SuperScript™ II Reverse Transcriptase (18,064,014 - ThermoFisher Scientific) according to the manufacturer protocol. Real-time PCR was performed with the SYBR Green Supermix (Biorad) on a Biorad’s thermal cycler with the following profile: 95 °C for 10 min; 40 cycles of amplification with successively 95 °C for 15 s, 60 °C for 10s, and 72 °C for 20s; one cycle for melting curve analysis with an acquisition every 0.5 °C from 65 °C to 95 °C to verify the presence of a single product. Each assay included four successive dilutions as standard to determine the reaction efficiencies and Ct values. Reactions were done in triplicate. All PCR amplicon have a length from 140 to 160 bp. Normalized was made using the reference gene Ci-actin (forward 5′ ATGTGCAAGGCCGGTTT 3′; reverse 5′ GACACGGAGTTCGTTGT 3′), already successfully used [66]. Genes of interest and their primer pair are: KH.C9.385 with forward 5′ ACTTCTGAAAAGAGCGGACG 3′ and reverse 5′ CATTGCAACAGACCATCTTGC 3′, KH.C7.314 with forward 5′ AGCCAACTACCGAATGGGA 3′ and reverse 5′ GTTCCACGTCTCCAACTCTC 3′, KH.C12.103 with forward 5′ TGTCTGAGATACGAAAGCGT 3′ and reverse 5′ CATTTGCTCAAGATCGGCG 3′, KH.C9.791 with forward 5′ GGATCTAAGCACTGTTCTGGAC 3′ and reverse 5′ CAAGGGCGTTGGTGTTCAGT 3′, KH.C9.170 with forward 5′ AATCCCCGCCCTTGAAGAAA 3′ and reverse 5′ GGGGTGGATATGAATAACATGG 3′, and KH.C1.910 with forward 5′ ACATGAAGATTTTGGTACCGAC 3′ and reverse 5′ CAAAACTTGGCCAAAAAGTTGG 3′. Differences was evaluated by Wilcoxon Mann-Whitney tests using R i386 4.0.3 software, and considered statistically significant for p-value < 0.05.

In situ hybridization combined with vasa immunnodetection

Larvae were fixed 2 h at 4 °C in MEM-PFA (4% paraformaldehyde, 0.1 M MOPS pH 7.4, 0.1 M NaCl, 1 mM EGTA, and 2 mM MgSO4, 0.05% Tween-20). After fixation larvae were washed three times and dehydrated through a graded series of ethanol/PBS baths and stocked in 100% ethanol at − 20 °C. Rehydratation was made by successive washes of ethanol/PBS solution to a full PBS final solution. Hybridization was made according to a previous protocol [67] except that methanol was replaced by ethanol. After hybridization, immunnodetection of VASA was conducted as we previously described [23].

Availability of data and materials

Any request can be addressed to the corresponding author. Raw data are available under BioProject PRJNA725676.


  1. Kerr JF, Wyllie AH, Currie AR. Apoptosis: a basic biological phenomenon with wide-ranging implications in tissue kinetics. Br J Cancer. 1972;26(4):239–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Hengartner MO. The biochemistry of apoptosis. Nature. 2000;407(6805):770–6.

    Article  CAS  PubMed  Google Scholar 

  3. Galluzzi L, Vitale I, Aaronson SA, Abrams JM, Adam D, Agostinis P, et al. Molecular mechanisms of cell death: recommendations of the nomenclature committee on cell death 2018. Cell Death Differ. 2018;25(3):486–541.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Lasi M, David CN, Böttger A. Apoptosis in pre-Bilaterians: Hydra as a model. Apoptosis. 2010;15(3):269–78.

    Article  CAS  PubMed  Google Scholar 

  5. Accordi F, Chimenti C. Programmed cell death in the pancreas of Bufo bufo during metamorphosis. J Anat. 2001;199(Pt 4):419–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Kiss T. Apoptosis and its functional significance in molluscs. Apoptosis. 2010;15(3):313–21.

    Article  PubMed  Google Scholar 

  7. Jiang C, Baehrecke EH, Thummel CS. Steroid regulated programmed cell death during Drosophila metamorphosis. Development. 1997;124(22):4673–83.

    Article  CAS  PubMed  Google Scholar 

  8. Jacobson MD, Weil M, Raff MC. Programmed cell death in animal development. Cell. 1997;88(3):347–54.

    Article  CAS  PubMed  Google Scholar 

  9. Wood W, Turmaine M, Weber R, Camp V, Maki RA, McKercher SR, et al. Mesenchymal cells engulf and clear apoptotic footplate cells in macrophageless PU.1 null mouse embryos. Development. 2000;127(24):5245–52.

    Article  CAS  PubMed  Google Scholar 

  10. Nishikawa A, Hayashi H. Spatial, temporal and hormonal regulation of programmed muscle cell death during metamorphosis of the frog Xenopus laevis. Differentiation. 1995;59(4):207–14.

    Article  CAS  PubMed  Google Scholar 

  11. Furuta I, Porkka-Heiskanen T, Scarbrough K, Tapanainen J, Turek FW, Hsueh AJ. Photoperiod regulates testis cell apoptosis in Djungarian hamsters. Biol Reprod. 1994;51(6):1315–21.

    Article  CAS  PubMed  Google Scholar 

  12. Fogarty CE, Bergmann A. The sound of silence: signaling by apoptotic cells. Curr Top Dev Biol. 2015;114:241–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Fuchs Y, Steller H. Live to die another way: modes of programmed cell death and the signals emanating from dying cells. Nat Rev Mol Cell Biol. 2015;16(6):329–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Chera S, Ghila L, Dobretz K, Wenger Y, Bauer C, Buzgariu W, et al. Apoptotic cells provide an unexpected source of Wnt3 signaling to drive hydra head regeneration. Dev Cell. 2009;17(2):279–89.

    Article  CAS  PubMed  Google Scholar 

  15. Li F, Huang Q, Chen J, Peng Y, Roop DR, Bedford JS, et al. Apoptotic cells activate the “phoenix rising” pathway to promote wound healing and tissue regeneration. Sci Signal. 2010;3:ra13.

    PubMed  PubMed Central  Google Scholar 

  16. Wittig K, Kasper J, Seipp S, Leitz T. Evidence for an instructive role of apoptosis during the metamorphosis of Hydractinia echinata (Hydrozoa). Zoology (Jena). 2011;114(1):11–22.

    Article  Google Scholar 

  17. Bilak A, Uyetake L, Su TT. Dying cells protect survivors from radiation-induced cell death in Drosophila. PLoS Genet. 2014;10(3):e1004220.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Lauber K, Bohn E, Kröber SM, Xiao Y, Blumenthal SG, Lindemann RK, et al. Apoptotic cells induce migration of phagocytes via caspase-3-mediated release of a lipid attraction signal. Cell. 2003;113(6):717–30.

    Article  CAS  PubMed  Google Scholar 

  19. Gude DR, Alvarez SE, Paugh SW, Mitra P, Yu J, Griffiths R, et al. Apoptosis induces expression of sphingosine kinase 1 to release sphingosine-1-phosphate as a “come-and-get-me” signal. FASEB J. 2008;22(8):2629–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Tseng A-S, Adams DS, Qiu D, Koustubhan P, Levin M. Apoptosis is required during early stages of tail regeneration in Xenopus laevis. Dev Biol. 2007;301(1):62–9.

    Article  CAS  PubMed  Google Scholar 

  21. Chambon J-P, Soule J, Pomies P, Fort P, Sahuquet A, Alexandre D, et al. Tail regression in Ciona intestinalis (Prochordate) involves a caspase-dependent apoptosis event associated with ERK activation. Development. 2002;129(13):3105–14.

    Article  CAS  PubMed  Google Scholar 

  22. Jeffery WR. Programmed cell death in the ascidian embryo: modulation by FoxA5 and Manx and roles in the evolution of larval development. Mech Dev. 2002;118(1-2):111–24.

    Article  CAS  PubMed  Google Scholar 

  23. Krasovec G, Robine K, Quéinnec E, Karaiskou A, Chambon JP. Ci-hox12 tail gradient precedes and participates in the control of the apoptotic-dependent tail regression during Ciona larva metamorphosis. Dev Biol. 2019;448(2):237–46.

    Article  CAS  PubMed  Google Scholar 

  24. Shirae-Kurabayashi M, Nishikata T, Takamura K, Tanaka KJ, Nakamoto C, Nakamura A. Dynamic redistribution of vasa homolog and exclusion of somatic cell determinants during germ cell specification in Ciona intestinalis. Development. 2006;133(14):2683–93.

    Article  CAS  PubMed  Google Scholar 

  25. Delsuc F, Brinkmann H, Chourrout D, Philippe H. Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature. 2006;439(7079):965–8.

    Article  CAS  PubMed  Google Scholar 

  26. Kawai N, Ogura Y, Ikuta T, Saiga H, Hamada M, Sakuma T, et al. Hox10-regulated endodermal cell migration is essential for development of the ascidian intestine. Dev Biol. 2015;403(1):43–56.

    Article  CAS  PubMed  Google Scholar 

  27. Nakayama A, Satou Y, Satoh N. Isolation and characterization of genes that are expressed during Ciona intestinalis metamorphosis. Dev Genes Evol. 2001;211(4):184–9.

    Article  CAS  PubMed  Google Scholar 

  28. Nakayama A, Satou Y, Satoh N. Further characterization of genes expressed during Ciona intestinalis metamorphosis. Differentiation. 2002;70(8):429–37.

    Article  CAS  PubMed  Google Scholar 

  29. Chambon J-P, Nakayama A, Takamura K, McDougall A, Satoh N. ERK- and JNK-signalling regulate gene networks that stimulate metamorphosis and apoptosis in tail tissues of ascidian tadpoles. Development. 2007;134(6):1203–19.

    Article  CAS  PubMed  Google Scholar 

  30. Terasaki AG, Hiruta J, Suzuki J, Sakamoto S, Nishioka T, Suzuki H, et al. A lasp family protein of Ciona intestinalis. Biochim Biophys Acta. 1779;2008(1):51–9.

    Article  CAS  Google Scholar 

  31. Kusakabe T, Yoshida R, Kawakami I, Kusakabe R, Mochizuki Y, Yamada L, et al. Gene expression profiles in tadpole larvae of Ciona intestinalis. Dev Biol. 2002;242(2):188–203.

    Article  PubMed  Google Scholar 

  32. Wang S, Lei Y, Cai Z, Ye X, Li L, Luo X, et al. Girdin regulates the proliferation and apoptosis of pancreatic cancer cells via the PI3K/Akt signalling pathway. Oncol Rep. 2018;40(2):599–608.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Hoshi H, Sawada T, Uchida M, Iijima H, Kimura K, Hirakawa K, et al. MUC5AC protects pancreatic cancer cells from TRAIL-induced death pathways. Int J Oncol. 2013;42(3):887–93.

    Article  CAS  PubMed  Google Scholar 

  34. Nonn L, Williams RR, Erickson RP, Powis G. The absence of mitochondrial thioredoxin 2 causes massive apoptosis, exencephaly, and early embryonic lethality in homozygous mice. Mol Cell Biol. 2003;23(3):916–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Tanaka T, Hosoi F, Yamaguchi-Iwai Y, Nakamura H, Masutani H, Ueda S, et al. Thioredoxin-2 (TRX-2) is an essential gene regulating mitochondria-dependent apoptosis. EMBO J. 2002;21(7):1695–703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Baker A, Santos BD, Powis G. Redox control of caspase-3 activity by thioredoxin and other reduced proteins. Biochem Biophys Res Commun. 2000;268(1):78–81.

    Article  CAS  PubMed  Google Scholar 

  37. Mitchell DA, Marletta MA. Thioredoxin catalyzes the S-nitrosation of the caspase-3 active site cysteine. Nat Chem Biol. 2005;1(3):154–8.

    Article  CAS  PubMed  Google Scholar 

  38. Reyland ME. Protein kinase C isoforms: Multi-functional regulators of cell life and death. Front Biosci (Landmark Ed). 2009;14:2386–99.

  39. Wang X, Ye X, Ji J, Wang J, Xu B, Zhang Q, et al. MicroRNA-155 targets myosin light chain kinase to inhibit the migration of human bone marrow-derived mesenchymal stem cells. Int J Mol Med. 2018;42(3):1585–92.

    Article  CAS  PubMed  Google Scholar 

  40. Schmid M, Prajczer S, Gruber LN, Bertocchi C, Gandini R, Pfaller W, et al. Uromodulin facilitates neutrophil migration across renal epithelial monolayers. Cell Physiol Biochem. 2010;26(3):311–8.

    Article  CAS  PubMed  Google Scholar 

  41. Grunewald TGP, Kammerer U, Schulze E, Schindler D, Honig A, Zimmer M, et al. Silencing of LASP-1 influences zyxin localization, inhibits proliferation and reduces migration in breast cancer cells. Exp Cell Res. 2006;312(7):974–82.

    Article  CAS  PubMed  Google Scholar 

  42. Lin YH, Park Z-Y, Lin D, Brahmbhatt AA, Rio M-C, Yates JR, et al. Regulation of cell migration and survival by focal adhesion targeting of Lasp-1. J Cell Biol. 2004;165(3):421–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Spiegel S, English D, Milstien S. Sphingosine 1-phosphate signaling: providing cells with a sense of direction. Trends Cell Biol. 2002;12(5):236–42.

    Article  CAS  PubMed  Google Scholar 

  44. Kassmer SH, Rodriguez D, Langenbacher AD, Bui C, De Tomaso AW. Migration of germline progenitor cells is directed by sphingosine-1-phosphate signalling in a basal chordate. Nat Commun. 2015;6(1):8565.

    Article  CAS  PubMed  Google Scholar 

  45. Jeffery WR, Gorički Š. Apoptosis is a generator of Wnt-dependent regeneration and homeostatic cell renewal in the ascidian Ciona. Biol Open. 2021;10:bio058526.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Jeffery WR. Distal regeneration involves the age dependent activity of branchial sac stem cells in the ascidian Ciona intestinalis. Regeneration (Oxf). 2015;2(1):1–18.

    Article  CAS  Google Scholar 

  47. Auger H, Sasakura Y, Joly J-S, Jeffery WR. Regeneration of oral siphon pigment organs in the ascidian Ciona intestinalis. Dev Biol. 2010;339(2):374–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Jeffery WR. Regeneration, stem cells, and aging in the tunicate Ciona: insights from the Oral siphon. Int Rev Cell Mol Biol. 2015;319:255–82.

    Article  CAS  PubMed  Google Scholar 

  49. Jeffery WR. Progenitor targeting by adult stem cells in Ciona homeostasis, injury, and regeneration. Dev Biol. 2019;448(2):279–90.

    Article  CAS  PubMed  Google Scholar 

  50. Tiozzo S, Ballarin L, Burighel P, Zaniolo G. Programmed cell death in vegetative development: apoptosis during the colonial life cycle of the ascidian Botryllus schlosseri. Tissue Cell. 2006;38(3):193–201.

    Article  CAS  PubMed  Google Scholar 

  51. Franchi N, Ballin F, Manni L, Schiavon F, Basso G, Ballarin L. Recurrent phagocytosis-induced apoptosis in the cyclical generation change of the compound ascidian Botryllus schlosseri. Dev Comp Immunol. 2016;62:8–16.

    Article  CAS  PubMed  Google Scholar 

  52. Ballarin L, Schiavon F, Manni L. Natural apoptosis during the blastogenetic cycle of the colonial ascidian Botryllus schlosseri: a morphological analysis. Zool Sci. 2010;27(2):96–102.

    Article  Google Scholar 

  53. Manni L, Zaniolo G, Cima F, Burighel P, Ballarin L. Botryllus schlosseri: a model ascidian for the study of asexual reproduction. Dev Dyn. 2007;236(2):335–52.

    Article  CAS  PubMed  Google Scholar 

  54. Voskoboynik A, Weissman IL. Botryllus schlosseri, an emerging model for the study of aging, stem cells, and mechanisms of regeneration. Invertebr Reprod Dev. 2015;59(sup1):33–8.

    Article  PubMed  Google Scholar 

  55. Zermati Y, Garrido C, Amsellem S, Fishelson S, Bouscary D, Valensi F, et al. Caspase activation is required for terminal erythroid differentiation. J Exp Med. 2001;193(2):247–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Ishizaki Y, Jacobson MD, Raff MC. A role for caspases in lens fiber differentiation. J Cell Biol. 1998;140(1):153–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Ohsawa S, Hamada S, Kuida K, Yoshida H, Igaki T, Miura M. Maturation of the olfactory sensory neurons by Apaf-1/caspase-9-mediated caspase activity. Proc Natl Acad Sci U S A. 2010;107(30):13366–71.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Denecker G, Ovaere P, Vandenabeele P, Declercq W. Caspase-14 reveals its secrets. J Cell Biol. 2008;180(3):451–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Kanuka H, Kuranaga E, Takemoto K, Hiratou T, Okano H, Miura M. Drosophila caspase transduces shaggy/GSK-3beta kinase activity in neural precursor development. EMBO J. 2005;24(21):3793–806.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Huang Q, Li F, Liu X, Li W, Shi W, Liu F-F, et al. Caspase 3-mediated stimulation of tumor cell repopulation during cancer radiotherapy. Nat Med. 2011;17(7):860–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012;28(8):1086–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.

    Article  CAS  PubMed  Google Scholar 

  66. Vizzini A, Pergolizzi M, Vazzana M, Salerno G, Di Sano C, Macaluso P, et al. FACIT collagen (1alpha-chain) is expressed by hemocytes and epidermis during the inflammatory response of the ascidian Ciona intestinalis. Dev Comp Immunol. 2008;32(6):682–92.

    Article  CAS  PubMed  Google Scholar 

  67. Christiaen L, Wagner E, Shi W, Levine M. Whole-mount in situ hybridization on sea squirt (Ciona intestinalis) embryos. Cold Spring Harb Protoc. 2009;2009:pdb.prot5348.

Download references


The authors strongly acknowledge Muriel Jager (Sorbonne Université, Paris, France), Thierry Jaffredo (Sorbonne Université, Paris, France), and the IBPS Proteomic facility (Sorbonne Université, Paris, France) for technical support. We thank Uri Frank (National University of Ireland Galway, Galway, Ireland) and Helen Horkan (National University of Ireland Galway, Galway, Ireland) for useful comments.


This work was supported by the Fondation ARC pour la recherche sur le cancer grant n° PJA 20141201892 to JP. The work of GK was supported by a Ph.D. fellowship from the French Ministry of Education, Research and Innovation.

Author information

Authors and Affiliations



JP managed the project and made in situ hybridization. GK performed Z-VAD-fmk experiments, RNA extractions, real-time PCR, analysis, created figures and wrote the manuscript. AK performed actinomycin and TUNEL experiments. All authors contributed to comments on the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Gabriel Krasovec.

Ethics declarations

Ethics approval and consent to participate

According to the French legislation (Décret n° 2013–118 du 1er février 2013 relatif à la protection des animaux utilisés à des fins scientifiques), experiments on Ciona intestinalis do not require an approval from Animal Ethics Commission because its not a “live vertebrate animals” and not a “live cephalopods”.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Fig. 1.

Experimental design. Six independent fertilisations were performed, and each of them was divided into two batches, in twelve petri dishes. After settlement, supernatant was discarded and replaced with filtered ASW containing DMSO (control) or Z-VAD-Fmk diluted in DSMO (treatment). RNA sampling was done at two time points according to the control: immediately after observation of tail regression beginning, and at mid-tail regression.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Krasovec, G., Karaiskou, A., Quéinnec, É. et al. Comparative transcriptomic analysis reveals gene regulation mediated by caspase activity in a chordate organism. BMC Mol and Cell Biol 22, 51 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: