|
|
||||||||
ENVIRONMENT, WELL-BEING, AND BEHAVIOR |

* Department of Biotechnology, Stephen F. Austin State University, P.O. Box 6132, Nacogdoches, TX 75965; and
USDA, Agricultural Research Service, Southern Plains Agricultural Research Center, Food and Feed Safety Research Unit, College Station, TX, 77845
1 Corresponding author: avankley{at}sfasu.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: denaturing gradient gel electrophoresis Eimeria coccidiosis species identification
| INTRODUCTION |
|---|
|
|
|---|
| MATERIALS AND METHODS |
|---|
|
|
|---|
Individual droppings from ten 30-d-old chickens from 10 poultry houses in east Texas were collected and stored at 4°C to prevent DNA degradation. The samples from each poultry house were pooled and 1 g of pooled sample was washed and vortexed with glass beads (710 to 1,180 µm, Sigma, Oakville, Ontario, Canada) for 15 min twice in 5 mL of water and then filtered through cheese cloth. After centrifugation, the samples were resuspended in stool lysis buffer. To extract the DNA we used the QIAamp DNA Stool Mini Kit (Qiagen Inc., Mississauga, Ontario, Canada) following the instructions of the manufacturer.
Primer Design
To design Eimeria genus-specific primers to amplify a fragment suitable for DGGE, 27 sequences of the nearly complete SSU rRNA gene were obtained from GenBank (Benson et al., 2006). The GenBank accession numbers and the species names of the sequences used in this study were as follows: E. acervulina (DQ136187.1 and DQ538351.1), Eimeria bovis (U77084.1 [GenBank] ), E. brunetti (U67116.1 [GenBank] ), E. maxima (DQ136186.1, DQ538350.1, EF210322.1, U67117.1 [GenBank] , DQ538348.1, DQ640012.1, DQ538349.1, and EF122251.2), E. mitis (U40262.1 [GenBank] ), E. mivati (U76748.1 [GenBank] ), E. necatrix (DQ136185.1 and U67119.1 [GenBank] ), E. praecox (U67120.1 [GenBank] ), and E. tenella (EF210325.1, DQ136181.1, AF026388.1, DQ640011.1, U40264.1 [GenBank] , DQ136183.1, DQ136184.1, DQ136180.1, DQ136179.1, and EF210326.1). Sequences were aligned with ClustalW (Thompson et al., 1994) and used to calculate the theoretical melting profiles with MELT94 software (Michikawa et al., 1997). The melting profiles and the multiple sequence alignment were overlaid to select a fragment suitable for DGGE analysis. The fragment was required to have a single melting domain and to be flanked by conserved regions to accommodate annealing of genus-specific DGGE primers. In addition, multiple sequence alignment was used to design another pair of Eimeria genus-specific SSU primers to amplify the nearly complete SSU rRNA gene for phylogenetic analysis.
PCR and Cloning of the Nearly Complete SSU rRNA Gene
To obtain PCR products suitable for cloning and sequencing, we carried out PCR with cloned Pfu DNA Polymerase AD (Stratagene Inc., La Jolla, CA) and Eimeria genus-specific SSU primers: forward primer, 5'-TTG TCT CAA AGA TTA AGC C-3'; reverse primer, 5'-AGC GAC GGG CGG TGT GTA CAA-3'. The extracted DNA (3 to 30 ng) with the SSU primer pair (0.5 µM) was used to amplify a long fragment (~1,600 bp) in a final volume of 25 µL. The PCR program included an initial denaturation cycle of 95°C for 3 min followed by 25 cycles of denaturation at 94°C for 35 s, annealing at 56°C for 25 s, and an extension at 72°C for 2 min. The final extension step was increased to 4 min at 72°C. The PCR product was cloned with a pPCR Script-Amp cloning kit (Stratagene Inc.) into a pPCR Script-Amp SK(+) vector according to the manufacturers protocol (Stratagene Inc.). The vector was transformed into Escherichia coli XL10 competent cells (Stratagene Inc.). Selection of transformants was done by blue-white screening. Colony PCR was performed for 350 clones to amplify the SSU rDNA insert. The screening yielded 58 positive transformants. Thirty randomly chosen clones were sent for sequencing. Two clones that failed to sequence and 5 clones that were only partially sequenced were not included in further analysis.
Sequence Analysis
Both strands of each of the clones were sequenced by Amplicon Express (Pullman, WA). The chromatograms obtained were assembled with the STADEN package (Staden et al., 2000). Sequences were checked for contamination using the VecScreen (http://www.ncbi.nlm.nih.gov/VecScreen/) and segments of vector origin were removed manually. The sequences were deposited in GenBank under accession numbers EU044765 [GenBank] –85. A search by BLAST (Altschul et al., 1990) was conducted against the GenBank database and the top 5 matching sequences for every clone were retrieved. The sequences were combined and after excluding redundant entries, aligned by ClustalW 1.81 (Thompson et al., 1994). Multiple alignment parameters such as gap opening and extension penalties available in ClustalW were modified to improve the quality of the alignment. The final alignment was manually examined for obvious errors. The alignment and additional information is available at http://biotech.sfasu.edu/bt/EimeriaID/. Phylogenetic trees were constructed by using maximum likelihood (ML) and Bayesian inference (BI) methods. The ML method was carried out in PHYML (Guindon and Gascuel, 2003) and Bayesian inference in MrBayes 3.0 (Huelsenbeck and Ronquist, 2001). For both analyses we used a general-time-reversible model for nucleotide substitution allowing both programs to estimate the proportion of invariable sites and the gamma distribution parameter. The BIONJ option was selected to generate an initial tree for the ML method and the BI was started with a random tree. The reliability of the ML tree was estimated by the approximate likelihood ratio test method (Anisimova and Gascuel, 2006). For the BI analysis, 4 Markov chains were run for 5,000,000 generations and sampled every 1,000 generations. After elimination of the first 25% of the trees as burn-in, 3,750 trees were used to generate a 50% majority-rule consensus-tree.
DGGE Identification Ladder Preparation
The 23 successfully sequenced plasmids were used to carry out PCR (final volume of 50 µL). One microliter of each purified recombinant plasmid was used as template to amplify an ~120-bp fragment with the genus-specific DGGE primers (0.5 µM) that were previously designed: forward primer containing GC-clamp 5'-GCC CGC CGC GCC CGC GCC CGT CCC GCC GCC CCC GCC CGG ATT AGA TAC AAA ACC AAC CC-3', and reverse primer 5'-GCT GAT AGG TCA GAA ACT TG-3'. The amplification process was performed using 25 µL of JumpStart REDTaq ReadyMix PCR Reaction Mix (Sigma Chemical Co., St. Louis, MO) in a Mastercycler (Eppendorf Scientific Inc., Westbury, NY) according to the following program: initial denaturation at 94°C for 3 min followed by 35 cycles of denaturation at 92°C for 30 s, annealing at 60°C for 35 s, and an extension at 72°C for 2 min. A final extension was performed at 72°C for 7 min. Three microliters of each PCR product was tested on a 2% agarose gel. The remaining 42 µL was mixed to obtain the DGGE identification ladder.
DGGE Analysis of Field Samples
Five microliters of DNA extracted from field samples was used to perform PCR with the pair of DGGE primers using the same PCR conditions described above. The electrophoresis was performed both with the PCR products (20 µL) and the DGGE identification ladder (20 µL) for 4.5 h at 250 V in 1 x Tris-acetate-EDTA (TAE) buffer at 60°C using a DGGE-2001 apparatus (CBS Scientific Co., Del Mar, CA). Polyacrylamide gels (8%, 0.75-mm thick) were prepared with a denaturing gradient of urea-formamide mix (Myers et al., 1987) ranging from 35 to 50%. After electrophoresis, the gels were stained with SYBR Green I (Invitrogen Inc., Carlsbad, CA) for 30 min and subsequently photographed using a BioRad Imager System equipped with a Gel Doc XR camera and Quantity-One software (Bio-Rad Inc., Hercules, CA).
| RESULTS |
|---|
|
|
|---|
To identify the species of Eimeria present in field samples we cloned and sequenced 22 nearly complete SSU rRNA genes. The BLAST queries against Gen-Bank identified several related sequences of Eimeria. However, we could not identify all sequences to species level. For example, a query of clone 93 returned an E. brunetti sequence with 98% identity (1,565/1,587) to our sequence as the first match and E. acervulina with 98% (1,559/1,591) as the second match. We attempted to clarify relationships among the sequences by using a character-based phylogenetic approach in which phylogenetic trees were inferred by ML and BI methods. All trees were rooted with the outgroup species of E. bovis (Figures 1
and 2
). As expected, the avian Eimeria species formed a monophyletic clade in both trees and addition of our sequences did not substantively change the previously published tree topology (Barta et al., 1997). In the first clade (Figures 1
and 2
), marked by Roman numeral I, there was a high likelihood for a close relationship between E. necatrix and E. tenella. Both analyses grouped E. tenella sequences with the sequence of E. necatrix; however, the branching patterns within the clade were less resolved. In the ML tree E. necatrix appeared to be more closely related to E. tenella (EF210325.1) by excluding the E. tenella (EF210326.1) sequence. The BI analysis did not resolve the branching pattern in this clade and collapsed to polytomy. Nonetheless, both phylogenetic trees showed affiliation of clone 69 with E. necatrix and clones 141, 144, and 206 with E. tenella (EF210326.1). In the second clade (II), E. maxima was strongly affiliated with E. praecox and E. mitis with E. mivati. However, the position of E. brunetti (U67116.1
[GenBank]
) according to ML was different from that resulting from BI analyses. The BI analysis supported identification of clones 112, 143, 165, and 325 by grouping them with the E. brunetti sequence, whereas other clones collapsed in an unresolved polytomy. The ML tree provided a different affiliation for E. brunetti; it was grouped with clones 30, 93, and 126, and clones 131 and 251 were associated with E. acervulina.
|
|
In conclusion, from sequences of 22 clones, 2 sequences were identified as chimeric, 11 were unambiguously identified to species, and the remaining 9 were identified differently by ML and BI methods. The summary of the results obtained from phylogenetic analysis is shown in Table 1
.
|
Based on multiple sequence alignment of 8 avian Eimeria SSU rDNA sequences from GenBank, a genus-specific primer set was designed to amplify a 114- to 118-bp fragment. The forward primer had a single nucleotide mismatch in its priming site with the sequences from E. mitis and E. maxima. In contrast, the reverse primer was completely homologous to all reference sequences. The fragment contained significant differences in the sequences of the species of avian Eimeria indicating the feasibility of using this fragment for species identification. The PCR amplification carried out with the set of primers was successful in all tested samples, which included cloned SSU rDNA as well as DNA extracted from fecal field samples. The expected size of the PCR-amplified fragments was verified on agarose gel (data not shown). To test the specificity of the primers we carried out a PCR with template DNA extracted from fecal samples obtained from coccidia-free 2-d-old chicks: no PCR product was obtained (data not shown). The PCR products obtained from amplification of clones were subjected to DGGE analysis. An example of a DGGE profile is shown in Figure 3
. For 22 sequenced clones we obtained 9 sequence-specific migrations on the 35 to 50% gradient gel. Clone 334 was disregarded, because both BI and ML analyses failed to identify its phylogenetic position. Clone 30 was also disregarded, because its position in the BI tree was unresolved and it was the only sequence that exhibited band position A. Sequences that could not be identified or were identified differently by BI and ML methods were disregarded for identification of ladder bands. The PCR products from other clones were mixed to obtain an identification ladder with 7 bands (Figure 4
; lanes 1 and 12) and were identified as follows: band B = E. tenella; C = E. necatrix; D = E. tenella; E = E. maxima/E. praecox/E. brunetti; F = E. maxima/E. praecox/E. brunetti; G = E. tenella/E. acervulina; H = E. acervulina.
|
|
| DISCUSSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
Received for publication March 2, 2008. Accepted for publication April 22, 2008.
| REFERENCES |
|---|
|
|
|---|
Anisimova, M., and O. Gascuel. 2006. Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst. Biol. 55:539–552.[CrossRef][Web of Science][Medline]
Barta, J. R., D. S. Martin, P. A. Liberator, M. Dashkevicz, J. W. Anderson, S. D. Feighner, A. Elbrecht, A. Perkins-Barrow, M. C. Jenkins, H. D. Danforth, M. D. Ruff, and H. Profous-Juchelka. 1997. Phylogenetic relationships among eight Eimeria species infecting domestic fowl inferred using complete small subunit ribosomal DNA sequences. J. Parasitol. 83:262–271.[CrossRef][Medline]
Benson, D. A., I. Karsch-Mizrachi, D. J. Lipman, J. Ostell, and D. L. Wheeler. 2006. GenBank. Nucleic Acids Res. 34:D16–D20.
Blake, D. P., A. L. Smith, and M. W. Shirley. 2003. Amplified fragment length polymorphism analyses of Eimeria spp.: An improved process for genetic studies of recombinant parasites. Parasitol. Res. 90:473–475.[CrossRef][Web of Science][Medline]
Cocolin, L., K. Rantsiou, L. Iacumin, C. Cantoni, and G. Comi. 2002. Direct identification in food samples of Listeria spp. and Listeria monocytogenes by molecular methods. Appl. Environ. Microbiol. 68:6273–6282.
Fernandez, S., A. H. Pagotto, M. M. Furtado, A. M. Katsuyama, A. M. Madeira, and A. Gruber. 2003. A multiplex PCR assay for the simultaneous detection and discrimination of the seven Eimeria species that infect domestic fowl. Parasitology 127:317–325.[CrossRef][Web of Science][Medline]
Gasser, R. B., R. Skinner, R. Fadavi, G. Richards, and G. Morris. 2005. High-throughput capillary electrophoresis for the identification and differentiation of seven species of Eimeria from chickens. Electrophoresis 26:3479–3485.[Medline]
Guindon, S., and O. Gascuel. 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52:696–704.
Huelsenbeck, J. P., and F. Ronquist. 2001. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17:754–755.
Lien, Y. Y., S. C. Sheu, H. J. Liu, S. C. Chen, M. Y. Tsai, S. C. Luo, K. C. Wu, S. S. Liu, and H. Y. Su. 2007. Cloning and nucleotide sequencing of the second internal transcribed spacer of ribosomal DNA for three species of Eimeria from chickens in Taiwan. Vet. J. 173:184–189.[CrossRef][Web of Science][Medline]
Long, P. L., and L. P. Joyner. 1984. Problems in the identification of species of Eimeria. J. Protozool. 31:535–541.[Medline]
Long, P. L., B. J. Millard, L. P. Joyner, and C. C. Norton. 1976. A guide to laboratory techniques used in the study and diagnosis of avian coccidiosis. Folia Vet. Lat. 6:201–217.[Medline]
Mathieu-Daude, F., J. Welsh, T. Vogt, and M. McClelland.1996. DNA rehybridization during PCR: The Cot effect and its consequences. Nucleic Acids Res. 24:2080–2086.
Michikawa, Y., G. Hofhaus, L. S. Lerman, and G. Attardi. 1997. Comprehensive, rapid and sensitive detection of sequence variants of human mitochondrial tRNA genes. Nucleic Acids Res. 25:2455–2463.
Myers, R. M., T. Maniatis, and L. S. Lerman. 1987. Detection and localization of single base changes by denaturing gradient gel electrophoresis. Methods Enzymol. 155:501–527.[Web of Science][Medline]
Satoh, M., and Y. Nakai. 2007. Discrimination of Cryptosporidium species by denaturing gradient gel electrophoresis. Parasitol. Res. 101:463–466.[CrossRef][Web of Science][Medline]
Shirley, M. W. 1975. Enzyme variation in Eimeria species of the chicken. Parasitology 71:369–376.[Web of Science][Medline]
Shirley, M. W., and N. Bumstead. 1994. Intra-specific variation within Eimeria tenella detected by the random amplification of polymorphic DNA. Parasitol. Res. 80:346–351.[CrossRef][Web of Science][Medline]
Staden, R., K. F. Beal, and J. K. Bonfield. 2000. The Staden package, 1998. Methods Mol. Biol. 132:115–130.[Medline]
Su, Y. C., A. C. Fei, and F. M. Tsai. 2003. Differential diagnosis of five avian Eimeria species by polymerase chain reaction using primers derived from the internal transcribed spacer 1 (ITS-1) sequence. Vet. Parasitol. 117:221–227.[CrossRef][Web of Science][Medline]
Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:4673–4680.
Walter, J., G. W. Tannock, A. Tilsala-Timisjarvi, S. Rodtong, D. M. Loach, K. Munro, and T. Alatossava. 2000. Detection and identification of gastrointestinal Lactobacillus species by using denaturing gradient gel electrophoresis and species-specific PCR primers. Appl. Environ. Microbiol. 66:297–303.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |