Raúl J. Bobes1†
Karel Estrada2†
Diana G. Rios-Valencia1†
Arturo Calderón-Gallegos1
Patricia de la Torre1Julio C. Carrero1
Alejandro Sanchez-Flores2*‡
Juan P. Laclette1*‡- 1Biomedical Research Institute, Universidad Nacional Autónoma de México, CDMX, Mexico
- 2Biotechnology Institute, Universidad Nacional Autónoma de México, Cuernavaca, Mexico
Human cysticercosis by Taenia solium is the major cause of neurological illness in countries of Africa, Southeast Asia, and the Americas. Publication of four cestode genomes (T. solium, Echinococcus multilocularis, E. granulosus and Hymenolepis microstoma) in the last decade, marked the advent of novel approaches on the study of the host-parasite molecular crosstalk for cestode parasites of importance for human and animal health. Taenia crassiceps is another cestode parasite, closely related to T. solium, which has been used in numerous studies as an animal model for human cysticercosis. Therefore, characterization of the T. crassiceps genome will also contribute to the understanding of the human infection. Here, we report the genome of T. crassiceps WFU strain, reconstructed to a noncontiguous finished resolution and performed a genomic and differential expression comparison analysis against ORF strain. Both strain genomes were sequenced using Oxford Nanopore (MinION) and Illumina technologies, achieving high quality assemblies of about 107 Mb for both strains. Dotplot comparison between WFU and ORF demonstrated that both genomes were extremely similar. Additionally, karyotyping results for both strains failed to demonstrate a difference in chromosome composition. Therefore, our results strongly support the concept that the absence of scolex in the ORF strain of T. crassiceps was not the result of a chromosomal loss as proposed elsewhere. Instead, it appears to be the result of subtle and extensive differences in the regulation of gene expression. Analysis of variants between the two strains identified 2,487 sites with changes distributed in 31 of 65 scaffolds. The differential expression analysis revealed that genes related to development and morphogenesis in the ORF strain might be involved in the lack of scolex formation.
Introduction
The tapeworm family Taeniidae has been extensively studied, particularly because of the importance of Taenia solium, Taenia asiatica (Wang et al., 2016; Gripper and Welburn, 2017), Echinococcus granulosus and Echinococcus multilocularis (Battelli, 2009; Deplazes et al., 2017) as public health threats. Other species also have veterinary importance such as Taenia saginata and Taenia pisiformis (Silva and Costa-Cruz, 2012; Yang et al., 2012).
T. solium is the cestode parasite that generates a more significant economic and health burden on humans and pigs, this parasite, is responsible of porcine cysticercosis, human teniasis and neurocysticercosis, considered the most frequent parasite disease of the central nervous system (Garcia and Modi, 2008; Fleury et al., 2010; Carpio et al., 2018), and a neglected tropical disease by the World Health Organization (World Health Organization, 2011). Humans acquire cysticercosis through fecal-oral contamination with T. solium eggs from human tapeworm carriers. The most common places for the establishment of the larval form in the intermediate host are the skeletal muscles, the ocular system and the central nervous system. Cysticercosis is of great economic and public health importance, causing morbidity and mortality in countries and regions with deficient sanitary conditions. Many of these developing countries are located in Southeast Asia, Africa and the Americas (Ito et al., 2003; Flisser et al., 2006; Lustigman et al., 2012; Zammarchi et al., 2013; Ito, 2015). In the last decades, cysticercosis has spread to Europe, the United States and Australia, mainly due to the increase in human migration from endemic regions (Gabriël et al., 2015; O’Neal and Flecker, 2015). T. solium has been globally classified as the main foodborne parasite by the United Nations Food and Agriculture Organization (FAO) in 2014 and the World Health Organization (Torgerson et al., 2015).
Taenia crassiceps, a close relative of T. solium, can be easily maintained in most animal facilities and the cysticerci can be manipulated in vivo and in vitro, providing a suitable model for the study of several aspects of human cysticercosis (Willms and Zurabian, 2010). Thus, the intraperitoneal murine infection by T. crassiceps has been widely used as a model of cysticercosis in studies of the host immune response, diagnosis, evaluation of anthelmintic drugs, vaccine development (Larralde et al., 1990; Toledo et al., 2001; Vargas-Villavicencio et al., 2005; Sciutto et al., 2011; Mahanty et al., 2013); sex and genetic influence on the susceptibility to infection (Sciutto et al., 1991), influence of the parasite on the hormonal status of the infected host (Morales-Montor et al., 2008), development of practical methods for transfection of exogenous genes (Moguel et al., 2015), among other aspects. Besides, T. crassiceps infections in humans are very rare, isolated cases have been reported in immunosuppressed individuals (Heldwein et al., 2006; Willms and Zurabian, 2010; Lescano and Zunt, 2013; Deplazes et al., 2019). Human infections by T. crassiceps include subcutaneous and intramuscular (François et al., 1998; Maillard et al., 1998; Heldwein et al., 2006; Goesseringer et al., 2011), intraocular (Shea et al., 1973; Arocker-Mettinger et al., 1992; Mougeot et al., 1996; Chuck et al., 1997) and cerebellum 97 locations (Ntoukas et al., 2013).
As an animal model, the metacestode stage of T. crassiceps is usually maintained through serial intraperitoneal passages of cysts from the peritoneal cavity of infected to recipient naïve mice, where they reproduce asexually by budding (Freeman, 1962; Esch and Smyth, 1976; Good and Miller, 1976; Willms and Zurabian, 2010). Definitive hosts for T. crassiceps are foxes, groundhogs, wolves, and dogs, intermediate hosts include field rats, mice and other small rodents, (Loos-Frank, 2000).
Several strains of T. crassiceps have been reported, such as WFU (Everhart et al., 2004), HYG (Sally et al., 1976) or KBS (Dorais and Esch, 1969), isolated from natural infections. The cysticercus possesses an invaginated scolex, which evaginates when ingested by the definitive host and develops into the adult worm. In contrast, the ORF strain of T. crassiceps, isolated in 1962 (Freeman et al., 1962), lacks an invaginated scolex and has been considered a mutant strain. It has been proposed that ORF suffered a chromosomal loss that impairs scolex formation (Smith et al., 1972), thus producing sterile larvae that cannot develop into adult worms in the wild (Mount, 1968; Dorais and Esch, 1969); can only be propagated in the laboratory and are considered safe for human manipulation (Freeman, 1962).
In terms of genomic resources, the firsts genomes reported were for Schistosoma japonicum (Zhou et al., 2009), S. mansoni (Berriman et al., 2009). Since then, several genomes have been published including trematodes of medical importance such as, and Clonorchis sinensis (Wang et al., 2011; Huang et al., 2013). In 2013, the genomes of four cestodes: Hymenolepis microstoma, Echinococcus granulosus, E. multilocularis and Taenia solium were reported (Tsai et al., 2013) and a free living flatworm Schmidtea mediterranea (Grohme et al., 2018). Afterwards, the International Helminth Genomes Consortium started a project for genome characterization whose database (WormbaseParasite) already provides 45 Platyhelminthes and 19 cestodes (International Helminth Genomes Consortium, 2019).
In the last five years the genomes of other taeniidae have been reported: T. asiatica, T. saginata (Wang et al., 2016), T. multiceps (Li et al., 2018), E. canadensis (Maldonado et al., 2017), E. oligarthrus (Maldonado et al., 2019), H. diminuta (Nowak et al., 2019), H. microstoma (Olson et al., 2020). Availability of these genomes has increased the knowledge about gene numbers and composition of these organisms, improving our understanding of their biology and pathogenic mechanisms.
Here we report the high resolution characterization for the genome of T. crassiceps WFU strain and a comparison analysis against ORF strain. Both genomes were compared in order to find differences that could be related to the lack of scolex development. Our results provide evidence that the absence of scolex in the ORF strain was not the result of a chromosomal loss proposed elsewhere. Instead, it appears to be the result of subtle differences in the regulation of gene expression or unidentified mutations in specific genes.
Materials and Methods
Parasite Culture
Cysticerci from T. crassiceps strains WFU and ORF (Supplementary Figure 1) were propagated in the peritoneum of 4-week-old female BALB/c mice, which were inoculated intraperitoneally with 10 cysticerci and after 90 days were sacrificed for the recovery of the larvae. The cysticerci were thoroughly washed with sterile Phosphate Buffer Saline (PBS), pH 7.4, and then frozen at -70°C until use. The use and humanitarian handling of animals in this project was authorized by the Institutional Committee for the Care and Use of Laboratory Animals (CICUAL), UNAM with the register number 6329.
Chromosome Spreading for Taenia crassiceps Karyotyping
The protocol for spreading chromosomes was adapted from (Guo et al., 2018 and Špakulová et al., 2011). Cysticerci recovered from the peritoneum of BALB/C mice with three months of infection were selected and placed in 16-well culture plates with RPMI-1640 media supplemented with 10% FBS and colchicine (0.25% w/v) during 5 h at 37°C. After incubation, cysticerci were washed twice with PBS and then placed in deionized water for 20 min and then fixed with a Methanol-Acetic acid solution (3:1) for 30 minutes. The chromosomes were spread by crushing the cysts between two slides with strong physical pressure. The slides were then frozen with liquid nitrogen and washed twice with PBS before being permeabilized with PBS-Triton X-100 (1%) for 5 min. Finally, the spread was stained with DAPI (1:500) for 20 minutes. Slides were mounted in Fluoroshield (Sigma) and observed under an inverted fluorescence microscope (Olympus IX71). Image analysis was made using FIJI software.
DNA and RNA Isolation and Sequencing
Fresh cysticerci (both WFU and ORF strains) obtained from the peritoneum of mice were used for DNA and RNA extraction after three times washing in PBS, pH 7.4. A 1 mL of supernatant from the PBS recovered cysticerci was centrifuged at 12,000 g. The pellet was resuspended in 200 µL of Biofluid & Solid Tissue Buffer, added with 20 µL of Proteinase K and 2 µL de RNAsa (10 mg/mL), and incubated at 55°C for one hour. The lysis reaction was transferred using a bore-tip to a tissue grinder for 10 minutes and 1 mL of phenol was added and after 2 minutes, the solution was centrifuged at 12,000 g for another 2 minutes. A 350 µL volume of the supernatant was transferred to a clean tube and 1mL of chloroform was added. Chloroform addition and centrifugation was repeated twice and 250 µL of the supernatant were recovered and placed in 150 µL of RNAse-free water. DNA was extracted using the Zymo Quick-DNA™ HMW MagBead Kit (Catalog No. D6060) according to the vendor’s protocol. The same protocol was adapted for RNA extraction, skipping the 55°C incubation and harsher disruption using high-speed vortexing.
Different Illumina sequencing libraries were prepared using the extracted DNA from both WFU and ORF strains. Two Nextera XT DNA libraries, one with ~150 bp insert size and another with a ~500 bp insert size. Also, two Illumina Mate Paired libraries with 3-4 kb and 4-6 kb insert sizes, respectively. All libraries were prepared following the vendor’s protocols. For Oxford Nanopore sequencing, the DNA libraries for WFU and ORF DNA samples were prepared using the SQK-LSK109 using 1µg of DNA for each sample and following the vendor’s protocol. Sequencing was performed using a whole minION flow cell (v9.4) for each sample. Basecalling was done by the program Guppy v4.0.14 with default parameters and a high accuracy error model (dna_r9.4.1_450bps_hac.cfg).
RNA sequencing libraries were prepared using the Illumina TruSeq RNA Sample Prep Kit v2 (Illumina) following the manufacturer’s protocol. High quality RNA samples from WFU and ORF strains, three biological replicates for each one with RIN numbers ranging from 8 to 10 were used for NGS library construction. The libraries were sequenced using the Illumina NextSeq 500 platform with a paired-end configuration using 150 cycles (2x75 reads).
Genome Assembly, Post-Assembly Improvement, Gene Prediction and Annotation
A de novo genome assembly for WFU strain was performed using a hybrid approach with two different sequencing technologies. The first assembly was obtained using the AllPaths-LG v.52448 (Butler et al., 2008) assembler using default parameters, with Illumina paired-end and mate pair libraries. Then, a second de novo assembly was obtained using Oxford Nanopore long reads and Canu v.1.9 (Koren et al., 2017) assembler. This assembly was corrected using two error correction tools, first Racon v1.3.1 (Vaser et al., 2002) with 3 iterations and later 3 iterations of Pilon v.1.23 (Walker et al., 2014), both used the illumina libraries employed in the first assembly. Finally, both assemblies were merged using the Ntjoin software v.1.0.8 (Coombe et al., 2020). The latest version of the genome can be found with the accession number JAKROA000000000, associated with the BioProject PRJNA807072.
For gene prediction, we first used RepeatModeler v1.0.11 (Smit et al., 2016) and RepeatMasker v4.08 (Smit et al., 2018) to softmasking the genome. We used Braker2 v2.1.6 (Brůna et al., 2021) which used 30,000 orthologous proteins that were obtained from T. solium, E. multilocularis and E. granulosus genomes and 10 Gb of RNAseq data from T. crassiceps to train AUGUSTUS v3.3.2 (Stanke et al., 2006) which run within the Braker pipeline. Briefly, we generated hints by mapping the RNAseq data using HISAT2 v2.2.1 (Kim et al., 2019) to generate a sorted BAM. Then, using the BAM file as hints, AUGUSTUS was trained along with the orthologous proteins to generate the gene predictions for T. crassiceps. The protein annotation was achieved using an in-house pipeline with the following steps: first we used Blastp v2.11.0 (Camacho et al., 2009) to align these against the curated protein database Swiss-Prot and assign the function of the protein with the greatest similarity. We used HmmScan v3.3.2 (Eddy, 2009) to search protein domains using the Pfam-A database, additionally we used SignalP v4.1 (Petersen et al., 2011) to identify annotated signal peptides. Finally, we enriched the annotation with the GO and KEGG database, assigning the GO terms associated with each protein and their relationship with metabolic pathways, respectively.
SNPs Analysis
For variant calling and annotation, we used BWA v0.7.17 (Li and Durbin, 2009) to map the Illumina DNA sequencing of the ORF strain vs the WFU genome. Then we used Freebayes v1.2.0 (Garrison and Gabor, 2012) to analyze the BAM resulting from the previous alignment and generate the VCF variants file.
Annotation of the variants was performed using SnpEff v4.3t (Cingolani et al., 2012); the program takes the WFU annotation file in GFF format and its genome as a database and the resulting VCF from the variant call.
Comparative Genomics
For the comparative genomic analysis, the genomes available at Genbank for E. multilocularis (GCA_000469725.3), E. granulosus (GCA_000524195.1), E. canadiensis (GCA_900004735.1), T. solium (GCA_001870725.1) and H. microstoma (GCA_000469805.3) were used. Also included was the T. crassiceps genome here described (JAKROA000000000). The Proteinortho v6.02 (Lechner et al., 2011) was used to obtain protein clusters of WFU genomes analyzed here with the following parameters: E-value for blast: 1e-05, minimum percent identity of best blast alignments: 25, minimum coverage of best blast alignments in percent: 50, minimum similarity for additional hits: 0.95. The results from the protein orthologous was also used to verify the annotation. Additionally, we used the Circos v0.69-8 (Krzywinski et al., 2009) package for circular data visualization.
Differential Expression Analysis
We used Bowtie2 v2.3.4.3 (Langmead and Salzberg, 2012) to map the sequences from WFU and ORF RNAseq libraries to the coding sequences (CDS) of WFU strain. Subsequently, the counting matrices were obtained with Express v1.5.1 (Roberts et al., 2011) and finally the EdgeR v3.24.3 (Robinson et al., 2010) package was used to obtain the results of the differential expression analysis. We used a FDR cutoff value equal or less than 0.05 and Log-Fold Change value equal or greater than 1 to consider a gene as differentially expressed.
Gene Ontology Analysis
We used the R package TopGO v2.46.0 (Alexa et al., 2006) to 256 perform enrichment analysis of the differentially expressed gene set and of those proteins that were unique to T. crassiceps in the comparison between orthologues.
Results
Genome Assembly and Curation
The complete genome for the T. crassiceps WFU strain was sequenced using genomic DNA obtained from cysts and reconstructed through a hybrid approach using Illumina and Nanopore technologies. The initial assembly produced 179 scaffolds for a total of 118.7 Mb. However, several of the shorter scaffolds presented repetitive patterns or high similarity to regions in larger scaffolds. After performing the comparison analysis at sequence level, we confirmed that 114 scaffolds belonged to telomeric or collapsed regions. Therefore, the curated final assembly we obtained included only 65 scaffolds for a total of 107.05 Mb. The statistics of the two genomic assemblies are shown in Table 1.
Table 1 T. crassiceps assemblies and genomic statistics.
Noncontiguous Finished Resolution for T. crassiceps Genome and Comparative Genomics
We compared the ten largest scaffolds in the curated assembly of T. crassiceps with data from genomes of two closely related species with chromosome resolution assemblies. A CIRCOS plot comparison between H. microstoma, E. multilocularis and T. crassiceps is shown in Figure 1A. Three T. crassiceps scaffolds presented a similar degree of completeness and synteny to E. multilocularis chromosomes: scaffolds 1-3 of T. crassiceps vs. scaffolds 2, 4 and 3 from E. multilocularis, respectively. This suggested that these scaffolds might correspond to whole chromosomes. However, none of the compared T. crassiceps scaffolds showed a clear correspondence to H. microstoma scaffolds demonstrating a larger phylogenetic distance.
Figure 1 Comparative genomics between strains of Taenia crassiceps and with different cestode species. (A) Comparative genomics of three cestode genomes [Hymenolepis microstoma (HMN); Echinococcus multilocularis (EmW) and Taenia crassiceps WFU (TcrWFU)]. Scaffold sequences used for this comparison are available in the NCBI and WormBase databases. Genomic fragment names are those from the original source annotation (B) Dotplot comparison of WFU vs ORF strains at genomic level demonstrated a clear synteny discarding large structural variations.
An extensive comparison of the identified proteins was carried out with five other species of cestodes to analyze their orthology with T. crassiceps. The proteomes of four taeniid species (T. solium, T. crassiceps, E. multilocularis, E. granulosus and E. canadiensis), as well as H. microstoma were compared (Supplementary File 1). A total of 4,665 orthologue groups were found in all species (cyclophyllidea core-genome); 5,221 were found in T. solium, T. crassiceps, E. multilocularis, E. granulosus and E. canadiensis (Taenidae core genome); 6,842 orthologue groups were shared between T. solium and T. crassiceps (Taenia core genome). Finally, 2,822 proteins were unique to T. crassiceps and did not cluster with other proteins in any other genome. Among these, 2,248 didn’t have a functional annotation.
Genomic Variation Between T. crassiceps WFU and ORF Strains
Using the curated assembly of T. crassiceps WFU genome as a template, we evaluated the differences of the ORF strain. A dotplot analysis indicated that both assemblies were extremely close (Figure 1B); no differences among them were found that might suggest that a catastrophic genomic event happened in the ORF strain, such as the loss of a complete chromosome. As for the differences in both strains, using short read mapping with the Illumina sequencing, we found 2,487 site variants in ORF, distributed in 31 of the 65 scaffolds in the WFU curated assembly. Of these, 155 variants affected 147 protein coding genes resulting in missense or frameshift mutations. Variants were mainly single-base substitutions, insertions or deletions. Since long read sequencing was also done for the ORF strain, we searched for larger structural changes but none was found. All the information regarding the variants between the two strains can be found in the Supplementary File 1.
WFU and ORF Karyotype Determination
It has been proposed that the lack of scolex in the ORF strain was due to an aneuploidy: a loss of a pair of chromosomes. Therefore, we also carried out a chromosome spreading to karyotype the metaphase nuclei of both strains. Interestingly the number 2n = 16 or 18 appeared with similar frequencies. This observation has been reported elsewhere (Špakulová et al., 2011). Whatever the number of chromosome pairs might be, we were unable to observe differences in both strains (Supplementary Figure 2). It is worth mentioning that the morphologic resolution we achieved did not allow us to observe if differences in the chromosome structure existed.
Differential Expression Analysis Between ORF and WFU Strains of T. crassiceps
To elucidate the molecular basis of the phenotypic differences between WFU and ORF in relation with the lack of scolex development, we performed a differential expression analysis. A total of 1080 genes were found differentially expressed in the ORF strain. Among them, 627 and 453 genes were up and down-regulated, respectively (Supplementary File 1). Using a GO term enrichment analysis of down-regulated genes with |logFc| >2 as threshold (having functional annotation in ORF), different protein families were identified. Interestingly, some of them were related to embryonic digestive tract morphogenesis, mesenchymal cell differentiation, maintenance of animal organ identity, establishment of protein localization, among others (Supplementary File 1). Differentially expressed genes were also characterized by GO revealing a particular enrichment for both up and down-regulated genes related to integral membrane proteins in the cellular component category (Figure 2). In the same category, a considerable number of genes related to cytoskeleton were down-regulated. For up-regulated genes (Figure 2), in the biological processes category, most of the differentially expressed genes were grouped within response to stimulus. In total, 45 UF genes were also down regulated in ORF (Supplementary File 1).
Figure 2 Gene ontology enrichment analysis in the T. crassiceps ORF strain with |logFc| >2 as threshold.
Differentially Expressed Genes and Genomic Variation
Analysis of genes with different expression levels whose coding sequence contained SNPs allowed the identification of 17 genes, including 5 of unknown function. Among the annotated genes, 3 were up-expressed and 7 were down-expressed in ORF, including some related with development, such as homeobox IRX-5 (TcrWFU_07074), rolling stone [rost] (TcrWFU_02139), as well as genes participating in the cell cycle: ubiquitin carboxyl terminal hydrolase 16 (TcrWFU_00103) and others involved in signal transduction pathways: Axin-1 (TcrWFU_07383) that participates in the Wnt/ß catenin pathway (Supplementary File 1).
Discussion
The present work describes the characterized genomes and RNAseq results of two T. crassiceps strains (WFU and ORF), that are widely used as animal models for the study of human cysticercosis. The use of T. crassiceps as an adequate animal model seems justified in the light of our results of orthologue proteins in both species shown here; 6,842 orthologue groups were shared with T. solium demonstrating a high evolutionary closeness with T. crassiceps.
As mentioned above, the ORF strain lost the capacity to develop a scolex, structure that gives rise to the adult worm. It has been proposed that the absence of the scolex in the ORF strain was due to the loss of a pair of chromosomes (Smith et al., 1972). We undertook a genomic and transcriptomic characterization of these two strains as a way to approach this proposal. Initially, a dotplot analysis using the two assemblies showed minimal differences between both strains. Moreover, our karyotyping results did not reveal differences between the chromosome numbers in both strains, such as the absence of a whole chromosome pair (aneuploidy) reported previously (Smith et al., 1972). Even without a full chromosome resolution, the genomic information from both strains revealed no differences or massive structural changes that might support the loss of two chromosomes. Only minor differences were found in the genome drafts of both strains, including gene number, GC content and SNPs number. Instead, we hypothesize that the loss of scolex is the result of changes in the expression of genes involved in the development of this crucial structure in the ORF strain. In order to identify the network of genes resulting in the lack of scolex, we decided to analyze differences in gene expression (up or down-regulated) in the ORF strain. Considering that the scolex is the structure that gives rise to the adult stage of the parasite, we explored a diversity of genes related to embryonic development, formation of the anteroposterior axis, cell division, establishment of mitotic spindle orientation, among others.
The genome of T. crassiceps WFU strain was reconstructed using long and short read sequencing technologies. This allowed us to generate a noncontiguous high resolution genome, which was compared to the ORF strain and other cestode genomes with similar or higher resolution. According to the BUSCO completeness analysis, the reconstructed genome had a 78.5% completeness which is similar to the E. multilocularis genome (GCA_000469725.3). It is expected that genomes from parasitic species have a lower completeness due to a lack of proteins that are present in other free living eukaryotes. Interestingly, the ten largest fragments have a size similar to chromosomes in other species and five of them had a similar arrangement to chromosomes in E. multilocularis, which is a manually curated genome and a closely related species to T. crassiceps; both are included within the same family (Taenidae). Three T. crassiceps chromosomes (scaffolds 1-3) presented a very similar synteny to E. multilocularis but some of the comparisons suggested that other large fragments of T. crassiceps can also be arranged on the structure of other chromosomes in E. multilocularis (Figure 1A), although rearrangements cannot be ruled out. Therefore, other approaches like optical mapping can be used to resolve the assembly contiguity in the near future. Unfortunately, the current fragmentation of the T. solium genome, didn’t allow us to include it in the comparison. Our orthology analysis among cyclophyllideans, taeniids and Taenia genomes is of particular interest; it clearly shows increasing orthology within closer clades of tapeworms, supporting the basis of current taxonomy for this group of parasites. T. solium and T. crassiceps shared 6,842 orthologue groups, whereas 2,822 proteins were unique proteins for T. crassiceps; including 2,248 without functional annotation, not significant similarity with sequences in public databases. These unannotated proteins deserve further exploration to elucidate their role in this parasite. Finally, 574 unique proteins in T. crassiceps genome were included in a GO term enriched analysis (Supplementary File 1). From those unique proteins for T. crassiceps, we found 60 that were differentially expressed in ORF with respect to WFU. Interestingly, proteins without orthologues in other cestodes were related to molecular functions/biological processes such as arginyltransferase activity/protein arginylation; peptidase inhibitor activity/negative regulation of peptidase activity and glutathione peroxidase activity/response to oxidative stress. These ontologies suggest that some of the proteins may be involved in protease degradation/blocking or oxidative stress response. However, these and other results require further investigation and experimental assessment to confirm the role of those unique proteins in T. crassiceps.
From the total of 1080 genes, 843 had annotated function; our RNA-seq results indicated that 627 were upregulated and 453 genes were downregulated in ORF, in comparison to the WFU strain, considered as the wild-type strain in this work. Those changes in expression included zinc finger proteins, transcription factors, homeoboxes, components of the Wnt pathway, etc., that had been identified in previous studies on other cestode species (Olson et al., 2018).
As an example, in the ORF strain, down-regulated genes included homeobox B7 (HoxB7) (TcrWFU_04386), which is mainly related with processes of cell differentiation and segmentation in the adult stage; since they are presented in protoscoleces and in adult worms of E. granulosus (Dezaki et al., 2016). HoxB7 420 participates in processes such as cell proliferation and differentiation including 421 anterior/posterior pattern formation (Carè et al., 1999). Moreover, an important role of HoxB7 is related with morphological diversification of various body structures during development (Wellik, 2009). Regarding HoxB7, a differential expression has been demonstrated at different stages of development and structures in E. granulosus (in vivo and in vitro). Down-regulation of HoxB7 in ORF could be related to the absence of the scolex.
Another gene of interest related to development is frizzeld. Our results showed that frizzled (TcrWFU_07188) was down-regulated in the ORF strain. The BLAST result shows 100 and 89% identity with frizzeld-10 and frizzeld-4 of E. granulosus, respectively. Regarding the Wnt ligand, the frizzeld-4 gene (em-fz4) is expressed in dispersed cells in the germinal layer, becomes strongly up-regulated during brood capsule development, and is always restricted to the posterior-most region of the protoscolex throughout development in H. microstoma and E. multilocularis (Koziol et al., 2016). This finding in closely related organisms could explain why the gene is down-regulated, impairing formation of the scolex in T. crassiceps ORF strain.
On the other hand, down-regulated in the ORF strain (in comparison to WFU) whose function could be related to its inability to form a scolex, were seven genes involved in embryonic development and morphogenesis (Supplementary File 1). Among them, transcription factor 21 (TcrWFU_01181), a member of a subfamily of basic-helix-loop-helix proteins (bHLH) showed a high differential expression (LogFC = -2.02). This factor has been reported to play an important role in embryonic development of mesodermal tissues (Quaggin et al., 1998). Since the mesoderm gives rise to the notochord and the neural tube, it is also involved in the definition of the anteroposterior axis (Corallo et al., 2015). Another down-regulated protein (LogFC = - 2.23) in ORF is the Tyrosine protein-kinase src-1 (TcrWFU_10053), which together with Wnt, plays an essential role in embryonic development contributing to spindle orientation and endoderm specification in early development (Bei et al., 2002). Noteworthy, this gene is exclusively expressed in some head neurons of the nematode Caenoharbditis elegans, but not in its tail (Hirose et al., 2003). Finally, the other down-regulated gene in ORF is the fibroblast growth factor receptor-like 1 [(LogFC = -2.08) (TcrWFU_04317)], involved in the development of organisms and expressed in the anterior region (Lamb and Harland, 1995; Teven et al., 2014). This gene has been identified as up-regulated in larval and adult H. microstoma (HmN_000423900) (Olson et al., 2018).
Regarding proteins that are up-regulated in the ORF strain, one was the G protein-signaling modulator 2 [(LogFC = 2.43) (TcrWFU_07618)] belonging to a group of proteins that have the ability to activate G proteins. They are proteins involved in different biological processes, for example, cell division, establishment and organization of mitotic spindle also playing an important role in asymmetric cell divisions via NuMA (Du et al., 2001; Zhu et al., 2011; Kiyomitsu and Cheeseman, 2012). Some reports focus on G protein regulators (GPRs) in asymmetric positioning of the mitotic spindle in the early C. elegans embryo (Manning, 2003). Also up-regulated in ORF were Dystrobrevin-1 (TcrWFU_10303) and aquaporin-4 (TcrWFU_02255), which have also been found in the Echinococcus encystment process, playing a role of vasopressor and water reabsorption (Fan et al., 2020). These mechanisms have been related to the transport of proteins, carbohydrates, and other substances (Park and Kwon, 2015).
Proteins related to maintenance of animal organ identity were also upregulated in ORF (Supplementary File 1); among them, Fibronectins type III domain [(TcrWFU_10230, TcrWFU_10661, TcrWFU_10672 and TcrWFU_08319) (LogFC = 7.03, 8.13, 4.09 and 7.31)] are conserved proteins widely distributed in different species. Three of them have been reported as oncospheral proteins, although our results here shown demonstrated that these genes were also expressed in the larva. Also upregulated in ORF was an oncospheral fibronectin type III domain which has been used to develop a successful vaccine against porcine cysticercosis (Flisser et al., 2004; Gonzalez et al., 2005; Gauci et al., 2012). Finally, ORF cysticerci have a quicker reproductive rate in the peritoneal cavity of mice (Everhart et al., 2004; Fragoso et al., 2008; Willms and Zurabian, 2010), suggesting that upregulation of proteins related to cell proliferation and differentiation are involved. These findings allow focusing future studies to determine the causative agent of the lack of scolex in ORF.
Genomic Variation Analysis of Differentially Expressed Genes
In addition to the genes that were found differentially expressed between the two strains, we identified a group of 12 other genes containing moderate and high impact SNPs, besides their differential expression state (Supplementary File 1). Three genes including high impact mutations that could affect their function, were found over-expressed in ORF. Protocadherin (TcrWFU_08480) which is necessary for cell-adhesion, might have an important role in the scolex formation. Other examples were the suppressor of tumorigenicity 7 protein (ST7) (TcrWFU_09775) and the adenylate cyclase type 9 (TcrWFU_09008), which are signaling proteins that have been well characterized in human development but are poorly studied in cestodes. If associated to scolex development, over expression of these regulators could be affected by mutations in those genes. For example, the gene that codes for Iroquois-class homeodomain protein IRX-5 (TcrWFU_07074), which is a transcription factor that has been associated with craniofacial development through modulating the migration of progenitor cells in the branchial arches of Xenopus laevis embryos; mutation of this gene prevents the correct craniofacial development causing dysmorphism (Bonnard et al., 2012). In ORF, we found that the SNP on IRX-5 contains a missense change (c.290T>C) with respect to WFU, suggesting that this mutation could affect its function and contribute to the lack in the development of the scolex.
Other gene with a high impact change (frameshift c.260delG) was the ubiquitin carboxyl-terminal hydrolase 16 [(UBP16) (TcrWFU_00103)], resulting in an amino acid deletion which might affect its important role in the deubiquitination of histone H2A, necessary for the progress of M phase in the cell cycle (Cai et al., 1999; Joo et al., 2007). On the other hand, UBP16 has been observed in X. laevis embryos, playing a role in the regulation of HOX genes through deubiquitination of histone H2A, which leads to the correct development of the anterior-posterior embryonic pattern (Joo et al., 2007).
Another downregulated gene in the ORF strain was Axin-1/DIX (TcrWFU_07383); which has been described as a protein with multiple functions, including the regulation of Wnt/β-catenin (cWnt) signaling pathway. Axin-1/DIX in ORF strain presented a SNP (c.1131delC) producing a frameshift of high-impact effect. Axins could interact with a β-catenin paralog and limit their accumulation; having highly segregated expression patterns along the anteroposterior axis in E. multilocularis and H. microstoma, this would indicate that β-catenin destruction complexes act during larval metamorphosis (Montagne et al., 2019). Dix domain is a family of proteins that form complexes with members of the disheveled protein family (Dvl) and Axin, acting as positive regulators of the Wnt signaling pathway (Shiomi et al., 2005). These complexes are involved in protein-protein interactions and are necessary for the ability of Dvl-1 and Axin to regulate the stability of β-catenin (Kishida et al., 1999).
Taken together, the above data suggest that the lack of scolex in the ORF strain of T. crassiceps could result from the action of multiple genes involved in the antero-posterior development of the parasite. Although a detailed scheme cannot be advanced at this point, something that can be ruled out is the proposal that the lack of scolex was due to the loss of complete chromosomes (Smith et al., 1972). Instead, our results suggest that it was due to a series of smaller changes, including single point mutations or differential expression of gene networks. Survival of this strain in natural conditions appears to be impractical. The role of 533 specific genes in the formation of the scolex could be approached through gene 534 silencing experiments on WFU cysts (Guerrero-Hernández et al., 2022).
Elucidation of the genomes of two strains of T. crassiceps will improve our knowledge about the genomic similarities and differences that exist between taeniid species, allowing novel approaches to a number of unsolved questions on tapeworm infections, including human neurocysticercosis caused by T. solium. We identified a core genome of conserved proteins for a group of cyclophyllidean cestodes, as well as some proteins that are unique to the Taenia genus and the T. crassiceps species. In this respect, improving the resolution of the T. solium genome is a must for future studies of comparative genomics within Taeniidae. Comparison between genomes might become a pipeline for the identification of potential targets and the design of new drugs, diagnostic tools or vaccines. Finally, availability of these high quality assemblies of T. crassiceps (WFU and ORF) can also contribute to defining the usefulness and limits of this animal model for the improvement of our understanding on human cysticercosis.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the accession number JAKROA000000000, associated with the BioProject PRJNA807072.
Ethics Statement
The animal study was reviewed and approved by Institutional Committee for the Care and Use of Laboratory Animals (CICUAL), Biomedical Research Institute-UNAM with the register number 6329.
Author Contributions
KE performed the hybrid genome assembly, gene prediction, RNAseq differential expression and GO/KEGG enrichment analysis; PT performed initial sequencing, PT, DR-V and RJB maintained and isolated the cysts in mice, RJB, KE, DR-V, AC-G and JCC carried out the gene annotation, differential protein expression analysis, AS-F and KE did the comparative genomics (SNP and synteny analyses). JPL and AS-F conceived the project. RJB, KE and DR-V prepared the first draft and AS-F and JPL developed the final version. JPL provided the budget from grants. All authors contributed to the article and approved the submitted version.
Funding
This report was supported in part by grants A1-5-11306 (CONACYT) and [IN 205820] PAPIIT-UNAM.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to R. Grande and J. Verleyen for sequencing and bioinformatic support and to “Unidad de Secuenciación Masiva y Bioinformática” (UUSMB) of the “Laboratorio Nacional de Apoyo Tecnológico a las Ciencias Genómicas,” (CONACyT) in the Biotechnology Institute (UNAM). Infected mice were maintained in the Unidad de Modelos Biológicos at the Biomedical Research Institute. We thank G. Díaz and D. Garzón for the technical support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.876839/full#supplementary-material
Supplementary Figure 1 | Whole mount localization of proliferative cells in T. crassiceps ORF and WFU strains cysticerci. Cysts were labelled with EdU (red) or DAPI (blue) to reveal proliferative cells (phase S) and cell nuclei, respectively. (A) ORF cysts, (B) WFU cysts, and (C) WFU scolex. Scale bars in (A, B) represent 200 μm; bars in (C) represent 100 μm.
Supplementary Figure 2 | Epifluorescence microscopy for karyotype determination in T. crassiceps WFU and ORF strains. Representative micrographs of the two different chromosome numbers observed in mitotic cytons from both strains of cysticerci: 2n=16, 2n=18.
References
Alexa, A., Rahnenführer, J., Lengauer, T. (2006). Improved Scoring of Functional Groups From Gene Expression Data by Decorrelating GO Graph Structure. Bioinformatics 22 (13), 1600–1607. doi: 10.1093/bioinformatics/btl140
Arocker-Mettinger, E., Huber-Spitzy, V., Auer, H., Grabner, G., Stur, M. (1992). Taenia Crassiceps in the Anterior Chamber of the Human Eye: A Case Report. Klin Monbl Augenheilkd 201 (1), 34–37. doi: 10.1055/s-2008-1045865
Battelli, G. (2009). Echinococcosis: Costs, Losses and Social Consequences of a Neglected Zoonosis. Vet. Res. Commun. Suppl 1, 47–52. doi: 10.1007/s11259-009-9247-y
Bei, Y., Hogan, J., Berkowitz, L. A., Soto, M., Rocheleau, C. E., Pang, K. M. (2002). SRC-1 and Wnt Signaling Act Together to Specify Endoderm and to Control Cleavage Orientation in Early C. Elegans Embryos. Dev. Cell 3 (1), 113–125. doi: 10.1016/s1534-5807(02)00185-5
Berriman, M., Haas, B. J., Loverde, P. T., Wilson, R. A., Dillon, G. P., Cerqueira, G. C., et al. (2009). The Genome of the Blood Fluke Schistosoma mansoni. Nature 460, 352–358. doi: 10.1038/nature08160
Bonnard, C., Strobl, A. C., Shboul, M., Lee, H., Merriman, B., Nelson, S. F., et al. (2012). Mutations in IRX5 Impair Craniofacial Development and Germ Cell Migration via SDF1. Nat. Genet. 44 (6), 709–13. doi: 10.1038/ng.2259
Brůna, T., Hoff, K. J., Lomsadze, A., Stanke, M., Borodovsky, M. (2021). BRAKER2: Automatic Eukaryotic Genome Annotation With GeneMark-EP+ and AUGUSTUS Supported by a Protein Database. NAR Genom Bioinform. 3 (1), lqaa108. doi: 10.1093/nargab/lqaa108
Butler, J., MacCallum, I., Kleber, M., Shlyakhter, I. A., Belmonte, M. K., Lander, E. S., et al. (2008). ALLPATHS: De Novo Assembly of Whole-Genome Shotgun Microreads. Genome Res. 18, 810–820. doi: 10.1101/gr.7337908
Cai, S. Y., Babbitt, R. W., Marchesi, V. T. (1999). A Mutant Deubiquitinating Enzyme (Ubp-M) Associates With Mitotic Chromosomes and Blocks Cell Division. Proc. Natl. Acad. Sci. 96, 2828–2833. doi: 10.1073/pnas.96.6.2828
Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: Architecture and Applications. BMC Bioinf. 10 (1), 1–9. doi: 10.1186/1471-2105-10-421
Carè, A., Valtieri, M., Mattia, G., Meccia, E., Masella, B., Luchetti, L., et al. (1999). Enforced Expression of HOXB7 Promotes Hematopoietic Stem Cell Proliferation and Myeloid-Restricted Progenitor Differentiation. Oncogene 18 (11), 1993–2001. doi: 10.1038/sj.onc.1202498
Carpio, A., Fleury, A., Matthew, L., Romo and Abraham, R. (2018). Neurocysticercosis: The Good, the Bad, and The Missing. Expert Rev. Neurother. 18 (4), 289–301. doi: 10.1080/14737175.2018.1451328
Chuck, R. S., Olk, R. J., Weil, G. J., Akduman, L., Benenson, I. L., Smith, M. E., et al. (1997). Surgical Removal of a Subretinal Proliferating Cysticercus of Taeniaeformis Crassiceps. Arch. Ophthalmol. 115 (4), 562–563. doi: 10.1001/archopht.1997.01100150564028
Cingolani, P., Platts, A., Wang le, L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A Program for Annotating and Predicting the Effects of Single Nucleotide Polymorphisms, SnpEff: SNPs in the Genome of Drosophila Melanogaster Strain W1118; Iso-2; Iso-3. Fly 6 (2), 80–92. doi: 10.4161/fly.19695
Coombe, L., Nikolić, V., Chu, J., Birol, I., Warren, R. L. (2020). Ntjoin: Fast and Lightweight Assembly-Guided Scaffolding Using Minimizer Graphs. Bioinformatics 36 (12), 3885–3887. doi: 10.1093/bioinformatics/btaa253
Corallo, D., Trapani, V., Bonaldo, P. (2015). The Notochord: Structure and Functions. Cell. Mol. Life Sci.: CMLS 72 (16), 2989–3008. doi: 10.1007/s00018-015-1897-z
Deplazes, P., Eichenberger, R. M., Grimm, F. (2019). Wildlife-Transmitted Taenia and Versteria Cysticercosis and Coenurosis in Humans and Other Primates. Int. J. Parasitol. Parasites Wildl. 11 (9), 342–358. doi: 10.1016/j.ijppaw.2019.03.013
Deplazes, P., Rinaldi, L., Alvarez-Rojas, C. A., Torgerson, P. R., Harandi, M. F., Romig, T., et al. (2017). Global Distribution of Alveolar and Cystic Echinococcosis. Adv. Parasitol. 95, 315–493. doi: 10.1016/bs.apar.2016.11.001
Dezaki, E. S., Yaghoobi, M. M., Taheri, E., Almani, P. G., Tohidi, F., Gottstein, B., et al. (2016). Differential Expression of Hox and Notch Genes in Larval and Adult Stages of Echinococcus Granulosus. Korean J. Parasitol. 54 (5), 653–658. doi: 10.3347/kjp.2016.54.5.653
Dorais, F. J., Esch, G. W. (1969). Growth Rate of Two Taenia Crassiceps Strains. Exp. Parasitol. 25, 395–398. doi: 10.1016/0014-4894(69)90086-1
Du, Q., Stukenberg, P. T., Macara, I. G. (2001). A Mammalian Partner of Inscuteable Binds NuMA and Regulates Mitotic Spindle Organization. Nat. Cell Biol. 3 (12), 1069–1075. doi: 10.1038/ncb1201-1069.s
Eddy, S. R. (2009). A New Generation of Homology Search Tools Based on Probabilistic Inference. Genome Inform. 23 (1), 205–211. doi: 10.1142/9781848165632_0019
Esch, G. W., Smyth, J. D. (1976). Studies on the In Vitro Culture of Taenia Crassiceps. Int. J. Parasitol. 6, 143–149. doi: 10.1016/0020-7519(76)90071-0
Everhart, M. E., Kuhn, R. E., Zelmer, D. A. (2004). Infrapopulation Dynamics of a Wild Strain of Taenia Crassiceps (WFU) (Cestoda: Taeniidae) in BALB/cJ Mice. J. Parasitol. 90 (1), 79–84. doi: 10.1645/GE-3255
Fan, J., Wu, H., Li, K., Liu, X., Tan, Q., Cao, W., et al. (2020). Transcriptomic Features of Echinococcus Granulosus Protoscolex During the Encystation Process. Korean J. Parasitol. 58, 287–299. doi: 10.3347/kjp.2020.58.3.287
Fleury, A., Moreno García, J., Valdez Aguerrebere, P., de Sayve Durán, M., Becerril Rodríguez, P., Larralde, C., et al. (2010). Neurocysticercosis, a Persisting Health Problem in Mexico. PloS Negl. Trop. Dis. 4 (8), e805. doi: 10.1371/journal.pntd.0000805
Flisser, A., Gauci, C. G., Zoli, A., Martinez-Ocana, J., Garza-Rodriguez, A., Dominguez-Alpizar, J. L., et al. (2004). Induction of Protection Against Porcine Cysticercosis by Vaccination With Recombinant Oncosphere Antigens. Infect. Immun. 72 (9), 5292–5297. doi: 10.1128/IAI.72.9.5292-5297.2004
Flisser, A., Rodríguez-Canul, R., Willingham, A. L. (2006). Control of the Taeniosis/Cysticercosis Complex: Future Developments. Vet. Parasitol. 139, 283–292. doi: 10.1016/j.vetpar.2006.04.019
Fragoso, G., Meneses, G., Sciutto, E., Fleury, A., Larralde, C. (2008). Preferential Growth of Taenia Crassiceps Cysticerci in Female Mice Holds Across Several Laboratory Mice Strains and Parasite Lines. J. Parasitol. 94 (2), 551–553. doi: 10.1645/GE-1287.1
François, A., Favennec, L., Cambon-Michot, C., Gueit, I., Biga, N., Tron, F., et al. (1998). Taenia Crassiceps Invasive Cysticercosis: A New Human Pathogen in Acquired Immunodeficiency Syndrome? Am. J. Surg. Pathol. 22 (4), 488–492. doi: 10.1097/00000478-199804000-00015
Freeman, R. S. (1962). Studies on the Biology of Taenia Crassiceps (ZEDER 1800) RudolphI 1810 (CESTODA). Can. J. Zool. 40, 969–990. doi: 10.1139/z62-086
Gabriël, S., Johansen, M. V., Pozio, E., Smit, G. S. A., Devleesschauwer, B., Allepuz, A., et al. (2015). Human Migration and Pig/Pork Import in the European Union: What are the Implications for Taenia Solium Infections? Vet. Parasitol. 213, 38–45. doi: 10.1016/j.vetpar.2015.03.006
Garcia, H. H., Modi, M. (2008). Helminthic Parasites and Seizures. Epilepsia 49, 25–32. doi: 10.1111/j.1528-1167.2008.01753.x
Garrison, E., Gabor, M. (2012). Haplotype-Based Variant Detection From Short-Read Sequencing. arXiv, arXiv:1207.3907.doi: 10.48550/arXiv.1207.3907
Gauci, C. G., Jayashi, C. M., Gonzalez, A. E., Lackenby, J., Lightowlers, M. W. (2012). Protection of Pigs Against Taenia Solium Cysticercosis by Immunization With Novel Recombinant Antigens. Vaccine 30 (26), 3824–3828. doi: 10.1016/j.vaccine.2012.04.019
Goesseringer, N., Lindenblatt, N., Mihic-Probst, D., Grimm, F., Giovanoli, P. (2011). Taenia Crassiceps Upper Limb Fasciitis in a Patient With Untreated Acquired Immunodeficiency Syndrome and Chronic Hepatitis C Infection–the Role of Surgical Debridement. J. Plast. Reconstr. Aesthet. Surg. 64 (7), e174-6. doi: 10.1016/j.bjps.2011.02.011
Gonzalez, A. E., Gauci, C. G., Barber, D., Gilman, R. H., Tsang, V. C., Garcia, H. H., et al. (2005). Vaccination of Pigs to Control Human Neurocysticercosis. Am. J. Trop. Med. Hyg. 72 (6), 837–839. doi: 10.4269/ajtmh.2005.72.837
Good, A. H., Miller, K. L. (1976). Depression of the Immune Response to Sheep Erythrocytes in Mice Infected With Taenia Crassiceps Larvae. Infect. Immun. 14, 449–456. doi: 10.1128/iai.14.2.449-456.1976
Gripper, L. B., Welburn, S. C. (2017). Neurocysticercosis Infection and Disease-A Review. Acta Trop. 166, 218–224. doi: 10.1016/j.actatropica.2016.11.015
Grohme, M. A., Schloissnig, S., Rozanski, A., Pippel, M., Young, G. R., Winkler, S., et al. (2018). The Genome of Schmidtea Mediterranea and the Evolution of Core Cellular Mechanisms. Nature 554 (7690), 56–61. doi: 10.1038/nature25473
Guerrero-Hernández, J., Bobes, R. J., García-Varela, M., Castellanos-Gonzalez, A., Laclette, J. P. (2022). Identification and Functional Characterization of the siRNA Pathway in Taenia Crassiceps by Silencing Enolase A. Acta Trop. 225, 106197. doi: 10.1016/j.actatropica.2021.106197
Guo, L., Accorsi, A., He, S., Guerrero-Hernández, C., Sivagnanam, S., McKinney, S., et al. (2018). An Adaptable Chromosome Preparation Methodology for Use in Invertebrate Research Organisms. BMC Biol. 16 (1), 25. doi: 10.1186/s12915-018-0497-4
Heldwein, K., Biedermann, H. G., Hamperl, W. D., Bretzel, G., Löscher, T., Laregina, D., et al. (2006). Subcutaneous Taenia Crassiceps Infection in a Patient With non-Hodgkin's Lymphoma. Am. J. Trop. Med. Hyg. 75 (1), 108–111. doi: 10.4269/ajtmh.2006.75.108
Hirose, T., Koga, M., Ohshima, Y., Okada, M. (2003). Distinct Roles of the Src Family Kinases, SRC-1 and KIN-22, That are Negatively Regulated by CSK-1 in C. Elegans. FEBS Lett. 534 (1-3), 133–138. doi: 10.1016/s0014-5793(02)03819-x
Huang, Y., Chen, W., Wang, X., Liu, H., Chen, Y., Guo, L., et al. (2013). The Carcinogenic Liver Fluke, Clonorchis Sinensis: New Assembly, Reannotation and Analysis of the Genome and Characterization of Tissue Transcriptomes. PloS One 8 (1), e54732. doi: 10.1371/journal.pone.0054732
International Helminth Genomes Consortium (2019). Comparative Genomics of the Major Parasitic Worms. Nat. Genet. 51 (1), 163–174. doi: 10.1038/s41588-018-0262-1
Ito, A. (2015). Recent Advances and Perspectives in Molecular Epidemiology of Taenia Solium Cysticercosis. Infect. Genet. Evol. 40, 357–367. doi: 10.1016/j.meegid.2015.06.022
Ito, A., Nakao, M., Wandra, T. (2003). Human Taeniasis and Cysticercosis in Asia. Lancet 362, 1918–1920. doi: 10.1016/S0140-6736(03)14965-3
Joo, H. Y., Zhai, L., Yang, C., Nie, S., Erdjument-Bromage, H., Tempst, P., et al. (2007). Regulation of Cell Cycle Progression and Gene Expression by H2A Deubiquitination. Nature 449, 1068–1072. doi: 10.1038/nature06256
Kim, D., Paggi, J. M., Park, C., Bennett, C., Salzberg, S. L. (2019). Graph-Based Genome Alignment and Genotyping With HISAT2 and HISAT-Genotype. Nat. Biotechnol. 37 (8), 907–915. doi: 10.1038/s41587-019-0201-4
Kishida, S., Yamamoto, H., Hino, S., Ikeda, S., Kishida, M., Kikuchi, A. (1999). DIX Domains of Dvl and Axin are Necessary for Protein Interactions and Their Ability to Regulate Beta-Catenin Stability. Mol. Cell Biol. 19 (6), 4414–4422. doi: 10.1128/MCB.19.6.4414
Kiyomitsu, T., Cheeseman, I. M. (2012). Chromosome- and Spindle-Pole-Derived Signals Generate an Intrinsic Code for Spindle Position and Orientation. Nat. Cell Biol. 14 (3), 311–317. doi: 10.1038/ncb2440
Koren, S., Walenz, B. P., Berlin, K., Miller, J. R., Bergman, N. H., Phillippy, A. M. (2017). Canu: Scalable and Accurate Long-Read Assembly via Adaptive K-Mer Weighting and Repeat Separation. Genome Res. 27 (5), 722–736. doi: 10.1101/gr.215087.116
Koziol, U., Jarero, F., Olson, P. D., Brehm, K. (2016). Comparative Analysis of Wnt Expression Identifies a Highly Conserved Developmental Transition in Flatworms. BMC Biol. 14, 10. doi: 10.1186/s12915-016-0233-x
Krzywinski, M., Schein, J., Birol, I., Connors, J., Gascoyne, R., Horsman, D., et al. (2009). Circos: An Information Aesthetic for Comparative Genomics. Genome Res. 19 (9), 1639–1645. doi: 10.1101/gr.092759.109
Lamb, T. M., Harland, R. M. (1995). Fibroblast Growth Factor is a Direct Neural Inducer, Which Combined With Noggin Generates Anterior-Posterior Neural Pattern. Development 121 (11), 3627–3636. doi: 10.1242/dev.121.11.3627
Langmead, B., Salzberg, S. L. (2012). Fast Gapped-Read Alignment With Bowtie 2. Nat. Methods 9 (4), 357–359. doi: 10.1038/nmeth.1923
Larralde, C., Sotelo, J., Montoya, R. M., Palencia, G., Padilla, A., Govezensky, T., et al. (1990). Immunodiagnosis of Human Cysticercosis in Cerebrospinal Fluid. Antigens From Murine Taenia Crassiceps Cysticerci Effectively Substitute Those From Porcine Taenia Solium. Arch. Pathol. Lab. Med. 114, 926–928.
Lechner, M., Findeiss, S., Steiner, L., Marz, M., Stadler, P. F., Prohaska, S. J., et al. (2011). Proteinortho: Detection of (Co-) Orthologs in Large-Scale Analysis. BMC Bioinf. 12, 124. doi: 10.1186/1471-2105-12-124
Lescano, A. G., Zunt, J. (2013). Other Cestodes: Sparganosis, Coenurosis and Taenia Crassiceps Cysticercosis. Handb. Clin. Neurol. 114, 335–345. doi: 10.1016/B978-0-444-53490-3.00027-3
Li, H., Durbin, R. (2009). Fast and Accurate Short Read Alignment With Burrows–Wheeler Transform. Bioinformatics 25 (14), 1754–1760. doi: 10.1093/bioinformatics/btp324
Li, W., Liu, B., Yang, Y., Ren, Y., Wang, S., Liu, C., et al. (2018). The Genome of Tapeworm Taenia Multiceps Sheds Light on Understanding Parasitic Mechanisms and Control of Coenurosis Disease. DNA Res. 25, 499–510. doi: 10.1093/dnares/dsy020
Loos-Frank, B. (2000). An Up-Date of Verster’s, (1969) “Taxonomic Revision of the Genus Taenia Linnaeus” (Cestoda) in Table Format. Syst. Parasitol. 45, 155–183. doi: 10.1023/a:1006219625792
Lustigman, S., Prichard, R. K., Gazzinelli, A., Grant, W. N., Boatin, B. A., McCarthy, J. S., et al. (2012). A Research Agenda for Helminth Diseases of Humans: The Problem of Helminthiases. PloS Negl. Trop. Dis. 6 (4), e1582. doi: 10.1371/journal.pntd.0001582
Mahanty, S., Madrid, E. M., Nash, T. E. (2013). Quantitative Screening for Anticestode Drugs Based on Changes in Baseline Enzyme Secretion by Taenia Crassiceps. Antimicrob. Agents Chemother. 57 (2), 990–995. doi: 10.1128/AAC.01022-12
Maillard, H., Marionneau, J., Prophette, B., Boyer, E., Célerier, P. (1998). Taenia Crassiceps Cysticercosis and AIDS. AIDS 12 (12), 1551–1552. doi: 10.1097/00002030-199812000-00019
Maldonado, L. L., Arrabal, J. P., Rosenzvit, M. C., Kamenetzky, L. (2019). Revisiting the Phylogenetic History of Helminths Through Genomics, the Case of the New. Echinococcus Oligarthrus. Genome 10. doi: 10.3389/fgene.2019.00708
Maldonado, L. L., Assis, J., Araújo, F. M. G., Salim, A. C. M., Macchiaroli, N., Cucher, M., et al. (2017). The Echinococcus Canadensis (G7) Genome: A Key Knowledge of Parasitic Platyhelminth Human Diseases. BMC Genomics 18 (1), 204. doi: 10.1186/s12864-017-3574-0
Manning, D. R. (2003). Evidence Mounts for Receptor-Independent Activation of Heterotrimeric G Proteins Normally In Vivo: Positioning of the Mitotic Spindle in C. Elegans. Sci. STKE 2003 (196), pe35. doi: 10.1126/stke.2003
Moguel, B., Moreno-Mendoza, N., Bobes, R. J., Carrero, J. C., Chimal-Monroy, J., Díaz-Hernández, M. E., et al. (2015). Transient Transgenesis of the Tapeworm Taenia Crassiceps. Springerplus 4, 496. doi: 10.1186/s40064-015-1278-y
Montagne, J., Preza, M., Castillo, E., Brehm, K., Koziol, U. (2019). Divergent Axin and GSK-3 Paralogs in the Beta-Catenin Destruction Complexes of Tapeworms. Dev. Genes Evol. 229 (4), 89–102. doi: 10.1007/s00427-019-00632-w
Morales-Montor, J., Escobedo, G., Vargas-Villavicencio, J. A., Larralde, C. (2008). The Neuroimmunoendocrine Network in the Complex Host-Parasite Relationship During Murine Cysticercosis. Curr. Top. Med. Chem. 8 (5), 400–407. doi: 10.2174/156802608783790866
Mougeot, G., Cambon, M., Dimeglio, V., Menerath, J. M. (1996). Cestodose Larvaire Intra-Oculaire Chez Un Jeune Garçon De 14 Ans En Auvergne [Intraocular Cestode Larva in a 14-Year-Old Boy in Auvergne]. Presse Med. 25 (25), 1168.
Mount, P. S. (1968). The Rostellar Morphogenesis of Two Strains of Larval Taenia Crassiceps. Unpublished Master's Thesis (Winston-Salem, North Carolina: Wake Forest University).
Nowak, R. M., Jastrzębski, J. P., Kuśmirek, W., Sałamatin, R., Rydzanicz, M., Sobczyk-Kopcioł, A., et al. (2019). Hybrid De Novo Whole-Genome Assembly and Annotation of the Model Tapeworm Hymenolepis Diminuta. Sci. Data 6 (1), 302. doi: 10.1038/s41597-019-0311-3
Ntoukas, V., Tappe, D., Pfütze, D., Simon, M., Holzmann, T. (2013). Cerebellar Cysticercosis Caused by Larval Taenia Crassiceps Tapeworm in Immunocompetent Woman, Germany. Emerg. Infect. Dis. 19 (12), 2008–2011. doi: 10.3201/eid1912.130284
Olson, P. D., Tracey, A., Baillie, A., James, K., Doyle, S. R., Buddenborg, S. K., et al. (2020). Complete Representation of a Tapeworm Genome Reveals Chromosomes Capped by Centromeres, Necessitating a Dual Role in Segregation and Protection. BMC Biol. 18 (1), 165. doi: 10.1186/s12915-020-00899-w
Olson, P. D., Zarowiecki, M., James, K., Baillie, A., Bartl, G., Burchell, P., et al. (2018). Genome − Wide Transcriptome Profiling and Spatial Expression Analyses Identify Signals and Switches of Development in Tapeworms. Evodevo 9, 21. doi: 10.1186/s13227-018-0110-5
O’Neal, S. E., Flecker, R. H. (2015). Hospitalization Frequency and Charges for Neurocysticercosis, United States. 2003–2012. Emerg. Infect. Dis. 21, 969–976. doi: 10.3201/eid2106.141324
Park, E. J., Kwon, T. H. (2015). A Minireview on Vasopressin-Regulated Aquaporin-2 in Kidney Collecting Duct Cells. E BP. 13 (1), 1–6. doi: 10.5049/EBP.2015.13.1.1
Petersen, T. N., Brunak, S., von Heijne, G., Nielsen, H.. (2011). SignalP 4.0: Discriminating Signal Peptides From Transmembrane Regions. Nat. Methods 8 (10), 785–86. doi: 10.1038/nmeth.1701
Quaggin, S. E., Vanden, H., Igarashi, P. (1998). Pod-1, a Mesoderm-Specific Basic-Helix-Loop-Helix Protein Expressed in Mesenchymal and Glomerular Epithelial Cells in the Developing Kidney. Mech. Dev. 71 (1-2), 37–48. doi: 10.1016/s0925-4773(97)00201-3
Roberts, A., Trapnell, C., Donaghey, J., Rinn, J. L., Pachter, L. (2011). Improving RNA-Seq Expression Estimates by Correcting for Fragment Bias. Genome Biol. 12 (3), R22. doi: 10.1186/gb-2011-12-3-r22
Robinson, M. D., McCarthy, D. J., Smyth, G. K. (2010). Edger: A Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26 (1), 139–140. doi: 10.1093/bioinformatics/btp616
Sally, C. Y., Chau, J., Freeman, R. S. (1976). Intraperitoneal Passage of Taenia Crassiceps in Rats. J. Parasitol. 62, 837–839. doi: 10.2307/3278978
Sciutto, E., Fragoso, G., Diaz, M. L., Valdez, F., Montoya, R. M., Govezensky, T., et al. (1991). Murine Taenia Crassiceps Cysticercosis: H-2 Complex and Sex Influence on Susceptibility. Parasitol. Res. 77 (3), 243–246. doi: 10.1007/BF00930866
Sciutto, E., Fragoso, G., Larralde, C. (2011). Taenia Crassiceps as a Model for Taenia Solium and the S3Pvac Vaccine. Parasite Immunol. 33 (1), 79–80. doi: 10.1111/j.1365-3024.2010.01257.x
Shea, M., Maberley, A. L., Walters, J., Freeman, R. S., Fallis, A. M. (1973). Intraocular Taenia Crassiceps (Cestoda). Trans. Am. Acad. Ophthalmol. Otolaryngol. 77 (6), OP778–OP783.
Shiomi, K., Kanemoto, M., Keino-Masu, K., Yoshida, S., Soma, K., Masu, M. (2005). Identification and Differential Expression of Multiple Isoforms of Mouse Coiled-Coil-DIX1 (Ccd1), a Positive Regulator of Wnt Signaling. Brain Res. Mol. Brain Res. 135, 169–180. doi: 10.1016/j.molbrainres.2004.12.002
Silva, C. V., Costa-Cruz, J. M. (2012). A Glance at Taenia Saginata Infection, Diagnosis, Vaccine, Biological Control and Treatment. Infect. Disord. - Drug Targets 10, 313–321. doi: 10.2174/187152610793180894
Smith, J. K., Esch, G. W., Kuhn, R. E. (1972). Growth and Development of Larval Taenia Crassiceps (Cestoda). I. Aneuploidy in the Anomalous ORF Strain. Int. J. Parasitol. 2 (2), 261–263. doi: 10.1016/0020-7519(72)90014-8
Smit, A. F. A., Hubley, R. (2016). RepeatMasker. Available at http://www.repeatmasker.org.
Smit, A. F. A., Hubley, R., Green, P. (2018). RepeatMasker. Available at http://www.repeatmasker.org.
Špakulová, M., Orosová, M., Mackiewicz, J. S. (2011). Cytogenetics and Chromosomes of Tapeworms (Platyhelminthes, Cestoda). Adv. Parasitol. 74, 177–230. doi: 10.1016/B978-0-12-385897-9.00003-3
Stanke, M., Keller, O., Gunduz, I., Hayes, A., Waack, S., Morgenstern, B. (2006). AUGUSTUS: Ab Initio Prediction of Alternative Transcripts. Nucleic Acids Res. 34 (suppl_2), W435–W439. doi: 10.1093/nar/gkl200
Teven, C. M., Farina, E. M., Rivas, J., Reid, R. R. (2014). Fibroblast Growth Factor (FGF) Signaling in Development and Skeletal Diseases. Genes Dis. 1 (2), 199–213. doi: 10.1016/j.gendis.2014.09.005
Toledo, A., Fragoso, G., Rosas, G., Hernández, M., Gevorkian, G., López-Casillas, F., et al. (2001). Two Epitopes Shared by Taenia Crassiceps and Taenia Solium Confer Protection Against Murine T. Crassiceps Cysticercosis Along With a Prominent T1 Response. Infect. Immun. 69 (3), 1766–1773. doi: 10.1128/IAI.69.3.1766-1773.2001
Torgerson, P. R., Devleesschauwer, B., Praet, N., Speybroeck, N., Willingham, A. L., Kasuga, F., et al. (2015). World Health Organization Estimates of the Global and Regional Disease Burden of 11 Foodborne Parasitic Diseases 2010: A Data Synthesis. PloS Med. 12 (12), e1001920. doi: 10.1371/journal.pmed.1001920
Tsai, I. J., Zarowiecki, M., Holroyd, N., Garciarrubio, A., Sanchez-Flores, A., Brooks, K. L., et al. (2013). The Genomes of Four Tapeworm Species Reveal Adaptations to Parasitism. Nature 496, 57–63. doi: 10.1038/nature12031
Vargas-Villavicencio, J. A., Larralde, C., De León-Nava, M. A., Morales-Montor, J. (2005). Regulation of the Immune Response to Cestode Infection by Progesterone is Due to its Metabolism to Estradiol. Microbes Infect. 3), 485–493. doi: 10.1016/j.micinf.2004.12.015
Vaser, R., Sovic, I., Nagarajan, N., Sikic, M. (2002). Racon-Rapid Consensus Module for Raw De Novo Genome Assembly of Long Uncorrected Reads. Bioinformatics 452-464, 6. doi: 10.1101/gr.214270.116
Walker, B. J., Abeel, T., Shea, T., Priest, M., Abouelliel, A., Sakthikumar, S., et al. (2014). Pilon: An Integrated Tool for Comprehensive Microbial Variant Detection and Genome Assembly Improvement. PloS One 9 (11), e112963. doi: 10.1371/journal.pone.0112963
Wang, X., Chen, W., Huang, Y., Sun, J., Men, J., Liu, H., et al. (2011). The Draft Genome of the Carcinogenic Human Liver Fluke Clonorchis sinensis. Genome Biol. 12 (10), R107. doi: 10.1186/gb-2011-12-10-r107
Wang, S., Wang, S., Luo, Y., Xiao, L., Luo, X., Gao, S., et al. (2016). Comparative Genomics Reveals Adaptive Evolution of Asian Tapeworm in Switching to a New Intermediate Host. Nat. Commun. 7, 12845. doi: 10.1038/ncomms12845
Wellik, D. M. (2009). Hox Genes and Vertebrate Axial Pattern. Curr. Top. Dev. Biol. 88, 257–278. doi: 10.1016/S0070-2153(09)88009-5
Willms, K., Zurabian, R. (2010). Taenia Crassiceps: In Vivo and In Vitro Models. Parasitology 137, 335–346. doi: 10.1017/S0031182009991442
World Health Organization (2011). Working to Overcome the Global Impact of Neglected Tropical Diseases. First WHO Report on Neglected Tropical Diseases Vol. 86 (Geneva, Switzerland:Wkly Epidemiol Rec), 113–120.
Yang, D., Fu, Y., Wu, X., Xie, Y., Nie, H., Chen, L., et al. (2012). Annotation of the Transcriptome From Taenia Pisiformis and its Comparative Analysis With Three Taeniidae Species. PloS One 7, e32283. doi: 10.1371/journal.pone.0032283
Zammarchi, L., Strohmeyer, M., Bartalesi, F., Bruno, E., Muñoz, J., Buonfrate, D. (2013). Epidemiology and Management of Cysticercosis and Taenia Solium Taeniasis in Europe, Systematic Review 1990–2011. PloS One 8, e69537. doi: 10.1371/journal.pone.0069537
Zhou, Y., Zheng, H., Chen, Y., Zhang, L., Wang, K., Guo, J., et al. (2009). The Schistosoma Japonicum Genome Reveals Features of Host-Parasite Interplay. Nature 460, 345–351. doi: 10.1038/nature08140
Keywords: Taenia crassiceps, ORF, WFU, comparative genomics, RNA-seq, differential expression, cysticercosis, animal model
Citation: Bobes RJ, Estrada K, Rios-Valencia DG, Calderón-Gallegos A, de la Torre P, Carrero JC, Sanchez-Flores A and Laclette JP (2022) The Genomes of Two Strains of Taenia crassiceps the Animal Model for the Study of Human Cysticercosis. Front. Cell. Infect. Microbiol. 12:876839. doi: 10.3389/fcimb.2022.876839
Received: 15 February 2022; Accepted: 12 April 2022;
Published: 10 May 2022.
Edited by:
Laurence A. Marchat, Instituto Politécnico Nacional, MexicoReviewed by:
Klaus Rüdiger Brehm, Julius Maximilian University of Würzburg, GermanyPablo Smircich, Instituto de Investigaciones Biológicas Clemente Estable (IIBCE), Uruguay
Copyright © 2022 Bobes, Estrada, Rios-Valencia, Calderón-Gallegos, de la Torre, Carrero, Sanchez-Flores and Laclette. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Juan P. Laclette, laclette@iibiomedicas.unam.mx; Alejandro Sanchez-Flores, alejandro.sanchez@ibt.unam.mx
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work