- Research article
- Open Access
Genetics of lineage diversification and the evolution of host usage in the economically important wheat curl mite, Aceria tosichella Keifer, 1969
BMC Evolutionary Biology volume 18, Article number: 122 (2018)
Understanding the mechanisms that underlie the diversification of herbivores through interactions with their hosts is important for their diversity assessment and identification of expansion events, particularly in a human-altered world where evolutionary processes can be exacerbated. We studied patterns of host usage and genetic structure in the wheat curl mite complex (WCM), Aceria tosichella, a major pest of the world’s grain industry, to identify the factors behind its extensive diversification.
We expanded on previous phylogenetic research, demonstrating deep lineage diversification within the taxon, a complex of distinctive host specialist and generalist lineages more diverse than previously assumed. Time-calibrated phylogenetic reconstruction inferred from mitochondrial DNA sequence data suggests that lineage diversification pre-dates the influence of agricultural practices, and lineages started to radiate in the mid Miocene when major radiations of C4 grasses is known to have occurred. Furthermore, we demonstrated that host specificity is not phylogenetically constrained, while host generalization appears to be a more derived trait coinciding with the expansion of the world’s grasslands. Demographic history of specialist lineages have been more stable when compared to generalists, and their expansion pre-dated all generalist lineages. The lack of host-associated genetic structure of generalists indicates gene flow between mite populations from different hosts.
Our analyses demonstrated that WCM is an unexpectedly diverse complex of genetic lineages and its differentiation is likely associated with the time of diversification and expansion of its hosts. Signatures of demographic histories and expansion of generalists are consistent with the observed proliferation of the globally most common lineages. The apparent lack of constrains on host use, coupled with a high colonization potential, hinders mite management, which may be further compromised by host range expansion. This study provides a significant contribution to the growing literature on host-association and diversification in herbivorous invertebrates.
Niche breadth varies significantly among species, with life history traits, physical tolerances, and evolutionary adaptations directly influencing the range of habitats which species can successfully colonize and thrive within [1, 2]. Consequently, generalist and specialist species differ in traits and adaptations that facilitate the inhabitation of wide and narrow environmental ranges, respectively [3,4,5]. Species that have obligatory associations with animal or plant hosts, which function as their local habitat, provide ideal model systems for studying niche breadth and host-associated genetic differentiation [6,7,8,9]. Long-term host affiliation in the absence of gene flow across populations occupying different hosts can lead to their genetic divergence due to host specialization [6,7,8, 10]. Yet, these processes are not static, and host ranges can change dynamically through time with environmental changes, specialists giving rise to generalists and vice versa. This nature of host range expansion and contraction is thought to be largely responsible for the astounding diversity of parasite and herbivorous invertebrate taxa that exists today [10,11,12].
Host-associations can change both markedly and rapidly under different scenarios . Populations of generalist species might restrict their host range as a result of diminishing host availability or costs of host adaptation . Conversely, specialist species might expand their host range or transit from one exclusive host to another (host shift) as a consequence of environmental change, including anthropogenic influences [13, 14]. Numerous studies of host-associated differentiation have contributed greatly to our understanding of host shifts and host range evolution [15,16,17,18,19,20], with the sympatric host race formation in the apple maggot fly Rhagoletis pomonella being the most well-known example . Yet there is still much to learn about mechanisms behind the diversification of herbivores and parasites through interactions with their hosts [22, 23].
The host specificity is among the most important life-history traits that influences population size and genetic structure [9, 24]. Studies on parasitic lice in birds  and nematodes in livestock  indicate that the degree of genetic structuring in host specialists will typically be higher and their genetic diversity will be reduced. With an exception of cases where hosts are highly prolific, gene flow in specialists is often more restricted due to reduced dispersal capacity, population size, or host availability . Yet these patterns have been established with grounds on observations based on laboratory models thus far, and therefore it is still unclear how broadly this concept applies to natural systems.
Understanding patterns of genetic diversity and host ranges is particularly important for identifying risks associated with parasites and herbivorous pests, especially those with high invasive potential [27, 28]. The wheat curl mite Aceria tosichella Keifer (WCM hereafter) is an obligate plant-feeding pest in wheat and many other cereal crops, with a high colonization potential driven by an ability for long distance dispersal and reproduction by arhenotokous parthenogenesis . Mite infestations can lead to yield losses through direct feeding and the transmission of several plant viruses [30, 31]. Wheat streak mosaic virus (WSMV) is a particularly damaging viral pathogen that primarily affects wheat in cropping regions of the Americas, Northern Africa, Asia, Europe and Australia, and can be responsible for annual yield losses of approximately 5%, with localized areas suffering complete yield loss [32,33,34]. Recent studies proved that the WCM is a complex of distinct mitochondrial lineages that differ in their host preference [35,36,37,38], with two globally distributed lineages (viz. MT-1 and MT-8) that appear to be responsible for global WSMV transmission [39, 40]. Recent sampling efforts suggest that lineage diversity within the WCM complex is currently understated, and further work is needed to characterize the true extent of genetic diversity and host ranges of the lineages [41, 42].
Of the WCM lineages described to date, all appear to be morphologically indistinguishable, but differ markedly in host-acceptance traits, including host specialists and generalists [37, 38, 42]. Thus, this is an ideal system to investigate host-association dynamics and associated evolutionary processes. In this study, we tested the hypotheses: (i) whether host specificity is a phylogenetically constrained trait, and if mite lineages are restricted to hosts of certain taxonomic groups, suggesting possible co-evolution; (ii) if there is an association between the timing of lineage diversification and the historical expansion of agriculture; and (iii) if there is host-associated structuring within generalist lineages, that may be indicative of incipient speciation. We drew on all available mitochondrial Cox1 and nuclear 28S D2 sequence data, representing mites from 25 host-plant species and undertook comprehensive phylogenetic and population genetic analyses to provide insight into the extent of lineage diversification within the WCM complex, and host ranges of respective lineages. We compared patterns of genetic structure between host specialists and generalists, and identified patterns of demographic expansion or stability of host specialists and generalists.
Sampling of mites
The cereal hosts (bread wheat Triticum aestivum L., triticale, ×Triticosecale Wittm. ex A. Camus, rye Secale cereale L., oat Avena sativa L., barley Hordeum vulgare L., maize Zea mays L.) and wild grass hosts (cockspur Echinochloa crus-galli (L.) Beauv, quackgrass Elymus repens (L.) Gould, tall oat-grass Arrhenatherum elatius (L.) P. Beauv. ex J. & C. Presl, soft brome Bromus hordeaceus L., smooth brome Bromus inermis Leyss., timothy-grass Phleum pretense L., and wall barley Hordeum murinum L.) of the wheat curl mite were collected in the field from a total of 316 sample sites in Poland (covering an area of 311,888 km2) between June and August 2012–2014 (see Additional file 1: Table S1). Plants were examined in the laboratory for the presence of WCM specimens under a stereo-microscope. Individual specimens were collected using insect pin and preserved in Eppendorf tubes with 180 μl of ATL extraction-buffer (Qiagen GmbH, Hilden, Germany) for subsequent DNA extraction and genetic analysis. Altogether 1187 mite specimens were collected (from 1 to 20 specimens per one Eppendorf tube) and subsequently analyzed.
DNA isolation, amplification and sequencing
DNA was isolated from specimens that had been stored in ATL buffer according to the non-destructive method of Dabert et al. . The exoskeletons of the digested mites were preserved in 70% ethyl alcohol, and later mounted on slides according to Monfreda et al.  for morphological WCM identification. A fragment of the mitochondrial cytochrome c oxidase subunit I (Cox1) gene (603 bp) was amplified by PCR using the degenerate primers bcdF01 and bcdR04 [43, 45]. Amplification of the ca. 600 bp D2 region in 28S rDNA was performed using the primers D1D2fw2  and 28SR0990 . Reactions steps and product handling followed protocols described by Szydło et al. . Products were sequenced with the same primers that were used for amplification and additionally D2 amplified fragments were sequenced with specific sequencing primers Er28SF05 and Er28SR05 . Trace files were checked and edited using MEGA 6 .
Mitochondrial and nuclear datasets
To be able to build an alignment consisting of sequences of comparable length (and therefore comprehensive datasets), we obtained sequences of similar length to these downloaded from Genbank. Newly generated 662 mitochondrial and 63 nuclear sequence data from WCM in Poland were combined with previously published 85 Cox1 and 46 28S rDNA D2 regions of WCM and outgroup Trisetacus species data [36, 39, 48, 50] for analysis (see Additional file 1: Table S1 for sample details and Genbank accession numbers). Trisetacus genus is a taxon belonging to the Eriophyoidae associated with coniferous plants, for which genetic differentiation according to host plants has been recorded . Sequences were aligned using the MAFFT algorithm (https://www.ebi.ac.uk/Tools/msa/mafft/). Altogether 747 Cox1 (662 + 85; 603 bp) and 109 28S D2 (63 + 46; 595 bp) sequences were analyzed and evaluated. They were obtained from mite specimens from 25 plant species (13 plant species collected in Poland in the course of this study plus 12 plant species from earlier sample collections) and eight countries (Additional file 1: Table S1). Sequences from both mitochondrial and nuclear datasets were collapsed to haplotypes using FABOX 1.41 DNA collapse tool  resulting in a batch of 291 unique Cox1 haplotypes, 39 unique D2 sequences, and 84 unique concatenated Cox1 + D2 sequences. Alignment reliability and product authenticity of Cox1 sequences was validated by translating aligned DNA sequences into amino acids and assessing the alignments for premature stop codons. A simplified Cox1 sequence dataset with a representative set of 158 sequences (for each WCM lineage a maximum of five sequences from the same country and same host) was used for the construction of phylogenetic trees. This selection was based on preliminary trees and aimed to simplify the phylogenetic tree as well as to keep the general tree structure. Cox1 and D2 sequences obtained from the same DNA isolate were concatenated resulting in a dataset consisting of 109 concatenated haplotypes that were collapsed to 84 unique concatenated sequences. Partition congruence was analyzed using the incongruence length difference test (ILD) with 100 partition homogeneity test replicates implemented in PAUP* 4.0a147 .
Phylogenetic analyses and divergence time estimation
Phylogenetic reconstructions were performed using Bayesian Inference (BI) methods implemented in BEAST 2.3.0 . General Time Reversible model  with gamma distribution of rates across sites (GTR + G) was selected as the best fit model of evolution for each of the mtDNA and nuclear genes, based on Akaike Information Criteria (AIC)  implemented in JMODELTEST v.0.1.1 . Operators were auto-optimized, and five independent Markov Chain Monte Carlo (MCMC) runs were performed using a Yule (speciation) tree-prior, each running for 5 × 106 generations, sampling every 5000 states. Log files were examined with TRACER v.1.5  to ensure that runs were sampling from the same posterior distribution, to determine appropriate burn-in, and to ensure that effective sample sizes (ESSs) of parameters of interest were greater than 1000. Tree files of independent runs were then combined with LOGCOMBINER v.2.1.3 , discarding the first 20% and re-sampling at lower frequency of 15,000. The maximum clade credibility (MCC) tree was recovered from a sample of 10,000 posterior trees, and branch support was annotated, using TREEANNOTATOR v.2.1.3 . All analyses started with a random starting tree and seed with no root specified. Sequence data from Trisetacus species was used to estimate the root of the mitochondrial gene tree. Additionally, Maximum Likelihood (ML) analysis was performed in PAUP* with starting trees obtained by Neighbor-joining, and Tree bisection-reconnection (TBR) as branch swapping algorithm. Bootstrap proportions  were obtained to access node robustness, using 1000 bootstrap replications.
In order to test the timing of diversification between WCM mitochondrial lineages, the mitochondrial gene tree was time calibrated, with divergence times of nodes being inferred from 95% highest posterior density (HPD) intervals. Substitution rates for the Cox1 locus in mite lineages have been shown to differ by up to four orders of magnitude, and at present there is no calibrated substitution rate for Acari specifically. Consequently, we chose to work with the universal invertebrate substitution rate 0.0115 substitutions/site/million years [60, 61] with relaxed clock log normal priors with standard deviations account allowing for uncertainty (four orders of magnitude) around the mean. Substitution rates were set in BEAUti 1.7.3 , and TRACER was then used to obtain parameter estimates for time to the most recent common ancestor (tMRCAs) for nodes within the gene tree.
DNA sequence-based species delineation
We implemented two DNA taxonomy approaches to evaluate the presence of additional putative mitochondrial lineages using the Cox1 dataset. (1) The general mixed Yule coalescent (GMYC) approach [62, 63] was applied to the ultrametric tree (produced by beast analyses) in R 2.15.3  with the ‘splits package’. The GMYC model is a process-based approach for detecting the threshold in a gene tree at which within-species processes (i.e., coalescence) shift to between-species processes (i.e., speciation and extinction). (2) We applied the combination of the Poisson Tree Processes model for species delimitation (PTP), and a Bayesian implementation of PTP (bPTP) to infer putative lineages boundaries on a given phylogenetic input tree . The PTP/bPTP model, unlike the GMYC model, requires a bifurcated phylogenetic tree, not an ultrametric tree. PTP/bPTP models speciation or branching events in terms of number of substitutions. We used the following parameters: MCMC 500,000 generations; Thinning, 100; Burnin, 0.1; Seed, 123, and assessed convergence in each case to ensure reliability of results. Additionally, the mean uncorrected p-distances with standard error estimates (obtained using a bootstrap procedure with 1000 replicates) within and between WCM lineages were calculated in MEGA6 .
Diversity within WCM lineages and tests for demographic history
Genetic diversity estimates including number of haplotypes (h), number of polymorphic sites (p), haplotype diversity (Hd) and nucleotide diversity (л) were calculated with ARLEQUIN v.220.127.116.11  for each major clade inferred from BI analyses. We investigated the historical demography of the WCM lineages with N ≥ 15 by calculating Tajima’s D  and Fu’s Fs . The significance of Fs statistic was determined with a simulated random distribution produced by 1000 permutations assuming neutrality and population equilibrium. We also analyzed the mismatch distribution of pairwise genetic differences . The fit to models of population expansion was tested with the sum of squared deviations between the observed and expected mismatch (SSD) and the raggedness index (HRag). The significance of HRag and SSD was determined with 1000 bootstrap replicates. Non-significant values for these statistics signify that the data do not deviate from that expected under the model of expansion. All above analyses were conducted with ARLEQUIN. Estimates of times for lineage expansion were calculated by t = τ/2 μ where τ is the number of generations and μ is the cumulative (across the sequence) probability of substitution, assuming a divergence rate of 0.0115 substitutions per nucleotide per million years and 18 generations per year.
Tests for host associated genetic structure in multi-host WCM lineages
For WCM lineages with multi-host association, we tested for patterns of host associated genetic structure. An Analysis of Molecular Variation (AMOVA) was performed to test for significant host-associated genetic structuring, indicating ecological isolation and potential incipient speciation. Partitioning of WCM mtDNA variation according to hosts was investigated by AMOVA using pairwise FST as the distance measure with 10,000 permutations, as implemented in ARLEQUIN. We also tested for host associated structuring by comparing relatedness estimates among haplotypes inferred from Parsimony Median Joining Network analyses  implemented in popart 1.7. FST values between host-associated populations (with n ≥ 10) were computed to disclose the level of sharing genetic diversity. Additionally, nucleotide diversity (л) and haplotype diversity (Hd) were calculated for each host population, within multi-host lineages, and were compared with the estimates for the whole lineage.
Divergence between WCM lineages
Phylogenetic analyses and molecular dating
The Bayesian Inference (BI) analysis WCM Cox1 sequences revealed a high level of genetic structuring with strong statistical support (posterior probability, PP > 0.95) for the monophyly of 20 major clades (hereafter lineages) and nine lineages represented by single sequences (Fig. 1a and Additional file 2: Figure S1). Trees produced by Maximum Likelihood analysis showed consistent topologies (Additional file 2: Figure S2). Uncorrected pairwise distances between WCM lineages ranged from 6.7 to 28.2% (Additional file 3: Table S2). Seven clades were specialists being associated with a single wild grass species, including lineages associated with smooth-brome (clades: MT-9, MT-10, MT-14, violet color on the Fig. 1a), wall barley (MT-7, light blue), tall oat-grass (MT-5, orange), quackgrass (MT-6, green) and timothy-grass (MT-16, brown) hosts. Two lineages were found on two host species but were mainly associated with one of them, hereafter named ‘semi-specialists’; these included MT-12 and MT-13 (black). Finally, six lineages were found on multiple host species (black, MT-1, MT-2, MT-3, MT-4, MT-8, MT-23), hereafter named as ‘generalists’. Note, that this classification was applied only for lineages with at least three records. The relationships among clades were largely resolved (PP > 0.95). Furthermore, evidence of paraphyly in mites which are specialized toward a particular host was observed, such as in mites associated with smooth brome (i.e. MT-9, MT-10 and MT-14). In contrast, the specialist lineages MT-5 and MT-7 each associated with a different host species (tall oat-grass and wall barley, respectively) formed a well-supported sister-relationship (PP > 0.95). There was also no support for genetic grouping for lineages associated with cereal or wild growing grass hosts, however lineages occurring on cereals were the most polyphagous. Additionally, WCM lineages associated with taxonomically related plant species, namely congeneric species, did not group together (Fig. 1a).
BI reconstructions of WCM interrelationships based on the D2 dataset (Fig. 2) supported high level of genetic structuring within WCM, however it failed to support some lineage sorting observed in mtDNA analysis. The D2 tree was only partly coinciding with the mitochondrial tree (ILD tests indicated phylogenetic incongruence between the Cox1 and D2, P = 0.01): some relationships were similar to those observed in the Cox1 gene tree (e.g. clustering of MT-6 and MT-16; MT-12 and MT-29 with PP > 0.95), but some others were not supported on D2 tree (e.g. between lineages MT-3, MT-13 and MT-24; MT-8 and MT-17). Branch lengths were considerably shorter compared to those within the Cox1 gene tree, but expected due to the comparatively slower evolutionary rate of the locus . The concatenated gene tree, that included a total of 84 unique haplotypes, recovered the major WCM clades established by the mitochondrial tree (PP > 0.95) (Additional file 2: Figure S3).
Our time calibrated Cox1 phylogeny indicates that the tMRCA for all genetic lineages described in this study is approximately 17.97 mya (16.27–19.60 mya). Two of the most genetically distinct WCM lineages (viz. MT-13 and MT-15) diverged approximately 16.78 mya (15.07–18.51 mya) and 13.77 mya (12.50–15.21 mya). Other WCM lineages appear to have started to radiate approximately 11.18 mya (10.15–12.14 mya) in mid-Miocene epoch and continued diversification after ca. 7.5 mya. Interestingly, two of the most polyphagous lineages, associated mostly with cereals (viz. MT-1 and MT-8), diverged approximately at the same time from different ancestors, i.e. ca. 3 mya (2.13–3.62 mya). Two specialist lineages (MT-5 and MT-7) shared a common ancestor with generalists MT-8 at 5.38 mya (4.60–6.19 mya), and they both diverged from their common ancestor at 4.33 mya (3.64–5.17 mya) (Fig. 3). Nevertheless, the above ages of radiation events within the WCM complex should be treated with caution due to the large 95% HPDs.
DNA sequence-based species delineation
Species delineation analyses (GYMC and PTP/bPTP methods) (Fig. 1b) were consistent with phylogenetic structuring by identifying each of the 20 major WCM clades as putative species. Additional delineation was detected within eight lineages (MT-1, MT-3, MT-6, MT-9, MT-12, MT-17, MT-18, and MT-26) including specialists, generalists and semi-specialists, in total identifying 45 (41–50) lineages (PTP/bPTP analysis) and 42 (95% CI: 39–50) lineages (GYMC analysis). Lineages assignments were highly consistent between the two methods; however, due to the varying degree of differentiation between lineages, only those identified by PTP/bPTP are shown in Fig. 1b as this is accepted as the most accurate measure when evolutionary distances between lineages are small . Interestingly, lineages within MT-26 and MT-18, both from Turkey, were associated with different grass species. Within the globally distributed MT-1 lineage, one distinct lineage was collected in Turkey, and the second one comprised mites from multiple countries including Poland, USA, Brazil, and Germany. Analyses also identified four distinct lineages within the specialist quackgrass-associated MT-6 lineage in Poland.
Diversity within WCM lineages and tests for demographic history
The majority of the WCM lineages had high Cox1 haplotype diversity (Hd ≥ 0.7), with exception of the MT-13 lineage which had Hd = 0.5256. Lineage MT-29 had a haplotype diversity of 0.000; however, this estimate is likely compromised by the small sample size (N = 2). Nucleotide Cox1 diversity (л) ranged from 0 to 0.0696. Haplotype and nucleotide D2 diversity was in most cases lower when compared to that of Cox1 (Additional file 3: Table S3).
Tajima’s D and Fu’s Fs analyses were negative and significant for the MT-3, MT-4, and MT-8 lineages, and Fu’s Fs was also significantly negative for the MT-5 lineage (Table 1), indicating deviation from neutrality. Mismatch distributions for the WCM lineages MT-1, MT-2, MT-3, MT-4, MT-5, and MT-8 were unimodal with a significant sudden expansion model fit and non-significant raggedness test values (Table 1 and Additional file 4: Figure S4). Whereas, the sudden expansion model was rejected for the MT-7 and MT-9 lineages (Table 1 and Additional file 4: Figure S4).
Genetic structure in multi-host WCM lineages
Haplotype networks for multi-host WCM lineages did not reveal clear host associated structure (Additional file 5: Figures S5-S9). Generalist WCM lineages MT-1, MT-3, and MT-8 had complex networks, showing star-like patterns. Common ancestral haplotypes of MT-1, MT-3 and MT-8 were polyphagous, being present in many hosts (Additional file 5: Figures S5, S7, S9). The majority of MT-1 (65.3%) and MT-8 (75.6%) haplotypes were associated with cereals, whereas the majority of MT-3 haplotypes (76.4%) were quackgrass-associated. Also, the most common ancestral haplotype of lineage MT-3 (h1, 28.5% frequency) was mainly associated with quackgrass (frequency 64.8%) (Additional file 5: Figure S7). The MT-2 network (Additional file 5: Figure S6) revealed a well resolved haplotype structure, lacking a common ancestral haplotype, and displaying two highly divergent haplotypes separated by more than 10 mutations from the main haplotype cluster. In the MT-4 lineage haplotype network (Additional file 5: Figure S8), the most common central haplotype (h2) was associated with five hosts.
Analysis of molecular variance partitioned by host (Table 2) revealed non-significant subdivision for WCM lineages MT-1 (− 5.01, P = 0.67644 + − 0.01490), MT-3 (1.69, P = 0.14761 + − 0.01375), and MT-4 (7.70, P = 0.29228 + − 0.01517). Variation among different host-associated groups was moderate but significant in MT-8 (6.90, P = 0.00587 + − 0.002600), and significant and high in MT-2 (44.13, P = 0.04301 + − 0.00576), however in the latter one results were based on small sample size.
Most pairwise FST values between host-associated groups within the generalist MT-1, MT-3, MT-8 lineages did not differ significantly from zero, indicating gene flow and shared genetic diversity between lineages. In contrast, FST between triticale- and maize-associated MT-8 lineages was significantly different from zero (FST = 0.53; P < 0.01) indicating an absence of gene flow (Additional file 3: Table S4).
Haplotype diversity of generalist lineages and host-associated groups within these lineages was high in the majority of cases (Additional file 3: Table S5), and the exceptions with lower Hd values refer to groups with very low Ns.
Comparing the nucleotide diversity of host-associated populations with the whole ancestral lineage, values were similar only for MT-2 and MT-8 lineages, with exception of the MT-8 oats-associated population which had lower nucleotide diversity. For MT-1 lineage only the wheat-associated population had comparable nucleotide diversity to that observed for the whole MT-1 lineage. Other cereal-associated populations had lower nucleotide diversity when compared to the whole MT-1 lineage, whereas, on the contrary, quackgrass-associated population had higher nucleotide diversity than MT-1 lineage. Within MT-3 lineage, barley- and smooth brome-associated populations had lower nucleotide diversity, while others (wheat-, triticale- and quackgrass-associated) had comparable values to that observed for the whole MT-3 lineage. Wheat-associated MT-4 population had much lower nucleotide diversity when compared to the values observed for the whole lineage and the two wild grass-associated populations (quackgrass and smooth brome) (Additional file 3: Table S5).
Results from this study provide new and important insights into the extent and timing of lineage diversification, and the diversity of host usage in the economically important mite taxon, which is recognized as one of the most prolific and damaging pests of cereal crops around the world . Being a complex of genetic lineages with diverse patterns of host usage, it is an ideal model system for studying the mechanisms underlying the diversification processes.
WCM complex differentiation
The phylogenetic reconstructions performed in this study confirmed deep lineage diversification within the wheat curl mite, indicating a complex of lineages more diverse than previously assumed [36,37,38]. We identified 29 divergent genetic lineages, with their mitochondrial Cox1 locus differing by 6.7 to 28.2% (uncorrected p-distance), which corresponds or exceeds levels of interspecific variation observed in other animal taxa [27, 73,74,75]. Moreover, within eight of these lineages additional genetic sub-structuring was detected, suggesting that the WCM could be a complex of more than 40 distinct lineages. It is likely that lineage diversity remains grossly underestimated, there is a need for more extensive global investigations of WCM species complex. This presents a significant challenge for biodiversity assessments and cereal crop protection, as distinct lineages are expected to be biologically different and vary in their ability to vector plant pathogens . Therefore, future studies and control measures will need to account for biological and ecological differences between WCM lineages.
Although the criticism toward mitochondrial COI as a barcode for identifications and as a sequence for measuring divergence of closely related species exists (e.g. [76, 77]) COI reference data availability enables distance comparisons between case studies. Some of the advantages of employing COI for our study are that: (i) it is comparable to data from other eriophyoid mite species and distinguishes them well (e.g. [78, 79]); (ii) it is relatively easy to obtain from even scarce DNA amount extracted out of as little as one specimen ; (iii) its resolution is high enough to work on genus, species and population level ; (iv) there are known closely related Aceria species that do not exhibit such a high variation in COI as WCM, i.e. A. tulipae , and this supports notion that the variation seen in WCM is related to their overall genetic divergence, and that WCM is not a single evolutionary unit.
We identified some discrepancies in lineage delineation between the mitochondrial and nuclear DNA datasets, a phenomenon frequently reported in many animal taxa [82,83,84,85,86]. The uniparental inheritance and smaller effective population size of mitochondrial DNA suggests that lineage sorting will occur faster in mitochondrial DNA than in nuclear DNA, rendering the locus more reliable for characterizing recent lineage divergences .
Diversity in host use in relation to phylogenetic structuring
Our research has characterized WCM as a complex of at least seven specialist, six generalist, and two semi-specialist taxa, with specificity levels yet to be assigned for several other distinct genetic lineages. We expect that semi-specialist lineages MT-12 and MT-13 are in fact host specialists, given they occur predominantly on single grass host species, and only sporadically on other hosts. This is in concordance with previous research on host-related variation in population density of WCM lineages, which has indicated that MT-12 and MT-13 exhibit preferences toward one host species (tall oat-grass and quackgrass, respectively), although the lineages have been found on the other hosts sporadically .
Our phylogenetic reconstructions also indicate that neither host generalization nor specialization is a phylogenetically constrained trait, as the tree topology provides clear evidence of their convergent evolution in unrelated taxa. Shifting patterns of host generalization to specialization between closely-related lineages suggest that host range expansion and contraction is likely to be dynamic through time. The phenomenon has been well documented and can be attributed to environmental changes, including anthropogenic influence [10, 12]. Although it is not clear which factors have influenced diversification of WCM, it can be hypothesized that diversity in host use within WCM complex had been linked to vegetation transitions during the Miocene and Pliocene (discussed in next section), and differences in dispersal and colonization abilities between taxa [87, 88].
Related phytophagous insect species tend to feed on plants from similar taxonomic groupings [89,90,91,92]. However, this does not appear to apply to WCM, as evidence of phylogenetic structuring of mite lineages associated with host specificity or host taxonomy was not detected. While this is an uncommon phenomenon, it has been observed previously in butterfly Vanessa cardui , leaf beetles  and monogeneans parasitizing fish . These findings provide new insights into the diversity of host usage, however greater sampling of WCM from a broader range of grass hosts and across multiple continents is needed to gain a more reliable appreciation of absolute host ranges for WCM lineages.
Timing of WCM lineages diversification
WCM is inherently associated with grass hosts which are intensively farmed across multiple continents. Therefore, the global distribution of WCM has been previously attributed to the extension of agriculture around the globe , and here we hypothesized that agricultural practices are also likely to have had an influence on lineage diversity. However, our time calibrated phylogeny suggests that WCM radiations pre-date the origin and expansion of global agriculture, which originated approximately 11,000 year ago in the Near East . Instead, we have shown that mite lineages diverged between 16.5 mya and 19 mya, which overlaps with the timing of major diversification events within Poaceae, including ecological shifts from heavy forestation to savannas during the mid-Miocene (Fig. 3) [96, 97]. The timing of additional divergence also coincides with important events in the history of the grass family, including the expansion of C4 grasses between 7.5 and 5.5 mya, and the worldwide spread of grasslands between ca. 4.0 and 2.5 mya [97,98,99]. The divergence of cereal-associated MT-1 and MT-8 generalist appears to correspond to the final phase of global grasslands spread around 3 mya, long before the development of agriculture. Therefore, we can conclude that while agricultural practices may have had a significant influence on WCM distributions, these activities have had no detectable impact on WCM mite lineage diversification, and instead this appears to be closely linked with the early diversification of grass hosts.
Patterns of WCM lineages demographic expansion
Our analyses suggest that some generalist WCM lineages have undergone major historical demographic expansion, but again the timing of these events precedes the potential influence of agriculture. Generalist lineages associated with wild grasses appear to have underwent much earlier demographic expansions during the Pleistocene 160,000–290,000 years ago. In contrast, the most widespread generalists MT-1 and MT-8, which nowadays predominate on cereal hosts, appear to have undergone expansion events more recently, approximately 27,000 and 86,000 years ago, respectively. It is likely that during the time of recent expansion MT-1 and MT-8 lineages infested wild grasses, since many cereal hosts are no older than 11,000 years old . It is possible that these lineages were associated with progenitors of cereals, e.g. quackgrass-related plants, and had some intrinsic characteristics for rapid colonization and proliferation .
In contrast, the demographic history of specialists MT-7 (associated with wall-barley) and MT-9 (associated with smooth brome) appears to have been more stable with evidence of expansion events predating all generalist lineages (approximately 390,000 and 950,000 years ago, respectively). Expansions seem to be more likely in generalists rather than in specialists, probably given their enhanced dispersal capacities and abilities to capitalize on host availability on greater geographic scales [24, 101].
Host-associated structuring within WCM generalist lineages
High host-specificity is a prevalent feeding strategy in eriophyoid mites , therefore we expected to observe host-associated structuring in generalist WCM lineages, indicating incipient speciation processes. However, we found no evidence of these within generalist lineages MT-1, MT-3, and MT-4. In contrast, some host structuring was detected within the MT-8 lineage, with an apparent differentiation between mites from triticale and maize hosts, suggesting limited gene flow only between these two host-associated populations. However, this general absence of host-associated genetic structure within host generalist lineages indicates ongoing gene flow assisted by high dispersal and rapid colonization potential [29, 103], and suggests that generalist WCM lineages may be capable of extending their host ranges. These results are consistent with previous observations that MT-1 and MT-8 lineages are the most pestiferous within the WCM complex . This has significant implications for future pest control in cereal crops, as over-summering grass hosts provide green bridge refuges for mite persistence and re-infestation of emerging crops.
Our research focuses on patterns of genetic diversity and host association in obligatory phytophagous wheat curl mite (WCM), which is a major agricultural pest in cropping regions of the world. We revealed that lineage diversity within the taxon is significantly greater than was assumed, and consists of at least 29 genetically divergent lineages with distinctive host-use patterns, ranges from host generalists to specialists. The timing of WCM diversification events pre-dates agriculture, instead appears to be associated with timing of global diversification events within Poaceae during the mid-Miocene and the Pliocene. These novel findings provide insight into co-evolutionary processes among herbivores and their hosts, and revoke hypothesis that WCM lineage diversity is primarily linked to human activity, similarly like it has been shown for whitefly Bemisia tabaci .
Furthermore, our results shed new light on the mechanisms associated with specialization and generalization, each appearing to have evolved convergently in unrelated WCM lineages on numerous occasions. Such a pattern supports the notion that host use strategy can be highly dynamic in herbivores [11, 105], although host specificity is also known to be phylogenetically constrained in phytophagous insects  and in animal parasites . The demographic histories of generalist mite lineages also appear to be dynamic, with evidence suggesting that many have undergone major demographic expansion events, with ongoing gene flow among populations from different hosts. This is consistent with our knowledge of generalist lineages being the most prolific and damaging in cereal cropping regions of the world. Overall, we demonstrate that WCM diversity is associated with diversity in host use, which is likely to be affected by a complex of factors, including invertebrate life history traits influencing dispersal and colonization abilities, and processes such as secondary contact, host plant diversification, and environmental changes.
Akaike Information Criteria
Analysis of molecular variation
Bayesian implementation of PTP
Cytochrome c oxidase subunit I
Deegres of freedom
Effective sample sizes
- FST :
General mixed Yule coalescent
- h :
Number of haplotypes
- H d :
Highest posterior density
- HRag :
Maximum clade credibility
Markov Chain Monte Carlo
- p :
Number of polymorphic sites
Poisson tree processes model for species delimitation
Observed and expected mismatch
Recent common ancestor
Wheat curl mite, Aceria tosichella
- л :
Pennak R. Comparative ecology of the interstitial fauna of fresh-water and marine beaches. Ann Biol. 1951;27:217–480.
Cook SA, Johnson MP. Adaptation to heterogeneous environments. I. Variation in heterophylly in Ranunculus flammula L. Evolution. 1967;22:496–516.
Futuyma DJ, Moreno G. The evolution of ecological specialization. Annu Rev Ecol Evol Syst. 1988;19:207–33.
Kawecki TJ. Accumulation of deleterious mutations and the evolutionary cost of being a generalist. Am Nat. 1994;144:833–8.
Kassen R. The experimental evolution of specialists, generalists, and the maintenance of diversity. J Evol Biol. 2002;15:173–90.
Feder JL, Chilcote CA, Bush GL. Genetic differentiation between sympatric host races of the apple maggot fly Rhagoletis pomonella. Nature. 1988;336:61–4.
Drés M, Mallet J. Host races in plant-feeding insects and their importance in sympatric speciation. Philos Trans R Soc Lond Ser B Biol Sci. 2002;357:471–92.
Coyne JA, Orr HA. Speciation. Sunderland, MA: Sinauer Associates; 2004.
Van der Mescht L, Matthee S, Matthee CA. Comparative phylogeography between two generalist flea species reveal a complex interaction between parasite life history and host vicariance: parasite-host association matters. BMC Evol Biol. 2015;15:105.
Nosil P. Transition rates between specialization and generalization in phytophagous insects. Evolution. 2002;56:1701–6.
Janz N, Nyblom K, Nylin S. Evolutionary dynamics of host-plant specialization: a case study of the tribe Nymphalini. Evolution. 2001;55:783–96.
Janz N, Nylin S. The oscillation hypothesis of host plant range and speciation. In: Tilmon KJ, editor. Specialization, speciation, and radiation: the evolutionary biology of herbivorous insects. Berkeley, CA: University of California Press; 2008. p. 203–15.
Singer MC, Thomas CD, Parmesan C. Rapid human-induced evolution of insect–host associations. Nature. 1993;366:681–3.
Singer MC, Wee B, Hawkins S, Butcher M. Rapid natural and anthropogenic diet evolution: three examples from checkerspot butterflies. In: Tilmon KJ, editor. The evolutionary ecology of herbivorous insects: speciation, specialization and radiation, vol. 2008. Berkeley, CA: University of California Press; 2008. p. 311–24.
Bush GL. Sympatric host race formation and speciation in frugivorous flies of the genus Rhagoletis (Diptera, Tephritidae). Evolution. 1969;23:237–51.
Stireman JO, Nason JD, Heard SB. Host-associated genetic differentiation in phytophagous insects: general phenomenon or isolated exceptions? Evidence from a goldenrod-insect community. Evolution. 2005;59:2573–87.
Magalhães S, Forbes MR, Skoracka A, Osakabe M, Chevillon C, McCoy CD. Host race formation in the Acari. Exp Appl Acarol. 2007;42:225–38.
Audusseau H, de la Paz C-MM, Janz N, Nylin S. Why stay in a bad relationship? The effect of local host phenology on a generalist butterfly feeding on a low-ranked host. BMC Evol Biol. 2016;16(1):144.
Augustyn WJ, Anderson B, van der Merwe JF, Ellis AG. Spatial turnover in host-plant availability drives host-associated divergence in a south African leafhopper (Cephalelus uncinatus). BMC Evol Biol. 2017;17(1):72.
Hall AAG, Steinbauer MJ, Taylor GS, Johnson SN, Cook JM, Riegler M. Unravelling mummies: cryptic diversity, host specificity, trophic and coevolutionary interactions in psyllid – parasitoid food webs. BMC Evol Biol. 2017;17:127.
Powell THQ, Forbes AA, Hood GR, Feder JL. Ecological adaptation and reproductive isolation in sympatry: genetic and phenotypic evidence for native host races of Rhagoletis pomonella. Mol Ecol. 2014;23:688–704.
Groen SC, Whiteman NK. Using Drosophila to study the evolution of herbivory and diet specialization. Curr Opin Insect Sci. 2016;14:66–72.
Vanhove MPM, Hablützel PI, Pariselle A, Šimková A, Huyse T, Raeymaekers JAM. Cichlids: a host of opportunities for evolutionary parasitology. Trends Parasitol. 2016;32:820–32.
Barrett LG, Thrall PH, Burdon JJ, Linde CC. Life history determines genetic structure and evolutionary potential of host–parasite interactions. Trends Ecol Evol. 2008;23:678–85.
Johnson KP, Williams BL, Drown DM, Adams RJ, Clayton DH. The population genetics of host specificity: genetic differentiation in dove lice (Insecta: Phthiraptera). Mol Ecol. 2002;11:25–38.
Archie EA, Ezenwa VO. Population genetic structure and history of a generalist parasite infecting multiple sympatric host species. Int J Parasitol. 2011;41:89–98.
Hill MP, Hoffmann AA, Umina PA, Cheng X, Miller AD. Genetic analysis along an invasion pathway reveals endemic cryptic taxa, but a single species with little population structure in the introduced range. Divers Distrib. 2016;22:57–72.
Fraser TA, Shao R, Fountain-Jones NM, Charleston M, Martin A, Whiteley P, Holme R, Carver S, Polkinghorne A. Mitochondrial genome sequencing reveals potential origins of the scabies mite Sarcoptes scabiei infesting two iconic Australian marsupials. BMC Evol Biol. 2017;17:233.
Miller AD, Umina P, Weeks AR, Hoffmann AA. Population genetics of the wheat curl mite (Aceria tosichella Keifer) in Australia – implications for the management of wheat pathogens. Bull Entomol Res. 2012;102:199–212.
Nault LR, Styer WE. The dispersal of Aceria tulipae and three other grass-infesting Eriophyid mites in Ohio. Ann Entomol Soc Am. 1969;62:1446–55.
Navia D, Ochoa R, Welbourn C, Ferragut F. Adventive eriophyoid mites: a global review of their impact, pathways, prevention and challenges. Exp App Acarol. 2010;51:225–55.
Christian ML, Willis WG. Survival of wheat streak mosaic virus in grass hosts in Kansas from wheat harvest to fall wheat emergence. Plant Dis. 1993;77:239–42.
French R, Stenger DC. Evolution of wheat streak mosaic virus: dynamics of population growth within plants may explain limited variation. Annu Rev Phytopathol. 2003;41:199–214.
Hadi BAR, Langham MAC, Osborne L, Tilmon KJ. Wheat streak mosaic virus on wheat: biology and management. J Integr Pest Manag. 2011;1(2):1–5.
Schiffer M, Umina P, Carew M, Hoffmann A, Rodoni B, Miller A. The distribution of wheat curl mite (Aceria tosichella) lineages in Australia and their potential to transmit wheat streak mosaic virus. Ann Appl Biol. 2009;155:371–9.
Skoracka A, Kuczyński L, de Mendonca R, Dabert M, Szydło W, Knihinicki D, Truol G, Navia D. Cryptic species within the wheat curl mite Aceria tosichella (Keifer) (Acari, Eriophyoidea) revealed by mitochondrial, nuclear and morphometric data. Invertebr Syst. 2012;26:417–33.
Miller AD, Skoracka A, Navia D, de Mendonça R, Szydło W, Schultz M, Smith CM, Truol G, Hoffmann AA. Phylogenetic analyses reveal extensive cryptic speciation and host specialization in an economically important mite taxon. Mol Phylogenet Evol. 2013;66(3):928–40.
Skoracka A, Kuczyński L, Szydło W, Rector B. The wheat curl mite Aceria tosichella (Acari: Eriophyoidea) is a complex of cryptic lineages with divergent host ranges: evidence from molecular and plant bioassay data. Biol J Linnean Soc. 2013;109:165–80.
Skoracka A, Rector B, Kuczyński L, Szydło W, Hein G, French R. Global spread of wheat curl mite by its most polyphagous and pestiferous lineages. Ann App Biol. 2014;165:222–35.
Wosula EN, McMechan AJ, Oliveira-Hofman C, Wegulo SN, Hein GL. Differential transmission of two isolates of wheat streak mosaic virus by five wheat curl mite populations. Plant Dis. 2016;100:154–8.
Skoracka A, Magalhães S, Rector B, Kuczyński L. Cryptic speciation in the Acari: a function of species lifestyles or our ability to separate species? Exp Appl Acarol. 2015;67:165–82.
Skoracka A, Lewandowski M, Rector BG, Szydło W, Kuczyński L. Spatial and host-related variation in prevalence and population density of wheat curl mite (Aceria tosichella) cryptic genotypes in agricultural landscapes. PLoS One. 2017;12(1):e0169874.
Dabert J, Ehrnsberger R, Dabert M. Glaucalgales tytonis sp. n. (Analgoidea, Xolalgidae) from the barn owl Tyto alba (Strigiformes, Tytonidea): compiling morphology with DNA barcode data for the taxon description in mites (Acari). Zootaxa. 2008;1719:41–52.
Monfreda R, Lekveishvili M, Petanovič R, Amrine JW Jr. Collection and detection of eriophyoid mites. Exp Appl Acarol. 2010;51:273–82.
Dabert M, Witaliński W, Kaźmierski A, Olszanowski Z, Dabert J. Molecular phylogeny of acariform mites (Acari, Arachnida): strong conflict between phylogenetic signal and long-branch attraction artifacts. Mol Phylogenet Evol. 2010;56:222–41.
Sonnenberg R, Nolte AW, Tautz D. An evaluation of LSU rDNA D1-D2 sequences for their use in species identification. Front Zool. 2007;4:6.
Mironov SV, Dabert J, Dabert M. A new feather mite species of the genus Proctophyllodes Robin, 1877 (Astigmata: Proctophyllodidae) from the long-tailed tit Aegithalos caudatus (Passeriformes: Aegithalidae) -morphological description with DNA barcode data. Zootaxa. 2012;3253:54–61.
Szydło W, Hein G, Denizhan E, Skoracka A. Exceptionally high levels of genetic diversity in wheat curl mite (Acari: Eriophyidae) populations from Turkey. J Econ Entomol. 2015;108:2030–9.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Lewandowski M, Skoracka A, Szydło W, Kozak M, Druciarek T, Griffiths D. Genetic and morphological diversity of Trisetacus species (Eriophyoidea: Phytoptidae) associated with coniferous trees in Poland: phylogeny, barcoding, host and habitat specialization. Exp Appl Acarol. 2014;63:497–520.
Villesen P. FaBox: an online toolbox for fasta sequences. Mol Ecol Notes. 2007;7:965–8.
Swofford DL. PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. Sunderland, MA: Sinauer Associates; 2003.
Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu CH, Xie D, Suchard MA, Rambaut A, Drummond AJ. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput Biol. 2014;10(4):e1003537.
Tavaré S. Some probabilistic and statistical problems in the analysis of DNA sequences. Lectures Math Life Sci. 1986;17:57–86.
Akaike H. Likelihood of a model and information criteria. J Econom. 1981;16:3–14.
Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.
Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214. https://doi.org/10.1186/1471-2148-7-214.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
Felsenstein J. Phylogenies and the comparative method. Am Nat. 1985;125:1–15.
Brower AV. Rapid morphological radiation and convergence among races of the butterfly Heliconius erato inferred from patterns of mitochondrial DNA evolution. Proc Natl Acad Sci U S A. 1994;91:6491–5.
Papadopoulou A, Anastasiou I, Vogler AP. Revisiting the insect mitochondrial molecular clock: the mid-Aegean trench calibration. Mol Biol Evol. 2010;27:1659–72.
Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler AP. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 2006;55:595–609.
Fujisawa T, Barraclough TG. Delimiting species using single-locus data and the generalized mixed Yule coalescent (GMYC) approach: a revised method and evaluation on simulated datasets. Syst Biol. 2013;62:707–24.
R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, 2008. http://www.R-project.org. Accessed 15 Jan 2017.
Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29:2869–76.
Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and windows. Mol Ecol Resour. 2010;10:564–7.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.
Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.
Rogers AR, Harpending H. Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992;9:552–69.
Bandelt H, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.
Yang Z. On the best evolutionary rate for phylogenetic analysis. Syst Biol. 1998;47:125–33.
Navia D, de Mendoça R, Skoracka A, Szydło W, Knihinicki D, Hein GL, da Silva Pereira PR, Truol G, Lau D. Wheat curl mite, Aceria tosichella, and transmitted viruses: an expanding pest complex affecting cereal crops. Exp Appl Acarol. 2013;59:95–143.
Hebert PDN, Ratnasingham S, de Waard JR. Barcoding animal life: cytochrome c oxidase subunit I divergences among closely related species. Proc R Soc Lond B Biol Sci. 2003;270:S96–9.
Lee Y, Lee W, Kanturski M, Foottit RG, Akimoto S-I, Lee S. Cryptic diversity of the subfamily Calaphidinae (Hemiptera: Aphididae) revealed by comprehensive DNA barcoding. PLoS One. 2017;12(4):e0176582.
Liu Y, Fend SV, Martinsson S, Erséus C. Extensive cryptic diversity in the cosmopolitan sludge worm Limnodrilus hoffmeisteri (Clitellata, Naididae). Org Divers Evol. 2017;17:477–95.
Moritz C, Cicero C. DNA barcoding: promise and pitfalls. PLoS Biol. 2004;2(10):e354.
Rubinoff D, Cameron S, Will K. A genomic perspective on the shortcomings of mitochondrial DNA for “barcoding” identification. J Hered. 2006;97:581–94.
Guo JF, Li HS, Wang B, Xue XF, Hong XY. DNA barcoding reveals the protogyne and deutogyne of Tegolophus celtis sp. nov. (Acari: Eriophyidae). Exp Appl Acarol. 2015;67:393–410.
Cvrković T, Chetverikov P, Vidović B, Petanović R. Cryptic speciation within Phytoptus avellanae s.L. (Eriophyoidea: Phytoptidae) revealed by molecular data and observations on molting Tegonotus-like nymphs. Exp Appl Acarol. 2016;68:83–96.
Navajas M, Navia D. DNA-based methods for eriophyoid mite studies: review, critical aspects, prospects and challenges. Exp Appl Acarol. 2010;51:257–71.
Kiedrowicz A, Rector BG, Lommen S, Kuczyński L, Szydło W, Skoracka A. Population growth rate of dry bulb mite, Aceria tulipae (Acariformes: Eriophyidae), on agriculturally important plants and implications for its taxonomic status. Exp Appl Acarol. 2017;73(1):1–10.
Funk DJ, Omland KE. Species level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu Rev Ecol Evol Syst. 2003;34:397–423.
Hudson RR, Turelli M. Stochasticity overrules the “three-times rule”: genetic drift, genetic draft, and coalescence times for nuclear loci versus mitochondrial DNA. Evolution. 2003;57:182–90.
Zink RM, Barrowclough GF. Mitochondrial DNA under siege in avian phylogeography. Mol Ecol. 2008;17:2107–21.
McKay BD, Zink RM. The causes of mitochondrial DNA gene tree papraphyly in birds. Mol Phylogenet Evol. 2010;54:647–50.
Thielsch A, Knell A, Mohammadyari A, Petrusek A, Schwenk K. Divergent clades or cryptic species? Mito-nuclear discordance in a Daphnia species complex. BMC Evol Biol. 2017;17:227.
Sabelis MW, Bruin J. Evolutionary ecology: life history patterns, food plant choice and dispersal. In: Lindquist EE, Sabelis MW, Bruin J, editors. Eriophyoid mites: their biology, natural enemies and control. Amsterdam: Elsevier Science BV; 1996. p. 329–66.
Dahirel M, Olivier E, Guiller A, Martin M-C, Madec L, Ansart A. Movement propensity and ability correlate with ecological specialization in European land snails: comparative analysis of a dispersal syndrome. J Anim Ecol. 2015;84:228–38.
Futuyma DJ. Some current approaches to the evolution of plant-herbivore interactions. Plant Species Biol. 2000;15:1–9.
Engelstädter J, Hurst GDD. The dynamics of parasite incidence across host species. Evol Ecol. 2006;20:603–16.
Poulin R. Decay of similarity with host phylogenetic distance in parasite faunas. Parasitology. 2010;137:733–41.
Cooper N, Griffin R, Franz M, Omotayo M, Nunn CL, Fryxell J. Phylogenetic host specificity and understanding parasite sharing in primates. Ecol Lett. 2012;15:1370–7.
Celorio-Mancera MP, Wheat CW, Huss M, Vezzi F, Neethiraj R, Reimegård J, Nylin S, Janz N. Evolutionary history of host use, rather than plant phylogeny, determines gene expression in a generalist butterfly. BMC Evol Biol 2016;16:59.
Mendlová M, Šimková A. Evolution of host specificity in monogeneans parasitizing African cichlid fish. Parasit Vectors. 2014;7:69.
Abbo S, Lev-Yadun S, Gopher A. Plant domestication and crop evolution in the near east: on events and processes. CRC Crit Rev Plant Sci. 2012;31:241–57.
Jacobs BF, Kingston JD, Jacobs LL. The origin of grass-dominated ecosystems. Ann Mo Bot Gard. 1999;86:590–643.
Kellogg E. Evolutionary history of the grasses. Plant Physiol. 2001;125:1198–205.
Edwards EJ, Osborne CP, Strömberg CAE, Smith SA. C 4 grasses consortium. The origins of C4 grasslands: integrating evolutionary and ecosystem science. Science. 2010;328:587–91.
Edwards EJ, Smith SA. Phylogenetic analyses reveal the shady history of C4 grasses. Proc Natl Acad Sci U S A. 2010;107:2532–7.
Preece C, Livarda A, Christin P-A, Wallace M, Martin G, Charles M, Jones G, Rees M, Osborne CP. How did the domestication of Fertile Crescent grain crops increase their yields? Funct Ecol. 2017;31:387–97.
Levins R. Theory of fitness in a heterogeneous environment. 1. Fitness set and adaptive function. Am Nat. 1962;96:361–73.
Skoracka A, Smith L, Oldfield G, Cristofarro M, Amrine JW. Host specificity and specialization in eriophyoid mites and their importance for the use of eriophyoid mites as biocontrol agents of weed. Exp Appl Acarol. 2010;51:93–113.
Kiedrowicz A, Kuczyński L, Lewandowski M, Proctor H, Skoracka A. Behavioural responses to potential dispersal cues in two economically important species of cereal-feeding eriophyid mites. Sci Rep. 2017;7:3890.
Boykin LM, Bell CD, Evans G, Small I, De Barro PJ. Is agriculture driving the diversification of the Bemisia tabaci species complex (Hemiptera: Sternorrhyncha: Aleyrodidae)?: dating, diversification and biogeographic evidence revealed. BMC Evol Biol. 2013;13:228.
Kohyama TI, Matsumoto K, Katakura H. Deep phylogeographical structure and parallel host range evolution in the leaf beetle Agelasa nigriceps. Mol Ecol. 2014;23:421–34.
Ehrlich PR, Raven PH. Butterflies and plants: a study in coevolution. Evolution. 1964;18:586–608.
Krasnov BR, Mouillot D, Shenbrot GI, Khokhlova IS, Poulin R. Geographical variation in host specificity of fleas (Siphonaptera) parasitic on small mammals: the influence of phylogeny and local environmental conditions. Ecography. 2004;27:787–97.
The authors thank to Jarosław Raubic for help in laboratory analyses.
The study was funded by the National Science Centre in Poland, grants no. 2011/03/B/NZ8/00129 given to Skoracka A and 2011/01/N/NZ8/04540 given to Szydło W. Lopes LF was funded by the FCT fellowship SFRH/BPD/84478/2012. Alves MJ and Lopes LF were financed by Portuguese National Funds, through FCT – Fundação para a Ciência e a Tecnologia, within the project UID/BIA/00329/2013. The work was supported by the KNOW RNA Research Centre in Poznań (No. 01/KNOW2/2014) (Funding for open access charge: Ministry of Science and Higher Education of the Republic of Poland, from the quality promoting subsidy, under the Leading National Research Centre (KNOW) program for the years 2014–2019).
The funding bodies provided funds for data collection (1 and 2), analyses, data interpretation (1–4), and writing the manuscript (all 5).
Availability of data and materials
The datasets supporting the conclusions of this article are included in this published article and its supplementary information files. The sequence data used in this study have been submitted to the GenBank databases under accession numbers: MG076809- MG076875, MG193928- MG194416, MG208676- MG208850. The frozen DNA isolates and DNA samples (mites stored in ATL buffer) of field samples collected in Poland are available in the Faculty of Biology, Adam Mickiewicz University, Poznań, Poland (unless not exhausted for repeated or other analyses).
Ethics approval and consent to participate
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.
Sequence and sampling information. (DOCX 99 kb)
Phylogenetic trees for wheat curl mite sequence datasets. (DOCX 776 kb)
Statistics for genetic diversity of WCM genetic lineages. (DOCX 738 kb)
Mismatch distributions of WCM genetic lineages. (DOCX 747 kb)
Haplotype networks for the mtDNA Cox1 wheat curl mite (WCM). (DOCX 950 kb)
About this article
Cite this article
Skoracka, A., Lopes, L.F., Alves, M.J. et al. Genetics of lineage diversification and the evolution of host usage in the economically important wheat curl mite, Aceria tosichella Keifer, 1969. BMC Evol Biol 18, 122 (2018). https://doi.org/10.1186/s12862-018-1234-x
- Aceria tosichella
- Demographic history
- Genetic diversity
- Lineage diversification
- Species delimitation