Comparative analysis of microRNA profiles between adult Ascaris lumbricoides and Ascaris suum
© Shao et al.; licensee BioMed Central Ltd. 2014
Received: 3 February 2014
Accepted: 17 April 2014
Published: 27 April 2014
The parasitic nematodes Ascaris lumbricoides and A. suum are of great public health and economic significance, and the two taxa were proposed to represent a single species. miRNAs are known with functions of gene regulations at post-transcriptional level.
We herein compared the miRNA profiles of A. lumbricoides and A. suum female adults by Solexa deep sequencing combined with bioinformatics analysis and stem-loop real-time PCR. Using the A. suum genome as the reference genome, we obtained 171 and 494 miRNA candidates from A. lumbricoides and A. suum, respectively. Among which, 74 miRNAs were shared between the two taxa, 97 and 420 miRNAs were A. lumbricoides and A. suum specific. Target and function prediction revealed a significant set of targets which are related to ovarian message protein, vitellogenin and chondroitin proteoglycan of the two nematodes. Enrichment analysis revealed that the percentages of most predicted functions of the miRNA targets were similar, with some taxon specific or taxon enhanced functions, such as different target numbers, specific functions (NADH dehydrogenase and electron carrier functions), etc.
This study characterized comparatively the miRNAs of adult A. lumbricoides and A. suum, and the findings provide additional evidence that A. lumbricoides and A. suum represent a single species. Due to the fast evolution nature of miRNAs and the different parasitic living conditions of humans and pigs, the phenomenon above might indicate a fast evolution of miRNAs of Ascaris in humans and pigs.
KeywordsMicroRNA (miRNA) Ascaris lumbricoides Ascaris suum Comparative analysis
Roundworms Ascaris lumbricoides and A. suum are important parasites of human and pig health and socio-economic significance, with global distribution [1, 2]. Ascaris infects approximately 1.2 billion people globally and has been associated with intestinal pathology, respiratory symptoms and malnutrition in children from endemic areas . It was reported that humans and pigs can both be infected with the two nematodes . There are very limited nucleotide differences in the sequences of internal transcribed spacers (ITS) of ribosomal DNA between the two ascarid species . A. lumbricoides and A. suum are now proposed to represent a single species, and A. suum is considered a synonym of A. lumbricoides . The genome and transcripts of A. suum were recently available, which provided valuable resources for better understanding and further studies of the biology of the Ascaris parasites [6–8].
MicroRNAs (miRNAs) are non-coding small RNA of 18–24 nt in length. They are considered as key regulators for gene expression at the post-transcriptional level [9–12]. Due to their key regulating functions in growth, metabolism, development and cell differentiation and their ability to respond to environmental and developmental signals, miRNA is essential for the complex life cycles of human and animal parasites [13–22]. miRNAs are also important for pathogen-host interactions [23–26]. Recently studies have indicated that miRNAs may represent potential biomarkers for various biomedical problems, such as the differentiation of different tumors and muscular tissues [27–29], and they may provide biomarkers for the characterization of different genotypes of Toxoplasma gondii .
There might be differences in gene expression and regulation given the different living environments of A. lumbricoides and A. suum, although they are considered to represent the same species. Therefore, it would be interesting to characterize the expression profiles of miRNAs in the two taxa. Herein the objective of the present study was to examine and compare the miRNA profiles of A. lumbricoides and A. suum using an integrative approach combining Solexa deep sequencing combined with bioinformatics analysis and stem-loop real-time PCR.
The present study was approved by the Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (Approval No. LVRIAEC2011-006), and the A. lumbricoides and A. suum samples were collected strictly according to the requirements of the Ethics Procedures and Guidelines of the People's Republic of China.
Female adults of A. suum roundworms were obtained from slaughtered pigs in Shenzhen, China. Female adults of A. lumbricoides were obtained from a patient with ascariasis after being treated with piperazine in Zhanjiang, Guangdong Province, China. Worms were incubated in physiological saline for 3 h at 37°C and then washed 3 times to get rid of contamination from hosts. Female adults were identified morphologically and their identity was further ascertained by PCR amplification and sequencing of the first internal transcribed spacer (ITS-1) of rDNA . The nematodes were stored at −80°C until use.
Total RNA isolation, small RNA preparation and deep sequencing
Total RNA was prepared from a whole single adult of A. lumbricoides and A. suum respectively, using Trizol reagent (Invitrogen) according to the manufacturer’s protocol. Ten μg total RNA and Novex 15% TBE-Urea gel (Invitrogen) were used for small RNA isolation. The RNA fragments of 18–30 bases long were added with 5’ and 3’ adaptors (Illumina) to the both ends, reverse transcripted with a RT-PCR kit (Invitrogen), and sequenced employing a Solexa sequencer at Huada Genomics Institute Co. Ltd, China.
Adaptors, low quality sequences and reads smaller than 18 nt were removed from the raw sequencing data. The reads were then searched against GenBank and Rfam database (http://www.sanger.ac.uk/software/Rfam) to identify non-coding RNA. The remaining reads were mapped onto the A. suum genome by using the SOAP  with the sequences of pre-miRNA meeting the three criteria: 1) there was a standard stem-loop structure of pre-miRNA; 2) mature miRNAs were present in one arm instead of the loop of hairpin precursors; and 3) the free energy hybridization was lower than −18 kcal/mol. The identified miRNA candidates were searched against the known miRNAs of A. suum deposited in the Sanger miRBase with Blast . Unmatched miRNA candidates were marked as novel miRNA. Targets of miRNA candidates were predicated with RNAhybrid  with the following extra parameters: A) the △△G was set as lower than −25 kcal/mol; B) P-value was set as ≤ 0.01. The Gene Ontology database (GO, http://www.geneontology.org/), Blast and Interproscan (http://www.ebi.ac.uk/Tools/pfa/iprscan/) were used for prediction of functions of predicted targets.
Analysis of novel miRNA expression
The representative novel miRNA in the two nematodes were certified using stem-loop real-time reverse transcription polymerase chain reaction (RT-PCR) with SYBR Green . The house keeping gene β-actin of A. suum (GenBank accession no. BI594141) was used as the endogenous control with primers as follows: forward primer (5′-CTCGAAACAAGAATACGATG-3′) and reverse primer (5′-ACATGTGCCGTTGTATGATG-3′). Primers were synthesized by Shenggong Co, Ltd., China. The cycle conditions were as follows: 94°C for 30s, 52°C for 30s, 72°C for 30s and finally with a single extension at 72°C for 10 min. The quantification of each miRNA relative to β-actin gene was calculated using the equation: N = 2-ΔCt, ΔCt = CtmiRNA-Ctacin .
Profile differences in short RNAs between the two Ascaristaxa
Common and taxon-specific reads of Ascaris lumbricoides and Ascaris suum
Unique sRNAs (%)
Total sRNAs (%)
A. lumbricoides specificb
A. suum specificc
miRNA profiles of the two taxa
Comparison of miRNA profiles in Ascaris lumbricoides (Alu) and A. suum (Asu)
Target prediction and the functional prediction of the predicted targets
A total of 57,359 mRNA and EST items of Ascaris deposited in NCBI were downloaded and used for target prediction. Under the stringent matching criteria, it was found that the target numbers of both nematodes ranged from one to thousands. For A. lumbricoides, the target number ranged from one (Alu-miR-novel-012-3p, Alu-miR-novel-066-5p, Alu-miR-novel-102-3p) to 3,562 (Alu-miR-novel-063-3p), 255 in average. For A. suum, the target number ranged between one (asu-miR-novel-027-3p, asu-miR-novel-044-5p, and asu-miR-novel-068-5p) and 3,343 (asu-miR-novel-039-3p), 168 in average.
Functional prediction revealed a significant set of targets which are related to ovarian message protein in A. lumbricoides (n = 9) and A. suum (n = 10), a set of vitellogenin (5 in A. lumbricoides and 18 in A. suum), and a set of chondroitin proteoglycan (5 in A. lumbricoides and 7 in A. suum). The same phenomenon is found in NADH dehydrogenase (n = 7 in A. suum, n = 3 in A. lumbricoides), 40s ribosomal protein (n = 14 in A. suum, n = 2 in A. lumbricoides), and 60s ribosomal protein (n = 11 in A. suum, n = 7 in A. lumbricoides). As a distinguished character, there were more targets related to movement in A. suum, such as actin (n = 5) and tubulin (n = 6), while such targets were very few in A. lumbricoides.
Quantification of A. lumbricoides and A. suummiRNAs
Both of the mature miRNAs of the 6 selected miRNAs can be successfully detected, except the alu-miR-novel-64-5p and asu-miR-novel-383-3p. Two miRNAs named asu-miR-novel-031 and alu-miR-novel-53 had significant expression differences at 5p and 3p. For alu-miR-novel-53, the expression levels were 1.00 ± 0.12 at 5p, and it was 1.68 ± 0.17 at 3p. And for asu-miR-novel-031, the relative expression level was 1.85 ± 0.24 at 5p, while it was 1.56 ± 0.30 at 3p. The alu-miR-novel-64-5p and asu-miR-novel-383-3p were not successfully amplified, which could be attributed to the very low expression levels, despite that the primers and reaction mixture were modified several times.
A previous study revealed that A. lumbricoides and A. suum have identical 5.8S and ITS-2 rDNA sequences . There were only six (1.3%) nucleotide differences in ITS-1. All these genetic evidences supported the proposal that A. lumbricoides and A. suum represent a single species [4, 5].
Due to their specific expression in different organisms, tissues and cells, miRNAs may provide potential novel biomarkers [27–29]. Our previous study indicated that different T. gondii genotypes have their unique miRNA profiles . In the present study, the two nematode taxa shared very high percentage of total reads (72.47%), while having very small percentage of unique reads (4.35%), and the total percentage of ncRNA was at similar level (46.11% and 39.97% in A. lumbricoides and A. suum, respectively), which indicated huge redundance of some ncRNA of the two nematodes. When parasitic environments change, the gene expression of parasites will be correspondingly changed, and regulators such as miRNAs will also be changed, therefore the death and new creation of miRNAs are very common . For the different parasitic environment of pigs and humans, adaptive modifications should have occurred. So although the sRNA/shared-miRNAs for fundamental metabolism were similar in the two taxa, there still are some specific miRNAs in each nematode taxon. We obtained 171 and 494 miRNA candidates as expression profiles, with 74 miRNAs being shared by the two nematode taxa, and 97 and 420 miRNAs being A. lumbricoides and A. suum specific. The same precursors of a miRNA can be found in the genome of A. lumbricoides or A. suum, however, we obtained fewer miRNAs from A. lumbricoides than from A. suum, which indicated that some miRNAs in A. lumbricoides were not expressed. This difference might be resulted from the different parasitic living environment in pigs and humans.
Function prediction and enrichment analysis showed that targets of the miRNAs of the two parasites had similar metabolism patterns, including cellular component, molecular function and biological process, with some specific differences, such as different target number, specific functions (NADH dehydrogenase and electron carrier functions), etc.
The present study characterized comparatively the miRNAs of adult A. lumbricoides and A. suum, and the findings support the recent proposal that A. lumbricoides and A. suum represent a single species . Due to the fast evolution nature of miRNAs and the different parasitic living conditions of humans and pigs, the phenomenon above might indicate a fast evolution of miRNAs of Ascaris.
This work was supported in part by the International Science & Technology Cooperation Program of China (Grant No. 2013DFA31840), the Science Fund for Creative Research Groups of Gansu Province (Grant No. 1210RJIA006) to XQZ, and the China Postdoctoral Science Foundation (Grant No. 201104363) to MJX.
- Dold C, Holland CV: Ascaris and ascariasis. Microbes Infect. 2011, 13: 632-637. 10.1016/j.micinf.2010.09.012.View ArticlePubMedGoogle Scholar
- Crompton DW: Ascaris and ascariasis. Adv Parasitol. 2001, 48: 285-375.View ArticlePubMedGoogle Scholar
- Hagel I, Giusti T: Ascaris lumbricoides: an overview of therapeutic targets. Infect Disord Drug Targets. 2010, 10: 349-367. 10.2174/187152610793180876.View ArticlePubMedGoogle Scholar
- Zhu X, Chilton NB, Jacobs DE, Boes J, Gasser RB: Characterisation of Ascaris from human and pig hosts by nuclear ribosomal DNA sequences. Int J Parasitol. 1999, 29: 469-478. 10.1016/S0020-7519(98)00226-4.View ArticlePubMedGoogle Scholar
- Leles D, Gardner SL, Reinhard K, Iniguez A, Araujo A: Are Ascaris lumbricoides and Ascaris suum a single species?. Parasit Vectors. 2012, 5: 42-10.1186/1756-3305-5-42.PubMed CentralView ArticlePubMedGoogle Scholar
- Jex AR, Liu S, Li B, Young ND, Hall RS, Li Y, Yang L, Zeng N, Xu X, Xiong Z, Chen F, Wu X, Zhang G, Fang X, Kang Y, Anderson GA, Harris TW, Campbell BE, Vlaminck J, Wang T, Cantacessi C, Schwarz EM, Ranganathan S, Geldhof P, Nejsum P, Sternberg PW, Yang H, Wang J, Wang J, Gasser RB: Ascaris suum draft genome. Nature. 2011, 479: 529-533. 10.1038/nature10553.View ArticlePubMedGoogle Scholar
- Cantacessi C, Zou FC, Hall RS, Zhong W, Jex AR, Campbell BE, Ranganathan S, Sternberg PW, Zhu XQ, Gasser RB: Bioinformatic analysis of abundant, gender-enriched transcripts of adult Ascaris suum (Nematoda) using a semi-automated workflow platform. Mol Cell Probes. 2009, 23: 205-217. 10.1016/j.mcp.2009.03.003.View ArticlePubMedGoogle Scholar
- Huang CQ, Gasser RB, Cantacessi C, Nisbet AJ, Zhong W, Sternberg PW, Loukas A, Mulvenna J, Lin RQ, Chen N, Zhu XQ: Genomic-bioinformatic analysis of transcripts enriched in the third-stage larva of the parasitic nematode Ascaris suum. PLoS Negl Trop Dis. 2008, 2: e246-10.1371/journal.pntd.0000246.PubMed CentralView ArticlePubMedGoogle Scholar
- Du T, Zamore PD: Beginning to understand microRNA function. Cell Res. 2007, 17: 661-663. 10.1038/cr.2007.67.View ArticlePubMedGoogle Scholar
- Lim LP, Glasner ME, Yekta S, Burge CB, Bartel DP: Vertebrate microRNA genes. Science. 2003, 299: 1540-10.1126/science.1080372.View ArticlePubMedGoogle Scholar
- Wienholds E, Plasterk RH: MicroRNA function in animal development. FEBS Lett. 2005, 579: 5911-5922. 10.1016/j.febslet.2005.07.070.View ArticlePubMedGoogle Scholar
- Finnegan EF, Pasquinelli AE: MicroRNA biogenesis: regulating the regulators. Crit Rev Biochem Mol Biol. 2013, 48: 51-68. 10.3109/10409238.2012.738643.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu Q, Tuo W, Gao H, Zhu XQ: MicroRNAs of parasites: current status and future perspectives. Parasitol Res. 2010, 107: 501-507. 10.1007/s00436-010-1927-6.View ArticlePubMedGoogle Scholar
- Lin WC, Li SC, Lin WC, Shin JW, Hu SN, Yu XM, Huang TY, Chen SC, Chen HC, Chen SJ, Huang PJ, Gan RR, Chiu CH, Tang P: Identification of microRNA in the protist Trichomonas vaginalis. Genomics. 2009, 93: 487-493. 10.1016/j.ygeno.2009.01.004.View ArticlePubMedGoogle Scholar
- Wang CR, Xu MJ, Fu JH, Nisbet AJ, Chang QC, Zhou DH, Huang SY, Zou FC, Zhu XQ: Characterization of microRNAs from Orientobilharzia turkestanicum, a neglected blood fluke of human and animal health significance. PLoS ONE. 2012, 7: e47001-10.1371/journal.pone.0047001.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu MJ, Ai L, Fu JH, Nisbet AJ, Liu QY, Chen MX, Zhou DH, Zhu XQ: Comparative characterization of microRNAs from the liver flukes Fasciola gigantica and F. hepatica. PLoS ONE. 2012, 7: e53387-10.1371/journal.pone.0053387.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu MJ, Fu JH, Nisbet AJ, Huang SY, Zhou DH, Lin RQ, Song HQ, Zhu XQ: Comparative profiling of microRNAs in male and female adults of Ascaris suum. Parasitol Res. 2013, 112: 1189-1195. 10.1007/s00436-012-3250-x.View ArticlePubMedGoogle Scholar
- Xu MJ, Liu Q, Nisbet AJ, Cai XQ, Yan C, Lin RQ, Yuan ZG, Song HQ, He XH, Zhu XQ: Identification and characterization of microRNAs in Clonorchis sinensis of human health significance. BMC Genomics. 2010, 11: 521-10.1186/1471-2164-11-521.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu MJ, Wang CR, Huang SY, Fu JH, Zhou DH, Chang QC, Zheng X, Zhu XQ: Identification and characterization of microRNAs in the pancreatic fluke Eurytrema pancreaticum. Parasit Vectors. 2013, 6: 25-10.1186/1756-3305-6-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Hakimi MA, Cannella D: Apicomplexan parasites and subversion of the host cell microRNA pathway. Trends Parasitol. 2011, 27: 481-486. 10.1016/j.pt.2011.07.001.View ArticlePubMedGoogle Scholar
- Zhao GH, Xu MJ, Zhu XQ: Identification and characterization of microRNAs in Baylisascaris schroederi of the giant panda. Parasit Vectors. 2013, 6: 216-10.1186/1756-3305-6-216.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang J, Liu X, Jia B, Lu H, Peng S, Piao X, Hou N, Cai P, Yin J, Jiang N, Chen Q: A comparative study of small RNAs in Toxoplasma gondii of distinct genotypes. Parasit Vectors. 2012, 5: 186-10.1186/1756-3305-5-186.PubMed CentralView ArticlePubMedGoogle Scholar
- TenOever BR: RNA viruses and the host microRNA machinery. Nat Rev Microbiol. 2013, 11: 169-180. 10.1038/nrmicro2971.View ArticlePubMedGoogle Scholar
- Zhai H, Fesler A, Ju J: MicroRNA: a third dimension in autophagy. Cell Cycle. 2013, 12: 246-250. 10.4161/cc.23273.PubMed CentralView ArticlePubMedGoogle Scholar
- Marshall ES, Elshekiha HM, Hakimi MA, Flynn RJ: Toxoplasma gondii peroxiredoxin promotes altered macrophage function, caspase-1-dependent IL-1beta secretion enhances parasite replication. Vet Res. 2011, 42: 80-10.1186/1297-9716-42-80.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu Y, Lan J, Wu X, Yang D, Zhang Z, Nie H, Hou R, Zhang R, Zheng W, Xie Y, Yan N, Yang Z, Wang C, Luo L, Liu L, Gu X, Wang S, Peng X, Yang G: Identification of Dirofilaria immitis miRNA using Illumina deep sequencing. Vet Res. 2013, 44: 3-10.1186/1297-9716-44-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Søkilde R, Vincent M, Møller AK, Hansen A, Høiby PE, Blondal T, Nielsen BS, Daugaard G, Møller S, Litman T: Efficient identification of miRNAs for classification of tumor origin. J Mol Diagn. 2014, 16: 106-115. 10.1016/j.jmoldx.2013.10.001.View ArticlePubMedGoogle Scholar
- Endo K, Weng H, Naito Y, Sasaoka T, Takahashi A, Fukushima Y, Iwai N: Classification of various muscular tissues using miRNA profiling. Biomed Res. 2013, 34: 289-299. 10.2220/biomedres.34.289.View ArticlePubMedGoogle Scholar
- Zhu C, Ren C, Han J, Ding Y, Du J, Dai N, Dai J, Ma H, Hu Z, Shen H, Xu Y, Jin G: A five-microRNA panel in plasma was identified as potential biomarker for early detection of gastric cancer. Br J Cancer. 2014, doi:10.1038/bjc.2014.119Google Scholar
- Xu MJ, Zhou DH, Huang SY, Zhao FR, Nisbet AJ, Lin RQ, Song HQ, Zhu XQ: Comparative characterization of microRNA profiles of different genotypes of Toxoplasma gondii. Parasitology. 2013, 140: 1111-1118. 10.1017/S003118201300053X.View ArticlePubMedGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25: 1966-1967. 10.1093/bioinformatics/btp336.View ArticlePubMedGoogle Scholar
- Mount DW: Using the Basic Local Alignment Search Tool (BLAST). CSH Protoc. 2007, 2007: p17.Google Scholar
- Kruger J, Rehmsmeier M: RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 2006, 34: W451-W454. 10.1093/nar/gkl243.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: e179-10.1093/nar/gni178.PubMed CentralView ArticlePubMedGoogle 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-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Wei Y, Chen S, Yang P, Ma Z, Kang L: Characterization and comparative profiling of the small RNA transcriptomes in two phases of locust. Genome Biol. 2009, 10: R6-10.1186/gb-2009-10-1-r6.PubMed CentralView 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.