the prevalence of p. skrjabini in iran has been reported 4.2% in goat , 5.43% in sheep ( 1 , 2 ) , and 0.8% in wild sheep ( 1 , 3 ) .
researchers could analyze the divergences and genetic distances among taxa in the rdna nucleotide sequences to perform phylogenetic relationships analysis and species identification researches .
the rdna secondary structures are predicted from the corresponding primary sequences according to base pairing , containing all the sequence information .
although , there are significant variations in rdna sequences across different taxa ( 4 ) , the corresponding secondary structures of the transcribed rrna are highly conserved during evolution ( 5 ) , perhaps due to the important role of the rrna folding in holding the structural rna functions ( 6 ) .
in addition , the secondary structures are more conserved than the primary sequences for the semi - compensatory or compensatory mutations , and therefore when the multiple sequence alignments look less reliable due to deletion or insertion , the structures can help to make more reliable assignment of nucleotide homology with important role in the phylogeny ( 7 ) .
in addition , some changes , like expansions and deletions , of a certain helix could be specific to a taxon to help in species identification .
 so , the secondary structures have drawn a lot of attention from phylogenetic scientists 
however , the study on secondary structures takes a slower step than that primary sequences because of the limited sequence data in genbank suitable for structure prediction . in this study , the secondary structures of 5.8srrna of p. skrjabini were comprehensively investigated and compared with secondary structure of a p. skrjabini ( eu375510.1 ) and two species of h. muscae ( ay251024.1 ) and h. microstoma ( ay251023.1 ) sequences , which are located in the same family retrieved from ncbi .
such case studies are suitable to provision of basic data , both for reconstructing molecular evolution in expansive phylogenetic contexts and for analyzing function in ribosome biogenesis .
however , only one study could be derived on parabronema at a molecular level ( 9 ) .
the aim of this study was to analyze the 5.8s rrnagene sequence and to study the possible effect of nucleotide substitutions on the topology of the secondary structure of the 5.8s rrna molecule in p. skrjabini and related habronematidae species .
consequently , we studied the 5.8s rrna gene sequence of the members of the family habronematidae and determined the probable secondary structure of the5.8s rrna molecule for the first time .
dna extraction from worms was performed using an extraction kit ( mbst , iran ) according to the manufacturer s instructions .
the rrna gene was amplified using the primer pairs based on the r dna genome sequence ( 9 ) .
the forward primer was ps - f : 5-gta ggt gaa cct gcg gaa gg -3 and reverse primer was ps - r : 5-ttagtttcttttcct ccgct -3. the pcr reaction was carried out in a total volume of 100l containing 1pcr buffer , 100 mmol mgcl 

 2 

 , 100 m dntp mix ( cinaclone , iran ) , 20 mol of each primer ( cinaclone co. ) , 5 unit/l taq dna polymerase ( cinaclone ) and 1l of template dna ( 100 ng dna ) in an automated thermocycler .
the pcr was performed using the following protocol : 5 min incubation at 94 c , 33 cycles of 45 s at 94 c , 45s at 59 c , and 45 s at 72 c , with an additional extension step for 5 min at 72 c .
the pcr products were analyzed on 1% agarose gels in 0.5 tbe buffer and visualized using sybersafe staining ( cinaclon , iran ) and a uv illuminator .
the pcr product was purified using a quick pcr product purification kit ( mbst , iran ) according to the manufacturer s instructions .
genomic dna sequencing using the sanger method was performed in both directions on the pcr product by the kawsar biotech co. iran .
the sequences were analyzed using the geneious 5.1.6 software and compared against genbank ( www.ncbi.nlm.nih.gov/ ) using the  basic local alignment search tool ( blast ) .
the probable secondary structure of the 5.8s rrna molecule was constructed using mfold software ( http://mfold.rit.albany.edu/ ) .
common structural elements were initially recognized with the help of mfold(10 , 11 ) by screening for thermodynamically optimal and suboptimal secondary structures ( default settings , with t = 37 c ) .
energy levels of the presumptive secondary structures were then calculated with mfold ( 10 , 11 ) .
alignment , analysis of the derived nucleotide sequences , and cluster analysis were performed using mega 6.0 software ( 12 ) .
dna extraction from worms was performed using an extraction kit ( mbst , iran ) according to the manufacturer s instructions .
the rrna gene was amplified using the primer pairs based on the r dna genome sequence ( 9 ) .
the forward primer was ps - f : 5-gta ggt gaa cct gcg gaa gg -3 and reverse primer was ps - r : 5-ttagtttcttttcct ccgct -3. the pcr reaction was carried out in a total volume of 100l containing 1pcr buffer , 100 mmol mgcl 

 2 

 , 100 m dntp mix ( cinaclone , iran ) , 20 mol of each primer ( cinaclone co. ) , 5 unit/l taq dna polymerase ( cinaclone ) and 1l of template dna ( 100 ng dna ) in an automated thermocycler .
the pcr was performed using the following protocol : 5 min incubation at 94 c , 33 cycles of 45 s at 94 c , 45s at 59 c , and 45 s at 72 c , with an additional extension step for 5 min at 72 c .
the pcr products were analyzed on 1% agarose gels in 0.5 tbe buffer and visualized using sybersafe staining ( cinaclon , iran ) and a uv illuminator .
the pcr product was purified using a quick pcr product purification kit ( mbst , iran ) according to the manufacturer s instructions .
genomic dna sequencing using the sanger method was performed in both directions on the pcr product by the kawsar biotech co. iran .
the sequences were analyzed using the geneious 5.1.6 software and compared against genbank ( www.ncbi.nlm.nih.gov/ ) using the  basic local alignment search tool ( blast ) .
the probable secondary structure of the 5.8s rrna molecule was constructed using mfold software ( http://mfold.rit.albany.edu/ ) .
common structural elements were initially recognized with the help of mfold(10 , 11 ) by screening for thermodynamically optimal and suboptimal secondary structures ( default settings , with t = 37 c ) .
energy levels of the presumptive secondary structures were then calculated with mfold ( 10 , 11 ) .
alignment , analysis of the derived nucleotide sequences , and cluster analysis were performed using mega 6.0 software ( 12 ) .
the length of 5.8s rrna sequence for p.skrjabini#1 , p.skrjabini#2 , h. microstoma and h. muscae was 158 , 156 , 127 and 127bp , and the dg required for the formation of the secondary structure was 70.50 , 56.40 , 41.50 and 41.40 kcal / mol , respectively .
the stem loop structures were folded using the mfold web server ( http://mfold.rna.albany.edu/ ) zuker 2003 ( fig .
the 5.8s rrna gene sequences in the habronematidae samples analyzed maximization of the hydrogen bonding formed solid stems , and the largest negative delta g value ( free energy ) .
we predicted the motifs for 5.8s rdna of p. skrjabini and three other sequences in genebank ( p. skrjabini#2 , h. muscae and h. microstoma ) .
three motifs were identified for 5.8s of p. skrjabini#1 , the first motive- aggggg ( 6bp ) , the second motive- taaaaa ( 6bp ) and the third motive- caaaga ( 6bp ) .
the first motive is repeated in two positions [ 10 and 148 ] . in addition , the second motive is repeated in two positions [ 39 and 61 ] . for p. skrjabini#2
three motifs were identified , the first motive  gataaatagtgcgaattgca ( 20bp ) , the second motive- gtggat ( 6bp ) and the third motive- ccatcggg ( 8bp ) .
the second motive is repeated in two positions [ 13 and 29 ] . for h. muscae
there are three motifs , the first motive  catcccgatggt ( 12bp ) , the second motive  gtcgat ( 6bp ) and the third motive  cagacg ( 6bp ) .
the second motive is repeated in two positions [ 3 and 122 ] . for h. microstoma
we identified three motifs , the first motive  cccgatggt ( 9bp ) , the second motive  agctgc ( 6bp ) and the third motive  cagacg ( 6bp ) .
the number of different loops , motifs and the dg for formation of the secondary structure is shown in table 1 , fig .
2 . statistical information of the predicted secondary structure of parabronema skrjabini rrna phylogenetic tree of habronematidae family
the secondary structures contain more information than the primary sequences and are the bases of rrna function ; have gained a lot of attention in phylogenetic analysis . in this study , the secondary structure of 5.8s rrna of p. skrjabini was predicted and the structure comparison performed by predicating the secondary structures for a p. skrjabini and two species of habronema ( h. muscae and h. microstoma ) sequences , which is located in the same family retrieved from ncbi .
5.8s ribosomal rna ( 5.8s rrna ) is a non - coding rna component of the large subunit of the eukaryotic ribosome that plays an important role in protein translation . in this study secondary structure with bulge , hairpins , helices , interior ,
external and multi loops of 5.8s rrna sequence of p. skrjabini was reconstructed under specific settings for base pairing and compared with secondary structure of a sequence of p. skrjabini , and two species of habronema ( h. muscae and h. microstoma ) retrieved from ncbi .
the predicted ssu rrna secondary structure in the present study was the first model for p. skrjabini .
however , the prediction of the secondary structure was impaired due to a few complete sequences of limited species of parabronema . in total ,
41 variable sites in the 5.8s rrna gene sequence were detected in the samples analyze ( 
 fig .
p. skrjabini#2 was more similar to h. muscae and h. microstoma and was more invariable than p. skrjabini#1 .
there was a 31-nucleotide absence at the beginning of the 5.8s rrna gene of h. muscae and h. microstoma .
however , this insertion is present in the p. skrjabini#1and this presence is 30 for p.skrjabini#2 .
the motive 2a and 2b for p. skrjabini#2 and the motive 2a and 2b for h. muscae are different in only one nucleotide ( g - c ) .
the motive 3a in h. muscae and 3a in h. microstoma is common and motive 1a in h. muscae and h. microstoma had only three - nucleotide difference .
we isolated p. skrjabini#1 in sheep and from sanandaj ( west of iran ) , but p. skrjabini#2 is in camel and china , h. muscae and h. microstoama are in horse and italy . the phylogenetic tree ( fig .
d : h. muscae ( ay251024.1).the position of motifs is shown on the secondary structures .
the nucleotide substitutions are shown by arrows moreover , some changes , such as expansions and deletions , of a certain helix could be specific to a taxon to help a lot in species identification . using mfold software , the probable secondary structure of the 5.8s rrna gene
when constructing the secondary structure of the 5.8s rrna gene sequence , the sequence of p. skrjabini was used as a consensus sequence .
the secondary structure of the 5.8s rrna gene has not been described for family habronematidae , specially the genus of p. skrjabini . in our study , for the first time we determined the secondary structure of 5.8s rrna in the p. skrjabini and related habronematidae species .
we made a multiple alignment and found structural differences among the analyzed samples , parabronema and habronema , which could be further used in the structure modeling across habronematidae .
 background : genomic dna was isolated from parabronema skrjabini . 
 rrna region was amplified and sequenced.methods:the rna secondary structure was predicted using mfold software ( 
 http://mfold.rit.albany.edu 
 ) . 
 the secondary structure with bulge , hairpins , helices , interior , external and multi loops was predicted for 5.8srdna of our sequence of p. skrjabini and a sequence of p. skrjabini and two species of habronema ( h. microstoma and h. muscae ) in genbank . 
 rna motifs were predicted by meme program version 4.10.2.results:the length of 5.8s rrna sequence for p. skrjabini#1 , p. skrjabini#2 , h. microstoma and h. muscae was 158 , 156 , 127 and 127bp , and the dg required for the formation of the secondary structure was 70.50 , 56.40 , 41.50 and 41.40 kcal / mol , respectively . 
 common structural elements were initially recognized with the help of mfold by screening for thermodynamically optimal and suboptimal secondary structures ( default settings , with t = 37 c ) . 
 the energy levels of the presumptive secondary structures were then calculated with mfold at the dna level . both motifs and 
 the sequence of p. skrjabini#1 were completely different from the other analyzed samples . 
 this difference might be due to the differences in host and geographical area.conclusion:this is the first molecular study of p. skrjabini in sheep , which could be further used in the structure modeling across habronematidae . 
Introduction
Material and Methods
DNA Extraction and PCR
Isolation of DNA, PCR and Sequencing
Results
Discussion
Conclusion