Ancient rapid radiations of insects: challenges for phylogenetic analysis

Novel phylogenomic inference and 'Out of Asia' biogeography of cobras, coral snakes and their allies

Estimation of evolutionary relationships among lineages that rapidly diversified can be challenging, and, in such instances, inaccurate or unresolved phylogenetic estimates can lead to erroneous conclusions regarding historical geographical ranges of lineages. One example underscoring this issue has been the historical challenge posed by untangling the biogeographic origin of elapoid snakes, which includes numerous dangerously venomous species as well as species not known to be dangerous to humans. The worldwide distribution of this lineage makes it an ideal group for testing hypotheses related to historical faunal exchanges among the many continents and other landmasses occupied by contemporary elapoid species. We developed a novel suite of genomic resources, included worldwide sampling, and inferred a robust estimate of evolutionary relationships, which we leveraged to quantitatively estimate geographical range evolution through the deep-time history of this remarkable radiation. Our phylogenetic and biogeographical estimates of historical ranges definitively reject a lingering former 'Out of Africa' hypothesis and support an 'Out of Asia' scenario involving multiple faunal exchanges between Asia, Africa, Australasia, the Americas and Europe.

A new nuclear phylogeny of the tea family (Theaceae) unravels rapid radiations in genus Camellia

Molecular analyses of rapidly radiating groups often reveal incongruence between gene trees. This mainly results from incomplete lineage sorting, introgression, and gene tree estimation error, which complicate the estimation of phylogenetic relationships. In this study, we reconstruct the phylogeny of Theaceae using 348 nuclear loci from 68 individuals and two outgroup taxa. Sequence data were obtained by target enrichment using the recently released Angiosperm 353 universal probe set applied to herbarium specimens. The robustness of the topologies to variation in data quality was established under a range of different filtering schemes, using both coalescent and concatenation approaches. Our results confirmed most of the previously hypothesized relationships among tribes and genera, while clarifying additional interspecific relationships within the rapidly radiating genus Camellia. We recovered a remarkably high degree of gene tree heterogeneity indicative of rapid radiation in the group and observed cytonuclear conflicts, especially within Camellia. This was especially pronounced around short branches, which we primarily associate with gene tree estimation error. Our analysis also indicates that incomplete lineage sorting (ILS) contributed to gene-tree conflicts and accounted for approximately 14 % of the explained variation, whereas inferred introgression levels were low. Our study advances the understanding of the evolution of this important plant family and provides guidance on the application of target capture methods and the evaluation of key processes that influence phylogenetic discordances.

Comparative Mitogenome of Phylogenetic Relationships and Divergence Time Analysis within Potamanthidae (Insecta: Ephemeroptera)

Potamanthidae belongs to the superfamily Ephemeroidea but has no complete mt genome released in the NCBI (except for two unchecked and one partial mt genome). Since the sister clade to Potamanthidae has always been controversial, we sequenced seven mt genomes of Potamanthidae (two species from Rhoenanthus and five species from Potamanthus) in order to rebuild the phylogenetic relationships of Potamanthidae in this study. The divergence time of Potamanthidae was also investigated by utilizing five fossil calibration points because of the indeterminate origin time. In addition, because Rhoenanthus coreanus and Potamanthus luteus are always in low-temperature environments, we aimed to explore whether these two species were under positive selection at the mt genome level. Amongst the 13 PCGs, CGA was used as the start codon in COX1, whereas other genes conformed to initiating with an ATN start codon. From this analysis, UUA (L), AUU (I), and UUU (F) had the highest usage. Furthermore, the DHU arm was absent in the secondary structure of S1 in all species. By combining the 13 PCGs and 2 rRNAs, we reconstructed the phylogenetic relationship of Potamanthidae within Ephemeroptera. The monophyly of Potamanthidae and the monophyly of Rhoenanthus and Potamanthus were supported in the results. The phylogenetic relationship of Potamanthidae + (Ephemeridae + Polymitarcyidae) was also recovered with a high prior probability. The divergence times of Potamanthidae were traced to be 90.44 Mya (95% HPD, 62.80-121.74 Mya), and the divergence times of Rhoenanthus and Potamanthus originated at approximately 64.77 Mya (95% HPD, 43.82-88.68 Mya), thus belonging to the late Pliocene Epoch or early Miocene Epoch. In addition, the data indicated that R. coreanus was under negative selection and that ATP8 and ND2 in Potamanthidae had a high evolutionary rate.

Phylogeny of North Americas largest cicada radiation redefines Tibicinoides and Okanagana (Hemiptera: Auchenorrhyncha: Cicadidae: Tibicininae)

Tibicinoides, with three small endemic California cicada species, has a confusing, intertwined systematic history with Okanagana that we unravel here. An ingroup including all species of Tibicinoides and the majority (84.7%) of Okanagana species were sampled for six gene regions, polarized with Clidophleps, Okanagodes, Subpsaltria, and Tibicina outgroups, and subjected to Bayesian phylogenetic analysis. Although the ingroup was monophyletic from all outgroups including Tibicina, Tibicinoides rendered Okanagana paraphyletic among two major ingroup clades. To bring classification into agreement with phylogeny, we redescribe and redefine Tibicinoides to include all Okanagana species with a hooked uncus in the male genitalia, all of which grouped with the type T. cupreosparsa (Uhler, 1889) in the first of these clades: T. boweni (Chatfield-Taylor & Cole, 2020) comb. n., T. catalina (Davis, 1936) comb. n., T. hesperia (Uhler, 1876) comb. n., T. mercedita (Davis, 1915), T. minuta (Davis, 1915), T. pallidula (Davis, 1917a) comb. n., T. pernix (Bliven, 1964) comb. n., T. rubrovenosa (Davis, 1915) comb. n., T. simulata (Davis, 1921) comb. n., T. striatipes (Haldeman, 1852) comb. n., T. uncinata (Van Duzee, 1915) comb. n., T. utahensis (Davis, 1919) comb. n., and T. vanduzeei (Distant, 1914) comb. n. Okanagana is redescribed and restricted to the species of the second major clade which contained the type O. rimosa (Say, 1830). We describe two new genera for morphologically distinct orphan lineages: Chlorocanta gen. nov. for C. viridis (Davis, 1918) comb. n. and Hewlettia gen. nov. for H. nigriviridis (Davis, 1921) comb. n. We recognize O. rubrobasalis Davis, 1926 stat. rev. as a species and relegate two former species to junior subjective synonyms: O. noveboracensis (Emmons, 1854) = O. canadensis (Provancher, 1889) and O. occidentalis (Walker in Lord, 1866) = O. lurida Davis, 1919. Tibicinoides and Okanagana together represent a rapid radiation that presents challenges to phylogenetic analysis including suboptimal outgroups and short internodes.

Phylogeny, biogeography, and character evolution of the genus Sophora s.l. (Fabaceae, Papilionoideae)

The papilionoid legume genus Sophora (Fabaceae) exhibits a worldwide distribution, but a phylogenetic framework to understand the evolution of this group is lacking to date. Previous studies have demonstrated that Sophora is not monophyletic and might include Ammodendron, Ammothamnus, and Echinosophora, but the relationships among these four genera (defined as Sophora s.l.) are unclear. Here we used a nuclear DNA dataset (ETS, ITS, SQD1) and a plastid DNA dataset (matK, rbcL, rpl32-trnL, trnL-F) of 654 accession sequences to reconstruct the phylogenetic relationships, estimate the divergence times and ancestral range of Sophora s.l., and infer the evolution of chromosome number and morphological characteristics. Our major aim was to reconstruct phylogenetic relationships to test monophyly and elucidate relationships within the genus. Our results indicated that Ammodendron, Ammothamnus, and Echinosophora are embedded within Sophora s.s. and that nine well-supported clades can be recognized within comprise Sophora s.l. Ancestral character state estimation revealed that the most recent common ancestor of Sophora s.l. was a deciduous shrub that lacks rhizome spines and has unwinged legumes. Divergence times estimation and ancestral area reconstruction showed that Sophora s.l. originated in Central Asia and/or adjacent Southeast China in the early Oligocene (ca. 31 Mya) and dispersed from these regions into East and South Asia's adjacent areas and North America via the Bering land bridge. The analyses also supported a South American origin for S. sect. Edwardsia, which experienced rapid radiation with its major lineages diversifying over a relatively narrow timescale (8 Mya).

Ultraconserved elements-based phylogenomic systematics of the snake superfamily Elapoidea, with the description of a new Afro-Asian family

The highly diverse snake superfamily Elapoidea is considered to be a classic example of ancient, rapid radiation. Such radiations are challenging to fully resolve phylogenetically, with the highly diverse Elapoidea a case in point. Previous attempts at inferring a phylogeny of elapoids produced highly incongruent estimates of their evolutionary relationships, often with very low statistical support. We sought to resolve this situation by sequencing over 4,500 ultraconserved element loci from multiple representatives of every elapoid family/subfamily level taxon and inferring their phylogenetic relationships with multiple methods. Concatenation and multispecies coalescent based species trees yielded largely congruent and well-supported topologies. Hypotheses of a hard polytomy were not retained for any deep branches. Our phylogenies recovered Cyclocoridae and Elapidae as diverging early within Elapoidea. The Afro-Malagasy radiation of elapoid snakes, classified as multiple subfamilies of an inclusive Lamprophiidae by some earlier authors, was found to be monophyletic in all analyses. The genus Micrelaps was consistently recovered as sister to Lamprophiidae. We establish a new family, Micrelapidae fam. nov., for Micrelaps and assign Brachyophis to this family based on cranial osteological synapomorphy. We estimate that Elapoidea originated in the early Eocene and rapidly diversified into all the major lineages during this epoch. Ecological opportunities presented by the post-Cretaceous-Paleogene mass extinction event may have promoted the explosive radiation of elapoid snakes.

Diversity, Form, and Postembryonic Development of Paleozoic Insects

While Mesozoic, Paleogene, and Neogene insect faunas greatly resemble the modern one, the Paleozoic fauna provides unique insights into key innovations in insect evolution, such as the origin of wings and modifications of postembryonic development including holometaboly. Deep-divergence estimates suggest that the majority of contemporary insect orders originated in the Late Paleozoic, but these estimates reflect divergences between stem groups of each lineage rather than the later appearance of the crown groups. The fossil record shows the initial radiations of the extant hyperdiverse clades during the Early Permian, as well as the specialized fauna present before the End Permian mass extinction. This review summarizes the recent discoveries related to the documented diversity of Paleozoic hexapods, as well as current knowledge about what has actually been verified from fossil evidence as it relates to postembryonic development and the morphology of different body parts.

Phylogeny and biogeography of the northern temperate genus Dracocephalum s.l. (Lamiaceae)

The northern temperate genus Dracocephalum consists of approximately 70 species mainly distributed in the steppe-desert biomes of Central and West Asia and the alpine region of the Qinghai-Tibetan Plateau (QTP). Previous work has shown that Dracocephalum is not monophyletic and might include Hyssopus and Lallemantia. This study attempts to clarify the phylogenetic relationships, diversification patterns, and the biogeographical history of the three genera (defined as Dracocephalum s.l.). Based on a sampling of 66 taxa comprising more than 80% from extant species of Dracocephalum s.l., morphological, phylogenetic (maximum parsimony, likelihood, and Bayesian inference based on nuclear ITS and ETS, plastid rpl32-trnL, trnL-trnF, ycf1, and ycf1-rps15, and two low-copy nuclear markers AT3G09060 and AT1G09680), molecular dating, diversification, and ancestral range estimation analyses were carried out. Our results demonstrate that both Hyssopus and Lallemantia are embedded within Dracocephalum and nine well-supported clades can be recognized within Dracocephalum s.l. Analyses of divergence times suggest that the genus experienced an early rapid radiation during the middle to late Miocene with major lineages diversifying within a relatively narrow timescale. Ancestral area reconstruction analyses indicate that Dracocephalum s.l. originated in Central and West Asia and southern Siberia, and dispersed from Central and West Asia into the QTP and adjacent areas twice independently during the Pliocene. The aridification of the Asian interior possibly promoted the rapid radiation of Dracocephalum within this region, and the uplift of the QTP appears to have triggered the dispersal and recent rapid diversification of the genus in the QTP and adjacent regions. Combining molecular phylogenetic and morphological evidence, a revised infrageneric classification of Dracocephalum s.l. is proposed, which recognizes nine sections within the genus.

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

Rapid radiation and rampant reticulation: Phylogenomics of South American Liolaemus lizards

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction-site associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.

Increasing 28 mitogenomes of Ephemeroptera, Odonata and Plecoptera support the Chiastomyaria hypothesis with three different outgroup combinations

BACKGROUND

The phylogenetic relationships of Odonata (dragonflies and damselflies) and Ephemeroptera (mayflies) remain unresolved. Different researchers have supported one of three hypotheses (Palaeoptera, Chiastomyaria or Metapterygota) based on data from different morphological characters and molecular markers, sometimes even re-assessing the same transcriptomes or mitochondrial genomes. The appropriate choice of outgroups and more taxon sampling is thought to eliminate artificial phylogenetic relationships and obtain an accurate phylogeny. Hence, in the current study, we sequenced 28 mt genomes from Ephemeroptera, Odonata and Plecoptera to further investigate phylogenetic relationships, the probability of each of the three hypotheses, and to examine mt gene arrangements in these species. We selected three different combinations of outgroups to analyze how outgroup choice affected the phylogenetic relationships of Odonata and Ephemeroptera.

METHODS

Mitochondrial genomes from 28 species of mayflies, dragonflies, damselflies and stoneflies were sequenced. We used Bayesian inference (BI) and Maximum likelihood (ML) analyses for each dataset to reconstruct an accurate phylogeny of these winged insect orders. The effect of outgroup choice was assessed by separate analyses using three outgroups combinations: (a) four bristletails and three silverfish as outgroups, (b) five bristletails and three silverfish as outgroups, or (c) five diplurans as outgroups.

RESULTS

Among these sequenced mitogenomes we found the gene arrangement IMQM in Heptageniidae (Ephemeroptera), and an inverted and translocated tRNA-Ile between the 12S RNA gene and the control region in Ephemerellidae (Ephemeroptera). The IMQM gene arrangement in Heptageniidae (Ephemeroptera) can be explained via the tandem-duplication and random loss model, and the transposition and inversion of tRNA-Ile genes in Ephemerellidae can be explained through the recombination and tandem duplication-random loss (TDRL) model. Our phylogenetic analysis strongly supported the Chiastomyaria hypothesis in three different outgroup combinations in BI analyses. The results also show that suitable outgroups are very important to determining phylogenetic relationships in the rapid evolution of insects especially among Ephemeroptera and Odonata. The mt genome is a suitable marker to investigate the phylogeny of inter-order and inter-family relationships of insects but outgroup choice is very important for deriving these relationships among winged insects. Hence, we must carefully choose the correct outgroup in order to discuss the relationships of Ephemeroptera and Odonata.

Revisiting the phylogenetic predicament of the genus Huia (Amphibia: Ranidae) using molecular data and tadpole morphology

Despite a considerable amount of research, the systematics of the ranid genus Huia have remained unresolved, mostly owing to insufficient sampling and morphological similarities. As currently circumscribed, Huia consists of five species, but multiple studies have consistently demonstrated that it is not a monophyletic genus. However, no study has approached the problem with adequate data and provided a systematically sound solution, leaving the genus to languish in a classification that is phylogenetically incoherent. We generated the most comprehensive sampling of Huia to date, based on extensive fieldwork in Java and Sumatra. Using an integrative taxonomy framework, we analysed four mitochondrial and two nuclear markers and, in conjunction with tadpole morphology, investigated the phylogenetics of Huia and its congeners. Corroborating previous studies, Huia is recovered as a paraphyletic group. Huia cavitympanum emerges as the sister taxon to Meristogenys. The remaining members of Huia form a monophyletic group, sister to the H. cavitympanum + Meristogenys clade. Our extensive geographical sampling in Sumatra and Java reveals multiple highly divergent lineages that potentially represent undescribed species. Using our expanded molecular and morphological dataset, we resolve the paraphyly of Huia by restricting the genus to its type species and propose a new genus to accommodate the strongly supported clade of Sumatran and Javan populations previously belonging to Huia.

Phylogenomics of Ichneumonoidea (Hymenoptera) and implications for evolution of mode of parasitism and viral endogenization

Ichneumonoidea is one of the most diverse lineages of animals on the planet with >48,000 described species and many more undescribed. Parasitoid wasps of this superfamily are mostly beneficial insects that attack and kill other arthropods and are important for understanding diversification and the evolution of life history strategies related to parasitoidism. Further, some lineages of parasitoids within Ichneumonoidea have acquired endogenous virus elements (EVEs) that are permanently a part of the wasp's genome and benefit the wasp through host immune disruption and behavioral control. Unfortunately, understanding the evolution of viral acquisition, parasitism strategies, diversification, and host immune disruption mechanisms, is deeply limited by the lack of a robust phylogenetic framework for Ichneumonoidea. Here we design probes targeting 541 genes across 91 taxa to test phylogenetic relationships, the evolution of parasitoid strategies, and the utility of probes to capture polydnavirus genes across a diverse array of taxa. Phylogenetic relationships among Ichneumonoidea were largely well resolved with most higher-level relationships maximally supported. We noted codon use biases between the outgroups, Braconidae, and Ichneumonidae and within Pimplinae, which were largely solved through analyses of amino acids rather than nucleotide data. These biases may impact phylogenetic reconstruction and caution for outgroup selection is recommended. Ancestral state reconstructions were variable for Braconidae across analyses, but consistent for reconstruction of idiobiosis/koinobiosis in Ichneumonidae. The data suggest many transitions between parasitoid life history traits across the whole superfamily. The two subfamilies within Ichneumonidae that have polydnaviruses are supported as distantly related, providing strong evidence for two independent acquisitions of ichnoviruses. Polydnavirus capture using our designed probes was only partially successful and suggests that more targeted approaches would be needed for this strategy to be effective for surveying taxa for these viral genes. In total, these data provide a robust framework for the evolution of Ichneumonoidea.

The Perfect Storm: Gene Tree Estimation Error, Incomplete Lineage Sorting, and Ancient Gene Flow Explain the Most Recalcitrant Ancient Angiosperm Clade, Malpighiales

The genomic revolution offers renewed hope of resolving rapid radiations in the Tree of Life. The development of the multispecies coalescent model and improved gene tree estimation methods can better accommodate gene tree heterogeneity caused by incomplete lineage sorting (ILS) and gene tree estimation error stemming from the short internal branches. However, the relative influence of these factors in species tree inference is not well understood. Using anchored hybrid enrichment, we generated a data set including 423 single-copy loci from 64 taxa representing 39 families to infer the species tree of the flowering plant order Malpighiales. This order includes 9 of the top 10 most unstable nodes in angiosperms, which have been hypothesized to arise from the rapid radiation during the Cretaceous. Here, we show that coalescent-based methods do not resolve the backbone of Malpighiales and concatenation methods yield inconsistent estimations, providing evidence that gene tree heterogeneity is high in this clade. Despite high levels of ILS and gene tree estimation error, our simulations demonstrate that these two factors alone are insufficient to explain the lack of resolution in this order. To explore this further, we examined triplet frequencies among empirical gene trees and discovered some of them deviated significantly from those attributed to ILS and estimation error, suggesting gene flow as an additional and previously unappreciated phenomenon promoting gene tree variation in Malpighiales. Finally, we applied a novel method to quantify the relative contribution of these three primary sources of gene tree heterogeneity and demonstrated that ILS, gene tree estimation error, and gene flow contributed to 10.0$\%$, 34.8$\%$, and 21.4$\%$ of the variation, respectively. Together, our results suggest that a perfect storm of factors likely influence this lack of resolution, and further indicate that recalcitrant phylogenetic relationships like the backbone of Malpighiales may be better represented as phylogenetic networks. Thus, reducing such groups solely to existing models that adhere strictly to bifurcating trees greatly oversimplifies reality, and obscures our ability to more clearly discern the process of evolution. [Coalescent; concatenation; flanking region; hybrid enrichment, introgression; phylogenomics; rapid radiation, triplet frequency.].

Larger, unfiltered datasets are more effective at resolving phylogenetic conflict: Introns, exons, and UCEs resolve ambiguities in Golden-backed frogs (Anura: Ranidae; genus Hylarana)

Using FrogCap, a recently-developed sequence-capture protocol, we obtained >12,000 highly informative exons, introns, and ultraconserved elements (UCEs), which we used to illustrate variation in evolutionary histories of these classes of markers, and to resolve long-standing systematic problems in Southeast Asian Golden-backed frogs of the genus-complex Hylarana. We also performed a comprehensive suite of analyses to assess the relative performance of different genetic markers, data filtering strategies, tree inference methods, and different measures of branch support. To reduce gene tree estimation error, we filtered the data using different thresholds of taxon completeness (missing data) and parsimony informative sites (PIS). We then estimated species trees using concatenated datasets and Maximum Likelihood (IQ-TREE) in addition to summary (ASTRAL-III), distance-based (ASTRID), and site-based (SVDQuartets) multispecies coalescent methods. Topological congruence and branch support were examined using traditional bootstrap, local posterior probabilities, gene concordance factors, quartet frequencies, and quartet scores. Our results did not yield a single concordant topology. Instead, introns, exons, and UCEs clearly possessed different phylogenetic signals, resulting in conflicting, yet strongly-supported phylogenetic estimates. However, a combined analysis comprising the most informative introns, exons, and UCEs converged on a similar topology across all analyses, with the exception of SVDQuartets. Bootstrap values were consistently high despite high levels of incongruence and high proportions of gene trees supporting conflicting topologies. Although low bootstrap values did indicate low heuristic support, high bootstrap support did not necessarily reflect congruence or support for the correct topology. This study reiterates findings of some previous studies, which demonstrated that traditional bootstrap values can produce positively misleading measures of support in large phylogenomic datasets. We also showed a remarkably strong positive relationship between branch length and topological congruence across all datasets, implying that very short internodes remain a challenge to resolve, even with orders of magnitude more data than ever before. Overall, our results demonstrate that more data from unfiltered or combined datasets produced superior results. Although data filtering reduced gene tree incongruence, decreased amounts of data also biased phylogenetic estimation. A point of diminishing returns was evident, at which higher congruence (from more stringent filtering) at the expense of amount of data led to topological error as assessed by comparison to more complete datasets across different genomic markers. Additionally, we showed that applying a parameter-rich model to a partitioned analysis of concatenated data produces better results compared to unpartitioned, or even partitioned analysis using model selection. Despite some lingering uncertainties, a combined analysis of our genomic data and sequences supplemented from GenBank (on the basis of a few gene regions) revealed highly supported novel systematic arrangements. Based on these new findings, we transfer Amnirana nicobariensis into the genus Indosylvirana; and I. milleti and Hylarana celebensis to the genus Papurana. We also provisionally place H. attigua in the genus Papurana pending verification from positively identified (voucher substantiated) samples.

The Mitochondrial Genome of the Phytopathogenic Fungus Bipolaris sorokiniana and the Utility of Mitochondrial Genome to Infer Phylogeny of Dothideomycetes

A number of species in Bipolaris are important plant pathogens. Due to a limited number of synapomorphic characters, it is difficult to perform species identification and to estimate phylogeny of Bipolaris based solely on morphology. In this study, we sequenced the complete mitochondrial genome of Bipolaris sorokiniana, and presented the detailed annotation of the genome. The B. sorokiniana mitochondrial genome is 137,775 bp long, and contains two ribosomal RNA genes, 12 core protein-coding genes, 38 tRNA genes. In addition, two ribosomal protein genes (rps3 gene and rps5 gene) and the fungal mitochondrial RNase P gene (rnpB) are identified. The large genome size is mostly determined by the presence of numerous intronic and intergenic regions. A total of 28 introns are inserted in eight core protein-coding genes. Together with the published mitochondrial genome sequences, we conducted a preliminary phylogenetic inference of Dothideomycetes under various datasets and substitution models. The monophyly of Capnodiales, Botryosphaeriales and Pleosporales are consistently supported in all analyses. The Venturiaceae forms an independent lineage, with a distant phylogenetic relationship to Pleosporales. At the family level, the Mycosphaerellaceae, Botryosphaeriaceae. Phaeosphaeriaceae, and Pleosporaceae are recognized in the majority of trees.

The Mitochondrial Genome of Amara aulica (Coleoptera, Carabidae, Harpalinae) and Insights into the Phylogeny of Ground Beetles

Carabidae are one of the most species-rich families of beetles, comprising more than 40,000 described species worldwide. Forty-three complete or partial mitochondrial genomes (mitogenomes) from this family have been published in GenBank to date. In this study, we sequenced a nearly complete mitogenome of Amara aulica (Carabidae), using a next-generation sequencing method. This mitogenome was 16,646 bp in length, which encoded the typical 13 mitochondrial protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a putative control region. Combining with the published mitogenomes of Carabidae and five outgroup species from Trachypachidae, Gyrinidae and Dytiscidae, we performed phylogenetic estimates under maximum likelihood and Bayesian inference criteria to investigate the phylogenetic relationships of carabid beetles. The results showed that the family Carabidae was a non-monophyletic assemblage. The subfamilies Cicindelinae, Elaphrinae, Carabinae, Trechinae and Harpalinae were recovered as monophyletic groups. Moreover, the clade (Trechinae + (Brachininae + Harpalinae)) was consistently recovered in all analyses.

Molecular phylogenetics disentangles the limits of a group of related Neotropical doryctine parasitoid wasp genera (Hymenoptera : Braconidae)

The Doryctinae is one of the most specious subfamilies of the mainly parasitoid wasp family Braconidae. Members of this subfamily are characterised by having a very heterogeneous external morphology, which has made establishing the limits of several of its genera extremely difficult. The Neotropical Callihormius Ashmead, Platydoryctes Barbalho & Penteado-Dias, and Aphelopsia Marsh are prime examples of this, since they show extensive interspecific morphological variation but lack consistent diagnostic features among them. We carried out a comprehensive phylogenetic study among several representative species belonging to the above genera as well as other morphologically similar doryctine taxa based on three nuclear and two mitochondrial (mt)DNA sequence markers. Based on the relationships recovered, we assessed the generic limits among the examined taxa and investigated the times of origin and diversification in this group. Our best estimate of phylogeny did not recover Callihormius as monophyletic, and thus we propose a taxonomic arrangement where we elevate the status of Ca. (Dmitriohormius) Belokobylskij, Zaldívar-Riverón & Coronado-Blanco to genus level, stat. nov., and also describe four new genera: Macrometasoma, gen. nov., Caputlenis, gen. nov., Platyhormius, gen. nov. and Caputrugosus, gen. nov. Platydoryctes is delimited morphologically, and the body flatness is proposed to have independently evolved more than once within the group of genera examined. Aphelopsia was not recovered as monophyletic, though the relationships involved were weakly supported and therefore we maintain its species composition. Panama Marsh, stat. rev., which was considered a junior synonym of Leluthia Cameron, is resurrected. Basal relationships in the phylogeny had very short branch lengths, which could be due to ancient rapid radiation events that occurred during the early evolution of the group. The estimated times of divergence showed that the clade containing Callihormius and its related genera originated between the late Oligocene and middle Miocene, whereas its early diversification events probably occurred between the middle to late Miocene.

Maximize Resolution or Minimize Error? Using Genotyping-By-Sequencing to Investigate the Recent Diversification of Helianthemum (Cistaceae)

A robust phylogenetic framework, in terms of extensive geographical and taxonomic sampling, well-resolved species relationships and high certainty of tree topologies and branch length estimations, is critical in the study of macroevolutionary patterns. Whereas Sanger sequencing-based methods usually recover insufficient phylogenetic signal, especially in recently diversified lineages, reduced-representation sequencing methods tend to provide well-supported phylogenetic relationships, but usually entail remarkable bioinformatic challenges due to the inherent trade-off between the number of SNPs and the magnitude of associated error rates. The genus Helianthemum (Cistaceae) is a species-rich and taxonomically complex Palearctic group of plants that diversified mainly since the Upper Miocene. It is a challenging case study since previous attempts using Sanger sequencing were unable to resolve the intrageneric phylogenetic relationships. Aiming to obtain a robust phylogenetic reconstruction based on genotyping-by-sequencing (GBS), we established a rigorous methodological workflow in which we i) explored how variable settings during dataset assembly have an impact on error rates and on the degree of resolution under concatenation and coalescent approaches, ii) assessed the effect of two extreme parameter configurations (minimizing error rates vs. maximizing phylogenetic resolution) on tree topology and branch lengths, and iii) evaluated the effects of these two configurations on estimates of divergence times and diversification rates. Our analyses produced highly supported topologically congruent phylogenetic trees for both configurations. However, minimizing error rates did produce more reliable branch lengths, critically affecting the accuracy of downstream analyses (i.e. divergence times and diversification rates). In addition to recommending a revision of intrageneric systematics, our results enabled us to identify three highly diversified lineages in Helianthemum in contrasting geographical areas and ecological conditions, which started radiating in the Upper Miocene.

Evolutionary history of Polyneoptera and its implications for our understanding of early winged insects

Polyneoptera represents one of the major lineages of winged insects, comprising around 40,000 extant species in 10 traditional orders, including grasshoppers, roaches, and stoneflies. Many important aspects of polyneopteran evolution, such as their phylogenetic relationships, changes in their external appearance, their habitat preferences, and social behavior, are unresolved and are a major enigma in entomology. These ambiguities also have direct consequences for our understanding of the evolution of winged insects in general; for example, with respect to the ancestral habitats of adults and juveniles. We addressed these issues with a large-scale phylogenomic analysis and used the reconstructed phylogenetic relationships to trace the evolution of 112 characters associated with the external appearance and the lifestyle of winged insects. Our inferences suggest that the last common ancestors of Polyneoptera and of the winged insects were terrestrial throughout their lives, implying that wings did not evolve in an aquatic environment. The appearance of the first polyneopteran insect was mainly characterized by ancestral traits such as long segmented abdominal appendages and biting mouthparts held below the head capsule. This ancestor lived in association with the ground, which led to various specializations including hardened forewings and unique tarsal attachment structures. However, within Polyneoptera, several groups switched separately to a life on plants. In contrast to a previous hypothesis, we found that social behavior was not part of the polyneopteran ground plan. In other traits, such as the biting mouthparts, Polyneoptera shows a high degree of evolutionary conservatism unique among the major lineages of winged insects.

Noise and biases in genomic data may underlie radically different hypotheses for the position of Iguania within Squamata

Squamate reptiles are a major component of vertebrate biodiversity whose crown-clade traces its origin to a narrow window of time in the Mesozoic during which the main subclades diverged in rapid succession. Deciphering phylogenetic relationships among these lineages has proven challenging given the conflicting signals provided by genomic and phenomic data. Most notably, the placement of Iguania has routinely differed between data sources, with morphological evidence supporting a sister relationship to the remaining squamates (Scleroglossa hypothesis) and molecular data favoring a highly nested position alongside snakes and anguimorphs (Toxicofera hypothesis). We provide novel insights by generating an expanded morphological dataset and exploring the presence of phylogenetic signal, noise, and biases in molecular data. Our analyses confirm the presence of strong conflicting signals for the position of Iguania between morphological and molecular datasets. However, we also find that molecular data behave highly erratically when inferring the deepest branches of the squamate tree, a consequence of limited phylogenetic signal to resolve this ancient radiation with confidence. This, in turn, seems to result from a rate of evolution that is too high for historical signals to survive to the present. Finally, we detect significant systematic biases, with iguanians and snakes sharing faster rates of molecular evolution and a similarly biased nucleotide composition. A combination of scant phylogenetic signal, high levels of noise, and the presence of systematic biases could result in the misplacement of Iguania. We regard this explanation to be at least as plausible as the complex scenario of convergence and reversals required for morphological data to be misleading. We further evaluate and discuss the utility of morphological data to resolve ancient radiations, as well as its impact in combined-evidence phylogenomic analyses, with results relevant for the assessment of evidence and conflict across the Tree of Life.

Multispecies coalescent analysis confirms standing phylogenetic instability in Hexapoda

The multispecies coalescent (MSC) has been increasingly used in phylogenomic analyses due to the accommodation of gene tree topological heterogeneity by taking into account population-level processes, such as incomplete lineage sorting. In this sense, the phylogeny of insect species, which are characterized by their large effective population sizes, is suitable for a coalescent-based analysis. Furthermore, studies so far recovered short internal branches at early divergences of the insect tree of life, indicating fast evolutionary radiations that increase the probability of incomplete lineage sorting in deep time. Here, we investigated the performance of the MSC for a phylogenomic data set of hexapods compiled by Misof et al. (2014, Science 346:763). Our analysis recovered the monophyly of most insect orders, and major phylogenetic relationships were in agreement with current insect systematics. We identified, however, some evolutionary associations that were consistently problematic. Most noticeable, Hexapod monophyly was disrupted by the sister group relationship between the remiped crustacean and Insecta. Additionally, the interordinal relationships within Polyneoptera and Neuropteroidea were found to be phylogenetically unstable. We show that these controversial phylogenetic arrangements were also poorly supported by previous analyses, and therefore, we evaluated their robustness to stochastic errors from sampling sites and terminals, confirming standing problems in hexapod phylogeny in the genomics age.

Histochemistry, immunohistochemistry and cytochemistry of the anterior midgut region of the stingless bee Melipona quadrifasciata and honey bee Apis mellifera (Hymenoptera: Apidae)

The anterior midgut region of stingless bees is anatomically differentiated with tall and narrow cells, whereas in other social and solitary bees this anatomical gut region is lacking. The objective of the present study was to describe the histochemistry, immunohistochemistry and cytochemistry of the anterior midgut region of the stingless bee Melipona quadrifasciata in comparison with the honey bee Apis mellifera. The anterior midgut region of both species was evaluated for identification of the enzymes β- galactosidase, glucose-6-phosphatase, acid phosphatase, and alkaline phosphatase, the membrane transporter aquaporin, the hormone FMRF-amide, and lysosomes. Histology of the anterior midgut region showed that this region in M. quadrifasciata workers did not present external folds of the wall, whereas the following midgut wall presented many. In A. mellifera, folds in the midgut wall occur starting from the fore- midgut transition region. Despite these morphological differences, the tests evaluated were similar in both species. β-galactosidase was not found in the anterior midgut cells. Glucose-6-phosphatase and acid phosphatase occurred in the apical region of the gut epithelium. Alkaline phosphatase occurred in vesicles in apical cytoplasm and in the basal plasma membrane infoldings of the epithelial cells. Aquaporin was found in the basal region of the midgut epithelium and in the associated visceral muscles. FMRF-amide was found only in nerve endings in the anterior midgut region. All cells in the anterior midgut region were rich in lysosomes. These results suggest that in both bee species, although they have anatomically different anterior midgut regions, these regions present high metabolic activity and function in cellular homeostasis, lipid absorption and are under neurohormone control.

Reanalyzing the Palaeoptera problem - The origin of insect flight remains obscure

The phylogenetic relationships of the winged insect lineages - mayflies (Ephemeroptera), damselflies and dragonflies (Odonata), and all other winged insects (Neoptera) - are still controversial with three hypotheses supported by different datasets: Palaeoptera, Metapterygota and Chiastomyaria. Here, we reanalyze available phylogenomic data with a focus on detecting confounding and alternative signal. In this context, we provide a framework to quantitatively evaluate and assess incongruent molecular phylogenetic signal inherent in phylogenomic datasets. Despite overall support for the Palaeoptera hypothesis, we also found considerable signal for Chiastomyaria, which is not easily detectable by standardized tree inference approaches. Analyses of the accumulation of signal across gene partitions showed that signal accumulates gradually. However, even in case signal only slightly supported one over the other hypothesis, topologies inferred from large datasets switch from statistically strongly supported Palaeoptera to strongly supported Chiastomyaria. From a morphological point of view, Palaeoptera currently appears to be the best-supported hypothesis; however, recent analyses were restricted to head characters. Phylogenetic approaches covering all organ systems including analyses of potential functional or developmental convergence are still pending so that the Palaeoptera problem has to be considered an open question in insect systematics.

Systematics, Biology, and Evolution of Microgastrine Parasitoid Wasps

The braconid parasitoid wasp subfamily Microgastrinae is perhaps the most species-rich subfamily of animals on Earth. Despite their small size, they are familiar to agriculturalists and field ecologists alike as one of the principal groups of natural enemies of caterpillars feeding on plants. Their abundance and nearly ubiquitous terrestrial distribution, their intricate interactions with host insects, and their historical association with mutualistic polydnaviruses have all contributed to Microgastrinae becoming a key group of organisms for studying parasitism, parasitoid genomics, and mating biology. However, these rich sources of data have not yet led to a robust genus-level classification of the group, and some taxonomic confusion persists as a result. We present the current status of understanding of the general biology, taxonomic history, diversity, geographical patterns, host relationships, and phylogeny of Microgastrinae as a stimulus and foundation for further study. Current progress in elucidating the biology and taxonomy of this important group is rapid and promises a revolution in the classification of these wasps in the near future.

Mitochondrial genomes of blister beetles (Coleoptera, Meloidae) and two large intergenic spacers in Hycleus genera

BACKGROUND

Insect mitochondrial genomes (mitogenomes) exhibit high diversity in some lineages. The gene rearrangement and large intergenic spacer (IGS) have been reported in several Coleopteran species, although very little is known about mitogenomes of Meloidae.

RESULTS

We determined complete or nearly complete mitogenomes of seven meloid species. The circular genomes encode 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and two ribosomal RNAs (rRNAs), and contain a control region, with gene arrangement identical to the ancestral type for insects. The evolutionary rates of all PCGs indicate that their evolution is based on purifying selection. The comparison of tRNA secondary structures indicates diverse substitution patterns in Meloidae. Remarkably, all mitogenomes of the three studied Hycleus species contain two large intergenic spacers (IGSs). IGS1 is located between trnW and trnC, including a 9 bp consensus motif. IGS2 is located between trnS2 (UCN) and nad1, containing discontinuous repeats of a pentanucleotide motif and two 18-bp repeat units in both ends. To date, IGS2 is found only in genera Hycleus across all published Coleopteran mitogenomes. The duplication/random loss model and slipped-strand mispairing are proposed as evolutionary mechanisms for the two IGSs (IGS1, IGS2). The phylogenetic analyses using MrBayes, RAxML, and PhyloBayes methods based on nucleotide and amino acid datasets of 13 PCGs from all published mitogenomes of Tenebrionoids, consistently recover the monophylies of Meloidae and Tenebrionidae. Within Meloidae, the genus Lytta clusters with Epicauta rather than with Mylabris. Although data collected thus far could not resolve the phylogenetic relationships within Meloidae, this study will assist in future mapping of the Meloidae phylogeny.

CONCLUSIONS

This study presents mitogenomes of seven meloid beetles. New mitogenomes retain the genomic architecture of the Coleopteran ancestor, but contain two IGSs in the three studied Hycleus species. Comparative analyses of two IGSs suggest that their evolutionary mechanisms are duplication/random loss model and slipped-strand mispairing.

Phylogenomic analysis of a rapid radiation of misfit fishes (Syngnathiformes) using ultraconserved elements

Phylogenetics is undergoing a revolution as large-scale molecular datasets reveal unexpected but repeatable rearrangements of clades that were previously thought to be disparate lineages. One of the most unusual clades of fishes that has been found using large-scale molecular datasets is an expanded Syngnathiformes including traditional long-snouted syngnathiform lineages (Aulostomidae, Centriscidae, Fistulariidae, Solenostomidae, Syngnathidae), as well as a diverse set of largely benthic-associated fishes (Callionymoidei, Dactylopteridae, Mullidae, Pegasidae) that were previously dispersed across three orders. The monophyly of this surprising clade of fishes has been upheld by recent studies utilizing both nuclear and mitogenomic data, but the relationships among major lineages within Syngnathiformes remain ambiguous; previous analyses have inconsistent topologies and are plagued by low support at deep divergences between the major lineages. In this study, we use a dataset of ultraconserved elements (UCEs) to conduct the first phylogenomic study of Syngnathiformes. UCEs have been effective markers for resolving deep phylogenetic relationships in fishes and, combined with increased taxon sampling, we expected UCEs to resolve problematic syngnathiform relationships. Overall, UCEs were effective at resolving relationships within Syngnathiformes at a range of evolutionary timescales. We find consistent support for the monophyly of traditional long-snouted syngnathiform lineages (Aulostomidae, Centriscidae, Fistulariidae, Solenostomidae, Syngnathidae), which better agrees with morphological hypotheses than previously published topologies from molecular data. This result was supported by all Bayesian and maximum likelihood analyses, was robust to differences in matrix completeness and potential sources of bias, and was highly supported in coalescent-based analyses in ASTRAL when matrices were filtered to contain the most phylogenetically informative loci. While Bayesian and maximum likelihood analyses found support for a benthic-associated clade (Callionymidae, Dactylopteridae, Mullidae, and Pegasidae) as sister to the long-snouted clade, this result was not replicated in the ASTRAL analyses. The base of our phylogeny is characterized by short internodes separating major syngnathiform lineages and is consistent with the hypothesis of an ancient rapid radiation at the base of Syngnathiformes. Syngnathiformes therefore present an exciting opportunity to study patterns of morphological variation and functional innovation arising from rapid but ancient radiation.

Evolutionary dynamics of triosephosphate isomerase gene intron location pattern in Metazoa: A new perspective on intron evolution in animals

Intron evolution, including its dynamics in the evolutionary transitions and diversification of eukaryotes, remains elusive. Inadequate taxon sampling due to data shortage, unclear phylogenetic framework, and inappropriate outgroup application might be among the causes. Besides, the integrity of all the introns within a gene was often neglected previously. Taking advantage of the ancient conserved triosephosphate isomerase gene (tim), the relatively robust phylogeny of Metazoa, and choanoflagellates as outgroup, the evolutionary dynamics of tim intron location pattern (ILP) in Metazoa was investigated. From 133 representative species of ten phyla, 30 types of ILPs were identified. A most common one, which harbors the maximum six intron positions, is deduced to be the common ancestral tim ILP of Metazoa, which almost had formed in their protozoan ancestor and was surprisingly retained and passed down till to each ancestors of metazoan phyla. In the subsequent animal diversification, it underwent different evolutionary trajectories: within Deuterostomia, it was almost completely retained only with changes in a few species with relatively recently fast-evolving histories, while within the rapidly radiating Protostomia, besides few but remarkable retention, it usually displayed extensive intron losses and a few gains. Therefore, a common ancestral exon-intron arrangement pattern of an animal gene is definitely discovered; besides the 'intron-rich view' of early animal genes being confirmed, the novel insight that high exon-intron re-arrangements of genes seem to be associated with the relatively recently rapid evolution of lineages/species/genomes but have no correlation with the ancient major evolutionary transitions in animal evolution, is revealed.

Computational biomechanics changes our view on insect head evolution

Despite large-scale molecular attempts, the relationships of the basal winged insect lineages dragonflies, mayflies and neopterans, are still unresolved. Other data sources, such as morphology, suffer from unclear functional dependencies of the structures considered, which might mislead phylogenetic inference. Here, we assess this problem by combining for the first time biomechanics with phylogenetics using two advanced engineering techniques, multibody dynamics analysis and finite-element analysis, to objectively identify functional linkages in insect head structures which have been used traditionally to argue basal winged insect relationships. With a biomechanical model of unprecedented detail, we are able to investigate the mechanics of morphological characters under biologically realistic load, i.e. biting. We show that a range of head characters, mainly ridges, endoskeletal elements and joints, are indeed mechanically linked to each other. An analysis of character state correlation in a morphological data matrix focused on head characters shows highly significant correlation of these mechanically linked structures. Phylogenetic tree reconstruction under different data exclusion schemes based on the correlation analysis unambiguously supports a sistergroup relationship of dragonflies and mayflies. The combination of biomechanics and phylogenetics as it is proposed here could be a promising approach to assess functional dependencies in many organisms to increase our understanding of phenotypic evolution.

Microtomography of the Baltic amber tick Ixodes succineus reveals affinities with the modern Asian disease vector Ixodes ovatus

BACKGROUND

Fossil ticks are extremely rare and Ixodes succineus Weidner, 1964 from Eocene (ca. 44-49 Ma) Baltic amber is one of the oldest examples of a living hard tick genus (Ixodida: Ixodidae). Previous work suggested it was most closely related to the modern and widespread European sheep tick Ixodes ricinus (Linneaus, 1758).

RESULTS

Restudy using phase contrast synchrotron x-ray tomography yielded images of exceptional quality. These confirm the fossil's referral to Ixodes Latreille, 1795, but the characters resolved here suggest instead affinities with the Asian subgenus Partipalpiger Hoogstraal et al., 1973 and its single living (and medically significant) species Ixodes ovatus Neumann, 1899. We redescribe the amber fossil here as Ixodes (Partipalpiger) succineus.

CONCLUSIONS

Our data suggest that Ixodes ricinus is unlikely to be directly derived from Weidner's amber species, but instead reveals that the Partipalpiger lineage was originally more widely distributed across the northern hemisphere. The closeness of Ixodes (P.) succineus to a living vector of a wide range of pathogens offers the potential to correlate its spatial and temporal position (northern Europe, nearly 50 million years ago) with the estimated origination dates of various tick-borne diseases.

Progress, pitfalls and parallel universes: a history of insect phylogenetics

The phylogeny of insects has been both extensively studied and vigorously debated for over a century. A relatively accurate deep phylogeny had been produced by 1904. It was not substantially improved in topology until recently when phylogenomics settled many long-standing controversies. Intervening advances came instead through methodological improvement. Early molecular phylogenetic studies (1985-2005), dominated by a few genes, provided datasets that were too small to resolve controversial phylogenetic problems. Adding to the lack of consensus, this period was characterized by a polarization of philosophies, with individuals belonging to either parsimony or maximum-likelihood camps; each largely ignoring the insights of the other. The result was an unfortunate detour in which the few perceived phylogenetic revolutions published by both sides of the philosophical divide were probably erroneous. The size of datasets has been growing exponentially since the mid-1980s accompanied by a wave of confidence that all relationships will soon be known. However, large datasets create new challenges, and a large number of genes does not guarantee reliable results. If history is a guide, then the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analysed.

Revisiting the Zingiberales: using multiplexed exon capture to resolve ancient and recent phylogenetic splits in a charismatic plant lineage

The Zingiberales are an iconic order of monocotyledonous plants comprising eight families with distinctive and diverse floral morphologies and representing an important ecological element of tropical and subtropical forests. While the eight families are demonstrated to be monophyletic, phylogenetic relationships among these families remain unresolved. Neither combined morphological and molecular studies nor recent attempts to resolve family relationships using sequence data from whole plastomes has resulted in a well-supported, family-level phylogenetic hypothesis of relationships. Here we approach this challenge by leveraging the complete genome of one member of the order, Musa acuminata, together with transcriptome information from each of the other seven families to design a set of nuclear loci that can be enriched from highly divergent taxa with a single array-based capture of indexed genomic DNA. A total of 494 exons from 418 nuclear genes were captured for 53 ingroup taxa. The entire plastid genome was also captured for the same 53 taxa. Of the total genes captured, 308 nuclear and 68 plastid genes were used for phylogenetic estimation. The concatenated plastid and nuclear dataset supports the position of Musaceae as sister to the remaining seven families. Moreover, the combined dataset recovers known intra- and inter-family phylogenetic relationships with generally high bootstrap support. This is a flexible and cost effective method that gives the broader plant biology community a tool for generating phylogenomic scale sequence data in non-model systems at varying evolutionary depths.

Phylogenomic methods outperform traditional multi-locus approaches in resolving deep evolutionary history: a case study of formicine ants

BACKGROUND

Ultraconserved elements (UCEs) have been successfully used in phylogenomics for a variety of taxa, but their power in phylogenetic inference has yet to be extensively compared with that of traditional Sanger sequencing data sets. Moreover, UCE data on invertebrates, including insects, are sparse. We compared the phylogenetic informativeness of 959 UCE loci with a multi-locus data set of ten nuclear markers obtained via Sanger sequencing, testing the ability of these two types of data to resolve and date the evolutionary history of the second most species-rich subfamily of ants in the world, the Formicinae.

RESULTS

Phylogenetic analyses show that UCEs are superior in resolving ancient and shallow relationships in formicine ants, demonstrated by increased node support and a more resolved phylogeny. Phylogenetic informativeness metrics indicate a twofold improvement relative to the 10-gene data matrix generated from the identical set of taxa. We were able to significantly improve formicine classification based on our comprehensive UCE phylogeny. Our divergence age estimations, using both UCE and Sanger data, indicate that crown-group Formicinae are older (104-117 Ma) than previously suggested. Biogeographic analyses infer that the diversification of the subfamily has occurred on all continents with no particular hub of cladogenesis.

CONCLUSIONS

We found UCEs to be far superior to the multi-locus data set in estimating formicine relationships. The early history of the clade remains uncertain due to ancient rapid divergence events that are unresolvable even with our genomic-scale data, although this might be largely an effect of several problematic taxa subtended by long branches. Our comparison of divergence ages from both Sanger and UCE data demonstrates the effectiveness of UCEs for dating analyses. This comparative study highlights both the promise and limitations of UCEs for insect phylogenomics, and will prove useful to the growing number of evolutionary biologists considering the transition from Sanger to next-generation sequencing approaches.

Elusive ditrysian phylogeny: an account of combining systematized morphology with molecular data (Lepidoptera)

BACKGROUND

Ditrysia comprise close to 99 % of all butterflies and moths. The evolutionary relationships among the ditrysian superfamilies have received considerable attention in phylogenetic studies based on DNA and transcriptomic data, but the deepest divergences remain for large parts unresolved or contradictory. To obtain complementary insight into the evolutionary history of the clade, and to test previous hypotheses on the subdivision of Ditrysia based on morphology, we examine the morphology of larvae, pupae and adult males and females of 318 taxa representing nearly all ditrysian superfamilies and families. We present the most comprehensive morphological dataset on Ditrysia to date, consisting of over 500 morphological characters. The data are analyzed alone and combined with sequence data (one mitochondrial and seven nuclear protein-coding gene regions, sequenced from 422 taxa). The full dataset consists of 473 exemplar species. Analyses are performed using maximum likelihood methods, and parsimony methods for the morphological dataset. We explore whether combining morphological data and DNA-data can stabilize taxa that are unstable in phylogenetic studies based on genetic data only.

RESULTS

Morphological characters are found phylogenetically informative in resolving apical nodes (superfamilies and families), but characters serving as evidence of relatedness of larger assemblages are few. Results include the recovery of a monophyletic Tineoidea, Sesioidea and Cossoidea, and a stable position for some unstable taxa (e.g. Epipyropidae, Cyclotornidae, Urodoidea + Schreckensteinioidea). Several such taxa, however, remain unstable even though morphological characters indicate a position in the tree (e.g. Immidae). Evidence supporting affinities between clades are suggested, e.g. a novel larval synapomorphy for Tineidae. We also propose the synonymy of Tineodidae with Alucitidae, syn. nov.

CONCLUSIONS

The large morphological dataset provides information on the diversity and distribution of morphological traits in Ditrysia, and can be used in future research on the evolution of these traits, in identification keys and in identification of fossil Lepidoptera. The "backbone" of the phylogeny for Ditrysia remains largely unresolved. As previously proposed as an explanation for the scarcity of molecular signal in resolving the deeper nodes, this may be due to the rapid radiation of Ditrysia in the Cretaceous.

Explosive radiation or uninformative genes? Origin and early diversification of tachinid flies (Diptera: Tachinidae)

Molecular phylogenetic studies at all taxonomic levels often infer rapid radiation events based on short, poorly resolved internodes. While such rapid episodes of diversification are an important and widespread evolutionary phenomenon, much of this poor phylogenetic resolution may be attributed to the continuing widespread use of "traditional" markers (mitochondrial, ribosomal, and some nuclear protein-coding genes) that are often poorly suited to resolve difficult, higher-level phylogenetic problems. Here we reconstruct phylogenetic relationships among a representative set of taxa of the parasitoid fly family Tachinidae and related outgroups of the superfamily Oestroidea. The Tachinidae are one of the most species rich, yet evolutionarily recent families of Diptera, providing an ideal case study for examining the differential performance of loci in resolving phylogenetic relationships and the benefits of adding more loci to phylogenetic analyses. We assess the phylogenetic utility of nine genes including both traditional genes (e.g., CO1 mtDNA, 28S rDNA) and nuclear protein-coding genes newly developed for phylogenetic analysis. Our phylogenetic findings, based on a limited set of taxa, include: a close relationship between Tachinidae and the calliphorid subfamily Polleninae, monophyly of Tachinidae and the subfamilies Exoristinae and Dexiinae, subfamily groupings of Dexiinae+Phasiinae and Tachininae+Exoristinae, and robust phylogenetic placement of the somewhat enigmatic genera Strongygaster, Euthera, and Ceracia. In contrast to poor resolution and phylogenetic incongruence of "traditional genes," we find that a more selective set of highly informative genes is able to more precisely identify regions of the phylogeny that experienced rapid radiation of lineages, while more accurately depicting their phylogenetic context. Although much expanded taxon sampling is necessary to effectively assess the monophyly of and relationships among major tachinid lineages and their relatives, we show that a small number of well-chosen nuclear protein-coding genes can successfully resolve even difficult phylogenetic problems.

Patterns of host plant utilization and diversification in the brush-footed butterflies

Herbivorous insects represent one of the most successful animal radiations known. They occupy a wide range of niches, feed on a great variety of plants, and are species rich; yet the factors that influence their diversification are poorly understood. Host breadth is often cited as a major factor influencing diversification, and, according to the Oscillation Hypothesis, shifts from generalist to specialist feeding states increase the diversification rate for a clade. We explored the relationship between host breadth and diversification within the Nymphalidae (Lepidoptera) and explicitly tested predictions of the Oscillation Hypothesis. We found strong evidence of diversification rate heterogeneity, but no difference in host breadth between clades with a higher diversification rate compared to their sisters. We also found some clades exhibited phylogenetic nonindependence in host breadth and these clades had lower host plant turnover than expected by chance, suggesting host breadth is evolutionarily constrained. Finally, we found that transitions among host breadth categories varied, but the likelihood of reductions in host breadth was greater than that of increases. Our results indicate host breadth is decoupled from diversification rate within the Nymphalidae, and that constraints on diet breadth might play an important role in the evolution of herbivorous insects.

The compact mitochondrial genome of Zorotypus medoensis provides insights into phylogenetic position of Zoraptera

BACKGROUND

Zoraptera, generally regarded as a member of Polyneoptera, represents one of the most enigmatic insect orders. Although phylogenetic analyses based on a wide array of morphological and/or nuclear data have been performed, the position of Zoraptera is still under debate. Mitochondrial genome (mitogenome) information is commonly considered to be preferable to reconstruct phylogenetic relationships, but no efforts have been made to incorporate it in Zorapteran phylogeny. To characterize Zoraptera mitogenome features and provide insights into its phylogenetic placement, here we sequenced, for the first time, one complete mitogenome of Zoraptera and reconstructed the phylogeny of Polyneoptera.

RESULTS

The mitogenome of Zorotypus medoensis with an A+T content of 72.50% is composed of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a noncoding A+T-rich region. The gene content and arrangement are identical to those considered ancestral for insects. This mitogenome shows a number of very unusual features. First, it is very compact, comprising 14,572 bp, and is the smallest among all known polyneopteran mitogenomes. Second, both noncoding sequences and coding genes exhibit a significant decrease in size compared with those of other polyneopterans. Third, Z. medoensis mitogenome has experienced an accelerated substitution rate. Fourth, truncated secondary structures of tRNA genes occur with loss of dihydrouridine (DHU) arm in trnC, trnR, and trnS(AGN) and loss of TΨC arm in trnH and trnT. The phylogenetic analyses based on the mitogenome sequence information indicate that Zoraptera, represented by Z. medoensis, is recovered as sister to Embioptera. However, both Zoraptera and Embioptera exhibit very long branches in phylogenetic trees.

CONCLUSIONS

Characterization of Z. medoensis mitogenome contributes to our understanding of the enigmatic Zoraptera. Mitogenome data demonstrate an overall strong resolution of deep-level phylogenies of Polyneoptera but not Insecta. It is preferable to expand taxon sampling of Zoraptera and other poorly represented orders in future to break up long branches.

Three-dimensional reconstruction and the phylogeny of extinct chelicerate orders

Arachnids are an important group of arthropods. They are: diverse and abundant; a major constituent of many terrestrial ecosystems; and possess a deep and extensive fossil record. In recent years a number of exceptionally preserved arachnid fossils have been investigated using tomography and associated techniques, providing valuable insights into their morphology. Here we use X-ray microtomography to reconstruct members of two extinct arachnid orders. In the Haptopoda, we demonstrate the presence of 'clasp-knife' chelicerae, and our novel redescription of a member of the Phalangiotarbida highlights leg details, but fails to resolve chelicerae in the group due to their small size. As a result of these reconstructions, tomographic studies of three-dimensionally preserved fossils now exist for three of the four extinct orders, and for fossil representatives of several extant ones. Such studies constitute a valuable source of high fidelity data for constructing phylogenies. To illustrate this, here we present a cladistic analysis of the chelicerates to accompany these reconstructions. This is based on a previously published matrix, expanded to include fossil taxa and relevant characters, and allows us to: cladistically place the extinct arachnid orders; explicitly test some earlier hypotheses from the literature; and demonstrate that the addition of fossils to phylogenetic analyses can have broad implications. Phylogenies based on chelicerate morphology-in contrast to molecular studies-have achieved elements of consensus in recent years. Our work suggests that these results are not robust to the addition of novel characters or fossil taxa. Hypotheses surrounding chelicerate phylogeny remain in a state of flux.