A large-scale, higher-level, molecular phylogenetic study of the insect order Lepidoptera (moths and butterflies)

Mitogenomic sequences effectively recover relationships within brush-footed butterflies (Lepidoptera: Nymphalidae)

BACKGROUND

Mitogenomic phylogenies have revealed well-supported relationships for many eukaryote groups. In the order Lepidoptera, 113 species mitogenomes had been sequenced (May 14, 2014). However, these data are restricted to ten of the forty-three recognised superfamilies, while it has been challenging to recover large numbers of mitogenomes due to the time and cost required for primer design and sequencing. Nuclear rather than mitochondrial genes have been preferred to reconstruct deep-level lepidopteran phylogenies, without seriously evaluating the potential of entire mitogenomes. Next-generation sequencing methods remove these limitations by providing efficiently massive amounts of sequence data. In the present study, we simultaneously obtained a large number of nymphalid butterfly mitogenomes to evaluate the utility of mitogenomic phylogenies by comparing reconstructions to the now quite well established phylogeny of Nymphalidae.

RESULTS

We newly obtained 30 nymphalid mitogenomes via pyrosequencing on the Roche 454 GS Junior system, and combined these sequences with publicly accessible data to provide a 70-taxa dataset covering 37 genes for a 15,495 bp alignment. Polymorphic sites were not homogeneously distributed across the gene. Two gene regions, nad6 and 3' end of nad5, were most variable, whereas the cox1 and 5' ends of rrnL were most conserved. Phylogenetic relationships inferred by two likelihood methods were congruent and strongly supported (>0.95 posterior probability; ML bootstrap >85%), across the majority of nodes for multiple partitioning strategies and substitution models. Bayes factor results showed that the most highly partitioned dataset is the preferred strategy among different partitioning schemes. The most striking phylogenetic findings were that the subfamily Danainae not Libytheinae was sister of the remaining brush-footed butterflies and that, within Limenitidini, the genus Athyma was clearly polyphyletic. None of the single-gene phylogenies recovered the highly supported topologies generated on the basis of the whole mitogenomic data.

CONCLUSIONS

Thirty mitogenomes were assembled with 89% completeness from the contigs of pyrosequencing-derived reads. Entire mitogenomes or higher-quality sequences could be obtained by increasing pyrosequencing read coverage or by additional Sanger sequencing. Our mitogenomic phylogenies provide robust nodal support at a range of levels, demonstrating that mitogenomes are both accurate and efficient molecular markers for inferring butterfly phylogeny.

Towards a mitogenomic phylogeny of Lepidoptera

The backbone phylogeny of Lepidoptera remains unresolved, despite strenuous recent morphological and molecular efforts. Molecular studies have focused on nuclear protein coding genes, sometimes adding a single mitochondrial gene. Recent advances in sequencing technology have, however, made acquisition of entire mitochondrial genomes both practical and economically viable. Prior phylogenetic studies utilised just eight of 43 currently recognised lepidopteran superfamilies. Here, we add 23 full and six partial mitochondrial genomes (comprising 22 superfamilies of which 16 are newly represented) to those publically available for a total of 24 superfamilies and ask whether such a sample can resolve deeper lepidopteran phylogeny. Using recoded datasets we obtain topologies that are highly congruent with prior nuclear and/or morphological studies. Our study shows support for an expanded Obtectomera including Gelechioidea, Thyridoidea, plume moths (Alucitoidea and Pterophoroidea; possibly along with Epermenioidea), Papilionoidea, Pyraloidea, Mimallonoidea and Macroheterocera. Regarding other controversially positioned higher taxa, Doidae is supported within the new concept of Drepanoidea and Mimallonidae sister to (or part of) Macroheterocera, while among Nymphalidae butterflies, Danainae and not Libytheinae are sister to the remainder of the family. At the deepest level, we suggest that a tRNA rearrangement occurred at a node between Adeloidea and Ditrysia+Palaephatidae+Tischeriidae.

A gateway for phylogenetic analysis powered by grid computing featuring GARLI 2.0

We introduce molecularevolution.org, a publicly available gateway for high-throughput, maximum-likelihood phylogenetic analysis powered by grid computing. The gateway features a garli 2.0 web service that enables a user to quickly and easily submit thousands of maximum likelihood tree searches or bootstrap searches that are executed in parallel on distributed computing resources. The garli web service allows one to easily specify partitioned substitution models using a graphical interface, and it performs sophisticated post-processing of phylogenetic results. Although the garli web service has been used by the research community for over three years, here we formally announce the availability of the service, describe its capabilities, highlight new features and recent improvements, and provide details about how the grid system efficiently delivers high-quality phylogenetic results. [garli, gateway, grid computing, maximum likelihood, molecular evolution portal, phylogenetics, web service.]

Does more sequence data improve estimates of galliform phylogeny? Analyses of a rapid radiation using a complete data matrix

The resolution of rapid evolutionary radiations or "bushes" in the tree of life has been one of the most difficult and interesting problems in phylogenetics. The avian order Galliformes appears to have undergone several rapid radiations that have limited the resolution of prior studies and obscured the position of taxa important both agriculturally and as model systems (chicken, turkey, Japanese quail). Here we present analyses of a multi-locus data matrix comprising over 15,000 sites, primarily from nuclear introns but also including three mitochondrial regions, from 46 galliform taxa with all gene regions sampled for all taxa. The increased sampling of unlinked nuclear genes provided strong bootstrap support for all but a small number of relationships. Coalescent-based methods to combine individual gene trees and analyses of datasets that are independent of published data indicated that this well-supported topology is likely to reflect the galliform species tree. The inclusion or exclusion of mitochondrial data had a limited impact upon analyses upon analyses using either concatenated data or multispecies coalescent methods. Some of the key phylogenetic findings include support for a second major clade within the core phasianids that includes the chicken and Japanese quail and clarification of the phylogenetic relationships of turkey. Jackknifed datasets suggested that there is an advantage to sampling many independent regions across the genome rather than obtaining long sequences for a small number of loci, possibly reflecting the differences among gene trees that differ due to incomplete lineage sorting. Despite the novel insights we obtained using this increased sampling of gene regions, some nodes remain unresolved, likely due to periods of rapid diversification. Resolving these remaining groups will likely require sequencing a very large number of gene regions, but our analyses now appear to support a robust backbone for this order.

Differences in nicotine metabolism of two Nicotiana attenuata herbivores render them differentially susceptible to a common native predator

BACKGROUND

Nicotiana attenuata is attacked by larvae of both specialist (Manduca sexta) and generalist (Spodoptera exigua) lepidopteran herbivores in its native habitat. Nicotine is one of N. attenuata's important defenses. M. sexta is highly nicotine tolerant; whether cytochrome P450 (CYP)-mediated oxidative detoxification and/or rapid excretion is responsible for its exceptional tolerance remains unknown despite five decades of study. Recently, we demonstrated that M. sexta uses its nicotine-induced CYP6B46 to efflux midgut-nicotine into the hemolymph, facilitating nicotine exhalation that deters predatory wolf spiders (Camptocosa parallela). S. exigua's nicotine metabolism is uninvestigated.

METHODOLOGY/PRINCIPAL FINDINGS

We compared the ability of these two herbivores to metabolize, tolerate and co-opt ingested nicotine for defense against the wolf spider. In addition, we analyzed the spider's excretion to gain insights into its nicotine metabolism. Contrary to previous reports, we found that M. sexta larvae neither accumulate the common nicotine oxides (cotinine, cotinine N-oxide and nicotine N-oxide) nor excrete them faster than nicotine. In M. sexta larvae, ingestion of nicotine as well as its oxides increases the accumulation of CYP6B46 transcripts. In contrast, S. exigua accumulates nicotine oxides and exhales less (66%) nicotine than does M. sexta. Spiders prefer nicotine-fed S. exigua over M. sexta, a preference reversed by topical or headspace nicotine supplementation, but not ingested or topically-coated nicotine oxides, suggesting that externalized nicotine but not the nicotine detoxification products deter spider predation. The spiders also do not accumulate nicotine oxides.

CONCLUSIONS

Nicotine oxidation reduces S. exigua's headspace-nicotine and renders it more susceptible to predation by spiders than M. sexta, which exhales unmetabolized nicotine. These results are consistent with the hypothesis that generalist herbivores incur costs of detoxification, which include the ecological costs of greater predation risks, in addition to the previously demonstrated energetic, physiological and metabolic costs.

Morphology-based phylogeny and biogeography of Wockia (Lepidoptera: Urodidae) with description of a new species from Japan and South Korea

Wockia, one of the six genera within the lepidopteran family Urodidae, currently includes nine species distributed in Holarctic, Oriental, and Neotropical regions. A new species of the genus, W. magna sp. nov., is described from Japan and South Korea. This is the first record of Urodidae from Japan. A cladistic analysis was constructed based on 27 morphological characters from five ingroup species and two outgroup taxa. A single most parsimonious tree was found (length = 38, CI = 71, RI = 70). The resulting tree failed to recover the monophyly of Wockia chewbacca and a clade Wockia sensu stricto, including four congeners from North America, temperate East Asia and Southeast Asia. Wockia sensu stricto was resolved as monophyletic and was divided into two subgroups, one including W. asperipunctella and W. magna, and the other including W. koreana and W. variata. The sister-group relationship of W. asperipunctella and W. magna was moderately supported. Synapomorphies of Wockia sensu stricto are provided from the resulting phylogeny. Systematic definitions of Wockia and other allied genera are revised. Optimal ancestral area reconstruction implemented in DIVA resulted in different hypotheses for Wockia, depending on how to constrain the maximum number of areas. The maximum area number set as two yielded a more likely scenario suggesting that Wockia sensu stricto originated in temperate East Asia and then dispersed into the Oriental region, North America, and Europe. A working hypothesis and other equally possible alternative explanations for the biogeography of Wockia are provided.

Statistical tests to identify appropriate types of nucleotide sequence recoding in molecular phylogenetics

Background

Under a Markov model of evolution, recoding, or lumping, of the four nucleotides into fewer groups may permit analysis under simpler conditions but may unfortunately yield misleading results unless the evolutionary process of the recoded groups remains Markovian. If a Markov process is lumpable, then the evolutionary process of the recoded groups is Markovian.

Results

We consider stationary, reversible, and homogeneous Markov processes on two taxa and compare three tests for lumpability: one using an ad hoc test statistic, which is based on an index that is evaluated using a bootstrap approximation of its distribution; one that is based on a test proposed specifically for Markov chains; and one using a likelihood-ratio test. We show that the likelihood-ratio test is more powerful than the index test, which is more powerful than that based on the Markov chain test statistic. We also show that for stationary processes on binary trees with more than two taxa, the tests can be applied to all pairs. Finally, we show that if the process is lumpable, then estimates obtained under the recoded model agree with estimates obtained under the original model, whereas, if the process is not lumpable, then these estimates can differ substantially. We apply the new likelihood-ratio test for lumpability to two primate data sets, one with a mitochondrial origin and one with a nuclear origin.

Conclusions

Recoding may result in biased phylogenetic estimates because the original evolutionary process is not lumpable. Accordingly, testing for lumpability should be done prior to phylogenetic analysis of recoded data.

CLIMBER: Climatic niche characteristics of the butterflies in Europe

Detailed information on species’ ecological niche characteristics that can be related to declines and extinctions is indispensable for a better understanding of the relationship between the occurrence and performance of wild species and their environment and, moreover, for an improved assessment of the impacts of global change. Knowledge on species characteristics such as habitat requirements is already available in the ecological literature for butterflies, but information about their climatic requirements is still lacking. Here we present a unique dataset on the climatic niche characteristics of 397 European butterflies representing 91% of the European species (see Appendix). These characteristics were obtained by combining detailed information on butterfly distributions in Europe (which also led to the ‘Distribution Atlas of Butterflies in Europe’) and the corresponding climatic conditions. The presented dataset comprises information for the position and breadth of the following climatic niche characteristics: mean annual temperature, range in annual temperature, growing degree days, annual precipitation sum, range in annual precipitation and soil water content. The climatic niche position is indicated by the median and mean value for each climate variable across a species’ range, accompanied by the 95% confidence interval for the mean and the number of grid cells used for calculations. Climatic niche breadth is indicated by the standard deviation and the minimum and maximum values for each climatic variable across a species’ range. Database compilation was based on high quality standards and the data are ready to use for a broad range of applications.

It is already evident that the information provided in this dataset is of great relevance for basic and applied ecology. Based on the species temperature index (STI, i.e. the mean temperature value per species), the community temperature index (CTI, i.e. the average STI value across the species in a community) was recently adopted as an indicator of climate change impact on biodiversity by the pan-European framework supporting the Convention on Biological Diversity (Streamlining European Biodiversity Indicators 2010) and has already been used in several scientific publications. The application potential of this database ranges from theoretical aspects such as assessments of past niche evolution or analyses of trait interdependencies to the very applied aspects of measuring, monitoring and projecting historical, ongoing and potential future responses to climate change using butterflies as an indicator.

Can RNA-Seq resolve the rapid radiation of advanced moths and butterflies (Hexapoda: Lepidoptera: Apoditrysia)? An exploratory study

Recent molecular phylogenetic studies of the insect order Lepidoptera have robustly resolved family-level divergences within most superfamilies, and most divergences among the relatively species-poor early-arising superfamilies. In sharp contrast, relationships among the superfamilies of more advanced moths and butterflies that comprise the mega-diverse clade Apoditrysia (ca. 145,000 spp.) remain mostly poorly supported. This uncertainty, in turn, limits our ability to discern the origins, ages and evolutionary consequences of traits hypothesized to promote the spectacular diversification of Apoditrysia. Low support along the apoditrysian "backbone" probably reflects rapid diversification. If so, it may be feasible to strengthen resolution by radically increasing the gene sample, but case studies have been few. We explored the potential of next-generation sequencing to conclusively resolve apoditrysian relationships. We used transcriptome RNA-Seq to generate 1579 putatively orthologous gene sequences across a broad sample of 40 apoditrysians plus four outgroups, to which we added two taxa from previously published data. Phylogenetic analysis of a 46-taxon, 741-gene matrix, resulting from a strict filter that eliminated ortholog groups containing any apparent paralogs, yielded dramatic overall increase in bootstrap support for deeper nodes within Apoditrysia as compared to results from previous and concurrent 19-gene analyses. High support was restricted mainly to the huge subclade Obtectomera broadly defined, in which 11 of 12 nodes subtending multiple superfamilies had bootstrap support of 100%. The strongly supported nodes showed little conflict with groupings from previous studies, and were little affected by changes in taxon sampling, suggesting that they reflect true signal rather than artifacts of massive gene sampling. In contrast, strong support was seen at only 2 of 11 deeper nodes among the "lower", non-obtectomeran apoditrysians. These represent a much harder phylogenetic problem, for which one path to resolution might include further increase in gene sampling, together with improved orthology assignments.

Timing and patterns in the taxonomic diversification of Lepidoptera (butterflies and moths)

The macroevolutionary history of the megadiverse insect order Lepidoptera remains little-known, yet coevolutionary dynamics with their angiospermous host plants are thought to have influenced their diversification significantly. We estimate the divergence times of all higher-level lineages of Lepidoptera, including most extant families. We find that the diversification of major lineages in Lepidoptera are approximately equal in age to the crown group of angiosperms and that there appear to have been three significant increases in diversification rates among Lepidoptera over evolutionary time: 1) at the origin of the crown group of Ditrysia about 150 million years ago (mya), 2) at the origin of the stem group of Apoditrysia about 120 mya and finally 3) a spectacular increase at the origin of the stem group of the quadrifid noctuoids about 70 mya. In addition, there appears to be a significant increase in diversification rate in multiple lineages around 90 mya, which is concordant with the radiation of angiosperms. Almost all extant families appear to have begun diversifying soon after the Cretaceous/Paleogene event 65.51 mya.

Notable plesiomorphies and notable specializations: head structure of the primitive "tongue moth" Acanthopteroctetes unifascia (Lepidoptera: Acanthopteroctetidae)

The Acanthopteroctetidae are one of the first-originated family-group lineages within "tongue moths" (Lepidoptera-Glossata). The purpose of this study is to provide a comprehensive account (based on whole mount preparations, serial sections, and Scanning electron microscopy) of the cephalic structure of an adult exemplar of the family, to supplement the sparse available information. Notable plesiomorphies include the retention of frontal retractors of the narrow labrum, a high supraocular index linked to strong development of cranio-mandibular ad- and abductors, and perhaps the unusually short but still coilable (just ca. 1.5 turns) galeal "tongue." Notable specializations (probably mostly family autapomorphies) include a complement of large sensilla placodea on the male antennae, an apical attachment of the long dorsal tentorial arm to the cranium, an extreme reduction of the single-segmented labial palps, a particularly strong subgenal bridge and a surface structure of near-parallel ridges on the ommatidial corneae. The presence of sizable saccular mandibular (type 1) glands opening into the adductor apodeme is unexpected, no counterparts being known from neighboring taxa. The same is true for ventral salivarium dilator muscles originating on the prelabium; and tentatively suggested to be homologues of the extrinsic palp flexors (the insertion shift being related to loss of original function due to palp reduction), rather than to the ventral salivarium muscles of more basal insects. A complete "deutocerebral loop"' may or may not be developed, as is true for a mutual appression of the optic lobe and circumoesophageal connective/suboesophageal ganglion, enclosing the anterior tentorial arm between them; a suboesophageal innervation of the retrocerebral complex was not observed. No characters bearing on the monophyly of the Coelolepida were identified. The scapo-pedicellar articulation with a scapal process and a smooth intercalary sclerite is reminiscent of conditions in Neopseustidae, but remains debatable as a synapomorphy of the two families.

Hawkmoths produce anti-bat ultrasound

Bats and moths have been engaged in aerial warfare for nearly 65 Myr. This arms race has produced a suite of counter-adaptations in moths, including bat-detecting ears. One set of defensive strategies involves the active production of sound; tiger moths' ultrasonic replies to bat attack have been shown to startle bats, warn the predators of bad taste and jam their biosonar. Here, we report that hawkmoths in the Choerocampina produce entirely ultrasonic sounds in response to tactile stimulation and the playback of biosonar attack sequences. Males do so by grating modified scraper scales on the outer surface of the genital valves against the inner margin of the last abdominal tergum. Preliminary data indicate that females also produce ultrasound to touch and playback of echolocation attack, but they do so with an entirely different mechanism. The anti-bat function of these sounds is unknown but might include startling, cross-family acoustic mimicry, warning of unprofitability or physical defence and/or jamming of echolocation. Hawkmoths present a novel and tractable system to study both the function and evolution of anti-bat defences.

Predicting novel trophic interactions in a non-native world

Humans are altering the global distributional ranges of plants, while their co-evolved herbivores are frequently left behind. Native herbivores often colonise non-native plants, potentially reducing invasion success or causing economic loss to introduced agricultural crops. We developed a predictive model to forecast novel interactions and verified it with a data set containing hundreds of observed novel plant-insect interactions. Using a food network of 900 native European butterfly and moth species and 1944 native plants, we built an herbivore host-use model. By extrapolating host use from the native herbivore-plant food network, we accurately forecasted the observed novel use of 459 non-native plant species by native herbivores. Patterns that governed herbivore host breadth on co-evolved native plants were equally important in determining non-native hosts. Our results make the forecasting of novel herbivore communities feasible in order to better understand the fate and impact of introduced plants.

Chromosomal evolution in tortricid moths: conserved karyotypes with diverged features

Moths of the family Tortricidae constitute one of the major microlepidopteran groups in terms of species richness and economic importance. Yet, despite their overall significance, our knowledge of their genome organization is very limited. In order to understand karyotype evolution in the family Tortricidae, we performed detailed cytogenetic analysis of Grapholita molesta, G. funebrana, Lobesia botrana, and Eupoecilia ambiguella, representatives of two main tortricid subfamilies, Olethreutinae and Tortricinae. Besides standard cytogenetic methods, we used fluorescence in situ hybridization for mapping of major rRNA and histone gene clusters and comparative genomic hybridization to determine the level of molecular differentiation of the W and Z sex chromosomes. Our results in combination with available data in the codling moth, Cydia pomonella, and other tortricids allow us a comprehensive reconstruction of chromosomal evolution across the family Tortricidae. The emerging picture is that the karyotype of a common ancestor of Tortricinae and Olethreutinae differentiated from the ancestral lepidopteran chromosome print of n = 31 by a sex chromosome-autosome fusion. This rearrangement resulted in a large neo-sex chromosome pair and a karyotype with n = 30 conserved in most Tortricinae species, which was further reduced to n = 28 observed in Olethreutinae. Comparison of the tortricid neo-W chromosomes showed differences in their structure and composition presumably reflecting stochasticity of molecular degeneration of the autosomal part of the neo-W chromosome. Our analysis also revealed conservative pattern of the histone distribution, which is in contrast with high rDNA mobility. Despite the dynamic evolution of rDNA, we can infer a single NOR-chromosome pair as an ancestral state not only in tortricids but probably in all Lepidoptera. The results greatly expand our knowledge of the genome architecture in tortricids, but also contribute to the understanding of chromosomal evolution in Lepidoptera in general.

Evolution of Manduca sexta hornworms and relatives: biogeographical analysis reveals an ancestral diversification in Central America

The hawkmoth genus Manduca is a diverse group of very large, conspicuous moths that has served as an important model across many biological disciplines. Two species in particular, the tobacco hornworm (Manduca sexta) and the tomato hornworm (Manduca quinquemaculatus) have been researched extensively. Studies across biological fields have referred to these two species as being closely related or even sister species, but the extent to which these two model organisms are related remains largely unknown. We conducted a comprehensive multi-gene phylogenetic analysis of Manduca, based on both an ML and Bayesian framework, which resulted in a monophyletic Manduca but only when two other genera, Dolba and Euryglottis are included. We tentatively conclude that the sister group to Manduca sexta comprises the Caribbean M. afflicta and M. johanni, and the sister lineage to this clade includes M. quinquemaculatus and the Hawaiian M. blackburni. Thus, M. sexta and M. quinquemaculatus are closely related, but are not sister species. Biogeographical analyses reveal an ancestral center of diversification in Central America, and Manduca appears to have subsequently colonized North and South America. Our phylogeny provides an important foundation for comparative studies of two model organisms and their relatives.