Li Wei Wu

The Complete Mitochondrial Genome of the Near-Threatened Swallowtail, Agehana maraho (Lepidoptera: Papilionidae): Evaluating Sequence Variability and Suitable Markers for Conservation Genetic Studies

ABSTRACT: Agehana maraho (Shiraki and Sonan, 1934) is a near-threatened swallowtail butterfly endemic to Taiwan. As a first step in evaluating the most variable molecular markers for further population genetic and conservation studies of this and other insects, the entire mitochondrial genome (mitogenome) was sequenced (16,094bp). The most distinctive structure of the Agehana mitogenome is the control region (CR; 1,270bp). This is the longest CR found so far in any lepidopteran, and it also represents the first known case of two units of macro repeats within a tandem region. In a comparison with another 12 lepidopteran mitogenomes, the genes atp8, atp6, and nad6 were found to be more variable than cox1, suggesting an undue focus on cox1 (COI) in identification and phylogeographic studies. A combination of these first three genes plus the CR, comprising micro as well as macro repeats, may thus provide more suitable markers for conservation genetic studies, not only of this near-threatened species, but also of many other insects.

Ancient Neotropical origin and recent recolonisation: Phylogeny, biogeography and diversification of the Riodinidae (Lepidoptera: Papilionoidea)

We present the first dated higher-level phylogenetic and biogeographic analysis of the butterfly family Riodinidae. This family is distributed worldwide, but more than 90% of the c. 1500 species are found in the Neotropics, while the c. 120 Old World species are concentrated in the Southeast Asian tropics, with minor Afrotropical and Australasian tropical radiations, and few temperate species. Morphologically based higher classification is partly unresolved, with genera incompletely assigned to tribes. Using 3666bp from one mitochondrial and four nuclear markers for each of 23 outgroups and 178 riodinid taxa representing all subfamilies, tribes and subtribes, and 98 out of 145 described genera of riodinids, we estimate that Riodinidae split from Lycaenidae about 96Mya in the mid-Cretaceous and started to diversify about 81Mya. The Riodinidae are monophyletic and originated in the Neotropics, most likely in lowland proto-Amazonia. Neither the subfamily Euselasiinae nor the Nemeobiinae are monophyletic as currently constituted. The enigmatic, monotypic Neotropical genera Styx and Corrachia (most recently treated in Euselasiinae: Corrachiini) are highly supported as derived taxa in the Old World Nemeobiinae, with dispersal most likely occurring across the Beringia land bridge during the Oligocene. Styx and Corrachia, together with all other nemeobiines, are the only exclusively Primulaceae-feeding riodinids. The steadily increasing proliferation of the Neotropical Riodininae subfamily contrasts with the decrease in diversification in the Old World, and may provide insights into factors influencing the diversification rate of this relatively ancient clade of Neotropical insects.

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

BackgroundMitogenomic 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.ResultsWe 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.ConclusionsThirty 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.

Phylogeny and Historical Biogeography of Asian Pterourus Butterflies (Lepidoptera: Papilionidae): A Case of Intercontinental Dispersal from North America to East Asia

The phylogenetic status of the well-known Asian butterflies often known as Agehana (a species group, often treated as a genus or a subgenus, within Papilio sensu lato) has long remained unresolved. Only two species are included, and one of them especially, Papilio maraho, is not only rare but near-threatened, being monophagous on its vulnerable hostplant, Sassafras randaiense (Lauraceae). Although the natural history and population conservation of “Agehana” has received much attention, the biogeographic origin of this group still remains enigmatic. To clarify these two questions, a total of 86 species representatives within Papilionidae were sampled, and four genes (concatenated length 3842 bp) were used to reconstruct their phylogenetic relationships and historical scenarios. Surprisingly, “Agehana” fell within the American Papilio subgenus Pterourus and not as previously suggested, phylogenetically close to the Asian Papilio subgenus Chilasa. We therefore formally synonymize Agehana with Pterourus. Dating and biogeographic analysis allow us to infer an intercontinental dispersal of an American ancestor of Asian Pterourus in the early Miocene, which was coincident with historical paleo-land bridge connections, resulting in the present “East Asia-America” disjunction distribution. We emphasize that species exchange between East Asia and America seems to be a quite frequent occurrence in butterflies during the Oligocene to Miocene climatic optima.

Microbial Community Composition and Functional Capacity in a Terrestrial Ferruginous, Sulfate-Depleted Mud Volcano

Terrestrial mud volcanoes (MVs) are an important natural source of methane emission. The role of microbial processes in methane cycling and organic transformation in such environments remains largely unexplored. In this study, we aim to uncover functional potentials and community assemblages across geochemical transitions in a ferruginous, sulfate-depleted MV of eastern Taiwan. Geochemical profiles combined with 16S rRNA gene abundances indicated that anaerobic oxidation of methane (AOM) mediated by ANME-2a group coincided with iron/manganese reduction by Desulfuromonadales at shallow depths deprived of sulfate. The activity of AOM was stimulated either by methane alone or by methane and a range of electron acceptors, such as sulfate, ferrihydrite, and artificial humic acid. Metagenomic analyses revealed that functional genes for AOM and metal reduction were more abundant at shallow intervals. In particular, genes encoding pili expression and electron transport through multi-heme cytochromes were prevalent, suggesting potential intercellular interactions for electron transport involved in AOM. For comparison, genes responsible for methanogenesis and degradation of chitin and plant-derived molecules were more abundant at depth. The gene distribution combined with the enhanced proportions of 16S rRNA genes related to methanogens and heterotrophs, and geochemical characteristics suggest that particulate organic matter was degraded into various organic entities that could further fuel in situ methanogenesis. Finally, genes responsible for aerobic methane oxidation were more abundant in the bubbling pool and near-surface sediments. These methane oxidizers account for the ultimate attenuation of methane discharge into the atmosphere. Overall, our results demonstrated that various community members were compartmentalized into stratified niches along geochemical gradients. These community members form a metabolic network that cascades the carbon transformation from the upstream degradation of recalcitrant organic carbon with fermentative production of labile organic entities and methane to downstream methane oxidation and metal reduction near the surface. Such a metabolic architecture enables effective methane removal under ferruginous, sulfate-depleted conditions in terrestrial MVs.

Unravelling relationships among the shared stripes of sailors: Mitogenomic phylogeny of Limenitidini butterflies (Lepidoptera, Nymphalidae, Limenitidinae), focusing on the genera Athyma and Limenitis

The phylogenetic relationships of the nymphalid butterfly tribe Limenitidini are best known for the genera Limenitis and Adelpha, model taxa for evolutionary processes such as Batesian mimicry and rapid adaptive radiations. Whereas these American limenitidines have received the most attention, phylogenetic relationships of their Asian relatives are still controversial and largely unexplored. Even one of the largest genera in Asia, Athyma, is polyphyletic. To clarify the phylogenetic relationships of these Asian Limenitidini, a total of 53 representatives were sampled; 37 have their mitogenomes sequenced for the first time. Our phylogenetic results confirm that mitogenomic data provides well-resolved relationships at most major levels of the phylogeny, even using different partition schemes or different inference methods. Interestingly, our results show that some Athyma taxa are embedded within the genus Limenitis, whereas the genus Tacola, previously considered to be a synonym of Athyma, needs to be recognized as a valid clade. Additionally, the other Limenitidini genera in Asia (namely Tarattia, Litinga, Sumalia, Pandita and Patsuia) are now grouped either within Athyma or Limenitis, so these genera need to be sunk. Importantly, we also show that the mainly Old World Limenitis and entirely New World Adelpha are sister groups, confirming the relevance of Asian lineages to global studies of Limenitis evolution.

A distinct species, Dodona formosana, detected in the dodona eugenes species complex: Clarification of the taxonomic status of the Punch butterfly in Taiwan

Abstract The Tailed Punch, Dodona eugenes, is widely distributed in East Asia with seven subspecies currently recognized. However, two of them, namely ssp. formosana and ssp. esakii found in Taiwan, are hard to distinguish from each other due to ambiguous diagnostic characters. In this study, their taxonomic status is clarified by comparing genitalia characters and phylogenetic relationships based on mitochondrial sequences, COI and COII (total 2211 bps). Our results show that there is no reliable feature to separate these two subspecies. Surprisingly we found that Dodona in Taiwan is more closely related to the Orange Punch, D. egeon, than to other subspecies of D. eugenes. Therefore, the following nomenclatural changes are proposed: Dodona eugenes formosana is revised to specific status as Dodona formosana Matsumura, 1919, stat. rev, and ssp. esakii is sunk to a junior synonym of Dodona formosana syn. n.

Higher DNA insert fragment sizes improve mitogenomic assemblies from metagenomic pyrosequencing datasets: an example using Limenitidinae butterflies (Lepidoptera, Nymphalidae)

Abstract A large number of diverse mitogenomic sequences can be obtained more easily and affordably via mitochondrial metagenomics, which generates high-throughput sequences directly from sheared DNA extractions and assembles mitogenomic sequences using a few bioinformatic processing steps. However, following de novo assembly analysis, the optimal DNA fragment insert size is unclear. In this study, four extracted Limenitidinae butterfly DNA samples were sonically fragmented, and two fragment size ranges (200–400 and 400–600 bp) of each sample were tagged with different barcodes, producing pyrosequencing datasets. The results show that the datasets generated from longer DNA insert fragments result in better coverage and more complete mitogenomic sequences, and the phylogenetic analysis shows high support at nodes, revealing that Athyma butterflies do not represent a monophyletic group. Therefore, we recommend using longer insert DNA fragment sizes to generate high-throughput datasets for obtaining complete mitogenomic sequences which can improve phylogenetic studies.