BACKGROUND
the predictability and abundance of host resources is often assumed to favour ecological specialisation in populations of parasites that have a large spectrum of possible hosts through the disruptive selection of host preference, leading gradually to reproductive isolation and speciation if assortative mating occurs  <cit> .
host-associated biotypes are well documented in plant-feeding insects, and many examples of host races have been described in phytophagous insects in the last ten years.
as reviewed by dres and mallet   <cit> , host races are defined as sets of populations of the same species that use different host taxa in the wild, exhibit host fidelity, coexist in sympatry and correspond to separated genetic clusters among which gene flow still occurs.
moreover, these populations usually display a correlation between host choice and mate choice, have a higher fitness on natal than on alternative hosts and produce hybrids that are less fit than the parental forms.
host race status has been especially discussed within the aphididae  <cit> .
in the case of aphids, host plants exert considerable selective pressure because the plants constitute the aphids' feeding, mating and oviposition sites.
therefore, host specialisation in aphids evolves through the selection of behavioural and chemical mechanisms of host-plant location and recognition as well as through metabolic adaptation to the phloem content of the host plant.
how these adaptive traits evolve in an aphid species depends on the complexity of the aphids' annual life cycle.
like many parasites, aphids have evolved complex life cycles, with successive generations of individuals often specialised to cope with different ecological conditions.
however, aphids are among the few organisms that show plasticity in their reproductive mode as an adaptive response to seasonal changes  <cit> .
the typical aphid annual life cycle, called holocycly or cyclical parthenogenesis, includes a single sexual generation in the autumn composed of males and mating females  that lay resting eggs.
in the spring, a wingless  female  hatches from the egg and produces other parthenogenetic females .
the number of successive viviparae generations in the spring and summer depends on environmental factors.
all of the descendants of a single fundatrix are genotypically identical and constitute a clone, although several phenotypes are expressed.
indeed, the viviparae may be alate or apterous in response to environmental conditions, with development of alate forms stimulated in crowded colonies.
finally, in the autumn, in response to predictive environmental conditions, the viviparae will produce sexual morphs.
approximately 90% of aphid species are monoecious, i.e., they remain on the same host species or closely related plant species throughout their life cycle  <cit> .
in the remaining 10%, which are heteroecious, the aphids alternate between the primary host, usually a woody plant, onto which they migrate during autumn for sexual reproduction, and secondary hosts, usually herbaceous plants, onto which they migrate during spring and develop numerous parthenogenetic colonies.
under specific circumstances, some aphids do not complete such a life cycle and reproduce parthenogenetically throughout the year on herbaceous plants.
this life cycle is termed anholocycly.
approximately 3% of all aphid species are entirely anholocyclic throughout their range, and it has been suggested that 30% to 50% of holocyclic aphid species display a variation in their reproductive mode between and even within populations  <cit> .
only 5% of aphid species are polyphagous  <cit> , comprising crop-pest species that frequently encompass host races or host-specialised populations with a narrow host range.
such patterns have been described for the pea aphid acyrthosiphon pisum  <cit>  and the grain aphid sitobion avenae  <cit> , which are monoecious with a host-plant range restricted to fabaceae and poaceae, respectively.
therefore, ecological specialisation can evolve through feeding behaviour that determines both the host choice - and thus the pool of potential mates - and performance on different host plants  <cit> .
conversely, the evolution of host races in heteroecious  aphid species may represent a paradox because lineages that differ in their preferences for secondary hosts become admixed at mating sites, i.e., on the primary host.
in this case, specialisation could evolve through a reinforcement process whereby selection against hybrids - due to host-associated genetic trade-offs - can drive the increase of assortative mating  <cit> .
host-specialised populations were detected in two heteroecious pest aphid species, the peach-potato aphid myzus persicae  <cit>  and the cotton-melon aphid aphis gossypii  <cit> .
these species have a large spectrum of secondary hosts that belong to numerous plant families, among which cultivated crop species represent a very abundant resource  <cit> .
in m. persicae, host specialisation was evidenced for anholocyclic and holocyclic lineages, and assortative mating seemed to be promoted by differences in the pattern of sexual activity of both the males and mating females  <cit> .
in a. gossypii, the specialised populations consisted of a few anholocyclic lineages  <cit> , suggesting that host-plant specialisation evolved through the loss of the sexual generation in the life cycle  <cit> .
the cotton-melon aphid a. gossypii is distributed worldwide and is found in both tropical and subtropical areas and in temperate zones.
this species has been detected on more than  <dig> plant species  <cit>  and in biotopes affected by human activity  <cit> .
throughout most of its distribution area, including europe, the species is considered to be anholocyclic  <cit> .
however, in areas with rigorous winters in japan, china, korea, india and the united states, a sexual generation has been described on a limited number of primary hosts, including hibiscus syriacus, catalpa bignoïdes and other ligneous plants .
sexual reproduction on h. syriacus has also been reported in italy  <cit> .
genetic studies of a. gossypii populations collected from large apterous colonies sampled in heavily infested crops all over the world revealed a structure in host races on cucurbitaceae, cotton, solanum and pepper plants.
each host race was dominated by a few asexual clones, probably because of host selection, clone competition and pest-management practices  <cit> .
plant-transfer experiments confirmed the existence of host-associated trade-offs in these specialised clones of a. gossypii.
moreover, both the genetic and experimental data suggested that plants of the genus hibiscus may be used as refuge by the specialised clones when the seasonal plant resource is no longer available  <cit> .
the aim of this field study was to investigate how winged spring-migrant populations may contribute to the evolution and maintenance of host specialisation in a. gossypii through host-plant selection.
we also assessed whether host-specialised genotypes corresponded exclusively to anholocyclic lineages in relation to the environmental conditions.
the sampling was completed in areas where the studied agrosystems were dominated or not dominated by cucurbitaceae  crops.
the host-plant selection and reproductive mode were evidenced from a comparison of the genetic structure of the alate spring-migrant populations that alighted on melon crops with the structure of the apterous populations that consisted of the offspring of the migrants.
RESULTS
genetic diversity of the aphis gossypii populations in melon crops
to investigate the genetic diversity of a. gossypii, we sampled alate and apterous aphids from  <dig> to  <dig> in four locations in the southeast  of france, one location in the southwest  and west  of france and one location in the lesser antilles  .
the aphids were genotyped using eight microsatellite markers.
among the  <dig> aphids analysed, we discriminated  <dig> multilocus genotypes .
the majority of the mlgs were present as a single copy, and only four of the mlgs had a frequency higher than 5% of the total number of individual genotypes .
we discriminated  <dig> mlgs among the  <dig> alate aphids;  <dig> to  <dig> of the alate mlgs were characterised in the populations from the different growing areas.
in the alate aphid samples, the allelic diversity ranged from  <dig> to  <dig> alleles per locus, with  <dig> alleles identified across all of the loci.
we discriminated  <dig> mlgs among the  <dig> apterous aphids; three to  <dig> of the apterous mlgs were characterised in the populations from the different growing areas.
in the apterous aphid samples, the allelic diversity ranged from  <dig> to  <dig> alleles per locus, with  <dig> alleles identified across all of the loci.
number of individuals , number of multilocus genotypes , number of multilocus genotypes when the frequencies were higher than 5%  and when the frequencies were lower than 5% , and the number of unique mlgs , in the alate  and apterous  a. gossypii sampled in france and the lesser antilles.
among the  <dig> mlgs discriminated,  <dig> were observed within only the alate samples,  <dig> were observed within only the apterous samples and  <dig> were observed within both of the morphs .
the distribution was different for the numbers of individuals: 8% of the individuals had an mlg that was observed within only the alate samples, 6% of the individuals had an mlg that was observed within only the apterous samples and 86% of individuals had an mlg that was observed within both the alate and apterous samples.
we observed  <dig> mlgs that were present in only one of the studied years, but more than 70% of the individuals had an mlg that was observed in two or more years .
we observed  <dig> mlgs that were present in only one of the areas, but approximately 50% of the individuals had an mlg that was observed in two or more areas .
none of the mlgs that were observed in the la were observed in the other areas.
genotypic diversity of the alate populations vs. the apterous populations
among each morph, the genotypic diversity of the samples was investigated at three levels: within an area, between years and between areas.
three diversity indexes were considered, reflecting the clonal richness, r ; the clonal diversity, d*; and the clonal evenness, v.
the four se locations had similar values for these indexes for both the alate populations and apterous populations .
for example, in the alate populations sampled in  <dig>  the simpson's evenness index  was  <dig>  ±  <dig>  in saint-andiol,  <dig>  ±  <dig>  in montfavet and  <dig>  ±  <dig>  in eyragues .
no strong differences were observed in these three indexes over the different years in each location.
the comparison between the alate and apterous populations showed a significant decrease in the clonal richness of the apterous populations, regardless of the location or the area .
we decided to pool the data from the four locations within the se area for further analyses.
the comparison of the four areas  revealed that the lowest diversity of the alate populations was observed within the populations from the la region and that the highest diversity was observed within the populations from the w region.
an intermediate level of diversity was observed in the populations from the se and sw areas.
the clonal richness  was significantly higher in the sw populations than in the se populations, while the simpson's evenness index  was significantly higher in the populations sampled from the se area than in the populations sampled from the sw area .
in all of the cases, a significant decrease in the clonal richness was observed from the alate to the apterous populations .
genetic structure of the a. gossypii populations
we investigated the aphid genetic structure using amova analyses, considering the effects of the morph , the year of sampling, the growing area and the location .
more than 75% of the variance was distributed within the populations for each of the factors.
the only factor that had a significant effect among the groups was the growing-area factor , which explained 20% of the variance.
the amova analysis confirmed the absence of a significant difference of genetic structure within the populations from the se region.
therefore, we pooled the samples from all of the se locations.
to search for clustering of the mlgs found in this study, we analysed the mlgs along with the  <dig> mlgs previously described and assigned to host races .
using a bayesian structuring program, the likeliest numbers of the genetic clusters were k =  <dig> and k =  <dig> .
a. gossypii populations are grouped according to the aphid morph, year of sampling, growing area or se location.
when k =  <dig>   <dig> mlgs of the  <dig> that were detected in our samples were assigned to cluster a and represented only 5% of the aphids, and  <dig> mlgs were assigned to cluster b, accounting for 95% of the aphids .
the fst value between the clusters was  <dig> .
none of the mlgs that were grouped into cluster a had been described before.
when k =  <dig>  cluster a still consisted of the  <dig> mlgs, but cluster b was divided into three sub-clusters .
cluster x consisted of  <dig> mlgs, representing only 4% of the individuals; cluster x also encompassed mlgs that had been previously collected on solanum spp., pepper, cotton, citrus, strawberry and hibiscus spp.
plants.
cluster y contained  <dig> mlgs, representing  <dig> % of the individuals sampled and  <dig> mlgs that were characteristic of aphids from the cucurbitaceae host race.
therefore, cluster y was called the cucurbit cluster.
cluster z, containing  <dig> mlgs representing 18% of the individuals sampled, contained nm <dig> and nm1-like mlgs that also characterised aphids from cucurbitaceae hosts along with mlgs that were previously identified on cotton  and hibiscus spp.
hosts.
therefore, cluster z was called the nm <dig> cluster.
we analysed the mlgs that were assigned to cluster b and those that were assigned to cluster a separately using the same parameters of the structure program.
the number of genetic sub-clusters that best fitted the data was k =  <dig> for cluster b , and the structure was identical to the structure described above for all of the mlgs when k =  <dig> .
the number of genetic sub-clusters within cluster a was k =  <dig> , and the  <dig> mlgs were quite uniformly distributed in the  <dig> sub-clusters .
we considered the distribution of the alate and apterous individuals according to the existence of the four genetic clusters a, x, y and z, within each of the four geographic areas .
in the se region, 70% of the alate individuals and 80% of the apterous individuals were assigned to cluster y .
in the sw region, 60% of the alate individuals and 80% of the apterous individuals were also assigned to the cucurbit cluster.
in both the se and sw areas,  <dig> to 25% of the individuals, either alate or apterous, belonged to cluster z .
mlgs characteristic of cluster a were present among the alate individuals in the sw region .
in the w region, 5% of the alate and 28% of the apterous individuals had an mlg assigned to the cucurbit cluster, and  <dig> and 30% of the alate and apterous individuals, respectively, had an mlg assigned to the nm <dig> cluster; approximately 50% of the individuals, regardless of the morph, were assigned to cluster a. in the la region, most of the alate and apterous individuals had an mlg assigned to the cucurbit cluster , and these mlgs were not observed in the other growing areas.
in all of the areas,  <dig> to 20% of the alate individuals but almost none of the apterous individuals exhibited an mlg from cluster x .
moreover, in all of the areas, the number of individuals with an mlg assigned to the cucurbit cluster or to the nm <dig> cluster  was significantly higher in the apterous populations than in the alate populations .
notably, although the percentage of apterous individuals belonging to the cucurbitaceae host race was equivalent in the se and sw areas, the diversity among the cucurbit mlgs was considerably higher  in the se region  than in the sw region  .
tracking the genetic fingerprints of sexual reproduction
we examined the genetic variability of the alate populations to track the genetic fingerprints of sexual reproduction.
the dataset was completed by the mlgs that were observed within only the apterous samples.
we calculated three parameters: the hardy-weinberg equilibrium , the linkage disequilibrium  and the fis index.
to reduce the wahlund effect, we split the populations a) according to their cluster membership  and their area   or according to the sub-clusters  for the individuals assigned to cluster a .
number of individuals , number of multilocus genotypes , number of loci in hardy-weinberg equilibrium , number of pairs of loci in linkage disequilibrium , and fis values across all of the loci in the alate a. gossypii populations according to the sampling areas in france and four genetic clusters
number of individuals , number of multilocus genotypes , number of loci in hardy-weinberg equilibrium , number of pairs of loci in linkage disequilibrium , and fis values across all of the loci in the alate a. gossypii populations according to the six genetic sub-clusters within cluster a
the fingerprints of obligate parthenogenesis were departure from the hw equilibrium due to excess heterozygotes, frequent ld and a negative fis value.
these characteristics were observed in the la region subpopulations in cluster x, in the se area subpopulation in cluster z  and in all of the subpopulations of cluster y  except for the w subpopulation.
the fingerprints of panmictic sexual reproduction were frequent hwe, infrequent ld and a positive fis value close to  <dig>  the only subpopulation with these characteristics was the w subpopulation in cluster y. the fingerprints of sexual reproduction by inbreeding were a departure from hwe due to heterozygote deficiency, frequent ld and a positive fis value close to  <dig>  these fingerprints were observed in the se area subpopulation in cluster x, the sw and w area subpopulations in cluster z and the w area subpopulation in cluster a. we observed a fourth situation with a frequent hwe, rare ld and fis values ranging from  <dig>  to  <dig>  for the sw and w area subpopulations assigned to cluster x and the se and sw area subpopulations assigned to cluster a. when cluster a was split into six sub-clusters, hwe was observed for several loci , and the number of significant ld was low .
for the other loci, we observed a heterozygote deficiency that was significant at a level of 5% but not of 1%.
moreover, the fis values across the loci varied from  <dig>  to  <dig> , suggesting the occurrence of sexual reproduction with a varying level of inbreeding in these subpopulations.
discussion
the spring populations of the cotton-melon aphids visiting newly planted melon crops exhibited an unexpectedly high level of genetic diversity that contrasted with the very low diversity characterising the host-specialised populations of this aphid species.
the host races identified in a. gossypii were dominated by a few asexual clones  <cit> , as a result of host selection, clone competition and pest-management practices  <cit> .
because of the use of a novel sampling strategy, the present study illustrated in natura host-plant-selection pressure by showing the great differences in genetic diversity between the spring-migrant populations  and the melon-infesting populations .
moreover, an analysis of the genetic composition of these alate and apterous populations in four geographic regions suggested differences in their life-history strategies and questioned the common assertion that a. gossypii is an anholocyclic species in most of its distribution area, including europe  <cit> .
a strikingly high genetic variability
using the same set of microsatellite markers, previous genetic analyses of more than  <dig> a. gossypii individuals discriminated only  <dig> mlgs.
each of these individuals corresponded to an apterous aphid taken from a different colony on vegetable crops of the cucurbitaceae, solanaceae, malvaceae, rutaceae and rosaceae families in different countries  <cit> .
almost every mlg was present in multiple copies, and 70% of the aphids were distributed in only  <dig> mlgs that were each associated with a particular crop.
the level of genetic diversity detected in the present study was more than ten times higher:  <dig> mlgs were identified from approximately  <dig> aphids sampled exclusively on melon crops in france and the lesser antilles .
this difference may be explained by the time of sampling.
in the present study, we collected the aphids as soon as they migrated into the melon parcels at the very beginning of the crop season, while in the previous studies, apterous aphids were collected in infested crops from colonies that had already undergone the selective filter of plant choice and acceptance.
the majority of the  <dig> mlgs  identified in the present study were present as single copies, and only  <dig> of the mlgs had a frequency greater than 5% .
this distribution was roughly the same in the two studied morphotypes, alate and apterous.
nevertheless, the percentage of individuals that exhibited unique or rare mlgs was only 22% among the alate aphids and 9% among the apterous aphids.
therefore, although a high genetic variability was found, this study showed that the number of clones that recurred over time and space were dominated by only a dozen mlgs .
the analysis of the  <dig> mlgs revealed four genetic clusters .
three of the clusters had already been identified in previous studies, although the genetic diversity observed within each of them was much higher in the present study.
these three clusters contained the mlgs that characterised the aphid clones that were specialised on cultivated plants.
cluster y contained the mlgs that defined the cucurbitaceae host race ; these mlgs were found in more than 70% of the whole dataset.
cluster z contained some mlgs that specialised on cucurbitaceae  and some that specialised on malvaceae; these mlgs were found in 18% of all of the individuals.
cluster x contained mlgs from the solanaceae and cotton host races; these mlgs were found in less than 5% of the aphids.
in contrast, cluster a contained  <dig> mlgs that had never been observed in previous studies.
however, these new mlgs represented less than 5% of the sampled aphids.
the hypothesis that the aphids bearing these mlgs may have originated from a broad array of wild host plants is supported by the huge host-plant spectrum of a. gossypii, which consists of more than  <dig> plant species over  <dig> botanical families  <cit> .
host-plant selection
previous studies demonstrated the existence of host races in the species a. gossypii.
these hyper-specialised clones are found all over the world on the same crops.
because these aphids specialise on annual crops, they are submitted to strong demographic bottlenecks, and they have to produce alate morphs at least twice during the annual cycle, first at the end of the crop season to disperse on refuge plants and again at the beginning of the next season to migrate back to the host crops.
therefore, the settlement of a new colony challenges alate aphids because they may have to make long flights to find a suitable host plant  <cit> .
alate aphids have been shown to use visual and olfactory cues to alight  <cit> .
only a few studies have been conducted on winged-aphid olfactory responses to plant volatiles  <cit> .
studies of flight behaviour suggested that landing frequency increased in traps that released host-plant volatiles compared to empty control traps.
however, it was not demonstrated that winged aphids landed significantly more frequently on their specific host plants.
moreover, the ability to recognise a suitable host plant from a distance, which should be selectively advantageous for a clone, has never been shown in natura.
in the present study, we observed alate individuals that were obviously not specialised on cucurbitaceae on melon crops, i.e., individuals belonging to the clusters a and x. this finding suggests that olfactory cues are unlikely to be decisive in natura for host recognition, which probably occurs only after landing and plant probing.
to test this hypothesis, we compared the diversity and the genetic structure of alate and apterous populations on melon crops at the very beginning of the crop season.
as expected, the clonal richness was significantly lower in the apterous samples than in the alate samples in all of the studied fields .
moreover, this decrease in diversity was mainly due to a loss of the mlgs belonging to the clusters x and a. in the south of france  and in the lesser antilles , these mlgs, which represented up to a quarter of the alate populations, almost completely disappeared in the apterous populations .
this tendency was confirmed in the west of france, where the proportion of the mlgs of the clusters a and x decreased from approximately three quarters in the alate population to one third in the apterous population.
this finding argues in favour of host-plant recognition by alate aphids after probing.
the clones that were not adapted to melons might lay a few larvae before taking off again.
indeed, in a. fabae, studies investigating the relationships between stylet penetration, i.e., probing, and parturition showed that parturition is stimulated by chemical cues that are detected before contact with the phloem  <cit> .
however, these few nymphs were not adapted to feed and therefore to reproduce on melon, as exemplified by the fact that none of them developed colonies.
all of the colonies that were observed corresponded to cucurbitaceae-specialised mlgs from the clusters y or z. our results clearly demonstrate that the melon plant acts as a selective filter against the reproduction of non-specialised individuals.
geographical structure
the genetic diversity of the populations on the melon crops was geographically structured, as revealed by the amova  and pairwise fst analysis .
this structure could be related to the characteristics of the four agrosystems under study.
cucurbitaceae are one of the main vegetable crops cultivated in france based on the surface area used for production, more than  <dig>  ha in mainland france, equally distributed among three regions , and  <dig>  ha in the lesser antilles .
the differentiation was especially strong between the aphid populations from the lesser antilles and those from france.
the genetic diversity in the lesser antilles appeared to be unique and very low; the majority of the individuals exhibited the same mlg belonging to the cucurbit cluster .
obviously, these particular mlgs do not result from a recent introduction of aphids through the high level of commercial exchanges between france and the french west indies.
the introduction might be related to the triangular trade between europe, africa and america in the 17th century, or even earlier, to the colonisation by native americans who travelled from the region of venezuela through the whole chain of the west indies approximately  <dig>  years ago up to the arrival of europeans.
as mentioned by zepeda-paulo et al.
<cit> , the evolutionary processes associated with aphid introductions into new geographical areas have not been extensively studied, but as suggested by normark   <cit> , a temperate aphid lineage transported to the tropics would likely switch to obligate parthenogenesis.
therefore, the very low genetic diversity observed in the melon aphid populations in the lesser antilles could result from the introduction of a small number of clones, of which only those able to adapt to the selective pressures of the new environmental conditions could rapidly multiply and disperse.
when considering the populations from france, a geographical structure of the genetic diversity was detected .
in the se region, we observed that 70% of the individuals were assigned to the cucurbit cluster .
these individuals exhibited  <dig> out of the  <dig> mlgs identified within this cluster, which was the highest clonal richness that we observed within the cucurbit cluster.
in the sw region, 60% of the individuals fell within the cucurbit cluster, and those individuals displayed only  <dig> of the mlgs in that cluster.
in the w region, only 5% of the individuals were assigned to the cucurbit cluster and shared  <dig> mlgs.
these genetic differences could be related to differences in the importance of cucurbitaceae crops, both across seasons and years in these three production areas.
from an historical point of view, the melon crop has been of significance for many decades in the se and sw regions, while in the w region, melons have been grown on a large scale for only ten years.
therefore, cucurbit-adapted mlgs have been favoured and selected for longer in the south than in the western region.
in addition, in the se region, cucurbitaceae crops are grown from january to the end of the fall, either under cover or in open fields, while in the sw region, the plants are restricted to open fields in the spring and summer.
the availability of the cucurbitaceae resource throughout the aphid annual cycle varies considerably among regions, and that variation could explain the differences in the level of genetic diversity of the melon-adapted populations.
in the se area, the time window after which the cucurbitaceae resource is no longer available is very small.
therefore, the cucurbitaceae-specialised populations do not experience a harsh bottleneck during the winter, allowing for the maintenance of slightly divergent mlgs produced by mutation.
the high genetic diversity observed in the cucurbitaceae-specialised populations in the se region probably results from both long-term selection and weak demographic bottlenecks during the winters.
another significant difference lies in the proportion of the alate aphids belonging to cluster a, which was ten times higher in the samples from the sw region  than in the se region .
this difference might be due to a difference in the uniformity of the agricultural environment between the two regions, resulting in a lower abundance of the overwintering colonies of specialised cucurbitaceae clones in the sw region, and therefore, reducing the alate spring population size.
reproductive mode
parthenogenesis is much more common among insect pests than related non-pest insect species, most likely because agricultural environments are stable with an abundance of resources  <cit> .
the same genotypes may be continuously favoured by selection, leading to a selective advantage of some parthenogenetic lineages over a sexual population.
almost  <dig> mlgs were detected in this study, of which 28% were observed more than once over the entire sample of approximately  <dig> aphids.
these repeated mlgs characterised 90% of the individuals, most of which belonged to the cucurbit cluster.
the values of the genetic indices shown in tables  <dig> and  <dig>  strongly suggest that these individuals correspond to obligate parthenogenetic lineages, in accordance with previous genetic studies  <cit> .
sexual-reproduction fingerprints were, however, detected in the cucurbitaceae-specialised population collected in the w region .
although this result needs to be confirmed in a larger sample, it is supported by biological studies that revealed that a few clones sampled from cucurbitaceae hosts were able to induce oviparous female and/or male morphs during conditions of shortened daylight periods  <cit> , carletto .
moreover, as mentioned earlier, the relative heterogeneity of the cucurbitaceae agrosystem in the w region could explain the maintenance of sexual reproduction in cucurbit-specialised populations.
another cucurbitaceae-specialised mlg, nm <dig>  was found in 76% of the individuals from cluster z. this mlg also corresponded to a parthenogenetic lineage.
to date, this mlg has been reported only in aphid samples from cucurbitaceae crops in france  <cit> .
the nm <dig> mlg was largely predominant in the se population where it accounted for the negative value of the fis index .
the frequency of this mlg in the sw and w regions was less than 50%.
in addition to nm <dig>  cluster z contained  <dig> unique mlgs and  <dig> mlgs present more than once .
a genetic analysis comparing the populations from the sw and w regions suggested that sexual reproduction occurred in these populations .
the positive values of fis argued in favour a high level of inbreeding but they could also be artificially inflated by pooling different reproductive units  or because of a mix of asexual and sexual lineages.
cluster x was composed of  <dig> mlgs, of which  <dig> were unique, and the  <dig> mlgs left existed in two or three copies .
as shown in figure  <dig>  there was no redundancy in the mlgs among the different areas, and the sum of the mlgs found in every region was equal to  <dig>  the genetic indices suggested the occurrence of sexual reproduction with inbreeding in populations from the se, sw and w regions.
the  <dig> mlgs belonging to cluster x were grouped with mlgs detected over several years in france and collected from cultivated plant hosts, such as solanum spp., pepper, citrus, strawberry and hibiscus spp.
.
the origin of the new mlgs remains unknown.
cluster a was composed of  <dig> mlgs from  <dig> individuals, indicating that almost every aphid had a different mlg .
most of these genotypes  were detected in populations from the w region  and they were quite uniformly distributed in the six sub-clusters of the a cluster .
the genotypes detected in populations from se and sw regions were also distributed in the six sub-clusters.
the genetic analyses within the populations belonging to the sub-clusters a suggested that sexual reproduction occurred with a high level of inbreeding in these populations .
as mentioned earlier, these aphids are assumed to colonise wild plant species.
the uncultivated plant compartment could be considered as a heterogeneous and disturbed environment, in which the maintenance of sexuality confers advantages, while crop plants correspond to simplified and uniform environments that favoured a few host-adapted asexual lineages  <cit> .
CONCLUSIONS
in conclusion, the genetic survey of the alate aphids revealed that sexual reproduction occurs in french populations of a. gossypii, challenging the commonly shared view that the melon aphid reproduces exclusively by obligate parthenogenesis even in temperate regions  <cit> .
approximately 75% of the alate aphids collected in france were characterised by a dozen multi-copy mlgs and corresponded to parthenogenetic lineages that are highly specialised on cucurbitaceae.
the availability and abundance of the host-plant resource likely selected for the parthenogenetic reproductive mode of plant-adapted genotypes  <cit> .
population genetic structure analyses conducted on two pests of cereal crops, sitobion avenae  <cit>  and rhopalosiphum padi  <cit> , revealed that the populations were distributed in two clusters according to their reproductive mode .
populations of r. padi were not genetically differentiated according to host-plant origin  <cit> , whereas in s. avanae, a much higher proportion of sexual lineages were found on uncultivated host plants  <cit> .
distinct clusters corresponding to populations with a different reproductive mode were also evidenced within the host race nicotianae of m. persicae  <cit> .
our study demonstrated that the differentiation of the a. gossypii populations depended on both the reproductive strategy and host plant.
the bayesian clustering analysis was used to first separate the populations into two groups.
one group consisted of cluster a, with individuals assumed to reproduce sexually and to develop on wild plants, and the other group consisted of cluster b, in which asexual lineages co-existed with a few sexual lineages.
the sub-clustering of b led to the identification of the cluster y, containing mainly parthenogenetic lineages specialised on cucurbitaceae, and the clusters x and z, which likely displayed a mixed reproductive strategy.
the high level of diversity detected in the spring alate aphids was generated by the reproduction of the sexual lineages.
however, this huge diversity was counter-selected by the melon plants, as evidenced by the much lower diversity found in the apterous offspring.
moreover, none of the non-specialised offspring was able to develop a colony.
we cannot exclude that some of the parthenogenetic lineages from the cucurbitaceae host race are androcyclic, producing males and parthenogenetic females or are intermediate, investing into both sexual and parthenogenetic reproduction.
their sexual morphs could migrate to the primary host to mate either with an unrelated mlg or with a cucurbit mlg .
indeed, the cucurbit-specialised clone c <dig>  which was found in high frequency in cucurbitaceae crops all over the world  <cit> , could be intermediate as it was shown to develop sexual morphs in the lab when submitted to induction conditions  <cit> .
new asexual lineages could result from the occasional mating between females from sexual lineages and males produced by asexual lineages that could transmit asexuality genes  <cit> .
therefore, the evolution of a. gossypii as a pest of melon crops would depend on the balance between the gain of genetic diversity and the loss of host-adapted gene combinations resulting from sexual reproduction.
the evolution of parthenogenetic lineages in the cucurbitaceae host race also relies on the selection of new mutants exhibiting a higher fitness in response to pest-management practices and to the clonal amplification of these new variants  <cit> .

 the predictability and abundance of host resources is often assumed to favour ecological specialisation in populations of parasites that have a large spectrum of possible hosts through the disruptive selection of host preference, leading gradually to reproductive isolation and speciation if assortative mating occurs. host-associated biotypes are well documented in plant-feeding insects, and many examples of host races have been described in phytophagous insects in the last ten years. as reviewed by dres and mallet, host races are defined as sets of populations of the same species that use different host taxa in the wild, exhibit host fidelity, coexist in sympatry and correspond to separated genetic clusters among which gene flow still occurs. moreover, these populations usually display a correlation between host choice and mate choice, have a higher fitness on natal than on alternative hosts and produce hybrids that are less fit than the parental forms. host race status has been especially discussed within the aphididae. in the case of aphids, host plants exert considerable selective pressure because the plants constitute the aphids' feeding, mating and oviposition sites. therefore, host specialisation in aphids evolves through the selection of behavioural and chemical mechanisms of host-plant location and recognition as well as through metabolic adaptation to the phloem content of the host plant. how these adaptive traits evolve in an aphid species depends on the complexity of the aphids' annual life cycle. like many parasites, aphids have evolved complex life cycles, with successive generations of individuals often specialised to cope with different ecological conditions. however, aphids are among the few organisms that show plasticity in their reproductive mode as an adaptive response to seasonal changes. the typical aphid annual life cycle, called holocycly or cyclical parthenogenesis, includes a single sexual generation in the autumn composed of males and mating females that lay resting eggs. in the spring, a wingless female hatches from the egg and produces other parthenogenetic females. the number of successive viviparae generations in the spring and summer depends on environmental factors. all of the descendants of a single fundatrix are genotypically identical and constitute a clone, although several phenotypes are expressed. indeed, the viviparae may be alate or apterous in response to environmental conditions, with development of alate forms stimulated in crowded colonies. finally, in the autumn, in response to predictive environmental conditions, the viviparae will produce sexual morphs. approximately 90% of aphid species are monoecious, i.e., they remain on the same host species or closely related plant species throughout their life cycle. in the remaining 10%, which are heteroecious, the aphids alternate between the primary host, usually a woody plant, onto which they migrate during autumn for sexual reproduction, and secondary hosts, usually herbaceous plants, onto which they migrate during spring and develop numerous parthenogenetic colonies. under specific circumstances, some aphids do not complete such a life cycle and reproduce parthenogenetically throughout the year on herbaceous plants. this life cycle is termed anholocycly. approximately 3% of all aphid species are entirely anholocyclic throughout their range, and it has been suggested that 30% to 50% of holocyclic aphid species display a variation in their reproductive mode between and even within populations. only 5% of aphid species are polyphagous, comprising crop-pest species that frequently encompass host races or host-specialised populations with a narrow host range. such patterns have been described for the pea aphid acyrthosiphon pisum and the grain aphid sitobion avenae, which are monoecious with a host-plant range restricted to fabaceae and poaceae, respectively. therefore, ecological specialisation can evolve through feeding behaviour that determines both the host choice - and thus the pool of potential mates - and performance on different host plants. conversely, the evolution of host races in heteroecious aphid species may represent a paradox because lineages that differ in their preferences for secondary hosts become admixed at mating sites, i.e., on the primary host. in this case, specialisation could evolve through a reinforcement process whereby selection against hybrids - due to host-associated genetic trade-offs - can drive the increase of assortative mating. host-specialised populations were detected in two heteroecious pest aphid species, the peach-potato aphid myzus persicae and the cotton-melon aphid aphis gossypii. these species have a large spectrum of secondary hosts that belong to numerous plant families, among which cultivated crop species represent a very abundant resource. in m. persicae, host specialisation was evidenced for anholocyclic and holocyclic lineages, and assortative mating seemed to be promoted by differences in the pattern of sexual activity of both the males and mating females. in a. gossypii, the specialised populations consisted of a few anholocyclic lineages, suggesting that host-plant specialisation evolved through the loss of the sexual generation in the life cycle. the cotton-melon aphid a. gossypii is distributed worldwide and is found in both tropical and subtropical areas and in temperate zones. this species has been detected on more than plant species and in biotopes affected by human activity. throughout most of its distribution area, including europe, the species is considered to be anholocyclic. however, in areas with rigorous winters in japan, china, korea, india and the united states, a sexual generation has been described on a limited number of primary hosts, including hibiscus syriacus, catalpa bignoïdes and other ligneous plants. sexual reproduction on h. syriacus has also been reported in italy. genetic studies of a. gossypii populations collected from large apterous colonies sampled in heavily infested crops all over the world revealed a structure in host races on cucurbitaceae, cotton, solanum and pepper plants. each host race was dominated by a few asexual clones, probably because of host selection, clone competition and pest-management practices. plant-transfer experiments confirmed the existence of host-associated trade-offs in these specialised clones of a. gossypii. moreover, both the genetic and experimental data suggested that plants of the genus hibiscus may be used as refuge by the specialised clones when the seasonal plant resource is no longer available. the aim of this field study was to investigate how winged spring-migrant populations may contribute to the evolution and maintenance of host specialisation in a. gossypii through host-plant selection. we also assessed whether host-specialised genotypes corresponded exclusively to anholocyclic lineages in relation to the environmental conditions. the sampling was completed in areas where the studied agrosystems were dominated or not dominated by cucurbitaceae crops. the host-plant selection and reproductive mode were evidenced from a comparison of the genetic structure of the alate spring-migrant populations that alighted on melon crops with the structure of the apterous populations that consisted of the offspring of the migrants.  genetic diversity of the aphis gossypii populations in melon crops to investigate the genetic diversity of a. gossypii, we sampled alate and apterous aphids from to in four locations in the southeast of france, one location in the southwest and west of france and one location in the lesser antilles. the aphids were genotyped using eight microsatellite markers. among the aphids analysed, we discriminated multilocus genotypes. the majority of the mlgs were present as a single copy, and only four of the mlgs had a frequency higher than 5% of the total number of individual genotypes. we discriminated mlgs among the alate aphids; to of the alate mlgs were characterised in the populations from the different growing areas. in the alate aphid samples, the allelic diversity ranged from to alleles per locus, with alleles identified across all of the loci. we discriminated mlgs among the apterous aphids; three to of the apterous mlgs were characterised in the populations from the different growing areas. in the apterous aphid samples, the allelic diversity ranged from to alleles per locus, with alleles identified across all of the loci. number of individuals, number of multilocus genotypes, number of multilocus genotypes when the frequencies were higher than 5% and when the frequencies were lower than 5%, and the number of unique mlgs, in the alate and apterous a. gossypii sampled in france and the lesser antilles. among the mlgs discriminated, were observed within only the alate samples, were observed within only the apterous samples and were observed within both of the morphs. the distribution was different for the numbers of individuals: 8% of the individuals had an mlg that was observed within only the alate samples, 6% of the individuals had an mlg that was observed within only the apterous samples and 86% of individuals had an mlg that was observed within both the alate and apterous samples. we observed mlgs that were present in only one of the studied years, but more than 70% of the individuals had an mlg that was observed in two or more years. we observed mlgs that were present in only one of the areas, but approximately 50% of the individuals had an mlg that was observed in two or more areas. none of the mlgs that were observed in the la were observed in the other areas. genotypic diversity of the alate populations vs. the apterous populations among each morph, the genotypic diversity of the samples was investigated at three levels: within an area, between years and between areas. three diversity indexes were considered, reflecting the clonal richness, r; the clonal diversity, d*; and the clonal evenness, v. the four se locations had similar values for these indexes for both the alate populations and apterous populations. for example, in the alate populations sampled in the simpson's evenness index was ± in saint-andiol, ± in montfavet and ± in eyragues. no strong differences were observed in these three indexes over the different years in each location. the comparison between the alate and apterous populations showed a significant decrease in the clonal richness of the apterous populations, regardless of the location or the area. we decided to pool the data from the four locations within the se area for further analyses. the comparison of the four areas revealed that the lowest diversity of the alate populations was observed within the populations from the la region and that the highest diversity was observed within the populations from the w region. an intermediate level of diversity was observed in the populations from the se and sw areas. the clonal richness was significantly higher in the sw populations than in the se populations, while the simpson's evenness index was significantly higher in the populations sampled from the se area than in the populations sampled from the sw area. in all of the cases, a significant decrease in the clonal richness was observed from the alate to the apterous populations. genetic structure of the a. gossypii populations we investigated the aphid genetic structure using amova analyses, considering the effects of the morph, the year of sampling, the growing area and the location. more than 75% of the variance was distributed within the populations for each of the factors. the only factor that had a significant effect among the groups was the growing-area factor, which explained 20% of the variance. the amova analysis confirmed the absence of a significant difference of genetic structure within the populations from the se region. therefore, we pooled the samples from all of the se locations. to search for clustering of the mlgs found in this study, we analysed the mlgs along with the mlgs previously described and assigned to host races. using a bayesian structuring program, the likeliest numbers of the genetic clusters were k = and k =. a. gossypii populations are grouped according to the aphid morph, year of sampling, growing area or se location. when k = mlgs of the that were detected in our samples were assigned to cluster a and represented only 5% of the aphids, and mlgs were assigned to cluster b, accounting for 95% of the aphids. the fst value between the clusters was. none of the mlgs that were grouped into cluster a had been described before. when k = cluster a still consisted of the mlgs, but cluster b was divided into three sub-clusters. cluster x consisted of mlgs, representing only 4% of the individuals; cluster x also encompassed mlgs that had been previously collected on solanum spp., pepper, cotton, citrus, strawberry and hibiscus spp. plants. cluster y contained mlgs, representing % of the individuals sampled and mlgs that were characteristic of aphids from the cucurbitaceae host race. therefore, cluster y was called the cucurbit cluster. cluster z, containing mlgs representing 18% of the individuals sampled, contained nm and nm1-like mlgs that also characterised aphids from cucurbitaceae hosts along with mlgs that were previously identified on cotton and hibiscus spp. hosts. therefore, cluster z was called the nm cluster. we analysed the mlgs that were assigned to cluster b and those that were assigned to cluster a separately using the same parameters of the structure program. the number of genetic sub-clusters that best fitted the data was k = for cluster b, and the structure was identical to the structure described above for all of the mlgs when k =. the number of genetic sub-clusters within cluster a was k =, and the mlgs were quite uniformly distributed in the sub-clusters. we considered the distribution of the alate and apterous individuals according to the existence of the four genetic clusters a, x, y and z, within each of the four geographic areas. in the se region, 70% of the alate individuals and 80% of the apterous individuals were assigned to cluster y. in the sw region, 60% of the alate individuals and 80% of the apterous individuals were also assigned to the cucurbit cluster. in both the se and sw areas, to 25% of the individuals, either alate or apterous, belonged to cluster z. mlgs characteristic of cluster a were present among the alate individuals in the sw region. in the w region, 5% of the alate and 28% of the apterous individuals had an mlg assigned to the cucurbit cluster, and and 30% of the alate and apterous individuals, respectively, had an mlg assigned to the nm cluster; approximately 50% of the individuals, regardless of the morph, were assigned to cluster a. in the la region, most of the alate and apterous individuals had an mlg assigned to the cucurbit cluster, and these mlgs were not observed in the other growing areas. in all of the areas, to 20% of the alate individuals but almost none of the apterous individuals exhibited an mlg from cluster x. moreover, in all of the areas, the number of individuals with an mlg assigned to the cucurbit cluster or to the nm cluster was significantly higher in the apterous populations than in the alate populations. notably, although the percentage of apterous individuals belonging to the cucurbitaceae host race was equivalent in the se and sw areas, the diversity among the cucurbit mlgs was considerably higher in the se region than in the sw region. tracking the genetic fingerprints of sexual reproduction we examined the genetic variability of the alate populations to track the genetic fingerprints of sexual reproduction. the dataset was completed by the mlgs that were observed within only the apterous samples. we calculated three parameters: the hardy-weinberg equilibrium, the linkage disequilibrium and the fis index. to reduce the wahlund effect, we split the populations a) according to their cluster membership and their area or according to the sub-clusters for the individuals assigned to cluster a. number of individuals, number of multilocus genotypes, number of loci in hardy-weinberg equilibrium, number of pairs of loci in linkage disequilibrium, and fis values across all of the loci in the alate a. gossypii populations according to the sampling areas in france and four genetic clusters number of individuals, number of multilocus genotypes, number of loci in hardy-weinberg equilibrium, number of pairs of loci in linkage disequilibrium, and fis values across all of the loci in the alate a. gossypii populations according to the six genetic sub-clusters within cluster a the fingerprints of obligate parthenogenesis were departure from the hw equilibrium due to excess heterozygotes, frequent ld and a negative fis value. these characteristics were observed in the la region subpopulations in cluster x, in the se area subpopulation in cluster z and in all of the subpopulations of cluster y except for the w subpopulation. the fingerprints of panmictic sexual reproduction were frequent hwe, infrequent ld and a positive fis value close to the only subpopulation with these characteristics was the w subpopulation in cluster y. the fingerprints of sexual reproduction by inbreeding were a departure from hwe due to heterozygote deficiency, frequent ld and a positive fis value close to these fingerprints were observed in the se area subpopulation in cluster x, the sw and w area subpopulations in cluster z and the w area subpopulation in cluster a. we observed a fourth situation with a frequent hwe, rare ld and fis values ranging from to for the sw and w area subpopulations assigned to cluster x and the se and sw area subpopulations assigned to cluster a. when cluster a was split into six sub-clusters, hwe was observed for several loci, and the number of significant ld was low. for the other loci, we observed a heterozygote deficiency that was significant at a level of 5% but not of 1%. moreover, the fis values across the loci varied from to, suggesting the occurrence of sexual reproduction with a varying level of inbreeding in these subpopulations. discussion the spring populations of the cotton-melon aphids visiting newly planted melon crops exhibited an unexpectedly high level of genetic diversity that contrasted with the very low diversity characterising the host-specialised populations of this aphid species. the host races identified in a. gossypii were dominated by a few asexual clones, as a result of host selection, clone competition and pest-management practices. because of the use of a novel sampling strategy, the present study illustrated in natura host-plant-selection pressure by showing the great differences in genetic diversity between the spring-migrant populations and the melon-infesting populations. moreover, an analysis of the genetic composition of these alate and apterous populations in four geographic regions suggested differences in their life-history strategies and questioned the common assertion that a. gossypii is an anholocyclic species in most of its distribution area, including europe. a strikingly high genetic variability using the same set of microsatellite markers, previous genetic analyses of more than a. gossypii individuals discriminated only mlgs. each of these individuals corresponded to an apterous aphid taken from a different colony on vegetable crops of the cucurbitaceae, solanaceae, malvaceae, rutaceae and rosaceae families in different countries. almost every mlg was present in multiple copies, and 70% of the aphids were distributed in only mlgs that were each associated with a particular crop. the level of genetic diversity detected in the present study was more than ten times higher: mlgs were identified from approximately aphids sampled exclusively on melon crops in france and the lesser antilles. this difference may be explained by the time of sampling. in the present study, we collected the aphids as soon as they migrated into the melon parcels at the very beginning of the crop season, while in the previous studies, apterous aphids were collected in infested crops from colonies that had already undergone the selective filter of plant choice and acceptance. the majority of the mlgs identified in the present study were present as single copies, and only of the mlgs had a frequency greater than 5%. this distribution was roughly the same in the two studied morphotypes, alate and apterous. nevertheless, the percentage of individuals that exhibited unique or rare mlgs was only 22% among the alate aphids and 9% among the apterous aphids. therefore, although a high genetic variability was found, this study showed that the number of clones that recurred over time and space were dominated by only a dozen mlgs. the analysis of the mlgs revealed four genetic clusters. three of the clusters had already been identified in previous studies, although the genetic diversity observed within each of them was much higher in the present study. these three clusters contained the mlgs that characterised the aphid clones that were specialised on cultivated plants. cluster y contained the mlgs that defined the cucurbitaceae host race; these mlgs were found in more than 70% of the whole dataset. cluster z contained some mlgs that specialised on cucurbitaceae and some that specialised on malvaceae; these mlgs were found in 18% of all of the individuals. cluster x contained mlgs from the solanaceae and cotton host races; these mlgs were found in less than 5% of the aphids. in contrast, cluster a contained mlgs that had never been observed in previous studies. however, these new mlgs represented less than 5% of the sampled aphids. the hypothesis that the aphids bearing these mlgs may have originated from a broad array of wild host plants is supported by the huge host-plant spectrum of a. gossypii, which consists of more than plant species over botanical families. host-plant selection previous studies demonstrated the existence of host races in the species a. gossypii. these hyper-specialised clones are found all over the world on the same crops. because these aphids specialise on annual crops, they are submitted to strong demographic bottlenecks, and they have to produce alate morphs at least twice during the annual cycle, first at the end of the crop season to disperse on refuge plants and again at the beginning of the next season to migrate back to the host crops. therefore, the settlement of a new colony challenges alate aphids because they may have to make long flights to find a suitable host plant. alate aphids have been shown to use visual and olfactory cues to alight. only a few studies have been conducted on winged-aphid olfactory responses to plant volatiles. studies of flight behaviour suggested that landing frequency increased in traps that released host-plant volatiles compared to empty control traps. however, it was not demonstrated that winged aphids landed significantly more frequently on their specific host plants. moreover, the ability to recognise a suitable host plant from a distance, which should be selectively advantageous for a clone, has never been shown in natura. in the present study, we observed alate individuals that were obviously not specialised on cucurbitaceae on melon crops, i.e., individuals belonging to the clusters a and x. this finding suggests that olfactory cues are unlikely to be decisive in natura for host recognition, which probably occurs only after landing and plant probing. to test this hypothesis, we compared the diversity and the genetic structure of alate and apterous populations on melon crops at the very beginning of the crop season. as expected, the clonal richness was significantly lower in the apterous samples than in the alate samples in all of the studied fields. moreover, this decrease in diversity was mainly due to a loss of the mlgs belonging to the clusters x and a. in the south of france and in the lesser antilles, these mlgs, which represented up to a quarter of the alate populations, almost completely disappeared in the apterous populations. this tendency was confirmed in the west of france, where the proportion of the mlgs of the clusters a and x decreased from approximately three quarters in the alate population to one third in the apterous population. this finding argues in favour of host-plant recognition by alate aphids after probing. the clones that were not adapted to melons might lay a few larvae before taking off again. indeed, in a. fabae, studies investigating the relationships between stylet penetration, i.e., probing, and parturition showed that parturition is stimulated by chemical cues that are detected before contact with the phloem. however, these few nymphs were not adapted to feed and therefore to reproduce on melon, as exemplified by the fact that none of them developed colonies. all of the colonies that were observed corresponded to cucurbitaceae-specialised mlgs from the clusters y or z. our results clearly demonstrate that the melon plant acts as a selective filter against the reproduction of non-specialised individuals. geographical structure the genetic diversity of the populations on the melon crops was geographically structured, as revealed by the amova and pairwise fst analysis. this structure could be related to the characteristics of the four agrosystems under study. cucurbitaceae are one of the main vegetable crops cultivated in france based on the surface area used for production, more than ha in mainland france, equally distributed among three regions, and ha in the lesser antilles. the differentiation was especially strong between the aphid populations from the lesser antilles and those from france. the genetic diversity in the lesser antilles appeared to be unique and very low; the majority of the individuals exhibited the same mlg belonging to the cucurbit cluster. obviously, these particular mlgs do not result from a recent introduction of aphids through the high level of commercial exchanges between france and the french west indies. the introduction might be related to the triangular trade between europe, africa and america in the 17th century, or even earlier, to the colonisation by native americans who travelled from the region of venezuela through the whole chain of the west indies approximately years ago up to the arrival of europeans. as mentioned by zepeda-paulo et al., the evolutionary processes associated with aphid introductions into new geographical areas have not been extensively studied, but as suggested by normark, a temperate aphid lineage transported to the tropics would likely switch to obligate parthenogenesis. therefore, the very low genetic diversity observed in the melon aphid populations in the lesser antilles could result from the introduction of a small number of clones, of which only those able to adapt to the selective pressures of the new environmental conditions could rapidly multiply and disperse. when considering the populations from france, a geographical structure of the genetic diversity was detected. in the se region, we observed that 70% of the individuals were assigned to the cucurbit cluster. these individuals exhibited out of the mlgs identified within this cluster, which was the highest clonal richness that we observed within the cucurbit cluster. in the sw region, 60% of the individuals fell within the cucurbit cluster, and those individuals displayed only of the mlgs in that cluster. in the w region, only 5% of the individuals were assigned to the cucurbit cluster and shared mlgs. these genetic differences could be related to differences in the importance of cucurbitaceae crops, both across seasons and years in these three production areas. from an historical point of view, the melon crop has been of significance for many decades in the se and sw regions, while in the w region, melons have been grown on a large scale for only ten years. therefore, cucurbit-adapted mlgs have been favoured and selected for longer in the south than in the western region. in addition, in the se region, cucurbitaceae crops are grown from january to the end of the fall, either under cover or in open fields, while in the sw region, the plants are restricted to open fields in the spring and summer. the availability of the cucurbitaceae resource throughout the aphid annual cycle varies considerably among regions, and that variation could explain the differences in the level of genetic diversity of the melon-adapted populations. in the se area, the time window after which the cucurbitaceae resource is no longer available is very small. therefore, the cucurbitaceae-specialised populations do not experience a harsh bottleneck during the winter, allowing for the maintenance of slightly divergent mlgs produced by mutation. the high genetic diversity observed in the cucurbitaceae-specialised populations in the se region probably results from both long-term selection and weak demographic bottlenecks during the winters. another significant difference lies in the proportion of the alate aphids belonging to cluster a, which was ten times higher in the samples from the sw region than in the se region. this difference might be due to a difference in the uniformity of the agricultural environment between the two regions, resulting in a lower abundance of the overwintering colonies of specialised cucurbitaceae clones in the sw region, and therefore, reducing the alate spring population size. reproductive mode parthenogenesis is much more common among insect pests than related non-pest insect species, most likely because agricultural environments are stable with an abundance of resources. the same genotypes may be continuously favoured by selection, leading to a selective advantage of some parthenogenetic lineages over a sexual population. almost mlgs were detected in this study, of which 28% were observed more than once over the entire sample of approximately aphids. these repeated mlgs characterised 90% of the individuals, most of which belonged to the cucurbit cluster. the values of the genetic indices shown in tables and strongly suggest that these individuals correspond to obligate parthenogenetic lineages, in accordance with previous genetic studies. sexual-reproduction fingerprints were, however, detected in the cucurbitaceae-specialised population collected in the w region. although this result needs to be confirmed in a larger sample, it is supported by biological studies that revealed that a few clones sampled from cucurbitaceae hosts were able to induce oviparous female and/or male morphs during conditions of shortened daylight periods, carletto. moreover, as mentioned earlier, the relative heterogeneity of the cucurbitaceae agrosystem in the w region could explain the maintenance of sexual reproduction in cucurbit-specialised populations. another cucurbitaceae-specialised mlg, nm was found in 76% of the individuals from cluster z. this mlg also corresponded to a parthenogenetic lineage. to date, this mlg has been reported only in aphid samples from cucurbitaceae crops in france. the nm mlg was largely predominant in the se population where it accounted for the negative value of the fis index. the frequency of this mlg in the sw and w regions was less than 50%. in addition to nm cluster z contained unique mlgs and mlgs present more than once. a genetic analysis comparing the populations from the sw and w regions suggested that sexual reproduction occurred in these populations. the positive values of fis argued in favour a high level of inbreeding but they could also be artificially inflated by pooling different reproductive units or because of a mix of asexual and sexual lineages. cluster x was composed of mlgs, of which were unique, and the mlgs left existed in two or three copies. as shown in figure there was no redundancy in the mlgs among the different areas, and the sum of the mlgs found in every region was equal to the genetic indices suggested the occurrence of sexual reproduction with inbreeding in populations from the se, sw and w regions. the mlgs belonging to cluster x were grouped with mlgs detected over several years in france and collected from cultivated plant hosts, such as solanum spp., pepper, citrus, strawberry and hibiscus spp.. the origin of the new mlgs remains unknown. cluster a was composed of mlgs from individuals, indicating that almost every aphid had a different mlg. most of these genotypes were detected in populations from the w region and they were quite uniformly distributed in the six sub-clusters of the a cluster. the genotypes detected in populations from se and sw regions were also distributed in the six sub-clusters. the genetic analyses within the populations belonging to the sub-clusters a suggested that sexual reproduction occurred with a high level of inbreeding in these populations. as mentioned earlier, these aphids are assumed to colonise wild plant species. the uncultivated plant compartment could be considered as a heterogeneous and disturbed environment, in which the maintenance of sexuality confers advantages, while crop plants correspond to simplified and uniform environments that favoured a few host-adapted asexual lineages.  in conclusion, the genetic survey of the alate aphids revealed that sexual reproduction occurs in french populations of a. gossypii, challenging the commonly shared view that the melon aphid reproduces exclusively by obligate parthenogenesis even in temperate regions. approximately 75% of the alate aphids collected in france were characterised by a dozen multi-copy mlgs and corresponded to parthenogenetic lineages that are highly specialised on cucurbitaceae. the availability and abundance of the host-plant resource likely selected for the parthenogenetic reproductive mode of plant-adapted genotypes. population genetic structure analyses conducted on two pests of cereal crops, sitobion avenae and rhopalosiphum padi, revealed that the populations were distributed in two clusters according to their reproductive mode. populations of r. padi were not genetically differentiated according to host-plant origin, whereas in s. avanae, a much higher proportion of sexual lineages were found on uncultivated host plants. distinct clusters corresponding to populations with a different reproductive mode were also evidenced within the host race nicotianae of m. persicae. our study demonstrated that the differentiation of the a. gossypii populations depended on both the reproductive strategy and host plant. the bayesian clustering analysis was used to first separate the populations into two groups. one group consisted of cluster a, with individuals assumed to reproduce sexually and to develop on wild plants, and the other group consisted of cluster b, in which asexual lineages co-existed with a few sexual lineages. the sub-clustering of b led to the identification of the cluster y, containing mainly parthenogenetic lineages specialised on cucurbitaceae, and the clusters x and z, which likely displayed a mixed reproductive strategy. the high level of diversity detected in the spring alate aphids was generated by the reproduction of the sexual lineages. however, this huge diversity was counter-selected by the melon plants, as evidenced by the much lower diversity found in the apterous offspring. moreover, none of the non-specialised offspring was able to develop a colony. we cannot exclude that some of the parthenogenetic lineages from the cucurbitaceae host race are androcyclic, producing males and parthenogenetic females or are intermediate, investing into both sexual and parthenogenetic reproduction. their sexual morphs could migrate to the primary host to mate either with an unrelated mlg or with a cucurbit mlg. indeed, the cucurbit-specialised clone c which was found in high frequency in cucurbitaceae crops all over the world, could be intermediate as it was shown to develop sexual morphs in the lab when submitted to induction conditions. new asexual lineages could result from the occasional mating between females from sexual lineages and males produced by asexual lineages that could transmit asexuality genes. therefore, the evolution of a. gossypii as a pest of melon crops would depend on the balance between the gain of genetic diversity and the loss of host-adapted gene combinations resulting from sexual reproduction. the evolution of parthenogenetic lineages in the cucurbitaceae host race also relies on the selection of new mutants exhibiting a higher fitness in response to pest-management practices and to the clonal amplification of these new variants. 
the purpose of this field study was to determine how winged spring-migrant populations contribute to the evolution and maintenance of host specialisation in aphis gossypii through host-plant choice and acceptance. how these adaptive traits evolve in an aphid species depends on the complexity of the annual life cycle of that species. moreover, an analysis of the genetic composition of these alate and apterous populations in four geographic regions suggested differences in life-history strategies, such as host choice and reproductive mode, and questioned the common assertion that a. gossypii is an anholocyclic species throughout its distribution area, including europe. the spring populations of cotton-melon aphids visiting newly planted melon crops exhibited an unexpectedly high level of genetic diversity that contrasted with the very low diversity characterising the host-specialised populations of this aphid species. this study illustrated in natura host-plant-selection pressure by showing the great differences in genetic diversity between the spring-migrant populations and the melon-infesting populations. our results clearly demonstrate that the melon plant acts as a selective filter against the reproduction of non-specialised individuals. therefore, host specialisation in aphids evolves through selection of the behavioural and chemical mechanisms of host-plant location and recognition, and through metabolic adaptation to the phloem content of the host plant. the agroecosystem structure and history of the four studied regions may have partially shaped the genetic structure of the spring-migrant populations of a. gossypii. we also assessed whether host-specialised genotypes corresponded exclusively to anholocyclic lineages regardless of the environmental conditions. cucurbitaceae-specialised genotypes corresponded exclusively to anholocyclic lineages, regardless of the environmental conditions. however, some genotypes that were genetically close to the host-specialised genotypes and some genotypes that probably originated from wild plants had never been previously sampled; both were holocylic. host plants exert considerable selective pressure on aphids because the plants constitute their feeding, mating and oviposition sites. we showed that olfactory cues are unlikely to be decisive in natura for host recognition by spring-migrant aphid populations that are not specialised on cucurbitaceae.
 host plants exert considerable selective pressure on aphids because the plants constitute their feeding, mating and oviposition sites. therefore, host specialisation in aphids evolves through selection of the behavioural and chemical mechanisms of host-plant location and recognition, and through metabolic adaptation to the phloem content of the host plant. how these adaptive traits evolve in an aphid species depends on the complexity of the annual life cycle of that species. the purpose of this field study was to determine how winged spring-migrant populations contribute to the evolution and maintenance of host specialisation in aphis gossypii through host-plant choice and acceptance. we also assessed whether host-specialised genotypes corresponded exclusively to anholocyclic lineages regardless of the environmental conditions.  the spring populations of cotton-melon aphids visiting newly planted melon crops exhibited an unexpectedly high level of genetic diversity that contrasted with the very low diversity characterising the host-specialised populations of this aphid species. this study illustrated in natura host-plant-selection pressure by showing the great differences in genetic diversity between the spring-migrant populations and the melon-infesting populations. moreover, an analysis of the genetic composition of these alate and apterous populations in four geographic regions suggested differences in life-history strategies, such as host choice and reproductive mode, and questioned the common assertion that a. gossypii is an anholocyclic species throughout its distribution area, including europe.  our results clearly demonstrate that the melon plant acts as a selective filter against the reproduction of non-specialised individuals. we showed that olfactory cues are unlikely to be decisive in natura for host recognition by spring-migrant aphid populations that are not specialised on cucurbitaceae. the agroecosystem structure and history of the four studied regions may have partially shaped the genetic structure of the spring-migrant populations of a. gossypii. cucurbitaceae-specialised genotypes corresponded exclusively to anholocyclic lineages, regardless of the environmental conditions. however, some genotypes that were genetically close to the host-specialised genotypes and some genotypes that probably originated from wild plants had never been previously sampled; both were holocylic. 