Hojun Song

Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified

Nuclear mitochondrial pseudogenes (numts) are nonfunctional copies of mtDNA in the nucleus that have been found in major clades of eukaryotic organisms. They can be easily coamplified with orthologous mtDNA by using conserved universal primers; however, this is especially problematic for DNA barcoding, which attempts to characterize all living organisms by using a short fragment of the mitochondrial cytochrome c oxidase I (COI) gene. Here, we study the effect of numts on DNA barcoding based on phylogenetic and barcoding analyses of numt and mtDNA sequences in two divergent lineages of arthropods: grasshoppers and crayfish. Single individuals from both organisms have numts of the COI gene, many of which are highly divergent from orthologous mtDNA sequences, and DNA barcoding analysis incorrectly overestimates the number of unique species based on the standard metric of 3% sequence divergence. Removal of numts based on a careful examination of sequence characteristics, including indels, in-frame stop codons, and nucleotide composition, drastically reduces the incorrect inferences of the number of unique species, but even such rigorous quality control measures fail to identify certain numts. We also show that the distribution of numts is lineage-specific and the presence of numts cannot be known a priori. Whereas DNA barcoding strives for rapid and inexpensive generation of molecular species tags, we demonstrate that the presence of COI numts makes this goal difficult to achieve when numts are prevalent and can introduce serious ambiguity into DNA barcoding.

A mitochondrial genome phylogeny of the Neuropterida (lace-wings, alderflies and snakeflies) and their relationship to the other holometabolous insect orders

We present a mitochondrial (mt) genome phylogeny inferring relationships within Neuropterida (lacewings, alderflies and camel flies) and between Neuropterida and other holometabolous insect orders. Whole mt genomes were sequenced for Sialis hamata (Megaloptera: Sialidae), Ditaxis latistyla (Neuroptera: Mantispidae), Mongoloraphidia harmandi (Raphidioptera: Raphidiidae), Macrogyrus oblongus (Coleoptera: Gyrinidae), Rhopaea magnicornis (Coleoptera: Scarabaeidae), and Mordella atrata (Coleoptera: Mordellidae) and compared against representatives of other holometabolous orders in phylogenetic analyses. Additionally, we test the sensitivity of phylogenetic inferences to four analytical approaches: inclusion vs. exclusion of RNA genes, manual vs. algorithmic alignments, arbitrary vs. algorithmic approaches to excluding variable gene regions and how each approach interacts with phylogenetic inference methods (parsimony vs. Bayesian inference). Of these factors, phylogenetic inference method had the most influence on interordinal relationships. Bayesian analyses inferred topologies largely congruent with morphologically‐based hypotheses of neuropterid relationships, a monophyletic Neuropterida whose sister group is Coleoptera. In contrast, parsimony analyses failed to support a monophyletic Neuropterida as Raphidioptera was the sister group of the entire Holometabola excluding Hymenoptera, and Neuroptera + Megaloptera is the sister group of Diptera, a relationship which has not previously been proposed based on either molecular or morphological data sets. These differences between analytical methods are due to the high among site rate heterogeneity found in insect mt genomes which is properly modelled by Bayesian methods but results in artifactual relationships under parsimony. Properly analysed, the mt genomic data set presented here is among the first molecular data to support traditional, morphology‐based interpretations of relationships between the three neuropterid orders and their grouping with Coleoptera.

300 million years of diversification: elucidating the patterns of orthopteran evolution based on comprehensive taxon and gene sampling

Orthoptera is the most diverse order among the polyneopteran groups and includes familiar insects, such as grasshoppers, crickets, katydids, and their kin. Due to a long history of conflicting classification schemes based on different interpretations of morphological characters, the phylogenetic relationships within Orthoptera are poorly understood and its higher classification has remained unstable. In this study, we establish a robust phylogeny of Orthoptera including 36 of 40 families representing all 15 currently recognized superfamilies and based on complete mitochondrial genomes and four nuclear loci, in order to test previous phylogenetic hypotheses and to provide a framework for a natural classification and a reference for studying the pattern of divergence and diversification. We find strong support for monophyletic suborders (Ensifera and Caelifera) as well as major superfamilies. Our results corroborate most of the higher‐level relationships previously proposed for Caelifera, but suggest some novel relationships for Ensifera. Using fossil calibrations, we provide divergence time estimates for major orthopteran lineages and show that the current diversity has been shaped by dynamic shifts of diversification rates at different geological times across different lineages. We also show that mitochondrial tRNA gene orders have been relatively stable throughout the evolutionary history of Orthoptera, but a major tRNA gene rearrangement occurred in the common ancestor of Tetrigoidea and Acridomorpha, thereby representing a robust molecular synapomorphy, which has persisted for 250 Myr.

A century of paraphyly: A molecular phylogeny of katydids (Orthoptera: Tettigoniidae) supports multiple origins of leaf-like wings

The phylogenetic relationships of Tettigoniidae (katydids and bush-crickets) were inferred using molecular sequence data. Six genes (18S rDNA, 28S rDNA, Cytochrome Oxidase II, Histone 3, Tubulin Alpha I, and Wingless) were sequenced for 135 ingroup taxa representing 16 of the 19 extant katydid subfamilies. Five subfamilies (Tettigoniinae, Pseudophyllinae, Mecopodinae, Meconematinae, and Listroscelidinae) were found to be paraphyletic under various tree reconstruction methods (Maximum Likelihood, Bayesisan Inference and Maximum Parsimony). Seven subfamilies - Conocephalinae, Hetrodinae, Hexacentrinae, Saginae, Phaneropterinae, Phyllophorinae, and Lipotactinae - were each recovered as well-supported monophyletic groups. We mapped the small and exposed thoracic auditory spiracle (a defining character of the subfamily Pseudophyllinae) and found it to be homoplasious. We also found the leaf-like wings of katydids have been derived independently in at least six lineages.

Mitochondrial genomics in Orthoptera using MOSAS

We present complete mitochondrial genomes (mitogenomes) for three orthopterans (Xyleus modestus, Physemacris variolosa, and Ellipes minuta) and describe MOSAS (manipulation, organization, storage, and analysis of sequences), software we developed to facilitate annotation and analysis. We analyze the base composition, start and stop codons, non-coding regions, and gene order among these and 18 other orthopteran mitogenomes from GenBank and reconstruct a phylogeny of Orthoptera. We propose a tetranucleotide start codon for cox1, and hypothesize that the tRNAAsp–tRNALys rearrangement is a synapomorphy for Acridomorpha, but not Caelifera. We further describe MOSAS, user-friendly software we used for this analysis. MOSAS streamlines sequence data storage, organization, annotation, and alignment, and provides convenient search tools for dataset construction and a robust annotation engine particularly suited to annotating mitogenomes (available at http://mosas.byu.edu).

Phylogeny of bird‐grasshopper subfamily Cyrtacanthacridinae (Orthoptera: Acrididae) and the evolution of locust phase polyphenism

Locust phase polyphenism is an extreme form of density‐dependent phenotypic plasticity in which solitary and cryptic grasshoppers can transform into gregarious and conspicuous locusts in response to an increase in local population density. We investigated the evolution of this complex phenotypic plasticity in a phylogenetic framework using a morphological phylogeny of Cyrtacanthacridinae, which contains some of the most important locust species, and a comprehensive literature review on the biology and ecology of all known members of the subfamily. A phylogenetic analysis based on 71 morphological characters yielded a well‐resolved tree and found that locust phase polyphenism evolved multiple times within the subfamily. The literature review demonstrated that many cyrtacanthacridine species, both locust and sedentary, are capable of expressing density‐dependent color plasticity. When this color plasticity was divided into two smaller components, background coloration and development of black pigmentation, and when these plastic traits were optimized on to the phylogeny, we found that the physiological mechanisms underlying this plasticity were plesiomorphic for the subfamily. We also found that different locust species in Cyrtacanthacridinae express both similarities and differences in their locust phase polyphenism. Because locust phase polyphenism is a complex syndrome consisting of numerous plastic traits, we treat it as a composite character and dissected it into smaller components. The similarities among locust species could be attributed to shared ancestry and the differences could be attributed to the certain components of locust phase polyphenism evolving at different rates.

Density-Dependent Phase Polyphenism in Nonmodel Locusts: A Minireview

Although the specific mechanisms of locust phase transformation are wellunderstood for model locust species such as the desert locust Schistocerca gregaria and the migratory locust Locusta migratoria, the expressions of density-dependent phase polyphenism in other nonmodel locust species are not wellknown. The present paper is an attempt to review and synthesize what we know about these nonmodel locusts. Based on all available data, I find that locust phase polyphenism is expressed in many different ways in different locust species and identify a pattern that locust species often belong to large taxonomic groups which contain mostly nonswarming grasshopper species. Although locust phase polyphenism has evolved multiple times within Acrididae, I argue that its evolution should be studied from a phylogenetic perspective because I find similar density-dependent phenotypic plasticity among closely related species. Finally, I emphasize the importance of comparative analyses in understanding the evolution of locust phase and propose a phylogeny-based research framework.

Rampant nuclear insertion of mtDNA across diverse lineages within Orthoptera (Insecta).

Nuclear mitochondrial pseudogenes (numts) are non-functional fragments of mtDNA inserted into the nuclear genome. Numts are prevalent across eukaryotes and a positive correlation is known to exist between the number of numts and the genome size. Most numt surveys have relied on model organisms with fully sequenced nuclear genomes, but such analyses have limited utilities for making a generalization about the patterns of numt accumulation for any given clade. Among insects, the order Orthoptera is known to have the largest nuclear genome and it is also reported to include several species with a large number of numts. In this study, we use Orthoptera as a case study to document the diversity and abundance of numts by generating numts of three mitochondrial loci across 28 orthopteran families, representing the phylogenetic diversity of the order. We discover that numts are rampant in all lineages, but there is no discernable and consistent pattern of numt accumulation among different lineages. Likewise, we do not find any evidence that a certain mitochondrial gene is more prone to nuclear insertion than others. We also find that numt insertion must have occurred continuously and frequently throughout the diversification of Orthoptera. Although most numts are the result of recent nuclear insertion, we find evidence of very ancient numt insertion shared by highly divergent families dating back to the Jurassic period. Finally, we discuss several factors contributing to the extreme prevalence of numts in Orthoptera and highlight the importance of exploring the utility of numts in evolutionary studies.

Molecular phylogeny of the Kosciuscola grasshoppers endemic to the Australian alpine and montane regions

Abstract. Diversity and speciation in Australia’s alpine biota are poorly understood. Here we present a molecular phylogeny of the Australian alpine grasshopper genus Kosciuscola (Sjösted) that currently includes five described species. These grasshoppers are of interest not only because of their alpine distribution but also for the extraordinary colour change exhibited by the species K. tristis, whose males turn turquoise when their body temperature exceeds 25°C. We reconstructed the phylogeny with two fragments of the mitochondrial genome using parsimony, maximum likelihood and Bayesian analyses and our data support the current taxonomy. Further, our data show little geographic structuring within some clades, which is puzzling since members of Kosciuscola are brachypterous. Finally, our data coupled with our observations on colouration provide evidence for a genetically distinct clade of K. tristis in the Victorian Alps. This is among the first molecular studies of an alpine invertebrate and one of a few on non-endangered, widespread Australian alpine species. More phylogenetic studies in the Australian Alps are required if we are to understand the evolution of alpine fauna and establish baseline data to monitor their response to climate change.

Non-swarming grasshoppers exhibit density-dependent phenotypic plasticity reminiscent of swarming locusts.

Locusts are well known for exhibiting an extreme form of density-dependent phenotypic plasticity known as locust phase polyphenism. At low density, locust nymphs are cryptically colored and shy, but at high density they transform into conspicuously colored and gregarious individuals. Most of what we know about locust phase polyphenism come from the study of the desert locust Schistocerca gregaria (Forskål), which is a devastating pest species affecting many countries in North Africa and the Middle East. The desert locust belongs to the grasshopper genus Schistocerca Stål, which includes mostly non-swarming, sedentary species. Recent phylogenetic studies suggest that the desert locust is the earliest branching lineage within Schistocerca, which raises a possibility that the presence of density-dependent phenotypic plasticity may be a plesiomorphic trait for the whole genus. In order to test this idea, we have quantified the effect of rearing density in terms of the resulting behavior, color, and morphology in two non-swarming Schistocerca species native to Florida. When reared in both isolated and crowded conditions, the two non-swarming species, Schistocerca americana (Drury) and Schistocerca serialis cubense (Saussure) clearly exhibited plastic reaction norms in all traits measured, which were reminiscent of the desert locust. Specifically, we found that both species were more active and more attracted to each other when reared in a crowded condition than in isolation. They were mainly bright green in color when isolated, but developed strong black patterns and conspicuous background colors when crowded. We found a strong effect of rearing density in terms of size. There were also more mechanoreceptor hairs on the outer face of the hind femora in the crowded nymphs in both species. Although both species responded similarly, there were some clear species-specific differences in terms of color and behavior. Furthermore, we compare and contrast our findings with those on the desert locust and other relevant studies. We attribute the presence of density-dependent phenotypic plasticity in the non-swarming Schistocerca species to phylogenetic conservatism, but there may be a possible role of local adaptation in further shaping the ultimate expressions of plasticity.