Higher level phylogeny of Satyrinae butterflies (Lepidoptera: Nymphalidae) based on DNA sequence data

Characterization of the complete mitochondrial genome of an endemic species in China, Aulocera merlina (Lepidoptera: Nymphalidae: Satyrinae) and phylogenetic analysis within Satyrinae

The mitochondrial genome (mitogenome) has been extensively used as molecular markers in determining the insect phylogenetic relationships. In order to resolve the relationships among tribes and subtribes of Satyrinae at the mitochondrial genomic level, we obtained the complete mitogenome of Aulocera merlina (Oberthür, 1890) (Lepidoptera: Nymphalidae: Satyrinae) with a size of 15,259 bp. The mitogenome consisted of 37 typical genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A + T-rich region. The gene organization and arrangement were similar to those of all other known Satyrinae mitogenomes. All PCGs were initiated with the canonical codon pattern ATN, except for the cox1 gene, which used an atypical CGA codon. Nine PCGs used the complete stop codon TAA, while the remaining PCGs (cox1, cox2, nad4, and nad5) were terminated with a single T nucleotide. The canonical cloverleaf secondary structures were found in all tRNAs, except for trnS1 which lacked a dihydrouridine arm. The 448 bp A + T-rich region was located between rrnS and trnM, and it included the motif ATAGA followed by a 19-bp poly-T stretch and a microsatellite-like (TA)6 element preceded by the ATTTA motif. The phylogenetic tree, inferred using Bayesian inference and maximum likelihood methods, generated similar tree topologies, revealing well-supported monophyletic groups at the tribe level and recovering the relationship ((Satyrini + Melanitini) + ((Amathusiini + Elymniini) + Zetherini)). The close relationship between Satyrina and Melanargiina within the Satyrini was widely accepted. Additionally, Lethina, Parargina, and Mycalesina were closely related and collectively formed a sister group to Coenonymphina. Moreover, A. merlina was closely related to Oeneis buddha within the Satyrina. These findings will provide valuable information for future studies aiming to elucidate the phylogenetic relationships of Satyrinae.

Chromosome-Level Assembly and Annotation of the Pearly Heath Coenonympha arcania Butterfly Genome

We present the first chromosome-level genome assembly and annotation of the pearly heath Coenonympha arcania, generated with a PacBio HiFi sequencing approach and complemented with Hi-C data. We additionally compare synteny, gene, and repeat content between C. arcania and other Lepidopteran genomes. This reference genome will enable future population genomics studies with Coenonympha butterflies, a species-rich genus that encompasses some of the most highly endangered butterfly taxa in Europe.

Wing pattern diversity in Eunica butterflies (Nymphalidae: Biblidinae): phylogenetic analysis implies decoupled adaptive trends in dorsal sexual dimorphism and ventral eyespot evolution

Butterfly eyespots are wing patterns reminiscent of vertebrate eyes, formed by concentric rings of contrastingly coloured scales. Eyespots are usually located close to the wing margin and often regarded as the single most conspicuous pattern element of butterfly wing colour displays. Recent efforts to understand the processes involved in the formation of eyespots have been driven mainly by evo-devo approaches focused on model species. However, patterns of change implied by phylogenetic relationships can also inform hypotheses about the underlying developmental mechanisms associated with the formation or disappearance of eyespots, and the limits of phenotypic diversity occurring in nature. Here we present a combined evidence phylogenetic hypothesis for the genus Eunica, a prominent member of diverse Neotropical butterfly communities, that features notable variation among species in eyespot patterns on the ventral hind wing surface. The data matrix consists of one mitochondrial gene region (COI), four nuclear gene regions (GAPDH, RPS5, EF1a and Wingless) and 68 morphological characters. A combined cladistic analysis with all the characters concatenated produced a single most parsimonious tree that, although fully resolved, includes many nodes with modest branch support. The phylogenetic hypothesis presented corroborates a previously proposed morphological trend leading to the loss of eyespots, together with an increase in the size of the conserved eyespots, relative to outgroup taxa. Furthermore, wing colour pattern dimorphism and the presence of androconia suggest that the most remarkable instances of sexual dimorphism are present in the species of Eunica with the most derived eyespot patterns, and are in most cases accompanied by autapomorphic combinations of scent scales and "hair pencils". We discuss natural and sexual selection as potential adaptive explanations for dorsal and ventral wing patterns.

A meta-analysis of butterfly structural colors: their color range, distribution and biological production

Butterfly scales are among the richest natural sources of optical nanostructures, which produce structural color and iridescence. Several recurring nanostructure types have been described, such as ridge multilayers, gyroids and lower lamina thin films. While the optical mechanisms of these nanostructure classes are known, their phylogenetic distributions and functional ranges have not been described in detail. In this Review, we examine a century of research on the biological production of structural colors, including their evolution, development and genetic regulation. We have also created a database of more than 300 optical nanostructures in butterflies and conducted a meta-analysis of the color range, abundance and phylogenetic distribution of each nanostructure class. Butterfly structural colors are ubiquitous in short wavelengths but extremely rare in long wavelengths, especially red. In particular, blue wavelengths (around 450 nm) occur in more clades and are produced by more kinds of nanostructures than other hues. Nanostructure categories differ in prevalence, phylogenetic distribution, color range and brightness. For example, lamina thin films are the least bright; perforated lumen multilayers occur most often but are almost entirely restricted to the family Lycaenidae; and 3D photonic crystals, including gyroids, have the narrowest wavelength range (from about 450 to 550 nm). We discuss the implications of these patterns in terms of nanostructure evolution, physical constraint and relationships to pigmentary color. Finally, we highlight opportunities for future research, such as analyses of subadult and Hesperid structural colors and the identification of genes that directly build the nanostructures, with relevance for biomimetic engineering.

Description of life cycle and immature stages of Yphthimoides celmis (Godart, [1824]) (Lepidoptera: Nymphalidae: Satyrinae): a contribution to the taxonomic status of the genus

This article describes the immature stages of the satyrid butterfly Yphthimoides celmis (Godart) for the first time. Morphology, head and body chaetotaxy, and coloration patterns are compared with those of other Yphthimoides species. Larvae behavioral aspects and oviposition behavior are also documented here. Overwintering larvae are observed for the first time for this genus. The larval stage has four or five instars in non-overwintering while six in overwintering larvae. The average duration of the life cycle is around 82 days for non-overwintering, while it is 227 days for overwintering larvae under laboratory conditions. The increase in instar stages and the development time duration of immature stages in overwintering larvae may be due to climatic conditions in the study area and also a strategy to survive the cold season. Eggs have irregular pentagonal and hexagonal cells. First instar larvae have a dark head capsule and a pattern of reddish stripes along their body, with clubbed body setae. Pupae are short and smooth with a non-uniform brown color, presenting a lighter clypeus and a darker ventral portion in the cremaster. We consider that the morphological information of immature stages can improve further phylogenetic studies.

Mitogenomes of Eight Nymphalidae Butterfly Species and Reconstructed Phylogeny of Nymphalidae (Nymphalidae: Lepidoptera)

The Nymphalidae family of cosmopolitan butterflies (Lepidoptera) comprises approximately 7200 species found on all continents and in all habitats. However, debate persists regarding the phylogenetic relationships within this family. In this study, we assembled and annotated eight mitogenomes of Nymphalidae, constituting the first report of complete mitogenomes for this family. Comparative analysis of 105 mitochondrial genomes revealed that the gene compositions and orders were identical to the ancestral insect mitogenome, except for Callerebia polyphemus trnV being before trnL and Limenitis homeyeri having two trnL genes. The results regarding length variation, AT bias, and codon usage were consistent with previous reports on butterfly mitogenomes. Our analysis indicated that the subfamilies Limenitinae, Nymphalinae, Apaturinae, Satyrinae, Charaxinae, Heliconiinae, and Danainae are monophyletic, while the subfamily the subfamily Cyrestinae is polyphyletic. Danainae is the base of the phylogenetic tree. At the tribe level, Euthaliini in Limenitinae; Melitaeini and Kallimini in Nymphalinae; Pseudergolini in Cyrestinae; Mycalesini, Coenonymphini, Ypthimini, Satyrini, and Melanitini in Satyrinae; and Charaxini in Charaxinae are regarded as monophyletic groups. However, the tribe Lethini in Satyrinae is paraphyletic, while the tribes Limenitini and Neptini in Limenitinae, Nymphalini and Hypolimni in Nymphalinae, and Danaini and Euploeini in Danainae are polyphyletic. This study is the first to report the gene features and phylogenetic relationships of the Nymphalidae family based on mitogenome analysis, providing a foundation for future studies of population genetics and phylogenetic relationships within this family.

Morphology and phylogenetic significance of the pelvic articular region in elasmobranchs (Chondrichthyes)

The morphology of paired fins is commonly overlooked in morphological studies, particularly the pelvic girdle and fins. Consequently, previous phylogenetic studies incorporating morphological data used few skeletal characters from this complex. In this paper, the phylogenetic significance of pelvic articular characters for elasmobranchs is discussed in light of the morphological variation observed in 130 species, the most comprehensive study exploring the morphology of the pelvic girdle done so far. The 10 morphological characters proposed herein for the pelvic articulation were incorporated into a molecular matrix of NADH2 sequences and submitted to an analysis of maximum parsimony employing extended implied weighting. The most stable tree was selected based on the distortion coefficients, SPR distances (subtree pruning and regrafting) and fit values. Some of the striking synapomorphies recovered within elasmobranchs include the presence of an articular surface for the first enlarged pelvic radial supporting Elasmobranchii and the pelvic articular region for the basipterygium extending from the posterolatral margin of the pelvic girdle over its lateral surface in Echinorhinus + Hexanchiformes. Additionally, the proposed characters and their distributions are discussed considering the relationships recovered and also compared with previous morphological and molecular phylogenetic hypotheses.

Biogeographical events, not cospeciation, might be the main drivers in the historical association between Noctiliostrebla species (Streblidae) and their bulldog bat hosts

The genus Noctiliostrebla Wenzel (Diptera: Streblidae) stands out for its high degree of specificity exhibited in relation to bat species of the genus Noctilio Linnaeus and provides an exciting system for understanding the history behind host–parasite associations. Here, we present a phylogeny of Noctiliostrebla based on an analysis of DNA sequences and morphological characters, along with cophylogenetic and biogeographical analyses. Our results strongly support the monophyly of Noctiliostrebla, but with uncertainties within the genus. With a low frequency of cospeciation events explaining the associations between hosts and parasites, cophylogenetic analyses did not show an overall congruence between the host and parasite phylogenies. Indeed, two parallel histories were recovered in the host–parasite associations, which might indicate that niche segregation is determined evolutionarily, facilitating the coexistence of parasites and promoting diversification. Biogeographical analysis showed a strong spatial congruence between disjunct distributions of Noctiliostrebla and major river basins in South America and with areas of higher elevation, which might be associated with the glacial periods throughout the Pliocene and Pleistocene. Overall, our findings suggest an agreement with the expectations of the ‘Stockholm paradigm’ framework, in which biogeographical events and ecological factors act as important components to explain the associations, instead of cospeciation events.

Exploring Cold Hardiness within a Butterfly Clade: Supercooling Ability and Polyol Profiles in European Satyrinae

The cold hardiness of overwintering stages affects the distribution of temperate and cold-zone insects. Studies on Erebia, a species-rich cold-zone butterfly genus, detected unexpected diversity of cold hardiness traits. We expanded our investigation to eight Satyrinae species of seven genera. We assessed Autumn and Winter supercooling points (SCPs) and concentrations of putatively cryoprotective sugars and polyols via gas chromatography–mass spectrometry. Aphantopus hyperantus and Hipparchia semele survived freezing of body fluids; Coenonympha arcania, C. gardetta, and Melanargia galathea died prior to freezing; Maniola jurtina, Chazara briseis, and Minois dryas displayed a mixed response. SCP varied from −22 to −9 °C among species. Total sugar and polyol concentrations (TSPC) varied sixfold (2 to 12 μg × mg−1) and eightfold including the Erebia spp. results. SCP and TSPC did not correlate. Alpine Erebia spp. contained high trehalose, threitol, and erythritol; C. briseis and C. gardetta contained high ribitol and trehalose; lowland species contained high saccharose, maltose, fructose, and sorbitol. SCP, TSPC, and glycerol concentrations were affected by phylogeny. Species of mountains or steppes tend to be freeze-avoidant, overwinter as young larvae, and contain high concentrations of trehalose, while those of mesic environments tend to be freeze-tolerant, overwinter as later instars, and rely on compounds such as maltose, saccharose, and fructose.

Phylogenetic systematics of Yphthimoides Forster, 1964 and related taxa, with notes on the biogeographical history of Yphthimoides species

Species losses are increasing and may have an impact on our understanding of patterns of evolutionary pathways and phylogenetic relationships among the groups being lost. The knowledge of such patterns can contribute to preventing future losses by identifying which lineages have higher or lower diversification rates, thus informing conservation strategies. Recent years have seen a significant growth in studies of butterfly systematics, allowing a better understanding of evolutionary relationships among most groups and revealing significant taxonomic chaos in several groups. One of the latter groups is the nymphalid subtribe Euptychiina (Satyrinae), which has been shown to include a number of non-monophyletic genera based on recent molecular phylogenetic analyses. Among others, these genera include Yphthimoides, which is widespread throughout the Neotropical region but particularly diverse in the southeastern Neotropics, and a pair of related genera, Pharneuptychia Forster, 1964 and Moneuptychia Forster, 1964. Using molecular data, this study scope and aims was to provide a phylogenetic hypothesis that corroborates Yphthimoides as presently conceived being non-monophyletic, a result reinforced by a comparative study of the male genitalic morphology. Our results also show that Pharneuptychia and Moneuptychia, plus a species misplaced elsewhere in the Euptychiina, Euptychoides castrensis (Schaus, 1902), form a well supported clade, and that the latter 'species' is a complex of cryptic species. We therefore propose a number of taxonomic rearrangements in the present work to resolve these issues: Yphthimoides eriphule (A. Butler, 1867) will be moved to a new genus; Y. affinis (A. Butler, 1867), Y. maepius (Godart, [1824]), Y. mimula (Hayward, 1954), Y. neomaenas (Hayward, 1967) and Y. mythra (Weymer, 1911) are being transferred to Malaveria Viloria & Benmesbah, 2021; Pharneuptychia innocentia (Godart, [1824]) will be moved to another genus to be described; and Euptychoides castrensis, Pharneuptychia romanina (Bryk, 1953) and Yphthimoides viviana (Romieux, 1927) are being moved to Moneuptychia. The dating of divergences points to a split between the ancestral lineage of Yphthimoides and its sister group, Carminda Ebert and Dias, inDias 1998, during the last half of the Miocene, around 11.86 Mya, and to the diversification of the Pharneuptychia during the same time 11.35 (± 3.52) Mya. Biogeographic analysis showed that the most recent common ancestor of Yphthimoides started to diversify either in the the Brazilian Cerrado savannas or in a combined area of Cerrado and South Atlantic Forest, with a possible change in the ancestral habitat of Carminda. Furthermore, ancestral character mapping favors a savanna origin hypothesis over a forest origin hypothesis.

Patterns of diversification and phylogenetic structure in the dorsolateral head musculature of Neotropical electric eels (Ostariophysi: Gymnotiformes), with a myological synonymy

Abstract The present study offers a broad comparative analysis of the dorsolateral head musculature in the Gymnotiformes, with detailed descriptions and illustrations of the dorsolateral head muscles of 83 species representing combined all valid genera. Results permit a detailed assessment of primary homologies and taxonomically-relevant variation across the order. This provides the basis for a myological synonymy, which organizes 33 previously proposed names for 15 recognized muscles. Morphological variation derived from dorsolateral head musculature was coded into 56 characters. When analyzed in isolation, that set of characters results in Gymnotidae as the sister group of remaining gymnotiforms, and all other currently recognized families as monophyletic groups. In a second analysis, myological characters were concatenated with other previously proposed characters into a phenotypic matrix. Results of that analysis reveal new myological synapomorphies for nearly all taxonomic categories within Gymnotiformes. A Partitioned Bremer Support (PBS) was used to asses the significance of comparative myology in elucidating phylogenetic relationships. PBS values show strongly non-uniform distributions on the tree, with positive scores skewed towards more inclusive taxa, and negative PBS values concentrated on less inclusive clades. Our results provide background for future studies on biomechanical constraints evolved in the early stages of gymnotiform evolution.

A Contribution Towards Resolving the Systematics of the High-Altitude Tropical Andean Satyrine Genus Altopedaliodes Forster, 1964 (Lepidoptera, Nymphalinae: Satyrinae)

The first phylogenetic hypotheses for the high Andean satyrine butterfly genus Altopedaliodes is proposed based on sequence data from mitochondrial (COI and COII) and nuclear (EF-1α) genes. Four species previously included in the genus were found not to be closely related to the clade containing the type species for Altopedaliodes, and these species are therefore removed and placed in the appropriate genus: Pedaliodes cocytia (C. Felder and R. Felder, 1867) n. comb., Pedaliodes nebris Thieme 1905rev. comb., Neopedaliodes tamaensis (Pyrcz and Viloria 2007) n. comb., and Altopedaliodes similis Henao, Páez and Rodríguez-M., 2017; the last name is here synonymized with Pedaliodes nebris (Thieme 1905) n. syn. A taxonomic revision of the twelve species in Altopedaliodes as newly circumscribed is presented and all the currently recognized taxa are revised. We describe two new species, A. llanganati Padrón, Pyrcz and Willmott n. sp. (Ecuador, Napo) and A. pilimbala Pyrcz and Boyer n. sp. (Colombia, Cauca), and three new subspecies, A. reissi papallacta Padrón, Pyrcz and Willmott n. ssp. (Ecuador, Napo), A. reissi dominica Pyrcz and Padrón n. ssp. (Colombia, Valle del Cauca) and A. halli cagnoni Boyer, Pyrcz and Padrón n. ssp. (Ecuador, Morona-Santiago). A new status is proposed for A. nucea Pyrcz and Viloria 1999, n. stat., we combine Altopedaliodes scydmaena (Hayward, 1968) n. comb. and Altopedaliodes belmira (Pyrcz and Rodríguez, 2004), n. comb., n. stat. into the genus and treat the latter taxon as a species, and we newly treat Altopedaliodes scydmaena kruegeri Pyrcz 1999n. stat. as a subspecies.

Miocene Climate and Habitat Change Drove Diversification in Bicyclus, Africa's Largest Radiation of Satyrine Butterflies

Compared to other regions, the drivers of diversification in Africa are poorly understood. We studied a radiation of insects with over 100 species occurring in a wide range of habitats across the Afrotropics to investigate the fundamental evolutionary processes and geological events that generate and maintain patterns of species richness on the continent. By investigating the evolutionary history of Bicyclus butterflies within a phylogenetic framework, we inferred the group's origin at the Oligo-Miocene boundary from ancestors in the Congolian rainforests of central Africa. Abrupt climatic fluctuations during the Miocene (ca. 19-17 Ma) likely fragmented ancestral populations, resulting in at least eight early-divergent lineages. Only one of these lineages appears to have diversified during the drastic climate and biome changes of the early Miocene, radiating into the largest group of extant species. The other seven lineages diversified in forest ecosystems during the late Miocene and Pleistocene when climatic conditions were more favourable-warmer and wetter. Our results suggest changing Neogene climate, uplift of eastern African orogens, and biotic interactions might have had different effects on the various subclades of Bicyclus, producing one of the most spectacular butterfly radiations in Africa.

A new species of the enigmatic genus Oressinoma Doubleday from the Andes of central Peru (Lepidoptera: Nymphalidae, Satyrinae)

A new, and only the third known species of the Neotropical montane genus Oressinoma Doubleday is describedO. sorina n. sp., from the Andes of central Peru. It is distinguishable immediately from the other two congeners by the shape of the hindwing underside submarginal orange band, and by the male genitalia. The systematics of Oressinoma are reviewed. A preliminary analysis is carried out based on COI barcode confirming the separate specific status of O. sorina n. sp. in relation to other two congeners. Both barcode and genital morphology data suggest that the widespread O. typhla Doubleday may be a complex of allopatric or, locally parapatric species. The genus Oressinoma is the only neotropical member of the predominantly Australian subtribe Coenonymphina, represented in the entire Holarctic by one genus onlyCoenonympha Hbner, considered as the putative sister-genus of Oressinoma. Their origins and relationships are briefly discussed.

Phylogenetic position and taxonomic rearrangement of Davidina (Lepidoptera: Nymphalidae), an enigmatic butterfly genus new for Europe and America

Davidina, an enigmatic butterfly genus described from China in the 19th century, has for a long time been considered a member of the family Pieridae due to its pierid-like wing pattern. In the 20th century, it was transferred to Satyridae (now subfamily Satyrinae of Nymphalidae) based on analysis of the structure of genitalia and placed next to the species-rich genus Oeneis (subtribe Satyrina), being separated from the latter by supposed differences in wing venation. We have conducted a phylogenetic and taxonomic study of the subtribe Satyrina using analysis of molecular and morphological characters. We show that the genus Oeneis is not monophyletic, and consists of two genetically diverged and morphologically differentiated groups that are not sister-groups (Oeneis s.s. and Protoeneis). We also demonstrate that Davidina is closely related to Protoeneis, but not to Oeneis s.s. To resolve this newly discovered non-monophyly and morphological heterogeneity, several species should be removed from Oeneis and transferred to the genus Davidina. As a consequence, we synonymize the name ProtoeneisGorbunov, 2001 with DavidinaOberthür, 1879. We conclude that Davidina is not a monotypic Chinese endemic genus, as has been previously supposed, but is composed of nine species that have a broad distribution area across the Holarctic region, extending to Europe and America.

A release from developmental bias accelerates morphological diversification in butterfly eyespots

Development can bias the independent evolution of traits sharing ontogenetic pathways, making certain evolutionary changes less likely. The eyespots commonly found on butterfly wings each have concentric rings of differing colors, and these serially repeated pattern elements have been a focus for evo-devo research. In the butterfly family Nymphalidae, eyespots have been shown to function in startling or deflecting predators and to be involved in sexual selection. Previous work on a model species of Mycalesina butterfly, Bicyclus anynana, has provided insights into the developmental control of the size and color composition of individual eyespots. Experimental evolution has also shown that the relative size of a pair of eyespots on the same wing surface is highly flexible, whereas they are resistant to diverging in color composition, presumably due to the underlying shared developmental process. This fixed color composition has been considered as a prime example of developmental bias with significant consequences for wing pattern evolution. Here, we test this proposal by surveying eyespots across the whole subtribe of Mycalesina butterflies and demonstrate that developmental bias shapes evolutionary diversification except in the genus Heteropsis which has gained independent control of eyespot color composition. Experimental manipulations of pupal wings reveal that the bias has been released through a novel regional response of the wing tissue to a conserved patterning signal. Our study demonstrates that development can bias the evolutionary independence of traits, but it also shows how bias can be released through developmental innovations, thus, allowing rapid morphological change, facilitating evolutionary diversification.

Fourteen complete mitochondrial genomes of butterflies from the genus Lethe (Lepidoptera, Nymphalidae, Satyrinae) with mitogenome-based phylogenetic analysis

The mitochondrial genome (mitogenome) can help us understand the phylogenetic relationships within the genus Lethe and the subfamily Satyrinae. In this study, we sequenced the complete mitogenomes of 14 Lethe species, which range in size from 15,225 to 15,271 bp, with both 37 genes (13 PCGs, 22 tRNAs, 2 rRNAs) and a noncoding A + T-rich region. The gene arrangement and orientation is similar to typical mitogenomes of Lepidoptera. The Ka/Ks ratio shows that cox1 has the slowest evolutionary rate. The secondary structure of trnN lacks the Pseudouracil loop (TψC loop) in most Lethe species. The inferred phylogenetic analyses show that Lethe is a well-supported monophyletic group, and reveal 2 major clades within the genus Lethe, which is consistent with previous morphological classifications.

A complete time-calibrated multi-gene phylogeny of the European butterflies

With the aim of supporting ecological analyses in butterflies, the third most species-rich superfamily of Lepidoptera, this paper presents the first time-calibrated phylogeny of all 496 extant butterfly species in Europe, including 18 very localised endemics for which no public DNA sequences had been available previously. It is based on a concatenated alignment of the mitochondrial gene COI and up to eleven nuclear gene fragments, using Bayesian inferences of phylogeny. To avoid analytical biases that could result from our region-focussed sampling, our European tree was grafted upon a global genus-level backbone butterfly phylogeny for analyses. In addition to a consensus tree, the posterior distribution of trees and the fully concatenated alignment are provided for future analyses. Altogether a complete phylogenetic framework of European butterflies for use by the ecological and evolutionary communities is presented.

Revisional Notes on the Cloud Forest Butterfly Genus Oxeoschistus Butler in Central America (Lepidoptera: Nymphalidae: Satyrinae)

So far, six species of Oxeoschistus Butler, including its junior synonym Dioriste Thieme, were listed from Central America, with five of them from Costa Rica alone, which appears to represent the highest regional diversity of this Neotropical montane butterfly genus. Our research based on field work, morphological studies and barcode analysis proved that one record is a misunderstanding perpetuated in scientific literature for over a century: Oxeoschistus cothonides Grose-Smith is identified here as an individual form of the female of O. cothon Salvin. The presence of Oxeoschistus tauropolis (Westwood) in Costa Rica, subject to some controversy, is confirmed, and a new local subspecies is described from Costa Rica, O. tauropolis mitsuko Pyrcz & Nakamura n. ssp. Specific status of O. euriphyle Butler is reinstated based on morphological and molecular data. A new subspecies O. hilara lempira Pyrcz n. ssp. is described from Honduras. O. puerta submaculatus Butler is reported for the first time from the Darién region on the Panama-Colombia border. Species relationships are preliminarily evaluated based on COI data concluding, among others, that O. hilara and O. euriphyle are less closely related than previously assumed. Altitudinal and distributional data are revised, and ecological and behavioural information of all the species of Central American Oxeoschistus is provided.

Considerations on the Taxonomy of the Genus Arhuaco Adams and Bernard 1977, and its Relationships with the Genus Pronophila Doubleday [1849] (Nymphalidae, Satyrinae)

Arhuaco Adams & Bernard (1977) is one of the least known genera of Neotropical Satyrinae. It comprises two species and presents an unusual disjunct distribution, with A. ica Adams & Bernard (1977), endemic to the isolated Colombian Sierra Nevada de Santa Marta, and A. dryadina (Schaus 1913) found in the mountains of Costa Rica and Panama. Here, the female of A. dryadina is described, and a new generic diagnosis is presented. Affinities with other genera of the subtribe Pronophilina, in particular the potential closest relatives, such as Pronophila Doubleday (1849), are investigated based on morphological, molecular, ecological, and behavioral data. Results from molecular and morphological sources are incongruent. Molecular data indicate that Arhuaco is paraphyletic, with A. dryadina segregating within the Pronophila clade. Morphological data, by contrast, indicate a closer affinity between the two species currently placed in Arhuaco, favoring the monophyly of the genus, and show no consistent synapomorphies for Arhuaco + Pronophila. A vicariance biogeographical scenario is evaluated.

Uncovered Diversity of a Predominantly Andean Butterfly Clade in the Brazilian Atlantic Forest: a Revision of the Genus Praepedaliodes Forster (Lepidoptera: Nymphalidae, Satyrinae, Satyrini)

The genus Praepedaliodes Forster, 1964, the only representative of the mega-diverse mostly Andean Pedaliodes complex lineage in the Brazilian Atlantic Forest, is revised. Prior to this study, four species were known, P. phanias (Hewitson, 1862), P. granulata (Butler, 1868), P. amussis (Thieme, 1905) and P. exul (Thieme, 1905). Here, a further six are described, all from SE Brazil, expanding to 10 the number of species in this genus. Lectotypes are designated for P. phanias, P. granulata and P. amussis. The genus is most diverse in the Serra da Mantiqueira (São Paulo, Rio de Janeiro, Minas Gerais) and in the Serra Geral (Paraná, Santa Catarina) with seven species occurring in both ranges. Praepedaliodes phanias is the most widespread species and the only one found in the western part of the Atlantic Forest; only this species and P. duartei Dias, Dolibaina & Pyrcz n. sp. occurring to near sea level. Other species, P. zaccae Dolibaina, Dias & Pyrcz n. sp., P. francinii Freitas & Pyrcz n. sp., P. sequeirae Pyrcz, Dias & Dolbaina n. sp., P. landryi Pyrcz & Freitas n. sp. and P. pawlaki Pyrcz & Boyer n. sp. are strictly montane and the highest species richness is reached at 1400-1800 m. One species, P. sequeirae n. sp., is a narrow endemic found only at timberline in the Agulhas Negras massif above 2300 m. Immature stages are described for two species, P. phanias and P. landryi n. sp. Molecular data (barcode region of cytochrome oxidase, subunit I) and adult morphology, including male and female genitalia, support the genus as monophyletic, belonging to a predominantly Andean clade of the Pedaliodes Butler, 1867 complex. Morphological evidences, in particular female genitalia comparative analysis, indicate the genera Physcopedaliodes Forster, 1964 and Panyapedaliodes Forster, 1964 as possibly the closest relatives to Praepedaliodes. Molecular data are inconclusive in this respect.

Molecular systematics of the subfamily Limenitidinae (Lepidoptera: Nymphalidae)

We studied the systematics of the subfamily Limenitidinae (Lepidoptera: Nymphalidae) using molecular methods to reconstruct a robust phylogenetic hypothesis. The molecular data matrix comprised 205 Limenitidinae species, four outgroups, and 11,327 aligned nucleotide sites using up to 18 genes per species of which seven genes (CycY, Exp1, Nex9, PolII, ProSup, PSb and UDPG6DH) have not previously been used in phylogenetic studies. We recovered the monophyly of the subfamily Limenitidinae and seven higher clades corresponding to four traditional tribes Parthenini, Adoliadini, Neptini, Limenitidini as well as three additional independent lineages. One contains the genera Harma + Cymothoe and likely a third, Bhagadatta, and the other two independent lineages lead to Pseudoneptis and to Pseudacraea. These independent lineages are circumscribed as new tribes. Parthenini was recovered as sister to rest of Limenitidinae, but the relationships of the remaining six lineages were ambiguous. A number of genera were found to be non-monophyletic, with Pantoporia, Euthalia, Athyma, and Parasarpa being polyphyletic, whereas Limenitis, Neptis, Bebearia, Euryphura, and Adelpha were paraphyletic.

Eocene lantern fruits from Gondwanan Patagonia and the early origins of Solanaceae

The nightshade family Solanaceae holds exceptional economic and cultural importance. The early diversification of Solanaceae is thought to have occurred in South America during its separation from Gondwana, but the family's sparse fossil record provides few insights. We report 52.2-million-year-old lantern fruits from terminal-Gondwanan Patagonia, featuring highly inflated, five-lobed calyces, as a newly identified species of the derived, diverse New World genus Physalis (e.g., groundcherries and tomatillos). The fossils are considerably older than corresponding molecular divergence dates and demonstrate an ancient history for the inflated calyx syndrome. The derived position of these early Eocene fossils shows that Solanaceae were well diversified long before final Gondwanan breakup.

An illustrated checklist of the genus Elymnias Hübner, 1818 (Nymphalidae, Satyrinae)

We review the genus Elymnias Hübner, 1818, a morphologically diverse satyrine butterfly clade involved in multifarious Batesian mimicry relationships throughout Asia and Africa. A variety of different model species are mimicked, and many Elymnias species are sexually dimorphic mimics, with males and females resembling different model species. We revise species and subspecies delimitations in light of an integrative taxonomic investigation using external morphology, male and female genital morphology, and a multi-locus molecular phylogeny. There is little interspecific genitalic variation among species in this group, and previous taxonomists therefore relied almost entirely on wing patterns. Our molecular phylogenetic analysis reveals several examples of polymorphism or wing pattern divergence within a single species currently classified as two or more different species. We also found examples of wing pattern convergence among disparate lineages that mimic the same widespread model species. Frequently, two or more phenotypically similar species were classified as a single species. This comprehensive checklist reviews all names associated with Elymnias to align its taxonomy with the evolutionary history of the group. All available information on nomenclature, type localities, repositories of type specimens, and geographical distributions is summarized, and images of adult specimens and genitalia are provided along with distribution maps of all species and selected subspecies. We identify 2 species incertae sedis, establish 15 monophyletic species groups (including 1 species unplaced in any species group), and make 49 taxonomic changes, including 35 new synonyms, 7 new combinations (2 of which have new status), 1 resurrected combination, 1 resurrected subspecies, and 7 status changes.

DNA barcodes and citizen science provoke a diversity reappraisal for the “ring” butterflies of Peninsular Malaysia (Ypthima: Satyrinae: Nymphalidae: Lepidoptera)

The “rings” belonging to the genus Ypthima are amongst the most common butterflies in Peninsular Malaysia. However, the species can be difficult to tell apart, with keys relying on minor and often non-discrete ring characters found on the hindwing. Seven species have been reported from Peninsular Malaysia, but this is thought to be an underestimate of diversity. DNA barcodes of 165 individuals, and wing and genital morphology, were examined to reappraise species diversity of this genus in Peninsular Malaysia. DNA barcodes collected during citizen science projects—School Butterfly Project and Peninsular Malaysia Butterfly Count—recently conducted in Peninsular Malaysia were included. The new DNA barcodes formed six groups with different Barcode Index Numbers (BINs) representing four species reported in Peninsular Malaysia. When combined with public DNA barcodes from the Barcode Of Life Datasystems, several taxonomic issues arose. We consider the taxon Y. newboldi, formerly treated as a subspecies of Y. baldus, as a distinct species. DNA barcodes also supported an earlier suggestion that Y. nebulosa is a synonym under Y. horsfieldii humei. Two BINs of the genus Ypthima comprising DNA barcodes collected during citizen science projects did not correspond to any species previously reported in Peninsular Malaysia.

The Origin of Large-Bodied Shrimp that Dominate Modern Global Aquaculture

Several shrimp species from the clade Penaeidae are farmed industrially for human consumption, and this farming has turned shrimp into the largest seafood commodity in the world. The species that are in demand for farming are an anomaly within their clade because they grow to much larger sizes than other members of Penaeidae. Here we trace the evolutionary history of the anomalous farmed shrimp using combined data phylogenetic analysis of living and fossil species. We show that exquisitely preserved fossils of †Antrimpos speciosus from the Late Jurassic Solnhofen limestone belong to the same clade as the species that dominate modern farming, dating the origin of this clade to at least 145 mya. This finding contradicts a much younger Late Cretaceous age (ca. 95 mya) previously estimated for this clade using molecular clocks. The species in the farmed shrimp clade defy a widespread tendency, by reaching relatively large body sizes despite their warm water lifestyles. Small body sizes have been shown to be physiologically favored in warm aquatic environments because satisfying oxygen demands is difficult for large organisms breathing in warm water. Our analysis shows that large-bodied, farmed shrimp have more gills than their smaller-bodied shallow-water relatives, suggesting that extra gills may have been key to the clade's ability to meet oxygen demands at a large size. Our combined data phylogenetic tree also suggests that, during penaeid evolution, the adoption of mangrove forests as habitats for young shrimp occurred multiple times independently.

Nutrition shapes life-history evolution across species

Nutrition is a key component of life-history theory, yet we know little about how diet quality shapes life-history evolution across species. Here, we test whether quantitative measures of nutrition are linked to life-history evolution across 96 species of butterflies representing over 50 independent diet shifts. We find that butterflies feeding on high nitrogen host plants as larvae are more fecund, but their eggs are smaller relative to their body size. Nitrogen and sodium content of host plants are also both positively related to eye size. Some of these relationships show pronounced lineage-specific effects. Testis size is not related to nutrition. Additionally, the evolutionary timing of diet shifts is not important, suggesting that nutrition affects life histories regardless of the length of time a species has been adapting to its diet. Our results suggest that, at least for some lineages, species with higher nutrient diets can invest in a range of fitness-related traits like fecundity and eye size while allocating less to each egg as offspring have access to a richer diet. These results have important implications for the evolution of life histories in the face of anthropogenic changes in nutrient availability.

Systematics and historical biogeography of the old world butterfly subtribe Mycalesina (Lepidoptera: Nymphalidae: Satyrinae)

BACKGROUND

Butterflies of the subtribe Mycalesina have radiated successfully in almost all habitat types in Africa, Madagascar, the Indian subcontinent, Indo-China and Australasia. Studies aimed at understanding the reasons behind the evolutionary success of this spectacular Old World butterfly radiation have been hampered by the lack of a stable phylogeny for the group. Here, we have reconstructed a robust phylogenetic framework for the subtribe using 10 genes from 195 exemplar taxa.

RESULTS

We recovered seven well supported clades within the subtribe corresponding to the five traditional genera (Lohora, Heteropsis, Hallelesis, Bicyclus, Mycalesis), one as recently revised (Mydosama) and one newly revised genus (Culapa). The phylogenetic relationships of these mycalesine genera have been robustly established for the first time. Within the proposed phylogenetic framework, we estimated the crown age of the subtribe to be 40 Million years ago (Mya) and inferred its ultimate origin to be in Asia. Our results reveal both vicariance and dispersal as factors responsible for the current widespread distribution of the group in the Old World tropics. We inferred that the African continent has been colonized at least twice by Asian mycalesines within the last 26 and 23 Mya. In one possible scenario, an Asian ancestor gave rise to Heteropsis on continental Africa, which later dispersed into Madagascar and most likely back colonised Asia. The second colonization of Africa by Asian ancestors resulted in Hallelesis and Bicyclus on continental Africa, the descendants of which did not colonise other regions but rather diversified only in continental Africa. The genera Lohora and Mydosama are derivatives of ancestors from continental Asia.

CONCLUSION

Our proposed time-calibrated phylogeny now provides a solid framework within which we can implement mechanistic studies aimed at unravelling the ecological and evolutionary processes that culminated in the spectacular radiation of mycalesines in the Old World tropics.

Morphological characters are compatible with mitogenomic data in resolving the phylogeny of nymphalid butterflies (lepidoptera: papilionoidea: nymphalidae)

Nymphalidae is the largest family of butterflies with their phylogenetic relationships not adequately approached to date. The mitochondrial genomes (mitogenomes) of 11 new nymphalid species were reported and a comparative mitogenomic analysis was conducted together with other 22 available nymphalid mitogenomes. A phylogenetic analysis of the 33 species from all 13 currently recognized nymphalid subfamilies was done based on the mitogenomic data set with three Lycaenidae species as the outgroups. The mitogenome comparison showed that the eleven new mitogenomes were similar with those of other butterflies in gene content and order. The reconstructed phylogenetic trees reveal that the nymphalids are made up of five major clades (the nymphaline, heliconiine, satyrine, danaine and libytheine clades), with sister relationship between subfamilies Cyrestinae and Biblidinae, and most likely between subfamilies Morphinae and Satyrinae. This whole mitogenome-based phylogeny is generally congruent with those of former studies based on nuclear-gene and mitogenomic analyses, but differs considerably from the result of morphological cladistic analysis, such as the basal position of Libytheinae in morpho-phylogeny is not confirmed in molecular studies. However, we found that the mitogenomic phylogeny established herein is compatible with selected morphological characters (including developmental and adult morpho-characters).

Diversity dynamics in Nymphalidae butterflies: effect of phylogenetic uncertainty on diversification rate shift estimates

The species rich butterfly family Nymphalidae has been used to study evolutionary interactions between plants and insects. Theories of insect-hostplant dynamics predict accelerated diversification due to key innovations. In evolutionary biology, analysis of maximum credibility trees in the software MEDUSA (modelling evolutionary diversity using stepwise AIC) is a popular method for estimation of shifts in diversification rates. We investigated whether phylogenetic uncertainty can produce different results by extending the method across a random sample of trees from the posterior distribution of a Bayesian run. Using the MultiMEDUSA approach, we found that phylogenetic uncertainty greatly affects diversification rate estimates. Different trees produced diversification rates ranging from high values to almost zero for the same clade, and both significant rate increase and decrease in some clades. Only four out of 18 significant shifts found on the maximum clade credibility tree were consistent across most of the sampled trees. Among these, we found accelerated diversification for Ithomiini butterflies. We used the binary speciation and extinction model (BiSSE) and found that a hostplant shift to Solanaceae is correlated with increased net diversification rates in Ithomiini, congruent with the diffuse cospeciation hypothesis. Our results show that taking phylogenetic uncertainty into account when estimating net diversification rate shifts is of great importance, as very different results can be obtained when using the maximum clade credibility tree and other trees from the posterior distribution.

Phylogenetic utility of ribosomal genes for reconstructing the phylogeny of five Chinese satyrine tribes (Lepidoptera, Nymphalidae)

Satyrinae is one of twelve subfamilies of the butterfly family Nymphalidae, which currently includes nine tribes. However, phylogenetic relationships among them remain largely unresolved, though different researches have been conducted based on both morphological and molecular data. However, ribosomal genes have never been used in tribe level phylogenetic analyses of Satyrinae. In this study we investigate for the first time the phylogenetic relationships among the tribes Elymniini, Amathusiini, Zetherini and Melanitini which are indicated to be a monophyletic group, and the Satyrini, using two ribosomal genes (28s rDNA and 16s rDNA) and four protein-coding genes (EF-1α, COI, COII and Cytb). We mainly aim to assess the phylogenetic informativeness of the ribosomal genes as well as clarify the relationships among different tribes. Our results show the two ribosomal genes generally have the same high phylogenetic informativeness compared with EF-1α; and we infer the 28s rDNA would show better informativeness if the 28s rDNA sequence data for each sampling taxon are obtained in this study. The placement of the monotypic genus Callarge Leech in Zetherini is confirmed for the first time based on molecular evidence. In addition, our maximum likelihood (ML) and Bayesian inference (BI) trees consistently show that the involved Satyrinae including the Amathusiini is monophyletic with high support values. Although the relationships among the five tribes are identical among ML and BI analyses and are mostly strongly-supported in BI analysis, those in ML analysis are lowly- or moderately- supported. Therefore, the relationships among the related five tribes recovered herein need further verification based on more sampling taxa.

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.

Adapting simultaneous analysis phylogenomic techniques to study complex disease gene relationships

The characterization of complex diseases remains a great challenge for biomedical researchers due to the myriad interactions of genetic and environmental factors. Network medicine approaches strive to accommodate these factors holistically. Phylogenomic techniques that can leverage available genomic data may provide an evolutionary perspective that may elucidate knowledge for gene networks of complex diseases and provide another source of information for network medicine approaches. Here, an automated method is presented that leverages publicly available genomic data and phylogenomic techniques, resulting in a gene network. The potential of approach is demonstrated based on a case study of nine genes associated with Alzheimer Disease, a complex neurodegenerative syndrome. The developed technique, which is incorporated into an update to a previously described Perl script called "ASAP," was implemented through a suite of Ruby scripts entitled "ASAP2," first compiles a list of sequence-similarity based orthologues using PSI-BLAST and a recursive NCBI BLAST+ search strategy, then constructs maximum parsimony phylogenetic trees for each set of nucleotide and protein sequences, and calculates phylogenetic metrics (Incongruence Length Difference between orthologue sets, partitioned Bremer support values, combined branch scores, and Robinson-Foulds distance) to provide an empirical assessment of evolutionary conservation within a given genetic network. In addition to the individual phylogenetic metrics, ASAP2 provides results in a way that can be used to generate a gene network that represents evolutionary similarity based on topological similarity (the Robinson-Foulds distance). The results of this study demonstrate the potential for using phylogenomic approaches that enable the study of multiple genes simultaneously to provide insights about potential gene relationships that can be studied within a network medicine framework that may not have been apparent using traditional, single-gene methods. Furthermore, the results provide an initial integrated evolutionary history of an Alzheimer Disease gene network and identify potentially important co-evolutionary clustering that may warrant further investigation.

Tracing the radiation of Maniola (Nymphalidae) butterflies: new insights from phylogeography hint at one single incompletely differentiated species complex

The use of DNA sequence data often leads to the recognition of cryptic species within putatively well-known taxa. The opposite case, detecting less diversity than originally described, has, however, far more rarely been documented. Maniola jurtina, the Meadow Brown butterfly, occurs all over Europe, whereas all other six species in the genus Maniola are restricted to the Mediterranean area. Among them, three are island endemics on Sardinia, Cyprus, and Chios, respectively. Maniola species are almost indistinguishable morphologically, and hybridization seems to occur occasionally. To clarify species boundaries and diversification history of the genus, we reconstructed the phylogeography and phylogeny of all seven species within Maniola analyzing 138 individuals from across its range using mitochondrial and nuclear genetic markers. Examination of variation in mitochondrial and nuclear DNA surprisingly revealed a case of taxonomic "oversplitting". The topology of the recovered phylogenetic tree is not consistent with accepted taxonomy, but rather reveals haplotype clades that are incongruent with nominal species boundaries: instead of seven species, we recognized only two major, yet incompletely segregated, lineages. Our results are consistent with the hypothesis that Maniola originated in Africa. We suggest that one lineage dispersed over the Strait of Gibraltar and the Iberian Peninsula to the west of Europe, while the other lineage spreads eastward through Asia Minor and over the Bosporus to Eastern Europe.

Community structure of skipper butterflies (Lepidoptera, Hesperiidae) along elevational gradients in Brazilian Atlantic forest reflects vegetation type rather than altitude

Species turnover across elevational gradients has matured into an important paradigm of community ecology. Here, we tested whether ecological and phylogenetic structure of skipper butterfly assemblages is more strongly structured according to altitude or vegetation type along three elevation gradients of moderate extent in Serra do Mar, Southern Brazil. Skippers were surveyed along three different mountain transects, and data on altitude and vegetation type of every collection site were recorded. NMDS ordination plots were used to assess community turnover and the influence of phylogenetic distance between species on apparent community patterns. Ordinations based on ecological similarity (Bray-Curtis index) were compared to those based on phylogenetic distance measures (MPD and MNTD) derived from a supertree. In the absence of a well-resolved phylogeny, various branch length transformation methods were applied together with four different null models, aiming to assess if results were confounded by low-resolution trees. Species composition as well as phylogenetic community structure of skipper butterflies were more prominently related to vegetation type instead of altitude per se. Phylogenetic distances reflected spatial community patterns less clearly than species composition, but revealed a more distinct fauna of monocot feeders associated with grassland habitats, implying that historical factors have played a fundamental role in shaping species composition across elevation gradients. Phylogenetic structure of community turned out to be a relevant additional tool which was even superior to identify faunal contrasts between forest and grassland habitats related to deep evolutionary splits. Since endemic skippers tend to occur in grassland habitats in the Serra do Mar, inclusion of phylogenetic diversity may also be important for conservation decisions.

The complete mitochondrial genome of Callerebia suroia (Lepidoptera: Nymphalidae: Satyrinae)

The complete mitochondrial genome (mitogenome) of Callerebia suroia (Lepidoptera: Nymphalidae: Satyrinae) was determined and analyzed in this paper. The circular genome is 15,208 bp long, including 37 typical mitochondrial genes and one non-coding AT-rich region. All protein-coding genes (PCGs) started with ATN, except for COI gene with CGA(R), which is often found in other butterflies; nine PCGs harbor the typical stop codon TAA, whereas COI, COII, ND5 and ND4 end with a single T. All tRNA genes display typical secondary clover-leaf structures, except for tRNA(Ser)(AGN), whose dihydrouridine (DHU) arm is replaced by a simple loop. The lrRNA and srRNA genes are 1,347 bp and 753 bp in length, with their AT contents of 84.4% and 85.4%, respectively. The 417 bp AT-rich region contains non repetitive sequences, but harbor several features common to the lepidopterans, including the motif ATAGA followed by a 19-bp poly-T stretch and a microsatellite-like (TA)8 element preceded by the ATTTA motif.

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.

An eight-gene molecular phylogeny of the Kickxellomycotina, including the first phylogenetic placement of Asellariales

Kickxellomycotina is a recently described subphylum encompassing four zygomycete orders (Asellariales, Dimargaritales, Harpellales, Kickxellales). These fungi are united by the formation of disciform septal pores containing lenticular plugs. Morphological diversification and life history evolution has made the relationships within and among the four orders difficult to resolve on those grounds alone. Here we infer the phylogeny of the Kickxellomycotina based on an eight-gene supermatrix including both ribosomal rDNA (18S, 28S, 5.8S) and protein sequences (MCM7, TSR1, RPB1, RPB2, β-tubulin). The results of this study demonstrate that Kickxellomycotina is monophyletic and related to members of the Zoopagomycotina. Eight unique clades are distinguished in the Kickxellomycotina, including the four defined orders (Asellariales, Dimargaritales, Harpellales, Kickxellales) as well as four genera previously placed within two of these orders (Barbatospora, Orphella, Ramicandelaber, Spiromyces). Dimargaritales and Ramicandelaber are the earliest diverging members of the subphylum, although the relationship between these taxa remains uncertain. The remaining six clades form a monophyletic group, with Barbatospora diverging first. The next split divides the remaining members of the subphylum into two subclades: (i) Asellariales and Harpellales and (ii) Kickxellales, Orphella and Spiromyces. Estimation of ancestral states for four potentially informative morphological and ecological characters reveals that arthropod endosymbiosis might have been an important factor in the early evolution of the Kickxellomycotina.