The genetics of obligate parthenogenesis in an aphid species and its consequences for the maintenance of alternative reproductive modes

Transcriptomic basis of sex loss in the pea aphid

BACKGROUND

Transitions from sexual to asexual reproduction are common in eukaryotes, but the underlying mechanisms remain poorly known. The pea aphid-Acyrthosiphon pisum-exhibits reproductive polymorphism, with cyclical parthenogenetic and obligate parthenogenetic lineages, offering an opportunity to decipher the genetic basis of sex loss. Previous work on this species identified a single 840 kb region controlling reproductive polymorphism and carrying 32 genes. With the aim of identifying the gene(s) responsible for sex loss and the resulting consequences on the genetic programs controlling sexual or asexual embryogenesis, we compared the transcriptomic response to photoperiod shortening-the main sex-inducing cue-of a sexual and an obligate asexual lineage of the pea aphid, focusing on heads (where the photoperiodic cue is detected) and embryos (the final target of the cue).

RESULTS

Our analyses revealed that four genes (one expressed in the head, and three in the embryos) of the region responded differently to photoperiod in the two lineages. We also found that the downstream genetic programs expressed during embryonic development of a future sexual female encompass ∼1600 genes, among which miRNAs, piRNAs and histone modification pathways are overrepresented. These genes mainly co-localize in two genomic regions enriched in transposable elements (TEs).

CONCLUSIONS

Our results suggest that the causal polymorphism(s) in the 840 kb region somehow impair downstream epigenetic and post-transcriptional regulations in obligate asexual lineages, thereby sustaining asexual reproduction.

Morph-specific fitness throughout the life cycle of the grain aphid, nonhost-alternating, holocyclic Sitobion avenae (Hemiptera: Aphididae)

Aphids exhibit seasonally alternating asexual and sexual reproductive modes. Different morphs are produced throughout the life cycle. To evaluate morph-specific fitness during reproductive switching, holocyclic Sitobion avenae were induced continuously under short light conditions, and development and reproduction were compared in each morph. Seven morphs, including apterous and alate virginoparae, apterous and alate sexuparae, oviparae, males, and fundatrices, were produced during the life cycle. The greatest proportions of sexuparae, oviparae, males, and virginoparae were in the G1, G2, G3, and G4 generations, respectively. Regardless of asexual or sexual morphs, alate morphs exhibited a marked delay in age at maturity compared with that of apterous morphs. Among the alate morphs, males had the longest age at maturity, followed by sexuparae and virginoparae. Among the apterous morphs, sexuparae were older at maturity than the fundatrices, virginoparae, and oviparae. The nymphs of each morph had equal survival potentials. For the same wing morphs, apterous sexuparae and oviparae exhibited substantial delays in the pre-reproductive period and considerable reductions in fecundity, compared with those of apterous virginoparae and fundatrices, whereas alate sexuparae and alate virginoparae had similar fecundity. The seven morphs exhibited Deevey I survivorship throughout the life cycle. These results suggest that sexual production, particularly in males, has short-term development and reproduction costs. The coexistence of sexual and asexual morphs in sexuparae offspring may be regarded as an adaptive strategy for limiting the risk of low fitness in winter.

Contrasting Evolutionary Patterns Between Sexual and Asexual Lineages in a Genomic Region Linked to Reproductive Mode Variation in the pea aphid

Although asexual lineages evolved from sexual lineages in many different taxa, the genetics of sex loss remains poorly understood. We addressed this issue in the pea aphid Acyrthosiphon pisum, whose natural populations encompass lineages performing cyclical parthenogenesis (CP) and producing one sexual generation per year, as well as obligate parthenogenetic (OP) lineages that can no longer produce sexual females but can still produce males. An SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from 6 populations) revealed that an X-linked region is associated with the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2,000 SNPs or short Indels showing a high degree of association with the phenotypic trait. In OP populations specifically, this region also shows reduced diversity and Tajima's D, consistent with the OP phenotype being a derived trait in aphids. Interestingly, the low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.5%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. These genetic exchanges, combined with the selection of the OP and CP reproductive modes under different climates, probably contribute to the long-term persistence of the cp and op alleles.

No recombination suppression in asexually produced males of Daphnia pulex

Obligate parthenogenesis (OP) is often thought to evolve by disruption of reductional meiosis and suppression of crossover recombination. In the crustacean Daphnia pulex, OP lineages, which have evolved from cyclical parthenogenetic (CP) ancestors, occasionally produce males that are capable of reductional meiosis. Here, by constructing high-density linkage maps, we find that these males show only slightly and nonsignificantly reduced recombination rates compared to CP males and females. Both meiosis disruption and recombination suppression are therefore sex-limited (or partly so), which speaks against the evolution of OP by disruption of a gene that is essential for meiosis or recombination in both sexes. The findings may be explained by female-limited action of genes that suppress recombination, but previously identified candidate genes are known to be expressed in both sexes. Alternatively, and equally consistent with the data, OP might have evolved through a reuse of the parthenogenesis pathways already present in CP and through their extension to all events of oogenesis. The causal mutations for the CP to OP transition may therefore include mutations in genes involved in oogenesis regulation and may not necessarily be restricted to genes of the "meiosis toolkit." More generally, our study emphasizes that there are many ways to achieve asexuality, and elucidating the possible mechanisms is key to ultimately identify the genes and traits involved.

Evolutionary consequences of loss of sexual reproduction on male-related traits in parthenogenetic lineages of the pea aphid

Transition from sexual reproduction to parthenogenesis constitutes a major life-history change with deep evolutionary consequences for sex-related traits, which are expected to decay. The pea aphid Acyrthosiphon pisum shows intraspecific reproductive polymorphism, with cold-resistant cyclically parthenogenetic (CP) lineages that alternate sexual and asexual generations and cold-sensitive obligately parthenogenetic (OP) lineages that produce only asexual females but still males. Here, the genotyping of 219 pea aphid lineages collected in cold-winter and mild-winter regions revealed contrasting population structures. Samples from cold-winter regions consisted mostly of distinct multilocus genotypes (MLGs) usually represented by a single sample (101 different MLGs for 111 samples) and were all phenotyped as CP. In contrast, fewer MLGs were found in mild-winter regions (28 MLGs for 108 samples), all but one being OP. Since the males produced by OP lineages are unlikely to pass on their genes (sexual females being rare in mild-winter regions), we tested the hypothesis that their traits could degenerate due to lack of selection by comparing male production and male reproductive success between OP and CP lineages. Male production was indeed reduced in OP lineages, but a less clear pattern was observed for male reproductive success: females mated with OP males laid fewer eggs (fertilized or not) but OP and CP males fertilized the same proportion of eggs. These differences may stem from the type of selective forces: male production may be counter-selected whereas male performances may evolve under the slower process of relaxed selection. The overall effective reproductive capacity of OP males could result from recent sex loss in OP lineages or underestimated reproductive opportunities.

Comparative Proteomic Analysis Provides New Insights into the Molecular Basis of Thermal-Induced Parthenogenesis in Silkworm (Bombyx mori)

Artificial parthenogenetic induction via thermal stimuli in silkworm is an important technique that has been used in sericultural production. However, the molecular mechanism underlying it remains largely unknown. We have created a fully parthenogenetic line (PL) with more than 85% occurrence and 80% hatching rate via hot water treatment and genetic selection, while the parent amphigenetic line (AL) has less than 30% pigmentation rate and less than 1% hatching rate when undergoing the same treatment. Here, isobaric tags for relative and absolute quantitation (iTRAQ)-based analysis were used to investigate the key proteins and pathways associated with silkworm parthenogenesis. We uncovered the unique proteomic features of unfertilized eggs in PL. In total, 274 increased abundance proteins and 211 decreased abundance proteins were identified relative to AL before thermal induction. Function analysis displayed an increased level of translation and metabolism in PL. After thermal induction, 97 increased abundance proteins and 187 decreased abundance proteins were identified. An increase in stress response-related proteins and decrease in energy metabolism suggested that PL has a more effective response to buffer the thermal stress than AL. Cell cycle-related proteins, including histones, and spindle-related proteins were decreased in PL, indicating an important role of this decrease in the process of ameiotic parthenogenesis.

A comprehensive review on the documented characteristics of four Reticulitermes termites (Rhinotermitidae, Blattodea) of China

Termites are known as social insects worldwide. Presently in China 473 species, 44 genera and 4 families of termites have been reported. Of them, 111 Reticulitermes species are widely spread in different zones of China. The dispersion flight season of these Chinese Reticulitermes species are usually started from February to June, but in some regions different species are distributed, sharing their boundaries and having overlapping flight seasons. These reasons become important sources of hybridization between two different heterospecific populations of termites. It was confirmed that the fertilized eggs and unfertilized eggs of some Reticulitermes termites have the capacity of cleavage. While the unfertilized eggs of R. aculabialis, R. chinensis and R. labralis cleaved normally and the only R. aculabialis unfertilized eggs develop in embryos. While, the R. flaviceps and R. chinensis were observed with their abnormal embryonic development, and not hatching of eggs parthenogenetically. They were reported more threatening to Chinese resources as they propagate with parthenogenesis, hybridization and sexual reproduction. Eggshell and macrophiles of eggs play important roles in species identification and control. Although, they are severe pests and cause a wide range of damages to wooden structures and products in homes, buildings, building materials, trees, crops, and forests in China's Mainland.

Transitions to asexuality and evolution of gene expression in Artemia brine shrimp

While sexual reproduction is widespread among many taxa, asexual lineages have repeatedly evolved from sexual ancestors. Despite extensive research on the evolution of sex, it is still unclear whether this switch represents a major transition requiring major molecular reorganization, and how convergent the changes involved are. In this study, we investigated the phylogenetic relationship and patterns of gene expression of sexual and asexual lineages of Eurasian Artemia brine shrimp, to assess how gene expression patterns are affected by the transition to asexuality. We find only a few genes that are consistently associated with the evolution of asexuality, suggesting that this shift may not require an extensive overhauling of the meiotic machinery. While genes with sex-biased expression have high rates of expression divergence within Eurasian Artemia, neither female- nor male-biased genes appear to show unusual evolutionary patterns after sexuality is lost, contrary to theoretical expectations.

Global patterns in genomic diversity underpinning the evolution of insecticide resistance in the aphid crop pest Myzus persicae

The aphid Myzus persicae is a destructive agricultural pest that displays an exceptional ability to develop resistance to both natural and synthetic insecticides. To investigate the evolution of resistance in this species we generated a chromosome-scale genome assembly and living panel of >110 fully sequenced globally sampled clonal lines. Our analyses reveal a remarkable diversity of resistance mutations segregating in global populations of M. persicae. We show that the emergence and spread of these mechanisms is influenced by host-plant associations, uncovering the widespread co-option of a host-plant adaptation that also offers resistance against synthetic insecticides. We identify both the repeated evolution of independent resistance mutations at the same locus, and multiple instances of the evolution of novel resistance mechanisms against key insecticides. Our findings provide fundamental insights into the genomic responses of global insect populations to strong selective forces, and hold practical relevance for the control of pests and parasites.

Insect alarm pheromones in response to predators: Ecological trade-offs and molecular mechanisms

Insect alarm pheromones are chemical substances that are synthesized and released in response to predators to reduce predation risk. Alarm pheromones can also be perceived by predators, who take advantage of alarm cues to locate prey. While selection favors evolution of alarm pheromone signals that are not easily detectable by predators, predator evolution selects for better prey detection ability. Here, we review the diversity of alarm signals, and consider the behavioral and ecological conditions under which they have evolved. We show that components of alarm pheromones are similar across many insects, although aphids exhibit different behavioral responses to alarm cues compared to social insects. The effects of alarm pheromones on prey behavior depend on factors such as the concentration of pheromones and the density of conspecifics. We also discuss the molecular mechanisms of alarm pheromone perception underlying the evolutionary arms race between predators and prey, and the function of olfactory proteins and receptors in particular. Our review provides a novel synthesis of the diversity and function of insect alarm pheromones, while suggesting avenues that might better allow researchers to exploit population-level responses to alarm signaling for the sustainable management of pests and vector-borne pathogens.

A functional genomics approach to dissect spotted alfalfa aphid resistance in Medicago truncatula

Aphids are virus-spreading insect pests affecting crops worldwide and their fast population build-up and insecticide resistance make them problematic to control. Here, we aim to understand the molecular basis of spotted alfalfa aphid (SAA) or Therioaphis trifolii f. maculata resistance in Medicago truncatula, a model organism for legume species. We compared susceptible and resistant near isogenic Medicago lines upon SAA feeding via transcriptome sequencing. Expression of genes involved in defense and stress responses, protein kinase activity and DNA binding were enriched in the resistant line. Potentially underlying some of these changes in gene expression was the finding that members of the MYB, NAC, AP2 domain and ERF transcription factor gene families were differentially expressed in the resistant versus susceptible lines. A TILLING population created in the resistant cultivar was screened using exome capture sequencing and served as a reverse genetics tool to functionally characterise genes involved in the aphid resistance response. This screening revealed three transcription factors (a NAC, AP2 domain and ERF) as important regulators in the defence response, as a premature stop-codon in the resistant background led to a delay in aphid mortality and enhanced plant susceptibility. This combined functional genomics approach will facilitate the future development of pest resistant crops by uncovering candidate target genes that can convey enhanced aphid resistance.

Gene Expression Profiling Indicated Diverse Functions and Characteristics of Core Genes in Pea Aphid

The pea aphid is a global insect pest, and variable phenotypes can be produced by pea aphids in the same genotype in response to changes in external environmental factors. However, detailed dynamic gene regulation networks and the core markers involved in different biological processes of pea aphids have not yet been reported. In this study, we obtained the published genomic and transcriptomic data, and performed transcriptome profiling of five pea aphid morphs (winged asexual female, wingless asexual female, wingless sexual female, winged male and wingless male) from each of three pea aphid genotypes, i.e., the transcriptomes from a total of 15 types of pea aphids were analyzed and the type-specific expression of genes in five different morphs was identified. The expression profiling was verified by quantitative real-time PCR (qPCR) analysis. Moreover, we determined the expression features and co-expression networks of highly variable genes. We also used the ARACNe method to obtain 263 core genes related to different biological pathways. Additionally, eight of the identified genes were aligned with transcription factor families, indicating that they act as transcription factors and regulate downstream genes. Furthermore, we found reliable markers using random forest methodology to distinguish different morphs of pea aphids. Our study provides a systematic and comprehensive approach for analyzing the core genes that may play important roles in a multitude of biological processes from the insect transcriptomes.

Genetic Divergence of Two Sitobion avenae Biotypes on Barley and Wheat in China

Host plant affinity and geographic distance can play critical roles in the genetic divergence of insect herbivores and evolution of insect biotypes, but their relative importance in the divergence of insect populations is still poorly understood. We used microsatellite markers to test the effects of host plant species and geographic distance on divergence of two biotypes of the English grain aphid, Sitobion avenae (Fabricius). We found that clones of S. avenae from western provinces (i.e., Xinjiang, Gansu, Qinghai and Shaanxi) had significantly higher genetic diversity than those from eastern provinces (i.e., Anhui, Henan, Hubei, Zhejiang and Jiangsu), suggesting their differentiation between both areas. Based on genetic diversity and distance estimates, biotype 1 clones of eastern provinces showed high genetic divergence from those of western provinces in many cases. Western clones of S. avenae also showed higher genetic divergence among themselves than eastern clones. The Mantel test identified a significant isolation-by-distance (IBD) effect among different geographic populations of S. avenae, providing additional evidence for a critical role of geography in the genetic structure of both S. avenae biotypes. Genetic differentiation (i.e., F_ST) between the two biotypes was low in all provinces except Shaanxi. Surprisingly, in our analyses of molecular variance, non-significant genetic differentiation between both biotypes or between barley and wheat clones of S. avenae was identified, showing little contribution of host-plant associated differentiation to the divergence of both biotypes in this aphid. Thus, it is highly likely that the divergence of the two S. avenae biotypes involved more geographic isolation and selection of some form than host plant affinity. Our study can provide insights into understanding of genetic structure of insect populations and the divergence of insect biotypes.

Parthenogenesis and developmental constraints

The absence of a paternal contribution in an unfertilized ovum presents two developmental constraints against the evolution of parthenogenesis. We discuss the constraint caused by the absence of a centrosome and the one caused by the missing set of chromosomes and how they have been broken in specific taxa. They are examples of only a few well-underpinned examples of developmental constraints acting at macro-evolutionary scales in animals. Breaking of the constraint of the missing chromosomes is the best understood and generally involves rare occasions of drastic changes of meiosis. These drastic changes can be best explained by having been induced, or at least facilitated, by sudden cytological events (e.g., repeated rounds of hybridization, endosymbiont infections, and contagious infections). Once the genetic and developmental machinery is in place for regular or obligate parthenogenesis, shifts to other types of parthenogenesis can apparently rather easily evolve, for example, from facultative to obligate parthenogenesis, or from pseudoarrhenotoky to haplodiploidy. We argue that the combination of the two developmental constraints forms a near-absolute barrier against the gradual evolution from sporadic to obligate or regular facultative parthenogenesis, which can probably explain why the occurrence of the highly advantageous mode of regular facultative parthenogenesis is so rare and entirely absent in vertebrates.

Differences in egg hatching time between cyclical and obligate parthenogenetic lineages of aphids

Many aphid species exhibit a variation in reproductive mode which is influenced by winter climate regimes, with cyclical parthenogenetic (CP) lines dominating in cold winter areas (because they produce cold-resistant eggs) and obligate parthenogenetic (OP) ones in mild winter regions (because of their parthenogenetic overwintering). Genetic studies on several aphid species have shown that the OP trait can be transmitted during sexual events involving the 2 types of lines. This genetic system could be considered as a local safeguarding mechanism for OP alleles in case severe frost would have killed all parthenogenetically overwintering individuals. However, this strategy would only be efficient in restoring local polymorphism in breeding systems if the newly hatched OP recombinants remain competitive over their CP counterparts. In this study we compared egg hatching sequences of CP and OP F1 clones from several crosses obtained for 2 cereal aphid species, Sitobion avenae (constant 5 °C, 8 h of light) and Rhopalosiphum padi (winter outdoor conditions). For S. avenae, we obtained F1 offspring from 6 crosses, involving 4 clones while in R. padi F1 were obtained from 11 crosses involving 14 clones. We showed that in both species proportions of OP clones were higher in the first half of the progeny relative to the second half. In addition, F1 OP clones hatched in the mean about a week earlier than their CP sibs, which gives them a demographic advantage at the start of the growth season. We then discussed the consequences of this fitness advantage for the maintenance and spread of the OP trait in aphid populations.

Biological and genetic features of introduced aphid populations in agroecosystems

In agroecosystems, introduced aphids that reproduce by obligate parthenogenesis (OP) show strong biased representation of a few genotypes (superclones), whereas species with cyclical parthenogenesis (CP) exhibit the opposite trend with many unique genotypes. We analyzed the biological and genetic features of 23 different aphid species introduced in different geographic areas and climates, finding putative superclones in about 60% of them. We have examined the proximal causes for aphid establishment and spread after their introduction, and found that OP, host availability, and phenotypic plasticity are among the main variables underpinning the ability of aphids to succeed in new geographic areas, which may explain the high potential for invasion in this group of pest insects.

Diapause and maintenance of facultative sexual reproductive strategies

Facultative sex combines sexual and asexual reproduction in the same individual (or clone) and allows for a large diversity of life-history patterns regarding the timing, frequency and intensity of sexual episodes. In addition, other life-history traits such as a diapause stage may become linked to sex. Here, we develop a matrix modelling framework for addressing the cost of sex in facultative sexuals, in constant, periodic and stochastically fluctuating environments. The model is parametrized using life-history data from Brachionus calyciflorus, a facultative sexual rotifer in which sex and diapause are linked. Sexual propensity was an important driver of costs in constant environments, in which high costs (always > onefold, and sometimes > twofold) indicated that asexuals should outcompete facultative sexuals. By contrast, stochastic environments with high temporal autocorrelation favoured facultative sex over obligate asex, in particular, if the penalty to fecundity in 'bad' environments was large. In such environments, obligate asexuals were constrained by their life cycle length (i.e. time from birth to last reproductive adult age class), which determined an upper limit to the number of consecutive bad periods they could tolerate. Nevertheless, when facultative asexuals with different sexual propensities competed simultaneously against each other and asex, the lowest sex propensity was the most successful in stochastic environments with positive autocorrelation. Our results suggest that a highly specific mechanism (i.e. diapause linked to sex) can alone stabilize facultative sex in these animals, and protect it from invasion of both asexual and pure sexual strategies.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Reproductive switching analysis of Daphnia similoides between sexual female and parthenogenetic female by transcriptome comparison

The water flea Daphnia are planktonic crustaceans commonly found in freshwater environment that can switch their reproduction mode from parthenogenesis to sexual reproduction to adapt to the external environment. As such, Daphnia are great model organisms to study the mechanism of reproductive switching, the underlying mechanism of reproduction and development in cladocerans and other animals. However, little is known about the Daphnia's reproductive behaviour at a molecular level. We constructed a genetic database of the genes expressed in a sexual female (SF) and a parthenogenetic female (PF) of D. similoides using Illumina HiSeq 2500. A total of 1,763 differentially expressed genes (865 up- and 898 down-regulated) were detected in SF. Of the top 30 up-regulated SF unigenes, the top 4 unigenes belonged to the Chitin_bind_4 family. In contrast, of the top down-regulated SF unigenes, the top 3 unigenes belonged to the Vitellogenin_N family. This is the first study to indicate genes that may have a crucial role in reproductive switching of D. similoides, which could be used as candidate genes for further functional studies. Thus, this study provides a rich resource for investigation and elucidation of reproductive switching in D. similoides.

Detection of Gene Flow from Sexual to Asexual Lineages in Thrips tabaci (Thysanoptera: Thripidae)

Populations of Thrips tabaci are known to have two sympatric but genetically isolated reproductive modes, arrhenotoky (sexual reproduction) and thelytoky (asexual reproduction). Herein, we report behavioral, ecological and genetic studies to determine whether there is gene flow between arrhenotokous and thelytokous T. tabaci. We did not detect significant preference by arrhenotokous males to mate with females of a particular reproductive mode, nor did we detect significant behavioral differences between arrhenotokous males mated with arrhenotokous or thelytokous females in their pre-copulation, copulation duration and mating frequency. Productive gene transfer resulting from the mating between the two modes was experimentally confirmed. Gene transfer from arrhenotokous T. tabaci to thelytokous T. tabaci was further validated by confirmation of the passage of the arrhenotokous male-originated nuclear gene (histone H3 gene) allele to the F₂ generation. These behavioral, ecological and genetic studies confirmed gene transfer from the sexual arrhenotokous mode to the asexual thelytokous mode of T. tabaci in the laboratory. These results demonstrate that asexual T. tabaci populations may acquire genetic variability from sexual populations, which could offset the long-term disadvantage of asexual reproduction.

Genetic control of contagious asexuality in the pea aphid

Although evolutionary transitions from sexual to asexual reproduction are frequent in eukaryotes, the genetic bases of such shifts toward asexuality remain largely unknown. We addressed this issue in an aphid species where both sexual and obligate asexual lineages coexist in natural populations. These sexual and asexual lineages may occasionally interbreed because some asexual lineages maintain a residual production of males potentially able to mate with the females produced by sexual lineages. Hence, this species is an ideal model to study the genetic basis of the loss of sexual reproduction with quantitative genetic and population genomic approaches. Our analysis of the co-segregation of ∼ 300 molecular markers and reproductive phenotype in experimental crosses pinpointed an X-linked region controlling obligate asexuality, this state of character being recessive. A population genetic analysis (>400-marker genome scan) on wild sexual and asexual genotypes from geographically distant populations under divergent selection for reproductive strategies detected a strong signature of divergent selection in the genomic region identified by the experimental crosses. These population genetic data confirm the implication of the candidate region in the control of reproductive mode in wild populations originating from 700 km apart. Patterns of genetic differentiation along chromosomes suggest bidirectional gene flow between populations with distinct reproductive modes, supporting contagious asexuality as a prevailing route to permanent parthenogenesis in pea aphids. This genetic system provides new insights into the mechanisms of coexistence of sexual and asexual aphid lineages.

Evolution of asexuality via different mechanisms in grass thrips (thysanoptera: Aptinothrips)

Asexual lineages can derive from sexual ancestors via different mechanisms and at variable rates, which affects the diversity of the asexual population and thereby its ecological success. We investigated the variation and evolution of reproductive systems in Aptinothrips, a genus of grass thrips comprising four species. Extensive population surveys and breeding experiments indicated sexual reproduction in A. elegans, asexuality in A. stylifer and A. karnyi, and both sexual and asexual lineages in A. rufus. Asexuality in A. stylifer and A. rufus coincides with a worldwide distribution, with sexual A. rufus lineages confined to a limited area. Inference of molecular phylogenies and antibiotic treatment revealed different causes of asexuality in different species. Asexuality in A. stylifer and A. karnyi has most likely genetic causes, while it is induced by endosymbionts in A. rufus. Endosymbiont-community characterization revealed presence of Wolbachia, and lack of other bacteria known to manipulate host reproduction. However, only 69% asexual A. rufus females are Wolbachia-infected, indicating that either an undescribed endosymbiont causes asexuality in this species or that Wolbachia was lost in several lineages that remained asexual. These results open new perspectives for studies on the maintenance of mixed sexual and asexual reproduction in natural populations.

Coexistence of sexual individuals and genetically isolated asexual counterparts in a thrips

Sex is a paradoxical phenomenon because it is less efficient compared with asexual reproduction. To resolve this paradox we need a direct comparison between sexual and asexual forms. In many organisms, however, sexual and asexual forms do not occur in the same habitat, or at the same time. In a few cases where sexual and asexual forms are found in a single population, some (though rare) genetic exchange is usually detected between the two forms. When genetic exchange occurs a direct comparison is impossible. Here we investigate a thrips exhibiting both sexual and asexual forms (lineages) that are morphologically indistinguishable. We examine if the two forms are genetically isolated. Phylogeny based on nuclear genes confirms that the sexual and asexual lineages are genetically differentiated. Thus we demonstrate that the current system has certain advantages over existing and previously used model systems in the evolution of sexual reproduction.

Williams' paradox and the role of phenotypic plasticity in sexual systems

As George Williams pointed out in 1975, although evolutionary explanations, based on selection acting on individuals, have been developed for the advantages of simultaneous hermaphroditism, sequential hermaphroditism and gonochorism, none of these evolutionary explanations adequately explains the current distribution of these sexual systems within the Metazoa (Williams' Paradox). As Williams further pointed out, the current distribution of sexual systems is explained largely by phylogeny. Since 1975, we have made a great deal of empirical and theoretical progress in understanding sexual systems. However, we still lack a theory that explains the current distribution of sexual systems in animals and we do not understand the evolutionary transitions between hermaphroditism and gonochorism. Empirical data, collected over the past 40 years, demonstrate that gender may have more phenotypic plasticity than was previously realized. We know that not only sequential hermaphrodites, but also simultaneous hermaphrodites have phenotypic plasticity that alters sex allocation in response to social and environmental conditions. A focus on phenotypic plasticity suggests that one sees a continuum in animals between genetically determined gonochorism on the one hand and simultaneous hermaphroditism on the other, with various types of sequential hermaphroditism and environmental sex determination as points along the spectrum. Here I suggest that perhaps the reason we have been unable to resolve Williams' Paradox is because the problem was not correctly framed. First, because, for example, simultaneous hermaphroditism provides reproductive assurance or dioecy ensures outcrossing does not mean that there are no other evolutionary paths that can provide adaptive responses to those selective pressures. Second, perhaps the question we need to ask is: What selective forces favor increased versus reduced phenotypic plasticity in gender expression? It is time to begin to look at the question of sexual system as one of understanding the timing and degree of phenotypic plasticity in gender expression in the life history in terms of selection acting on a continuum, rather than on a set of discrete sexual systems.