A Novel Terminal-Repeat Retrotransposon in Miniature (TRIM) Is Massively Expressed in Echinococcus multilocularis Stem Cells

Helenus and Ajax, Two Groups of Non-Autonomous LTR Retrotransposons, Represent a New Type of Small RNA Gene-Derived Mobile Elements

Terminal repeat retrotransposons in miniature (TRIMs) are short non-autonomous long terminal repeat (LTR) retrotransposons found from various eukaryotes. Cassandra is a unique TRIM lineage which contains a 5S rRNA-derived sequence in its LTRs. Here, two new groups of TRIMs, designated Helenus and Ajax, are reported based on bioinformatics analysis and the usage of Repbase. Helenus is found from fungi, animals, and plants, and its LTRs contain a tRNA-like sequence. It includes two LTRs and between them, a primer-binding site (PBS) and polypurine tract (PPT) exist. Fungal and plant Helenus generate 5 bp target site duplications (TSDs) upon integration, while animal Helenus generates 4 bp TSDs. Ajax includes a 5S rRNA-derived sequence in its LTR and is found from two nemertean genomes. Ajax generates 5 bp TSDs upon integration. These results suggest that despite their unique promoters, Helenus and Ajax are TRIMs whose transposition is dependent on autonomous LTR retrotransposon. These TRIMs can originate through an insertion of SINE in an LTR of TRIM. The discovery of Helenus and Ajax suggests the presence of TRIMs with a promoter for RNA polymerase III derived from a small RNA gene, which is here collectively termed TRIMp3.

Genome-wide transcriptome analysis of Echinococcus multilocularis larvae and germinative cell cultures reveals genes involved in parasite stem cell function

The lethal zoonosis alveolar echinococcosis is caused by tumour-like growth of the metacestode stage of the tapeworm Echinococcus multilocularis within host organs. We previously demonstrated that metacestode proliferation is exclusively driven by somatic stem cells (germinative cells), which are the only mitotically active parasite cells that give rise to all differentiated cell types. The Echinococcus gene repertoire required for germinative cell maintenance and differentiation has not been characterised so far. We herein carried out Illumina sequencing on cDNA from Echinococcus metacestode vesicles, from metacestode tissue depleted of germinative cells, and from Echinococcus primary cell cultures. We identified a set of ~1,180 genes associated with germinative cells, which contained numerous known stem cell markers alongside genes involved in replication, cell cycle regulation, mitosis, meiosis, epigenetic modification, and nucleotide metabolism. Interestingly, we also identified 44 stem cell associated transcription factors that are likely involved in regulating germinative cell differentiation and/or pluripotency. By in situ hybridization and pulse-chase experiments, we also found a new general Echinococcus stem cell marker, EmCIP2Ah, and we provide evidence implying the presence of a slow cycling stem cell sub-population expressing the extracellular matrix factor Emkal1. RNA-Seq analyses on primary cell cultures revealed that metacestode-derived Echinococcus stem cells display an expanded differentiation capability and do not only form differentiated cell types of the metacestode, but also cells expressing genes specific for protoscoleces, adult worms, and oncospheres, including an ortholog of the schistosome praziquantel target, EmTRPMPZQ. Finally, we show that primary cell cultures contain a cell population expressing an ortholog of the tumour necrosis factor α receptor family and that mammalian TNFα accelerates the development of metacestode vesicles from germinative cells. Taken together, our analyses provide a robust and comprehensive characterization of the Echinococcus germinative cell transcriptome, demonstrate expanded differentiation capability of metacestode derived stem cells, and underscore the potential of primary germinative cell cultures to investigate developmental processes of the parasite. These data are relevant for studies into the role of Echinococcus stem cells in parasite development and will facilitate the design of anti-parasitic drugs that specifically act on the parasite germinative cell compartment.

Cell repertoire and proliferation of germinative cells of the model cestode Mesocestoides corti

The phylum Platyhelminthes shares a unique population of undifferentiated cells responsible for the proliferation capacity needed for cell renewal, growth, tissue repair and regeneration. These cells have been extensively studied in free-living flatworms, whereas in cestodes the presence of a set of undifferentiated cells, known as germinative cells, has been demonstrated in classical morphology studies, but poorly characterized with molecular biology approaches. Furthermore, several genes have been identified as neoblast markers in free-living flatworms that deserve study in cestode models. Here, different cell types of the model cestode Mesocestoides corti were characterized, identifying differentiated and germinative cells. Muscle cells, tegumental cells, calcareous corpuscle precursor cells and excretory system cells were identified, all of which are non-proliferative, differentiated cell types. Besides those, germinative cells were identified as a population of small cells with proliferative capacity in vivo. Primary cell culture experiments in Dulbecco's Modified Eagle Medium (DMEM), Echinococcus hydatid fluid and hepatocyte conditioned media in non-reductive or reductive conditions confirmed that the germinative cells were the only ones with proliferative capacity. Since several genes have been identified as markers of undifferentiated neoblast cells in free-living flatworms, the expression of pumilio and pL10 genes was analysed by qPCR and in situ hybridization, showing that the expression of these genes was stronger in germinative cells but not restricted to this cell type. This study provides the first tools to analyse and further characterise undifferentiated cells in a model cestode.

Genome-wide transcriptome analysis of the early developmental stages of Echinococcus granulosus protoscoleces reveals extensive alternative splicing events in the spliceosome pathway

BACKGROUND

The complex life cycle of Echinococcus granulosus involves numerous environmental conditions within different intermediate and definitive hosts. This requires adaptation at different levels of transcript regulation. Alternative splicing (AS) and the related cellular functions as one of the major fields of post-genomics has been poorly studied in tapeworms. In the present study, we investigated AS events and their potential biological effects in E. granulosus.

METHODS

Whole transcriptome sequencing data of four groups of protoscoleces were prepared for RNA-seq library construction. Fresh protoscoleces were either used as non-induced controls (NT group) or incubated for 15 min with pepsin (PEP group) and cultivated in a biphasic medium for 12 and 24 h (12 and 24 h groups). The frequency and different types of AS events were identified using rMATS software. Functional annotations and gene ontology of differential AS (DAS) genes were performed using Blast2GO software. AS events were experimentally validated by PCR on the protoscolex cDNAs using specific primers for each gene.

RESULTS

At least one AS event was found in 38.1% of the genes (3904 out of 10,245) in the protoscoleces during early strobilar development. The genes were associated primarily with cellular and metabolic processes and binding and catalytic activity. KEGG pathway analysis of DAS events revealed a number of genes belonging to different components of the spliceosome complex. These genes tended to belong to common SR proteins, U1-related factors, U2-related factors, complex A-specific factors and other splicing-related proteins.

CONCLUSIONS

The high number of AS events in the transcriptome regulatory mechanisms indicates the essential rapid molecular processes required by the parasite for adaptation in different environments.

Identification and culture of proliferative cells in abnormal Taenia solium larvae: Role in the development of racemose neurocysticercosis

Racemose neurocysticercosis is an aggressive disease caused by the aberrant expansion of the cyst form of Taenia solium within the subarachnoid spaces of the human brain and spinal cord resulting in a mass effect and chronic inflammation. Although expansion is likely caused by the proliferation and growth of the parasite bladder wall, there is little direct evidence of the mechanisms that underlie these processes. Since the development and growth of cysts in related cestodes involves totipotential germinative cells, we hypothesized that the expansive growth of the racemose larvae is organized and maintained by germinative cells. Here, we identified proliferative cells expressing the serine/threonine-protein kinase plk1 by in situ hybridization. Proliferative cells were present within the bladder wall of racemose form and absent from the homologous tissue surrounding the vesicular form. Cyst proliferation in the related model species Taenia crassiceps (ORF strain) occurs normally by budding from the cyst bladder wall and proliferative cells were concentrated within the growth buds. Cells isolated from bladder wall of racemose larvae were established in primary cell culture and insulin stimulated their proliferation in a dose-dependent manner. These findings indicate that the growth of racemose larvae is likely due to abnormal cell proliferation. The different distribution of proliferative cells in the racemose larvae and their sensitivity to insulin may reflect significant changes at the cellular and molecular levels involved in their tumor-like growth. Parasite cell cultures offer a powerful tool to characterize the nature and formation of the racemose form, understand the developmental biology of T. solium, and to identify new effective drugs for treatment.

RNA interference dynamics in juvenile Fasciola hepatica are altered during in vitro growth and development

For over a decade RNA interference (RNAi) has been an important molecular tool for functional genomics studies in parasitic flatworms. Despite this, our understanding of RNAi dynamics in many flatworm parasites, such as the temperate liver fluke (Fasciola hepatica), remains rudimentary. The ability to maintain developing juvenile fluke in vitro provides the opportunity to perform functional studies during development of the key pathogenic life stage. Here, we investigate the RNAi competence of developing juvenile liver fluke. Firstly, all life stages examined possess, and express, core candidate RNAi effectors encouraging the hypothesis that all life stages of F. hepatica are RNAi competent. RNAi effector analyses supported growing evidence that parasitic flatworms have evolved a separate clade of RNAi effectors with unknown function. Secondly, we assessed the impact of growth/development during in vitro culture on RNAi in F. hepatica juveniles and found that during the first week post-excystment liver fluke juveniles exhibit quantitatively lower RNAi mediated transcript knockdown when maintained in growth inducing media. This did not appear to occur in older in vitro juveniles, suggesting that rapidly shifting transcript dynamics over the first week following excystment alters RNAi efficacy after a single 24 h exposure to double stranded (ds)RNA. Finally, RNAi efficiency was found to be improved through use of a repeated dsRNA exposure methodology that has facilitated silencing of genes in a range of tissues, thereby increasing the utility of RNAi as a functional genomics tool in F. hepatica.

Terminal-Repeat Retrotransposons in Miniature (TRIMs) in bivalves

Terminal repeat retrotransposons in miniature (TRIMs) are small non-autonomous LTR retrotransposons consisting of two terminal direct repeats surrounding a short internal domain. The detection and characterization of these elements has been mainly limited to plants. Here we present the first finding of a TRIM element in bivalves, and among the first known in the kingdom Animalia. Class Bivalvia has high ecological and commercial importance in marine ecosystems and aquaculture, and, in recent years, an increasing number of genomic studies has addressed to these organisms. We have identified biv-TRIM in several bivalve species: Donax trunculus, Ruditapes decussatus, R. philippinarum, Venerupis corrugata, Polititapes rhomboides, Venus verrucosa, Dosinia exoleta, Glycymeris glycymeris, Cerastoderma edule, Magallana gigas, Mytilus galloprovincialis. biv-TRIM has several characteristics typical for this group of elements, exhibiting different variations. In addition to canonically structured elements, solo-TDRs and tandem repeats were detected. The presence of this element in the genome of each species is <1%. The phylogenetic analysis showed a complex clustering pattern of biv-TRIM elements, and indicates the involvement of horizontal transfer in the spreading of this element.

Evidence for densovirus integrations into tapeworm genomes

BACKGROUND

Tapeworms lack a canonical piRNA-pathway, raising the question of how they can silence existing mobile genetic elements (MGE). Investigation towards the underlying mechanisms requires information on tapeworm transposons which is, however, presently scarce.

METHODS

The presence of densovirus-related sequences in tapeworm genomes was studied by bioinformatic approaches. Available RNA-Seq datasets were mapped against the Echinococcus multilocularis genome to calculate expression levels of densovirus-related genes. Transcription of densovirus loci was further analyzed by sequencing and RT-qPCR.

RESULTS

We herein provide evidence for the presence of densovirus-related elements in a variety of tapeworm genomes. In the high-quality genome of E. multilocularis we identified more than 20 individual densovirus integration loci which contain the information for non-structural and structural virus proteins. The majority of densovirus loci are present as head-to-tail concatemers in isolated repeat containing regions of the genome. In some cases, unique densovirus loci have integrated close to histone gene clusters. We show that some of the densovirus loci of E. multilocularis are actively transcribed, whereas the majority are transcriptionally silent. RT-qPCR data further indicate that densovirus expression mainly occurs in the E. multilocularis stem cell population, which probably forms the germline of this organism. Sequences similar to the non-structural densovirus genes present in E. multilocularis were also identified in the genomes of E. canadensis, E. granulosus, Hydatigera taeniaeformis, Hymenolepis diminuta, Hymenolepis microstoma, Hymenolepis nana, Taenia asiatica, Taenia multiceps, Taenia saginata and Taenia solium.

CONCLUSIONS

Our data indicate that densovirus integration has occurred in many tapeworm species. This is the first report on widespread integration of DNA viruses into cestode genomes. Since only few densovirus integration sites were transcriptionally active in E. multilocularis, our data are relevant for future studies into gene silencing mechanisms in tapeworms. Furthermore, they indicate that densovirus-based vectors might be suitable tools for genetic manipulation of cestodes.

Unexpected encounter of the parasitic kind

Both parasitology and stem cell research are important disciplines in their own right. Parasites are a real threat to human health causing a broad spectrum of diseases and significant annual rates morbidity and mortality globally. Stem cell research, on the other hand, focuses on the potential for regenerative medicine for a range of diseases including cancer and regenerative therapies. Though these two topics might appear distant, there are some “unexpected encounters”. In this review, we summarise the various links between parasites and stem cells. First, we discuss how parasites’ own stem cells represent interesting models of regeneration that can be translated to human stem cell regeneration. Second, we explore the interactions between parasites and host stem cells during the course of infection. Third, we investigate from a clinical perspective, how stem cell regeneration can be exploited to help circumvent the damage induced by parasitic infection and its potential to serve as treatment options for parasitic diseases in the future. Finally, we discuss the importance of screening for pathogens during organ transplantation by presenting some clinical cases of parasitic infection following stem cell therapy.

Region-specific regulation of stem cell-driven regeneration in tapeworms

Tapeworms grow at rates rivaling the fastest-growing metazoan tissues. To propagate they shed large parts of their body; to replace these lost tissues they regenerate proglottids (segments) as part of normal homeostasis. Their remarkable growth and regeneration are fueled by adult somatic stem cells that have yet to be characterized molecularly. Using the rat intestinal tapeworm, Hymenolepis diminuta, we find that regenerative potential is regionally limited to the neck, where head-dependent extrinsic signals create a permissive microenvironment for stem cell-driven regeneration. Using transcriptomic analyses and RNA interference, we characterize and functionally validate regulators of tapeworm growth and regeneration. We find no evidence that stem cells are restricted to the regeneration-competent neck. Instead, lethally irradiated tapeworms can be rescued when cells from either regeneration-competent or regeneration-incompetent regions are transplanted into the neck. Together, the head and neck tissues provide extrinsic cues that regulate stem cells, enabling region-specific regeneration in this parasite.

Many worms have remarkable abilities to regrow and repair their bodies. The parasitic tapeworms, for example, can reach lengths of several meters and grow much more quickly than tissues in humans and other complex animals. This growth allows tapeworms to counteract the continual loss of the segments that make up their bodies, known as proglottids – a process that happens throughout their lives.

The capacity to regenerate thousands of lost body segments and maintain an overall body length suggests that tapeworms have groups of stem cells in their body which can grow and divide to produce the new body parts. Yet, regeneration in tapeworms has not been closely studied.

Rozario et al. have now examined Hymenolepsis diminuta, the rat tapeworm, and identified the neck of the tapeworm as crucial for its ability to regrow lost body segments. Further analysis identified two genes, zmym3 and pogzl, that are essential for cell division during tapeworm growth. However, Rozario et al. showed that these genes are active elsewhere in the worm’s body and that it is the conditions found specifically in the tapeworm’s neck that create the right environment for stem cells to enable regeneration of new segments.

Tapeworms provide a valuable example for studying the growth of stem cells and these findings highlight the important role that the cells’ surroundings play in driving stem cell activity. These findings could also lead to new insights into how stem cells behave in other animals and could potentially lead to new approaches to prevent or treat tapeworm infections.

The role of fibroblast growth factor signalling in Echinococcus multilocularis development and host-parasite interaction

BACKGROUND

Alveolar echinococcosis (AE) is a lethal zoonosis caused by the metacestode larva of the tapeworm Echinococcus multilocularis. The infection is characterized by tumour-like growth of the metacestode within the host liver, leading to extensive fibrosis and organ-failure. The molecular mechanisms of parasite organ tropism towards the liver and influences of liver cytokines and hormones on parasite development are little studied to date.

METHODOLOGY/PRINCIPAL FINDINGS

We show that the E. multilocularis larval stage expresses three members of the fibroblast growth factor (FGF) receptor family with homology to human FGF receptors. Using the Xenopus expression system we demonstrate that all three Echinococcus FGF receptors are activated in response to human acidic and basic FGF, which are present in the liver. In all three cases, activation could be prevented by addition of the tyrosine kinase (TK) inhibitor BIBF 1120, which is used to treat human cancer. At physiological concentrations, acidic and basic FGF significantly stimulated the formation of metacestode vesicles from parasite stem cells in vitro and supported metacestode growth. Furthermore, the parasite's mitogen activated protein kinase signalling system was stimulated upon addition of human FGF. The survival of metacestode vesicles and parasite stem cells were drastically affected in vitro in the presence of BIBF 1120.

CONCLUSIONS/SIGNIFICANCE

Our data indicate that mammalian FGF, which is present in the liver and upregulated during fibrosis, supports the establishment of the Echinococcus metacestode during AE by acting on an evolutionarily conserved parasite FGF signalling system. These data are valuable for understanding molecular mechanisms of organ tropism and host-parasite interaction in AE. Furthermore, our data indicate that the parasite's FGF signalling systems are promising targets for the development of novel drugs against AE.

Identification of EmSOX2, a member of the Sox family of transcription factors, as a potential regulator of Echinococcus multilocularis germinative cells

Larvae of the tapeworm Echinococcus multilocularis cause alveolar echinococcosis (AE), one of the most lethal helminthic infections in humans. The germinative cells, a population of stem cell-like cells, are considered to drive the continuous growth of the metacestodes within the host. The mechanisms and relative molecules controlling the behavior of germinative cells are poorly understood. Sox transcription factors play important roles in maintenance and regulation of stem/progenitor cells. We here describe the identification of a Sox family member in E. multilocularis, EmSOX2, as a potential regulator of germinative cells. Replacement of mouse Sox2 with EmSox2 could derive induced pluripotent stem cells (iPSCs) from somatic cells, suggesting that EmSOX2 is functionally related to mammalian SOX2. EmSOX2 is actively expressed in the proliferating germinative cells in E. multilocularis, and is significantly downregulated upon specific depletion of the germinative cell population by hydroxyurea treatment. These findings suggest that EmSOX2 may play a critical role in regulating E. multilocularis germinative cells.

A TRIM insertion in the promoter of Ms2 causes male sterility in wheat

The male-sterile ms2 mutant has been known for 40 years and has become extremely important in the commercial production of wheat. However, the gene responsible for this phenotype has remained unknown. Here we report the map-based cloning of the Ms2 gene. The Ms2 locus is remarkable in several ways that have implications in basic biology. Beyond having no functional annotation, barely detectable transcription in fertile wild-type wheat plants, and accumulated destructive mutations in Ms2 orthologs, the Ms2 allele in the ms2 mutant has acquired a terminal-repeat retrotransposon in miniature (TRIM) element in its promoter. This TRIM element is responsible for the anther-specific Ms2 activation that confers male sterility. The identification of Ms2 not only unravels the genetic basis of a historically important breeding trait, but also shows an example of how a TRIM element insertion near a gene can contribute to genetic novelty and phenotypic plasticity.

The Echinococcus canadensis (G7) genome: a key knowledge of parasitic platyhelminth human diseases

BACKGROUND

The parasite Echinococcus canadensis (G7) (phylum Platyhelminthes, class Cestoda) is one of the causative agents of echinococcosis. Echinococcosis is a worldwide chronic zoonosis affecting humans as well as domestic and wild mammals, which has been reported as a prioritized neglected disease by the World Health Organisation. No genomic data, comparative genomic analyses or efficient therapeutic and diagnostic tools are available for this severe disease. The information presented in this study will help to understand the peculiar biological characters and to design species-specific control tools.

RESULTS

We sequenced, assembled and annotated the 115-Mb genome of E. canadensis (G7). Comparative genomic analyses using whole genome data of three Echinococcus species not only confirmed the status of E. canadensis (G7) as a separate species but also demonstrated a high nucleotide sequences divergence in relation to E. granulosus (G1). The E. canadensis (G7) genome contains 11,449 genes with a core set of 881 orthologs shared among five cestode species. Comparative genomics revealed that there are more single nucleotide polymorphisms (SNPs) between E. canadensis (G7) and E. granulosus (G1) than between E. canadensis (G7) and E. multilocularis. This result was unexpected since E. canadensis (G7) and E. granulosus (G1) were considered to belong to the species complex E. granulosus sensu lato. We described SNPs in known drug targets and metabolism genes in the E. canadensis (G7) genome. Regarding gene regulation, we analysed three particular features: CpG island distribution along the three Echinococcus genomes, DNA methylation system and small RNA pathway. The results suggest the occurrence of yet unknown gene regulation mechanisms in Echinococcus.

CONCLUSIONS

This is the first work that addresses Echinococcus comparative genomics. The resources presented here will promote the study of mechanisms of parasite development as well as new tools for drug discovery. The availability of a high-quality genome assembly is critical for fully exploring the biology of a pathogenic organism. The E. canadensis (G7) genome presented in this study provides a unique opportunity to address the genetic diversity among the genus Echinococcus and its particular developmental features. At present, there is no unequivocal taxonomic classification of Echinococcus species; however, the genome-wide SNPs analysis performed here revealed the phylogenetic distance among these three Echinococcus species. Additional cestode genomes need to be sequenced to be able to resolve their phylogeny.

Echinococcus-Host Interactions at Cellular and Molecular Levels

The potentially lethal zoonotic diseases alveolar and cystic echinococcosis are caused by the metacestode larval stages of the tapeworms Echinococcus multilocularis and Echinococcus granulosus, respectively. In both cases, metacestode growth and proliferation occurs within the inner organs of mammalian hosts, which is associated with complex molecular host-parasite interactions that regulate nutrient uptake by the parasite as well as metacestode persistence and development. Using in vitro cultivation systems for parasite larvae, and informed by recently released, comprehensive genome and transcriptome data for both parasites, these molecular host-parasite interactions have been subject to significant research during recent years. In this review, we discuss progress in this field, with emphasis on parasite development and proliferation. We review host-parasite interaction mechanisms that occur early during an infection, when the invading oncosphere stage undergoes a metamorphosis towards the metacestode, and outline the decisive role of parasite stem cells during this process. We also discuss special features of metacestode morphology, and how this parasite stage takes up nutrients from the host, utilizing newly evolved or expanded gene families. We comprehensively review mechanisms of host-parasite cross-communication via evolutionarily conserved signalling systems and how the parasite signalling systems might be exploited for the development of novel chemotherapeutics. Finally, we point to an urgent need for the development of functional genomic techniques in this parasite, which will be imperative for hypothesis-driven analyses into Echinococcus stem cell biology, developmental mechanisms and immunomodulatory activities, which are all highly relevant for the development of anti-infective measures.

Evolutionary developmental biology (evo-devo) of cestodes

Cestodes (tapeworms) have complex adaptations to their obligatory parasitic life-style. Among these adaptations, they show many evolutionary innovations in their development, including complex life-cycles with multiple hosts and life-stages, several independent origins of asexual reproduction, and the evolution of segmentation as a mean to generate massive reproductive output. Therefore, cestodes offer many opportunities for the investigation of the evolutionary origins of developmental novelties (evo-devo). However, cestodes have not been exploited as major models for evo-devo research due to the considerable technical difficulties involved in their study. In this review, a panoramic view is given of classical aspects, methods and hypothesis of cestode development, together with recent advances in phylogenetics, genomics, culture methods, and comparative analysis of cestode gene expression. Together with the availability of powerful models for related free-living flatworms, these developments should encourage the incorporation of these fascinating parasites into the first-line of evo-devo research.

Evolutionary characterization of Ty3/gypsy-like LTR retrotransposons in the parasitic cestode Echinococcus granulosus

Cyclophyllidean cestodes including Echinococcus granulosus have a smaller genome and show characteristics such as loss of the gut, a segmented body plan, and accelerated growth rate in hosts compared with other tissue-invading helminths. In an effort to address the molecular mechanism relevant to genome shrinkage, the evolutionary status of long-terminal-repeat (LTR) retrotransposons, which are known as the most potent genomic modulators, was investigated in the E. granulosus draft genome. A majority of the E. granulosus LTR retrotransposons were classified into a novel characteristic clade, named Saci-2, of the Ty3/gypsy family, while the remaining elements belonged to the CsRn1 clade of identical family. Their nucleotide sequences were heavily corrupted by frequent base substitutions and segmental losses. The ceased mobile activity of the major retrotransposons and the following intrinsic DNA loss in their inactive progenies might have contributed to decrease in genome size. Apart from the degenerate copies, a gag gene originating from a CsRn1-like element exhibited substantial evidences suggesting its domestication including a preserved coding profile and transcriptional activity, the presence of syntenic orthologues in cestodes, and selective pressure acting on the gene. To my knowledge, the endogenized gag gene is reported for the first time in invertebrates, though its biological function remains elusive.