Distinct patterns of gene expression during regeneration and asexual reproduction in the annelid Pristina leidyi

Received: 9 December 2021 | Revised: 3 April 2022 | Accepted: 4 May 2022 DOI: 10.1002/jez.b.23143 RESEARCH ARTICLE Distinct patterns of gene expression during regeneration and asexual reproduction in the annelid Pristina leidyi Irene del Olmo 1 |Aida Verdes 2 |Patricia Álvarez‐Campos 1 1 Department of Biology (Zoology), Universidad Autónoma de Madrid, Madrid, Spain 2 Department of Biodiversity and Evolutionary Biology, Museo Nacional de Ciencias Naturales de Madrid, Madrid, Spain Correspondence Patricia Álvarez‐Campos, Department of Biology (Zoology), Universidad Autónoma de Madrid, Madrid, Spain, Email: [email protected] Funding information H2020 Marie Sklodowska‐Curie Actions, Grant/Award Number: 841576; Comunidad de Madrid, Grant/Award Number: SI1/PJI/ 2019‐00532; Universidad Autónoma de Madrid, Grant/Award Number: 3931‐326F‐ 726DP6A38‐4F77; European Molecular Biology Organization, Grant/Award Number: ALTF‐217‐2018 Abstract Regeneration, the ability to replace lost body parts, is a widespread phenomenon in the animal kingdom often connected to asexual reproduction or fission, since the only difference between the two appears to be the stimulus that triggers them. Both developmental processes have largely been characterized; however, the molecular toolkit and genetic mechanisms underlying these events remain poorly unexplored. Annelids, in particular the oligochaete Pristina leidyi, provide a good model system to investigate these processes as they show diverse ways to regenerate, and can reproduce asexually through fission under laboratory conditions. Here, we used a comparative transcriptomics approach based on RNA‐sequencing and differential gene expression analyses to understand the molecular mechanisms involved in anterior regeneration and asexual reproduction. We found 291 genes upregulated during anterior regeneration, including several regeneration‐related genes previously reported in other annelids such as frizzled,paics, and vdra. On the other hand, during asexual reproduction, 130 genes were found upregulated, and unexpectedly, many of them were related to germline development during sexual reproduction. We also found important differences between anterior regeneration and asexual reproduction, with the latter showing a gene expression profile more similar to that of control individuals. Nevertheless, we identified 35 genes that were upregulated in both conditions, many of them related to cell pluripotency, stem cells, and cell proliferation. Overall, our results shed light on the molecular mechanisms that control anterior regeneration and asexual reproduction in annelids and reveal similarities with other animals, suggesting that the genetic machinery controlling these processes is conserved across metazoans. KEYWORDS Annelida, differential gene expression, paratomic fission, Pristina leidyi, regeneration, reproduction J Exp Zool B Mol Dev Evol. 2022;1–16. wileyonlinelibrary.com/journal/jezb | 1 This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes. © 2022 The Authors. Journal of Experimental Zoology Part B: Molecular and Developmental Evolution published by Wiley Periodicals LLC. 1|INTRODUCTION Regeneration and asexual reproduction are two postembryonic developmental processes that have attracted the interest of biologists for centuries (Bely & Wray, 2001; Sanchez Alvarado, 2000). Regeneration, the ability to replace lost body parts, is a widespread phenomenon in the animal kingdom, present in several metazoan lineages (Bely & Nyberg, 2010) including annelids or segmented worms (Annelida). The ability to regenerate lost body parts is widely distributed across annelids and typically occurs through epimorphosis, which involves cell proliferation in a newly formed specialized structure called blastema that differentiates to restore the missing parts (e.g., Bely & Nyberg, 2010; de Jong & Seaver, 2018; Kostyuchenko & Kozin, 2021; Sanchez Alvarado, 2000). On the other hand, asexual reproduction or fission implies offspring production without the involvement of germ cells or gametes, and it is also a process that has evolved independently in numerous lineages throughout the tree of life (Campagna et al., 2016; Dolmatov et al., 2018; Zattara & Bely, 2016). Within the phylum Annelida, asexual reproduction can be achieved through two main ways: (1) architomy or fragmentation, when the worm first splits apart and then each fragment regenerates the missing structures; or (2) paratomy or paratomic fission, when a new head and tail is developed in a specific segment with high proliferative activity (i.e., fission zone) and subsequently the worm splits apart (Bely & Wray, 2001). Both abilities (regeneration and asexual reproduction) have often been connected in many organisms, since under a morphological perspective, both processes show great similarities, and thus the main difference between the two events appears to be the stimulus triggering them. While regeneration occurs following a deleterious incident such as an injury, an internal stimulus triggered by environmental factors such as food availability or temperature is necessary for asexual reproduction to happen (Sanchez Alvarado, 2000). In fact, it has been proposed that fission derives from regeneration, since agametic reproduction has evolved mainly in groups with extensive regeneration abilities, and all animals capable to reproduce asexually are also able to regenerate (Bely & Wray, 2001). Specifically in annelids, it has been suggested that asexual reproduction is a novel developmental process evolved through the cooption of the anterior regeneration process (Bely & Wray, 2001; Zattara & Bely, 2011,2016), and several studies have tried to shed light into the molecular framework underlying these two developmental processes. However, most previous studies focused on the characterization of a few candidate genes (Kostyuchenko & Kozin, 2021 and the references herein), or samples from different stages of regeneration and fission pooled together (Nyberg et al., 2012), hindering the possibility of differentiating the specific mechanisms that regulate each event. In our study, we have analyzed the two processes independently to be able to describe the genetic toolkit underlying each of these developmental mechanisms. However, due to limited resources, we have only been able to include one stage of differentiation, which might offer limited but highly valuable information since this is the first study of its kind. In future studies, we will describe more stages to understand the changes in gene expression throughout the whole duration of these developmental processes. Additionally, the source of cell production for the proliferative activity entailing the generation of new tissue in annelids is still unclear. It has been suggested that, similarly to other invertebrates with huge regenerative abilities such as planarian worms (e.g., Brockes & Kumar, 2008), pluripotent cells replace the missing structures during annelid regeneration. However, in contrast to the typical self‐renewal capacities of planarian stem cells (i.e., neoblasts), it seems that in annelids the pluripotent cells migrate from pre‐ existing segments to the wound site, proliferate in the blastema, and redifferentiate in the missing tissues (Bilello & Potswald, 1974;de Jong & Seaver, 2018; Sugio et al., 2012; Zattara & Turlington, & Bely, 2016). In the case of asexual reproduction, it has been suggested that cells of already differentiated tissues dedifferentiate and migrate from neighboring areas to the fission zone where active cell proliferation occurs (Kostyuchenko et al., 2016; Özpolat & Bely, 2015). The freshwater oligochaete Pristina leidyi (Clitellata, Naididae) represents an outstanding annelid model to study regeneration and asexual reproduction, due to its remarkable anterior and posterior regenerating abilities, asexual reproduction by paratomic fission, and continuous growth under favorable laboratory conditions (Bely & Wray, 2001; Özpolat & Bely, 2015; Smith, 1896; Zattara & Bely, 2011). Previous studies performed by Zattara and Bely (2011) have characterized morphologically the regeneration (Figure 1a) and fission (Figure 1b) processes of P. leidyi recognizing different stages easily identified using a stereomicroscope. Therefore, the phenotypic features and morphological changes that take place during regeneration and asexual reproduction in P. leidyi are well characterized. However, our knowledge about the underlying molecular mechanisms in this species in particular, and in annelids in general, is very limited. Thus, the present study aims to further investigate the molecular machinery involved in regeneration and asexual reproduction in the species P. leidyi, using a transcriptomics approach to study the genes potentially involved in each process, with special attention to pluripotency and stem cells gene markers. To unravel the mechanisms underlying these developmental processes and to evaluate the hypothesized link between the two, we have characterized gene expression patterns during anterior regeneration and paratomic fission, identifying candidate genes that may regulate these processes. Moreover, we have compared both events and identified genes involved in both anterior regeneration and asexual reproduction, to assess whether these two processes are evolutionarily linked. Finally, we have identified candidate genes related to cell pluripotency and/or multipotency that might be involved in these two developmental processes. Moreover, the abundant transcriptome data generated here can provide new insights into annelid development, setting the base for future comparative studies that will enable us to understand how these events, and the cellular types involved in them, have evolved within Annelida and across metazoans. 2 | DEL OLMO ET AL. 2|RESULTS AND DISCUSSION 2.1 |General characterization of de novo assembled transcriptome of P. leidyi Three libraries were prepared for each of the three experimental conditions (i.e., anteriorly regenerating, asexually reproducing, and control individuals). Nonetheless, one anterior regeneration library failed the quality controls, and therefore only eight of the nine libraries were sequenced (see Materials and Methods section). Therefore, sequencing of eight complementary DNA (cDNA) libraries (anterior regeneration, n= 2; asexual reproduction, n= 3; control, n= 3) resulted in approximately 223.6 million raw reads (31.6 GB), of which a total of 218,666,879 (~98%) were retained after trimming low‐quality portions and reads shorter than 30 bp (Supporting Information: Table S1). The filtered reads were assembled de novo using Trinity software generating a reference transcriptome for P. leidyi with 360,928 transcripts, an average GC content of 44.45% and N50 = 3545 bp (Supporting Information: Table S2). Of the total assembled transcripts, 73,959 (20.49%) had a BLAST hit in the UniProt Knowledgebase (UniProtKB) database and were functionally annotated. More detailed information on raw sequencing data, assembly statistics, and functional annotation of the reference transcriptome can be found in Supporting Information: Tables S1 and S2. Our assembled transcriptome was highly complete, with a BUSCO score of 97.8%. The analysis recovered 943 of the total 954 metazoan single‐copy orthologs. These corresponded to 933 complete BUSCOs (160 single‐copy BUSCOs, and 773 duplicated BUSCOs) and 10 fragmented BUSCOs (1% of the recovered BUSCOs). Only 1.2% (11 BUSCOs) of the total metazoan BUSCOs were missing in our transcriptome. These high levels of completeness are similar to those found in other annelid species (Álvarez‐Campos et al., 2019). A summary of the transcriptome completeness statistics can be found in Supporting Information: Table S2. 2.2 |Gene expression patterns during anterior regeneration and asexual reproduction are largely different Over the last few decades, several authors have pointed out developmental similarities between regeneration and asexual reproduction (e.g., Bely & Wray, 2001; Zattara & Bely, 2011). In fact, a FIGURE 1 Anterior regeneration and asexual reproduction in Pristina leidyi. Brief description of major events, schematic representation and bright‐field microscopy images of the different stages defined for (a) anterior regeneration and (b) asexual reproduction through paratomic fission in the annelid P. leidyi. Scale bars = 100 μm (based on Zattara & Bely, 2011). DEL OLMO ET AL. | 3 recent study that examined the presence of regeneration and fission capabilities across Annelida proposed that fission is a process derived from regeneration and that anterior regeneration specifically must be present for fission to evolve (Zattara & Bely, 2016). This would suggest that gene expression patterns of the two processes at a comparable time point should be somewhat similar. Indeed, expression of some genes such as orthodenticle (otx) and engrailed (en), as well as a lumbriculid Lan 3‐2 glycoepitope, has been reported in both regeneration and asexual reproduction in several annelid species (Bely & Wray, 2001; Martinez et al., 2005). We detected a total of 1496 differentially expressed genes (up‐and downregulated) in the pairwise comparisons performed (Supporting Information: Figure S1). Of those, 1285 genes were differentially expressed in the comparison of anteriorly regenerating and control individuals (291 upregulated and 994 downregulated in anterior regeneration), while 211 genes were found differentially expressed in the comparison of asexually reproducing and control individuals (130 upregulated and 81 downregulated in asexual reproduction) (Supporting Information: Table S3). Hierarchical clustering of samples, according to gene expression patterns, showed that asexually reproducing individuals have a genetic profile more similar to control samples (nonregenerating, nonreproducing individuals), compared to anteriorly regenerating individuals (Figure 2a). This better correlation between asexual reproducing and control worms is to be expected since both conditions show complete and entire worms. Conversely, anteriorly regenerating worms lack most of the anterior structures and therefore the genes herein expressed may explain the high number of downregulated genes in this regeneration condition, and thus it would be reasonable to find several anterior identities and head markers downregulated in this condition. However, we have found no differential expression of several anterior neural markers previously reported in adult and larval annelids and other invertebrates as aristaless‐related homeobox, dachshund,empty spiracles homeobox,forebrain zinc‐finger,orthopedia, paired box 3/6,retinal homeobox,six 3/6 homeobox,orsynaptotagmin 1 (Gąsiorowski et al., 2021; Kerbl et al., 2016; Santagata et al., 2012). Even though these genes have been reported in our transcriptome, none of them shows a significatively distinctive expression with the defined thresholds, which could suggest that brain structures are developing in the 72 hpa anterior blastema. In any case, the anterior head gene expression within the downregulated genes in the anterior regenerating comparison needs further research. Moreover, it is worth noting that individuals of P. leidyi maintained under optimal conditions are continuously undergoing fission, and therefore it is possible that the molecular toolkit necessary for asexual reproduction is triggered even before the associated morphological changes can be observed. This could contribute to even larger similarities between the two conditions. In fact, recruitment of the developmental programs involved in asexual reproduction has already been observed in the earthworm Lumbriculus sp. 1 week before the process starts (Martinez et al., 2005). It should be noted that this great difference in the number of differentially expressed genes during regeneration and asexual reproduction when compared to control individuals (1285 genes were differentially expressed during anterior regeneration, whereas only 211 genes were differentially expressed during asexual reproduction; Supporting Information: Table S3) could be suggesting blastema formation during anterior regeneration, which implies the reestablishment of anterior identity and de novo renewal of anterior‐ specific structures, and may demand a more complex molecular toolkit than the development of a new individual from preexisting tissues during asexual reproduction. Nonetheless, further research of the downregulated genes in this anteriorly regenerating comparison, and especially those anterior identity markers, would be needed to explain the great differences found here. Overall, these results indicate that although regeneration and asexual reproduction are substantially similar morphologically, the molecular toolkits controlling them seem to be rather different. These results agree with recent studies suggesting that even though these two processes might share a common origin, there are marked differences that might be due to their divergent evolution (Zattara & Bely, 2011). On the other hand, several common features between the two processes have been noted. By comparing both upregulated subsets, that is, anterior regeneration upregulated genes and asexual reproduction upregulated genes, we have found a statistically significant overlap of 35 genes that are upregulated during anterior regeneration and asexual reproduction processes when compared to control individuals (Figure 2a,b and Supporting Information: Tables S4,S5, and S6A). Among these 35 genes, some of them have been reported in other invertebrates and mammals to show a relevant role during regeneration, such as vitamin D3 receptor A (vdra; Chen et al., 2021), and germline development, including ubiquitin‐like modifier‐activating enzyme 1 (uba1; Yi et al., 2010). Nevertheless, the importance of these overlapped genes relies on those stem‐cell‐ related markers that seem to play a role during anterior regeneration and asexual reproduction, including connector enhancer of kinase suppressor of Ras 2 (cnkr2), cyclin‐Y‐like (ccnyl1), HECT, UBA, and WWE domain‐containing E3 ubiquitin protein ligase 1 (huwe1), and pumilio homolog 1 (pum1) orthologs (Bose et al., 2017; Henderson et al., 2015; Juliano et al., 2010; Vizziano‐Cantonnet et al., 2018; Zeng et al., 2016). The upregulation of these stem cell genes may suggest the presence of pluripotent‐like cell populations mediating both processes. Enrichment analyses showed that all shared upregulated genes were associated with Gene Ontology (GO) terms in the Molecular Function (MF) category, including nucleoside phosphate binding, nucleotide binding, adenyl nucleotide binding, or purine ribonucleotide binding (Figure 2c and Supporting Information: Table S6B). Interestingly, these metabolic pathways related to the biosynthesis of purine nucleotides and DNA/RNA are also enriched in rapidly proliferating cancer cells (Ngoka, 2008; Weber, 1983). Thus, the high demand for nucleic acids during regeneration and asexual reproduction in P. leidyi might indicate an elevated proliferative activity during the development of the blastema and the fission zone, respectively. This is consistent with previous findings showing epidermal cell proliferation and blastema/fission zone growth as the major shared events occurring during these stages (Zattara & 4 | DEL OLMO ET AL. Bely, 2011). The GO term motor activity was also found enriched in both processes, illustrating the importance of cellular reorganization. Other GO terms enriched during these stages of anterior regeneration and asexual reproduction included ubiquitin‐activating enzyme activity, due to the upregulation of uba1 ortholog, which is involved in maintaining genomic integrity during cell proliferation, as well as in different neuronal processes such as neuron differentiation, growth, and development (Lambert‐Smith et al., 2020). This suggests that nervous system development is already an important process during these stages of postembryonic development. Remarkably, we also found a B9 domain‐containing protein 1 (b9d1) ortholog upregulated during both regeneration and asexual reproduction. This gene is responsible for cilia biogenesis and is involved in the hedgehog signaling pathway (hedgehog receptor activity, MF category), showing an important role during regeneration, segmentation, and fission in several organisms (e.g., Dolmatov et al., 2018; Dray et al., 2010; Schnapp et al., 2005). 2.3 |Gene expression patterns and functional enrichment during anterior regeneration We detected a total of 1285 genes differentially expressed during the second stage of anterior regeneration, of which 291 (~23%) were upregulated and 994 (~77%) downregulated (Figure 3a,b and Supporting Information: Table S3). Of these, 210 (~72%) of the upregulated genes and 529 (~53%) of the downregulated genes had a blast hit in the UniProtKB database, and thus were functionally (a) (b) (c) FIGURE 2 Gene expression patterns during anterior regeneration and asexual reproduction. (a) Hierarchically clustered heatmap of the 40 most upregulated, annotated, and nonredundant genes from pairwise comparisons of anteriorly regenerating, asexually reproducing and control individuals. Upregulated genes in both anteriorly regenerating and fissioning comparisons are shown in bold. (b) Scaled Venn diagram showing the overlapped annotated genes upregulated in both regeneration and asexual reproduction conditions. (c) Enriched Gene Ontology terms are associated with the 35 upregulated genes shared between anterior regeneration and asexual reproduction (molecular function category). DEL OLMO ET AL. | 5 annotated (Supporting Information: Table S4). Among the 210 upregulated genes, 42 have been previously cited in different regeneration processes in other species (Supporting Information: Table S7) and 7 have been specifically reported during annelid regeneration. For instance, we found frizzled class receptor 1 (fzd1) upregulated, one of the major receptors in the wnt signaling pathway, which is involved in several developmental processes including regeneration. This gene has already been reported during fission and regeneration in P. leidyi and in several other organisms such as sea cucumbers and newts (e.g., Girich et al., 2017; Nyberg et al., 2012; Singh et al., 2018). Our results also showed upregulation of brain tumor (brat), glutamate dehydrogenase 1 (GLUD1), multifunctional protein ADE2 (PUR6, also known as paics), ABC transporter G family member 21 (AB21G), vdra, and some genes from the Sox family (a) (b) (c) (d) FIGURE 3 Gene expression patterns during anterior regeneration. (a) Volcano plot displaying the –log 10 pvalue (false discovery rate [FDR]) as a function of fold change in the regenerating and control individuals. Labeled genes are discussed in the text. (b) Hierarchically clustered heatmap of the most important upregulated annotated and nonredundant genes in this comparison, and their categories according to the main function discussed. (c) Gene Ontology enrichment analysis of the upregulated annotated genes. (d) Gene ontology enrichment analysis of the downregulated annotated genes. 6 | DEL OLMO ET AL. (i.e., Sox‐14 and Sox‐21), which have been previously found upregulated during regeneration in other metazoans, including several annelid species such as Enchytraeus japonensis,Eisenia fetida, Syllis gracilis,Sphaerosyllis hystrix, and Lumbriculus variegatus (Bhambri et al., 2018; Gómez et al., 2018; Myohara et al., 2006; Ribeiro et al., 2019; J. Tao et al., 2019; Tellez‐Garcia et al., 2021). Several regeneration‐related genes previously reported in nonannelid species, including fish, amphibians, mammals, or planarians, were also found upregulated during anterior regeneration in P. leidyi. Some of these genes include pim3,splA/ryanodine receptor domain‐ containing SOCS box 3 (spsb3), or vdra, which have been reported in limb, fin, or caudal regeneration of amphibians and zebrafish (Baddar et al., 2021; Chen et al., 2021; Ivanova et al., 2018). We also found upregulation of cyclin‐F(ccnf), ccnyl1, and growth arrest‐specific 6 (gas6), which have been also reported during mammals' liver regeneration (Couchie et al., 2005; L. Huang et al., 2015; Pibiri, 2018). In addition, we also found upregulated genes involved in different aspects of muscle and nerve development in planarians, zebrafish, amphibians, and mammals, including bystin,innexin Unc‐7(unc‐7), kif19,GLUD1, and gas6 (Gibbs et al., 2011; Güiza et al., 2018;T.C. Huang et al., 2020; Sheng et al., 2004; Stratton et al., 2018). The upregulation of these genes highlights the importance of nervous system development during this stage of regeneration, as suggested in previous studies (Zattara & Bely, 2011). Furthermore, we also found upregulation of glutamine/glutamate metabolism genes (glutamine synthetase, GLUD, glutaminase), which recent transcriptomic studies in other annelid species have revealed to be important in the early stages of regeneration (e.g., Ribeiro et al., 2019; J. Tao et al., 2019; Tellez‐Garcia et al., 2021). Interestingly, we have also found several upregulated genes that have been previously linked to neoblast and stem cell regulation and development. For instance, rad54B,huwe1, and DEAD‐box helicase 23 (ddx23) are neoblast‐specific genes required for different aspects of neoblast cell cycle progression and regulation in planarias (Galloni, 2012; Henderson et al., 2015; Lei et al., 2016). In addition, our results showed upregulation of bmi1, an epigenetic regulator of cell cycle and self‐renewal capacity of mammalian stem cells, required for neoblast function in planarians (Supporting Information: Table S7) (e.g., Önal et al., 2012; Robson et al., 2011), tissue regeneration in several mammals (e.g., Fukuda et al., 2012), and recently reported during early regeneration of the annelid L. variegatus (Tellez‐Garcia et al., 2021). In addition, we also found upregulation of the translational repressor pum1, whose expression in planarians seems to be restricted to neoblasts and is related to cell self‐renewal in mammals (Salvetti et al., 2005; Spassov & Jurecic, 2003). Additional markers of pluripotency from other cell types have also been found upregulated in regenerating P. leidyi individuals, including adaptor‐related protein complex 2 (ap2a2),alpha 2 subunit,CCR4‐NOT transcription complex subunit 6‐like‐B(CN6LB), and fam184A (Elmén et al., 2020; Kokkaliaris et al., 2012; Tauran et al., 2019). The presence of genes related to stem cell maintenance and pluripotency, and more specifically those markers exclusive to neoblasts, suggests the presence of stem cells and their relevance during P. leidyi regeneration. These results provide additional evidence supporting previous studies that indicate the presence of pluripotent cell populations able to migrate and participate in regeneration processes in annelids (Supporting Information: Table S7) (e.g., Bilello & Potswald, 1974; de Jong & Seaver, 2018; Sugio et al., 2012). However, the presence of these genes does not imply that stem cells are the only source of cellular production during tissue renewal, and it might also be possible that dedifferentiated tissues make a significant contribution during annelid regeneration as it has been previously suggested (Bely, 2014). Surprisingly, we also found an ovochymase‐1(ovch1) ortholog upregulated in anteriorly regenerating P. leidyi individuals, a gene typically involved in oogenesis in other invertebrates including the annelid Syllis magdalena, where it may assist maturation of oocytes and prevent self‐fertilization (Álvarez‐Campos et al., 2019). In addition, bmi1 and pum1, required for germ stem cell maintenance and germline development and maturation in several species (Juliano et al., 2010; Komai et al.,2014), were also found upregulated. Germ cell precursors have been identified during anterior regeneration in other annelids (Tadokoro et al., 2006), and thus the upregulation of ovch1,bmi1, and pum1 may indicate the restoration of lost gonadal tissue that takes place during anterior regeneration in P. leidyi. GO analysis of the upregulated genes once again shows a significant enrichment of purine metabolism (Figure 3b and Supporting Information: Table S8), which is required for nucleotide biosynthesis (Ngoka, 2008; Weber, 1983) and is also necessary for the high proliferative activity that takes place in the blastema during annelid regeneration (Bely, 2014). In addition, other enriched GO terms in the Cellular Component (CC) category, such as cyclin/CDK‐ positive transcription elongation factor complex, chromatin, and chromosome, confirm the importance of this proliferative activity through the regulation of the cell cycle (Malumbres & Barbacid, 2009). Other upregulated GO terms in the MF category, including catalytic activity and kinase activity and cell junction in the CC category, have also been previously reported during the regeneration of several annelid species (Paul et al., 2021; Ribeiro et al., 2019). The enrichment of GO terms such as nervous system development in the Biological Process (BP) category, as well as other related GO terms in the MF category (e.g., axonal growth cone) point out the importance of neural development during this stage of P. leidyi regeneration. In fact, 72 h after amputation, the horizontal nerves become evident and the central nervous system starts to develop in the new head (Zattara & Bely, 2011). As anteriorly regenerating worms lack most of the anterior tissues and structures, as well as the anterior brain, the downregulation of different anterior identity markers is to be expected. However, none of the genes previously reported (Gąsiorowski et al., 2021; Kerbl et al., 2016; Santagata et al., 2012) has been found to be significantly differentially expressed in this condition. However, among the 529 genes downregulated during anterior regeneration (Supporting Information: Table S5), 14 of them were identified as regeneration‐related genes previously reported in other species (reviewed by Zhao et al., 2016). For instance, Willebrand DEL OLMO ET AL. | 7 factor D and EGF domains (vwde), important for blastema formation in some vertebrates (Leigh et al., 2020), was also found downregulated during early regeneration in the annelid L. variegatus (Tellez‐Garcia et al., 2021), although its role in this process is still unknown. The downregulation of other genes, such as casp3,casp8,orepidermal grow factor receptor (egfr), which modulate apoptosis in annelids and regulate blastema growth in planarias (Campagna et al., 2016; Fraguas et al., 2011; Shao et al., 2020), may indicate the important role of apoptotic processes in the control of cell proliferation (Kostyuchenko & Kozin, 2021) during P. leidyi regeneration. Actually, the downregulation of these pro‐and antiapoptotic genes and other regeneration‐related genes, such as neuropilin and receptor‐interacting serine‐threonine kinase 3, led to the enrichment of several GO terms (BP category) related to the regulation of cell death, which have been shown to be important during anterior regeneration in cnidarians (Chera et al., 2011) (Figure 3c and Supporting Information: Table S8). In addition, we found downregulation of bone morphogenetic protein (bmp) and neuronal nicotinic acetylcholine receptor (chrna), two genes involved in signaling pathways controlling differentiation and development of the nervous system. Contrary to our results, both genes have been previously reported upregulated in other regenerating metazoan species, including annelids (Bhambri et al., 2018; Cho et al., 2009; Nyberg et al., 2012; Gómez et al., 2018), the former in late blastema developing stages during patterning (Bandyopadhyay et al., 2006), and the latter especially upregulated within the first 24 h after amputation (Cho et al., 2009). This suggests that these genes may show a dynamic expression that decreases during the second stage of P. leidyi regeneration, that is, when the blastema is completely formed and still elongating. Other downregulated GO terms in the BP category such as cellular hormone metabolic processes were also found enriched. It is well known that the annelid brain secretes hormones that promote and regulate posterior regeneration (e.g., de Jong & Seaver, 2018). However, during anterior regeneration the brain is absent and thus is incapable of secreting hormones, which might explain the downregulation of these hormone metabolic pathways. The GO term fatty acid metabolic process (BP category) was also enriched due to the downregulation of fatty acid desaturase 2 (fads2), among others, which maintains the structure and function of cell membranes (Los & Murata, 1998) (Figure 3c and Supporting Information: Table S8). This suggests the importance of cell reorganization during P. leidyi regeneration, probably due to the high proliferative activity. 2.4 |Gene expression patterns and functional enrichment during asexual reproduction We detected 211 differentially expressed genes during the second stage of asexual reproduction, of which 130 (~62%) were upregulated, whereas 81 (~38%) were downregulated (Figure 4a and Supporting Information: Table S3). Among these, 100 (~77%) of the upregulated genes and 59 (~73%) of the downregulated genes had a blast hit in the UniProtKB database and thus were functionally annotated (Supporting Information: Table S5). Little is known about asexual reproduction in P. leidyi in particular, and annelids in general, with only a few molecular studies carried out in the phylum (Bely & Wray, 2001; Martinez et al., 2005; Nyberg et al., 2012; Özpolat & Bely, 2015). In this study, we provide a first overview of the molecular toolkit involved in annelid asexual reproduction and greatly expand the knowledge regarding the genetic machinery underlying this process, as only four of the upregulated genes we identified have been previously reported during asexual reproduction in other metazoans and only one of them in annelids. For instance, fibroblast growth factor receptor 2 (fgfr2) was the only gene previously found highly expressed during P. leidyi asexual reproduction (Nyberg et al., 2012). However, its expression cannot be exclusively attributed to asexual reproduction given that it was also found in regenerative worms. In fact, fgfr2 seems to be essential in the initial formation and outgrowth of the vertebrate's limb and fin regeneration blastema (Brockes & Kumar, 2008). A member of the cadherin family (protocadherin‐11 X‐linked,pcdh11x) and tuberin (TSC complex subunit 2,tsc2), a tumor suppressor gene related to stem cell behavior, were found upregulated in this condition and have also been reported during the asexual cycle of the chordate Botryllus schlosseri (Campagna et al., 2016;Ricci et al., 2016). Finally, Pbx/knotted 1 homeobox 2 (pknox2)wasalso found upregulated, a gene whose upregulation has also been reported in fissioning individuals of the sea cucumber Cladolabes schmeltzii (Dolmatov et al., 2018). Although the expression of these genes has been reported during asexual reproduction, their specific role in the process is still unknown. Agrin, a structural component of the extracellular matrix, which also shows a relevant role in the nervous system (Daniels, 2012), was found upregulated during P. leidyi asexual reproduction. Conversely, this gene has been reported to be upregulated in C. schmeltzii not‐fissioning individuals (Dolmatov et al., 2018), which might suggest that in this species, unlike in P. leidyi, nervous system development continues after fissioning has finished. Other genes found upregulated during asexual reproduction were also involved in neurogenesis, nervous system development, neural differentiation, and muscle formation. These genes include ectoderm‐ neural cortex 1 (enc1), arginine‐glutamic acid dipeptide repeats (rere), and triple functional domain (trio), among others (B. J. Kim & Scott, 2014; S. G. Kim et al., 2009; T. Tao et al., 2020). Once again, the expression of these genes may indicate that nervous system development occurs early during asexual reproduction, as suggested in previous studies (Zattara & Bely, 2011). Additionally, we also found several upregulated genes in asexually reproducing individuals that play a relevant role not only during germline maintenance but also during germ cell development and differentiation in sexual reproduction. For instance, DEAD‐box helicase 3 X‐linked (ddx3x) is expressed during both oogenesis and spermatogenesis processes in mice and fish (Matsumura et al., 2019; L. Sun et al., 2020). Genes preferentially expressed during spermatogenesis in different species such as mammals, fish, and invertebrates include some previously reported upregulated, 8 | DEL OLMO ET AL. anteriorly regenerating P. leidyi individuals including pum1,eukaryotic translation initiation factor 4 gamma 3 (eIF4G3), uba1, and bromodomain adjacent to zinc finger domain protein 1A (baz1a) (Dowdle et al., 2013; F. Sun et al., 2010; Yi et al., 2010; Zheng et al., 2005). Besides, other genes exclusively upregulated during this condition such as sperm surface protein Sp17 and intraflagellar transport protein 56 (Dias et al., 2020; Teves et al., 2020) are also required for spermatogenesis. In addition, genes such as pum1 and socs7 have been reported during female germ cell development (oogenesis) (Assou et al., 2006; Juliano et al., 2010; Virant‐Klun et al., 2013). Finally, other stem‐cell‐related genes (e.g., cnkr2,huwe1,ccnyl1) involved in the regulation, differentiation, and/or maintenance of germline precursors (Bose et al., 2017; Vizziano‐Cantonnet et al., 2018; Zeng et al., 2016) were also upregulated during (a) (b) (c) (d) FIGURE 4 Gene expression patterns during asexual reproduction by paratomic fission. (a) Volcano plot displaying the –log 10 pvalue (false discovery rate [FDR]) as a function of fold change in fissioning and control individuals. Labeled genes are discussed in the text. (b) Hierarchically clustered heatmap of the most important upregulated annotated and nonredundant genes in this comparison, and their categories according to the main function discussed. (c) Gene Ontology enrichment analysis of upregulated annotated genes. (d) Gene Ontology enrichment analysis of downregulated annotated genes. DEL OLMO ET AL. | 9 and Engineering: C,110, 110671. https://doi.org/10.1016/j.msec. 2020.110671 Zattara, E. E., & Bely, A. E. (2011). Evolution of a novel developmental trajectory: Fission is distinct from regeneration in the annelid Pristina leidyi.Evolution & Development,13,80–95. https://doi.org/10.1111/ j.1525-142X.2010.00458.x Zattara, E. E., & Bely, A. E. (2016). Phylogenetic distribution of regeneration and asexual reproduction in Annelida: Regeneration is ancestral and fission evolves in regenerative clades. Invertebrate Biology,135,400–414. https://doi.org/10.1111/ ivb.12151 Zattara, E. E., Turlington, K. W., & Bely, A. E. (2016). Long‐term time‐lapse live imaging reveals extensive cell migration during annelid regeneration. BMC Developmental Biology,16,1–21. https://doi. org/10.1186/s12861-016-0104-2 Zeng, L., Cai, C., Li, S., Wang, W., Li, Y., Chen, J., Zhu, X., & Zeng, Y. A. (2016). Essential roles of cyclin Y‐like 1 and cyclin Y in dividing Wnt‐ responsive mammary stem/progenitor cells. PLoS Genetics,12, e1006055. https://doi.org/10.1371/journal.pgen.1006055 Zhao, M., Rotgans, B., Wang, T., & Cummins, S. F. (2016). REGene: A literature‐based knowledgebase of animal regeneration that bridge tissue regeneration and cancer. Scientific Reports,6,1–11. https:// doi.org/10.1038/srep23167 Zheng, P., Patel, B., McMenamin, M., Moran, E., Paprocki, A. M., Kihara, M., Schramm, R. D., & Latham, K. E. (2005). Effects of follicle size and oocyte maturation conditions on maternal messenger RNA regulation and gene expression in rhesus monkey oocytes and embryos. Biology of Reproduction,72, 890–897. https://doi.org/ 10.1095/biolreprod.104.035881 SUPPORTING INFORMATION Additional supporting information can be found online in the Supporting Information section at the end of this article. How to cite this article: del Olmo, I., Verdes, A., & Álvarez‐ Campos, P. (2022). Distinct patterns of gene expression during regeneration and asexual reproduction in the annelid Pristina leidyi.Journal of Experimental Zoology Part B: Molecular and Developmental Evolution,1–16. https://doi.org/10.1002/jez.b.23143 16 | DEL OLMO ET AL.