to test whether grass species represent habitats of variable quality for thrips , we measured reproductive success on different grasses for females of the sexual species a. elegans , a sexual a. rufus lineage , and asexual a. rufus .
we grew five grass species in 50 ml falcon tubes , with 2466 replicates per aptinothrips species and grass species : molinia caerulea ( moorgrass ) , bromus hordeaceus ( brome ) , festuca rubra ( fescue ) , triticum aestivum ( wheat ) , and secale cereale ( rye ) .
one fieldcollected , adult female was added to each tube and allowed to reproduce for 24 days at 23c and 10:14 night / day cycles .
reproductive success of each female was then determined by counting all offspring present in a tube and these offspring numbers were compared between treatments using a glm with a poissondistributed error and an empirically estimated scale parameter .
whether thrips species , grass species , and their interaction explained a significant amount of deviance in the glm was evaluated using a test . based on the results from the previous experiment ,
we chose three easytogrow grass species that varied in their quality for thrips : wheat , fescue , and rye . with these grasses , we generated four different types of mesocosms : one with a heterogeneous habitat ( the three species mixed in equal proportions ) and three with homogeneous habitats consisting only of one grass species .
mesocosms were generated in transparent , cylindrical cages ( 15 cm diameter , 25 cm high ) , randomly placed on watering tables in a greenhouse under a light regime of l / d = 14:10 superimposed over a natural light regime .
these mesocosms were used to test whether habitat heterogeneity affects the relative success of competing sexual and asexual lineages in two different assays , each running over approximately 12 thrips generations ( thrips have overlapping generations ) . in the first assay , we used two aptinothrips
lineages that frequently cooccur in natural populations but belong to different species : sexual a. elegans and asexual a. rufus .
females for this assay were sampled in a hay meadow in switzerland ( n46.657906 , w6.533641 ) , where the two species are sympatric with a ratio of 49:51 ( n > 3000 ) and that is mainly composed of the grasses arrhenatherum elatius , dactylis glomerata , and bromus hordeaceus .
aptinothrips elegans and a. rufus females were sorted by species under a dissecting scope at 40 magnification , after co2 anesthetization .
the two species are very similar but can be distinguished by the shape of antennal segments and position of hairs on the abdomen ( palmer 1975 ) . in the second competition assay , we used a. rufus sexual and asexual strains from laboratory stock populations to account for eventual species differences in the first assay .
the stock populations have been maintained in the laboratory for over 3 years on wheat ( > 30 generations ) . under the laboratory conditions used ( 23c and 10:14 night / day cycles ) , the sexual a. rufus populations were characterized by an approximately constant proportion of 0.32 males .
for each competition assay , we prepared a total of 160 mesocosms ( 40 heterogeneous ones and 40 homogeneous ones for each of the three grass species ) .
six asexual and six mated sexual females were then introduced into each of the 160 mesocosms .
this is below the carrying capacity of these mesocosms , which means that competition occurred after a colonisation episode .
furthermore , it is important to note that given the very high clonal diversities in asexual a. rufus populations , the six asexual females per mesocosm most likely comprised several different clones . to quantify the frequency of sexual and asexual lineages present ,
we removed randomly without replacement ten heterogeneous mesocosms and 30 homogeneous ones ( 10 per grass species ) at three successive intervals of 40  2 days ( approximately two thrips generations ) . thrips in the remaining ten heterogeneous mesocosms and 30 homogeneous ones were counted after a final interval of 120 days ( 240  6 days from the start of the experiment ) . in the sexual a. elegans versus asexual a. rufus assay , we determined the frequency of sexuals and asexuals by counting the number of adult females of each species . in the sexual versus asexual a. rufus assay , we counted the proportion of males among adults because sexual and asexual a. rufus females can not be distinguished morphologically . because the mesocosm conditions differed from the conditions in the field ( notably because of different grasses ) , we may have selected for adaptation to new host plants . to test for evidence of adaptation , we collected individuals alive from the 40 mesocosms of the first interval of competing a. elegans versus asexual a. rufus .
individual a. elegans and a. rufus females were then transferred to 50 ml tubes with the same habitat as the mesocosm they were extracted from and allowed to reproduce for 24 days as described above . reproductive output for each thrips species
was then compared to the reproductive output under the same conditions prior to the experiment .
thrips from all the remaining mesocosms were counted after extraction into alcohol ( adults and larvae ) via heat gradients on berlese funnels . to test the effect of habitat heterogeneity on the success of sexual and asexual grassthrips in natural populations , we took advantage of experimental hay meadows in the inner european alps ( valais , sw switzerland ) , where sexual a. elegans cooccur with asexual a. rufus , a. stylifer , and a. karnyi .
these meadows have been extensively managed for a long time before the start of the experiment . within each meadow ,
grass diversity was manipulated via six different irrigation and fertilization treatments ( andrey et  al .
these meadows are strongly water and nutrient limited and are occupied by xerophilic species , such that low or intermediate management intensities were expected to increase plant diversity in general and grass diversity in particular ( peter et  al .
precise details on site selection and experimental setups are presented elsewhere ( andrey et  al . 2014 ) .
briefly , six circular plots ( 20 m diameter ) were delineated at each of 11 study sites ( 66 different plots ; see fig .
in the first treatment (  standard  ) , hay meadows were managed as prior to the start of the experiment with grass cut once per year and no other manipulation .
the second treatment consisted in watering the plots once per week , the third treatment in fertilizing them .
for the forth to sixth treatments , plots were simultaneously watered and fertilized in an increasing intensity ( see table s1 ) .
treatments were initiated in spring 2010 and continued over five years , inducing a gradual transition from the original to new plant communities ( andrey et  al .
2014 ) over 4050 thrips generations . to estimate the final grass diversity in each plot , we used coverage data for grass species from exhaustive plant inventories in rectangles of 2  4 m within each plot , conducted in may  june 2014 ; thrips were collected in june 2014 ( see supporting information ) .
wiener diversity indices from grass species coverage in experimental plots using the r package vegan ( oksanen et  al .
we tested the effects of the treatments on habitat heterogeneity ( i.e. , grass diversity ) using a linear mixed model ( lmm ) .
the effect of the treatment and of habitat heterogeneity on the weighted counts of sexual aptinothrips was tested using generalized lmms ( glmms ) implemented in the r package lme4 ( bates et  al .
2015 ) . because the residuals in some of our models were overdispersed , we corroborated each glmm result with an anova using the proportion of sexual aptinothrips rather than their weighted counts ( see supporting information ) . because our analyses revealed a negative effect of grass diversity on the success of sexual thrips ( see results ) , we investigated whether specific types of grass communities could favor sexual aptinothrips .
we first generated a grass community distance dendrogram using braycurtis distances , as implemented in the r package vegan .
we then used this dendrogram to test whether similar communities were characterized by similar sexual thrips abundances by calculating pagel 's  with the r package phytools ( revell 2012 ) .
pagel 's  is typically used to assess whether phylogenetically related species share similar values of a quantitative trait . in our analyses
pagel 's  thus indicates whether similar grass communities comprise more similar abundances of sexual thrips than expected by chance . for visual evaluation
, we further performed a principal component analysis ( pca ) as implemented in the r package vegan , using a hellinger transformation for scaling grass cover data ( legendre and gallagher 2001 ) .
to test whether grass species represent habitats of variable quality for thrips , we measured reproductive success on different grasses for females of the sexual species a. elegans , a sexual a. rufus lineage , and asexual a. rufus .
we grew five grass species in 50 ml falcon tubes , with 2466 replicates per aptinothrips species and grass species : molinia caerulea ( moorgrass ) , bromus hordeaceus ( brome ) , festuca rubra ( fescue ) , triticum aestivum ( wheat ) , and secale cereale ( rye ) .
one fieldcollected , adult female was added to each tube and allowed to reproduce for 24 days at 23c and 10:14 night / day cycles .
reproductive success of each female was then determined by counting all offspring present in a tube and these offspring numbers were compared between treatments using a glm with a poissondistributed error and an empirically estimated scale parameter .
whether thrips species , grass species , and their interaction explained a significant amount of deviance in the glm was evaluated using a test .
based on the results from the previous experiment , we chose three easytogrow grass species that varied in their quality for thrips : wheat , fescue , and rye . with these grasses , we generated four different types of mesocosms : one with a heterogeneous habitat ( the three species mixed in equal proportions ) and three with homogeneous habitats consisting only of one grass species .
mesocosms were generated in transparent , cylindrical cages ( 15 cm diameter , 25 cm high ) , randomly placed on watering tables in a greenhouse under a light regime of l / d = 14:10 superimposed over a natural light regime .
these mesocosms were used to test whether habitat heterogeneity affects the relative success of competing sexual and asexual lineages in two different assays , each running over approximately 12 thrips generations ( thrips have overlapping generations ) . in the first assay ,
we used two aptinothrips lineages that frequently cooccur in natural populations but belong to different species : sexual a. elegans and asexual a. rufus .
females for this assay were sampled in a hay meadow in switzerland ( n46.657906 , w6.533641 ) , where the two species are sympatric with a ratio of 49:51 ( n > 3000 ) and that is mainly composed of the grasses arrhenatherum elatius , dactylis glomerata , and bromus hordeaceus .
aptinothrips elegans and a. rufus females were sorted by species under a dissecting scope at 40 magnification , after co2 anesthetization .
the two species are very similar but can be distinguished by the shape of antennal segments and position of hairs on the abdomen ( palmer 1975 ) . in the second competition assay
, we used a. rufus sexual and asexual strains from laboratory stock populations to account for eventual species differences in the first assay .
the stock populations have been maintained in the laboratory for over 3 years on wheat ( > 30 generations ) . under the laboratory conditions used ( 23c and 10:14 night / day cycles ) ,
the sexual a. rufus populations were characterized by an approximately constant proportion of 0.32 males . for each competition assay , we prepared a total of 160 mesocosms ( 40 heterogeneous ones and 40 homogeneous ones for each of the three grass species ) .
six asexual and six mated sexual females were then introduced into each of the 160 mesocosms .
this is below the carrying capacity of these mesocosms , which means that competition occurred after a colonisation episode .
furthermore , it is important to note that given the very high clonal diversities in asexual a. rufus populations , the six asexual females per mesocosm most likely comprised several different clones . to quantify the frequency of sexual and asexual lineages present , we removed randomly without replacement ten heterogeneous mesocosms and 30 homogeneous ones ( 10 per grass species ) at three successive intervals of 40  2 days ( approximately two thrips generations ) .
thrips in the remaining ten heterogeneous mesocosms and 30 homogeneous ones were counted after a final interval of 120 days ( 240  6 days from the start of the experiment ) . in the sexual a. elegans versus asexual a. rufus assay , we determined the frequency of sexuals and asexuals by counting the number of adult females of each species . in the sexual versus asexual a. rufus assay
, we counted the proportion of males among adults because sexual and asexual a. rufus females can not be distinguished morphologically . because the mesocosm conditions differed from the conditions in the field ( notably because of different grasses ) , we may have selected for adaptation to new host plants . to test for evidence of adaptation , we collected individuals alive from the 40 mesocosms of the first interval of competing a. elegans versus asexual a. rufus .
individual a. elegans and a. rufus females were then transferred to 50 ml tubes with the same habitat as the mesocosm they were extracted from and allowed to reproduce for 24 days as described above . reproductive output for each thrips species
was then compared to the reproductive output under the same conditions prior to the experiment .
thrips from all the remaining mesocosms were counted after extraction into alcohol ( adults and larvae ) via heat gradients on berlese funnels .
to test the effect of habitat heterogeneity on the success of sexual and asexual grassthrips in natural populations , we took advantage of experimental hay meadows in the inner european alps ( valais , sw switzerland ) , where sexual a. elegans cooccur with asexual a. rufus , a. stylifer , and a. karnyi .
these meadows have been extensively managed for a long time before the start of the experiment . within each meadow ,
grass diversity was manipulated via six different irrigation and fertilization treatments ( andrey et  al .
these meadows are strongly water and nutrient limited and are occupied by xerophilic species , such that low or intermediate management intensities were expected to increase plant diversity in general and grass diversity in particular ( peter et  al .
precise details on site selection and experimental setups are presented elsewhere ( andrey et  al . 2014 ) .
briefly , six circular plots ( 20 m diameter ) were delineated at each of 11 study sites ( 66 different plots ; see fig .
in the first treatment (  standard  ) , hay meadows were managed as prior to the start of the experiment with grass cut once per year and no other manipulation .
the second treatment consisted in watering the plots once per week , the third treatment in fertilizing them . for
the forth to sixth treatments , plots were simultaneously watered and fertilized in an increasing intensity ( see table s1 ) .
treatments were initiated in spring 2010 and continued over five years , inducing a gradual transition from the original to new plant communities ( andrey et  al .
2014 ) over 4050 thrips generations . to estimate the final grass diversity in each plot , we used coverage data for grass species from exhaustive plant inventories in rectangles of 2  4 m within each plot , conducted in may  june 2014 ; thrips were collected in june 2014 ( see supporting information ) .
we computed shannon  wiener diversity indices from grass species coverage in experimental plots using the r package vegan ( oksanen et  al .
we tested the effects of the treatments on habitat heterogeneity ( i.e. , grass diversity ) using a linear mixed model ( lmm ) .
the effect of the treatment and of habitat heterogeneity on the weighted counts of sexual aptinothrips was tested using generalized lmms ( glmms ) implemented in the r package lme4 ( bates et  al .
because the residuals in some of our models were overdispersed , we corroborated each glmm result with an anova using the proportion of sexual aptinothrips rather than their weighted counts ( see supporting information ) . because our analyses revealed a negative effect of grass diversity on the success of sexual thrips ( see results ) , we investigated whether specific types of grass communities could favor sexual aptinothrips .
we first generated a grass community distance dendrogram using braycurtis distances , as implemented in the r package vegan .
we then used this dendrogram to test whether similar communities were characterized by similar sexual thrips abundances by calculating pagel 's  with the r package phytools ( revell 2012 ) .
pagel 's  is typically used to assess whether phylogenetically related species share similar values of a quantitative trait . in our analyses
pagel 's  thus indicates whether similar grass communities comprise more similar abundances of sexual thrips than expected by chance . for visual evaluation
, we further performed a principal component analysis ( pca ) as implemented in the r package vegan , using a hellinger transformation for scaling grass cover data ( legendre and gallagher 2001 ) .
grass species vary significantly in quality for thrips as revealed by the number of progeny females produced on the five species wheat , brome , fescue , rye , and moorgrass ( thrips species effect : df = 2 , p < 0.0001 ; grass species effect : df = 4 , p < 0.0001 ; interaction : df = 8 , p = 0.11 ; fig .
females of the three tested thrips species / lineages were unable to reproduce on moorgrass and we found up to sevenfold differences in reproductive output among the four remaining grasses . given this broad quality variation among grasses , grass diversity is a suitable proxy for the level of habitat heterogeneity in a thrips population .
grassthrips females produce different numbers of offspring depending on the grass species they breed on , indicating that different grasses correspond to habitats of different quality .
n = 2466 females per thrips and grass species ; residual offspring numbers after accounting for variation among thrips species are presented .
we then used three of the five grass species to generate heterogeneous ( three species mixed ) and homogeneous ( one species ) mesocosms in which we followed the relative proportion of sexual and asexual thrips over time . within 80 days of the experiment ( approximately four thrips generations ) ,
each mesocosm comprised close to 300 individuals , a number that remained approximately constant until the 12th generation ( fig .
2 ) . independently of the mesocosm type ( heterogeneous or different homogeneous ones ) and the sexual lineage we used ( sexual a. rufus or a. elegans ) , sexual females were always completely replaced by asexual ones after 12 generations ( fig .
we therefore found no positive effect of habitat heterogeneity on the maintenance of sexual thrips .
( a ) proportion of sexual a. elegans females and ( b ) proportion of a. rufus males in mesocosms across different generations ; the dotted line in panel ( b ) indicates the proportion of males in laboratory populations of sexual a. rufus ( i.e. , the proportion of males corresponding to the fixation of sexual lineages ) .
generation 0 started with six asexual females and six mated sexual females in all mesocosms , n
the proportions of males , respectively sexual and asexual lineages , were not determined for generations 810 .
different colours in the homogeneous mesocosms indicate different grass species used in the homogeneous mesocosms ( dark green : fescue , intermediate : wheat , light green : rye ) .
right axis : population size , as measured by the total number of aptinothrips in the mesocosm .
the sexual thrips species a. elegans frequently cooccurs with asexual a. rufus , with individuals used in the mesocosm experiments stemming from a location with approximately equal frequencies of both species . nevertheless , the rate of population increase of asexual a. rufus was much higher than of sexual a. elegans in all mesocosms , leading to complete disappearance of a. elegans after 412 generations .
after approximately six generations , a. elegans was completely extinct in 10 homogeneous rye mesocosms , and in seven respectively eight of the homogeneous fescue and wheat ones ( without a single a. elegans individual among over 7000 thrips ) .
aptinothrips elegans females were still present in the 15 remaining mesocosms , but at a very low frequency in 14 of the 15 ( one to five individuals among several hundred asexual females ) .
only one homogeneous wheat mesocosm still comprised 12 a. elegans females ( among 369 asexual females ) after six generations .
however , these residual individuals from sexual lineages were completely replaced by asexual females by generation 12 , at which point we did not find a single a. elegans in any of the 40 mescosms ( combined over 11,000 females ) .
thus , although the final extinction of sexuals was somewhat delayed in heterogeneous mesocosms ( after six generations : complete extinction in 25 out of 30 homogeneous habitats vs. 0 out of 10 heterogeneous ones ; fisher 's exact test p = 0.01 ) , sexuals were unable to persist when in competition with asexuals , whatever the conditions .
part of the success of asexual a. rufus , when in competition with a. elegans under the used rearing conditions stems from an approximately twofold fecundity advantage of a. rufus that is unrelated to reproductive mode but characterizes different species ( see below ; fig .
however , we also found complete extinction of the sexual lineages after 12 generations for the mesocosms where asexual a. rufus were in competition with sexual a. rufus lineages .
after generation 12 , we found not a single male among over 12,000 females ( fig .
in contrast to the first competition experiment where sexual lineages were largely extinct by generation six , the proportion of males present at generation six ( 00.11 per mesocosm , fig .
consequently , we did not find any difference between homogeneous and heterogeneous mesocosms in the residual persistence of sexuals ( number of mesocosms without males after six generations : 4 out of 30 homogeneous ones vs. 1 out of 10 heterogeneous ones , p > 0.99 ) .
fitness ( number of offspring produced in 24 days ) of asexual a. rufus ( a ) and of sexual a. elegans females ( b ) prior to the experiment ( generation 0 ) and after approximately two and 12 ( for a. rufus ) generations of adaptation to the new conditions .
asexual a. rufus and sexual a. elegans used in the first mesocosm experiment were collected from a meadow comprising the grass species arrhenatherum elatius , dactylis glomerata , and bromus hordeaceus and without any of the three species used in the mesocosm experiment .
we therefore expected to observe some adaptation to the new grass species in the mesocosms .
we tested for evidence of adaptation by measuring reproductive output of individual females before the start of the experiment and after two generations in the mesocosms .
for asexual a. rufus in homogeneous mesocosms , we found , as expected , that the average number of offspring females produced within 24 days increased after ca .
two generations of selection ( grass species : df = 2 , p < 0.0001 ; generation : df = 2 , p < 0.0001 ; fig .
the reproductive output increase after two generations was accompanied by a significant decrease in variation among females in offspring number on wheat and fescue ( onetailed ansari  bradley tests , ab = 1478.5 , p = 0.02 and ab = 1451 , p = 0.05 , respectively ) as expected if the frequency of genotypes poorly performing on a given plant decreased as a consequence of directional selection .
a significant variation decrease was not observed on rye ( ab = 1749 , p = 0.35 ) .
in contrast to the improved reproductive output after two generations of selection in homogeneous habitats , we found no evidence for adaptation in the heterogeneous mesocosms ( generation effect : df = 2 , p = 0.98 , fig .
3a ) as expected if adaptation to the new conditions was hampered by the presence of multiple plants ( i.e. , spatially fluctuating selection ) .
however , in contrast to asexual a. rufus , we found no evidence for adaptation to mesocosm conditions in sexual a. elegans .
the reproductive output of sexual females did not improve significantly after two generations in the mesocosms , neither for homogeneous ( plant effect : df = 2 , p < 0.0001 ; generation effect : df = 1 , p = 0.166 ) or heterogeneous habitats ( generation effect : df = 1 , p = 0.293 ; fig .
no measurement of fitness was conducted past the second generation in a. elegans because we were unable to collect enough a. elegans females from the mesocosms at that point . in summary ,
our results from the mesocosm experiments show that sexual lineages ( a. elegans and a. rufus ) are unable to persist when in competition with asexual a. rufus lineages in homogeneous habitats consisting of a single grass species and in heterogeneous habitats consisting of a mix of three species .
however , although a mix of three grasses corresponds to a more complex habitat for thrips than a single grass species , a mix of three grasses may not be representative of natural populations , where typically more than three different grass species cooccur within a field .
therefore , the most relevant insights into the effect of habitat heterogeneity on the success of sexual thrips stem from diversity variation in natural populations . in the natural setting represented by 66 plots in 11 hay meadows ,
29 different grass species were inventoried , with an average of 8.6 ( range 213 ) species per plot .
the fertilization and irrigation treatments applied to the plots significantly changed grass diversity ( lmm with meadow location as random effect : p = 0.04 ; fig .
the treatments also influenced the proportion of sexual thrips in each plot ( glmm ,  = 26.62 , df = 5 , p = 0.00007 ; lmm : p = 0.0002 ; fig .
contrary to the predictions , the frequency of sexual thrips was negatively affected by treatments increasing grass diversity ( i.e. , the level of habitat heterogeneity ) .
this effect was most obvious when directly using grass diversity independently of fertilization and irrigation treatments , with a significant negative effect of the grass species diversity on the proportion of sexuals ( glmm  = 4.96 , df = 1 , p = 0.026 ; lmm : p = 0.011 , see also supporting information ) .
this negative correlation was also observed within standard plots only ( r 
  = 0.46 , p = 0.02 ) , indicating it was not solely driven by the shifts in grass communities following the irrigation and fertilization treatments .
( a ) habitat heterogeneity ( measured as grass diversity indices ) among plots in experimental hay meadows exposed to fertilization and irrigation treatments . combined irrigating and
fertilizing treatments ( irrigated + fertilized ) were applied at three different intensity levels ( low , intermediate , and high ) .
( b ) proportion of individuals from sexual aptinothrips species with respect to plot treatment .
( c ) correlation between grass diversity and the proportion of sexual thrips independently of plot treatments . given the unexpected finding that the frequency of sexual thrips was negatively affected by increased grass diversities in the experimental plots
, we examined the plots in which a high frequency of sexual thrips was maintained .
specifically , we examined whether plots with a high frequency of sexual thrips comprised specific grass communities , which would indicate that sexual thrips only occur in a specific type of niche . to this end , we grouped plots by their grass species composition in a dendrogram , using bray 
curtis distances to quantify the differences between plots . mapping the proportions of sexual thrips on this dendrogram revealed that plots with similar frequencies of sexual thrips were also characterized by similar grass communities ( pagel 's  = 0.55 , p < 0.0001 ; fig .
similarly , plots with high frequencies of sexual thrips tended to be clustered when the grass communities of the 66 plots were displayed along the first two axes obtained from a pca ( fig .
( a ) similarities of grass communities in the experimental plots displayed as a dendrogram .
plot labels are shaded proportionally to the frequency of sexual aptinothrips species . in combination ,
the grass community analyses support the idea that sexual thrips occur in a specific type of niche . as a first step to characterize this niche
, we identified the grass species that were the most strongly correlated with sexual thrips abundances .
only the grass species bromus erectus was consistently associated with sexual a. elegans ( see supplementary methods and results and tables s2 and s3 in the supporting information ) , with a strong correlation between the coverage of b. erectus in the plots and the proportion of sexual thrips ( r 
  = 0.58 , p < 0.0001 ) .
explaining the overwhelming success of sexual reproduction has been a challenge in evolutionary biology for decades .
recent theoretical approaches have identified plausible mechanisms that could provide benefits to sex and therefore help explain its maintenance .
for some of these mechanisms , experimental evolution has provided a  proof of principle  that the predicted mechanisms can work in real organisms given the appropriate conditions ( colegrave 2002 ; goddard et  al .
however , these conditions may not be realized in natural populations such that these approaches provide little insights into the maintenance of sex in nature .
indeed , it is impossible to know whether any benefit to sex detected under artificial conditions ( that may include controlled migration rates , specific population densities or sizes ) would outweigh its immediate costs expressed under natural conditions . to fill this gap ,
we tested whether one of the classical conditions predicted to provide benefits to sex , habitat heterogeneity , can contribute to the maintenance of sex under natural conditions . by analyzing the frequency of sexual and asexual grassthrips in experimental hay meadows
, we found that an increase of habitat heterogeneity ( measured via grass diversity ) results in a decrease of sexual grassthrips and favors asexual ones .
the frequency of sexuals and grass diversity was also negatively correlated in unmanipulated grass communities . in seminatural mesocosm conditions
, habitat heterogeneity did not affect the outcome of competition between sexual and asexual grassthrips ; asexual grassthrips completely replaced their sexual relatives within 412 generations . in combination , these experiments show that the potential benefits that sex confers under habitat heterogeneity are negligible compared to its costs when sexual lineages compete with asexuals under natural or seminatural conditions .
experimental evolution studies have also shown that sex can facilitate adaptation to new environments ( colegrave 2002 ; goddard et  al .
the treatments applied to the experimental hay meadows were initiated in spring 2010 and induced a gradual transition from the original to new plant communities ( andrey et  al .
the mesocosm habitats also represented new conditions for the thrips . yet in both competition assays , the one involving sexual a. elegans and asexual a. rufus collected from the field , and the one involving sexual and asexual a. rufus from laboratory populations , asexuals completely replaced their sexual counterparts .
part of the success of asexuals in the mesocosm experiments could stem from a colonisation advantage at the beginning of the experiment when only 12 individuals were present .
however most mesocosms still contained sexuals when the carrying capacity was reached ( by generation four ) , indicating that a colonization advantage of asexuals is unlikely to fully explain the results .
for asexual a. rufus collected in the field , we further observed an increase of average fitness in two of the three homogeneous mesocosms , and a decrease in fitness variation among females , as expected under adaptation following directional selection .
a possible explanation for this pattern is that the response to selection in sexuals is slow , perhaps because of epistasis effects among loci involved in adaptation to different grasses , while genotype sorting in asexuals may happen more rapidly
. however , the most likely explanation for these unexpected findings is that natural populations of asexual aptinothrips are characterized by more standing genetic variation for adaptation to different host plants and that the sexual species a. elegans is specialized on one or a few grass species .
the same explanation could account for the unexpected success of asexual thrips in experimental plots with high grass diversity .
as asexual thrips would be able to efficiently exploit many different species of grasses , and sexual thrips only few of them , the niche space available to asexual thrips would increase as grass diversities increase .
additional studies are required to infer the fundamental niche of the different sexual and asexual grassthrips lineages and their level of specialism vs generalism with respect to host plant use .
the proposition that increased genetic variation improves a lineage 's ability to exploit alternative niches is at the root of tangled bank related theories ( reviewed in bell 1982 ; song et  al .
however , in contrast to our findings , tangled bank theories assume that sexual species , rather than asexual ones , would feature extensive genetic variation for exploiting alternative niches .
thus , these theories in fact propose that there is an advantage for genetically based diversity in niche usage ; they do not directly show an advantage for sex . in systems where such diversity is larger in asexuals than sexuals ,
as appears to be the case in grassthrips , tangled bank mechanisms predict an advantage for asexual species .
such nichespecific persistence might however be transient , pending colonization of that niche by a new set of asexual clones .
asexual a. rufus populations are distributed throughout the world ( van der kooi and schwander 2014 ; fontcuberta et  al .
2016 ) , and the consequently enormous population sizes may be sufficient for mutation , selection , and migration to generate clones adapted to different grass species . furthermore , asexual a. rufus populations consist of a diverse mixture of clones , including clones that independently derived from sexual ancestors ( fontcuberta et  al .
such high levels of clonal diversity may allow different clones to be selected on alternative plant species . more generally , given that populations of asexual  species  are often a mixture of different and independently derived clones ( e.g. , neiman et  al . 2005 ; janko et  al .
2011 ) , the scope for  genetic diversity  based advantages of sex may be more limited in natural populations than typically assumed ( e.g. , judson 1997 ; ladle et  al .
1993 ; lively 2010 ; neiman and schwander 2011 ) . in summary , we maintained competing sexual and asexual thrips in spatially ( and temporally ) heterogeneous habitat conditions , conditions that have been shown to provide benefits to sex . in spite of this , sexual species decreased in frequency and , in the mesocosm experiments , went extinct .
these findings indicate that habitat heterogeneity is unlikely to maintain sex in natural populations given the direct costs sex generates .
the question , of course , remains ; whether sex in natural grassthrips populations is on the verge of extinction or whether sexual lineages are maintained by parasiteimposed selection or other factors unrelated to variation in grass composition over space or time .
additional supporting information may be found in the online version of this article at the publisher 's website : 
 figure s1 .
correlation between grass diversity and the proportion of sexual thrips separately for each hay meadow . 

details on the treatments applied to the hay meadows used in the field experiment . 

contribution of the different inventoried grass species to the clustering of plant communities click here for additional data file .
 explaining the overwhelming success of sex among eukaryotes is difficult given the obvious costs of sex relative to asexuality . 
 different studies have shown that sex can provide benefits in spatially heterogeneous environments under specific conditions , but whether spatial heterogeneity commonly contributes to the maintenance of sex in natural populations remains unknown . 
 we experimentally manipulated habitat heterogeneity for sexual and asexual thrips lineages in natural populations and under seminatural mesocosm conditions by varying the number of hostplants available to these herbivorous insects . 
 asexual lineages rapidly replaced the sexual ones , independently of the level of habitat heterogeneity in mesocosms . in natural populations , the success of sexual thrips decreased with increasing habitat heterogeneity , with sexual thrips apparently only persisting in certain types of hostplant communities . 
 our results illustrate how genetic diversitybased mechanisms can favor asexuality instead of sex when sexual lineages cooccur with genetically variable asexual lineages . 
Methods
THRIPS FITNESS ON DIFFERENT GRASS SPECIES
MESOCOSM EXPERIMENT
FIELD EXPERIMENT
STATISTICAL ANALYSES
Results
Discussion
Supporting information