- Research article
- Open Access
Phylogenetic signals in pest abundance and distribution range of spider mites
BMC Evolutionary Biology volume 19, Article number: 223 (2019)
Attributes of pest species like host range are frequently reported as being evolutionarily constrained and showing phylogenetic signal. Because these attributes in turn could influence the abundance and impact of species, phylogenetic information could be useful in predicting the likely status of pests. In this study, we used regional (China) and global datasets to investigate phylogenetic patterns in occurrence patterns and host ranges of spider mites, which constitute a pest group of many cropping systems worldwide.
We found significant phylogenetic signal in relative abundance and distribution range both at the regional and global scales. Relative abundance and range size of spider mites were positively correlated with host range, although these correlations became weaker after controlling for phylogeny.
The results suggest that pest impacts are evolutionarily constrained. Information that is easily obtainable – including the number of known hosts and phylogenetic position of the mites – could therefore be useful in predicting future pest risk of species.
The human modification of natural environments including expansion of agricultural production areas has been a primary driver of terrestrial biodiversity loss . Although hundreds of species have been documented as dramatically declining under habitat modification, some species are thriving , including agricultural pests and pathogens, which in turn have led to additional stresses on non-pest species . Understanding why some species fare poorly whereas others do well has been a key issue of concern to biologists, ecologists, agriculturists and policymakers [46, 61], and is an important consideration when assessing future risks of species extinctions as well as pest outbreaks .
Species extinction risk is often not randomly spread across phylogeny , indicating that phylogeny could be useful in predicting the fate of species [19, 52]. In risk assessments, phylogenetic information has also been used to predict which plant species are likely to be susceptible to a particular pest [23, 55], because closely-related plants tend to have similar traits (e.g. plant defensive chemicals) and host similar pests when compared to evolutionarily distant plant species . However, the host plant records for many novel pests are incomplete, and the severity of pest outbreaks may not be closely linked to this factor; for instance, information of pest host range and host phylogeny was insufficient to determine whether pests on a given host or novel region were severe or benign .
Predicting potential risks posed by a pest or pathogen requires an understanding of a range of biological and ecological characteristics for adapting to particular hosts and agricultural contexts [11, 16], as well as an assessment of the degree to which these are constrained within phylogeny [35, 66]. Species’ performance is often determined by traits that show a strong tendency to take similar values among closely related species , including host plant range  and thermal resistance . These traits in turn are likely to alter the demographic characteristics of species and link to species distribution and abundance ([25, 32]) which in turn are expected to be phylogenetically structured. Phylogenetic distance between poorly-known novel pests and well-studied pest species with known occurrence pattern may therefore help predict whether a novel pest is likely to have severe effects.
To test these ideas, we examined phylogenetic patterns of pest occurrence in a phytophagous species group, the spider mite family (Acari: Tetranychidae). This family includes more than 1300 species (around 100 of which are considered pests) that share similar morphological characters but vary in host breadth (extremely polyphagous vs. highly host-specific) and distribution range (widespread vs. narrow distributions) . These features make spider mites a useful group to investigate host- and distribution-related ecological and evolutionary questions [63, 64], although most previous work has focused on the intra-specific level rather than the comparative level. We therefore aim to test for phylogenetic signal in relative abundance and host range and link findings to pest outbreak patterns. This work can guide phytosanitary risk analysis of pests and their potential impact before pests arise in a region [25, 54].
Here we used long-term survey data from 2008 to 2017 in China and species information from a global dataset to test for phylogenetic signatures in spider mites. Our analysis addresses three main questions. First, are the host range and relative abundance of spider mites non-random within phylogeny? Second, does host range correlate with species occurrence patterns? And third, how does phylogeny influence the host range - abundance relationship?
Phylogenetic signal of pest occurrence
In our field survey, twelve spider mite taxa were found from 318 populations at 180 sites (Fig. 1, also see Additional file 1: Table S1), which included 7596 samples. Ten species of spider mites belonged to Tetranychus, one species belonged to Panonychus and one species belonged to Amphitetranychus (Fig. 1, also see Additional file 1: Table S3. The two measures of mite abundance (NOC and AF) were strongly and positively correlated (Pearson r = 0.996, p < 0.001). High correlation coefficient (Pearson r = 0.959, p < 0.001) also was found between two measures of host range (host species and host family). We therefore only considered NOC and host family in the analyses.
Ten of 12 species occurring in this survey belonged to the genus Tetranychus, and the molecular distance between these species was compared to the pattern of occurrence and host range in China. Phylogenetic trees used for phylogenetic signal tests were reconstructed with the BI and ML methods (Additional file 1: Figure S2 in ESM). The BI and ML trees of the combined three DNA fragments (COI, 18S and 28S) resulted in identical topologies for both China dataset (Additional file 1: Figure S1) and global dataset (Additional file 1: Figure S2, S3). In both cases, the topology is mostly well supported, with a bootstrap support value and posterior probability value in major nodes. The phylogenies were corresponding with prior phylogenies . No significant correlation between relative abundance and genetic distance was detected using a Mantel test (Additional file 1: Table S4, p = 0.119). However, the phylogenetic signal for species abundance (across all 10 species) [log10 (NOC)] was higher than random expectations and larger than expected under a Brownian motion model of character evolution (K = 1.032, p = 0.033) and a significant phylogenetic signal was also detected using Abouheif’s test (p = 0.013). The species distribution range size (latitudinal span) showed significant phylogenetic signal on all three measures (Table 1). For the global dataset, weak phylogenetic signal was detected in occurrence patterns (historical records number and distributed country number) with all three measures (Table 1).
For Tetranychus species in China, the relative abundance (Fig. 2a, r = 0.943, p < 0.001) and distribution range (Fig. 2b, r = 0.924, p < 0.001) of species declined significantly with increasing phylogenetic distance to focal species based on a correlation analysis. Such a pattern also existed in the global dataset (Fig. 2c, d). For the global dataset, the correlation coefficients and their significance were different among the genera (Additional file 1: Table S5). Three of four tested genera (Eotetranychus, Oligonychus and Panonychus) showed significant correlations between occurrence patterns (historical records number and distributed country number) and genetic distance to the focal species for the global dataset.
Relationships between host range and species occurrence
For species in China, the relative abundance and latitudinal range of spider mites were significantly associated with host range (Fig. 3a, b). To investigate whether phylogeny influences the host range – species occurrence relationship, we performed phylogenetically corrected correlations between host range and pest occurrence. The relationships between host range and pest occurrence tended to become weaker with lower coefficients after PIC and PGLS correction for phylogeny (Table 2). Similar patterns were also found in the global dataset (Fig. 3c, d). Although significant correlations between host range and pest occurrence remained, the strength of all relationships was reduced by phylogenic correction (Table 2).
The ancestral trait reconstruction showed different patterns of host range evolution among clades at the genus and subgenus levels (Fig. 4). This analysis suggested a monophagous origin of spider mites. The ancestral state of narrow host range seems persist in other clades within the evolutionary history of spider mites. But host range expanded rapidly in the clade Tetranychus, and the evolution of host range expansion was mostly restricted to this group. Several species in other groups also had a wide host range (e.g. Oligonychus coffeae), yet most maintained a narrow host range as for the ancestral form.
We found strong phylogenetic signal in pest abundance and distribution of Tetranychus mites at a regional scale (China), and also detected a phylogenetic signature for species occurrence when analyzing 88 spider mite species using a global dataset. Pest occurrence (relative abundance and distribution range size) declined predictably with increasing genetic distance from the most abundant pest species. These results suggest that species occurrence can be partly predicted by evolutionary relationships in the spider mite group.
Several studies have introduced phylogenetic information into pest risk assessments [22, 23]. However, such information was not sufficient to evaluate whether pest damage on a given host or in a novel region is severe or benign . Here we tested another hypothesis, namely that the phylogeny of pests themselves can be used to predict which pest species are likely to be abundant.
The evolutionary history of pests can significantly affect their capacity to adapt to new host plants or novel environments [13, 50], leading to a potential relationship between pest severity and phylogeny. This was confirmed in spider mites; there were significant phylogenetic signal in both relative abundance and distribution range size. Many species belonging to the genus Tetranychus were relatively common and also were serious pests locally (e.g. T. truncatus in China)  and with the potential to become global pests. For example, T. urticae and T. evansi have expanded their distribution and become serious pests in many regions [7, 26, 62]. However a categorical metric of pest risk decision process, (e.g., all pest within a genus are risky, and others are not) is not ideal for risk analysis of novel pests and pathogens . Because we found that the relative abundance and distribution range size declined as a function of phylogenetic distance between congeneric spider mites, we suspect that species relatedness data within genera may be useful in pest risk assessments in the absence of other empirical information.
In this study, we detected strong and positive correlations between host range and relative abundance at both the country and global scales. Species niche breadth is often considered to reflect an evolutionary trade-off between a species’ ability to exploit a wide range of resources and the effectiveness of exploitation [9, 27], resulting in a lower abundance of species exploting a broader host range [8, 65]. However, there are other theories [e.g. the hierarchical theory posed by Passy ] arguing that species with the highest maximum abundance and regional prevalence possess the broadest niches, especially under a stressful environment (e.g. short resource supply, human impact). The applicability of these hypotheses may reflect the degree of disturbance in the environment, perhaps caused by human-associated changes (e.g. agriculture and urban expanding, polluting or global warming) [4, 5]. This may help explain why generalists seem to benefit from global change more than specialists .
In this study, there was some evidence for a lower abundance of specialists. When we analyzed each host individually, there were diverse relationships (positive, negative or none) between abundance and specialization. When we considered all paired values for the two datasets (12 species in China and the global dataset), there was no significant association for the China dataset (see Additional file 1: Figure S4a). However, for the global dataset, we found negative correlations between host specialization and species occurrence on specific hosts (see Additional file 1: Figure S4b). This runs contrary to the expectation of a trade-off between species niche breadth and performance on a particular host. In insect herbivores, a global scale study showed that more diverse lineages of plants support assemblages of relatively more specialized herbivores . This suggests a lower abundance of specialized herbivores could relate to reduced plant diversity. However, positive host range –abundance relationships in spider mites may reflect the fact mites with a wider niche breadth can reproduce and persist in agricultural ecosystems on a greater range of crops and therefore build up across time.
The moderate to high phylogenetic signal in host range at both regional and global scales suggests non-random evolution of host range in spider mites. Two phylogenetically corrected correlation analyses (PIC and PGLS) showed that relationships tended to become weaker with lower coefficients after correction for phylogeny. This finding was supported by the ancestral trait reconstruction analysis (Fig. 4) showing that the evolutionary pattern of host range was different among clades. The majority of clades showed relative conservative patterns in host range, whereas host range in some clades have rapidly expanded after an early split with others. In particular, the frequency of evolutionary expansion in host range appears to increase dramatically at the Tetranychus group. Such patterns indicate evolutionary history is important to understanding species’ status in community [30, 67].
Compared to local scale, we found host range and occurrence showed lower phylogenetic signals (Table 1). The association between species occurrence and genetic distance also tend to be weaker at global scale (Fig. 2). However, the association between species occurrence and host range tend to be stronger at global scale than at local scale (Table 2). These results suggested host range maybe more relevant than phylogenetic signals on predicting pest risk at global scale.
False positive correlations can be produced across species comparisons, including scale selection , sampling effects  and statistics . The different strength and significance level of phylogenetic signal between regional and global scale suggests phylogenetic patterns could be influenced by sampling issues . In this study, only 88 of 1300 species were used for phylogenetic analysis, although the most comprehensive tree was developed based on available data . Further analysis of additional species (and particularly rare species) may provide insights into how species abundance and niches are distributed across subgeneric-level phylogenies. For tests of the host range – pest occurrence relationship, sampling may generate a positive relationship . For our data from China, we suspect that these issues are not likely to obscure patterns. The occurrence data were derived from a long-term survey, and the survey locations covered all major regions of China.
In summary, we found pest abundance/distribution and host range showed significant phylogenetic signal. Relative abundance and geographic range size of spider mites were positively correlated with host range. These results suggest that phylogenetic information could help to understand the community assembly of this pest group from an evolutionary perspective. Information that is easily obtainable – including the number of known hosts and phylogenetic position of the mites –may contribute to risk analyses of pest outbreaks.
Distribution, relative abundance and host range for species from China
We collected spider mites during the summers of 2008 to 2017 across major regions of China (Fig. 1, also see Additional file 1: Table S1 in ESM). Since the number of spider mites at a collection site is affected by many local factors including pesticide application, host type and sample period , the total abundance of mites is expected to differ even among nearby sites. In contrast, the species composition of spider mites at a larger scale is relatively stable . We therefore focused on surveys of multiple sites and estimated relative abundance based on occurrences across sites rather than resampling each site multiple times . At each site (around 3000 m2), our strategy was to collect a maximum of three mites per plant, with plants separated by a minimum of 1 m. Overall, 318 geographic or host-associated populations were collected from 180 sites that spanned the native range of spider mites, from Northeastern China to Southwestern China (Fig. 1). As a metric of range size, we calculated the latitudinal span covered by each species . Relative abundance for each species was represented by the total number of occurrences (NOC) in our survey of 318 populations. To minimize any bias associated with intensive sampling in one site, average frequency of occurrence (AF) across different sites was calculated as a second index of mites abundance.
In general, monophagous species tend to feed on a single plant species, oligophagous species tend to feed on one genus and polyphagous species feed on at least one family of plants . Some spider mite species (e.g., Tetranychus evansi) can feed on more than 300 host species, but most hosts belong to one family . Other mites may have very few host species, but the hosts belong to more than one host family . We therefore used both host plant species records and host plant family records to represent host range. The host records for each species were obtained from the survey.
Distribution, occurrence and host range for species from a global dataset
Spider mites web (http://www1.montpellier.inra.fr/CBGP/spmweb/) is a global database which includes host records, distribution countries and historical records for more than 1300 spider mites species . The frequency of occurrence for each species at a global scale was indirectly counted as the number of historical records that had clear host information on Spider mites web . This method may overestimate or underestimate the abundance of a mite species in the wild because of a likely focus on economically important species, but it is still likely to reflect the relative abundance of species within agricultural settings. Host range at a global scale was represented by host family number, which was derived from Spider mites web . Most of the species lacked detailed location information, and the distribution range for each species was estimated from the number of countries where each species was found.
Phylogeny and phylogenetic signal analyses
The 18S gene, 5′ end of the 28S rRNA gene and mitochondrial COI gene were used for phylogenetic analyses  (for GenBank accession ID, see Additional file 1: Table S2 in ESM). Phylogenetic trees were constructed using maximum likelihood (ML) and Bayesian inference (BI) methods following the protocols described by Xue et al. . BI analyses were performed with MrBayes 3.2.2 , and two independent runs were conducted, each with four Markov Chains (one cold chain and three heated chains). GTR + I + G was the model chosen by jModelTest 2.1.1 . ML analyses were performed using the GTRGAMMAI model in raxmlHPC-PTHREADS  implemented in raxmlGUI1.3 . Genetic distances between species were calculated in MEGA7  applying the Kimura 2-parameter model , with 1000 bootstrap replicates.
To test whether species occurrence (relative abundance and distribution) and host range were non-randomly associated with genetic similarity between the species, we used Mantel tests to compare these characteristics with a genetic distance matrix  in R-3.4.4 for Windows (R ). To quantify the phylogenetic signal of species characters, we computed Blomberg’s K  in the package ‘picante’  and Abouheif’s test  in package ‘adephylo’ . Both tests were performed in R-3.4.4 for Windows . Blomberg’s K quantifies the amount of phylogenetic signal in the data relative to a Brownian motion model of trait evolution. K = 1 corresponds to a Brownian motion pattern and K = 0 corresponds to a random distribution of the trait across the phylogeny. The higher the K statistic, the more phylogenetic signal in a trait. Traits with PIC.variance probabilities < 0.05 have significant phylogenetic signal. PIC.variance probabilities is the quantile of the observed phylogenetically independent contrast variance versus the null distribution, which can be used as a one-tailed p-value to test for greater phylogenetic signal than expected . Abouheif’s test for serial independence is based on the sum of the successive squared differences between trait values of neighboring species . Traits with probabilities < 0.05 were considered phylogenetically structured.
To test whether we could use phylogenic distance to predict pest occurrence for each genus, we first identified the most abundant species as focal species – the species for which the measured response was strongest . Then, we calculated the genetic distance between the focal species and other spider mites, respectively. The relationships between species occurrence (abundance and distribution) and phylogenetic distance were investigated using Pearson correlation analyses.
Host range and pest occurrence relationships
To test for relationships between host range and species occurrence, we computed correlations between these variables (Pearson’s r). The PGLS (Phylogenetic generalized linear models)  function in the ‘caper’ package  and PIC (phylogenetically independent contrasts function)  in the ‘ape’ package  were then used to calculate phylogenetically-corrected correlation coefficients for host range and abundance accounting for variable levels of phylogenetic signal. Both programs provide a phylogenetically corrected r value giving an estimation of the association between the host range and abundance variables following correction for phylogeny. To illustrate how host range evolved within the evolutionary history of spider mites, we reconstruct ancestral states for host range using a maximum likelihood approach based on a BM model in the ‘geiger’ package . Data (Relative abundance, host range and global distribution countries number) were log-transformed to meet requirements for normality in all analyses.
Availability of data and materials
The phylogenetic trees generated during the current study are available in the Treebase: http://purl.org/phylo/treebase/phylows/study/TB2:S25372. Sequences can be retrieved from Genbank with accession numbers provided in Additional file 1: Table S2. The other data supporting findings are contained within the manuscript and in the supplemental files. If any additional information is necessary, please contact the corresponding author.
Average frequency of occurrence
Number of distributed countries
Number of host families
Total number of occurrences
Phylogenetic generalized linear models
Phylogenetically independent contrasts function
Abouheif A. A method for testing the assumption of phylogenetic independance in comparative data. Evol Ecol Res. 1999;1:895–909.
Angert AL, Crozier LG, Rissler LJ, et al. Do species’ traits predict recent shifts at expanding range edges? Ecol Lett. 2011;14:677–89. https://doi.org/10.1111/j.1461-0248.2011.01620.x.
Arbetman MP, Gleiser G, Morales CL, et al. Global decline of bumblebees is phylogenetically structured and inversely related to species range size and pathogen incidence. Proceeding R Soc B Biol Sci. 2017;284:20170204. https://doi.org/10.1098/rspb.2017.0204.
Aronson MFJ, La Sorte FA, Nilon CH, et al. A global analysis of the impacts of urbanization on bird and plant diversity reveals key anthropogenic drivers. Proc R Soc B Biol Sci. 2014;281:20133330. https://doi.org/10.1098/rspb.2013.3330.
Blair RB. Land use and avian species diversity along an urban gradient. Ecol Appl. 1996;6:506–19. https://doi.org/10.2307/2269387.
Blomberg SP, Garland T, Ives AR. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution. 2003;57:717–45. https://doi.org/10.1111/j.0014-3820.2003.tb00285.x.
Boubou A, Migeon A, Roderick GK, Navajas M. Recent emergence and worldwide spread of the red tomato spider mite, Tetranychus evansi: genetic variation and multiple cryptic invasions. Biol Invasions. 2010;13:81–92. https://doi.org/10.1007/s10530-010-9791-y.
Boulangeat I, Lavergne S, Van Es J, et al. Niche breadth, rarity and ecological characteristics within a regional flora spanning large environmental gradients. J Biogeogr. 2012;39:204–14. https://doi.org/10.1111/j.1365-2699.2011.02581.x.
Broennimann O, Thuiller W, Hughes G, et al. Do geographic distribution, niche property and life form explain plants’ vulnerability to global change? Glob Chang Biol. 2006;12:1079–93. https://doi.org/10.1111/j.1365-2486.2006.01157.x.
Ceballos G, Ehrlich PR, Barnosky AD, et al. Accelerated modern human – induced species losses: entering the sixth mass extinction. Sci Adv. 2015;1:1–5. https://doi.org/10.1126/sciadv.1400253.
Cho S, Mitchell A, Mitter C, et al. Molecular phylogenetics of heliothine moths (Lepidoptera: Noctuidae: Heliothinae), with comments on the evolution of host range and pest status. Syst Entomol. 2008;33:581–94. https://doi.org/10.1111/j.1365-3113.2008.00427.x.
Darriba D, Taboada GL, Doallo R, Posada D. JModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772. https://doi.org/10.1038/nmeth.2109.
Davis EB. Comparison of climate space and phylogeny of Marmota (Mammalia: Rodentia) indicates a connection between evolutionary history and climate preference. Proc R Soc B Biol Sci. 2005;272:519–26. https://doi.org/10.1098/rspb.2004.2979.
Felsenstein J. Phylogenies and the comparative method. Am Nat. 1985;125:1–15. https://doi.org/10.1086/284325.
Forister ML, Novotny V, Panorska AK, et al. The global distribution of diet breadth in insect herbivores. Proc Natl Acad Sci. 2015;112:442–7. https://doi.org/10.1073/pnas.1423042112.
Franklin J. Species distribution models in conservation biogeography: developments and challenges. Divers Distrib. 2013;19:1217–23. https://doi.org/10.1111/ddi.12125.
Freckleton RP, Harvey PH, Pagel M. Phylogenetic analysis and comparative data: a test and review of evidence. Am Nat. 2002;160:712–26. https://doi.org/10.1086/343873.
Fritz SA, Purvis A. Selectivity in mammalian extinction risk and threat types: a new measure of phylogenetic signal strength in binary traits. Conserv Biol. 2010;24:1042–51. https://doi.org/10.1111/j.1523-1739.2010.01455.x.
Gallagher AJ, Hammerschlag N, Cooke SJ, et al. Evolutionary theory as a tool for predicting extinction risk. Trends Ecol Evol. 2015;30:61–5. https://doi.org/10.1016/j.tree.2014.12.001.
Ganaha-kikumura T, Ohno S, Kijima K, et al. Species composition of thrips (Thysanoptera: Thripidae) and spider mites (Acari: Tetranychidae) on cultivated chrysanthemum (Asteraceae) in Okinawa, southwestern Japan. Entomol Sci. 2012;15:232–7. https://doi.org/10.1111/j.1479-8298.2011.00501.x.
Gaston KJ, Blackburn TM, Lawton JH. Interspecific abundance-range size relationships: an appraisal of mechanisms. J Anim Ecol. 1997;66:579–601. https://doi.org/10.2307/5951.
Gilbert GS, Briggs HM, Magarey R. The impact of plant enemies shows a phylogenetic signal. PLoS One. 2015;10:1–11. https://doi.org/10.1371/journal.pone.0123758.
Gilbert GS, Magarey R, Suiter K, Webb CO. Evolutionary tools for phytosanitary risk analysis: phylogenetic signal as a predictor of host range of plant pests and pathogens. Evol Appl. 2012;5:869–78. https://doi.org/10.1111/j.1752-4571.2012.00265.x.
Gilbert GS, Webb CO. Phylogenetic signal in plant pathogen-host range. Proc Natl Acad Sci. 2007;104:4979–83. https://doi.org/10.1073/pnas.0607968104.
Godefroid M, Rasplus J-Y, Rossi J-P. Is phylogeography helpful for invasive species risk assessment? The case study of the bark beetle genus Dendroctonus. Ecography. 2016;39:1197–209. https://doi.org/10.1111/ecog.01474.
Grbić M, Van Leeuwen T, Clark RM, et al. The genome of Tetranychus urticae reveals herbivorous pest adaptations. Nature. 2011;479:487–92. https://doi.org/10.1038/nature10640.
Gregory RD, van Strien A, Vorisek P, et al. Developing indicators for European birds. Philos Trans R Soc B Biol Sci. 2005;360:269–88. https://doi.org/10.1098/rstb.2004.1602.
Halstead NT, Hoover CM, Arakala A, et al. Agrochemicals increase risk of human schistosomiasis by supporting higher densities of intermediate hosts. Nat Commun. 2018;9:837. https://doi.org/10.1038/s41467-018-03189-w.
Harmon LJ, Weir JT, Brock CD, et al. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24:129–31. https://doi.org/10.1093/bioinformatics/btm538.
Hawkins BA, Devries PJ. Tropical niche conservatism and the species richness gradient of north American butterflies. J Biogeogr. 2009;36:1698–711. https://doi.org/10.1111/j.1365-2699.2009.02119.x.
Hong X-Y. Agricultural acarology. Beijing: China Agriculture Press; 2012.
Jarzyna MA, Jetz W. Detecting the multiple facets of biodiversity. Trends Ecol Evol. 2016;31:527–38. https://doi.org/10.1016/j.tree.2016.04.002.
Jin P-Y, Tian L, Chen L, Hong X-Y. Spider mites of agricultural importance in China, with focus on species composition during the last decade (2008–2017). Syst Appl Acarol. 2018;23:2087–98. https://doi.org/10.11158/saa.23.11.1.
Jombart T, Dray S. Adephylo: exploratory analyses for the phylogenetic comparative method. Bioinformatics. 2010;26:1–21. https://doi.org/10.1093/bioinformatics/btq292.
Kellermann V, Loeschcke V, Hoffmann AA, et al. Phylogenetic constraints in key functional traits behind species’ climate niches: patterns of desiccation and cold resistance across 95 Drosophila species. Evolution. 2012;66:3377–89. https://doi.org/10.1111/j.1558-5646.2012.01685.x.
Kembel SW, Cowan PD, Helmus MR, et al. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26:1463–4. https://doi.org/10.1093/bioinformatics/btq166.
Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20. https://doi.org/10.1007/BF01731581.
Kolb A, Barsch F, Diekmann M. Determinants of local abundance and range size in forest vascular plants. Glob Ecol Biogeogr. 2006;15:237–47. https://doi.org/10.1111/j.1466-822X.2005.00210.x.
Krasnov BR, Poulin R, Mouillot D. Scale-dependence of phylogenetic signal in ecological traits of ectoparasites. Ecography. 2011;34:114–22. https://doi.org/10.1111/j.1600-0587.2010.06502.x.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4. https://doi.org/10.1093/molbev/msw054.
Le Viol I, Jiguet F, Brotons L, et al. More and more generalists: two decades of changes in the European avifauna. Biol Lett. 2012;8:780–2. https://doi.org/10.1098/rsbl.2012.0496.
Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27:209–20. https://doi.org/10.1088/1751-8113/44/8/085201.
Matsuda T, Morishita M, Hinomoto N, Gotoh T. Phylogenetic analysis of the spider mite sub-family tetranychinae (Acari: Tetranychidae) based on the mitochondrial COI gene and the 18S and the 5′ end of the 28S rRNA genes indicates that several genera are polyphyletic. PLoS One. 2014;9:e108672. https://doi.org/10.1371/journal.pone.0108672.
McKinney ML, Lockwood JL. Biotic homogenization: a few winners replacing many losers in the next mass extinction. Trends Ecol Evol. 1999;14:450–3. https://doi.org/10.1016/S0169-5347(99)01679-1.
Migeon A, Nouguier E, Dorkeld F. Spider mites web: a comprehensive database for the Tetranychidae. In: Trends in Acarology; 2011. p. 557–60.
Miller JR. Biodiversity conservation and the extinction of experience. Trends Ecol Evol. 2005;20:430–4. https://doi.org/10.1016/j.tree.2005.05.013.
Navajas M, de Moraes GJ, Auger P, Migeon A. Review of the invasion of Tetranychus evansi: biology, colonization pathways, potential expansion and prospects for biological control. Exp Appl Acarol. 2013;59:43–65. https://doi.org/10.1007/s10493-012-9590-5.
Orme D, Freckleton R, Thomas G et al (2013) Caper: comparative analyses of Phylogenetics and evolution in R. R package (version 0.5.2).
Passy SI. A hierarchical theory of macroecology. Ecol Lett. 2012;15:923–34. https://doi.org/10.1111/j.1461-0248.2012.01809.x.
Petit C, Dupas S, Thiéry D, et al. Do the mechanisms modulating host preference in holometabolous phytophagous insects depend on their host plant specialization? A quantitative literature analysis. J Pest Sci. 2017;90:797–805. https://doi.org/10.1007/s10340-017-0833-4.
Popescu AA, Huber KT, Paradis E. Ape 3.0: new tools for distance-based phylogenetics and evolutionary analysis in R. Bioinformatics. 2012;28:1536–7. https://doi.org/10.1093/bioinformatics/bts184.
Purvis A. Phylogenetic approaches to the study of extinction. Annu Rev Ecol Evol Syst. 2008;39:301–19. https://doi.org/10.1146/annurev-ecolsys-063008-1.
R Development Core Team. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.
Raje KR, Ferris VR, Holland JD. Phylogenetic signal and potential for invasiveness. Agric For Entomol. 2016;18:260–9. https://doi.org/10.1111/afe.12158.
Robles-Fernández ÁL, Lira-Noriega A. Combining phylogenetic and occurrence information for risk assessment of pest and pathogen interactions with host plants. Front Appl Math Stat. 2017;3:1–9. https://doi.org/10.3389/fams.2017.00017.
Ronquist F, Teslenko M, Van Der Mark P, et al. Mrbayes 3.2: efficient bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42. https://doi.org/10.1093/sysbio/sys029.
Sekar S. A meta-analysis of the traits affecting dispersal ability in butterflies: can wingspan be used as a proxy? J Anim Ecol. 2012;81:174–84. https://doi.org/10.1111/j.1365-2656.2011.01909.x.
Sih A, Ferrari MCO, Harris DJ. Evolution and behavioural responses to human-induced rapid environmental change. Evol Appl. 2011;4:367–87. https://doi.org/10.1111/j.1752-4571.2010.00166.x.
Silvestro D, Michalak I. raxmlGUI : a graphical front-end for RaxML. Org Divers Evol. 2011;12:335–7. https://doi.org/10.1007/s13127-011-0056-0.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90. https://doi.org/10.1093/bioinformatics/btl446.
Sullivan AP, Bird DW, Perry GH. Human behaviour as a long-term ecological driver of non-human evolution. Nat Ecol Evol. 2017;1:0065. https://doi.org/10.1038/s41559-016-0065.
Sun J-T, Lian C, Navajas M, Hong X-Y. Microsatellites reveal a strong subdivision of genetic structure in Chinese populations of the mite Tetranychus urticae Koch (Acari: Tetranychidae). BMC Genet. 2012;13:8. https://doi.org/10.1186/1471-2156-13-8.
Van Petegem K, Boeye J, Stoks R, Bonte D. Spatial selection and local adaptation jointly shape life-history evolution during range expansion. Am Nat. 2016;188:485–98. https://doi.org/10.1086/688666.
Van Petegem K, Moerman F, Dahirel M, et al. Kin competition accelerates experimental range expansion in an arthropod herbivore. Ecol Lett. 2018;21:225–34. https://doi.org/10.1111/ele.12887.
Verberk WCEP, van der Velde G, Esselink H. Explaining abundance-occupancy relationships in specialists and generalists: a case study on aquatic macroinvertebrates in standing waters. J Anim Ecol. 2010;79:589–601. https://doi.org/10.1111/j.1365-2656.2010.01660.x.
Williams SE, Shoo LP, Isaac JL, et al. Towards an integrated framework for assessing the vulnerability of species to climate change. PLoS Biol. 2008;6:e325. https://doi.org/10.1371/journal.pbio.0060325.
Willis CG, Ruhfel B, Primack RB, et al. Phylogenetic patterns of species loss in Thoreau’s woods are driven by climate change. Proc Natl Acad Sci U S A. 2008;105:17029–33. https://doi.org/10.1073/pnas.0806446105.
Wilson PD. The pervasive influence of sampling and methodological artefacts on a macroecological pattern: the abundance-occupancy relationship. Glob Ecol Biogeogr. 2008;17:457–64. https://doi.org/10.1111/j.1466-8238.2008.00385.x.
Wright DH. Correlations between incidence and abundance are expected by chance. J Biogeogr. 1991;18:463–6. https://doi.org/10.2307/2845487.
Xue XF, Dong Y, Deng W, et al. The phylogenetic position of eriophyoid mites (superfamily Eriophyoidea) in Acariformes inferred from the sequences of mitochondrial genomes and nuclear small subunit (18S) rRNA gene. Mol Phylogenet Evol. 2017;109:271–82. https://doi.org/10.1016/j.ympev.2017.01.009.
We sincerely thank Hao-Sen Li, Si-Xiang Yang, Tong-Pu Li, Xing-Zhi Duan, Xu Zhang, Xue Xia, Yan-Kai Zhang and Yu-Xi Zhu of Nanjing Agricultural University, China for their help with the sample collection and species identification. We also thank Dr. Xiao-Li Bing of Nanjing Agricultural University, China for his helpful comments on the manuscript. We also wish to thank the Associate Editor for helpful comments and the Anonymous Reviewers who helped us improve the manuscript.
This study was supported in part by a grant-in-aid for Scientific Research (31672035, 31871976) from the National Natural Science Foundation of China (NSFC). Opinions, findings, conclusions and recommendations expressed in this publication is that of the authors, and the NSFC accepts no liability whatsoever in this regard.
Ethics approval and consent to participate
This article does not contain any studies with human participants or animals performed by any of the authors.
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.
Additional file 1: Table S1. Sample information for samples from China used in this study. Table S2. GenBank accession of sequences used in phylogenetic reconstruction at global scale. Table S3. Relative abundance, host range and distribution of each spider mite species in China. Table S4. Mantel tests of association between genetic distance and species abundance, distribution range and host range. Table S5. Pearson correlations between species occurrence and genetic distance to the focal species of different genera. Figure S1. Bayesian tree with posterior probabilities (a) and RAxML ML tree with bootstrap proportions from 1000 rapid bootstrap replicates (b) for Tetranychus species in China. Figure S2. Bayesian tree with posterior probabilities of 88 species. Figure S3. ML trees based on 1000 rapid bootstrap replicates of 88 species. Figure S4. Correlations between host specialization and local occurrence in the China dataset (a) and the global dataset (b).
About this article
Cite this article
Jin, P., Sun, J., Hoffmann, A. et al. Phylogenetic signals in pest abundance and distribution range of spider mites. BMC Evol Biol 19, 223 (2019) doi:10.1186/s12862-019-1548-3
- Pest occurrence
- Phylogenetic signal
- Host range
- Spider mite