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
Skip to main content
© 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.
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 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.
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.
Common and taxon-specific reads of Ascaris lumbricoides and Ascaris suum
Unique sRNAs (%)
Total sRNAs (%)
A. lumbricoides specificb
A. suum specificc
Comparison of miRNA profiles in Ascaris lumbricoides (Alu) and A. suum (Asu)
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.
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.
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.