The evolution of geographic parthenogenesis in Timema walking-sticks

Does ecology shape geographical parthenogenesis? Evidence from the facultatively parthenogenetic stick insect Megacrania batesii

Closely related sexual and parthenogenetic species often show distinct distribution patterns, known as geographical parthenogenesis. Similar patterns, characterized by the existence of separate sexual and parthenogenetic populations across their natural range, can also be found in facultative parthenogens - species in which every female is capable of both sexual and parthenogenetic reproduction. The underlying mechanisms driving this phenomenon in nature remain unclear. Features of the habitat, such as differences in host-plant phenotypes or niche breadth, could favour sexual or asexual reproductive modes and thus help to explain geographical parthenogenesis in natural insect populations. Megacrania batesii is a facultatively parthenogenetic stick insect that displays geographical parthenogenesis in the wild. We aimed to explore whether sexual and parthenogenetic populations of M. batesii displayed niche differentiation or variations in niche breadth that could explain the separation of the two population types. To do this, we sampled host plants from across the range of M. batesii and quantified phenotypic traits that might affect palatability or accessibility for M. batesii, including leaf thickness, toughness, spike size and density, plant height, and chemical composition. We also quantified host-plant density, which could affect M. batesii dispersal. We found little evidence of phenotypic differences between host plants supporting sexual versus asexual M. batesii populations, and no difference in host-plant density or niche breadth between the two population types. Our results suggest that habitat parameters do not play a substantial role in shaping patterns of geographical parthenogenesis in wild populations of M. batesii. Instead, population sex ratio variation could result from interactions between the sexes or dispersal dynamics.

Genetic and Phenotypic Consequences of Local Transitions between Sexual and Parthenogenetic Reproduction in the Wild

AbstractTransitions from sexual to asexual reproduction have occurred in numerous lineages, but it remains unclear why asexual populations rarely persist. In facultatively parthenogenetic animals, all-female populations can arise when males are absent or become extinct, and such populations could help to understand the genetic and phenotypic changes that occur in the initial stages of transitions to asexuality. We investigated a naturally occurring spatial mosaic of mixed-sex and all-female populations of the facultatively parthenogenetic Australian phasmid Megacrania batesii. Analysis of single-nucleotide polymorphisms indicated multiple independent transitions between reproductive modes. All-female populations had much lower heterozygosity and allelic diversity than mixed-sex populations, but we found few consistent differences in fitness-related traits between population types. All-female populations exhibited more frequent and severe deformities in their (flight-incapable) wings but did not show higher rates of appendage loss. All-female populations also harbored more ectoparasites in swamp (but not beach) habitats. Reproductive mode explained little variation in female body size, fecundity, or egg hatch rate. Our results suggest that transitions to parthenogenetic reproduction can lead to dramatic genetic changes with little immediate effect on performance. All-female M. batesii populations appear to consist of high-fitness genotypes that might be able to thrive for many generations in relatively constant and benign environments but could be vulnerable to environmental challenges, such as increased parasite abundance.

Evidence for cryptic sex in parthenogenetic stick insects of the genus Timema

Obligately parthenogenetic species are expected to be short lived since the lack of sex and recombination should translate into a slower adaptation rate and increased accumulation of deleterious alleles. Some, however, are thought to have been reproducing without males for millions of years. It is not clear how these old parthenogens can escape the predicted long-term costs of parthenogenesis, but an obvious explanation is cryptic sex. In this study, we screen for signatures of cryptic sex in eight populations of four parthenogenetic species of Timema stick insects, some estimated to be older than 1 Myr. Low genotype diversity, homozygosity of individuals and high linkage disequilibrium (LD) unaffected by marker distances support exclusively parthenogenetic reproduction in six populations. However, in two populations (namely, of the species Timema douglasi and T. monikensis) we find strong evidence for cryptic sex, most likely mediated by rare males. These populations had comparatively high genotype diversities, lower LD, and a clear LD decay with genetic distance. Rare sex in species that are otherwise largely parthenogenetic could help explain the unusual success of parthenogenesis in the Timema genus and raises the question whether episodes of rare sex are in fact the simplest explanation for the persistence of many old parthenogens in nature.

Descriptions of two new stick insect species of Cnipsomorpha Hennemann, Conle, Zhang & Liu (Phasmatodea) from China based on integrative taxonomy

Accurate taxonomical identification is an extremely important basis for stick insect research, including evolutionary biology but also applied biology such as pest control. In addition, genetic methods are a valuable identification auxiliary technology at present. Therefore, this paper used morphological and molecular data to investigate five stick insect specimens from the genus Cnipsomorpha in Yunnan, successfully identifying two new species: Cnipsomorphayunnanensis Xu, Jiang & Yang, sp. nov. and C.yuxiensis Xu, Jiang & Yang, sp. nov. A phylogenetic tree was constructed through their 28S and COI genes in order to infer the phylogenetic position of the two new species. Photographs of the new species and a key to all known Cnipsomorpha species are provided.

Climatic similarity and genomic background shape the extent of parallel adaptation in Timema stick insects

Evolution can repeat itself, resulting in parallel adaptations in independent lineages occupying similar environments. Moreover, parallel evolution sometimes, but not always, uses the same genes. Two main hypotheses have been put forth to explain the probability and extent of parallel evolution. First, parallel evolution is more likely when shared ecologies result in similar patterns of natural selection in different taxa. Second, parallelism is more likely when genomes are similar because of shared standing variation and similar mutational effects in closely related genomes. Here we combine ecological, genomic, experimental and phenotypic data with Bayesian modelling and randomization tests to quantify the degree of parallelism and its relationship with ecology and genetics. Our results show that the extent to which genomic regions associated with climate are parallel among species of Timema stick insects is shaped collectively by shared ecology and genomic background. Specifically, the extent of genomic parallelism decays with divergence in climatic conditions (that is, habitat or ecological similarity) and genomic similarity. Moreover, we find that climate-associated loci are likely subject to selection in a field experiment, overlap with genetic regions associated with cuticular hydrocarbon traits and are not strongly shaped by introgression between species. Our findings shed light on when evolution is most expected to repeat itself.

Phylogeography of an endemic California silkmoth genus suggests the importance of an unheralded central California province in generating regional endemic biodiversity

The California Floristic province is a biodiversity hotspot. Understanding the phylogeographic patterns that exist in this diverse region is essential to understanding its evolution and for guiding conservation efforts. Calosaturnia is a charismatic silkmoth genus endemic to large portions of the region with three described species, C. mendocino, C. walterorum, and C. albofasciata. We sampled all three species from across their ranges, sequenced 1463 bp of mitochondrial COI and 1941 bp of nuclear DNA from three genes, and reconstructed phylogenetic relationships and estimated divergence times within the lineages. All three species show pronounced evidence of isolation and, in two cases, secondary reconnection. An unexpected monophyletic mtDNA lineage was found in the Central Coast region, in a region thought to represent an intergrade between C. mendocino and C. walterorum. Our genetic data also significantly revise previous hypotheses as to which species occur in which regions, suggesting that historical ecological changes around four Ma ago isolated some lineages, and a secondary isolation event two Ma ago led to isolation of populations both in the Central Coast region and the southern Sierra Nevada. Our results add to a currently under-appreciated pattern suggesting that coastal Central California is not a transition zone between Northern and Southern California Floristic Province faunas but rather its own unique, periodically isolated, biogeographic region. They also suggest cryptic diversity may be present in many other currently unrecognized groups. Additional research should focus on this central California region because many species may be highly restricted in range and in need of conservation attention.

Integrated phylogenomic and fossil evidence of stick and leaf insects (Phasmatodea) reveal a Permian-Triassic co-origination with insectivores

Stick and leaf insects (Phasmatodea) are a distinctive insect order whose members are characterized by mimicking various plant tissues such as twigs, foliage and bark. Unfortunately, the phylogenetic relationships among phasmatodean subfamilies and the timescale of their evolution remain uncertain. Recent molecular clock analyses have suggested a Cretaceous-Palaeogene origin of crown Phasmatodea and a subsequent Cenozoic radiation, contrasting with fossil evidence. Here, we analysed transcriptomic data from a broad diversity of phasmatodeans and, combined with the assembly of a new suite of fossil calibrations, we elucidate the evolutionary history of stick and leaf insects. Our results differ from recent studies in the position of the leaf insects (Phylliinae), which are recovered as sister to a clade comprising Clitumninae, Lancerocercata, Lonchodinae, Necrosciinae and Xenophasmina. We recover a Permian to Triassic origin of crown Phasmatodea coinciding with the radiation of early insectivorous parareptiles, amphibians and synapsids. Aschiphasmatinae and Neophasmatodea diverged in the Jurassic-Early Cretaceous. A second spur in origination occurred in the Late Cretaceous, coinciding with the Cretaceous Terrestrial Revolution, and was probably driven by visual predators such as stem birds (Enantiornithes) and the radiation of angiosperms.

Large-scale mutation in the evolution of a gene complex for cryptic coloration

The types of mutations affecting adaptation in the wild are only beginning to be understood. In particular, whether structural changes shape adaptation by suppressing recombination or by creating new mutations is unresolved. Here, we show that multiple linked but recombining loci underlie cryptic color morphs of Timema chumash stick insects. In a related species, these loci are found in a region of suppressed recombination, forming a supergene. However, in seven species of Timema, we found that a megabase-size "supermutation" has deleted color loci in green morphs. Moreover, we found that balancing selection likely contributes more to maintaining this mutation than does introgression. Our results show how suppressed recombination and large-scale mutation can help to package gene complexes into discrete units of diversity such as morphs, ecotypes, or species.

Adaptive zones shape the magnitude of premating reproductive isolation in Timema stick insects

Simpson's fossil-record inspired model of 'adaptive zones' proposes that evolution is dominated by small fluctuations within adaptive zones, occasionally punctuated by larger shifts between zones. This model can help explain why the process of population divergence often results in weak or moderate reproductive isolation (RI), rather than strong RI and distinct species. Applied to the speciation process, the adaptive zones hypothesis makes two inter-related predictions: (i) large shifts between zones are relatively rare, (ii) when large shifts do occur they generate stronger RI than shifts within zones. Here, we use ecological, phylogenetic and behavioural data to test these predictions in Timema stick insects. We show that host use in Timema is dominated by moderate shifts within the systematic divisions of flowering plants and conifers, with only a few extreme shifts between these divisions. However, when extreme shifts occur, they generate greater RI than do more moderate shifts. Our results support the adaptive zones model, and suggest that the net contribution of ecological shifts to diversification is dependent on both their magnitude and frequency. We discuss the generality of our findings in the light of emerging evidence from diverse taxa that the evolution of RI is not always the only factor determining the origin of species diversity. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Loss and gain of sexual reproduction in the same stick insect

The outcome of competition between different reproductive strategies within a single species can be used to infer selective advantage of the winning strategy. Where multiple populations have independently lost or gained sexual reproduction it is possible to investigate whether the advantage is contingent on local conditions. In the New Zealand stick insect Clitarchus hookeri, three populations are distinguished by recent change in reproductive strategy and we determine their likely origins. One parthenogenetic population has established in the United Kingdom and we provide evidence that sexual reproduction has been lost in this population. We identify the sexual population from which the parthenogenetic population was derived, but show that the UK females have a post‐mating barrier to fertilisation. We also demonstrate that two sexual populations have recently arisen in New Zealand within the natural range of the mtDNA lineage that otherwise characterizes parthenogenesis in this species. We infer independent origins of males at these two locations using microsatellite genotypes. In one population, a mixture of local and nonlocal alleles suggested males were the result of invasion. Males in another population were most probably the result of loss of an X chromosome that produced a male phenotype in situ. Two successful switches in reproductive strategy suggest local competitive advantage for outcrossing over parthenogenetic reproduction. Clitarchus hookeri provides remarkable evidence of repeated and rapid changes in reproductive strategy, with competitive outcomes dependent on local conditions.

Parthenogenetic Females of the Stick Insect Clitarchus hookeri Maintain Sexual Traits

The New Zealand stick insect Clitarchus hookeri has both sexual and parthenogenetic (all-female) populations. Sexual populations exhibit a scramble competition mating system with distinctive sex roles, where females are signalers and males are searchers, which may lead to differences in the chemical and morphological traits between sexes. Evidence from a range of insect species has shown a decay of sexual traits is common in parthenogenetic lineages, especially those traits related to mate attraction and location, presumably due to their high cost. However, in some cases, sexual traits remain functional, either due to the recent evolution of the parthenogenetic lineage, low cost of maintenance, or because there might be an advantage in maintaining them. We measured morphological and chemical traits of C. hookeri to identify differences between males and females and between females from sexual and parthenogenetic populations. We also tested the ability of males to discriminate between sexual and parthenogenetic females in a laboratory bioassay. Our results show that male C. hookeri has morphological traits that facilitate mobility (smaller body with disproportionately longer legs) and mate detection (disproportionately longer antennae), and adult females release significantly higher amounts of volatile organic compounds than males when this species is sexually active, in accordance with their distinctive sex roles. Although some differences were detected between sexual and parthenogenetic females, the latter appear to maintain copulatory behaviors and chemical signaling. Males were unable to distinguish between sexual and parthenogenetic females, suggesting that there has been little decay in the sexual traits in the parthenogenetic lineage of C. hookeri.

Consequences of Asexuality in Natural Populations: Insights from Stick Insects

Recombination is a fundamental process with significant impacts on genome evolution. Predicted consequences of the loss of recombination include a reduced effectiveness of selection, changes in the amount of neutral polymorphisms segregating in populations, and an arrest of GC-biased gene conversion. Although these consequences are empirically well documented for nonrecombining genome portions, it remains largely unknown if they extend to the whole genome scale in asexual organisms. We identify the consequences of asexuality using de novo transcriptomes of five independently derived, obligately asexual lineages of stick insects, and their sexual sister-species. We find strong evidence for higher rates of deleterious mutation accumulation, lower levels of segregating polymorphisms and arrested GC-biased gene conversion in asexuals as compared with sexuals. Taken together, our study conclusively shows that predicted consequences of genome evolution under asexuality can indeed be found in natural populations.

The geography of sex: sexual conflict, environmental gradients and local loss of sex in facultatively parthenogenetic animals

Obligately asexual organisms tend to occur at higher altitudes or latitudes and occupy larger ranges than their obligately sexual relatives-a phenomenon called geographical parthenogenesis. Some facultatively parthenogenetic organisms that reproduce both sexually and asexually also exhibit spatial variation in reproductive mode. Theory suggests that sexual conflict and mate limitation can determine the relative frequency of sex in facultative parthenogens, but the effect of these dynamics on spatial distributions is unknown. Here, we use individual-based models to investigate whether these dynamics can generate local differences in the reproductive mode in a facultatively parthenogenetic metapopulation occupying an environmental gradient. We find that selection for resistance and high fecundity creates positive epistasis in virgin females between a mutant allele for parthenogenesis and alleles for resistance, resulting in female-biased sex ratios and higher resistance and coercion towards the productive 'core' of the metapopulation. However, steeper environmental gradients, which lead to lower density and less mating at the 'edge', generate female bias without promoting coercion or resistance. Our analysis shows that local adaptation of facultatively parthenogenetic populations subject to sexual conflict and productivity gradients can generate striking spatial variation suggesting new patterns for empirical investigation. These findings could also help to explain the rarity of facultative parthenogenesis in animals.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Population genomics and geographical parthenogenesis in Japanese harvestmen (Opiliones, Sclerosomatidae, Leiobunum)

Naturally occurring population variation in reproductive mode presents an opportunity for researchers to test hypotheses regarding the evolution of sex. Asexual reproduction frequently assumes a geographical pattern, in which parthenogenesis-dominated populations are more broadly dispersed than their sexual conspecifics. We evaluate the geographical distribution of genomic signatures associated with parthenogenesis using nuclear and mitochondrial DNA sequence data from two Japanese harvestman sister taxa, Leiobunum manubriatum and Leiobunum globosum. Asexual reproduction is putatively facultative in these species, and female-biased localities are common in habitat margins. Past karyotypic and current cytometric work indicates L. globosum is entirely tetraploid, while L. manubriatum may be either diploid or tetraploid. We estimated species phylogeny, genetic differentiation, diversity, and mitonuclear discordance in females collected across the species range in order to identify range expansion toward marginal habitat, potential for hybrid origin, and persistence of asexual lineages. Our results point to northward expansion of a tetraploid ancestor of L. manubriatum and L. globosum, coupled with support for greater male gene flow in southern L. manubriatum localities. Specimens from localities in the Tohoku and Hokkaido regions were indistinct, particularly those of L. globosum, potentially due to little mitochondrial differentiation or haplotypic variation. Although L. manubriatum overlaps with L. globosum across its entire range, L. globosum was reconstructed as monophyletic with strong support using mtDNA, and marginal support with nuclear loci. Ultimately, we find evidence for continued sexual reproduction in both species and describe opportunities to clarify the rate and mechanism of parthenogenesis.

General environmental heterogeneity as the explanation of sexuality? Comparative study shows that ancient asexual taxa are associated with both biotically and abiotically homogeneous environments

Ecological theories of sexual reproduction assume that sexuality is advantageous in certain conditions, for example, in biotically or abiotically more heterogeneous environments. Such theories thus could be tested by comparative studies. However, the published results of these studies are rather unconvincing. Here, we present the results of a new comparative study based exclusively on the ancient asexual clades. The association with biotically or abiotically homogeneous environments in these asexual clades was compared with the same association in their sister, or closely related, sexual clades. Using the conservative definition of ancient asexuals (i.e., age >1 million years), we found eight pairs of taxa of sexual and asexual species, six differing in the heterogeneity of their inhabited environment on the basis of available data. The difference between the environmental type associated with the sexual and asexual species was then compared in an exact binomial test. The results showed that the majority of ancient asexual clades tend to be associated with biotically, abiotically, or both biotically and abiotically more homogeneous environments than their sexual controls. In the exploratory part of the study, we found that the ancient asexuals often have durable resting stages, enabling life in subjectively homogeneous environments, live in the absence of intense biotic interactions, and are very often sedentary, inhabiting benthos, and soil. The consequences of these findings for the ecological theories of sexual reproduction are discussed.

Transitions between phases of genomic differentiation during stick-insect speciation

Speciation can involve a transition from a few genetic loci that are resistant to gene flow to genome-wide differentiation. However, only limited data exist concerning this transition and the factors promoting it. Here, we study phases of speciation using data from >100 populations of 11 species of Timema stick insects. Consistent with early phases of genic speciation, adaptive colour-pattern loci reside in localized genetic regions of accentuated differentiation between populations experiencing gene flow. Transitions to genome-wide differentiation are also observed with gene flow, in association with differentiation in polygenic chemical traits affecting mate choice. Thus, intermediate phases of speciation are associated with genome-wide differentiation and mate choice, but not growth of a few genomic islands. We also find a gap in genomic differentiation between sympatric taxa that still exchange genes and those that do not, highlighting the association between differentiation and complete reproductive isolation. Our results suggest that substantial progress towards speciation may involve the alignment of multi-faceted aspects of differentiation.

Responses of holocyclic and anholocyclic Rhopalosiphum padi populations to low-temperature and short-photoperiod induction

The different life cycles of aphid species make these organisms good models for studying the short‐term consequences of sex. The bird cherry‐oat aphid Rhopalosiphum padi has a wide geographic distribution and correspondingly different life cycles. In this study, the life cycles of R. padi collected from six different regions in China were characterized experimentally by comparing the responses of holocyclic and anholocyclic populations to low‐temperature and short‐photoperiod induction. Clones collected from Chuzhou, Taian, and Taigu consistently reproduced via obligate parthenogenesis, whereas clones from Hami and Baicheng were holocyclic in their response, and those from Lanzhou were both holocyclic and anholocyclic. Prolonged exposure to low temperature and a short photoperiod (LS) had negative effects on the offspring of anholocyclic aphids with regard to adult lifespan, total longevity, and fecundity compared with aphids maintained at a normal temperature and a long photoperiod (NL). Holocyclic LS R. padi had longer developmental times at all nymph stages, a shorter adult lifespan, shorter total longevity, and a lower fecundity than NL counterparts. The adult prereproduction period of gynoparae was significantly longer than that of virginoparae, and the total longevity of gynoparae was significantly shorter than that of virginoparae. Moreover, the net reproductive and gross reproduction rates, as well as the total fecundity, were roughly fivefold higher in virginoparae than in gynoparae, indicating that there is the short‐term cost of sex. When maintained on their secondary host (Triticum aestivum), gynoparae, males, and oviparae produced by holocyclic populations could survive, and gynoparae produced oviparae. However, under NL conditions, oviparae could not produce overwintering eggs on the secondary host, whereas a few overwintering eggs were generated by oviparae under LS conditions. Taken together, these results illuminate the complexity of insect responses and contribute to a complete understanding of the aphid life cycle and its evolution.

Nematode endoparasites do not codiversify with their stick insect hosts

Host-parasite coevolution stems from reciprocal selection on host resistance and parasite infectivity, and can generate some of the strongest selective pressures known in nature. It is widely seen as a major driver of diversification, the most extreme case being parallel speciation in hosts and their associated parasites. Here, we report on endoparasitic nematodes, most likely members of the mermithid family, infecting different Timema stick insect species throughout California. The nematodes develop in the hemolymph of their insect host and kill it upon emergence, completely impeding host reproduction. Given the direct exposure of the endoparasites to the host's immune system in the hemolymph, and the consequences of infection on host fitness, we predicted that divergence among hosts may drive parallel divergence in the endoparasites. Our phylogenetic analyses suggested the presence of two differentiated endoparasite lineages. However, independently of whether the two lineages were considered separately or jointly, we found a complete lack of codivergence between the endoparasitic nematodes and their hosts in spite of extensive genetic variation among hosts and among parasites. Instead, there was strong isolation by distance among the endoparasitic nematodes, indicating that geography plays a more important role than host-related adaptations in driving parasite diversification in this system. The accumulating evidence for lack of codiversification between parasites and their hosts at macroevolutionary scales contrasts with the overwhelming evidence for coevolution within populations, and calls for studies linking micro- versus macroevolutionary dynamics in host-parasite interactions.

Horizontal Gene Transfer of Pectinases from Bacteria Preceded the Diversification of Stick and Leaf Insects

Genes acquired by horizontal transfer are increasingly being found in animal genomes. Understanding their origin and evolution requires knowledge about the phylogenetic relationships from both source and recipient organisms. We used RNASeq data and respective assembled transcript libraries to trace the evolutionary history of polygalacturonase (pectinase) genes in stick insects (Phasmatodea). By mapping the distribution of pectinase genes on a Polyneoptera phylogeny, we identified the transfer of pectinase genes from known phasmatodean gut microbes into the genome of an early euphasmatodean ancestor that took place between 60 and 100 million years ago. This transfer preceded the rapid diversification of the suborder, enabling symbiont-free pectinase production that would increase the insects’ digestive efficiency and reduce dependence on microbes. Bacteria-to-insect gene transfer was thought to be uncommon, however the increasing availability of large-scale genomic data may change this prevailing notion.

Sticky Genomes: Using NGS Evidence to Test Hybrid Speciation Hypotheses

Hypotheses of hybrid origin are common. Here we use next generation sequencing to test a hybrid hypothesis for a non-model insect with a large genome. We compared a putative hybrid triploid stick insect species (Acanthoxyla geisovii) with its putative paternal diploid taxon (Clitarchus hookeri), a relationship that provides clear predictions for the relative genetic diversity within each genome. The parental taxon is expected to have comparatively low allelic diversity that is nested within the diversity of the hybrid daughter genome. The scale of genome sequencing required was conveniently achieved by extracting mRNA and sequencing cDNA to examine expressed allelic diversity. This allowed us to test hybrid-progenitor relationships among non-model organisms with large genomes and different ploidy levels. Examination of thousands of independent loci avoids potential problems produced by the silencing of parts of one or other of the parental genomes, a phenomenon sometimes associated with the process of stabilisation of a hybrid genome. Transcript assembles were assessed for evidence of paralogs and/or alternative splice variants before proceeding. Comparison of transcript assemblies was not an appropriate measure of genetic variability, but by mapping reads back to clusters derived from each species we determined levels of allelic diversity. We found greater cDNA sequence diversity among alleles in the putative hybrid species (Acanthoxyla geisovii) than the non-hybrid. The allelic diversity within the putative paternal species (Clitachus hookeri) nested within the hybrid-daughter genome, supports the current view of a hybrid-progenitor relationship for these stick insect species. Next generation sequencing technology provides opportunities for testing evolutionary hypotheses with non-model organisms, including, as here, genomes that are large due to polyploidy.

Extensive variation in chromosome number and genome size in sexual and parthenogenetic species of the jumping-bristletail genus Machilis (Archaeognatha)

Parthenogenesis in animals is often associated with polyploidy and restriction to extreme habitats or recently deglaciated areas. It has been hypothesized that benefits conferred by asexual reproduction and polyploidy are essential for colonizing these habitats. However, while evolutionary routes to parthenogenesis are manifold, study systems including polyploids are scarce in arthropods. The jumping-bristletail genus Machilis (Insecta: Archaeognatha) includes both sexual and parthenogenetic species, and recently, the occurrence of polyploidy has been postulated. Here, we applied flow cytometry, karyotyping, and mitochondrial DNA sequencing to three sexual and five putatively parthenogenetic Eastern-Alpine Machilis species to investigate whether (1) parthenogenesis originated once or multiply and (2) whether parthenogenesis is strictly associated with polyploidy. The mitochondrial phylogeny revealed that parthenogenesis evolved at least five times independently among Eastern-Alpine representatives of this genus. One parthenogenetic species was exclusively triploid, while a second consisted of both diploid and triploid populations. The three other parthenogenetic species and all sexual species were diploid. Our results thus indicate that polyploidy can co-occur with parthenogenesis, but that it was not mandatory for the emergence of parthenogenesis in Machilis. Overall, we found a weak negative correlation of monoploid genome size (Cx) and chromosome base number (x), and this connection is stronger among parthenogenetic species alone. Likewise, monoploid genome size decreased with elevation, and we therefore hypothesize that genome downsizing could have been crucial for the persistence of alpine Machilis species. Finally, we discuss the evolutionary consequences of intraspecific chromosomal rearrangements and the presence of B chromosomes. In doing so, we highlight the potential of Alpine Machilis species for research on chromosomal and genome-size alterations during speciation.

Alteration of sexual reproduction and genetic diversity in the kelp species Laminaria digitata at the southern limit of its range

Adaptation to marginal habitats at species range-limits has often been associated with parthenogenetic reproduction in terrestrial animals and plants. Laboratory observations have shown that brown algae exhibit a high propensity for parthenogenesis by various mechanisms. The kelp Laminaria digitata is an important component of the ecosystem in Northern European rocky intertidal habitats. We studied four L. digitata populations for the effects of marginality on genetic diversity and sexual reproduction. Two populations were marginal: One (Locquirec, in Northern Brittany) was well within the geographic range, but was genetically isolated from other populations by large stretches of sandy beaches. Another population was at the range limits of the species (Quiberon, in Southern Brittany) and was exposed to much higher seasonal temperature changes. Microsatellite analyses confirmed that these populations showed decreased genetic and allelic diversity, consistent with marginality and genetic isolation. Sporophytes from both marginal populations showed greatly diminished spore-production compared to central populations, but only the southern-limit population (Quiberon) showed a high propensity for producing unreduced (2N) spores. Unreduced 2N spores formed phenotypically normal gametophytes with nuclear area consistent with ≥2N DNA contents, and microsatellite studies suggested these were produced at least in part by automixis. However, despite this being the dominant path of spore production in Quiberon sporophyte individuals, the genetic evidence indicated the population was maintained mostly by sexual reproduction. Thus, although spore production and development showed the expected tendency of geographical parthenogenesis in marginal populations, this appeared to be a consequence of maladaptation, rather than an adaptation to, life in a marginal habitat.

Hydrocarbon divergence and reproductive isolation in Timema stick insects

Background

Individuals commonly prefer certain trait values over others when choosing their mates. If such preferences diverge between populations, they can generate behavioral reproductive isolation and thereby contribute to speciation. Reproductive isolation in insects often involves chemical communication, and cuticular hydrocarbons, in particular, serve as mate recognition signals in many species. We combined data on female cuticular hydrocarbons, interspecific mating propensity, and phylogenetics to evaluate the role of cuticular hydrocarbons in diversification of Timema walking-sticks.

Results

Hydrocarbon profiles differed substantially among the nine analyzed species, as well as between partially reproductively-isolated T. cristinae populations adapted to different host plants. In no-choice trials, mating was more likely between species with similar than divergent hydrocarbon profiles, even after correcting for genetic divergences. The macroevolution of hydrocarbon profiles, along a Timema species phylogeny, fits best with a punctuated model of phenotypic change concentrated around speciation events, consistent with change driven by selection during the evolution of reproductive isolation.

Conclusion

Altogether, our data indicate that cuticular hydrocarbon profiles vary among Timema species and populations, and that most evolutionary change in hydrocarbon profiles occurs in association with speciation events. Similarities in hydrocarbon profiles between species are correlated with interspecific mating propensities, suggesting a role for cuticular hydrocarbon profiles in mate choice and speciation in the genus Timema.

Neutral and selection-driven decay of sexual traits in asexual stick insects

Environmental shifts and lifestyle changes may result in formerly adaptive traits becoming non-functional or maladaptive. The subsequent decay of such traits highlights the importance of natural selection for adaptations, yet its causes have rarely been investigated. To study the fate of formerly adaptive traits after lifestyle changes, we evaluated sexual traits in five independently derived asexual lineages, including traits that are specific to males and therefore not exposed to selection. At least four of the asexual lineages retained the capacity to produce males that display normal courtship behaviours and are able to fertilize eggs of females from related sexual species. The maintenance of male traits may stem from pleiotropy, or from these traits only regressing via drift, which may require millions of years to generate phenotypic effects. By contrast, we found parallel decay of sexual traits in females. Asexual females produced altered airborne and contact signals, had modified sperm storage organs, and lost the ability to fertilize their eggs, impeding reversals to sexual reproduction. Female sexual traits were decayed even in recently derived asexuals, suggesting that trait changes following the evolution of asexuality, when they occur, proceed rapidly and are driven by selective processes rather than drift.

Can resource costs of polyploidy provide an advantage to sex?

The predominance of sexual reproduction despite its costs indicates that sex provides substantial benefits, which are usually thought to derive from the direct genetic consequences of recombination and syngamy. While genetic benefits of sex are certainly important, sexual and asexual individuals, lineages, or populations may also differ in physiological and life history traits that could influence outcomes of competition between sexuals and asexuals across environmental gradients. Here, we address possible phenotypic costs of a very common correlate of asexuality, polyploidy. We suggest that polyploidy could confer resource costs related to the dietary phosphorus demands of nucleic acid production; such costs could facilitate the persistence of sex in situations where asexual taxa are of higher ploidy level and phosphorus availability limits important traits like growth and reproduction. We outline predictions regarding the distribution of diploid sexual and polyploid asexual taxa across biogeochemical gradients and provide suggestions for study systems and empirical approaches for testing elements of our hypothesis.

Geographic parthenogenesis and plant-enemy interactions in the common dandelion

BACKGROUND

Many species with sexual and asexual variants show a pattern of geographic parthenogenesis where asexuals have broader and higher-latitude distribution than sexuals. Because sexual reproduction is often considered a costly evolutionary strategy that is advantageous in the face of selection by coevolving pests and pathogens, one possible explanation for geographic parthenogenesis is that populations at higher latitudes are exposed to fewer pests and pathogens. We tested this hypothesis in the common dandelion (Taraxacum officinale), a species with well-established geographic parthenogenesis, by screening prevalence and effects of several specialized pests and pathogens in natural dandelion populations.

RESULTS

We did a population survey of 18 dandelion populations along a geographic transect that ranged from the area where sexual and asexual dandelions co-occur northward into the area where only asexuals occur. In addition we used four southern and four northern populations in a 8x8 cross-inoculation greenhouse experiment in which plants were exposed experimentally to each other's natural field soil microbial communities. The cross-inoculation experiment indicated a higher pathogenicity of soil microbial communities from the southern, mostly sexual, populations compared to soil microbial communities from the northern asexual populations. Northern dandelion populations also showed reduced infestation by a specialized seed-eating weevil. A similar trend of reduced rust fungus infection in northern populations was observed but this trend was not statistically significant.

CONCLUSIONS

The prevalence of pests and pathogens decreased along the south-to-north axis of geographic parthenogenesis. This highlights the potential of biotic interactions in shaping patterns of geographic parthenogenesis.

Sexually antagonistic coevolution for sexual harassment can act as a barrier to further invasions by parthenogenesis

The assumption of a twofold cost of sex not only complicates the maintenance of sex but also sets conditions for sexual conflict: in organisms with the twofold cost, males often sexually harass females. Sexual harassment is detrimental to female fitness and thus might help maintain sexual populations if male harassment inflicts a harsher cost on parthenogens than on sexual females (asymmetric harassment cost). However, the generality of this concept is now considered doubtful because selective harassment of parthenogens results in loss of mating opportunities for males. Using three mathematical models, I show here that sexual harassment still can impose the asymmetric cost on parthenogens. First, I apply the Lotka-Volterra model to show the degree of asymmetric harassment cost that permits sex to be maintained stably in the population. Second, using adaptive dynamics, I examine whether sexually antagonistic coevolution for sexual harassment, which occurs only in sexual populations, can promote the asymmetric harassment cost. Finally, an individual-based model, which assumes a spatial structure unlike that in the other two, demonstrates that the asymmetric evolution of harassment cost prevents further invasions of parthenogens from different patches into sexual lineages; these mechanisms may account for allopatric distributions of sexual and parthenogenetic lineages as well as the maintenance of sex.

Deleterious mutation accumulation in asexual Timema stick insects

Sexual reproduction is extremely widespread in spite of its presumed costs relative to asexual reproduction, indicating that it must provide significant advantages. One postulated benefit of sex and recombination is that they facilitate the purging of mildly deleterious mutations, which would accumulate in asexual lineages and contribute to their short evolutionary life span. To test this prediction, we estimated the accumulation rate of coding (nonsynonymous) mutations, which are expected to be deleterious, in parts of one mitochondrial (COI) and two nuclear (Actin and Hsp70) genes in six independently derived asexual lineages and related sexual species of Timema stick insects. We found signatures of increased coding mutation accumulation in all six asexual Timema and for each of the three analyzed genes, with 3.6- to 13.4-fold higher rates in the asexuals as compared with the sexuals. In addition, because coding mutations in the asexuals often resulted in considerable hydrophobicity changes at the concerned amino acid positions, coding mutations in the asexuals are likely associated with more strongly deleterious effects than in the sexuals. Our results demonstrate that deleterious mutation accumulation can differentially affect sexual and asexual lineages and support the idea that deleterious mutation accumulation plays an important role in limiting the long-term persistence of all-female lineages.

Diversity in the reproductive modes of European Daphnia pulicaria deviates from the geographical parthenogenesis

BACKGROUND

Multiple transitions to obligate parthenogenesis have occurred in the Daphnia pulex complex in North America. These newly formed asexual lineages are differentially distributed being found predominantly at high latitudes. This conforms to the rule of geographical parthenogenesis postulating prevalence of asexuals at high latitudes and altitudes. While the reproductive mode of high-latitude populations is relatively well studied, little is known about the reproduction mode in high altitudes. This study aimed to assess the reproductive mode of Daphnia pulicaria, a species of the D. pulex complex, from high altitude lakes in Europe.

METHODOLOGY/PRINCIPAL FINDINGS

Variation at eight microsatellite loci revealed that D. pulicaria from the High Tatra Mountains (HTM) had low genotype richness and showed excess of heterozygotes and significant deviations from Hardy-Weinberg expectations, and was thus congruent with reproduction by obligate parthenogenesis. By contrast, populations from the Pyrenees (Pyr) were generally in Hardy-Weinberg equilibrium and had higher genotypic richness, suggesting that they are cyclic parthenogens. Four lakes from lowland areas (LLaP) had populations with an uncertain or mixed breeding mode. All D. pulicaria had mtDNA ND5 haplotypes of the European D. pulicaria lineage. Pyr were distinct from LLaP and HTM at the ND5 gene. By contrast, HTM shared two haplotypes with LLaP and one with Pyr. Principal Coordinate Analysis of the microsatellite data revealed clear genetic differentiation into three groups. HTM isolates were intermediate to Pyr and LLaP, congruent with a hybrid origin.

CONCLUSION/SIGNIFICANCE

Inferred transitions to obligate parthenogenesis have occurred only in HTM, most likely as a result of hybridizations. In contrast to North American populations, these transitions do not appear to involve meiosis suppressor genes and have not been accompanied by polyploidy. The absence of obligate parthenogenesis in Pyr, an environment highly similar to the HTM, may be due to the lack of opportunities for hybridization.

Geographic parthenogenesis in a consumer-resource model for sexual reproduction

The phenomenon of geographic parthenogenesis is closely tied to the question of why sexual reproduction is the dominant mode of reproduction in animals and plants. Geographic parthenogenesis describes the fact that many species reproduce asexually at the boundaries of their range. We present a mathematical model that derives the dominance of sexuals at the center and the dominance of asexuals at the boundary of a species' range from exactly the same mechanism. Our model is based on a set of resources that regrow slowly and that can be consumed only by those individuals that have a suitable genotype. Genotype is implemented by a multilocus model with two alleles at each locus, and with free recombination during production of sexual offspring. The model is tailored to seasonal species with intermittent mixis and low survival of offspring, such as Daphnia and aphids. Several patches of resources are arranged in a row, with a gradient of those parameters that typically vary through the range of species. By letting sexually and asexually reproducing populations compete, we obtain the typical patterns of geographic parthenogenesis.

Evolutionary origin and phylogeography of the diploid obligate parthenogen Artemia parthenogenetica (Branchiopoda: Anostraca)

BACKGROUND

Understanding the evolutionary origin and the phylogeographic patterns of asexual taxa can shed light on the origin and maintenance of sexual reproduction. We assessed the geographic origin, genetic diversity, and phylogeographic history of obligate parthenogen diploid Artemia parthenogenetica populations, a widespread halophilic crustacean.

METHODOLOGY/PRINCIPAL FINDINGS

We analysed a partial sequence of the Cytochrome c Oxidase Subunit I mitochondrial gene from an extensive set of localities (including Eurasia, Africa, and Australia), and examined their phylogeographic patterns and the phylogenetic relationships of diploid A. parthenogenetica and its closest sexual relatives. Populations displayed an extremely low level of mitochondrial genetic diversity, with one widespread haplotype shared by over 79% of individuals analysed. Phylogenetic and phylogeographic analyses indicated a multiple and recent evolutionary origin of diploid A. parthenogenetica, and strongly suggested that the geographic origin of parthenogenesis in Artemia was in Central Asia. Our results indicate that the maternal sexual ancestors of diploid A. parthenogenetica were an undescribed species from Kazakhstan and A. urmiana.

CONCLUSIONS/SIGNIFICANCE

We found evidence for multiple origin of parthenogenesis in Central Asia. Our results indicated that, shortly after its origin, diploid A. parthenogenetica populations underwent a rapid range expansion from Central Asia towards the Mediterranean region, and probably to the rest of its current geographic distribution. This contrasts with the restricted geographic distribution, strong genetic structure, and regional endemism of sexual Artemia lineages and other passively dispersed sexual continental aquatic invertebrates. We hypothesize that diploid parthenogens might have reached their current distribution in historical times, with a range expansion possibly facilitated by an increased availability of suitable habitat provided by anthropogenic activities, such as the spread of solar saltworks, aided by their natural dispersal vectors (i.e., waterbirds).

Effects of plant sex on range distributions and allocation to reproduction

Despite an abundance of theory, few empirical studies have explored the ecological and evolutionary consequences of sex. We used a comparative phylogenetic approach to examine whether transitions between sexual and asexual reproduction are associated with changes in the size and distribution of species' geographical ranges, and their investment in reproduction. Here, we reconstructed the phylogeny of the genus Oenothera sections Oenothera and Calylophus (Onagraceae), which contain 35 sexual and 30 functionally asexual species. From each species, we collected data on the geographical distribution and variation in plant traits related to reproduction. Functionally asexual species occurred at higher latitudes, but did not differ in range size, compared with sexual species. Transitions to asexuality were associated with decreased investment in floral structures, including the length of petals, floral tubes and styles. Decreased anther size and increased seed size within asexual species also suggest altered allocation to male and female fitness. The observed range shifts are consistent with superior colonization of environments by asexual species following glaciation, and the observed changes in reproductive allocation support predictions made by models relating to the evolution of selfing. Our results suggest that the evolutionary consequences of asexual reproduction might be less restrictive than previously thought.

Cytotype stability, facultative apomixis and geographical parthenogenesis in Ranunculus kuepferi (Ranunculaceae)

BACKGROUND AND AIMS

Asexual organisms are more widespread in previously glaciated areas than their sexual relatives ('geographical parthenogenesis'). In plants, this pattern is probably dependent on reproductive isolation and stability of cytotypes within their respective distribution areas. Both partial apomixis and introgressive hybridization potentially destabilize the spatial separation of sexual and apomictic populations. The wide distribution of apomicts may be further enhanced by uniparental reproduction which is advantageous for colonization. These factors are studied in the alpine species Ranunculus kuepferi.

METHODS

Geographical distribution, diversity and mode of reproduction of cytotypes were assessed using flow cytometry and flow cytometric seed screening on samples from 59 natural populations of Ranunculus kuepferi. Seed set of cytotypes was compared in the wild.

KEY RESULTS

Diploid sexuals are confined to the south-western parts of the Alps, while tetraploid apomicts dominate in previously glaciated and in geographically isolated areas despite a significantly lower fertility. Other cytotypes (3x, 5x and 6x) occur mainly in the sympatric zone, but without establishing populations. The tetraploids are predominantly apomictic, but also show a partial apomixis via an uncoupling of apomeiosis and parthenogenesis in the seed material. Both pseudogamy and autonomous endosperm formation are observed which may enhance uniparental reproduction.

CONCLUSIONS

Diploids occupy a glacial relic area and resist introgression of apomixis, probably because of a significantly higher seed set. Among the polyploids, only apomictic tetraploids form stable populations; the other cytotypes arising from partial apomixis fail to establish, probably because of minority cytotype disadvantages. Tetraploid apomicts colonize previously devastated and also distant areas via long-distance dispersal, confirming Baker's law of an advantage of uniparental reproduction. It is concluded that stability of cytotypes and of modes of reproduction are important factors for establishing a pattern of geographical parthenogenesis.