Analysis and characterization of the genes encoding the Dicer and Argonaute proteins of Schistosoma japonicum
© Luo et al; licensee BioMed Central Ltd. 2010
Received: 11 August 2010
Accepted: 17 September 2010
Published: 17 September 2010
The Dicer and Argonaute(AGO) proteins within the small RNA regulatory pathways (SRRPs) play an indispensable role in regulation of gene expression. In this study, we analyzed two genes, Dicer and Argonaute, from Schistosoma japonicum, along with their expression through a combination of bioinformatics and experimental approaches.
Our results indicate that one Dicer and four Argonaute genes exist in Schistosoma japonicum, termed SjDicer and SjAGO1, 2, 3, and 4, respectively. SjDicer encodes 2590 amino acid residues that contains 5 conserved domains, including one amino-terminal helicase domain, one PAZ (Piwi-Argonaut-Zwille) domain, two RNAse III domains, and one dsRNA-binding domain. SjAGO1, 2, and 3 encode 1009, 945, and 904 amino acid residues, respectively, all of which contain PAZ and PIWI domains. In addition, we analyzed the expression profiles of SjDicer and SjAGO1 genes by qRT-PCR in eggs, miracidium, cercariae, schistosomula, and adult worms. Results showed consistent expression of both SjDicer and SjAGO1 in different stages; however, their expression levels were stage-dependent, with the highest being in the miracidium stage.
This study provided the sequence of the Dicer and Ago genes of S. japonicum and their expression profiles which are essential for further investigation of functions of miRNA in Schistosoma japonicum.
Small RNA-mediated gene silencing pathways play important and diverse roles in development and differentiation of organisms[1–5]. Dicer and Argonaute (AGO) are the two core proteins involved in this pathway . Dicer processes linear dsRNA into small interfering RNA (siRNA) duplexes and also excises mature miRNAs from pre-miRNAs in the cytoplasm . Mature miRNAs and siRNAs bind to distinct members of the Argonaute protein family to form ribonucleoprotein complexes that recognize partially, or nearly perfect, complementary nucleic acid targets, and then mediate a variety of regulatory processes, including transcriptional and post-transcriptional gene silencing [8–10].
Dicer belongs to the Ribonuclease III family and has many subtypes, all of which contain one PAZ (Piwi-Argonaut-Zwille) domain, two RNAse III domains, and one RNA-binding domain (RBD) . Besides the PAZ domain, all Argonaute proteins share the structural features of PIWI domains. Both Dicer and Argonaute proteins are highly conserved between species, and many organisms encode multiple members of the family. Although those subtypes are similar to each other in amino acids sequence and structure, they participate in different small RNA regulatory pathways (SRRPs)[6, 14]. In general, in miRNA pathway Dicer-1 cleaves pre-miRNAs to form mature miRNAs, which are transferred to Ago1 for silencing of target mRNAs[3, 15, 16]. In siRNA pathway, long double-stranded RNA is cleaved by Dicer-2 to form mature siRNA and then complexes with a RISC-like structure for silencing of target mRNA.
Human schistosomiasis is one of the most prevalent and serious parasitic diseases in tropical and subtropical regions[17–19]. Schistosomiasis is mainly caused by species of Schistosoma including S. japonicum, S. mansoni, and S. haematobium. Recently, Gomes and Krautz-peterson and Skelly reported the discovery of Dicer and the Ago family using bioinformatics and experimental approaches in S. mansoni. We have previously identified small regulatory RNA in S. japonicum, including miRNAs, by traditional cloning and high-throughput sequencing approaches[22, 23]. To further investigate potential functions of the small RNAs and their role in the regulatory networks in S. japonicum, it is meaningful and important to analyse and characterize the important proteins, Dicer and Ago. In this study, we obtained the sequence of the Dicer and Ago genes in S. japonicum through a combination of bioinformatics and experimental methods and analyzed their expression profiles in different stages of the parasite.
Materials and methods
Sequence retrieval of Dicer and Argonaute families
S. japonicum genome and transcriptome sequence data were downloaded from SCBIT (http://www.scbit.org/pages/index.do), as well as the Sanger Institute (http://www.sanger.ac.uk/) and NCBI (http://www.ncbi.nlm.nih.gov/). Amino acid sequences of D. melanogaster and S. mansoni orthologs were used as queries. The BLASTp algorithm, underpinned by the Pfam and CDD databases was used for searches of conserved protein domains or motifs.
Multiple alignments and phylogenetic analyses
Multiple alignments of protein sequence were performed by ClustalX and phylogenetic analyses were conducted in MEGA 4. Phylogenetic trees of these sequences were inferred using the Neighbor-Joining method. The bootstrap consensus tree inferred from 1000 replicates was used to represent the evolutionary history of the taxa analyzed. Branches corresponding to partitions reproduced in less than 50% bootstrap replicates are collapsed. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches . The tree was drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. All positions containing gaps and missing data were eliminated from the dataset.
Parasite culturing was performed as described previously. S. japonicum was maintained by routine passage through Oncomelania hupensis snails and BALB/c mice. The infected snails were induced to shed cercariae under light exposure for 2 h, and the cercariae were recovered by sedimentation on ice. Schistosomula or adult worms were obtained by liver perfusion of mice after 3 or 6 wks of infection. Eggs were obtained from liver of mice 6 wks after infected by S. japonicum. Eggs were then incubated to produce miracidium at room temperature and normal light conditions. After collection, all freshly isolated samples were washed three times with 1 × phosphate buffered saline (PBS), pH 7.4, and were immediately used for extraction of total RNA or stored in liquid nitrogen. All procedures performed on animals within this study were conducted in accordance with and by approval of the Tongji University Committee on Use and Care of Animals.
PCR amplification and RACE (Rapid Amplification of cDNA Ends)
Total RNA was extracted using Trizol reagent. To make sure there was no contamination by S. Japonicum DNA, total RNA was treated with RNase-free DNase I. And then total RNA was quantified using a spectrophotometer, and 1 μg of RNA was reverse transcribed using an oligo-dT primer. Reverse-transcribed cDNA samples were used as templates for PCR amplification. The 5'UTR and 3'UTR of genes were amplified via RACE (rapid amplification of cDNA ends) using the Full RACE kit from TaKaRa, and target sequences were then cloned and sequenced.
Real-time PCR analysis
Eggs, cercariae, schistosomula, and adult worms were collected for real-time PCR to quantify developmental expression of SjDicer and SjAGO. Reverse-transcribed cDNA was prepared as above-mentioned and was used as templates for real-time PCR amplification using SYBR Green (TaKaRa) and a 7300 Real-Time PCR System (Applied Biosystems). Specific primers for S. japonicum GAPDH were used as a control. Relative expression was calculated according to the 2-ΔΔCt method.
Data mining of Dicer- and Argonaute-related sequences in S. japonicum
Homology analyses at the amino acid level revealed that the members of small RNA regulatory pathways (SRRPs) are well conserved among diverse organisms of D. melanogaster, C. elegans, Homo sapiens, Mus musculus and S. Mansoni. Due to the closer relationship between S. japonicum and S. mansoni, we chose SmDicer and SmAGO1/2/3/4 as reference sequences to retrieve their orthologs in S. japonicum.
Transcriptome analysis via BLAST of S. japonicum revealed one EST (AY810618) of 538 bp on supercontig CCON0000000957 that exhibited sequence similarity to proteins belonging to the Dicer family. Based on data mining of the transcriptome database of S. japonicum, we found two mRNA sequences (IDs: AY814744 and AY808628) with similarity to AGO1, three mRNA sequences (IDs: FN330861, AY809895, and AY809756) with similarity to AGO2, one mRNA sequence (ID: Sjc_0103990) with similarity to AGO3, and one protein sequence (ID: sjp_0117000) with similarity to AGO4. Therefore, we suggest that S. japonicum contains four members of the Argonaute protein family. These genes were named SjAGO1, SjAGO2, SjAGO3, SjAGO4, and SjDicer, respectively; however, all of them are partial coding sequences but SjAGO3. The complete CDS of SjAGO3 gene is 2715 bp coding a 904 aa protein in size. Therefore, we further obtained the remaining sequences of these genes through experimental approaches.
Exploring the entire coding sequences of SjDicer in the genome
Oligonucleotides used to characterize SjDicer and SjAGO1/2
The length of PCR productions
Real time F
Real time R
Real time F
Real time R
Real time F
Real time R
Analysis of conserved protein domains in SjDicer
The BLASTp algorithm, underpinned by the Pfam (v22.0) and CDD databases, was used for searches of conserved protein domains or motifs from SjDicer sequences. Results revealed that SjDicer contains an amino-terminal helicase domain, a DUF283 domain, a PAZ domain, two tandom RNAse III domains, and a dsRNA-binding domain (Fig. 1).
Signature motif in the RNAse III and four catalytic residues of SjDicer protein
ERLE TIGD......D CVE
QRLE FLGD.....D IFE
ERME TIGD......D CVE
QRLE FLGD.....D IFE
ERFE TIGD......D AVE
QRLE FLGD.....D IFE
ERLE TIGD......D CVE
QRLE FLGD.....D VFE
ERLE TIGD......D CVE
QRLE FLGD.....D VFE
ERLE MLGD......D CVE
QRLE FLGD.....D IFE
Exploring the entire coding sequences of SjAGO1 and SjAGO2 in the genome
AGO1 and AGO2 play a pivotal role in small RNA regulatory pathways (SRRPs). Therefore, our study focused mainly on SjAGO1 and SjAGO2. Transcriptome analysis of S. japonicum demonstrated that two ESTs, AY814744 and AY808628, were highly similar to the upstream and downstream regions of the SmAGO1 gene. At the same time, these analyses easily identified the 3' end of SjAGO1, since the extreme 3' end of SmAGO1 matched the sequence of AY814744. However, the extreme 5' end of the SjAGO1 cDNA is not well conserved and could not be identified by analysis of the genomic contig sequences. Moreover, there was no overlap between AY814744 and AY808628. Therefore, we designed primers based on the two known sequences, in order to amplify overlapping fragments of the predicted SjAGO1 coding gene. The 5' end of the SjAGO1 was cloned by RACE(Table 1: prmier 7, 8 and 9). Finally, the full SjAGO1 gene sequence was obtained after spliced sequences were concatenated. All exons combine to form a 3266 bp ORF, potentially encoding a 1009-aa protein with a predicted molecular weight of 111,647 Da.
Analysis of conserved protein domains in SjAGO
Analysis of expression profiles of SjDicer and SjAGO in different stages of the parasite by Real-Time PCR
In S. japonicum, miRNAs have been identified by traditional approaches or high-throughput sequencing[22, 23]. Meanwhile, endo-siRNAs were also found in S. japonicum via Solexa technology. All of these imply the presence of constituents of the small RNA regulatory pathways in S. japonicum. Through a combination of phylogenetic tree, similarity between domains and mining of the S. japonium database we were able to show putative sequences with conserved domains of the Small RNA-mediated gene silencing pathways in the parasite. The results revealed that there was a single Dicer and four members of the AGO family in S. japonicum, named SjDicer and SjAGO1, 2, 3, and 4, respectively. However, all of these are partial coding sequences. Our study mainly focused on SjDicer and SjAGO1/2 in an effort to elucidate their full-length gene sequences through experimental approaches.
Our data indicate that the SjDicer gene comprises 28 exons that potentially encode a 2590 amino acid protein. Like SmDicer, SjDicer protein contains all domains that are characteristic of metazoan dicers including an amino terminal helicase domain, DUF283, a PAZ domain, two RNAse III domains and an RNA binding domain[20, 21]. The identity between the SmDicer and SjDicer is as much as 71.8% in protein level. Among these domains, RNase III domains specifically cleave one strand of dsRNA. The first RNAse III domain is predicted to cut the RNA strand bearing a 3'-hydroxy group approximately 21 nucleotides from the end. In conjunction with the PAZ domain, the RNAse III domain is thought to be responsible for determining the distance from the terminus of the RNA to the cleavage site. Dicer proteins are key components of the biogenesis of small RNA such as miRNA and siRNA. For example, Lee et al. showed that Drosophila Dcr-1 mutants were defective in processing miRNA precursors, whereas Dcr-2 mutants had reduced siRNA levels and a complete RNAi defect in the eye.
Different Ago proteins have been found in several organisms, ranging from a single member in Schizosaccharomyces pombe, to more than twenty in C. Elegans[13, 31]. Alignments between SjAGO1 and SmAGO1(SjAGO2 and SmAGO2) proteins revealed 81.1%(81.2%) identity. The role of Ago family in Drosophila and in C. elegans development is well documented for D. melanogaster and C. elegans. It was observed that Ago1/2 are primary effector proteins responsible for the small RNA silencing pathway in fly and in C. Elegans. For example, Okamura et al. showed that Drosophila AGO1 is essential for mature miRNA production, whereas AGO2 is responsible for siRNA-directed RNAi. The PAZ domain of Argonaute is a small RNA-binding domain, while the PIWI domain is an RNase-H-type domain that relies on divalent cation binding to facilitate dsRNA-guided cleavage of ssRNA .
Most small-RNA regulatory pathways are expressed in a developmental or organ-specific manner, or both, which provides information regarding their functions. As mentioned, SjDicer and SjAGO1 are crucial proteins in the miRNA regulatory pathway. Therefore, we analyzed the expression of SjDicer and SjAGO1 by relative qRT-PCR in different developmental stages of S. japonicum. The expression patterns of SjDicer and SjAGO1 during the life cycle of S. japonicum indicate that the miRNA regulatory pathway might take part in the transformation and development of S. japonicum.
In summary, This study demonstrated existence of Dicer and Argonaute genes in S. japonicum and the 2590 amino acid residues of SjDicer containing 5 conserved domains, and the 1009, 945, and 904 amino acid residues of the SjAGO1, 2, and 3, respectively. In addition, expression profiles of SjDicer and SjAGO1 genes showed consistent expression of both genes in different stages, but difference in their expression levels. The study provides a theoretical platform for exploration of the functions of siRNA and miRNA in S. japonicum. Further investigations should be considered to explore the functions of the SjDicer and SjAGO families through experimental approaches.
This work is supported by grants from the National Basic Research Program (973 Program) in China (2007CB513100).
- Keene KM, Foy BD, Sanchez-Vargas I, Beaty BJ, Blair CD, Olson KE: RNA interference acts as a natural antiviral response to O'nyong-nyong virus (Alphavirus; Togaviridae) infection of Anopheles gambiae. Proc Natl Acad Sci USA. 2004, 101 (49): 17240-17245. 10.1073/pnas.0406983101.PubMed CentralView ArticlePubMedGoogle Scholar
- Nakahara K, Kim K, Sciulli C, Dowd SR, Minden JS, Carthew RW: Targets of microRNA regulation in the Drosophila oocyte proteome. Proc Natl Acad Sci USA. 2005, 102 (34): 12023-12028. 10.1073/pnas.0500053102.PubMed CentralView ArticlePubMedGoogle Scholar
- Olsen PH, Ambros V: The lin-4 regulatory RNA controls developmental timing in Caenorhabditis elegans by blocking LIN-14 protein synthesis after the initiation of translation. Dev Biol. 1999, 216 (2): 671-680. 10.1006/dbio.1999.9523.View ArticlePubMedGoogle Scholar
- Saito K, Nishida KM, Mori T, Kawamura Y, Miyoshi K, Nagami T, Siomi H, Siomi MC: Specific association of Piwi with rasiRNAs derived from retrotransposon and heterochromatic regions in the Drosophila genome. Genes Dev. 2006, 20 (16): 2214-2222. 10.1101/gad.1454806.PubMed CentralView ArticlePubMedGoogle Scholar
- Vagin VV, Sigova A, Li C, Seitz H, Gvozdev V, Zamore PD: A distinct small RNA pathway silences selfish genetic elements in the germline. Science. 2006, 313 (5785): 320-324. 10.1126/science.1129333.View ArticlePubMedGoogle Scholar
- Campbell Corey: Comparative genomics of small RNA regulatory pathway components in vector mosquitoes. BMC Genomics. 2008, 9: 425-10.1186/1471-2164-9-425.PubMed CentralView ArticlePubMedGoogle Scholar
- Nelson Peter, Kiriakidou Marianthi, Sharma Anup, Elsa Maniataki, Zissimos Mourelatos: The microRNA world: small is mighty. TRENDS in Biochemical Sciences. 2003, 28 (10): 10.1016/j.tibs.2003.08.005.Google Scholar
- Song J, Joshua-Tor L: Argonaute and RNA -- getting into the groove. Curr Opin Struck Biol. 2006, 16: 5-11. 10.1016/j.sbi.2006.01.010.View ArticleGoogle Scholar
- Liu J, Valencia-Sanchez MA, Hannon GJ, Parker R: MicroRNA dependent localization of targeted mRNAs to mammalian P-bodies. Nat Cell Biol. 2005, 7 (7): 719-723. 10.1038/ncb1274.PubMed CentralView ArticlePubMedGoogle Scholar
- Souret FF, Kastenmayer JP, Green PJ: AtXRN4 degrades mRNA in Arabidopsis and its substrates include selected miRNA targets. Mol Cell. 2004, 15 (2): 173-183. 10.1016/j.molcel.2004.06.006.View ArticlePubMedGoogle Scholar
- Macrae IJ, Zhou K, Li F, Repic A, Brooks AN, Cande WZ, Adams PD, Doudna JA: Structural basis for double-stranded RNA processing by Dicer. Science. 2006, 311 (5758): 195-198. 10.1126/science.1121638.View ArticlePubMedGoogle Scholar
- Cerutti L, Mian N, Bateman A: Domains in gene silencing and cell differentiation proteins: the novel PAZ domain and redefinition of the Piwi domain. Trends Biochem Sci. 2000, 25 (10): 481-482. 10.1016/S0968-0004(00)01641-8.View ArticlePubMedGoogle Scholar
- Carmell MA, Xuan Z, Zhang MQ, Hannon GJ: The Argonaute family: tentacles that reach into RNAi, developmental control, stem cell maintenance, and tumorigenesis. Genes Dev. 2002, 16: 2733-2742. 10.1101/gad.1026102.View ArticlePubMedGoogle Scholar
- He L, Hannon GJ: MicroRNAs: Small RNAs with a big role in gene regulation. Nature Rev Genet. 2004, 5: 522-531. 10.1038/nrg1379.View ArticlePubMedGoogle Scholar
- Behm-Ansmant I, Rehwinkel J, Doerks T, Stark A, Bork P, Izaurralde E: mRNA degradation by miRNAs and GW182 requires both CCR4:NOT deadenylase and DCP1:DCP2 decapping complexes. Genes Dev. 2006, 20 (14): 1885-1898. 10.1101/gad.1424106.PubMed CentralView ArticlePubMedGoogle Scholar
- Ketting RF, Fischer SE, Bernstein E, Sijen T, Hannon GJ, Plasterk RH: Dicer functions in RNA interference and in synthesis of small RNA involved in developmental timing in C. elegans. Genes Dev. 2001, 15 (20): 2654-2659. 10.1101/gad.927801.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang L, Utzinger J, Zhou XN: Schistosomiasis control: experiences and lessons from China. Lancet. 2008, 372 (9652): 1793-1795. 10.1016/S0140-6736(08)61358-6.View ArticlePubMedGoogle Scholar
- King CH, Dickman K, Tisch DJ: Reassessment of the cost of chronic helmintic infection: a meta-analysis of disability-related outcomes in endemic schistosomiasis. Lancet. 2005, 365 (9470): 1561-1569. 10.1016/S0140-6736(05)66457-4.View ArticlePubMedGoogle Scholar
- van der Werf MJ, de Vlas SJ, Brooker S, Looman CW, Nagelkerke NJ: Quantification of clinical morbidity associated with schistosome infection in sub-Saharan Africa. Acta Trop. 2003, 86: 125-10.1016/S0001-706X(03)00029-9.View ArticlePubMedGoogle Scholar
- Gomes MS: Preliminary analysis of miRNA pathway in Schistosoma mansoni. Parasitol Int. 2008Google Scholar
- Krautz-Peterson G, Skelly PJ: Schistosoma mansoni: The dicer gene and its expression. Exp Parasitol. 2007Google Scholar
- Xue X, Sun J, Zhang Q, Wang Z, Huang Y: Identification and Characterization of Novel MicroRNAs from Schistosoma japonicum. PLoS ONE. 2008, 3 (12): e4034-10.1371/journal.pone.0004034.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Z, Xue X, Sun J, Luo R, Xu X: An ''In-Depth'' Description of the Small Non-coding RNA Population of Schistosoma japonicum Schistosomulum. PLoS Negl Trop Dis. 2010, 4 (2): e596-10.1371/journal.pntd.0000596.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-8. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-25.PubMedGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-91. 10.2307/2408678.View ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001, 25: 402-8. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Cordin O, Banroques J, Tanner Nk: The DEAD-box protein family of RNA helicase. Gene. 2006, 367: 17-37. 10.1016/j.gene.2005.10.019.View ArticlePubMedGoogle Scholar
- Zhang H, Kolb FA, Jaskiewicz L, Westhof E, Filipowicz W: Single processing center models for human Dicer and bacterial RNase III. Cell. 2004, 118: 57-68. 10.1016/j.cell.2004.06.017.View ArticlePubMedGoogle Scholar
- Lee YS, Nakahara K, Pham JW, Kim K, He Z, Sontheimer EJ, Carthew RW: Distinct roles for Drosophila Dicer-1 andDicer -2 in the siRNA/miRNA silencing pathways. Cell. 2004, 117: 69-81. 10.1016/S0092-8674(04)00261-2.View ArticlePubMedGoogle Scholar
- Yigit E, Batista PJ, Bei Y, Pang KM, Chen CG, Tolia NH: Analysis of the C. elegans Argonaute family reveals that distinct Argonautes act sequentially during RNAi. Cell. 2006, 127: 747-57. 10.1016/j.cell.2006.09.033.View ArticlePubMedGoogle Scholar
- Farazi Thalia, Juranek Stefan, Tuschl Thomas: The growing catalog of small RNAs and their association with distinct Argonaute/Piwi family members. Development. 2008, 135: 1201-1214. 10.1242/dev.005629.View ArticlePubMedGoogle Scholar
- Okamura K, Ishizuka A, Siomi H, Siomi MC: Distinct roles for Argonaute proteins in small RNA-directed RNA cleavage pathways. Genes Dev. 2004, 18: 1655-1666. 10.1101/gad.1210204.PubMed CentralView ArticlePubMedGoogle Scholar
- Joshua-Tor L: The Argonautes. Cold Spring Harb Symp Quant Biol. 2006, 71: 67-72. 10.1101/sqb.2006.71.048.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.