The mitochondrial cytochrome oxidase II gene in Bacillus stick insects: ancestry of hybrids, androgenesis, and phylogenetic relationships

Parthenogenetic Stick Insects Exhibit Signatures of Preservation in the Molecular Architecture of Male Reproduction

After the loss of a trait, theory predicts that the molecular machinery underlying its phenotypic expression should decay. Yet, empirical evidence is contrasting. Here, we test the hypotheses that (i) the molecular ground plan of a lost trait could persist due to pleiotropic effects on other traits and (ii) that gene co-expression network architecture could constrain individual gene expression. Our testing ground has been the Bacillus stick insect species complex, which contains close relatives that are either bisexual or parthenogenetic. After the identification of genes expressed in male reproductive tissues in a bisexual species, we investigated their gene co-expression network structure in two parthenogenetic species. We found that gene co-expression within the male gonads was partially preserved in parthenogens. Furthermore, parthenogens did not show relaxed selection on genes upregulated in male gonads in the bisexual species. As these genes were mostly expressed in female gonads, this preservation could be driven by pleiotropic interactions and an ongoing role in female reproduction. Connectivity within the network also played a key role, with highly connected-and more pleiotropic-genes within male gonad also having a gonad-biased expression in parthenogens. Our findings provide novel insight into the mechanisms which could underlie the production of rare males in parthenogenetic lineages; more generally, they provide an example of the cryptic persistence of a lost trait molecular architecture, driven by gene pleiotropy on other traits and within their co-expression network.

Gene transcriptional profiles in gonads of Bacillus taxa (Phasmida) with different cytological mechanisms of automictic parthenogenesis

The evolution of automixis - i.e., meiotic parthenogenesis - requires several features, including ploidy restoration after meiosis and maintenance of fertility. Characterizing the relative contribution of novel versus pre-existing genes and the similarities in their expression and sequence evolution is fundamental to understand the evolution of reproductive novelties. Here we identify gonads-biased genes in two Bacillus automictic stick-insects and compare their expression profile and sequence evolution with a bisexual congeneric species. The two parthenogens restore ploidy through different cytological mechanisms: in Bacillus atticus, nuclei derived from the first meiotic division fuse to restore a diploid egg nucleus, while in Bacillus rossius, diploidization occurs in some cells of the haploid blastula through anaphase restitution. Parthenogens' gonads transcriptional program is found to be largely assembled from genes that were already present before the establishment of automixis. The three species transcriptional profiles largely reflect their phyletic relationships, yet we identify a shared core of genes with gonad-biased patterns of expression in parthenogens which are either male gonads-biased in the sexual species or are not differentially expressed there. At the sequence level, just a handful of gonads-biased genes were inferred to have undergone instances of positive selection exclusively in the parthenogen species. This work is the first to explore the molecular underpinnings of automixis in a comparative framework: it delineates how reproductive novelties can be sustained by genes whose origin precedes the establishment of the novelty itself and shows that different meiotic mechanisms of reproduction can be associated with a shared molecular ground plan.

A second view on the evolution of flight in stick and leaf insects (Phasmatodea)

BACKGROUND

The re-evolution of complex characters is generally considered impossible, yet, studies of recent years have provided several examples of phenotypic reversals shown to violate Dollo's law. Along these lines, the regain of wings in stick and leaf insects (Phasmatodea) was hypothesised to have occurred several times independently after an ancestral loss, a scenario controversially discussed among evolutionary biologists due to overestimation of the potential for trait reacquisition as well as to the lack of taxonomic data.

RESULTS

We revisited the recovery of wings by reconstructing a phylogeny based on a comprehensive taxon sample of over 500 representative phasmatodean species to infer the evolutionary history of wings. We additionally explored the presence of ocelli, the photoreceptive organs used for flight stabilisation in winged insects, which might provide further information for interpreting flight evolution. Our findings support an ancestral loss of wings and that the ancestors of most major lineages were wingless. While the evolution of ocelli was estimated to be dependent on the presence of (fully-developed) wings, ocelli are nevertheless absent in the majority of all examined winged species and only appear in the members of few subordinate clades, albeit winged and volant taxa are found in every euphasmatodean lineage.

CONCLUSION

In this study, we explored the evolutionary history of wings in Phasmatodea and demonstrate that the disjunct distribution of ocelli substantiates the hypothesis on their regain and thus on trait reacquisition in general. Evidence from the fossil record as well as future studies focussing on the underlying genetic mechanisms are needed to validate our findings and to further assess the evolutionary process of phenotypic reversals.

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.

Androgenesis: where males hijack eggs to clone themselves

Androgenesis is a form of quasi-sexual reproduction in which a male is the sole source of the nuclear genetic material in the embryo. Two types of androgenesis occur in nature. Under the first type, females produce eggs without a nucleus and the embryo develops from the male gamete following fertilization. Evolution of this type of androgenesis is poorly understood as the parent responsible for androgenesis (the mother) gains no benefit from it. Ultimate factors driving the evolution of the second type of androgenesis are better understood. In this case, a zygote is formed between a male and a female gamete, but the female genome is eliminated. When rare, androgenesis with genome elimination is favoured because an androgenesis-determining allele has twice the reproductive success of an allele that determines sexual reproduction. Paradoxically, except in hermaphrodites, a successful androgenetic strain can drive such a male-biased sex ratio that the population goes extinct. This likely explains why androgenesis with genome elimination appears to be rarer than androgenesis via non-nucleate eggs, although both forms are either very rare or remain largely undetected in nature. Nonetheless, some highly invasive species including ants and freshwater clams are androgenetic, for reasons that are largely unexplained.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Obligatory parthenogenesis and TE load: Bacillus stick insects and the R2 non-LTR retrotransposon

Transposable elements (TEs) are selfish genetic elements whose self-replication is contrasted by the host genome. In this context, host reproductive strategies are predicted to impact on both TEs load and activity. The presence and insertion distribution of the non-LTR retrotransposon R2 was here studied in populations of the strictly bisexual Bacillus grandii maretimi and of the obligatory parthenogenetic Bacillus atticus atticus. Furthermore, data were also obtained from the offspring of selected B. a. atticus females. At the population level, the gonochoric B. g. maretimi showed a significantly higher R2 load than the obligatory parthenogenetic B. a. atticus. The comparison with bisexual and unisexual Bacillus rossius populations showed that their values were higher than those recorded for B. a. atticus and similar, or even higher, than those of B. g. maretimi. Consistently, an R2 load reduction is scored in B. a. atticus offspring even if with a great variance. On the whole, data here produced indicate that in the obligatory unisexual B. a. atticus R2 is active and that mechanisms of molecular turnover are effective. Furthermore, progeny analyses show that, at variance of the facultative parthenogenetic B. rossius, the R2 activity is held at a lower rate. Modeling parental-offspring inheritance, suggests that in B. a. atticus recombination plays a major role in eliminating insertions rather than selection, as previously suggested for unisexual B. rossius progeny, even if in both cases a high variance is observed. In addition to this, mechanisms of R2 silencing or chances of clonal selection cannot be ruled out.

Hybridogenesis and a potential case of R2 non-LTR retrotransposon horizontal transmission in Bacillus stick insects (Insecta Phasmida)

Horizontal transfer (HT) is an event in which the genetic material is transferred from one species to another, even if distantly related, and it has been demonstrated as a possible essential part of the lifecycle of transposable elements (TEs). However, previous studies on the non-LTR R2 retrotransposon, a metazoan-wide distributed element, indicated its vertical transmission since the Radiata-Bilateria split. Here we present the first possible instances of R2 HT in stick insects of the genus Bacillus (Phasmida). Six R2 elements were characterized in the strictly bisexual subspecies B. grandii grandii, B. grandii benazzii and B. grandii maretimi and in the obligatory parthenogenetic taxon B. atticus. These elements were compared with those previously retrieved in the facultative parthenogenetic species B. rossius. Phylogenetic inconsistencies between element and host taxa, and age versus divergence analyses agree and support at least two HT events. These HT events can be explained by taking into consideration the complex Bacillus reproductive biology, which includes also hybridogenesis, gynogenesis and androgenesis. Through these non-canonical reproductive modes, R2 elements may have been transferred between Bacillus genomes. Our data suggest, therefore, a possible role of hybridization for TEs survival and the consequent reshaping of involved genomes.

Dead element replicating: degenerate R2 element replication and rDNA genomic turnover in the Bacillus rossius stick insect (Insecta: Phasmida)

R2 is an extensively investigated non-LTR retrotransposon that specifically inserts into the 28S rRNA gene sequences of a wide range of metazoans, disrupting its functionality. During R2 integration, first strand synthesis can be incomplete so that 5’ end deleted copies are occasionally inserted. While active R2 copies repopulate the locus by retrotransposing, the non-functional truncated elements should frequently be eliminated by molecular drive processes leading to the concerted evolution of the rDNA array(s). Although, multiple R2 lineages have been discovered in the genome of many animals, the rDNA of the stick insect Bacillus rossius exhibits a peculiar situation: it harbors both a canonical, functional R2 element (R2Br^fun) as well as a full-length but degenerate element (R2Br^deg). An intensive sequencing survey in the present study reveals that all truncated variants in stick insects are present in multiple copies suggesting they were duplicated by unequal recombination. Sequencing results also demonstrate that all R2Br^deg copies are full-length, i. e. they have no associated 5' end deletions, and functional assays indicate they have lost the active ribozyme necessary for R2 RNA maturation. Although it cannot be completely ruled out, it seems unlikely that the degenerate elements replicate via reverse transcription, exploiting the R2Br^fun element enzymatic machinery, but rather via genomic amplification of inserted 28S by unequal recombination. That inactive copies (both R2Br^deg or 5'-truncated elements) are not eliminated in a short term in stick insects contrasts with findings for the Drosophila R2, suggesting a widely different management of rDNA loci and a lower efficiency of the molecular drive while achieving the concerted evolution.

Evolutionary dynamics of R2 retroelement and insertion inheritance in the genome of bisexual and parthenogenetic Bacillus rossius populations (Insecta Phasmida)

Theoretical and empirical studies have shown differential management of transposable elements in organisms with different reproductive strategies. To investigate this issue, we analysed the R2 retroelement structure and variability in parthenogenetic and bisexual populations of Bacillus rossius stick insects, as well as insertions inheritance in the offspring of parthenogenetic isolates and of crosses. The B. rossius genome hosts a functional (R2Br(fun) ) and a degenerate (R2Br(deg) ) element, their presence correlating with neither reproductive strategies nor population distribution. The median-joining network method indicated that R2Br(fun) duplicates through a multiple source model, while R2Br(deg) is apparently still duplicating via a master gene model. Offspring analyses showed that unisexual and bisexual offspring have a similar number of R2Br-occupied sites. Multiple or recent shifts from gonochoric to parthenogenetic reproduction may explain the observed data. Moreover, insertion frequency spectra show that higher-frequency insertions in unisexual offspring significantly outnumber those in bisexual offspring. This suggests that unisexual offspring eliminate insertions with lower efficiency. A comparison with simulated insertion frequencies shows that inherited insertions in unisexual and bisexual offspring are significantly different from the expectation. On the whole, different mechanisms of R2 elimination in unisexual vs bisexual offspring and a complex interplay between recombination effectiveness, natural selection and time can explain the observed data.

Random DNA libraries from three species of the stick insect genus Bacillus (Insecta: Phasmida): repetitive DNA characterization and first observation of polyneopteran MITEs

The repetitive DNA content of the stick insect species Bacillus rossius (facultative parthenogenetic), Bacillus grandii (gonochoric), and Bacillus atticus (obligate parthenogenetic) was analyzed through the survey of random genomic libraries roughly corresponding to 0.006% of the genome. By repeat masking, 19 families of transposable elements were identified (two LTR and six non-LTR retrotransposons; 11 DNA transposons). Moreover, a de novo analysis revealed, among the three libraries, the first MITE family observed in polyneopteran genomes. On the whole, transposable element abundance represented 23.3% of the genome in B. rossius, 22.9% in B. atticus, and 18% in B. grandii. Tandem repeat content in the three libraries is much lower: 1.32%, 0.64%, and 1.86% in B. rossius, B. grandii, and B. atticus, respectively. Microsatellites are the most abundant in all species. Minisatellites were only found in B. rossius and B. atticus, and five monomers belonging to the Bag320 satellite family were detected in B. atticus. Assuming the survey provides adequate representation of the relative genome, the obligate parthenogenetic species (B. atticus), compared with the other two species analyzed, does not show a lower transposable element content, as expected from some theoretical and empirical studies.

Do clones degenerate over time? Explaining the genetic variability of asexuals through population genetic models

Background

Quest for understanding the nature of mechanisms governing the life span of clonal organisms lasts for several decades. Phylogenetic evidence for recent origins of most clones is usually interpreted as proof that clones suffer from gradual age-dependent fitness decay (e.g. Muller's ratchet). However, we have shown that a neutral drift can also qualitatively explain the observed distribution of clonal ages. This finding was followed by several attempts to distinguish the effects of neutral and non-neutral processes. Most recently, Neiman et al. 2009 (Ann N Y Acad Sci.:1168:185-200.) reviewed the distribution of asexual lineage ages estimated from a diverse array of taxa and concluded that neutral processes alone may not explain the observed data. Moreover, the authors inferred that similar types of mechanisms determine maximum asexual lineage ages in all asexual taxa. In this paper we review recent methods for distinguishing the effects of neutral and non-neutral processes and point at methodological problems related with them.

Results and Discussion

We found that contemporary analyses based on phylogenetic data are inadequate to provide any clear-cut answer about the nature and generality of processes affecting evolution of clones. As an alternative approach, we demonstrate that sequence variability in asexual populations is suitable to detect age-dependent selection against clonal lineages. We found that asexual taxa with relatively old clonal lineages are characterised by progressively stronger deviations from neutrality.

Conclusions

Our results demonstrate that some type of age-dependent selection against clones is generally operational in asexual animals, which cover a wide taxonomic range spanning from flatworms to vertebrates. However, we also found a notable difference between the data distribution predicted by available models of sequence evolution and those observed in empirical data. These findings point at the possibility that processes affecting clonal evolution differ from those described in recent studies, suggesting that theoretical models of asexual populations must evolve to address this problem in detail.

Reviewers

This article was reviewed by Isa Schön (nominated by John Logsdon), Arcady Mushegian and Timothy G. Barraclough (nominated by Laurence Hurst).

The mitochondrial genome of Bacillus stick insects (Phasmatodea) and the phylogeny of orthopteroid insects

The Order Phasmatodea (stick and leaf insects) includes many well-known species of cryptic phytophagous insects. In this work, we sequenced the almost complete mitochondrial genomes of two stick insect species of the genus Bacillus. Phasmatodea pertain to the Polyneoptera, and represent one of the major clades of heterometabolous insects. Orthopteroid insect lineages arose through rapid evolutionary radiation events, which likely blurred the phylogenetic reconstructions obtained so far; we therefore performed a phylogenetic analysis to resolve and date all major splits of orthopteroid phylogeny, including the relationships between Phasmatodea and other polyneopterans. We explored several molecular models, with special reference to data partitioning, to correctly detect any phylogenetic signal lying in rough data. Phylogenetic Informativeness analysis showed that the maximum resolving power on the orthopteroid mtDNA dataset is expected for the Upper Cretaceous, about 80millionyears ago (Mya), but at least 70% of the maximum informativeness is also expected for the 150-200 Mya timespan, which makes mtDNA a suitable marker to study orthopteroid splits. A complete chronological calibration has also been computed following a Penalized Likelihood method. In summary, our analysis confirmed the monophyly of Phasmatodea, Dictyoptera and Orthoptera, and retrieved Mantophasmatodea as sister group of Phasmatodea. The origin of orthopteroid insects was also estimated to be in the Middle Triassic, while the order Phasmatodea seems to appear in the Upper Jurassic. The obtained results evidenced that mtDNA is a suitable marker to unravel the ancient splits leading to the orthopteroid orders, given a proper methodological approach.

Reticulate evolution in stick insects: the case of Clonopsis (Insecta Phasmida)

Background

Phasmids show noteworthy abilities to overcome species-specific reproductive isolation mechanisms, including hybridization, polyploidy, parthenogenesis, hybridogenesis and androgenesis. From an evolutionary standpoint, such tangled reproductive interactions lead to the complex phyletic relationships known as "reticulate evolution". Moroccan stick insects of the genus Clonopsis include one bisexual (C. felicitatis) and two closely related parthenogenetic forms (C. gallica, C. soumiae), which represent a polyploid series in chromosome number, but with apparent diploid karyotypes. Moreover, two Clonopsis strains of ameiotic males have been described, C. androgenes-35 and C. androgenes-53. As a consequence, Clonopsis stick insects may have experienced complex micro-evolutionary events, which we try to disentangle in this study.

Results

Mitochondrial cox2 analysis supports a recent divergence of Clonopsis, while AFLPs evidence genetic differentiation not linked to karyotypes, so that parthenogenetic C. gallica and C. soumiae appear to be a mix of strains of polyphyletic origin rather than single parthenogenetic species. Moreover, an admixed hybrid origin seems to be confirmed for C. androgenes.

Conclusion

On the whole, Clonopsis is an intriguing case of reticulate evolution. Actually, complex cladogenetic events should be taken into account to explain the observed genetic structure, including diploidization of polyploid karyotypes, possibly coupled with hybridization and androgenesis. We also proposed a "working hypothesis" to account for the observed data, which deserves further studies, but fits the observed data very well.

Thelytokous parthenogenesis, male clonality and genetic caste determination in the little fire ant: new evidence and insights from the lab

Previous studies indicate that some populations of the little fire ant, Wasmannia auropunctata, display an unusual reproduction system polymorphism. Although some populations have a classical haplodiploid reproduction system, in other populations queens are produced by thelytokous parthenogenesis, males are produced by a male clonality system and workers are produced sexually. An atypical genetic caste determination system was also suggested. However, these conclusions were indirectly inferred from genetic studies on field population samples. Here we set up experimental laboratory nests that allow the control of the parental relationships between individuals. The queens heading those nests originated from either putatively clonal or sexual populations. We characterized the male, queen and worker offspring they produced at 12 microsatellite loci. Our results unambiguously confirm the unique reproduction system polymorphism mentioned above and that male clonality is strictly associated with thelytokous parthenogenesis. We also observed direct evidence of the rare production of sexual gynes and arrhenotokous males in clonal populations. Finally, we obtained evidence of a genetic basis for caste determination. The evolutionary significance of the reproduction system polymorphism and genetic caste determination as well as future research opportunities are discussed.

Twigs on the tree of life? Neutral and selective models for integrating macroevolutionary patterns with microevolutionary processes in the analysis of asexuality

Neutral models characterize evolutionary or ecological patterns expected in the absence of specific causal processes, such as natural selection or ecological interactions. In this study, we describe and evaluate three neutral models that can, in principle, help to explain the apparent 'twigginess' of asexual lineages on phylogenetic trees without involving the negative consequences predicted for the absence of recombination and genetic exchange between individuals. Previously, such phylogenetic twiggyness of asexual lineages has been uncritically interpreted as evidence that asexuality is associated with elevated extinction rates and thus represents an evolutionary dead end. Our first model uses simple phylogenetic simulations to illustrate that, with sexual reproduction as the ancestral state, low transition rates to stable asexuality, or low rates of ascertained 'speciation' in asexuals, can generate twiggy distributions of asexuality, in the absence of high extinction rates for asexual lineages. The second model, developed by Janko et al. (2008), shows that a dynamic equilibrium between origins and neutral losses of asexuals can, under some conditions, generate a relatively low mean age of asexual lineages. The third model posits that the risk of extinction for asexual lineages may be higher than that of sexuals simply because asexuals inhabit higher latitudes or altitudes, and not due to effects of their reproductive systems. Such neutral models are useful in that they allow quantitative evaluation of whether empirical data, such as phylogenetic and phylogeographic patterns of sex and asexuality, indeed support the idea that asexually reproducing lineages persist over shorter evolutionary periods than sexual lineages, due to such processes as mutation accumulation, slower rates of adaptive evolution, or relatively lower levels of genetic variability.

The IGS-ETS in Bacillus (Insecta Phasmida): molecular characterization and the relevance of sex in ribosomal DNA evolution

Background

DNA encoding for ribosomal RNA (rDNA) is arranged in tandemly-repeated subunits, each containing ribosomal genes and non-coding spacers. Because tandemly-repeated, rDNA evolves under a balanced influence of selection and "concerted evolution", which homogenizes rDNA variants over the genome (through genomic turnover mechanisms) and the population (through sexuality).

Results

In this paper we analyzed the IGS-ETS of the automictic parthenogen Bacillus atticus and the bisexual B. grandii, two closely related stick-insect species. Both species share the same IGS-ETS structure and sequence, including a peculiar head-to-tail array of putative transcription enhancers, here named Bag530. Sequence variability of both IGS-ETS and Bag530 evidenced a neat geographic and subspecific clustering in B. grandii, while B. atticus shows a little but evident geographic structure. This was an unexpected result, since the parthenogen B. atticus should lack sequence fixation through sexuality. In B. atticus a new variant might spread in a given geographic area through colonization by an all-female clone, but we cannot discard the hypothesis that B. atticus was actually a bisexual taxon in that area at the time the new variant appeared. Moreover, a gene conversion event between two Bag530 variants of B. grandii benazzii and B. grandii maretimi suggested that rRNA might evolve according to the so-called "library hypothesis" model, through differential amplification of rDNA variants in different taxa.

Conclusion

On the whole, Bacillus rDNA evolution appears to be under a complex array of interacting mechanisms: homogenization may be achieved through genomic turnover that stabilizes DNA-binding protein interactions but, simultaneously, new sequence variants can be adopted, either by direct appearance of newly mutated repeats, or by competition among repeats, so that both DNA-binding proteins and repeat variants drive each other's evolution. All this, coupled with chromosome reshuffling due to sexuality (when present), might drive a quick fixation of new rDNA variants in the populations.

TheLeptynia hispanicaspecies complex (Insecta Phasmida): polyploidy, parthenogenesis, hybridization and more

The Leptynia hispanica stick insect species complex includes bisexuals, triploid and tetraploid parthenogenetic populations, suggesting that polyploidy has played a central role in the evolution of this complex. An analysis of karyotype, mitochondrial DNA (cox2) and nuclear DNA (ef1-alpha) markers was carried out to clarify phylogenetic relationships and microevolutionary/phylogeographical patterns of the L. hispanica complex. Our analyses suggested a subdivision of bisexual populations into four groups, tentatively proposed as incipient species. Moreover, triploids and tetraploids showed two independent origins, the latter being more ancient than the former. From ef1-alpha analysis, triploids showed hybrid constitution, while the hybrid constitution of tetraploids is likely, but more data are needed. We suggest that L. hispanica is a case of 'geographical parthenogenesis' with parthenogenetic strains colonizing large peripheral ranges, and bisexuals confined to glacial refuge areas. Moreover, the age, wide distribution and competitive advantage of polyploids over diploids, demonstrate their significance in the evolution of the L. hispanica species complex.

Sex and Clonality in the Little Fire Ant

Reproduction systems are controlling the creation of new genetic variants as well as how natural selection can operate on these variants. Therefore, they had historically been one of the main foci of evolutionary biology studies. The little fire ant, Wasmannia auropunctata, has been found to display an extraordinary reproduction system, in which both males and female queens are produced clonally. So far, native sexual populations of W. auropunctata have not been identified. Our goals were to identify such sexual populations and investigate the origins of female parthenogenesis and male clonality. Using mitochondrial DNA and microsatellite markers in 17 native populations, we found that traditional sexual populations occurred in W. auropunctata and are likely the recent source of neighboring clonal populations. Queen parthenogenesis has probably evolved several times through mutational events. Male clonality is tightly linked to queen parthenogenesis and thus appears to be female controlled. Its origin could be accounted for by 2 mutually exclusive hypotheses: either by the expected coevolution of the 2 sexes (i.e., a variant of the maternal genome elimination hypothesis) or by a shared mechanistic origin (i.e., by the production of anucleate ovules by parthenogenetic queens). Our results also show that W. auropunctata males and females do not form separate evolutionary units and are unlikely to be engaged in an all-out battle of sexes. This work opens up new perspectives for studies on the adaptive significance and evolutionary stability of mixed sexual and clonal reproduction systems in living organisms.

The first fossil leaf insect: 47 million years of specialized cryptic morphology and behavior

Stick and leaf insects (insect order Phasmatodea) are represented primarily by twig-imitating slender forms. Only a small percentage ( approximately 1%) of extant phasmids belong to the leaf insects (Phylliinae), which exhibit an extreme form of morphological and behavioral leaf mimicry. Fossils of phasmid insects are extremely rare worldwide. Here we report the first fossil leaf insect, Eophyllium messelensis gen. et sp. nov., from 47-million-year-old deposits at Messel in Germany. The new specimen, a male, is exquisitely preserved and displays the same foliaceous appearance as extant male leaf insects. Clearly, an advanced form of extant angiosperm leaf mimicry had already evolved early in the Eocene. We infer that this trait was combined with a special behavior, catalepsy or "adaptive stillness," enabling Eophyllium to deceive visually oriented predators. Potential predators reported from the Eocene are birds, early primates, and bats. The combination of primitive and derived characters revealed by Eophyllium allows the determination of its exact phylogenetic position and illuminates the evolution of leaf mimicry for this insect group. It provides direct evidence that Phylliinae originated at least 47 Mya. Eophyllium enlarges the known geographical range of Phylliinae, currently restricted to southeast Asia, which is apparently a relict distribution. This fossil leaf insect bears considerable resemblance to extant individuals in size and cryptic morphology, indicating minimal change in 47 million years. This absence of evolutionary change is an outstanding example of morphological and, probably, behavioral stasis.

Molecular phylogeny, historical biogeography, and divergence time estimates for swallowtail butterflies of the genus Papilio (Lepidoptera: Papilionidae)

Swallowtail butterflies are recognized as model organisms in ecology, evolutionary biology, genetics, and conservation biology but present numerous unresolved phylogenetic problems. We inferred phylogenetic relationships for 51 of about 205 species of the genus Papilio (sensu lato) from 3.3-Kilobase (kb) sequences of mitochondrial and nuclear DNA (2.3 kb of cytochrome oxidases I and II and 1.0 kb of elongation factor 1 alpha). Congruent phylogenetic trees were recovered within Papilio from analyses of combined data using maximum likelihood, Bayesian analysis, and maximum parsimony bootstrap consensus. Several disagreements with the traditional classification of Papilio were found. Five major previously hypothesized subdivisions within Papilio were well supported: Heraclides, Pterourus, Chilasa, Papilio (sensu stricto), and Eleppone. Further studies are required to clarify relationships within traditional "Princeps," which was paraphyletic. Several biologically interesting characteristics of Papilio appear to have polyphyletic origins, including mimetic adults, larval host associations, and larval morphology. Early diversification within Papilio is estimated at 55-65 million years ago based on a combination of biogeographic time constraints rather than fossils. This divergence time suggests that Papilio has slower apparent substitution rates than do Drosophila and fig-pollinating wasps and/or divergences corrected using best-fit substitution models are still being consistently underestimated. The amount of sequence divergence between Papilio subdivisions is equivalent to divergences between genera in other tribes of the Papilionidae, and between genera of moths of the noctuid subfamily Heliothinae.

Phylogeny and karyotype evolution of the Iberian Leptynia attenuata species complex (Insecta Phasmatodea)

An in-depth analysis of the Leptynia attenuata species complex has been performed by cytochrome oxidase subunit 2 (cox2) gene sequencing as well as karyotype and allozyme analysis. The whole set of data allows to largely resolve the taxonomy of the group and suggests an overall trend of chromosomal repatternings through a progressive reduction of the chromosome number. A previously suggested new species has been also confirmed on a genetic basis. Data are discussed in order to depict a phylogenetic and phylogeographic scenario fitting the observed genetic relationships between the different species of the group. Chromosome rearrangements are proposed as the major speciation driving force within the group and androgenetic reproduction is suggested as a shortcut to overcome the problem of fixing chromosomal rearrangements that are strongly underdominant in heterozygotes.

Polymerase chain reaction amplification of the Bag320 satellite family reveals the ancestral library and past gene conversion events in Bacillus rossius (Insecta Phasmatodea)

Polymerase chain reaction amplifications of genomic DNA in 17 individuals of bisexual and parthenogenetic populations of three subspecies of Bacillus rossius (Insecta Phasmatodea) revealed that the species still harbours the whole variability of the ancestral Bag320 satellite family, since monomers of all non-hybrid Bacillus taxa plus private sequences occur in it. Bag320 monomers had not been rescued as a major satellite component in B. rossius, but possibly represent the remnant of a set of diverging sequences present in the Bacillus ancestor. Following the library hypothesis, these monomer variants have been differently amplified along the evolutionary pathways leading to present taxa in agreement with the mitochondrial phylogeny of the genus. The putative converted tracts observed are explained as the results of past gene conversion events.

Phylogenetic evidence for hybrid origins of asexual lineages in an aphid species

Understanding the mode of origin of asexuality is central to ongoing debates concerning the evolution and maintenance of sexual reproduction in eukaryotes. This is because it has profound consequences for patterns of genetic diversity and ecological adaptability of asexual lineages, hence on the outcome of competition with sexual relatives both in short and longer terms. Among the possible routes to asexuality, hybridization is a very common mechanism in animals and plants. Aphids present frequent transitions from their ancestral reproductive mode (cyclical parthenogenesis) to permanent asexuality, but the mode of origin of asexual lineages is generally not known because it has never been thoroughly investigated with appropriate molecular tools. Rhopalosiphum padi is an aphid species with coexisting sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages that are genetically distinct. Previous studies have shown that asexual lineages of R. padi are heterozygous at most nuclear loci, suggesting either that they have undergone long-term asexuality (under which heterozygosity tends to increase) or that they have hybrid origins. To discriminate between these alternatives, we conducted an extensive molecular survey combining the sequence analysis of alleles of two nuclear DNA markers and mitochondrial DNA haplotypes in sexual and asexual lineages of R. padi. Both nuclear and cytoplasmic markers clearly showed that many asexual lineages have hybrid origins, the first such demonstration in aphids. Our results also indicated that asexuals result from multiple events of hybridization between R. padi and an unknown sibling species, and are of recent origin (contradicting previous estimates that asexual R. padi lineages were of moderate longevity). This study constitutes another example that putatively ancient asexual lineages are actually of much more recent origin than previously thought. It also presents a robust approach for testing whether hybrid origin of asexuality is also a common phenomenon in aphids.

Recent and ancient asexuality in Timema walkingsticks

Determining the evolutionary age of asexual lineages should help in inferring the temporal scale under which asexuality and sex evolve and assessing selective factors involved in the evolution of asexuality. We used 416 bp of the mitochondrial COI gene to infer phylogenetic relationships of virtually all known Timema walkingstick species, including extensive intraspecific sampling for all five of the asexuals and their close sexual relatives. The asexuals T. douglasi and T. shepardii were very closely related to each other and evolutionarily young (less than 0.5 million years old). For the asexuals T. monikensis and T. tahoe, evidence for antiquity was weak since only one population of each was sampled, intraspecific divergences were low, and genetic distances to related sexuals were high: maximum-likelihood molecular-clock age estimates ranged from 0.26 to 2.39 million years in T. monikensis and from 0.29-1.06 million years in T. tahoe. By contrast, T. genevieve was inferred to be an ancient asexual, with an age of 0.81 to 1.42 million years. The main correlate of the age of asexual lineages was their geographic position, with younger asexuals being found further north.

The evolution of alternative genetic systems in insects

There are three major classes of insect genetic systems: those with diploid males (diplodiploidy), those with effectively haploid males (haplodiploidy), and those without males (thelytoky). Mixed systems, involving cyclic or facultative switching between thelytoky and either of the other systems, also occur. I present a classification of the genetic systems of insects and estimate the number of evolutionary transitions between them that have occurred. Obligate thelytoky has arisen from each of the other systems, and there is evidence that over 900 such origins have occurred. The number of origins of facultative thelytoky and the number of reversions from obligate thelytoky to facultative and cyclic thelytoky are difficult to estimate. The other transitions are few in number: five origins of cyclic thelytoky, eight origins of obligate haplodiploidy (including paternal genome elimination), the strange case of Micromalthus, and the two reversions from haplodiploidy to diplodiploidy in scale insects. Available evidence tends to support W.D. Hamilton's hypothesis that maternally transmitted endosymbionts have been involved in the origins of haplodiploidy. Bizarre systems of extrazygotic inheritance in Sternorrhyncha are not easily accommodated into any existing classification of genetic systems.