- Open Access
Structural and functional characterization of a novel gene, Hc-daf-22, from the strongylid nematode Haemonchus contortus
Parasites & Vectorsvolume 9, Article number: 422 (2016)
The strongylid nematode Haemonchus contortus is a parasite of major concern for modern livestock husbandry because hostile environmental conditions may induce diapause in the early fourth-stage larvae.
A new gene Hc-daf-22 was identified which is the homologue of Ce-daf-22 and human SCPx. Genome walking and RACE were performed to obtain the whole cDNA and genomic sequence of this gene. Using qRT-PCR with all developmental stages as templates to explore the transcription level and micro-injection was applied to confirm the promoter activity of the 5′-flanking region. Overexpression, rescue and RNA interference experiments were performed in N2, daf-22 mutant (ok 693) strains of C. elegans to study the gene function of Hc-daf-22.
The full length gene of Hc-daf-22 (6,939 bp) contained 16 exons separated by 15 introns, and encoded a cDNA of 1,602 bp (533 amino acids, estimated at about 59.3 kDa) with a peak in L3 and L4 in transcriptional level. The Hc-DAF-22 protein was consisted of a 3-oxoacyl-CoA thiolase domain and a SCP2 domain and evolutionarily conserved. The 1,548 bp fragment upstream of the 5′-flanking region was confirmed to have promoter activity compared with 5′-flanking region of Ce-daf-22. The rescue experiment by micro-injection of daf-22 (ok693) mutant strain showed significant increase in body size and brood size in the rescued worms with significantly reduced or completely absent fat granules confirmed by Oil red O staining, indicating that Hc-daf-22 could partially rescue the function of Ce-daf-22. Furthermore, RNAi with Hc-daf-22 could partially silence the endogenous Ce-daf-22 in N2 worms and mimic the phenotype of daf-22 (ok693) mutants.
The gene Hc-daf-22 was isolated and its function identified using C. elegans as a model organism. Our results indicate that Hc-daf-22 shared similar characteristics and function with Ce-daf-22 and may play an important role in peroxisomal β-oxidation and the development in H. contortus.
The free-living nematode Caenorhabditis elegans normally develops from egg to adult through L1 to L4 larval stages in suitable environment. However, when worm density is too high or in starvation, L3 larva can transform into a dauer stage, in which worms would accumulate large amount of fat and have more spindly pharynx with thicker cuticle . These physiological and metabolic changes can ensure a longer survival time for worms in harsh environments; normal worms can live for about two weeks while dauer worms can live for months . In addition, these changes can be reversed when conditions become suitable for nematode growth, and worms can develop into regular reproductive adult stage . Caenorhabditis elegans can sense population density through a group of small molecules named dauer pheromone which is secreted into environment . The dauer pheromone consists of ascaroside derivatives with short, fatty-acid like side chains .
Ce-daf-22 is homologous to the N-terminal portion of the invertebrate sterol carrier protein x (SCPx) , and UNC-24 is homologous to the C-terminal portion of the vertebrate SCPx  which is one of these lip/sterol intra cellular transport proteins  belonging to the well-characterized SCP2 gene family  and containing a sterol-binding domain (SCP2 domain). The vertebrate SCPx genes encode bipartite proteins containing 3-ketoacyl-CoA thiolase and sterol carrier protein-2 (SCP2) domains, located at the N- and C-terminus, respectively [8–10]. SCPx is found exclusively in peroxisomes . It has been demonstrated that SCPx plays a unique role in the peroxisomal β-oxidation of branched-chain fatty acids  and in bile acid formation [4, 12]. In C. elegans, the daf-22 mutant was originally identified as a mutant failed to induce dauer formation and was defective in peroxisomal β-oxidation and accumulated massive amounts of very-long-chain fatty acid (VLCFAs) and very long acyl-CoAs which are toxic . The massive accumulation of fatty acids and their acyl-CoAs causes severe developmental defects. Researchers found that Ce-daf-22 was also required for the biosynthesis of the short-chain fatty acid-derived side chains of the dauer pheromone , which indicated that there was certain link between dauer pheromone production and the general fatty acid metabolism . Putative SCPx genes have been identified in parasitic nematodes Ascaris suum (ADY43045.1) and Brugia malayi  based on sequence similarities. However, the biological importance of SCPx in parasitic nematodes has not been yet reported.
Haemonchus contortus could also reach a special developmental stage in harsh environments such as low oxygen tension, short photoperiod, low temperature and host immune response, called diapause to avoid being expelled by hosts . Little is known about details of the fatty acid β-oxidation pathway in this parasite and whether this pathway is also involved in the formation of diapause. Although significant progress in the molecular characterization of animal parasitic nematodes has been made during the last decade, the complexity of their life-cycles, the difficulties in collection of larvae of every stage, and the fact that adults cannot be cultivated in vitro, still make such studies difficult. Caenorhabditis elegans has been used as a heterologous expression system for defining the gene function of parasitic nematodes for many years. In this study, the complete cDNA of Ce-daf-22 orthologue was first identified in H. contortus, with the name of Hc-daf-22. The 5′-flanking region of Hc-daf-22 was examined for its promoter activity in C. elegans and the coding sequence containing the thiolase domain of Hc-daf-22 was expressed in the C. elegans daf-22 (ok693) strain by micro-injection in order to see whether it can rescue peroxisomal β-oxidation and accumulation of fat droplets in the intestine of daf-22 loss-of-function mutant or not.
Adults of H. contortus (ZJ strain) were collected from sheep abomasa (sheep abomasa were obtained from the Hu Zhou Slaughter house), washed by phosphate buffer saline several times and stored in liquid nitrogen for later use. The C. elegans strains were maintained on NGM agar plates seeded with E. coli (OP50 strain) at 20 °C unless otherwise stated. The strains employed here were Bristol N2 and RB859 daf-22 Y57A10C.6 (ok693) II. The daf-22 (ok693) mutant used in this study was obtained from Caenorhabditis Genetics Center (Minneapolis, MN, USA). This mutant strain has a 790 bp deletion of Ce-daf-22 gene and shows different ascaroside biosynthesis which leads to a phenotype of dauer pheromone production defect and fat associated bodies increase.
Isolation of the full-length genomic cDNA of Hc-daf-22
The translated amino acids of Ce-daf-22 were used to search the European Bioinformatics Institute’s parasite genomes database using tBLASTN algorithm. An EST of H. contortus (BF662749) with significant similarity to Ce-daf-22 was identified. Full-length cDNA was obtained from total RNA drawn from H. contortus adults using 5′ and 3′ RACE kit (Takara, Dalian, China) according to the manufacturer’s instructions with gene-specific primers (Table 1). The products were cloned into pMD18-T vector and sequenced. To verify the full-length sequence of Hc-daf-22, primers were designed based on the results of 5′ and 3′ RACE to perform RT-PCR. The product was linked into pMD18-T vector and sequenced to confirm the whole sequence of this gene.
Acquisition of full-length genomic DNA of Hc-daf-22
Adult H. contortus worms, thawed from liquid nitrogen, were washed three times in phosphate buffer saline and homogenized in an ice bath. Total genomic DNA was extracted according to the protocol of Genomic DNA extraction kit (Takara, Dalian, China) with gene-specific primers (Table 1) designed based on the cDNA sequences and used to amplify the entire Hc-daf-22 by long-PCR from total genomic DNA of adult H. contortus using LA Taq. The cycling conditions were: 94 °C for 1 min; 94 °C for 45 s; 60 °C for 30 s; 72 °C for 3 min for 30 cycles, with a final extension at 72 °C for 10 min. The PCR products were cloned into pMD18-T vector and assembled manually after carefully sequenced. The complete cDNA sequence and genomic DNA sequence were analyzed to determine intron/exon boundaries, following the AG-GT rule.
Amplification of the 5′-flanking region of Hc-daf-22 gene
The 5′-flanking region of Hc-daf-22 was amplified by using a genome walking kit (Takara) with primers designed according to the genomic sequence acquired as described above. The primary PCR was carried out using an adapter primer AP1 and specific primer Hc-daf-22-SP1 with ~20 ng of total genomic DNA from H. contortus adults. The cycling conditions were: 1 min at 94 °C, 1 min at 98 °C, 30 s at 94 °C, 1 min at 65 °C, 2 min at 72 °C for 5 cycles; then 30 s at 94 °C, 3 min at 25 °C, 2 min at 72 °C, 30 s at 94 °C, 1 min at 65 °C, 2 min at 72 °C, 30 s at 94 °C, 1 min at 44 °C, 2 min at 72 °C for 15 cycles; followed by 10 min at 72 °C for the final extension. Primers AP1 and Hc-daf-22-SP2 were used at the second round of PCR using 1 μl first round PCR product as a template undergoing the following conditions: 15 cycles of 30 s at 94 °C, 1 min at 65 °C, 2 min at 72 °C, 30 s at 94 °C, 1 min at 44 °C, and 2 min at 72 °C; followed by 10 min at 72 °C for the final extension. Then 1 μl of the second round PCR product was used as template for the third round of PCR with primers AP1 and Hc-daf-22-SP3. The reaction conditions were the same as for the second round. PCR products of the third round were purified and cloned into pMD18-T vector for further sequencing.
Sequence analysis and phylogenetic analysis
Homology of Hc-daf-22 gene analysis was carried out using blastp at the National Center for Biotechnology Information . Sequence similarities between the amino acid sequence of Hc-daf-22 and other related species (C. elegans, Brugia malayi and Ascaris suum) were aligned using Clustal W method. Protein motifs were identified by searching the databases of PROSITE  and Pfam . Phylogenetic analyses were performed using neighbor-joining (NJ) and maximum likelihood (ML) methods in Molecular Evolutionary Genetic Analysis (MEGA) .
Determination of Hc-daf-22 mRNA transcription levels in different developmental stages
Quantitative reverse transcription PCR (qRT-PCR) with specific primers was performed using different developmental stages of H. contortus as templates to determine the mRNA transcriptional levels (Table 1). Briefly, the total RNA of each sample, extracted by Trizol (Invitrogen, Shanghai, China), was used to synthesize the first-strand cDNA by random priming using ReverTra Ace-α (Cat. No. FSK-100, TOYOBO, Shanghai, China). The qRT-PCR (25 μl) using the THUNDERBIRD SYBR qPCR Mix (TOYOBO, Shanghai, China) in an ABI 7300 thermal cycler was performed under the following conditions: 50 °C for 2 min and 95 °C for 1 min for the first cycle followed by 95 °C for 15 s, 60 °C for 15 s and 72 °C for 31 s for 40 cycles. The dissociation curve was generated under the following conditions: 95 °C for 15 s, 60 °C for 1 min, 95 °C for 15 s and 60 °C for 15 s. Each sample was applied in triplicate using β-tubulin 8–9 (Accession: M76493) gene as a normaliser using specific primers (Table 1). The whole experiment was repeated three times. The mean threshold cycle (CT) values were used for further analysis.
Hc-daf-22 promoter amplification and transformation of C. elegans
The affliction of Hc-daf-22 gene promoter was carried out using specific primers (Table 1) which were designed based on the 5′-flanking region of Hc-daf-22. The PCR product was purified and cloned into the upstream of gfp region of the pPD 95.77 expression vector to verify the ability of driving the green fluorescent protein (GFP) expression in C. elegans. This recombinant plasmid, designated as pPD 95.77Hc-daf-22-prom, was microinjected into the gonad of young adults of the wild type (N2) hermaphrodites as described together with another plasmid pRF4, which bears a dominant mutant allele of the rol-6 gene, each at a final concentration of 50 μg/ml. The progeny worms showing the roller phenotype were collected to determine the expression patterns of GFP using a compound fluorescence microscope (Olympus IX71).
Expression of Hc-daf-22 in N2 and daf-22 (ok693) strains of C. elegans
The coding sequence containing the thiolase domain of Hc-daf-22 was amplified by PCR with primers listed in Table 1 with SmaI/KpnI sites underlined. The product was purified and cloned into pPD 95.77 vector between the 2061 bp Ce-daf-22 promoter region and the gfp region. The promoter region was amplified by PCR from N2 strain genomic DNA using the primers listed in Table 1 with restriction sites XbaI and SmaI underlined. The recombinant plasmid was microinjected into N2 and daf-22 (ok693) strains as described above. Transgenic worms displaying the roller phenotype were further analyzed for examination of GFP activity as described above. CeP-pPD95.77-Hc-daf-22d (pPD95.77 vector with Ce-daf-22 gene promoter region and Hc-daf-22 gene functional domain region) was also used as a rescue plasmid, which was microinjected into daf-22 (ok693) mutants as described above. Transformants were selected based on GFP expression and rescued homozygous mutants were identified by fat staining, body size , growth rate  and brood size  measurements.
Oil red O staining was performed as described . Briefly, 200–300 young adult worms were washed three times with 1× PBS and settled by gravity. Worms were resuspended with 120 μl of PBS, together with an equal volume of 2× MRWB buffer containing 2 % paraformaldehyde (PFA) to permeabilize the cuticle. The 2× MRWB buffer contained the following: 160 mM KCl, 40 mM NaCl, 14 mM Na2EGTA, 1 mM spermidine-HCl, 0.4 mM spermine, 30 mM Na-PIPES pH 7.4, 0.2 % β-mercaptoethanol. Worms were gently rocked for 1 h at room temperature, and allowed to settle by gravity. Worms were then washed with 1× PBS to remove PFA one time, and resuspended in 60 % isopropanol for 15 min at room temperature to dehydrate. Isopropanol was removed after the worms had settled and worms were incubated overnight with rocking after adding 1 ml of Oil red O staining solution. On the following day, the staining solution was removed and 200 μl of 1 × PBS 0.01 % Triton X-100 added. Worms were photographed with an Olympus color camera equipped with DIC optics.
Hc-daf-22d RNA interference
Hc-daf-22d was 857 bp, from 237–1,093 bp which contained a functional thiolase domain of Hc-daf-22 gene. This fragment was amplified with specific primers (Table 1) containing the restriction sites Hind III and Nhe. The product was purified and cloned into the L4440 vector. The recombinant plasmid was transformed into HT115 (DE3) cells, an RNase III-deficient Escherichia coli strain with isopropyl-β-D-thiogalac-topyranoside-inducible T7 polymerase activity. The empty L4440 vector was also transformed into HT115 cells as control. Experiments were carried out using the standard bacterial feeding protocol  and repeated at least three times. Worms were observed under a microscope and any change in fat storage was confirmed by Oil red O staining.
Statistical analysis for Hc-daf-22 mRNA transcription levels and parameters of C. elegans were carried out using one-way ANOVA in SPSS 16.0. P-values < 0.05 were considered statistically significant. Graphs were made by GraghPad Prism 6.
Characterization of Hc-daf-22 gene from H. contortus
A 1,777 bp transcript was obtained by overlapping 3′ and 5′ RACE results, including a 1,602 bp open reading frame (GenBank: HQ738470.1), a 58 bp 5′UTR and a 117 bp 3′ UTR. The whole gene (from the ATG initiation codon to the TAA stop codon) was 6,939 bp and consisted of 16 exons separated by 15 introns (Fig. 1). The 1,548 bp 5′-flanking region of Hc-daf-22 was obtained by genome walking. Further analysis indicated that the sequence contained a TATA box, GATA-1, CCAAT binding factor, Unc-86, AP-1, AP-2, and NFγB sites (Fig. 2). The Hc-daf-22 ORF encoded a 533-amino-acid protein (GenBank: AEO14647.1) with a mass of 59.2 kDa as predicted in the PROSITE database, which included the thiolases acyl-enzyme intermediate signature and thiolases signature 2 (Fig. 3a). The sequence was aligned with the DAF-22 protein of C. elegans, the sterol carrier protein of Brugia malayi, the non-specific lipid-transfer protein of Ascaris suum, sterol carrier proteins of chicken, Drosophila melanogaster, Homo sapiens and Mus musculus which showed similarities as follows: 80.3 %, 61.9 %, 73.0 %, 47.5 %, 45.6 %, 47.8 % and 48.2 %, respectively (Fig. 3a). Hc-DAF-22 showed a significantly higher similarity to parasitic nematodes and the Ce-DAF-22 in amino acids (Fig. 3b).
Transcriptional level of Hc-daf-22 throughout the lifespan of H. contortus
The transcription level of Hc-daf-22 was determined at different stages throughout the lifespan of H. contortus by qRT-PCR. Hc-daf-22 was transcribed at detectable levels at L3, L4, adult male and female stages, with a peak in L3 and L4 stages. There was no significant difference in expression levels between male and female adults (t-test: t = 0.2052, df = 4, P = 0.8474) (Fig. 4).
Promoter activity analysis of the 5′-flanking region of Hc-daf-22 in C. elegans
The N2 strain of C. elegans was transformed by microinjection with the reconstructed plasmid as described above. The Ce-daf-22 promoter was microinjected as control. Transgenic lines exhibiting the roller phenotype were selected. Fluorescence microscopy showed that GFP was mainly localized in the distal, middle and anterior part of the intestine, pharynx and excretory cells (Fig. 5) in the transgenic Hc-daf-22-promoter::gfp worms. However, compared to transgenetic Ce-daf-22-promoter::gfp worms, which could drive GFP expression in the hypodermis and throughout the intestine (Fig. 5), the level of GFP expression was clearly lower as we can see in the picture.
Expression of Hc-daf-22 in N2 strains
In order to investigate the gene function of Hc-daf-22 in vivo, the Hc-daf-22 coding region was overexpressed in N2 transgenetic line, using Ce-daf-22 5′-flanking region as promoter. As the results showed, the expression of Hc-daf-22::gfp fusion protein could be observed in all life-cycle stages (L1, L2, L3, L4 and adult). GFP was mainly localized in the whole intestine, and more strongly, in the intestinal nuclei (Fig. 6).
Rescue experiment of daf-22 (ok693) strain with Hc-daf-22
Based on existing data, Ce-daf-22 mutant tends to accumulate enlarged spherical intracellular structures which were confirmed to be fat droplets as a distinct phenotype and to cause defect in post-embryonic development . To testify whether Hc-daf-22 could rescue the loss function of Ce-daf-22 in daf-22 (ok693) strain, daf-22 (ok693) strain was transformed with Hc-daf-22::gfp. GFP could be detected in all life stages mainly in the intestine and tail parts of the worms (Fig. 7). The rescued transformed lines of the mutant strain exhibited reduced fat storage, especially significant reduction or total disappearance of huge fat granules (Fig. 8). In addition, there was a significant increase in the body length and width of young adults of the transformed lines compared to the daf-22 (ok693) mutants (t-test: t = 10.79, df =17, P < 0.0001; t = 8.857, df = 17. P < 0.0001). However, these measurements remained significantly lower than those of the N2 strain (t-test: t = 5.023, df = 17. P < 0.0001; t = 5.550, df = 17. P < 0.0001) (Fig. 9c, d); the same pattern was observed for the brood size (t-test: t = 8.605, df = 8. P < 0.0001; t = 3.333, df = 7. P = 0.0125) (Fig. 9a). No significant difference was observed in the growth rate between daf-22 (ok693) and rescued worms (t-test: t = 0.378, df = 4, P = 0.7247) and both were significantly lower in contrast with that of N2 strain (t-test: t = 166.3, df = 4. P < 0.0001; t = 144.5, df = 4; P < 0.0001). Furthermore, no worms could grow into adult stage either in daf-22 (ok693) mutant or rescued worms (Fig. 9b).
RNAi in C. elegans N2 strain
RNAi experiment was carried out as described above to verify whether the Hc-daf-22 gene could silence the endogenous mRNA of Ce-daf-22 in C. elegans N2 strain. The relative quantification of mRNA levels of N2 and Hc-daf-22 RNAi strain showed a significant reduction in gene transcription after RNAi with Hc-daf-22 gene (t-test: t = 3.563, df = 4, P = 0.0235) (Fig. 10). Worms grown on Hc-daf-22 (RNAi) plates (from egg to adult) at 20 °C were observed. There were no significant changes in body length, body width and growth rate compared with N2 worms (data not shown). However, the brood size of Hc-daf-22 RNAi line was significantly decreased (t-test: t = 4.308, df = 4, P = 0.0126) and worms also displayed a phenotype that was richer in fat storage and accumulating more fat granules than the control worms (Fig. 11), as seen in the daf-22 (ok693) mutants.
Recently, the whole genome of H. contortus was sequenced ; however, the precision and function annotation of the genes still need to be improved. In the present study, a new gene of H. contortus, Hc-daf-22, was isolated and characterized both in relation to gene structure and function. Hc-daf-22 is composed by 16 exons separated by 15 introns, and this composition is identical to human  and mouse  SCPx genes. Comparison of amino acid sequences of Hc-DAF-22 among different species revealed that the thiolase domain and SCP2 domain were conserved in nematodes and mammals. The C-terminal of Hc-DAF-22 contained a tripeptide (Ala-Lys-Ile), which has already been proven to be functional in porcine species , in contrast to the peroxisomal-targeting signal 1 (PTS1; typically Ser-Lys-Leu [27, 28] that are conserved among yeasts, plants, insects and mammals. Moreover, the C-terminal of Hc-DAF-22 was also found containing the amino acids that are important for SCP2 sterol transfer function (Fig. 3). The SCPx domain in the N-terminal of Hc-DAF-22 was more close to C. elegans compared to other species; however, thiolase and SCP-2 domains in C. elegans are coded by two separate genes and not fused together as one protein . These indicate that the SCPx domain is conserved through most species, while the fusion of SCP2 domain in the Hc-daf-22 gene indicates a more advanced evolutionary position of the parasitic H. contortus in relation to the free-living C. elegans. Due to the similarities in structure and function with Ce-daf-22 and other SCPx genes, Hc-daf-22 may act as a peroxisomal protein and play an important role in worm development.
We obtained the 5′-flanking region of Hc-daf-22 by genome walking to analyze sequence character and verify its promoter activity with the 5′-flanking region of Ce-daf-22 as a control. Results from our analysis showed that the 5′-flanking region of Hc-daf-22 contained a TATA box, GATA-1, CCAAT binding factor, CAC-binding protein, Unc-86, AP-1, AP-2, and NFγB sites, some of which can also be found in the promoter region of the SCPx gene in humans . The micro-injection results indicated that this region shows promoter activity which could drive the GFP expression in the transformed lines in N2 strain mostly in the intestine area. However, compared to the control lines, the level of GFP expression was slightly low. Since the promoter region was forced to function in an alien environment of C. elegans, the driving ability could be affected and causing the relatively low GFP levels. In addition, the difference of gene structure may also lead to the difference between two promoters. Hc-daf-22 contained two functional domains that are SCPx and SCP2 domains, while the Ce-daf-22 only has one SCPx domain; therefore, it was reasonable to infer that Hc-daf-22 function is more complicated than Ce-daf-22 which may lead to a different driving pattern of the promoter. There are other possible explanations such as the first exon of Hc-daf-22 has a distinct signal which may lead to different driving pattern of the promoter just like the first exon of human SCPx gene which has influence on the promoter region .
To demonstrate the gene function of Hc-daf-22 in vivo, transformations of N2 and daf-22 (ok693) strains with cp-Hc-daf-22::gfp were performed by micro-injection. The detection of GFP revealed that the recombinant Hc-DAF-22 could be expressed in the intestine region in a similar pattern to Ce-daf-22 in C. elegans . Besides, Hc-DAF-22 could partially restore the phenotype of the daf-22 mutant (ok693) as the fat granules were significantly reduced or even completely absent (Fig. 8). The subsequent analysis of the post-embryonic development showed that Hc-DAF-22 also partially rescued the body size and brood size of the daf-22 (ok693) mutants (Fig. 9a, c, d), although no significant change was seen in growth rates. Similar results were observed by other researchers when expressing Hc-hsp-90 in C. elegans daf-21 mutant. Hc-HSP-90 could only partially rescue the phenotype of the mutant . Nevertheless, other studies showed a full restore of transgenetic expression of parasitic genes in mutant strains of C. elegans. The expression of A. caninum slo-1 and C. oncophora slo-1 genes in C. elegans slo-1(js379) mutant was able to rescue the phenotype of worm locomotion and phenotypic behavior . In another study, the glutamate-gated chloride channel (GluCl) subunit of H. contortus was expressed in a C. elegans mutant (DA1316), which was able to rescue the ivermectin sensitivity of mutant C. elegans . These data indicate that the expression of parasitic genes in mutants of C. elegans could provide a first step in understanding gene functions in vivo, some genes could fully rescue the functions of C. elegans mutant while some could not. The similarity of parasitic genes with the homologues in C. elegans may play a key role in the successful rate of rescue. In our case, the partial rescue of daf-22 (ok693) mutant with Hc-daf-22 may be caused by the difference of gene structure between Hc-daf-22 and Ce-daf-22 as described above.
RNA interference (RNAi), inducing gene silencing, has been applied successfully to C. elegans and provides a functional genomic platform in a range of organisms including parasitic nematodes. Although RNAi experiments in parasitic nematodes such as in Ostertagia ostertagi , Heligmosomoides polygyrus , H. contortus  have been successful, there is still a limited number of genes which can be successfully silenced with RNAi due to many complicated reasons, including the lack of appropriate methods of dsRNA delivery and in vitro culture systems for parasitic nematodes, differences in RNAi effector protein functionality and in the complement of RNAi effectors between nematodes [33, 35, 37, 38]. In this study, we used RNAi to silence the Ce-daf-22 of C. elegans with Hc-daf-22 by feeding. Based on our observation, the daf-22 (RNAi) phenotype showed more fat storage and larger fat granules in a manner similar to daf-22 (ok693) mutant strain. Furthermore, the relative quantification of mRNA levels showed that Hc-daf-22 could successfully partially silence the Ce-daf-22, which further confirmed the similarity of gene structure and possible functions between these two genes.
Based on the expression pattern in different stages of H. contortus, it was noteworthy that Hc-daf-22 reached its peak in L3 and L4 stages. L3 stage worms have sealed mouth and live without feeding, so it is very important for them to take full advantage of the inner sources of energy like fat to survive, while L4 worms are just starting to feed on blood; this gene may play important roles in non-feeding stages of H. contortus.
In this study, we demonstrated that Hc-daf-22 has similar gene structure with the SCPx protein genes, which all act as peroximal β-oxidation enzymes associated in the fatty acid metabolism. We also confirmed that Hc-daf-22 could partially recued the function of Ce-daf-22 in daf-22 mutant (ok693). It is safe to infer that Hc-daf-22 may play a similar role in the same pathway in H. contortus. Further studies are needed to confirm its enzyme activity as thiolase, its position in the peroximal β-oxidation pathway and its relationship with diapause formation in H. contortus in developmental level. The results of our study primary demonstrated the possible role of this gene and showed a promising future in the research of parasitic diapause formation mechanisms.
In this study, a new gene Hc-daf-22 was identified which is the homologue of Ce-daf-22 and human SCPx. Genome walking and RACE confirmed the full length gene of Hc-daf-22 (6,939 bp) which contained 16 exons separated by 15 introns, and encoded a cDNA of 1,602 bp. The 1,548 bp fragment upstream of the 5′-flanking region was confirmed to have promoter activity compared with 5′-flanking region of Ce-daf-22 by micro-injection. The rescue experiment indicated that Hc-daf-22 could partially rescue the function of Ce-daf-22. Furthermore, RNAi with Hc-daf-22 could partially silence the endogenous Ce-daf-22 in N2 worms and mimic the phenotype of daf-22 (ok693) mutants. Hc-daf-22 shared similar characteristics and function with Ce-daf-22 and may play an important role in peroxisomal β-oxidation and development in Haemonchus contortus.
Golden JW, Riddle DL. The Caenorhabditis elegans dauer larva: developmental effects of pheromone, food, and temperature. Dev Biol. 1984;102(2):368–78.
Sommerville RI, Davey KG. Diapause in parasitic nematodes: a review. Can J Zoo. 2002;80(11):1817–40.
Butcher RA, Ragains JR, Li W, Ruvkun G, Clardy J, Mak HY. Biosynthesis of the Caenorhabditis elegans dauer pheromone. Proc Natl Acad Sci U S A. 2009;106(6):1875–9.
Bun-ya M, Maebuchi M, Kamiryo T, Kurosawa T, Sato M, Tohma M, et al. Thiolase involved in bile acid formation. J Biochem-Tokyo. 1998;123(2):347–52.
Barnes TM, Jin Y, Horvitz HR, Ruvkun G, Hekimi S. The Caenorhabditis elegans behavioral gene unc-24 encodes a novel bipartite protein similar to both erythrocyte band 7.2 (stomatin) and nonspecific lipid transfer protein. J Neurochem. 1996;67(1):46–57.
Wirtz KWA. Phospholipid transfer proteins. Annu R Biochem. 1991;60(1):73–99.
Gallegos AM, Atshaves BP, Storey SM, Starodub O, Petrescu AD, Huang H, et al. Gene structure, intracellular localization, and functional roles of sterol carrier protein-2. Prog Lipid Res. 2001;40(6):498–563.
Ossendorp BC, Heusden GPH, Beer ALJ, Bos K, Schouten GL, Wirtz KWA. Identification of the cDNA clone which encodes the 58-kDa protein containing the amino acid sequence of rat liver non-specific lipid-transfer protein (sterol-carrier protein 2). Eur J Biochem. 1991;201(1):233–9.
Mannaerts GP, Van Veldhoven PP, Casteels M. Peroxisomal lipid degradation via β-and α-oxidation in mammals. Cell Biochem. 2000;32(1–3):73–87.
Kitamura T, Kobayashi S, Okada M. Regional expression of the transcript encoding sterol carrier protein x-related thiolase and its regulation by homeotic genes in the midgut of Drosophila embryos. Dev Growth Differ. 1996;38(4):373–81.
Wanders RJA, Denis S, Wouters F, Wirtz KWA, Seedorf U. Sterol carrier protein X (SCPx) is a peroxisomal branched-chain β-ketothiolase specifically reacting with 3-oxo-pristanoyl-CoA: a new, unique role for SCPx in branched-chain fatty acid metabolism in peroxisomes. Biochem Bioph Res Co. 1997;236(3):565–9.
Wanders RJA, Denis S, Van Berkel E, Wouters F, Wirtz KWA, Seedorf U. Identification of the newly discovered 58 kDa peroxisomal thiolase SCPx as the main thiolase involved in both pristanic acid and trihydroxycholestanoic acid oxidation: Implications for peroxisomal β-oxidation disorders. J Inherit Metab Dis. 1998;21(3):302–5.
Ghedin E, Wang S, Spiro D, Caler E, Zhao Q, Crabtree J, et al. Draft genome of the filarial nematode parasite Brugia malayi. Science. 2007;317(5845):1756–60.
Bairoch A. The PROSITE dictionary of sites and patterns in proteins, its current status. Nucleic Acids Res. 1993;21(13):3097.
Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths‐Jones S, et al. The Pfam protein families database. Nucleic Acids Res. 2004;32 Suppl 1:D138–41.
Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Mörck C, Pilon M. C. elegans feeding defective mutants have shorter body lengths and increased autophagy. BMC Dev Biol. 2006;6(1):1.
Joo H-J, Yim Y-H, Jeong P-Y, Jin Y-X, Lee J-E, Kim H, et al. Caenorhabditis elegans utilizes dauer pheromone biosynthesis to dispose of toxic peroxisomal fatty acids for cellular homoeostasis. Biochem J. 2009;422(1):61–71.
Wong A, Boutis P, Hekimi S. Mutations in the clk-1 gene of Caenorhabditis elegans affect developmental and behavioral timing. Genetics. 1995;139(3):1247–59.
O’Rourke EJ, Soukas AA, Carr CE, Ruvkun G. C. elegans major fats are stored in vesicles distinct from lysosome-related organelles. Cell Metab. 2009;10(5):430–5.
Timmons L, Court DL, Fire A. Ingestion of bacterially expressed dsRNAs can produce specific and potent genetic interference in Caenorhabditis elegans. Gene. 2001;263(1):103–12.
Zhang SO, Box AC, Xu N, Le Men J, Yu J, Guo F, et al. Genetic and dietary regulation of lipid droplet expansion in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 2010;107(10):4640–5.
Erich MS, Pasi KK, Bronwyn EC, et al. The genome and developmental transcriptome of the strongylid nematode Haemonchus contortus. Genome Biol. 2013;14:R89.
Ohba T, Rennert H, Pfeifer SM, He Z, Yamamoto R, Holt JA, et al. The structure of the human sterol carrier protein X/sterol carrier protein 2 gene (SCP2). Genomics. 1994;24(2):370–4.
Seedorf U, Raabe M, Ellinghaus P, Kannenberg F, Fobker M, Engel T, et al. Defective peroxisomal catabolism of branched fatty acyl coenzyme A in mice lacking the sterol carrier protein-2/sterol carrier protein-x gene function. Genes Dev. 1998;12(8):1189–201.
Möller G, Lüders J, Markus M, Husen B, Van Veldhoven PP, Adamski J. Peroxisome targeting of porcine 17β-hydroxysteroid dehydrogenase type IV/D-specific multifunctional protein 2 is mediated by its C-terminal tripeptide AKI. J Cell Biochem. 1999;73(1):70–8.
Gould SJ, Keller G-A, Hosken N, Wilkinson J, Subramani S. A conserved tripeptide sorts proteins to peroxisomes. J Cell Biol. 1989;108(5):1657–64.
Gould SJ, Keller G-A, Schneider M, Howell SH, Garrard LJ, Goodman JM, et al. Peroxisomal protein import is conserved between yeast, plants, insects and mammals. EMBO J. 1990;9(1):85.
Pfeifer SM, Sakuragi N, Ryan A, Johnson AL, Deeley RG, Billheimer JT, et al. Chicken sterol carrier protein m2/sterol carrier protein x: cDNA cloning reveals evolutionary conservation of structure and regulated expression. Arch Biochem Biophys. 1993;304(1):287–93.
Ohba T, Holt JA, Billheimer JT, Strauss Iii JF. Human sterol carrier protein x/sterol carrier protein 2 gene has two promoters. Biochemistry. 1995;34(33):10660–8.
Gillan V, Maitland K, McCormack G, Him NAIIN, Devaney E. Functional genomics of hsp-90 in parasitic and free-living nematodes. Int J Parasitol. 2009;39(10):1071–81.
Welz C, Krüger N, Schniederjans M, Miltsch SM, Krücken J, Guest M, et al. SLO-1-channels of parasitic nematodes reconstitute locomotor behaviour and emodepside sensitivity in Caenorhabditis elegans slo-1 loss of function mutants. PLoS Pathog. 2011;7(4):e1001330.
Glendinning SK, Buckingham SD, Sattelle DB, Wonnacott S, Wolstenholme AJ. Glutamate-gated chloride channels of Haemonchus contortus restore drug sensitivity to ivermectin resistant Caenorhabditis elegans. PLoS One. 2011;6(7):e22390.
Visser A, Geldhof P, De Maere V, Knox DP, Vercruysse J, Claerebout E. Efficacy and specificity of RNA interference in larval life-stages of Ostertagia ostertagi. Parasitology. 2006;133(06):777–83.
Lendner M, Doligalska M, Lucius R, Hartmann S. Attempts to establish RNA interference in the parasitic nematode Heligmosomoides polygyrus. Mol Biochem Parasitol. 2008;161(1):21–31.
Samarasinghe B, Knox DP, Britton C. Factors affecting susceptibility to RNA interference in Haemonchus contortus and in vivo silencing of an H11 aminopeptidase gene. Int J Parasitol. 2011;41(1):51–9.
Knox DP, Geldhof P, Visser A, Britton C. RNA interference in parasitic nematodes of animals: a reality check? Trends Parasitol. 2007;23(3):105–7.
Viney ME, Thompson FJ. Two hypotheses to explain why RNA interference does not work in animal parasitic nematodes. Int J Parasitol. 2008;38(1):43–7.
We would like to thank Dr. Yang CL (Institute of Genetics and Developmental Biology, Chinese Academy of Sciences), Dr. Wang XC (the National Institute of Biological Science, Beijing) and Dr. Ren CH (Academy of Military Medical Science, Beijing) for their assistance in transgenic technique in C. elegans and for donating of vectors we needed. The C. elegans daf-22 (ok693) strain was provided by the Caenorhabditis Genetics Center, University of Minnesota.
This project was supported by grants from the National Key Basic Research Program (973 Program) of China (No. 2015CB150300), the National Natural Science Foundation of China (Nos. 31101807 and 31372426), the Natural Science Foundation of Zhejiang Province (Y201327077), and the Fundamental Research Funds for the Central Universities.
Availability of data and material
The datasets supporting the conclusions of this article are included within the article. The complete cDNA sequence of Hc-daf-22 is submitted in the GenBank database under accession number HQ738470.1 and the amino acid sequence is also submitted in the GenBank under accession number AEO14647.1.
XG and HZ finished the isolation of Hc-daf-22, genome walking and RACE experiments, meanwhile, completed the functional rescue test of daf-22(ok693) mutants, these two authors had equal contribution to this research. XZ performed the relative quantification PCR of RNAi. QZ was responsible for gene structural and phylogenetic analysis. YY did the fat staining and promoter activity analysis. XC performed the statistical analysis of all data, AD was the corresponding author. All authors contributed to this manuscript and read and approved the submitted version.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate