Phylogenomics resolves major relationships and reveals significant diversification rate shifts in the evolution of silk moths and relatives.

BACKGROUND: Silkmoths and their relatives constitute the ecologically and taxonomically diverse superfamily Bombycoidea, which includes some of the most charismatic species of Lepidoptera. Despite displaying spectacular forms and diverse ecological traits, relatively little attention has been given to understanding their evolution and drivers of their diversity. To begin to address this problem, we created a new Bombycoidea-specific Anchored Hybrid Enrichment (AHE) probe set and sampled up to 571 loci for 117 taxa across all major lineages of the Bombycoidea, with a newly developed DNA extraction protocol that allows Lepidoptera specimens to be readily sequenced from pinned natural history collections. RESULTS: The well-supported tree was overall consistent with prior morphological and molecular studies, although some taxa were misplaced. The bombycid Arotros Schaus was formally transferred to Apatelodidae. We identified important evolutionary patterns (e.g., morphology, biogeography, and differences in speciation and extinction), and our analysis of diversification rates highlights the stark increases that exist within the Sphingidae (hawkmoths) and Saturniidae (wild silkmoths). CONCLUSIONS: Our study establishes a backbone for future evolutionary, comparative, and taxonomic studies of Bombycoidea. We postulate that the rate shifts identified are due to the well-documented bat-moth "arms race". Our research highlights the flexibility of AHE to generate genomic data from a wide range of museum specimens, both age and preservation method, and will allow researchers to tap into the wealth of biological data residing in natural history collections around the globe.

Independent evolution of sexual dimorphism and female-limited mimicry in swallowtail butterflies (Papilio dardanus and Papilio phorcas).

Several species of swallowtail butterflies (genus Papilio) are Batesian mimics that express multiple mimetic female forms, while the males are monomorphic and nonmimetic. The evolution of such sex-limited mimicry may involve sexual dimorphism arising first and mimicry subsequently. Such a stepwise scenario through a nonmimetic, sexually dimorphic stage has been proposed for two closely related sexually dimorphic species: Papilio phorcas, a nonmimetic species with two female forms, and Papilio dardanus, a female-limited polymorphic mimetic species. Their close relationship indicates that female-limited polymorphism could be a shared derived character of the two species. Here, we present a phylogenomic analysis of the dardanus group using 3964 nuclear loci and whole mitochondrial genomes, showing that they are not sister species and thus that the sexually dimorphic state has arisen independently in the two species. Nonhomology of the female polymorphism in both species is supported by population genetic analysis of engrailed, the presumed mimicry switch locus in P. dardanus. McDonald-Kreitman tests performed on SNPs in engrailed showed the signature of balancing selection in a polymorphic population of P. dardanus, but not in monomorphic populations, nor in the nonmimetic P. phorcas. Hence, the wing polymorphism does not balance polymorphisms in engrailed in P. phorcas. Equally, unlike in P. dardanus, none of the SNPs in P. phorcas engrailed were associated with either female morph. We conclude that sexual dimorphism due to female polymorphism evolved independently in both species from monomorphic, nonmimetic states. While sexual selection may drive male-female dimorphism in nonmimetic species, in mimetic Papilios, natural selection for protection from predators in females is an alternative route to sexual dimorphism.

Mitogenomics of 'Old World Acraea' butterflies reveals a highly divergent 'Bematistes'.

Afrotropical Acraeini butterflies provide a fascinating potential model system to contrast with the Neotropical Heliconiini, yet their phylogeny remains largely unexplored by molecular methods and their generic level nomenclature is still contentious. To test the potential of mitogenomes in a simultaneous analysis of the radiation, we sequenced the full mitochondrial genomes of 19 African species. Analyses show the potential of mitogenomic phylogeny reconstruction in this group. Inferred relationships are largely congruent with a previous multilocus study. We confirm a monophyletic Telchinia to include the Asiatic Pareba with a complicated paraphylum, traditional (sub)genus Acraea, toward the base. The results suggest that several proposed subgenera and some species groups within Telchinia are not monophyletic, while two other (sub)genera could possibly be combined. Telchinia was recovered without strong support as sister to the potentially interesting system of distasteful model butterflies known as Bematistes, a name that is suppressed in some treatments. Surprisingly, we find that this taxon has remarkably divergent mitogenomes and unexpected synapomorphic tRNA rearrangements. These gene order changes, combined with evidence for deviating dN/dS ratios and evidence for episodal diversifying selection, suggest that the ancestral Bematistes mitogenome has had a turbulent past. Our study adds genetic support for treating this clade as a distinct genus, while the alternative option, adopted by some authors, of Acraea being equivalent to Acraeini merely promotes redundancy. We pave the way for more detailed mitogenomic and multi-locus molecular analyses which can determine how many genera are needed (possibly at least six) to divide Acraeini into monophyletic groups that also facilitate communication about their biology.

Interplay of robustness and plasticity of life-history traits drives ecotypic differentiation in thermally distinct habitats.

Phenotypic plasticity describes the ability of an individual to alter its phenotype in response to the environment and is potentially adaptive when dealing with environmental variation. However, robustness in the face of a changing environment may often be beneficial for traits that are tightly linked to fitness. We hypothesized that robustness of some traits may depend on specific patterns of plasticity within and among other traits. We used a reaction norm approach to study robustness and phenotypic plasticity of three life-history traits of the collembolan Orchesella cincta in environments with different thermal regimes. We measured adult mass, age at maturity and growth rate of males and females from heath and forest habitats at two temperatures (12 and 22 °C). We found evidence for ecotype-specific robustness of female adult mass to temperature, with a higher level of robustness in the heath ecotype. This robustness is facilitated by plastic adjustments of growth rate and age at maturity. Furthermore, female fecundity is strongly influenced by female adult mass, explaining the importance of realizing a high mass across temperatures for females. These findings indicate that different predicted outcomes of life-history theory can be combined within one species' ontogeny and that models describing life-history strategies should not assume that traits like growth rate are maximized under all conditions. On a methodological note, we report a systematic inflation of variation when standard deviations and correlation coefficients are calculated from family means as opposed to individual data within a family structure.

Why barcode? High-throughput multiplex sequencing of mitochondrial genomes for molecular systematics.

Mitochondrial genome sequences are important markers for phylogenetics but taxon sampling remains sporadic because of the great effort and cost required to acquire full-length sequences. Here, we demonstrate a simple, cost-effective way to sequence the full complement of protein coding mitochondrial genes from pooled samples using the 454/Roche platform. Multiplexing was achieved without the need for expensive indexing tags ('barcodes'). The method was trialled with a set of long-range polymerase chain reaction (PCR) fragments from 30 species of Coleoptera (beetles) sequenced in a 1/16th sector of a sequencing plate. Long contigs were produced from the pooled sequences with sequencing depths ranging from ∼10 to 100× per contig. Species identity of individual contigs was established via three 'bait' sequences matching disparate parts of the mitochondrial genome obtained by conventional PCR and Sanger sequencing. This proved that assembly of contigs from the sequencing pool was correct. Our study produced sequences for 21 nearly complete and seven partial sets of protein coding mitochondrial genes. Combined with existing sequences for 25 taxa, an improved estimate of basal relationships in Coleoptera was obtained. The procedure could be employed routinely for mitochondrial genome sequencing at the species level, to provide improved species 'barcodes' that currently use the cox1 gene only.

Sugar sweet springtails: on the transcriptional response of Folsomia candida (Collembola) to desiccation stress.

Several species of Collembola survive stressful desiccating conditions by absorbing water vapour from the environment. To obtain insight into the transcriptomic responses underlying this 'water vapour absorption' mechanism we subjected Folsomia candida (Collembola) to transcriptome profiling. We show that ecologically relevant desiccation stress leads to strong time-dependent transcriptomic changes. Exposure of F. candida to 98.2% relative humidity over an interval of 174 h resulted in a high number of gene transcripts being differentially expressed (up to 41%; P-value < 0.05). Additional Gene Ontology analyses suggest that carbohydrate transport, sugar catabolism and cuticle maintenance are biological processes involved in combating desiccation. However, many additional pathways seem to be affected; additional experiments are needed to elucidate which responses are primarily linked to desiccation resistance.

Revealing pancrustacean relationships: phylogenetic analysis of ribosomal protein genes places Collembola (springtails) in a monophyletic Hexapoda and reinforces the discrepancy between mitochondrial and nuclear DNA markers.

BACKGROUND: In recent years, several new hypotheses on phylogenetic relations among arthropods have been proposed on the basis of DNA sequences. One of the challenged hypotheses is the monophyly of hexapods. This discussion originated from analyses based on mitochondrial DNA datasets that, due to an unusual positioning of Collembola, suggested that the hexapod body plan evolved at least twice. Here, we re-evaluate the position of Collembola using ribosomal protein gene sequences. RESULTS: In total 48 ribosomal proteins were obtained for the collembolan Folsomia candida. These 48 sequences were aligned with sequence data on 35 other ecdysozoans. Each ribosomal protein gene was available for 25% to 86% of the taxa. However, the total sequence information was unequally distributed over the taxa and ranged between 4% and 100%. A concatenated dataset was constructed (5034 inferred amino acids in length), of which ~66% of the positions were filled. Phylogenetic tree reconstructions, using Maximum Likelihood, Maximum Parsimony, and Bayesian methods, resulted in a topology that supports monophyly of Hexapoda. CONCLUSION: Although ribosomal proteins in general may not evolve independently, they once more appear highly valuable for phylogenetic reconstruction. Our analyses clearly suggest that Hexapoda is monophyletic. This underpins the inconsistency between nuclear and mitochondrial datasets when analyzing pancrustacean relationships. Caution is needed when applying mitochondrial markers in deep phylogeny.

Allelic diversity of metallothionein in Orchesella cincta (L.): traces of natural selection by environmental pollution.

The advances made in statistical methods to detect selection from DNA sequence variation has resulted in an enormous increase in the number of studies reporting positive selection. However, a disadvantage of such statistical tests is that often no insight into the actual source of selection is obtained. Finer understanding of evolution can be obtained when those statistical tests are combined with field observations on allele frequencies. We assessed whether the metallothionein (mt) gene of Orchesella cincta (Collembola), which codes for a metal-binding protein, is subject to selection, by investigating alleles and allele frequencies among European metal-stressed and reference populations. Eight highly divergent alleles were resolved in Northwest Europe. At the nucleotide level, a total of 51 polymorphic sites (five of them implying amino-acid changes) were observed. Although statistical tests applied to the sequences alone showed no indication of selection, a G-test rejected the null hypothesis that alleles are homogeneously distributed over metal-stressed and reference populations. Analysis of molecular variance assigned a small, but significant amount of the total variance to differences between metal-stressed and non-stressed populations. In addition, it was shown that metal-stressed populations tend to be more genetically diversified at this locus than non-stressed ones. These results suggest that the mt gene and its surrounding DNA region are affected by environmental metal contamination. This study illustrates that, in addition to statistical tests, field observations on allele frequencies are needed to gain understanding of selection and adaptive evolution.

Genetic structure in Orchesella cincta (Collembola): strong subdivision of European populations inferred from mtDNA and AFLP markers.

Population genetic structure is determined both by current processes and historical events. Current processes include gene flow, which is largely influenced by the migration capacity of a species. Historical events are, for example, glaciation periods, which have had a major impact on the distribution of many species. Species with a low capacity or tendency to move about or disperse often exhibit clear spatial genetic structures, whereas mobile species mostly show less spatial genetic differentiation. In this paper we report on the genetic structure of a small, wingless arthropod species (Orchesella cincta: Collembola) in Europe. For this purpose we used mtDNA COII sequences and AFLP markers. We show that large genetic differences exist between populations of O. cincta, as expected from O. cincta's winglessness and sedentary lifestyle. Despite the fact that most variability was observed within populations (59%), a highly significant amount of AFLP variation (25%) was observed between populations from northwestern Europe, central Europe and Italy. This suggests that gene flow among regions is extremely low, which is additionally supported by the lack of shared mtDNA alleles between regions. Based on the genetic variation and sequence differences observed we conclude that the subdivision occurred long before the last glaciation periods. Although the populations still interbreed in the lab, we assume that in the long term the genetic isolation of these regions may lead to speciation processes.

Infectious parthenogenesis.

Parthenogenesis-inducing Wolbachia bacteria are reproductive parasites that cause infected female wasps to produce daughters without mating. This manipulation of the host's reproduction enhances the transmission of Wolbachia to future generations because the bacteria are passed on vertically only from mothers to daughters. Males are dead ends for cytoplasmically inherited bacteria: they do not pass them on to their offspring. Vertical transmission of Wolbachia has been previously considered to be the main mode of transmission. Here we report frequent horizontal transmission from infected to uninfected wasp larvae sharing a common food source. The transferred Wolbachia are then vertically transmitted to the new host's offspring. This natural and unexpectedly frequent horizontal transfer of parthenogensis-inducing Wolbachia intraspecifically has important implications for the co-evolution of Wolbachia and their host.