Skip to main content
  • Research article
  • Open access
  • Published:

Genetics of lineage diversification and the evolution of host usage in the economically important wheat curl mite, Aceria tosichella Keifer, 1969

Abstract

Background

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.

Results

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.

Conclusions

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.

Background

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 [11]. Populations of generalist species might restrict their host range as a result of diminishing host availability or costs of host adaptation [5]. 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 [21]. 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 [25] and nematodes in livestock [26] 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 [24]. 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 [29]. 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.

Methods

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. [43]. The exoskeletons of the digested mites were preserved in 70% ethyl alcohol, and later mounted on slides according to Monfreda et al. [44] 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 [46] and 28SR0990 [47]. Reactions steps and product handling followed protocols described by Szydło et al. [48]. 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 [48]. Trace files were checked and edited using MEGA 6 [49].

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 [50]. 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 [51] 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 [52].

Phylogenetic analyses and divergence time estimation

Phylogenetic reconstructions were performed using Bayesian Inference (BI) methods implemented in BEAST 2.3.0 [53]. General Time Reversible model [54] 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) [55] implemented in JMODELTEST v.0.1.1 [56]. 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 [57] 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 [58], 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 [58]. 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 [59] 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 [58], 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 [64] 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 [65]. 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 [49].

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.3.5.2.2 [66] 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 [67] and Fu’s Fs [68]. 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 [69]. 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 [70] 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.

Results

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).

Fig. 1
figure 1

a Bayesian inference (BI) tree constructed using the GTR + G model for the cytochrome c oxidase subunit 1 (Cox1) sequences of the wheat curl mite (WCM) species complex and outgroup species. WCM individuals on the tree are colored according to their grass host species; color boxes correspond to cultivated grasses (cereals) and color labels correspond to non-cultivated grasses. Sequence labels correspond to data in Additional file 1: Table S1and include information about the host species and geographic locality. Vertical bars indicate particular WCM lineages, and these of strict specialists (s) are colored according to their specific host plants. Lineages associated with several hosts are designated as generalists (g), and lineages found on two hosts with higher prevalence on one of them are designated as semi-specialists (ss). Numbers above branches are Bayesian posterior probabilities; only values > 0.6 are shown. b Results of combined Poisson Tree Processes model for species delimitation (PTP) and a Bayesian implementation of PTP (bPTP) identifying additional putative species groups within MT-1, MT-3, MT-6, MT-9, MT-12, MT-17 MT-18, and MT-26 lineages. Blue lines indicate the unique species groups and red clades indicate species groups contained more than one haplotypes. Numbers above the branches are Bayesian posterior probabilities

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 [71]. 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).

Fig. 2
figure 2

Bayesian inference (BI) tree constructed using the GTR + G model for the haplotypes of 28S rDNA D2 region of the wheat curl mite (WCM) species complex and outgroup species. Numbers above branches are Bayesian posterior probabilities; only values > 0.6 are shown

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.

Fig. 3
figure 3

Time calibrated mitochondrial Cox1 gene tree showing mean values in million years (mya) above branches and the 95% HPD bars for divergence time on nodes. Color boxes indicate events that may explain the diversification and divergence of WCM lineages

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 [65]. 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).

Table 1 Demographic analyses and recent expansion time for WCM lineages (for those with N ≥ 15). The estimates of neutrality tests for each lineage were determined by Tajima’s D and Fu’s Fs statistics. The validity of sudden expansion model was tested by mismatch distributions: SSD, the sum of square deviations; HRag, the raggedness index

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.

Table 2 Analyses of molecular variance (AMOVA) for generalist WCM lineages with host as a grouping variable

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).

Discussion

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 [72]. 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 [37]. 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 [80]; (iii) its resolution is high enough to work on genus, species and population level [80]; (iv) there are known closely related Aceria species that do not exhibit such a high variation in COI as WCM, i.e. A. tulipae [81], 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 [84].

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 [42].

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 [93], leaf beetles [89] and monogeneans parasitizing fish [94]. 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 [39], 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 [95]. 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 [100]. 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 [29].

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 [102], 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 [39]. 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.

Conclusions

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 [104].

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 [106] and in animal parasites [107]. 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.

Abbreviations

AIC:

Akaike Information Criteria

AMOVA:

Analysis of molecular variation

BI:

Bayesian inference

bPTP:

Bayesian implementation of PTP

Cox1:

Cytochrome c oxidase subunit I

d.f.:

Deegres of freedom

ESSs:

Effective sample sizes

FST :

Fixation index

GMYC:

General mixed Yule coalescent

h :

Number of haplotypes

H d :

Haplotype diversity

HPD:

Highest posterior density

HRag :

Raggedness index

MCC:

Maximum clade credibility

MCMC:

Markov Chain Monte Carlo

ML:

Maximum likelihood

p :

Number of polymorphic sites

PTP:

Poisson tree processes model for species delimitation

SSD:

Observed and expected mismatch

TBR:

Tree bisection-reconnection

tMRCAs:

Recent common ancestor

WCM:

Wheat curl mite, Aceria tosichella

л :

Nucleotide diversity

References

  1. Pennak R. Comparative ecology of the interstitial fauna of fresh-water and marine beaches. Ann Biol. 1951;27:217–480.

    Google Scholar 

  2. Cook SA, Johnson MP. Adaptation to heterogeneous environments. I. Variation in heterophylly in Ranunculus flammula L. Evolution. 1967;22:496–516.

    Google Scholar 

  3. Futuyma DJ, Moreno G. The evolution of ecological specialization. Annu Rev Ecol Evol Syst. 1988;19:207–33.

    Article  Google Scholar 

  4. Kawecki TJ. Accumulation of deleterious mutations and the evolutionary cost of being a generalist. Am Nat. 1994;144:833–8.

    Article  Google Scholar 

  5. Kassen R. The experimental evolution of specialists, generalists, and the maintenance of diversity. J Evol Biol. 2002;15:173–90.

    Article  Google Scholar 

  6. Feder JL, Chilcote CA, Bush GL. Genetic differentiation between sympatric host races of the apple maggot fly Rhagoletis pomonella. Nature. 1988;336:61–4.

    Article  Google Scholar 

  7. 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.

    Article  Google Scholar 

  8. Coyne JA, Orr HA. Speciation. Sunderland, MA: Sinauer Associates; 2004.

    Google Scholar 

  9. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Nosil P. Transition rates between specialization and generalization in phytophagous insects. Evolution. 2002;56:1701–6.

    Article  PubMed  CAS  Google Scholar 

  11. Janz N, Nyblom K, Nylin S. Evolutionary dynamics of host-plant specialization: a case study of the tribe Nymphalini. Evolution. 2001;55:783–96.

    Article  PubMed  CAS  Google Scholar 

  12. 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.

    Google Scholar 

  13. Singer MC, Thomas CD, Parmesan C. Rapid human-induced evolution of insect–host associations. Nature. 1993;366:681–3.

    Article  Google Scholar 

  14. 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.

    Google Scholar 

  15. Bush GL. Sympatric host race formation and speciation in frugivorous flies of the genus Rhagoletis (Diptera, Tephritidae). Evolution. 1969;23:237–51.

    Article  PubMed  Google Scholar 

  16. 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.

    Article  PubMed  CAS  Google Scholar 

  17. 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.

    Article  PubMed  Google Scholar 

  18. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. 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.

    Article  PubMed  Google Scholar 

  22. Groen SC, Whiteman NK. Using Drosophila to study the evolution of herbivory and diet specialization. Curr Opin Insect Sci. 2016;14:66–72.

    Article  PubMed  Google Scholar 

  23. 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.

    Article  PubMed  Google Scholar 

  24. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 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.

    Article  PubMed  CAS  Google Scholar 

  26. 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.

    Article  PubMed  Google Scholar 

  27. 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.

    Article  Google Scholar 

  28. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 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.

    Article  PubMed  CAS  Google Scholar 

  30. 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.

    Article  Google Scholar 

  31. 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.

    Article  Google Scholar 

  32. 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.

    Article  Google Scholar 

  33. 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.

    Article  PubMed  CAS  Google Scholar 

  34. 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.

    Article  Google Scholar 

  35. 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.

    Article  Google Scholar 

  36. 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.

    Article  Google Scholar 

  37. 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.

    Article  PubMed  Google Scholar 

  38. 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.

    Article  Google Scholar 

  39. 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.

    Article  Google Scholar 

  40. 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.

    Article  CAS  PubMed  Google Scholar 

  41. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  42. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  43. 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.

    Google Scholar 

  44. Monfreda R, Lekveishvili M, Petanovič R, Amrine JW Jr. Collection and detection of eriophyoid mites. Exp Appl Acarol. 2010;51:273–82.

    Article  PubMed  CAS  Google Scholar 

  45. 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.

    Article  PubMed  Google Scholar 

  46. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. 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.

    Google Scholar 

  48. 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.

    Article  PubMed  Google Scholar 

  49. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Villesen P. FaBox: an online toolbox for fasta sequences. Mol Ecol Notes. 2007;7:965–8.

    Article  CAS  Google Scholar 

  52. Swofford DL. PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. Sunderland, MA: Sinauer Associates; 2003.

    Google Scholar 

  53. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  54. Tavaré S. Some probabilistic and statistical problems in the analysis of DNA sequences. Lectures Math Life Sci. 1986;17:57–86.

    Google Scholar 

  55. Akaike H. Likelihood of a model and information criteria. J Econom. 1981;16:3–14.

    Article  Google Scholar 

  56. Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.

    Article  PubMed  CAS  Google Scholar 

  57. Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2148-7-214.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. Felsenstein J. Phylogenies and the comparative method. Am Nat. 1985;125:1–15.

    Article  Google Scholar 

  60. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. Papadopoulou A, Anastasiou I, Vogler AP. Revisiting the insect mitochondrial molecular clock: the mid-Aegean trench calibration. Mol Biol Evol. 2010;27:1659–72.

    Article  PubMed  CAS  Google Scholar 

  62. 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.

    Article  PubMed  Google Scholar 

  63. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  64. 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.

  65. Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29:2869–76.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. 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.

    Article  PubMed  Google Scholar 

  67. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.

    PubMed  PubMed Central  CAS  Google Scholar 

  68. Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.

    PubMed  PubMed Central  CAS  Google Scholar 

  69. Rogers AR, Harpending H. Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992;9:552–69.

    PubMed  CAS  Google Scholar 

  70. Bandelt H, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.

    Article  PubMed  CAS  Google Scholar 

  71. Yang Z. On the best evolutionary rate for phylogenetic analysis. Syst Biol. 1998;47:125–33.

    Article  PubMed  CAS  Google Scholar 

  72. 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.

    Article  PubMed  Google Scholar 

  73. 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.

    Article  CAS  Google Scholar 

  74. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  75. 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.

    Article  Google Scholar 

  76. Moritz C, Cicero C. DNA barcoding: promise and pitfalls. PLoS Biol. 2004;2(10):e354.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  77. Rubinoff D, Cameron S, Will K. A genomic perspective on the shortcomings of mitochondrial DNA for “barcoding” identification. J Hered. 2006;97:581–94.

    Article  PubMed  CAS  Google Scholar 

  78. 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.

    Article  PubMed  CAS  Google Scholar 

  79. 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.

    Article  PubMed  Google Scholar 

  80. Navajas M, Navia D. DNA-based methods for eriophyoid mite studies: review, critical aspects, prospects and challenges. Exp Appl Acarol. 2010;51:257–71.

    Article  PubMed  CAS  Google Scholar 

  81. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  82. 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.

    Article  Google Scholar 

  83. 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.

    PubMed  Google Scholar 

  84. Zink RM, Barrowclough GF. Mitochondrial DNA under siege in avian phylogeography. Mol Ecol. 2008;17:2107–21.

    Article  PubMed  CAS  Google Scholar 

  85. McKay BD, Zink RM. The causes of mitochondrial DNA gene tree papraphyly in birds. Mol Phylogenet Evol. 2010;54:647–50.

    Article  PubMed  Google Scholar 

  86. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  87. 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.

    Chapter  Google Scholar 

  88. 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.

    Article  PubMed  Google Scholar 

  89. Futuyma DJ. Some current approaches to the evolution of plant-herbivore interactions. Plant Species Biol. 2000;15:1–9.

    Article  Google Scholar 

  90. Engelstädter J, Hurst GDD. The dynamics of parasite incidence across host species. Evol Ecol. 2006;20:603–16.

    Article  Google Scholar 

  91. Poulin R. Decay of similarity with host phylogenetic distance in parasite faunas. Parasitology. 2010;137:733–41.

    Article  PubMed  CAS  Google Scholar 

  92. 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.

    Article  PubMed  Google Scholar 

  93. 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.

    Article  PubMed Central  CAS  Google Scholar 

  94. Mendlová M, Šimková A. Evolution of host specificity in monogeneans parasitizing African cichlid fish. Parasit Vectors. 2014;7:69.

    Article  PubMed  PubMed Central  Google Scholar 

  95. 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.

    Article  Google Scholar 

  96. Jacobs BF, Kingston JD, Jacobs LL. The origin of grass-dominated ecosystems. Ann Mo Bot Gard. 1999;86:590–643.

    Article  Google Scholar 

  97. Kellogg E. Evolutionary history of the grasses. Plant Physiol. 2001;125:1198–205.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  98. 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.

    Article  PubMed  CAS  Google Scholar 

  99. Edwards EJ, Smith SA. Phylogenetic analyses reveal the shady history of C4 grasses. Proc Natl Acad Sci U S A. 2010;107:2532–7.

    Article  PubMed  PubMed Central  Google Scholar 

  100. 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.

    Article  PubMed  Google Scholar 

  101. Levins R. Theory of fitness in a heterogeneous environment. 1. Fitness set and adaptive function. Am Nat. 1962;96:361–73.

    Article  Google Scholar 

  102. 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.

    Article  PubMed  Google Scholar 

  103. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  104. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  105. 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.

    Article  PubMed  CAS  Google Scholar 

  106. Ehrlich PR, Raven PH. Butterflies and plants: a study in coevolution. Evolution. 1964;18:586–608.

    Article  Google Scholar 

  107. 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.

    Article  Google Scholar 

Download references

Acknowledgements

The authors thank to Jarosław Raubic for help in laboratory analyses.

Funding

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).

Author information

Authors and Affiliations

Authors

Contributions

AS, LK designed the study. AS, LK, ML, WS, A Majer performed the field work and collected material. ML, WS, A Majer, ER done the laboratory work and obtained sequences. LFL, MJA, A Miller, AS analyzed the data. AS, A Miller, ML drafted the manuscripts. The manuscript was revised by all of the authors. All authors read, approved and contributed to the final manuscript.

Corresponding author

Correspondence to Anna Skoracka.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Sequence and sampling information. (DOCX 99 kb)

Additional file 2:

Phylogenetic trees for wheat curl mite sequence datasets. (DOCX 776 kb)

Additional file 3:

Statistics for genetic diversity of WCM genetic lineages. (DOCX 738 kb)

Additional file 4:

Mismatch distributions of WCM genetic lineages. (DOCX 747 kb)

Additional file 5:

Haplotype networks for the mtDNA Cox1 wheat curl mite (WCM). (DOCX 950 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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://0-doi-org.brum.beds.ac.uk/10.1186/s12862-018-1234-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12862-018-1234-x

Keywords