Great molecular variation within the species Phytoseius finitimus (Acari: Phytoseiidae): implications for diagnosis decision within the mite family Phytoseiidae

Molecular markers are increasingly used for species identification and new taxa description. However, rules to determine frontiers between populations and species are not clear depending on taxa considered. For mites, few studies deal with molecular diagnoses, making rules for associated decision difficult. The present study focuses on a species of the predatory mite family Phytoseiidae (Phytoseius finitimus), considered for biological control of mites and small insect pests in fruit orchards and vineyards in the Mediterranean basin. This paper aims to elucidate the causes of great molecular variations and questions the occurrence of cryptic species. Molecular (12S rRNA, CytB mtDNA, ITSS) and morphological analyses were performed on four populations collected in Corsica and Italy in crops (vine and kiwi) and in an uncultivated environment (Viburnum lantana). Different methods for identifying species have been used (tree approaches, distances and ABGD algorithms). A reference database of distances within and between Phytoseiidae species has been elaborated to inform the present question and to assist with further diagnosis within Acari. Mitochondrial DNA analyses show that specimens from V. lantana were well separated from the three other populations with high genetic distances, suggesting the existence of a cryptic species. Molecular ITSS analyses coupled with morphological features show however that the four populations seem to belong to the same species. The great mitochondrial polymorphism is discussed in regards to: (i) genetic distances reported for Phytoseiidae species and (ii) potential biological differences between populations (cultivated versus uncultivated areas). This study clearly emphasizes the necessity of integrative taxonomy approaches for diagnosis decisions. Furthermore, based on the polymorphism herein detected, maximal intraspecific distances are proposed (9, 23 and 2.8 % for 12S rRNA, CytB mtDNA and ITSS) for diagnosis decisions within Phytoseiidae. Further statistical analyses are however clearly required to determine statistical error for general and reliable decision making.


INTRODUCTION
Molecular markers for differentiating microorganism, plant and animal species are currently used in taxonomy. The rapid technical developments (i.e. DNA extraction, PCR and Sanger sequencing) during these last decades as well barcoding approaches and associated conceptual definitions have led to the development of molecular diagnosis (i.e. Hebert et al. 2003a,b;2004a,b;DeSalle et al. 2005;Goldstein and DeSalle 2011). Molecular markers were traditionally used when doubts on species identity arise from morphological comparisons and sometimes from biological features between populations. However, cryptic species are more and more revealed via comparison of DNA sequences of apparently morphologically similar specimens (i.e. Pons et al. 2006;Witt et al. 2006;Bickford et al. 2007;Valentini et al. 2009;Foottit et al. 2009;Pauls et al. 2010;Jorger et al. 2012;Kekkonen and Hebert, 2014;Kekkonen et al. 2015). The present study focuses on a species of the mite family Phytoseiidae: Phytoseius finitimus Ribaga. This family, containing 2,479 valid species (Demite et al. 2016), is one of the most well-known within the Mesostigmata, due to the predatory abilities of numerous species to control mite and small insect pests (McMurtry and Croft 1997;Gerson et al. 2003;McMurtry et al. 2013). Correct identification of these predators is clearly of huge importance for ensuring the success of biological control, as different species have different biological features (i.e. Tixier et al. 2006;McMurtry et al. 2013). Phytoseiidae species identification is essentially based on female morphological characters. However, the reliability of those characters as well as the low number of morphological features potentially available (due to the mite small size, 500 µm maximum in length) have led to the development of molecular studies to resolve questions concerning synonymies (i.e. Jeyaprakash and Hoy 2002;Noronha et al. 2003;Tixier et al. 2006Tixier et al. , 2008aTixier et al. , 2010aTixier et al. , 2011bTixier et al. , 2012Tixier et al. , 2014Okassa et al. 2009Okassa et al. , 2010Okassa et al. , 2011Okassa et al. , 2012Kanouh et al. 2010a,b;Bowman and Hoy 2012;Chao et al. 2012;Navia et al. 2014;Rezende et al. 2015). For ensuring molecular diagnosis, a database containing DNA sequences associated to the correct species name is essential (i.e. DeSalle et al. 2005;Collins and Cruickshank 2013). When elaborating such a database for Phytoseiidae species, high molecular variations were observed within specimens identified as P. finitimus during samplings carried out in Italy and France, two areas where this species is most commonly found.
Phytoseius finitimus is a generalist predator and is quite common in Mediterranean vineyards (i.e. Ragusa and Ciulla 1991;Duso and Vettorazzo 1999;Duso and Fontana 2002;Peverieri et al. 2009;Miñarro and Kreiter 2012). Some studies have shown its ability to control phytophagous mites, whiteflies and thrips and to develop when fed with pollen (Pappas et al. 2013) and Ahmad et al. (2015) recently studied its interspecific competition features. The present study aims to elucidate the great molecular distances observed between four P. finitimus populations and to determine how these great variations are related or not to cryptic species occurrence. Answering this question is clearly important for biological control applications, as implementations and advice for agro-environmental management would be different if one or two species are named under P. finitimus. To determine if one or several species were present, an integrative taxonomic approach based on three molecular fragments and morphological analyses has been carried out. The number of analytical developments for assessing species within a DNA sequence dataset is increasing for ensuring diagnosis decision (i.e. Kekkonen et al. 2015). Traditional cluster tree-based and distance threshold-based approach as well as the overlapping ABGD (Automatic Barcode Gap Discovery) algorithm have been herein used to determine their respective performance and congruence. Furthermore, a compilation of genetic distances between Phytoseiidae species retrieved from already published works is also provided as a reference comparison point. This review aims to provide additional information that can be used in informing the decision and more broadly aims to assist further analyses on DNA molecular diagnosis within the family Phytoseiidae for delimiting intra-and interspecific frontiers.

Species and populations studied
Phytoseius finitimus, even if reported many times in the literature, is not naturally observed in great densities and is quite difficult to efficiently retrieve. It is known to be frequent in Italy and in Corsica (France), especially in vineyards (Kreiter et al. 2000;Duso and Fontana 2002;Demite et al. 2014). Four populations of P. finitimus were considered, with the double aims to study populations (i) from the two main areas where this species is most abundant (Italy and France) and (ii) from cultivated (vine and kiwi) and uncultivated areas (colonized by Viburnum lantana L.). Two populations were collected in Italy in the same locality on vines (Vitis vinifera L.) and on uncultivated V. lantana plants located besides the vine plot. Two populations were collected in France -Corsica, on vine and kiwi (Actinidia deliciosa C.F. Liang & A.R. Ferguson.) on plots separated by approximately 15 kms ( Figure 1). The characteristics of the populations studied and accession numbers of DNA sequences in the Genbank database are presented in Table 1. The species Typhlodromus (Typhlodromus) phialatus Athias-Henriot, belonging to the sub-family Typhlodrominae, was used as an out-group for tree construction. Sequences of this latter species were obtained by present authors and already referenced in the Genbank database (12S rRNA: HM635274, ITSS: HQ404829, CytB mtDNA: HM635300) (Okassa et al. 2012;Tsolakis et al., 2012).

DNA extraction.
Total genomic DNA of a single female was extracted using a Qiagen DNeasy tissue kit (Qiagen, Hilden, Germany), according to the DNA extraction protocol « Purification of Total DNA from Animal Blood or Cells » (Spin-Column Protocol) adapted for extracting total DNA from mites (Kanouh et al. 2010b). Between six and thirteen females were considered per population (Table 1). After extraction, specimens were retrieved and mounted on slides according the method developed by Tixier et al. (2010 b).

Data analyses.
A preliminary analysis was conducted on the coding sequences (CytB mtDNA) to check for the absence of stop codons. Three molecular diagnosis approaches were used to assess the occurrence of cryptic species within the four populations of P. finitimus.

Tree-based cluster approach.
Neighbour-joining trees were constructed using the Kimura 2-parameter model and node support was determined using 1,000 bootstrap replicates calculated in Mega 6.0.6® (Tamura et al. 2013). The K2P model was used to make our results comparable with other studies on Phytoseiidae. Furthermore to determine how other tree models perform for resolving the question on cryptic species, Parsimony and Bayesian analyses were carried out for each fragment and combined dataset (applying simple partition procedure with one partition corresponding to each DNA marker). As for some specimens, it was impossible to obtain PCR products (especially for ITS), the combined data set was analysed following Wiens (1998a,b), Wiens et al. (2005) and Zheng and Wiens (2016), who demonstrated that missing data do not affect the phylogenetic accuracy.
The Incongruence Length Difference (ILD) P value (P= 0.07) suggests no conflicting signals between markers (Mickevich and Farris 1981;Sullivan 1996;Farris et al. 1994;Cunningham 1997).For Parsimony analyses a heuristic search procedure repeated 1,000 times, was performed with randomized taxa additions and branch-swapping algorithm (TBR). To reduce misleading effects of homoplasious characters, an a posteriori re-weighting was applied according to the rescaled consistency index (RC) after each tree search, until the number of trees stabilised (Farris 1969(Farris , 1989 (PAUP*, v.4.0b.10 / Swofford, 2002). Bayesian analysis was carried out using MrBayes 3.1 (Ronquist and Huelsenbeck, 2003); the best-fit-substitution models were determined by Modeltest 3.07 in PAUP (Posada and Crandall 1998) through hierarchical likelihood-ratio As tree clades do not necessary reflect species clusters but can represent populations of a same species (i.e. Moritz and Cicero 2004;Collins and Cruickshank 2013), the classical threshold-based approach was applied and distance matrices (using the Kimura 2-parameter) model were elaborated. With this approach, ideally no overlap should occur between intra and interspecific distances (i.e. Meyer and Paulay 2005). Even if the no-overlap model is sometimes point of contention especially because of sampling completion, knowledge of species considered and associated taxonomical features, distance comparison is currently used to distinguish between species for practical assignments. Two approaches for comparing distances were applied.

The overlap distance criteria.
First, distances within and between the four populations were manually compared to determine gap occurrence. In addition, the ABGD (Automatic Barcode Gap Discovery) algorithm was applied.
It aims at automatically determining the occurrence of different species within a dataset based on gap between batches of sequences pairs (Puillandre et al. 2012a), detecting the first barcode gap and using it to partition the dataset into candidate species. Analyses were performed in November 2016 on the web interface (http://wwwabi.snv.jussieu.fr/public/abgd/). They were conducted using three metrics (Jukes-Cantor (JC), Kimura 2 parameter (K2P) and simple p-distances) and default parameter values.

Empirical comparisons to references.
Distances usually constitute a useful point of comparison when the taxonomic status of the taxa compared is sure (Robinson et al. 2009). Thus, the genetic distances herein obtained were compared to those mentioned in 15 publications for 24 Phytoseiidae species. Only distances issued from published works and not calculated from all Phytoseiidae sequences deposited in Genbank were considered, because of great incertitude on the true identity of the sequences deposited in this world database (i.e. Tixier et al. 2011a). Two distances were compiled: (i) one corresponding to the maximal value observed within a same species (maximal intraspecific variation) and (ii) one corresponding the minimal value observed between two species of a same genus [minimal interspecific variation assumed to be the "Nearest Neighbour" (Kekkonen et al. 2015)].

Characters considered.
As female dorsal seta lengths are currently used to discriminate between Phytoseiidae species (Chant and McMurtry 2007), the seventeen dorsal idiosomal setae were measured for 5 to 10 females issued from the same populations as sequenced specimens: j1, j3, j4, j5, j6, J2, J5, z2, z3, z5, z4, Z4, Z5, s4, s6, r3 and R1. The lengths of the seta JV5 and the macroseta on the basitarsus of the leg IV (St IV) were also considered. Finally, peritreme, spermatheca, ventrianal shield and chelicera shapes were observed. Terminology for seta notation follows that of Lindquist and Evans (1965) as adapted by Rowell et al. (1978) for Phytoseiidae. Data analyses. The data corresponding to the 17 continuous characters were normally distributed and variances were equal (Bartlett test). ANOVA analyses were performed for each character to determine differences between the four populations. A principal component analysis was also carried out to determine whether the combination of the morphological characters differentiates specimens considered. All statistical analyses were carried out with Statistica (Statsoft France 2010).
To determine the variation of each seta length, the rCL 95 corresponding to the percentage of variation around the mean was calculated according to the formula provided in Tixier (2012). Furthermore, min and max values of the interval assumed to include 95 % of specimens of a same species for each setal length were modelled using the abacus provided in Tixier (2012) to determine whether the observed variation would be included in the modelled "normal" intraspecific variation.

Molecular analyses.
437, 418 and 690 base pairs (bp) were aligned for the 12S rRNA, Cytb mtDNA, and ITSS genes, respectively. Quite similar and constant rates of nucleotide substitutions were observed for all the populations studied. A BLAST search in Genbank showed that the sequences aligned with those of Phytoseiidae. The overall genetic distances between the specimens of P. finitimus for 12S rRNA (mean = 2.9 %, min-max = 0-9 %) and CytB mtDNA (mean = 7.6 %, min-max = 0-23 %) are high in regards to those obtained in other studies (Table 3) questioning the occurrence of cryptic species.

Diagnosis decision and tree cluster criteria.
The same tree topologies were observed with the three algorithms (Neighbour Joining, Parsimony, Bayesian analyses) for the three markers considered. Three clades were observed on the trees obtained with the two mitochondrial markers (12S rRNA, CytB mtDNA): one containing specimens collected in Corsica on vine, one containing specimens collected in Italy on vine and in Corsica on kiwi and   Figures 2, 3). One specimen collected on vine in Corsica was however located in the clade containing specimens from vine-Italy and kiwi-Corsica. These clades could be three species or three populations of a same species. Specimens included in the two former clades are more closer to each other than to specimens included in the third. Two not well-sustained clades are observed on the ITSS tree ( Figure 4), one containing specimens from V. lantana and one containing the specimens of the three other populations. Furthermore, the same clades as those obtained with mitochondrial markers were also obtained when concatenated data were considered for both Parsimony and Bayesian analyses ( Figure 5).

Diagnosis decision and the overlap criteria.
Figure 6 presents the four intra-populations distances and the six inter-population distances. For 12S rRNA, overlap is observed between intrapopulation distances and distances between specimens collected on vine in Italy, vine and kiwi in Corsica. However, no overlap is observed between intra-population distances and distances between specimens collected on V. lantana in Italy and those from the three other populations ( Figure 6). For CytB mtDNA, overlap is observed between intrapopulation and all inter-population distances, especially because of the high intra-population distances observed within the vine-France population. Finally, for ITSS marker, overlap is observed between all intra and inter-population distances.
With the ABGD analyses, the same groups were obtained for the Jukes-Cantor (JC69) and Kimura (K80) TS/TV distances but different clusters were observed for the P distance (Table 3). For the 12S rRNA and CytB mtDNA markers, four and five groups were observed with the K80 and JC69 distances, respectively. Globally, the same groups as those observed on the phylogenetic trees were obtained, but more clustering was observed within the population collected on vine in Corsica (two clusters for 12S rRNA and three clusters for CytB mtDNA). Using P distance, only two groups were observed with the two mitochondrial markers, one containing specimens collected on V. lantana and one containing all the others. For the ITSS marker, one group was observed whatever the distance considered. The ABGD method provides thus different information and groups depending on (i) mitochondrial DNA vs nuclear markers and (ii) the distance considered.

Diagnosis decision and comparison to references in literature.
The 12S molecular distances between clades I and II (0-5 %) correspond to values currently retrieved in literature for specimens belonging to the same species (Table 4). Furthermore, such distances are lower than the minimal distances observed until now between species of a same genus. In one case corresponding to distances between Phytoseiulus persimilis Athias-Henriot and P. macropilis (Banks), the presently observed distances are higher than the minimal distance observed between these two species (4 %). The maximal 12S rRNA distances between V. lantana specimens and specimens of clades I and II (9 %) is higher than the maximal intraspecific distances observed until now within a Phytoseiidae species (8 %) (Tables 2, 4). This value even if close to minimal distances observed between Neoseiulus californicus and N. fallacis (Garman) [10 % in Tixier et al. (2014) and 9.5 % in Jeyaprakash and Hoy (2002)] is however lower than distances observed until now between two Phytoseiidae species of the same genus.
The mean CytB distances between specimens from clades I and II are usually lower than genetic distances observed within a same Phytoseiidae species (Tables 2, 4). However, maximal distances (16 %) observed between specimens from vine in Italy and France, and kiwi from France as well those within vine-Corsica population (13 %), are higher or equal to the maximal intraspecific distances observed until now for Phytoseiidae (13 % for N. aceri). The maximal CytB distances between specimens from V. lantana and the three other populations (23%) are higher than those already observed within a same Phytoseiidae species. This value is just below the minimal interspecific distances observed between species of some genus (23.4 % between T. (T.) phialatus and T. (T.) exhilaratus and 24 % between N. californicus and N. picanus in Okassa et al. (2011Okassa et al. ( , 2012).
The highest ITSS distances are observed between specimens collected on V. lantana and the specimens of the other populations. The distances between clades I and II are very low and included in the variation range of intraspecific distances retrieved in literature (Table 3). The mean distances between these two clades and specimens from V. lantana are usually lower than those observed between different species of a same Phytoseiidae genus. The maximal distance is however similar to the distances obtained between T. ( T.) phialatus and T. (T.) exhilaratus (two morphological close species) by Tixier et al. (2006) whereas Navia et al. (2014) reported a distance of 4.75 % between these two same latter species (Table 4)

Morphological analyses.
Some statistical differences are observed for twelve continuous characters (Table 5). However, the statistical groups show that measurements of specimens from V. lantana are usually not different from specimens from the other three populations. The most frequent significant differences are observed between specimens from Corsica (vine and kiwi) and from Italy (vine and V. lantana). On the two axes of the multifactorial analysis explaining 47 % of the variation (Figure 7), specimens collected on V. lantana are close to those collected on vine-Italy, these two populations being quite well separated from the specimens collected in Corsica, both on vine and kiwi. The variation within and between populations from Italy is lower than that observed within and between populations from Corsica. Mean differences in seta length even if statically significant are however usually low. The most important difference is observed in s6 length, this setae being significantly shorter for specimens collected on V. lantana (mean = 81 µm) than for specimens of the three other populations (88, 89, and 94 µmfor mites collected on vine in Italy, on kiwi in Corsica and on vine in Corsica, respectively).

Viburnum lantana-Italy Vine-Italy
Kiwi-Corsica Vine -Corsica However the absolute value between these seta lengths, ranging between 7 and 13 µm, is according to the modelling of intraspecific seta length range proposed by Tixier (2012Tixier ( , 2013 included in intraspecific variation (Table 4). The rCL 95 calculated according to Tixier (2012) shows that most setae have a variation around the mean lower than 20 %, except for the smallest setae. This result is in accordance with results of Tixier (2012) when modelling the intraspecific variation rate around the mean for setae length of Phytoseiidae dorsal shield. These modelled variation intervals around the means (Tixier 2012) include the variation range of specimens herein studied for all seta considered.
Furthermore observations of shapes of spermatheca ( Figure 8) as well as ventrianal shield chaetotaxy and shape and peritremal extremity (reaching z2) do not reveal differences between the populations considered.

DISCUSSION
This study shows (i) the difficulty in interpreting DNA variation for species diagnosis and (ii) the necessity of an integrative taxonomic approach. Clear conclusion can be proposed for two clades (clades I and II) among the three observed on the mitochondrial DNA trees. Because (i) these two clades were not separated using the nuclear ITSS markers, (ii) distance overlap was observed for all markers (nuclear and mitochondrial) between intra and inter-population distances, (iii) these populations belong to a same group with the ABGD method (using P distance), one specimen of vine-Corsica is included in the clade containing specimens from vine-Italy and kiwi-Corsica and (iv) the higher interpopulation distances were much lower than the already observed ones between two species of the same genus, the convergent decision was that these two clades belong to the same species. However, for the status of the third clade containing specimens collected on V. lantana, the decision was more difficult. Because of (i) the well-separated clade including those specimens with both mitochondrial and nuclear markers, (ii) the high distances between those specimens and specimens from the three other populations, (iii) no clear overlap between sequences of these two groups for the 12S rRNA fragment and (iv) results obtained based on mtDNA fragments with the ABGD approach, two hypotheses could be proposed: (i) specimens of the four populations considered belong to a same species characterised by high molecular polymorphisms, or (ii) specimens from V. lantana belong to a species very close to the specimens of the three other populations. The final decision was made based on the observation of morphological characters. No difference between seta lengths of the four populations is observed, suggesting that they belong to the same species.
According to the markers and analytic concepts considered, diagnosis conclusion can be different and decision hard to dress.

Molecular markers and pitfalls.
Mitochondrial markers showed much more variation than was observed with the nuclear marker. Such mitochondrial polymorphism can be due to the nature of molecule evolution (i.e. mtDNA is more sensitive to population bottlenecks, the absence of recombination in mtDNA) and to ecological features (i.e. dispersal, fecundity) that lead to population separations. Furthermore, Phytoseiidae mites are haplodiploid. Even if there are few studies on the impact of such features on mt DNA variation outside of Hymenoptera, Lohse and Ross (2015) suggest that this characteristic can imply a higher introgressions for mitochondrial compared to nuclear markers ; they report current incongruence between nuclear and mitochondrial trees in haploidiploid organisms, including spiders.
As a consequence, some artificial groups can be obtained, underlying the necessity to use an integrative taxonomical approach, combining morpho-logical and molecular characters as well as different DNA fragments and associated conceptual analyses (i.e. Dayrat 2005).

Which analytic methods to use?
This study first stresses the difficult in interpreting phylogenetic trees. The same clades were obtained for the three algorithms used (Neighbour-Joining, Parsimony and Bayesian analyses), suggesting that at low taxonomic level the Neighbour-Joining tree provides the same clustering as the other model-based analyses (Holder and Lewis 2003;Yang and Rannala 2012). In addition the concatenated analysis of the three DNA markers does not provide a clearer answer as the same clusters were obtained as those with the two mitochondrial markers. The concatenated tree thus can hide some information as the signal obtained clearly depends on the resolution level of each DNA marker considered.
This study also stresses the necessity of having clear distance references for comparison. The ABGD aims to construct its own rules based on a data set (without a priori hypothesis on taxon status) (Puillandre et al. 2012a). This method overcomes common threshold applications, reported to be variable according to taxa considered (i.e. Collins and Cruickshank 2013). For the present dataset, more than one species are emphasized using the two mitochondrial markers. Four and five groups are even observed using the JC69 and K80 distances, suggesting (i) the high sensitivity of this clustering method and (ii) the impact of the number of sequences analysed (Puillandre et al. 2012b). Furthermore, even with the ABGD analysis based on P distance, the specimens from V. lantana were always assumed to belong to a different species. However using the ITSS marker, the results converge with other analyses, especially morphology suggesting that all specimens belong to a same species. The number of populations and sequences were quite low which can lead to artificial groups significantly detected by ABGD approach (Lhose 2009;Puillandre et al. 2012b). To our knowledge, this latter method has never been applied for the markers 12S rRNA, CytB mtDNA and only one time for ITS2 (Schwarzfeld and Sperling 2015); it has been essentially used for the traditional barcode fragment (COI mtDNA) and in some rare cases for 16S rRNA (i.e. Puillandre et al. 2012b;Jörger et al. 2012;Guarnizo et al. 2015;Han et al. 2016). Results herein obtained can be due to the fact that other DNA fragments were considered and prior parameters not adapted. Furthermore, it was also carried out only once on mites, on feather mites using COI (Doña et al. 2015) making conclusions of ABGD relevance for mite species delimitation difficult.
This study contributes to the characterization of decision rules for molecular diagnoses within Phytoseiidae and more generally speaking within mites. It is the first time that molecular variation of a species of the subfamily Phytoseiinae is studied. Its variation is higher that already observed within species of the other two subfamilies (Typhlodrominae and Amblyseiinae). Testing additional species of the sub-family Phytoseiinae would be interesting for understanding how subfamily and associated biological/ecological features can impact intraspecific variations (i.e. low dispersal, small population size). The present genetic distances are the highest ever observed within a Phytoseiidae species. Yet, it seems that genetic distances higher than 9 %, 23 % and 2.8 % would correspond to interspecific distances for the markers 12S rRNA, CytB mtDNA and ITSS, respectively. Those results are quite difficult to compare with those obtained for other groups, as the threshold can be different according to taxonomic entities (i.e. Collins and Cruickshank 2013). Furthermore, CytB mt DNA and 12S rRNA markers have been rarely used for species diagnosis within other mite groups. A recent study on ticks using the 12S rRNA markers proposes a minimal interspecific distance of 13.2% and an average intraspecific distance of 1.8%; no data is provided on maximal intraspecific distance (Lv et al. 2014).
These simple comparisons need to be completed with statistical analyses to determine the validity (and associated error) of a barcoding decision depending on the distance obtained between two sequences and to their proximity with the lower limit of intraspecific variation. Furthermore, when comparing with literature, these rules seem to be different depending on the genus as minimal genetic dis-tance of 4 % for the 12S rRNA is observed between two close species (P. persimilis and P. macropilis). However, in this latter case, it is worthwhile to test if: (i) this low intraspecific distance is due to particular biological features of species of the atypical genus Phytoseiulus within Phytoseiidae, or (ii) P. macropilis and P. persimilis are synonyms and a 12S rRNA distance of 4 % would correspond to an intraspecific value, as for all the Phytoseiidae species tested until now.

Implications for biological control?
Finally, for biological control matters, further experiments should be planned to determine whether the present mitochondrial clades show different behaviour and biological features. Tixier et al. (2010a) observed that mitochondrial clades within the species Phytoseiulus longipes Evans were associated with different feeding abilities (prey preferences). The closest populations with molecular markers correspond to specimens collected in cultivated areas (vine-Italy and -Corsica and kiwi-Corsica) but in geographically distinct areas. The well-differentiated population was collected from uncultivated plants (V. lantana) nearby the vine-Italy plot. It is interesting to note that morphological analyses do not show the same similarities, as morphological similarity seems more associated with geographical distances than with plants and cultivated situation. Further investigations would be to determine how selection pressures associated with agricultural practices are linked to such a differentiation for assessing the consequences on colonisation of specimens of P. finitimus dispersing from uncultivated vegetation towards crops.