1. Introduction {#sec0005}
===============

Apicomplexan parasites can cause serious human and animal disease (e.g. [@bib0105], [@bib0195], [@bib0170]). Control strategies vary between genera but effective or improved vaccines are a major goal for all. Experimental recombinant or subunit vaccines have been described for many apicomplexans with varying levels of efficacy obtained under controlled laboratory or experimental animal conditions. However, commercial translation has commonly been hindered by a range of factors including naturally occurring genetic diversity in field parasites, resulting in insufficient immunological protection in the vaccinated host and selection of resistant populations ([@bib0085], [@bib0040]). If vaccines are to be successful and remain effective in the long term it is essential to understand the impacts that parasite genetic (antigenic) diversity and population structures have on the selection of field populations capable of vaccine escape. Studies into naturally occurring diversity have yielded notable insights for parasites such as *Plasmodium falciparum*, revealing varied levels of polymorphism with evidence of clonal, as well as panmictic structures, depending on geographic location and rates of transmission ([@bib0010], [@bib0130]). Locus-specific analyses have also detected considerable diversity for *P. falciparum* vaccine candidates such as apical membrane antigen-1 and merozoite surface protein-1 (AMA1 and MSP1, ([@bib0115], [@bib0200])). In contrast, *Toxoplasma gondii* has been defined as having a clonal population structure across much of the world apart from South America, where genetic diversity and inter-clonal interbreeding is common ([@bib0155]). For most other apicomplexans details on diversity are scarce; for example for *Eimeria*, the cause of the enteric disease coccidiosis, rather little is known at the molecular level. Recently, populations of *Eimeria tenella* parasites isolated from chicken facilities in Egypt, Libya, India and Nigeria were genotyped using a Sequenom MassARRAY single nucleotide polymorphism (SNP) tool, revealing notable variation in haplotype diversity and population structure with a North/South regional divide ([@bib0040]). *E. tenella* populations sampled below the 30th parallel north (30°N latitude) were defined by extensive haplotype diversity with no significant signatures of selection. In contrast, *E. tenella* populations sampled north of this latitude presented limited haplotype diversity with significant disequilibrium ([@bib0040]). More recent analysis of internal transcribed spacer (ITS) sequence diversity for *E. tenella* has confirmed this apparent variation in population structure, with notably limited diversity detected in the United States of America (USA) ([@bib0060]). While these data are informative, expanding research to *E. tenella* populations in other regions and laboratories requires a more accessible technique.

Genetic approaches for differentiating between parasite isolates have improved rapidly over the last 20 years. Reviewed previously ([@bib0025]), examples including random amplification of polymorphic DNA (RAPD), single and multi-locus sequence typing, amplified fragment length polymorphism (AFLP) and variable number tandem repeats (VNTR) analyses have all been employed with success for *Eimeria* ([@bib0095], [@bib0185], [@bib0035], [@bib0160], [@bib0140]). Another widely used genotyping technology is polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). PCR-RFLP combines the sensitivity and specificity of PCR with genetic discrimination arising from substitutions, insertions and/or deletions located within recognition sites for restriction endonucleases ([@bib0165]). The approach is straightforward and applicable in many laboratories given provision of a thermal cycler and equipment for gel electrophoresis. Among apicomplexan parasites PCR-RFLP has been employed most notably to genotype *T. gondii*, including examples of parasites derived from multiple host species ([@bib0205], [@bib0080], [@bib0045]) and regions ([@bib0070], [@bib0075], [@bib0210]), and underpinning some of the most fundamental population genetics studies ([@bib0205]). The availability of more accessible and cost-effective population genetics tools will provide opportunities for complimentary, and importantly comparable, analyses of *E. tenella* from other populations.

2. Materials and methods {#sec0010}
========================

2.1. Reference parasites---production, processing and genomic DNA extraction {#sec0015}
----------------------------------------------------------------------------

The *E. tenella* Houghton, Nippon-2, Weybridge and Wisconsin reference isolates ([@bib0190], [@bib0175]) and a Nigerian field isolate ([@bib0120]) were maintained, purified, and propagated as described previously using specific pathogen-free White Leghorn chickens ([@bib0145], [@bib0150]). Genomic DNA was isolated and purified using a QIAamp DNA Stool mini kit as recommended by the manufacturer (Qiagen, Hilden, Germany), including an additional initial mechanical disruption step ([@bib0125]).

2.2. Field sample collection---processing and genomic DNA extraction {#sec0020}
--------------------------------------------------------------------

Twenty seven field samples collected from commercial broiler farms in South/South-East England and the Republic of Ireland between 2013 and 2015 and found by previous diagnostic PCR analysis to include oocysts of *E. tenella* were available for use in this study (n = 24 and 3 respectively, diagnostic analysis not shown). Each sample was collected as described previously ([@bib0125]) and consisted of pooled faecal samples collected from the floor of a poultry house, representing multiple individuals within a single flock. All farms sampled used commercial indoor broiler systems, including either Ross 308 or Cobb500 stock between four and six weeks of age (Supplementary file 1). Routine anticoccidial control was based on chemoprophylaxis on all farms sampled, with no recorded use of anticoccidial vaccination. Total oocysts (including *E. tenella*, plus any other *Eimeria* species which may have been present) were purified by saturated salt flotation, allowed to sporulate and then stored in 2.5% (w/v) potassium dichromate at +4°C. Total genomic DNA was recovered from each sample using a QIAamp mini stool DNA kit (Qiagen) as described previously ([@bib0125]).

2.3. Sequenom marker development for PCR-RFLP {#sec0025}
---------------------------------------------

Each of the 52 *E. tenella* species-specific SNPs included in a Sequenom genotyping array ([@bib0040]) and found previously to yield robust profiles which supported discrimination between field isolates were assessed for development as PCR-RFLPs. The online tool NEBcutter (version 2.0, New England Biolabs; ([@bib0215])) was used to identify restriction endonuclease recognition sites within genomic DNA sequences representing 15 bp centred on each candidate SNP. All SNP-types were tested at each SNP locus. Restriction endonucleases whose recognition site was disrupted by one, but not both SNP-types were recorded for each locus, generating a series of cut/non-cut genetic markers. Where more than one restriction endonuclease was available selection was made based upon (i) activity at 37 °C, (ii) use of the generic Cut Smart buffer (New England Biolabs) and (iii) application in one or more other PCR-RFLPs to minimise complexity of the protocol. Subsequently, 1200 bp genomic sequence centred on each candidate SNP was used as template for primer design with Primer3 ([@bib0180]). Target amplicon size was set to 400--1000 bp with an optimal predicted melting temperature (Tm) centred on 60 °C. A uniform Tm was required in order to simplify application of the eventual PCR-RFLP panel. All primers were synthesised by Sigma-Aldrich (UK).

2.4. Polymerase chain reaction {#sec0030}
------------------------------

Standard PCR amplification was completed using MyTaq DNA polymerase (Bioline). Briefly, each PCR reaction contained 1 μl template DNA (commonly containing between 10 and 50 ng), 1 μl of each of the relevant forward and reverse primers (10 μM stock) and 12.5 μl of MyTaq ×2 mastermix, made up to a final volume of 25 μl using molecular grade water (Sigma). Negative controls included molecular grade water as template. Positive controls included reference strain genomic DNA representative of each SNP type. Thermal cycler parameters were: 1× initial denaturation at 95 °C for 1 min, followed by 35× (denaturation 15 s at 95 °C, annealing 15 s at 56 °C, extension 1 min at 72 °C), followed by a final extension phase of 72 °C for 7 min. All primers used were synthesized by Sigma-Genosys and are shown in [Table 1](#tbl0005){ref-type="table"}. PCR amplicons were sequenced on both strands using the same primers employed in their original amplification (GATC Biotech). Sequence assembly, annotation, and interrogation were undertaken with CLC Main Workbench v6.0.2 (CLC Bio; ([@bib0050])) using BLASTn against the reference *E. tenella* Houghton genome sequence assembly to confirm identity.

2.5. Restriction fragment length polymorphism (RFLP) {#sec0035}
----------------------------------------------------

After thermal cycling each PCR reaction underwent restriction endonuclease digestion. Briefly, 10 μl PCR product was combined with 5 U restriction enzyme ([Table 1](#tbl0005){ref-type="table"}), 3 μl 10× Cut Smart buffer (New England Biolabs) and molecular biology grade water to a final volume of 30 μl, prior to incubation for 60 min at the temperature required for each enzyme (as shown in [Table 1](#tbl0005){ref-type="table"}). Restriction endonucleases were purchased from New England Biolabs, using high fidelity (-HF) versions where necessary to permit use of the same Cut Smart^®^ buffer in all reactions. Subsequently, each digested PCR reaction was resolved by agarose gel electrophoresis using a 2% (w/v) UltraPure agarose gel in 1× Tris-borate-EDTA buffer (TBE; all Sigma-Aldrich), including 0.01% (v/v) SafeView nucleic acid stain (NBS Biologicals). The results of electrophoresis were visualised using a U:Genius Gel Documentation System (Syngene).

2.6. Genetic analysis {#sec0040}
---------------------

Following PCR-RFLP a SNP type was assigned to all eleven loci for every sample by the presence of one or two bands following electrophoresis. The presence of a single band indicated occurrence of the SNP associated with restriction endonuclease recognition site disruption. The presence of two bands indicated occurrence of the SNP associated with maintenance of the site ([Table 1](#tbl0005){ref-type="table"}). Detection of both band profiles in a single sample indicated a polyclonal parasite population with genetic diversity at the target SNP (occurring between different oocysts, or the product of cross-fertilisation in hybrid oocysts). Where polyclonal patterns were detected the dominant (brightest) pattern was retained for genetic analysis. Band profiles were counted and compiled into a Microsoft Excel spreadsheet ([@bib0090]). Base specificity at each SNP was identified and combined for each sample to create a SNP haplotype. Haplotypes derived here from the UK and Irish field samples were combined with the equivalent SNPs extracted from the larger Sequenom SNP panel produced prior to this study (from 244 samples collected in Egypt, India, Libya and Nigeria ([@bib0040])) and converted to Phylip format using CLC Main Workbench v6.0.2.

The program DnaSP (DNA Sequence Polymorphism version 5.10.01 ([@bib0135])) was used to calculate the number of SNP haplotypes in the full and regional haplotype data sets. The standardised index of association (*I*~A~^S^) and associated statistical significance were calculated using LIAN v3.6 ([@bib0110]) as a measure of linkage disequilibrium, where values of 0 indicated no association and values significantly higher than 0 indicated association. The Monte Carlo test option was chosen with 10,000 iterations using p \< 0.05 as a measure of significance. Mean genetic diversity (Het) was also calculated using LIAN. A median joining phylogenetic network was produced for the PCR-RFLP SNP haplotypes using NETWORK version 4.5.1.1 ([@bib0020]).

2.7. Ethics approval {#sec0045}
--------------------

This study was carried out in strict accordance with the Animals (Scientific Procedures) Act 1986, an Act of Parliament of the United Kingdom. All protocols were approved by the Royal Veterinary College Ethical Review Committee. All animal studies were approved by the United Kingdom Government Home Office.

3. Results {#sec0050}
==========

3.1. Identification and validation of PCR-RFLP markers for *E. tenella* {#sec0055}
-----------------------------------------------------------------------

Fifty-two SNPs developed previously for use in a Sequenom MassARRAY SNP typing assay for *E. tenella* were assessed here for development as PCR-RFLPs. Eleven SNPs were prioritised following *in silico* identification of an appropriate restriction endonuclease and accompanying locus specific primer pair ([Table 1](#tbl0005){ref-type="table"}). Direct sequencing of each PCR amplicon derived using genomic DNA extracted from the *E. tenella* Houghton, Nippon-2, Weybridge or Wisconsin reference strains, and a Nigerian field isolate, indicated amplification of a single, correct target for each PCR-RFLP when compared with the reference Houghton strain genome sequence assembly and confirmed SNP identity. Restriction endonuclease digestion confirmed the uncut/cut electrophoresis profiles for each PCR-RFLP ([Fig. 1](#fig0005){ref-type="fig"}, [Table 2](#tbl0010){ref-type="table"}).

3.2. PCR-RFLP application to UK and Irish field isolates {#sec0060}
--------------------------------------------------------

Twenty seven samples collected from UK and Irish broiler farms (n = 24 and 3 respectively) and found previously by PCR to contain *E. tenella* were genotyped using the panel of 11 PCR-RFLPs. All samples were successfully genotyped at all 11 loci with no evidence of polyclonal populations (Supplementary file 1). Analysis of the PCR-RFLP profiles revealed four distinct haplotypes in the UK and Irish sample set, three of which were shared between the UK and Ireland ([Fig. 2](#fig0010){ref-type="fig"}). Mean genetic diversity was low, in agreement with the limited number of haplotypes detected, although the *I*~A~^S^ provided significant evidence of disequilibrium ([Table 3](#tbl0015){ref-type="table"}).

3.3. Comparison of genetic analyses using PCR-RFLP or sequenom SNP markers {#sec0065}
--------------------------------------------------------------------------

Previous analysis of *E. tenella* genetic diversity has been based on a panel of 52 SNP markers, assessed using Sequenom MassARRAY genotyping ([@bib0040]). In recognition of the difficulty many laboratories may experience employing the same technology we have compared results achieved previously using the full Sequenom panel with the subset of SNPs included in the 11 PCR-RFLPs from the same samples. The variation in *E. tenella* haplotype diversity and population structure described previously for North Africa and northern India, compared to Nigeria and southern India remained distinct, indicated by the total number of haplotypes, mean genetic diversity and *I*~A~^S^ ([Table 3](#tbl0015){ref-type="table"}). Populations exhibiting statistically significant disequilibrium were confirmed, despite an ∼five-fold reduction in markers. However, the regional specificity was less obvious, with greater haplotype overlap between regions. NETWORK analysis of the PCR-RFLP SNP subset of Sequenom markers illustrated these findings (compare [Fig. 2](#fig0010){ref-type="fig"}A and B). The previous full Sequenom panel revealed non-overlapping *E. tenella* populations defined by the frequent occurrence of a small number of haplotypes in North Africa and northern India, complemented by many diverse, often unique haplotypes from Nigeria and southern India ([@bib0040] and [Fig. 2](#fig0010){ref-type="fig"}A). Analysis of the PCR RFLP SNPs presented a similar overview of haplotype occurrence but the NETWORK remained unresolved, presenting multiple high-dimensional cubes and cycles ([Fig. 2](#fig0010){ref-type="fig"}B) which persisted even when using Maximum Parsimony. Only removal of the unique haplotypes permitted successful NETWORK construction at the expense of a complete portrayal of diversity. *E. tenella* PCR-RFLP SNP profiles from the UK and Ireland were notably distinct from those described in Africa and India ([Fig. 2](#fig0010){ref-type="fig"}B).

3.4. Identification of optimal PCR-RFLP assays by region {#sec0070}
--------------------------------------------------------

Comparison of PCR-RFLP SNP profiles between regions revealed considerable variation in discriminatory power. Every PCR-RFLP was found to be polymorphic in at least one region, although no single SNP was able to differentiate between genotypes in every region sampled ([Table 4](#tbl0020){ref-type="table"}). Only two (M3208 and C2F235220) provided discrimination between genotypes in the UK and Ireland. Markers considered most discriminatory within each region (assessed as the closest to 50% representation of each SNP-type) were less discriminatory elsewhere.

4. Discussion {#sec0075}
=============

Understanding variation in population structure and genetic diversity for pathogens such as *Eimeria* is important. As we seek to combat drug resistance and develop new, cost-effective anticoccidial vaccines for use with chickens, defining the epidemiology of parasite occurrence and disease in different regions and under different production systems has become essential. The apparently obligate requirement for sexual reproduction in eimerian life cycles ([@bib0220]), together with the vast number of immunologically naïve definitive hosts (chickens) produced every year, offers the potential for co-infection, resegregation/recombination and the generation of diversity among *Eimeria* strains, although the extent to which this occurs remains largely unknown. The development of a PCR-RFLP panel for the related coccidian *T. gondii* provides an excellent example of how having accessible population genetics tools provides opportunities for comparable analyses of parasites from around the world. There are more than 170 peer reviewed manuscripts available in PubMed under the keywords '*Toxoplasma'* and 'PCR-RFLP' (accessed 12th July; 2016); including surveillance of multiple host species in many different regions (e.g. [@bib0205], [@bib0070], [@bib0075], [@bib0080], [@bib0210], [@bib0045]). Development of PCR-RFLP for *Eimeria*; starting with *E. tenella*; will allow expansion of genetics studies with these parasites. Sequencing the PCR-RFLP loci from a panel of reference isolates confirmed the presence of the target SNPs; and provides a robust platform for application of the assay to field isolates.

Recent comparison of *E. tenella* populations in Africa and India found that samples collected in the more northern regions (Egypt, Libya and northern India) displayed less variation in genome-wide genetic diversity and population structure than those from southern regions (Nigeria and southern India) ([@bib0040]). Extension of this survey to include *E. tenella* populations collected from broilers reared in the UK and Ireland now reveals restricted genetic diversity with significant evidence of linkage disequilibrium comparable to the more northern regions sampled previously. It is important to note that just three samples from Ireland were characterised in this study, indicating that additional uncaptured diversity may occur. However, all three Irish haplotypes were also detected in the UK suggesting that *E. tenella* from Ireland and the UK represent one genetic population. Further sampling will be required to confirm this view. The absence of polyclonal *E. tenella* populations in this study may be a consequence of low levels of genetic diversity within the sampled regions, although this can also not be confirmed here. The reasons for reduced diversity among northern *E. tenella* such as those sampled in the UK and Ireland, compared with more southern populations, remain unclear. It is possible that regional variations may exist in biosecurity and the application of anticoccidial control measures, although direct evidence of this is limited at present. Chemoprophylaxis remains the dominant form of control in all regions tested with a variety of ionophore and chemical products ([@bib0030]), although it has been noted that proportional inclusion of the active ingredients can be highly variable ([@bib0225], [@bib0005]). Anticoccidial vaccination using live *Eimeria* parasites is widespread in European breeder and layer stock, but remains uncommon in the majority broiler sector. Vaccines are also used in Egypt, Nigeria and India, but with a lower overall uptake than chemoprophylaxis and with a similar bias against broilers. Regional differences in climate may also be important, with the humid tropical climate in the southern regions sampled possibly favouring oocyst sporulation and persistence ([@bib0065], [@bib0015]), and thus encouraging co-infection. The expansion of surveillance to other northern and southern regions is required to test these hypotheses, although it is noteworthy that the UK and US reference strains tested here differed at just two PCR-RFLP SNPs, both of which were found to be variant within UK and Irish broilers. Interestingly, comparison with the UK reference strains genotyped during validation highlighted notably different profiles with three more SNPs found to be polymorphic than in the recent field samples (compare [Table 2](#tbl0010){ref-type="table"}, [Table 4](#tbl0020){ref-type="table"}, also Supplementary file 1). The *E. tenella* Houghton and Weybridge strains were both isolated more than 60 years ago ([@bib0055]), providing considerable scope for mutation and genome evolution.

Comparison of PCR-RFLP SNP profiles highlighted varied discriminatory power in different regions. In routine or pilot studies knowledge of the most informative SNPs offers the opportunity to screen using a streamlined, region-specific panel to reduce reagent costs and operator time. For genetic analysis additional markers are likely to be required. Comparison of the 52 SNP Sequenom and 11 SNP PCR-RFLP panels revealed equivalent haplotype numbers, mean genetic diversity and significance of association within each population, albeit with less focus on regional specificity. While NETWORK analysis was not fully resolved, the results achieved here support the validity of the full PCR-RFLP panel. Nonetheless, confirmation of the limited diversity discovered among UK and Irish *E. tenella* would benefit from adoption of additional markers given that just two of the eleven developed here proved informative.

5. Conclusion {#sec0080}
=============

As next generation sequencing technologies develop and become more cost-effective it is anticipated that metagenomics approaches focused on complex field samples will eventually become available for *Eimeria*. However, such tools remain distant for many veterinary pathogens. For now PCR-RFLP can begin to fill an important knowledge gap, exploring aspects of genetic diversity and population structure with direct relevance to anticoccidial control. Application of the technology to *Eimeria* from different countries and production systems can provide an important dataset whose value will be enhanced given the capture of key variables such as system type, chicken breed and age, date of sampling and farm choice of anticoccidial control.

Conflict of interest statement {#sec0085}
==============================

The authors declare that they have no competing interests.

Appendix A. Supplementary data {#sec0095}
==============================

The following is Supplementary data to this article:Supplementary Table S1Detailed PCR-RFLP results and accompanying sample data. Collection details for each sample analysed in these studies and the PCR-RFLP profiles derived.

The work described here was funded by the Biotechnology and Biological Sciences Research Council and the Department for International Development \[grant numbers BB/H009337 and BB/L00478X\]. The funding bodies played no role in the study design; collection, analysis and interpretation of data; writing of the report; and in the decision to submit the article for publication. The authors thank National Academy of Sciences of the United States of America for permission to reproduce a figure and table published previously for the purpose of review and comparison ([@bib0040]). This manuscript has been assigned the reference PPB_01375 by the RVC.

Supplementary data associated with this article can be found, in the online version, at [http://dx.doi.org/10.1016/j.vetpar.2016.09.018](10.1016/j.vetpar.2016.09.018){#intr0005}.

![Eleven genotyping PCR-RFLPs for *Eimeria tenella*, showing uncut/cut genotypes and negative control for each. The marker name and locus (in parentheses) is shown above each set.](gr1){#fig0005}

![Median-joining phylogenetic NETWORKs illustrating *Eimeria tenella* SNP haplotype diversity assessed using (A) 52 SNP Sequenom MassARRAY or (B) 11 SNP PCR-RFLP genotyping protocols. Both illustrate a regional bias in haplotype occurrence with a small number of highly represented haplotypes detected in North Africa and northern India, and a much larger number of mostly unique haplotypes in southern India and Nigeria. Haplotypes from the UK and Ireland were only available using PCR-RFLP and appear distinct from the other regions. Note the incomplete NETWORK resolution achieved using the smaller number of PCR-RFLPs, indicated by multiple cycles or cubes. Node size indicates the frequency of haplotype occurrence. Panel A reproduced from ([@bib0040]).](gr2){#fig0010}

###### 

Polymerase chain reaction-restriction fragment length (PCR-RFLP) markers for *Eimeria tenella*.

Table 1

  --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Marker      Target locus                                 Genome contig[a](#tblfn0005){ref-type="table-fn"}   Primers (5′ − 3′)          Amplicon (bp)   Restriction Enzyme (reaction temperature)   Target site \[SNP\]\   Digest fragments (bp)
                                                                                                                                                                                                      Underlined = cuts      
  ----------- -------------------------------------------- --------------------------------------------------- -------------------------- --------------- ------------------------------------------- ---------------------- -----------------------
  M3208       Microneme protein 3 (MIC3)                   HG675067                                            F: CACTTGAGTCCACTGCTCCA\   504             *Dde*I\                                     \[C/T\]TNAG            293/211
                                                                                                               R: TCATTGACAGCGACAAAAGC                    (37 °C)                                                            

  M3330       Microneme protein 3 (MIC3)                   HG675067                                            As M3208                   As M3208        *Mbo*II\                                    \[G/A\]AAGA            416/88
                                                                                                                                                          (37 °C)                                                            

  SG580       Surface antigen SAG_B05                      HG673747                                            F: TGAACCAGTTCAACGCAGAC\   935             *Bts*CI\                                    GGAT\[G/A\]            558/377
                                                                                                               R: CCTGATGGCCACTGAGAAAT                    (50 °C)                                                            

  RJ1409      Rhoptry protein J (RopJ)                     HG677946                                            F: GACTTGGTACGTTGGCCACT\   454             *Eco*RV-HF\                                 GA\[T/A\]ATC           315/139
                                                                                                               R: TTCCAATGTCCTTGCCTTTC                    (37 °C)                                                            

  M42850      Microneme protein 4 (MIC4)                   HG675525                                            F: GACATCGACGAATGTGCAAG\   612             *Hind*III-HF\                               A\[A/G\]GCTT           519/93
                                                                                                               R: GTACGAACCCGCAGTGTTTT                    (37 °C)                                                            

  M42986      Microneme protein 4 (MIC4)                   HG675525                                            As M42850                  As M42850       *Fau*I\                                     CCCG\[C/G\]            383/229
                                                                                                                                                          (55 °C)                                                            

  M15391      Microneme protein 1 (MIC1)                   HG673835                                            F: CGAACAGGACAAATGGTGTG\   617             *Eco*RI-HF\                                 GAATT\[C/A\]           386/231
                                                                                                               R: GTTGGGAGTCTGCACAGTGA                    (37 °C)                                                            

  BT1959      β-tubulin                                    HG673764                                            F: ATGCTTCCCCCTGAATCTTT\   642             *Bsr*BI\                                    CCGCT\[C/G\]           598/44
                                                                                                               R: TATTCTTCGCGGACTTTGCT                    (37 °C)                                                            

  C128593     Chromosome 1, feature-poor region 1 (C1P1)   HG675753                                            F: AATTAAAAGAAGGCGCAGCA\   889             *Bts*IMutI\                                 CA\[G/C\]TG            592/297
                                                                                                               R: CATACAAGCACCCAATGACG                    (55 °C)                                                            

  C1510864    Chromosome 1, feature-rich region 1 (C1R1)   HG675718                                            F: CTGCAGTTTGCTGCTTCTTG\   599             *Bts*IMutI\                                 CA\[G/A\]TG            544/55
                                                                                                               R: GAGGCCCCATTATAGCTTCC                    (55 °C)                                                            

  C2F235220   Chromosome 2, feature-poor region 1 (C2P1)   HG673763                                            F: GCACCAGCTGAACGTTTGTA\   520             *Bts*IMutI\                                 CA\[G/A\]TG            444/76
                                                                                                               R: GAGAGCACGAACAACAACGA                    (55 °C)                                                            
  --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Eimeria tenella Houghton strain genome, version 2013-11-05 (accessed through ToxoDB [@bib0100]).

###### 

PCR-RFLP profiles defining five reference and field *Eimeria tenella* isolates. I = PCR-RFLP cut profile, II = PCR-RFLP uncut profile.

Table 2

  Origin:     UK   UK   USA   Nigeria   Japan
  ----------- ---- ---- ----- --------- -------
  M3208       I    I    II    I         I
  M3330       I    I    I     II        I
  SG580       I    I    I     I         II
  RJ1409      II   II   II    II        I
  M42850      II   II   II    I         I
  M42986      I    I    I     II        II
  M15391      I    I    I     I         II
  BT1959      I    I    I     II        I
  C128593     II   II   II    II        I
  C1510864    II   II   II    I         I
  C2F235220   II   II   I     I         I

###### 

Summary of *Eimeria tenella* genome-wide genetic data calculated using previous Sequenom MassARRAY genotyping (including 52 SNPs per sample; top figure in each cell as reported previously ([@bib0040])) or PCR-RFLP (11 SNPs per sample; bottom *italicised figure* in each cell).

Table 3

                                       SNP haplotypes                                        
  --------------- ------------ ------- ---------------- ------ ------------------ ---------- ------------
  All (no UK)     Sequenom     244     93               --     0.361 (0.019)      0.124      \<0.001
                  *PCR-RFLP*   *244*   *48*             *--*   *0.351 (0.046)*    *0.164*    *\<0.001*
                                                                                             
  All (Inc. UK)   Sequenom     nd      nd               nd     nd                 nd         nd
                  *PCR-RFLP*   *271*   *52*             *--*   *0.388 (0.031)*    *0.140*    *\<0.001*
                                                                                             
  N. India        Sequenom     86      8                7      0.333 (0.025)      0.115      \<0.001
                  *PCR-RFLP*   *86*    *8*              *2*    *0.276 (0.069)*    *0.203*    *\<0.001*
                                                                                             
  S. India        Sequenom     53      50               49     0.340 (0.021)      0.009      ns
                  *PCR-RFLP*   *53*    *36*             *27*   *0.314 (0.040)*    *0.022*    *ns*
                                                                                             
  Egypt           Sequenom     40      21               13     0.192 (0.032)      0.058      \<0.001
                  *PCR-RFLP*   *40*    *5*              *1*    *0.135 (0.070)*    *0.072*    *\<0.001*
                                                                                             
  Libya           Sequenom     51      11               3      0.161 (0.025)      0.053      \<0.001
                  *PCR-RFLP*   *51*    *8*              *3*    *0.138 (0.057)*    *0.117*    *\<0.001*
                                                                                             
  Egypt + Libya   Sequenom     91      25               24     0.179 (0.027)      0.053      \<0.001
                  *PCR-RFLP*   91      9                5      0.139 (0.061)      *0.084*    *\<0.001*
                                                                                             
  Nigeria         Sequenom     14      14               13     0.400 (0.0224)     0.001      ns
                  *PCR-RFLP*   *14*    *11*             *1*    *0.3916 (0.052)*   *0.013*    *ns*
                                                                                             
  UK + Ireland    Sequenom     nd      nd               nd     nd                 nd         nd
                  *PCR-RFLP*   *27*    *4*              *4*    *0.068 (0.048)*    *0.0021*   *\<0.0001*

N = the number of samples tested. The number of SNP haplotypes was calculated using DnaSP v5.10.01. Het = mean genetic diversity and I~A~^S^ = standardised index of association, with statistical significance indicated as p, calculated using LIAN v3.6. nd = not done, ns = not statistically significant.

###### 

Summary of PCR-RFLP genotypes from field samples found to contain *Eimeria tenella*. The dominant SNP type detected is shown for each region, together with the percentage of samples defined by that SNP type.

Table 4

              Dominant SNP type (% samples)                                                                                                                                                                 
  ----------- ---------------------------------------------- ---------------------------------------------- ---------------------------------------------- ------------------------------------------------ ----------------------------------------------
  M3208       C (58.1)                                       T (75.5)                                       T (98.9)                                       C/T (50.0)[a](#tblfn0010){ref-type="table-fn"}   T (85.2)
  M3330       G (68.6)                                       G (60.3)[a](#tblfn0010){ref-type="table-fn"}   G (92.3)                                       G (78.6)                                         G (100.0)
  SG580       A (61.6)                                       G (64.1)                                       G (98.9)                                       A/G (50.0)[a](#tblfn0010){ref-type="table-fn"}   A (100.0)
  RJ1409      C (100.0)                                      C (75.5)                                       C (100.0)                                      C (92.9)                                         T (100.0)
  M42850      G (59.3)                                       G (88.7)                                       G (100.0)                                      G (71.4)                                         G (100.0)
  M42986      C (73.3)                                       G (92.4)                                       G (100.0)                                      C/G (50.0)[a](#tblfn0010){ref-type="table-fn"}   G (100.0)
  M15391      A (86.0)                                       C (73.6)                                       C (100.0)                                      A (71.4)                                         A (100.0)
  BT1959      G (100.0)                                      G (96.2)                                       G (100.0)                                      G (100.0)                                        C (100.0)
  C128593     C (52.3)[a](#tblfn0010){ref-type="table-fn"}   G (83.0)                                       G (64.8)                                       G (64.3)                                         C (100.0)
  C1510864    A (100.0)                                      A (75.5)                                       G (56.0)[a](#tblfn0010){ref-type="table-fn"}   A (78.6)                                         G (100.0)
  C2F235220   G (100.0)                                      G (84.9)                                       A (74.7)                                       G (71.4)                                         A (63.0)[a](#tblfn0010){ref-type="table-fn"}

Most discriminatory marker within each region (column).

[^1]: Present address: Imperial College London, South Kensington Campus, London, SW7 2AZ, UK.

[^2]: Present address: Roslin Institute, University of Edinburgh, Midlothian, UK.
