Development of genome-wide polymorphic microsatellite markers for Trichinella spiralis
Parasites & Vectors volume 13, Article number: 58 (2020)
Trichinella nematodes are globally distributed food-borne pathogens, in which Trichinella spiralis is the most common species in China. Microsatellites are a powerful tool in population genetics and phylogeographic analysis. However, only a few microsatellite markers were reported in T. spiralis. Thus, there is a need to develop and validate genome-wide microsatellite markers for T. spiralis.
Microsatellites were selected from shotgun genomic sequences using MIcroSAtellite identification tool (MISA). The identified markers were validated in 12 isolates of T. spiralis in China.
A total of 93,140 microsatellites were identified by MISA from 9267 contigs in T. spiralis genome sequences, in which 16 polymorphic loci were selected for validation by PCR with single larvae from 12 isolates of T. spiralis in China. There were 7–19 alleles per locus (average 11.25 alleles per locus). The observed heterozygosity (HO) and expected heterozygosity (HE) ranged from 0.325 to 0.750 and 0.737 to 0.918, respectively. The polymorphism information content (PIC) ranged from 0.719 to 0.978 (average 0.826). Among the 16 loci, markers for 10 loci could be amplified from all 12 international standard strains of Trichinella spp.
Sixteen highly polymorphic markers were selected and validated for T. spiralis. Primary phylogenetic analysis showed that these markers might serve as a useful tool for genetic studies of Trichinella parasites.
Human trichinellosis is caused by eating raw or undercooked meat infected with Trichinella parasites . Trichinella parasites have a broad geographical distribution on all continents except Antarctica, and can infect > 150 animal species, including mammals, birds and reptiles . The genus Trichinella contains nine species and three genotypes that can be separated into two clades by the ability to form encapsulated and non-encapsulated larvae [3,4,5]. There are genetic variations in Trichinella spp. based on geographical distributions and host species [6, 7]. In China, Trichinella spp. have been reported in a range of animals, including foxes, bears, wild boar, weasels, raccoon dogs, rats, bamboo rats and civets . Only two Trichinella species (i.e. T. spiralis and T. nativa) have been identified in China [8,9,10,11,12]. However, little is known about the genetic variations among the Trichinella species in China.
Genetic variability in T. spiralis was first reported in 1992, with three allozyme patterns at the loci of glucose 6-phosphate dehydrogenase and glucose phosphate isomerase detected in 61 isolates of T. spiralis from zoogeographical regions . Genetic polymorphisms in T. spiralis were also studied using different molecular tools, such as restriction fragment length polymorphism and single-strand conformational polymorphism (RFLP-SSCP) [13, 14], non-isotopic single-strand conformation polymorphism (‘cold’ SSCP) , and deep resequencing of the mitochondrial genomes . Compared with other molecular markers, microsatellites exist throughout the genome. In addition, microsatellites are relatively easy to score, since their gel band patterns could provide unambiguous results. Thus, they have been widely used in genetic diversity, population genetic structure, genome mapping, parentage analysis, population genetics and phylogeography studies [17,18,19]. However, only a few microsatellites have been reported in T. spiralis [12, 20,21,22]. The present study was aimed to identify and characterize microsatellites in T. spiralis and to obtain polymorphic microsatellite markers for further study.
Twelve isolates of T. spiralis were obtained from seven regions in China: five from Tianjin city, two from Yunnan Province, and one each from Heilongjiang, Henan, Hubei, Shaanxi and Tibet, respectively (Fig. 1). All isolates were confirmed as T. spiralis using multiplex PCR method according to Zarlenga et al. . The following 15 international standard Trichinella strains were acquired from the International Trichinella Reference Centre (ITRC; Rome, Italy): T. spiralis (T1, ISS534 and ISS4); T. nativa (T2, ISS70); T. britovi (T3, ISS100); T. pseudospiralis (T4, ISS13, ISS141 and ISS470); T. murrelli (T5, ISS415); Trichinella T6 (ISS34); T. nelsoni (T7, ISS37); Trichinella T8 (ISS124); Trichinella T9 (ISS408); T. papuae (T10, ISS572); T. zimbabwensis (T11, ISS1029); and T. patagoniensis (T12, ISS1826). All isolates and strains were maintained by serial passages in ICR mice. Larvae were recovered from the muscle tissues of infected mice on day 35 post-infection by an artificial digestion method , and stored at − 80 °C until use.
Microsatellite identification and primer design
All 9267 contigs of T. spiralis were retrieved from GenBank database (https://www.ncbi.nlm.nih.gov/nuccore/ABIR00000000) and used to search for microsatellite sequences by MIcroSAtellite Identification Tool (MISA) that was configured with strict minimum motif repeat requirements . The criteria of motifs were that mono- to hexanucleotide repeats with a minimum of 12 bp and a minimum of two repeat units. The maximum length of sequence between two simple sequence repeats (SSRs) to register as compound SSR was 100 bp . The number of microsatellites, motif, number of repeats, length of the repeat sequence, repeat type, start and end position of the repeat sequence, and microsatellite sequence, were analyzed using MISA.
Primers flanking the putative microsatellite sequences were designed at the PRIMER3 online server (http://primer3.ut.ee) , using following parameters: optimal primer length = 20 bp (between 18–22 bp); optimal primer GC content = 50% (between 40–60%); optimal primer melting temperature = 58 °C (between 55.9–60.1 °C); and product size ranged from 150 to 300 bp. The melting temperatures between a pair of primer had < 1 °C difference. The specificity of primer sequences was determined by BLAST searches against the genome of T. spiralis (https://www.ncbi.nlm.nih.gov/tools/primer-blast/).
Screening of microsatellites by PCR
A total of 1000 SSR primer pairs were selected for preliminary screening by PCR using DNA from a pool of ~ 4000 muscle larvae (~ 350 larvae from each of the 12 T. spiralis isolates in China). For isolating DNA, all larvae were homogenized in 500 μl extraction buffer containing 500 mM NaCl, 10 mM Tris-Cl (pH 8.0), 50 mM EDTA (pH 8.0), 2% (w/v) SDS and 10 mM β-mercaptoethanol, followed by incubation with 5 μl of proteinase K (20 mg/ml) at 60 °C for 0.5–2 h, phenol-chloroform extraction (50:50%, v/v), precipitation with 70% ethanol, and resuspension in 30–50 μl of sterile water. DNA samples were stored at − 20 °C. PCR reactions were carried out in a final volume of 20 μl, consisting of ~ 50 ng of DNA, 2 μl of 10× Ex Taq buffer (20 mM Mg2+ Plus; TaKaRa, Kusatsu, Japan), 1.6 μl of dNTP mixture (2.5 mM each), 0.2 μl of Ex Taq DNA polymerase (5 U/μl) (TaKaRa), and 0.4 μl of each primer (10 pmol/μl). PCR amplifications were performed in a thermal cycler (Applied Biosystems, California, USA) using following program: 98 °C for 5 min; followed by 35 cycles of 98 °C for 10 s, a specified annealing temperature for each primer pair for 30 s, 72 °C for 30 s; and a final extension step at 72 °C for 7 min. PCR products were electrophoresed on 1% agarose gels, stained with ethidium bromide and visualized under UV illumination. Microsatellite markers producing single bands were selected as candidate loci for further validation.
Verification of microsatellite polymorphism
Each of the selected primers was validated with 40 single larvae of T. spiralis from seven regions in China. Single larva was digested with proteinase K for DNA extraction using a Tissue and Hair Extraction Kit and a DNA IQ™ System Extraction Kit (Promega, Madison, USA) with magnetic beads following manufacturer’s instructions. DNA was eluted in 25 μl of elution buffer. Whole genome amplification was performed using an Illustra™ Ready-To-Go™ GenomiPhi V3 DNA Amplification Kit (GE Healthcare, Pittsburgh, USA) to increase the quantity of DNA. Concentrations of DNA were measured in a NanoDrop 2000 photometer (Thermo Fisher Scientific, Waltham, USA).
PCR amplifications were performed in a 20 μl reaction using a primer mixture which contained three primers: a sequence-specific forward primer with M13-tail at its 5′-end, a sequence-specific reverse primer, and the universal fluorescent-labeled M13 primer (FAM-M13 primer) . A 20 μl reaction contained 0.05 μM forward primer, 0.25 μM reverse primer, 0.2 μM FAM-M13 primer, 0.16 mM dNTP, 1 U of Ex Taq DNA polymerase (TaKaRa), and ~ 50 ng of DNA from a single larva . The PCR program was run as follows: 98 °C for 5 min; 32 cycles of 98 °C for 10 s, an annealing temperature specified for a primer pair for 30 s, and 72 °C for 30 s; eight additional cycles of 98 °C for 10 s, 53 °C for 30 s and 72 °C for 30 s; a final extension at 72 °C for 7 min. PCR products were subjected to capillary electrophoresis analysis (CEA) with a 96-capillary 3730XL DNA Analyzer (Applied Biosystems). Data were analyzed with GeneMapper 4.0 (Applied Biosystems). A negative control with sterile water was included in each PCR run.
Finally, the microsatellite loci with high polymorphism were selected for further validation by PCR using DNA samples isolated from individual larvae from 12 isolates of T. spiralis in China (10 larvae per isolate; total 120 samples). PCR amplification and analysis followed the protocols described above.
For each locus, the number of alleles (Na), the effective number of alleles (Ne), the expected heterozygosity (HE) and the observed heterozygosity (HO) per locus were estimated using GENEPOP version 4.2 (http://genepop.curtin.edu.au/) . This same software was used to test the polymorphism information content (PIC) and possible deviations from Hardy–Weinberg equilibrium (HWE) with Bonferroni correction .
DNA samples were isolated from the 12 Trichinella international standard strains as described in section “Screening of microsatellites by PCR” above. Cross-amplifications at selected polymorphic loci were performed and analyzed by a capillary electrophoresis using the same PCR protocols as described in section “Verification of microsatellite polymorphism” above.
The PCR products amplified from 15 international standard strains at the TsMs03 locus were analyzed by 8% denaturing urea-polyacrylamide gel electrophoresis. The homozygous individuals were selected for sequencing. Multiple sequence alignments of nucleotide sequences at the TsMs03 locus were performed using Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clustalo/) . The phylogenetic tree was inferred by MEGA X using the Neighbor-Joining method with 1000 bootstrap replicates [31, 32].
Abundance and microsatellite characteristics
A total of 93,140 microsatellites were identified from 9267 contigs of the T. spiralis genome by MISA (Table 1). The microsatellite density was 1591 loci per Mb. Among motifs containing mono- to hexanucleotide repeats, the most abundant was hexanucleotides that accounted for 49.51% of the total, followed by trinucleotide (19.61%) and tetranucleotide (17.44%). The di-, penta-, and mononucleotide motifs accounted for 8.77%, 3.69%, and 0.98% of the total motifs, respectively. The significant decrease in abundance of microsatellites was accompanied by the increase in the number of motif repeats. The number of repeating nucleotide sets was two times in 97.81% of hexanucleotide repeats. Meanwhile the number was three times in 1.81% of hexanucleotide repeats. For the pentanucleotide repeats, 68.29% consisted of three repeats, 19.12% consisted of four repeats, 8.18% consisted of five repeats, and 1.63% consisted of six repeats (Fig. 2). The top 20 most frequently classified repeat types were listed in Fig. 3. The most common motifs in each type of repeats were A/T (59.43%), AT/AT (61.84%), AAT/ATT (39.28%), AAAT/ATTT (37.30%), AAAAT/ATTTT (18.07%) and AAAAAT/ATTTTT (10.87%). The longest repeat was (TATAA)98 which belonged to the pentanucleotide group (Table 2).
Polymorphic microsatellite screening
Among the 1000 microsatellite loci selected for primary screening, 676 loci generated PCR products at expected sizes. A total of 120 loci producing single bright band in gel electrophoresis were selected as candidate loci. Among them, 47 microsatellite loci were homozygotes, while 57 loci showed low polymorphism. Finally, we selected 16 loci that produced distinct bands among individual larvae originated from different regions in China with high polymorphism for further analysis (Table 3).
Na varied from 7 to 19, and Ne ranged from 5.655 to 14.452 (average 8.820) per locus. HO and HE ranged from 0.325 to 0.750 and 0.737 to 0.918, respectively. PIC ranged from 0.719 to 0.978 (average of 0.826). The final set of 16 microsatellite markers were all highly informative (PIC > 0.50), and four of the 16 loci showed significant deviations from HWE after Bonferroni correction (Table 4).
Among the final 16 loci, 10 produced PCR amplicons for all tested Trichinella spp. Four (i.e. TsMs01, TsMs04, TsMs10 and TsMs14) obtained PCR products only from the Trichinella spp. with encapsulated larvae. Most of these loci were homozygous in the T. britovi (encapsulated larvae) and species with non-encapsulated larvae (Table 5). In addition, the TsMs07 and TsMs08 loci were amplified from species with encapsulated and non-encapsulated larvae, except for T. pseudospiralis. The average number of amplified alleles in each of the Trichinella spp. ranged from 1.300 (T. papuae and T. zimbabwensis) to 2.938 (Trichinella T9). A maximum of six alleles was observed in Trichinella T9 strain at the TsMs03 locus. Allelic size varied among taxa at a given locus, and one allele was shared by two or three taxa commonly. Trichinella T9 had specific alleles at three loci (i.e. TsMs12, TsMs14 and TsMs16) that were different in allelic size from other Trichinella taxa. None of the alleles at a given locus were shared by all Trichinella spp.
Primary phylogenetic analysis showed that all Trichinella spp. clustered into two clades: encapsulated larvae and non-encapsulated larvae group (Fig. 4). Sister relationship was observed for T. spiralis and T. nelsoni in comparison to other species with encapsulated larvae. Trichinella papuae and T. zimbabwensis were more closely related to each other than to T. pseudospiralis.
Microsatellites have been used in genetic diversity and genetic mapping studies in various organisms [33,34,35], partly because of their high polymorphism and the ability to detect alleles at a given locus in individual organisms [36, 37]. In previous studies, most of microsatellites in T. spiralis were designed based on expressed sequence tag (EST) databases [20,21,22]. The present study identified 93,140 microsatellites in the T. spiralis genomes using MISA, which accounted for 2.25% of the total genome sequence. The relative abundance of microsatellite sequences was estimated at 1.591 loci per kb of the T. spiralis genomes.
Generally, microsatellites decrease in abundance with increasing repeat length [38, 39], and this trend has been observed in many organisms . Previous comparative studies of microsatellites from eukaryotic genomes have found that the composition characteristics and distribution patterns significantly varied by species [39, 41]. Caenorhabditis elegans has a low frequency of microsatellites in its genome, even lower than Saccharomyces cerevisiae and other fungi [19, 42, 43]. In general, eukaryotic genomes are characterized by the prevalence of mononucleotide repeat motifs [19, 44]. For instance, mononucleotide repeats are the most abundant class of microsatellites in C. elegans  and Meloidogyne incognita . However, dinucleotide repeats are the most abundant type of motif in rodents  and most dicot plant species . Moreover, trinucleotide repeats are dominant in some algae and fungi species [44, 47], potentially indicating their genomic structural similarity with prokaryotes . In contrast, tetra- to hexanucleotide repeats are less abundant in eukaryotic genomes [49, 50]. Intriguingly, our results suggested a different distribution pattern for T. spiralis: hexa- > tri- > tetra- > di- > penta- > mononucleotide repeats. The repeat frequency of hexanucleotides (49.51%) was higher than other repeat classes. This may be a characteristic that is unique to T. spiralis. It is also possible that the abundance of repeats is influenced by secondary structures and DNA replication .
Among mononucleotide repeats, the motif (A/T)n is predominant, while (C/G)n repeats are rare [45, 48]. Our results for the most dominant motif type in mono- to hexanucleotide repeat classes of T. spiralis showed similar (A+T)-rich motif patterns, where A/T, AT, AAT, AAAT, AAAAT, and AAAAAT were the predominant repeats. The possible reasons for this (A + T)-rich motif pattern may be as follows: (A + T)-rich motifs can decrease the annealing temperature and accelerate strand separation, and the AT content increases through DNA replication and slippage . Secondly, DNA methylation can generate regions with high mutagenic rates, where the cytidine monophosphate becomes transformed into thymine. This type of mutation results from the deamination of methylation sites, leading to a combination of (A + T)-rich repeats. DNA methylation has been confirmed in the three life-cycle stages of T. spiralis, making it the only nematode species known to date with epigenetic modification of its genome . In addition, these repeats may be favored because the order of bases can directly influence chromatin structure, protein coding and gene function .
Previous studies have shown that Trichinella spp. are considered to have low intraspecific genetic diversity and genetic differentiation between populations [6, 21, 52,53,54,55,56,57,58]. The unique life-cycle of Trichinella species can often promote sibling inbreeding and reduced population size . Therefore, successful selection of microsatellite markers with relatively high abundance and polymorphism might be very difficult. Although the microsatellites of T. spiralis were detected in 12% of the 1000 EST sequences by La Rosa et al. , only seven microsatellite markers were suitable for genetic subgroup analysis. In the present study, 16 microsatellite markers with high polymorphism were selected and identified from 1000 candidate microsatellite loci.
To verify microsatellite markers with high polymorphism, we ranked the informativeness of markers using 120 individuals into highly (PIC > 0.50), reasonably (PIC of 0.25–0.50) and slightly informative (PIC < 0.25), as proposed by Botstein et al. . Sixteen markers with high PIC were selected in 12 isolates of T. spiralis in China. The number of alleles per locus were positively correlated with the length of the repeat region, such as the locus TsMs03, which had the highest number of alleles and the longest repeat sequence (TAATT)17. Previous studies have shown that long loci have higher mutation rates than short loci [36, 60]. The HWE describes how allele and genotype frequencies are related. Deviations often occur in the presence of small sample size, inbreeding, or the effects of population subdivision . Unfortunately, however, four microsatellite sites in tested populations deviated significantly from HWE after Bonferroni correction (P < 0.003) . In addition, HO was much lower than HE in these 16 loci, which led to the observation of limited polymorphism to some extent.
Zarlenga et al.  found that T. spiralis diverged early in the genus Trichinella. An analysis of population variability used nine microsatellite markers and observed more allelic richness among eight isolates originating in Asia compared to the remaining isolates from Europe, North Africa, and North and South America, suggested that T. spiralis populations are more diverse in East Asia, where pigs were first domesticated . Hence, in this study, we developed microsatellite loci and selected the ones with high polymorphism in 12 isolates of T. spiralis in China. The flanking sequences of the selected loci were relatively conserved in other Trichinella spp. Thus, ten of the 16 loci were amplified successfully in all 12 Trichinella spp. Therefore, the microsatellite loci developed in this study are good candidate loci to study the genetic variation and structure of Trichinella spp. beyond T. spiralis. Two loci, TsMs07 and TsMs08, were successfully amplified from all Trichinella spp., except for T. pseudospiralis. Recent studies have indicated that all five geographical isolates of T. pseudospiralis had one geographical origin that might diverge from T. papuae and T. zimbabwensis. Taken together, our results were consistent with other studies that T. papuae and T. zimbabwensis appeared to be basal in the group of species with non-encapsulated larvae and T. pseudospiralis the most recently evolved. The microsatellite analyses confirmed relationships among Trichinella spp. with non-encapsulated larvae, showing the utility of the new markers for investigating distantly related species within the genus .
We reported the identification of microsatellite sequences from the genome sequence data of T. spiralis with MISA. Among them, 16 microsatellites with high polymorphisms among 12 isolates of T. spiralis from various geographical regions in China were identified, and 10 microsatellites could be amplified successfully from all 12 Trichinella spp. The primary phylogenetic analysis suggested that the newly selected microsatellite markers could be applied to the analysis of genetic relationship of Trichinella spp. These microsatellite markers might serve as an important resource for the further study of Trichinella spp.
Availability of data and materials
The datasets supporting the findings of this article are included within the article.
MIcroSAtellite Identification Tool
- H O :
- H E :
Polymorphism information content
Restriction fragment length polymorphism
Single-strand conformation polymorphism
Simple sequence repeat
Capillary electrophoresis analysis
- N a :
Number of alleles
- N e :
Effective number of alleles
Expressed sequence tag
Pozio E, Zarlenga DS. New pieces of the Trichinella puzzle. Int J Parasitol. 2013;43:983–97.
Gottstein B, Pozio E, Nockler K. Epidemiology, diagnosis, treatment, and control of trichinellosis. Clin Microbiol Rev. 2009;22:127–45.
Pozio E, La Rosa G, Murrell KD, Lichtenfels JR. Taxonomic revision of the genus Trichinella. J Parasitol. 1992;78:654–9.
Murrell KD, Lichtenfels RJ, Zarlenga DS, Pozio E. The systematics of the genus Trichinella with a key to species. Vet Parasitol. 2000;93:293–307.
Pozio E, Zarlenga DS. Recent advances on the taxonomy, systematics and epidemiology of Trichinella. Int J Parasitol. 2005;35:1191–204.
La Rosa G, Pozio E, Rossi P, Murrell KD. Allozyme analysis of Trichinella isolates from various host species and geographical regions. J Parasitol. 1992;78:641–6.
Pozio E. Factors affecting the flow among domestic, synanthropic and sylvatic cycles of Trichinella. Vet Parasitol. 2000;93:241–62.
Wang ZQ, Cui J, Shen LJ. The epidemiology of animal trichinellosis in China. Vet J. 2007;173:391–8.
Liu MY, Zhu XP, Xu KC, Lu Q, Boireau P. Biological and genetic characteristics of two Trichinella isolates in China; comparison with European species. Parasite. 2001;8(Suppl. 2):S34–8.
Gasser RB, Zhu XQ, Monti JR, Dou L, Cai X, Pozio E. PCR-SSCP of rDNA for the identification of Trichinella isolates from mainland China. Mol Cell Probes. 1998;12:27–34.
Liu M, Boireau P. Trichinellosis in China: epidemiology and control. Trends Parasitol. 2002;18:553–6.
Fu BQ, Liu MY, Yao CY, Li WH, Li YG, Wang YH, et al. Species identification of Trichinella isolates from China. Vet Parasitol. 2009;159:214–7.
Wu Z, Nagano I, Pozio E, Takahashi Y. Polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) for the identification of Trichinella isolates. Parasitology. 1999;118:211–8.
Wu Z, Nakada T, Nagano I, Takahashi Y. DNA fingerprints of Trichinella as revealed by restriction fragment length polymorphism and single-strand conformational polymorphism (RFLP-SSCP). Mol Cell Probes. 2000;14:291–7.
Gasser RB, Hu M, El-Osta YA, Zarlenga DS, Pozio E. Genetic analysis of Trichinella populations by ‘cold’ single-strand conformation polymorphism analysis. Vet Parasitol. 2005;132:23–6.
Webb KM, Rosenthal BM. Deep resequencing of Trichinella spiralis reveals previously un-described single nucleotide polymorphisms and intra-isolate variation within the mitochondrial genome. Infect Genet Evol. 2010;10:304–10.
Guichoux E, Lagache L, Wagner S, Chaumeil P, Leger P, Lepais O, et al. Current trends in microsatellite genotyping. Mol Ecol Resour. 2011;11:591–611.
Hodel RG, Segovia-Salcedo MC, Landis JB, Crowl AA, Sun M, Liu X, et al. The report of my death was an exaggeration: a review for researchers using microsatellites in the 21st century. Appl Plant Sci. 2016;4:1600025.
Toth G, Gaspari Z, Jurka J. Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res. 2000;10:967–81.
Rosenthal BM, LaRosa G, Zarlenga D, Dunams D, Chunyu Y, Mingyuan L, et al. Human dispersal of Trichinella spiralis in domesticated pigs. Infect Genet Evol. 2008;8:799–805.
La Rosa G, Marucci G, Rosenthal BM, Pozio E. Development of a single larva microsatellite analysis to investigate the population structure of Trichinella spiralis. Infect Genet Evol. 2012;12:369–76.
Wang ZQ, Li LZ, Jiang P, Liu LN, Cui J. Molecular identification and phylogenetic analysis of Trichinella isolates from different provinces in mainland China. Parasitol Res. 2012;110:753–7.
Zarlenga DS, Chute MB, Martin A, Kapel CM. A multiplex PCR for unequivocal differentiation of all encapsulated and non-encapsulated genotypes of Trichinella. Int J Parasitol. 1999;29:1859–67.
Wu XP, Fu BQ, Wang XL, Yu L, Yu SY, Deng HK, et al. Identification of antigenic genes in Trichinella spiralis by immunoscreening of cDNA libraries. Vet Parasitol. 2009;159:272–5.
Thiel T, Michalek W, Varshney RK, Graner A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003;106:411–22.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115.
Schuelke M. An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000;18:233–4.
Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89:583–90.
Rice WR. Analyzing tables of statistical tests. Evolution. 1989;43:223–5.
Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.
Kumar S, Nei M, Dudley J, Tamura K. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008;9:299–306.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–9.
Gower CM, Gabrielli AF, Sacko M, Dembele R, Golan R, Emery AM, et al. Population genetics of Schistosoma haematobium: development of novel microsatellite markers and their application to schistosomiasis control in Mali. Parasitology. 2011;138:978–94.
Laoprom N, Sithithaworn P, Ando K, Sithithaworn J, Wongkham S, Laha T, et al. Microsatellite loci in the carcinogenic liver fluke, Opisthorchis viverrini and their application as population genetic markers. Infect Genet Evol. 2010;10:146–53.
Shrivastava J, Barker GC, Johansen MV, Xiaonong Z, Aligui GD, McGarvey ST, et al. Isolation and characterization of polymorphic DNA microsatellite markers from Schistosoma japonicum. Mol Ecol Notes. 2003;3:406–8.
Ellegren H. Heterogeneous mutation processes in human microsatellite DNA sequences. Nat Genet. 2000;24:400–2.
Selkoe KA, Toonen RJ. Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers. Ecol Lett. 2006;9:615–29.
Temnykh S, DeClerck G, Lukashova A, Lipovich L, Cartinhour S, McCouch S. Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): frequency, length variation, transposon associations, and genetic marker potential. Genome Res. 2001;11:1441–52.
Grover A, Aishwarya V, Sharma PC. Biased distribution of microsatellite motifs in the rice genome. Mol Genet Genomics. 2007;277:469–80.
Karaoglu H, Lee CM, Meyer W. Survey of simple sequence repeats in completed fungal genomes. Mol Biol Evol. 2005;22:639–49.
Sonah H, Deshmukh RK, Sharma A, Singh VP, Gupta DK, Gacche RN, et al. Genome-wide distribution and organization of microsatellites in plants: an insight into marker development in Brachypodium. PLoS ONE. 2011;6:e21298.
Katti MV, Ranjekar PK, Gupta VS. Differential distribution of simple sequence repeats in eukaryotic genome sequences. Mol Biol Evol. 2001;18:1161–7.
Astolfi P, Bellizzi D, Sgaramella V. Frequency and coverage of trinucleotide repeats in eukaryotes. Gene. 2003;317:117–25.
Sharma PC, Grover A, Kahl G. Mining microsatellites in eukaryotic genomes. Trends Biotechnol. 2007;25:490–8.
Castagnone-Sereno P, Danchin EG, Deleury E, Guillemaud T, Malausa T, Abad P. Genome-wide survey and analysis of microsatellites in nematodes, with a focus on the plant-parasitic species Meloidogyne incognita. BMC Genomics. 2010;11:598.
Kumpatla SP, Mukhopadhyay S. Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species. Genome. 2005;48:985–98.
Kim TS, Booth JG, Gauch HG Jr, Sun Q, Park J, Lee YH, et al. Simple sequence repeats in Neurospora crassa: distribution, polymorphism and evolutionary inference. BMC Genomics. 2008;9:31.
Qi WH, Jiang XM, Du LM, Xiao GS, Hu TZ, Yue BS, et al. Genome-wide survey and analysis of microsatellite sequences in bovid species. PLoS ONE. 2015;10:e0133667.
Subramanian S, Mishra RK, Singh L. Genome-wide analysis of microsatellite repeats in humans: their abundance and density in specific genomic regions. Genome Biol. 2003;4:R13.
Xu Y, Li W, Hu Z, Zeng T, Shen Y, Liu S, et al. Genome-wide mining of perfect microsatellites and tetranucleotide orthologous microsatellites estimates in six primate species. Gene. 2018;643:124–32.
Gao F, Liu XL, Wu XP, Wang XL, Gong DS, Lu HL, et al. Differential DNA methylation in discrete developmental stages of the parasitic nematode Trichinella spiralis. Genome Biol. 2012;13:R100.
Dunams-Morel DB, Reichard MV, Torretti L, Zarlenga DS, Rosenthal BM. Discernible but limited introgression has occurred where Trichinella nativa and the T6 genotype occur in sympatry. Infect Genet Evol. 2012;12:530–8.
Dame JB, Murrell KD, Worley DE, Schad GA. Trichinella spiralis: genetic evidence for synanthropic subspecies in sylvatic hosts. Exp Parasitol. 1987;64:195–203.
La Rosa G, Pozio E, Rossi P. Biochemical resolution of European and African isolates of Trichinella nelsoni Britov and Boev, 1972. Parasitol Res. 1991;77:173–6.
La Rosa G, Pozio E. Molecular investigation of African isolates of Trichinella reveals genetic polymorphism in Trichinella nelsoni. Int J Parasitol. 2000;30(5):663–7.
La Rosa G, Marucci G, Zarlenga DS, Pozio E. Trichinella pseudospiralis populations of the Palearctic region and their relationship with populations of the Nearctic and Australian regions. Int J Parasitol. 2001;31:297–305.
La Rosa G, Marucci G, Zarlenga DS, Casulli A, Zarnke RL, Pozio E. Molecular identification of natural hybrids between Trichinella nativa and Trichinella T6 provides evidence of gene flow and ongoing genetic divergence. Int J Parasitol. 2003;33:209–16.
Cole R, Viney M. The population genetics of parasitic nematodes of wild animals. Parasit Vectors. 2018;11:590.
Botstein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet. 1980;32:314–31.
Schlotterer C. Evolutionary dynamics of microsatellite DNA. Chromosoma. 2000;109:365–71.
Lachance J, Tishkoff SA. SNP ascertainment bias in population genetic analyses: why it is important, and how to correct it. Bioessays. 2013;35:780–6.
Groenewald CW, Born C, Chown SL, van Vuuren BJ. Development of a microsatellite library for the flightless moth Pringleophaga marioni Viette (Lepidoptera: Tineidae). Conserv Genet Resour. 2011;3:291–4.
Zarlenga DS, Rosenthal BM, La Rosa G, Pozio E, Hoberg EP. Post-Miocene expansion, colonization, and host switching drove speciation among extant nematodes of the archaic genus Trichinella. Proc Natl Acad Sci USA. 2006;103:7354–9.
Korhonen PK, Pozio E, La Rosa G, Chang BC, Koehler AV, Hoberg EP, et al. Phylogenomic and biogeographic reconstruction of the Trichinella complex. Nat Commun. 2016;7:10513.
This study was supported by The National Key Research and Development Program of China (2017YFD0501300), the National Natural Science Foundation of China (31520103916, 31872467), Guangdong Innovative and Entrepreneurial Research Team Program (No. 2014ZT05S123) and Program for JLU Science and Technology Innovative Research Team.
Ethics approval and consent to participate
All mice were handled strictly in accordance with the Animal Ethics Procedures and Guidelines of the People’s Republic of China. The protocol was approved by the Institutional Animal Care and Use Committee of Jilin University (Protocol # 20170318).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, TT., Tang, B., Bai, X. et al. Development of genome-wide polymorphic microsatellite markers for Trichinella spiralis. Parasites Vectors 13, 58 (2020). https://doi.org/10.1186/s13071-020-3929-2