A highly conserved nuclear gene for low-level phylogenetics: elongation factor-1 alpha recovers morphology-based tree for heliothine moths

Specialization on Ficus Supported by Genetic Divergence and Morphometrics in Sympatric Host-Populations of the Camellia Aphid, Aphis aurantii

Adaptation to different host plants is considered to be an important driver of the divergence and speciation of herbivorous insects. The application of molecular data and integrated taxonomic practices in recent years may contribute to our understanding of population divergence and speciation, especially for herbivorous insects considered to be polyphagous. Aphis aurantii is an important agricultural and forestry pest with a broad range of host plants. In this study, samples of A. aurantii feeding on different host plants in the same geographical area were collected, and their population genetic divergence and morphological difference were analyzed. Phylogenetic analysis and haplotype network analysis based on five genes revealed that the population on Ficus exhibited significantly genetic divergence from populations on other host plants, which was also supported by the statistical analysis based on measurements of 38 morphological characters. Our results suggest that A. aurantii has undergone specialized evolution on Ficus, and the Ficus population may represent a lineage that is experiencing ongoing sympatric speciation.

A single point mutation causes one-way alteration of pheromone receptor function in two Heliothis species

The sex pheromone processing system of moths has been a major focus of research on olfaction and speciation, as it is highly specific and closely related to reproductive isolation. The two noctuid moths Heliothis virescens and Heliothis subflexa have been used as a model for deciphering the mechanisms underlying differentiation in pheromone communication, but no information exist regarding the functions of the pheromone receptors (PRs) of H. subflexa. Here, we functionally characterized all candidate PRs of H. subflexa, and found that only the response profile of OR6 differed between the two species. Through domain swapping and site-directed mutation followed by functional characterization, we identified a critical amino acid in OR6 caused a one-way alteration in specificity. This result suggests HsubOR6 evolved from an ancestral OR6 gene with a HvirOR6-like function and implies that the evolutionary direction of the receptor specificity was from the H. virescens-like pattern to H. subflexa-like pattern.

Mutation Associated with Orange Fruit Color Increases Concentrations of β-Carotene in a Sweet Pepper Variety (Capsicum annuum L.)

Pepper is the second most important vegetable crop in Bulgarian agriculture and has become the subject of extensive breeding programs that frequently employ induced mutagenesis. The success of breeding programs can be enhanced by the efficient and integral application of different biochemical and molecular methods to characterize specific mutant alleles. On the other hand, identifying new cost-effective methods is important under a limited-resources environment. In this paper we compare the levels of five health-related carotenoid compounds of fruits (α-carotene, β-carotene, lutein, β-cryptoxanthin, zeaxanthin) between a mutant variety Oranzheva kapia (possessing high ß-carotene concentration) and a corresponding initial pepper variety Pazardzhishka kapia 794. Both varieties are intended for fresh consumption. Pepper is a major natural source of β-carotene. It was observed that fruit at both commercial and botanical maturity from mutant variety had greater α-carotene and β-carotene concentrations to the initial variety (7.49 and 1.94 times higher, respectively) meaning that the mutant was superior in fruit quality to the initial genotype. Two hydroxylase enzymes, converting α- and β-carotene to lutein and zeaxanthin, respectively, are known to exist in pepper and are encoded by two genes on chromosomes 3 and 6-CrtZ_chr03 and CrtZ_chr06. The molecular characterization of the mutant variety through locus-specific Polymerase chain reaction amplification, gene cloning and sequencing as well as expression was performed. Our results suggest that the increased ß-carotene accumulation in the mutant variety Oranzheva kapia results from a biosynthetic pathway breakdown due to deletion of CrtZ_chr03 gene.

Evolutionary trade-offs between male secondary sexual traits revealed by a phylogeny of the hyperdiverse tribe Eumaeini (Lepidoptera: Lycaenidae)

Male butterflies in the hyperdiverse tribe Eumaeini possess an unusually complex and diverse repertoire of secondary sexual characteristics involved in pheromone production and dissemination. Maintaining multiple sexually selected traits is likely to be metabolically costly, potentially resulting in trade-offs in the evolution of male signals. However, a phylogenetic framework to test hypotheses regarding the evolution and maintenance of male sexual traits in Eumaeini has been lacking. Here, we infer a comprehensive, time-calibrated phylogeny from 379 loci for 187 species representing 91% of the 87 described genera. Eumaeini is a monophyletic group that originated in the late Oligocene and underwent rapid radiation in the Neotropics. We examined specimens of 818 of the 1096 described species (75%) and found that secondary sexual traits are present in males of 91% of the surveyed species. Scent pads and scent patches on the wings and brush organs associated with the genitalia were probably present in the common ancestor of Eumaeini and are widespread throughout the tribe. Brush organs and scent pads are negatively correlated across the phylogeny, exhibiting a trade-off in which lineages with brush organs are unlikely to regain scent pads and vice versa. In contrast, scent patches seem to facilitate the evolution of scent pads, although they are readily lost once scent pads have evolved. Our results illustrate the complex interplay between natural and sexual selection in the origin and maintenance of multiple male secondary sexual characteristics and highlight the potential role of sexual selection spurring diversification in this lineage.

Using the combined gene approach and multiple analytical methods to improve the phylogeny and classification of Bombus (Hymenoptera, Apidae) in China

Bumble bees are vital to our agro-ecological system, with approximately 250 species reported around the world in the single genus Bombus. However, the health of bumble bees is threatened by multiple factors: habitat loss, climate change, pesticide use, and disease caused by pathogens and parasites. It is therefore vitally important to have a fully developed phylogeny for bumble bee species as part of our conservation efforts. The purpose of this study was to explore the phylogenetic relationships of the dominant bumble bees on the Tibetan plateau and in northern China as well as their placement and classification within the genus Bombus. The study used combined gene analysis consisting of sequence fragments from six genes, 16S rRNA, COI, EF-1α, Argk, Opsin and PEPCK, and the phylogenetic relationships of 209 Bombus species were explored. Twenty-six species, including 152 gene sequences, were collected from different regions throughout China, and 1037 gene sequences representing 183 species were obtained from GenBank or BOLD. The results suggest that the 209 analyzed species belong to fifteen subgenera and that most of the subgenera in Bombus are monophyletic, which is in accordance with conventional morphology-based classifications. The phylogenetic trees also show that nearly all subgenera easily fall into two distinct clades: short-faced and long-faced. The study is the first to investigate the phylogenetic placement of Bombus turneri (Richards), Bombus opulentus Smith, Bombus pyrosoma Morawitz, Bombus longipennis Friese, Bombus minshanensis Bischoff, and Bombus lantschouensis Vogt, all of which are widely distributed throughout different regions of China. The knowledge and understanding gained from the findings can provide a molecular basis to accurately classify Bombus in China and to define strategies to conserve biodiversity and promote pollinator populations.

Multiple invasions of a generalist herbivore-Secondary contact between two divergent lineages of Nezara viridula Linnaeus in Australia

The presence of distinct evolutionary lineages within herbivorous pest insect taxa requires close attention. Scientific understanding, biosecurity planning and practice, and pest management decision-making each suffer when such situations remain poorly understood. The pest bug Nezara viridula Linnaeus has been recorded from numerous host plants and has two globally distributed mitochondrial (mtDNA) lineages. These mtDNA lineages co-occur in few locations globally, and the consequences of their divergence and recent secondary contact have not been assessed. We present evidence that both mtDNA lineages of N. viridula are present in Australia and their haplotype groups have a mostly separate distribution from one another. The north-western population has only Asian mtDNA haplotypes, and the population with an eastern distribution is characterized mostly by European mtDNA haplotypes. Haplotypes of both lineages were detected together at only one site in the north of eastern Australia, and microsatellite data indicate that this secondary contact has resulted in mating across the lineages. Admixture and the movement of mtDNA haplotypes outside of this limited area of overlap has not, however, been extensive. Some degree of mating incompatibility or differences in the climatic requirements and tolerances of the two lineages, and perhaps a combination of these influences, might limit introgression and the movement of individuals, but this needs to be tested. This work provides the foundation for further ecological investigation of the lineages of N. viridula, particularly the consequences of admixture on the ecology of this widespread pest. We propose that for now, the Asian and European lineages of N. viridula would best be investigated as subspecies, so that "pure" and admixed populations of this bug can each be considered directly with respect to management and research priorities.

Phylogeny and Biogeography of Spinicaudata (Crustacea: Branchiopoda)

Spinicaudata (spiny clam shrimp) is a taxon of Branchiopoda occurring since the Devonian and today it occurs nearly globally in temporary water bodies. We present the most species-rich phylogenetic analyses of this taxon based on four molecular loci: COI, 16S rRNA, EF1α and 28S rRNA. Our results support previous findings that Cyzicidae sensu lato is paraphyletic. To render Cyzicidae monophyletic we establish a fourth extant spinicaudatan family to accommodate Eocyzicus. Within Cyzicidae, none of the genera Cyzicus, Caenestheria or Caenestheriella are monophyletic, and the morphological characters used to define these genera (condyle length and rostrum shape) are not associated with well-delimited clades within Cyzicidae. There is insufficient resolution to elucidate the relationships within Leptestheriidae. However, there is sufficient evidence to show that the leptestheriid genera Eoleptestheria and Leptestheria are non-monophyletic, and there is no support for the genus Leptestheriella. Molecular clock analyses suggest that the wide geographic distribution of many spinicaudatan taxa across multiple continents is largely based on vicariance associated with the break-up of Pangea and Gondwana. Trans-oceanic dispersal has occurred in some taxa (e.g., Eulimnadia and within Leptestheriidae) but has been relatively rare. Our results highlight the need to revise the taxonomy of Cyzicidae and Leptestheriidae and provide evidence that the global spinicaudatan diversity may be underestimated due to the presence of numerous cryptic species. We establish Eocyzicidae fam. nov. to accommodate the genus Eocyzicus. Consequently, Cyzicidae comprises only two genera -Cyzicus and Ozestheria. Ozestheria occurs also in Africa and Asia and Ozestheria pilosa new comb. is assigned to this genus.

Evolution and losses of spines in slug caterpillars (Lepidoptera: Limacodidae)

Larvae of the cosmopolitan family Limacodidae, commonly known as "slug" caterpillars, are well known because of the widespread occurrence of spines with urticating properties, a morpho-chemical adaptive trait that has been demonstrated to protect the larvae from natural enemies. However, while most species are armed with rows of spines ("nettle" caterpillars), slug caterpillars are morphologically diverse with some species lacking spines and thus are nonstinging. It has been demonstrated that the evolution of spines in slug caterpillars may have a single origin and that this trait is possibly derived from nonstinging slug caterpillars, but these conclusions were based on limited sampling of mainly New World taxa; thus, the evolution of spines and other traits within the family remains unresolved. Here, we analyze morphological variation in slug caterpillars within an evolutionary framework to determine character evolution of spines with samples from Asia, Australia, North America, and South America. The phylogeny of the Limacodidae was reconstructed based on a multigene dataset comprising five molecular markers (5.6 Kbp: COI, 28S, 18S, EF-1α, and wingless) representing 45 species from 40 genera and eight outgroups. Based on this phylogeny, we infer that limacodids evolved from a common ancestor in which the larval type possessed spines, and then slug caterpillars without spines evolved independently multiple times in different continents. While larvae with spines are well adapted to avoiding generalist predators, our results imply that larvae without spines may be suited to different ecological niches. Systematic relationships of our dataset indicate six major lineages, several of which have not previously been identified.

Untangling the Hypogeococcus pungens species complex (Hemiptera: Pseudococcidae) for Argentina, Australia, and Puerto Rico based on host plant associations and genetic evidence

Hypogeococcus pungens, a mealybug native of southern South America, is devastating native cacti in Puerto Rico and threatening cactus diversity in the Caribbean, and potentially in Central and North America. The taxonomic status of H. pungens is controversial since it has been reported feeding not only on Cactaceae but also on other plant families throughout its distribution range. However, in Australia, where the species had been exported from Argentina to control weedy American cacti, it was never found on host plants other than Cactaceae. These conflicting pieces of evidence not only cast doubt on the species identity that invaded Puerto Rico, but also have a negative impact on the search for natural enemies to be used in biological control programs against this pest. Here we present reproductive incompatibility and phylogenetic evidences that give support to the hypothesis that H. pungens is a species complex in which divergence appears to be driven by the host plants. The nuclear EF1α and 18S and the mitochondrial COI genes were used as markers to evaluate the phylogenetic relationships among H. pungens populations collected in Argentina, Australia and Puerto Rico feeding on Cactaceae and/or Amaranthaceae. Additionally, we conducted reciprocal crosses between mealybugs from both hosts. Species delimitation analysis revealed two well-supported putative species within H. pungens, one including mealybugs feeding on Amaranthaceae (H. pungens sensu stricto), and a new undescribed species using Cactaceae as hosts. Additionally, we found asymmetric reproductive incompatibility between these putative species suggesting recent reproductive isolation. The Bayesian species delimitation also suggested that the Australian mealybug population may derive from another undescribed species. Overall, the patterns of genetic differentiation may be interpreted as the result of recent speciation events prompted by host plant shifts. Finally, the finding of a single haplotype in the Puerto Rico population suggests only one invasive event. We still need to identify the geographical origin of the pest in order to enable the use of biological control to reduce the threat to cacti diversity in the Caribbean.

Phylogeny, host use, and diversification in the moth family Momphidae (Lepidoptera: Gelechioidea)

Insect herbivores and their hostplants constitute much of Earth's described biological diversity, but how these often-specialized associations diversify is not fully understood. We combined detailed hostplant data and comparative phylogenetic analyses of the lepidopteran family Momphidae to explore how shifts in the use of hostplant resources, not just hostplant taxon, contribute to the diversification of a phytophagous insect lineage. We inferred two phylogenetic hypotheses emphasizing relationships among species in the nominate genus, Mompha Hübner. A six-gene phylogeny was constructed with reared exemplars and collections from hostplants in the family Onagraceae from western and southwestern USA, and a cytochrome c oxidase subunit 1 (COI) phylogeny was inferred from collections and publicly available accessions in the Barcode of Life Data System. Species delimitation analyses combined with morphological data revealed ca. 56 undescribed species-level taxa, many of which are hostplant specialists on Onagraceae in the southwestern USA. Our phylogenetic reconstructions divided Momphidae into six major clades: 1) an Onagraceae flower- and fruit-boring clade, 2) a Melastomataceae-galling clade, 3) a leafmining clade A, 4) a leafmining clade B, 5) a Zapyrastra Meyrick clade, and 6) a monobasic lineage represented by Mompha eloisella (Clemens). Ancestral trait reconstructions using the COI phylogeny identified leafmining on Onagraceae as the ancestral state for Momphidae. Our study finds that shifts along three hostplant resource axes (plant taxon, plant tissue type, and larval feeding mode) have contributed to the evolutionary success and diversification of momphids.

A Revision of North American Lactura (Lepidoptera, Zygaenoidea, Lacturidae)

The Lactura Walker, 1854 fauna north of Mexico is revised. Six species are documented, one new species Lacturanalli Matson & Wagner, sp. n. is described, and two new synonymies are proposed: Lacturapsammitis (Zeller, 1872), syn. n. and L.rhodocentra (Meyrick, 1913), syn. n. One new subspecies Lacturasubfervenssapeloensis Matson & Wagner, ssp. n. is also described. Adult and larval stages, male and female genitalia, are illustrated, a preliminary phylogeny is presented based on nuclear and mitochondrial data, distribution records provided for verified specimens, and the biology and life history for each species is briefly characterized. Phylogenetic analyses, larval phenotypes, and life history information reveal that much of the historic taxonomic confusion rampant across this group in North America traces to the phenotypic variation in just one species, L.subfervens (Walker, 1854).

The Klingon batbugs: Morphological adaptations in the primitive bat bugs, Bucimex chilensis and Primicimex cavernis, including updated phylogeny of Cimicidae

The Cimicidae is a family of blood-dependent ectoparasites in which dispersion capacity is greatly associated with host movements. Bats are the ancestral and most prevalent hosts for cimicids. Cimicids have a worldwide distribution matching that of their hosts, but the global classification is incomplete, especially for species outside the most common Cimicidae taxa. In this study, we place a little-studied cimicid species, Bucimex chilensis, within a comprehensive molecular phylogeny of Cimicidae by sequencing the genomic regions of this and other closely related species. For this study, we collected B. chilensis females from Myotis chiloensis in Tierra del Fuego, 1,300 km further south than previously known southernmost distribution boundary. We also sequenced COI regions from Primicimex cavernis, a species which together with B. chilensis comprise the entire subfamily Primiciminae. Using Bayesian posterior probability and maximum-likelihood approaches, we found that B. chilensis and P. cavernis clustered close to each other in the molecular analyses, receiving support from similar morphological features, agreeing with the morphology-based taxonomic placement of the two species within the subfamily Primiciminae. We also describe a previously unrecognized morphological adaptation of the tarsal structure, which allows the austral bat ectoparasite, B. chilensis, to cling on to the pelage of its known host, the Chilean myotis (Myotis chiloensis). Through a morphological study and behavioral observation, we elucidate how this tarsal structure operates, and we hypothesize that by clinging in the host pelage, B. chilensis is able to disperse effectively to new areas despite low host density. This is a unique feature shared by P. cavernis, the only other species in Primiciminae.

The genetic architecture of ecological adaptation: intraspecific variation in host plant use by the lepidopteran crop pest Chloridea virescens

Intraspecific variation in ecologically important traits is a cornerstone of Darwin's theory of evolution by natural selection. The evolution and maintenance of this variation depends on genetic architecture, which in turn determines responses to natural selection. Some models suggest that traits with complex architectures are less likely to respond to selection than those with simple architectures, yet rapid divergence has been observed in such traits. The simultaneous evolutionary lability and genetic complexity of host plant use in the Lepidopteran subfamily Heliothinae suggest that architecture may not constrain ecological adaptation in this group. Here we investigate the response of Chloridea virescens, a generalist that feeds on diverse plant species, to selection for performance on a novel host, Physalis angulata (Solanaceae). P. angulata is the preferred host of Chloridea subflexa, a narrow specialist on the genus Physalis. In previous experiments, we found that the performance of C. subflexa on P. angulata depends on many loci of small effect distributed throughout the genome, but whether the same architecture would be involved in the generalist's adoption of P. angulata was unknown. Here we report a rapid response to selection in C. virescens for performance on P. angulata, and establish that the genetic architecture of intraspecific variation is quite similar to that of the interspecific differences in terms of the number, distribution, and effect sizes of the QTL involved. We discuss the impact of genetic architecture on the ability of Heliothine moths to respond to varying ecological selection pressures.

Phlebotomus (Paraphlebotomus) chabaudi and Phlebotomus riouxi: closely related species or synonyms?

Phlebotomus riouxi Depaquit, Killick-Kendrick & Léger 1998 was described as a species closely related to Phlebotomus chabaudi Croset, Abonnenc & Rioux 1970, differing mainly by the size and number of setae of the coxite basal lobe. Molecular studies carried out on several populations from Algeria and Tunisia and based on mitochondrial genes cytochrome b (Cytb) and cytochrome oxidase I (COI) supported the typological validity of these two species. Recently, specimens from a single population in southern Tunisia were morphologically identified as Ph. riouxi, Ph. chabaudi and intermediates, but were clustered in the same clade according to their Cytb and nuclear gene elongation factor-1 α (EF-1α) sequences. These species were thus synonymized. To further explore this synonymy, we carried out a molecular study on specimens from Algeria and Tunisia using the same molecular markers and a part of 28S rDNA. We did not find any morphologically intermediate specimens in our sampling. We highlighted differences between the genetic divergence rates within and between the two species for the three markers and we identified new haplotypes. The sequence analysis did not reveal any signature of introgression in allopatric nor in sympatric populations such as in the Ghomrassen population. Phylogenetic analyses based on our specimens revealed that the two main clades are Ph. chabaudi and Ph. riouxi, in agreement with the morphological identification. These results support the validity of Ph. riouxi and Ph. chabaudi as typological species.

Coevolutionary arms race versus host defense chase in a tropical herbivore-plant system

Coevolutionary models suggest that herbivores drive diversification and community composition in plants. For herbivores, many questions remain regarding how plant defenses shape host choice and community structure. We addressed these questions using the tree genus Inga and its lepidopteran herbivores in the Amazon. We constructed phylogenies for both plants and insects and quantified host associations and plant defenses. We found that similarity in herbivore assemblages between Inga species was correlated with similarity in defenses. There was no correlation with phylogeny, a result consistent with our observations that the expression of defenses in Inga is independent of phylogeny. Furthermore, host defensive traits explained 40% of herbivore community similarity. Analyses at finer taxonomic scales showed that different lepidopteran clades select hosts based on different defenses, suggesting taxon-specific histories of herbivore-host plant interactions. Finally, we compared the phylogeny and defenses of Inga to phylogenies for the major lepidopteran clades. We found that closely related herbivores fed on Inga with similar defenses rather than on closely related plants. Together, these results suggest that plant defenses might be more evolutionarily labile than the herbivore traits related to host association. Hence, there is an apparent asymmetry in the evolutionary interactions between Inga and its herbivores. Although plants may evolve under selection by herbivores, we hypothesize that herbivores may not show coevolutionary adaptations, but instead "chase" hosts based on the herbivore's own traits at the time that they encounter a new host, a pattern more consistent with resource tracking than with the arms race model of coevolution.

Dated phylogeny and dispersal history of the butterfly subfamily Nymphalinae (Lepidoptera: Nymphalidae)

The origin and dispersal history of the large butterfly subfamily Nymphalinae are not fully understood, due to internal phylogenetic and time calibration issues. We conducted phylogenetic and dating analyses using mitochondrial and nuclear genes of biogeographically diverse groups of the Nymphalinae in order to resolve some controversial relationships and the paleobiogeographic pattern of the subfamily. Our results support the sister relationship of Vanessa (Tribe Nymphalini) and the Nymphalis-group, and the grouping of the three old-world genera (Rhinopalpa, Kallimoides and Vanessula) within Tribe Victorinini. Molecular dating analyses invoking two additional calibrations under the butterfly-host plant coevolutionary scenarios result in a relatively deeper divergence of the subfamily's two major clades (Nymphalini and the Kallimoids), compatible with the Cretaceous floral turnover scenario during the so-called Cretaceous Terrestrial Revolution. Phylobiogeographic analyses reveal that the Oriental region is probably the center of early divergences for Nymphalinae after the Cretaceous-Paleogene (K-Pg) mass extinction, followed by repeated dispersals into the rest of the Old World and the New World during various periods beginning in Eocene. The biogeographic history indicates that temperature changes and host-plant diversification may have facilitated the dispersals of this butterfly subfamily, with accelerated global colonization during the middle to late Miocene.

Species delimitation in asexual insects of economic importance: The case of black scale (Parasaissetia nigra), a cosmopolitan parthenogenetic pest scale insect

Asexual lineages provide a challenge to species delimitation because species concepts either have little biological meaning for them or are arbitrary, since every individual is monophyletic and reproductively isolated from all other individuals. However, recognition and naming of asexual species is important to conservation and economic applications. Some scale insects are widespread and polyphagous pests of plants, and several species have been found to comprise cryptic species complexes. Parasaissetia nigra (Nietner, 1861) (Hemiptera: Coccidae) is a parthenogenetic, cosmopolitan and polyphagous pest that feeds on plant species from more than 80 families. Here, we implement multiple approaches to assess the species status of P. nigra, including coalescence-based analyses of mitochondrial and nuclear genes, and ecological niche modelling. Our results indicate that the sampled specimens of P. nigra should be considered to comprise at least two ecotypes (or "species") that are ecologically differentiated, particularly in relation to temperature and moisture. The presence of more than one ecotype under the current concept of P. nigra has implications for biosecurity because the geographic extent of each type is not fully known: some countries may currently have only one of the biotypes. Introduction of additional lineages could expand the geographic extent of damage by the pest in some countries.

Genomic Mining of Phylogenetically Informative Nuclear Markers in Bark and Ambrosia Beetles

Deep level insect relationships are generally difficult to resolve, especially within taxa of the most diverse and species rich holometabolous orders. In beetles, the major diversity occurs in the Phytophaga, including charismatic groups such as leaf beetles, longhorn beetles and weevils. Bark and ambrosia beetles are wood boring weevils that contribute 12 percent of the diversity encountered in Curculionidae, one of the largest families of beetles with more than 50000 described species. Phylogenetic resolution in groups of Cretaceous age has proven particularly difficult and requires large quantity of data. In this study, we investigated 100 nuclear genes in order to select a number of markers with low evolutionary rates and high phylogenetic signal. A PCR screening using degenerate primers was applied to 26 different weevil species. We obtained sequences from 57 of the 100 targeted genes. Sequences from each nuclear marker were aligned and examined for detecting multiple copies, pseudogenes and introns. Phylogenetic informativeness (PI) and the capacity for reconstruction of previously established phylogenetic relationships were used as proxies for selecting a subset of the 57 amplified genes. Finally, we selected 16 markers suitable for large-scale phylogenetics of Scolytinae and related weevil taxa.

DNA Barcoding the Heliothinae (Lepidoptera: Noctuidae) of Australia and Utility of DNA Barcodes for Pest Identification in Helicoverpa and Relatives

Helicoverpa and Heliothis species include some of the world’s most significant crop pests, causing billions of dollars of losses globally. As such, a number are regulated quarantine species. For quarantine agencies, the most crucial issue is distinguishing native species from exotics, yet even this task is often not feasible because of poorly known local faunas and the difficulties of identifying closely related species, especially the immature stages. DNA barcoding is a scalable molecular diagnostic method that could provide the solution to this problem, however there has been no large-scale test of the efficacy of DNA barcodes for identifying the Heliothinae of any region of the world to date. This study fills that gap by DNA barcoding the entire heliothine moth fauna of Australia, bar one rare species, and comparing results with existing public domain resources. We find that DNA barcodes provide robust discrimination of all of the major pest species sampled, but poor discrimination of Australian Heliocheilus species, and we discuss ways to improve the use of DNA barcodes for identification of pests.

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

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

Multilocus phylogeny and ecological differentiation of the "Eupelmus urozonus species group" (Hymenoptera, Eupelmidae) in the West-Palaearctic

BACKGROUND

The ecological differentiation of insects with parasitic life-style is a complex process that may involve phylogenetic constraints as well as morphological and/or behavioural adaptations. In most cases, the relative importance of these driving forces remains unexplored. We investigate here this question for the "Eupelmus urozonus species group" which encompasses parasitoid wasps of potential interest in biological control. This was achieved using seven molecular markers, reliable records on 91 host species and a proxy of the ovipositor length.

RESULTS

After using an adequate partitioning scheme, Maximum likelihood and Bayesian approaches provide a well-resolved phylogeny supporting the monophyly of this species group and highlighting its subdivision into three sub-groups. Great variations of both the ovipositor length and the host range (specialist versus generalist) were observed at this scale, with these two features being not significantly constrained by the phylogeny. Ovipositor length was not shown as a significant predictor of the parasitoid host range.

CONCLUSIONS

This study provides firstly the first evidence for the strong lability of both the ovipositor's length and the realised host range in a set of phylogenetically related and sympatric species. In both cases, strong contrasts were observed between sister species. Moreover, no significant correlation was found between these two features. Alternative drivers of the ecological differentiation such as interspecific interactions are proposed and the consequences on the recruitment of these parasitoids on native and exotic pests are discussed.

A phylogenetic assessment of the polyphyletic nature and intraspecific color polymorphism in the Bactrocera dorsalis complex (Diptera, Tephritidae)

The Bactroceradorsalis complex (Tephritidae) comprises 85 species of fruit flies, including five highly destructive polyphagous fruit pests. Despite significant work on a few key pest species within the complex, little has been published on the majority of non-economic species in the complex, other than basic descriptions and illustrations of single specimens regarded as typical representatives. To elucidate the species relationships within the Bactroceradorsalis complex, we used 159 sequences from one mitochondrial (COI) and two nuclear (elongation factor-1α and period) genes to construct a phylogeny containing 20 described species from within the complex, four additional species that may be new to science, and 26 other species from Bactrocera and its sister genus Dacus. The resulting concatenated phylogeny revealed that most of the species placed in the complex appear to be unrelated, emerging across numerous clades. This suggests that they were placed in the Bactroceradorsalis complex based on the similarity of convergent characters, which does not appear to be diagnostic. Variations in scutum and abdomen color patterns within each of the non-economic species are presented and demonstrate that distantly-related, cryptic species overlap greatly in traditional morphological color patterns used to separate them in keys. Some of these species may not be distinguishable with confidence by means other than DNA data.

Specific olfactory neurons and glomeruli are associated to differences in behavioral responses to pheromone components between two Helicoverpa species

Sex pheromone communication of moths helps to understand the mechanisms underlying reproductive isolation and speciation. Helicoverpa armigera and Helicoverpa assulta use (Z)-11-hexadecenal (Z11-16:Ald) and (Z)-9-hexadecenal (Z9-16:Ald) as pheromone components in reversed ratios, 97:3 and 5:95, respectively. H. armigera also produces trace amount of (Z)-9-tetradecenal (Z9-14:Ald) in the sex pheromone gland, but H. assulta does not. Wind tunnel studies revealed that the addition of small amounts (0.3%) of Z9-14:Ald to the main pheromone blend of H. armigera increased the males' attraction, but at higher doses (1%, 10%) the same compound acted as an inhibitor. In H. assulta, Z9-14:Ald reduced male attraction when presented as 1% to the pheromone blend, but was ineffective at lower concentrations (0.3%). Three types (A–C) of sensilla trichodea in antennae were identified by single sensillum recording, responding to Z11-16:Ald, Z9-14:Ald, and both Z9-16:Ald and Z9-14:Ald, respectively. Calcium imaging in the antennal lobes (ALs) revealed that the input information of the three chemicals was transmitted to three units of the macroglomerular complex (MGC) in ALs in both species: a large glomerulus for the major pheromone components, a small one for the minor pheromone components, and a third one for the behavioral antagonists. The type A and C neurons tuned to Z11-16:Ald and Z9-16:Ald had a reversed target in the MGC between the two species. In H. armigera, low doses (1, 10 μg) of Z9-14:Ald dominantly activated the glomerulus which processes the minor pheromone component, while a higher dose (100 μg) also evoked an equal activity in the antagonistic glomerulus. In H. assulta, instead, Z9-14:Ald always strongly activated the antagonistic glomerulus. These results suggest that Z9-14:Ald plays different roles in the sexual communication of two Helicoverpa species through activation of functionally different olfactory pathways.

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.

Phylogeny and evolution of pharmacophagy in tiger moths (Lepidoptera: Erebidae: Arctiinae)

The focus of this study was to reconstruct a phylogenetic hypothesis for the moth subfamily Arctiinae (tiger moths, woolly bears) to investigate the evolution of larval and adult pharmacophagy of pyrrolizidine alkaloids (PAs) and the pathway to PA chemical specialization in Arctiinae. Pharmacophagy, collection of chemicals for non-nutritive purposes, is well documented in many species, including the model species Utetheisa ornatrix L. A total of 86 exemplar ingroup species representing tiger moth tribes and subtribes (68 genera) and nine outgroup species were selected. Ingroup species included the most species-rich generic groups to represent the diversity of host-plant associations and pharmacophagous behaviors found throughout Arctiinae. Up to nine genetic markers were sequenced: one mitochondrial (COI barcode region), one nuclear rRNA (D2 region, 28S rRNA), and seven nuclear protein-coding gene fragments: elongation factor 1-α protein, wingless, ribosomal protein subunit S5, carbamoylphosphate synthase domain regions, glyceraldehyde-3-phosphate dehydrogenase, isocitrate dehydrogenase and cytosolic malate dehydrogenase. A total of 6984 bp was obtained for most species. These data were analyzed using model-based phylogenetic methods: maximum likelihood (ML) and Bayesian inference (BI). Ancestral pharmacophagous behaviors and obligate PA associations were reconstructed using the resulting Bayes topology and Reconstructing Ancestral States in Phylogenies (RASP) software. Our results corroborate earlier studies on the evolution of adult pharmacophagous behaviors, suggesting that this behavior arose multiple times and is concentrated in the phaegopterine-euchromiine-ctenuchine clade (PEC). Our results suggest that PA specialization may have arisen early in the phylogeny of the subfamily and that facultative larval pharmacophagous behaviors are the derived condition.

Mitochondrial and nuclear phylogenetic analysis with Sanger and next-generation sequencing shows that, in Área de Conservación Guanacaste, northwestern Costa Rica, the skipper butterfly named Urbanus belli (family Hesperiidae) comprises three morphologically cryptic species

Background

Skipper butterflies (Hesperiidae) are a relatively well-studied family of Lepidoptera. However, a combination of DNA barcodes, morphology, and natural history data has revealed several cryptic species complexes within them. Here, we investigate three DNA barcode lineages of what has been identified as Urbanus belli (Hesperiidae, Eudaminae) in Área de Conservación Guanacaste (ACG), northwestern Costa Rica.

Results

Although no morphological traits appear to distinguish among the three, congruent nuclear and mitochondrial lineage patterns show that “Urbanus belli” in ACG is a complex of three sympatric species. A single strain of Wolbachia present in two of the three cryptic species indicates that Urbanus segnestami Burns (formerly Urbanus belliDHJ01), Urbanus bernikerni Burns (formerly Urbanus belliDHJ02), and Urbanus ehakernae Burns (formerly Urbanus belliDHJ03) may be biologically separated by Wolbachia, as well as by their genetics. Use of parallel sequencing through 454-pyrosequencing improved the utility of ITS2 as a phylogenetic marker and permitted examination of the intra- and interlineage relationships of ITS2 variants within the species complex. Interlineage, intralineage and intragenomic compensatory base pair changes were discovered in the secondary structure of ITS2.

Conclusion

These findings corroborate the existence of three cryptic species. Our confirmation of a novel cryptic species complex, initially suggested by DNA barcode lineages, argues for using a multi-marker approach coupled with next-generation sequencing for exploration of other suspected species complexes.

Phylogeny, systematics and biogeography of the genus panolis (lepidoptera: noctuidae) based on morphological and molecular evidence

The genus Panolis is a small group of noctuid moths with six recognized species distributed from Europe to East Asia, and best known for containing the widespread Palearctic pest species P. flammea, the pine beauty moth. However, a reliable classification and robust phylogenetic framework for this group of potentially economic importance are currently lacking. Here, we use morphological and molecular data (mitochondrial genes cytochrome c oxidase subunit I and 16S ribosomal RNA, nuclear gene elongation factor-1 alpha) to reconstruct the phylogeny of this genus, with a comprehensive systematic revision of all recognized species and a new one, P. ningshan sp. nov. The analysis results of maximum parsimony, maximum likelihood and Bayesian inferring methods for the combined morphological and molecular data sets are highly congruent, resulting in a robust phylogeny and identification of two clear species groups, i.e., the P. flammea species group and the P. exquisita species group. We also estimate the divergence times of Panolis moths using two conventional mutation rates for the arthropod mitochondrial COI gene with a comparison of two molecular clock models, as well as reconstruct their ancestral areas. Our results suggest that 1) Panolis is a young clade, originating from the Oriental region in China in the Late Miocene (6-10Mya), with an ancestral species in the P. flammea group extending northward to the Palearctic region some 3-6 Mya; 2) there is a clear possibility for a representative of the Palearctic clade to become established as an invasive species in the Nearctic taiga.

A multilocus phylogeny of the world Sycoecinae fig wasps (Chalcidoidea: Pteromalidae)

The Sycoecinae is one of five chalcid subfamilies of fig wasps that are mostly dependent on Ficus inflorescences for reproduction. Here, we analysed two mitochondrial (COI, Cytb) and four nuclear genes (ITS2, EF-1α, RpL27a, mago nashi) from a worldwide sample of 56 sycoecine species. Various alignment and partitioning strategies were used to test the stability of major clades. All topologies estimated using maximum likelihood and Bayesian methods were similar and well resolved but did not support the existing classification. A high degree of morphological convergence was highlighted and several species appeared best described as species complexes. We therefore proposed a new classification for the subfamily. Our analyses revealed several cases of probable speciation on the same host trees (up to 8 closely related species on one single tree of F. sumatrana), which raises the question of how resource partitioning occurs to avoid competitive exclusion. Comparisons of our results with fig phylogenies showed that, despite sycoecines being internally ovipositing wasps host-switches are common incidents in their evolutionary history. Finally, by studying the evolutionary properties of the markers we used and profiling their phylogenetic informativeness, we predicted their utility for resolving phylogenetic relationships of Chalcidoidea at various taxonomic levels.

Peripheral coding of sex pheromone blends with reverse ratios in two helicoverpa species

The relative proportions of components in a pheromone blend play a major role in sexual recognition in moths. Two sympatric species, Helicoverpa armigera and Helicoverpa assulta, use (Z)-11-hexadecenal (Z11-16: Ald) and (Z)-9-hexadecenal (Z9-16: Ald) as essential sex pheromone components but in very different ratios, 97∶3 and 7∶93 respectively. Using wind tunnel tests, single sensillum recording and in vivo calcium imaging, we comparatively studied behavioral responses and physiological activities at the level of antennal sensilla and antennal lobe (AL) in males of the two species to blends of the two pheromone components in different ratios (100∶0, 97∶3, 50∶50, 7∶93, 0∶100). Z11-16: Ald and Z9-16: Ald were recognized by two populations of olfactory sensory neurons (OSNs) in different trichoid sensilla on antennae of both species. The ratios of OSNs responding to Z11-16:Ald and Z9-16:Ald OSNs were 100∶28.9 and 21.9∶100 in H. armigera and H. assulta, respectively. The Z11-16:Ald OSNs in H. armigera exhibited higher sensitivity and efficacy than those in H. assulta, while the Z9-16:Ald OSNs in H. armigera had the same sensitivity but lower efficacy than those in H. assulta. At the dosage of 10 µg, Z11-16: Ald and Z9-16: Ald evoked calcium activity in 8.5% and 3.0% of the AL surface in H. armigera, while 5.4% and 8.6% of AL in H. assulta, respectively. The calcium activities in the AL reflected the peripheral input signals of the binary pheromone mixtures and correlated with the behavioral output. These results demonstrate that the binary pheromone blends were precisely coded by the firing frequency of individual OSNs tuned to Z11-16: Ald or Z9-16: Ald, as well as their population sizes. Such information was then accurately reported to ALs of H. armigera and H. assulta, eventually producing different behaviors.

Identification of Petriella setifera LH and characterization of its crude carboxymethyl cellulase for application in denim biostoning

The phylogenetic tree of the partial elongation factor-1 alpha gene fits better than the partial 18S rDNA for generic classification. From the results of the molecular tree and analysis of morphological characters, Petriella setifera LH was identified. It can be induced to produce carboxymethyl cellulase (CMCase). The crude CMCase only shows a 44.1-kDa band by activity staining after SDS-PAGE. It is optimally active at 55°C and pH 6.0, and is stable from pH 5.0-8.0 and at 45°C or below. The crude CMCase, which is not affected by Co(2+), is strongly activated in the presence of 10 mM Na(+), K(+), Ca(2+), Mg(2+), EDTA, and Mn(2+). It is strongly inhibited by 10 mM Fe(2+), Pb(2+), Al(3+), Zn(2+), Ag(+), Fe(3+), and Cu(2+). When compared with denim treatment by Novoprime A800 (a commercial neutral cellulase), crude CMCase exhibits a similar fabric weight loss and indigo dye removal. These results indicate that crude CMCase has potential application in denim biostoning.

The genetic architecture of a complex ecological trait: host plant use in the specialist moth, Heliothis subflexa

We used genetic mapping to examine the genetic architecture of differences in host plant use between two species of noctuid moths, Heliothis subflexa, a specialist on Physalis spp., and its close relative, the broad generalist H. virescens. We introgressed H. subflexa chromosomes into the H. virescens background and analyzed 1462 backcross insects. The effects of H. subflexa-origin chromosomes were small when measured as the percent variation explained in backcross populations (0.2-5%), but were larger when considered in relation to the interspecific difference explained (1.5-165%). Most significant chromosomes had effects on more than one trait, and their effects varied between years, sexes, and genetic backgrounds. Different chromosomes could produce similar phenotypes, suggesting that the same trait might be controlled by different chromosomes in different backcross populations. It appears that many loci of small effect contribute to the use of Physalis by H. subflexa. We hypothesize that behavioral changes may have paved the way for physiological adaptation to Physalis by the generalist ancestor of H. subflexa and H. virescens.

Increased gene sampling strengthens support for higher-level groups within leaf-mining moths and relatives (Lepidoptera: Gracillariidae)

BACKGROUND

Researchers conducting molecular phylogenetic studies are frequently faced with the decision of what to do when weak branch support is obtained for key nodes of importance. As one solution, the researcher may choose to sequence additional orthologous genes of appropriate evolutionary rate for the taxa in the study. However, generating large, complete data matrices can become increasingly difficult as the number of characters increases. A few empirical studies have shown that augmenting genes even for a subset of taxa can improve branch support. However, because each study differs in the number of characters and taxa, there is still a need for additional studies that examine whether incomplete sampling designs are likely to aid at increasing deep node resolution. We target Gracillariidae, a Cretaceous-age (~100 Ma) group of leaf-mining moths to test whether the strategy of adding genes for a subset of taxa can improve branch support for deep nodes. We initially sequenced ten genes (8,418 bp) for 57 taxa that represent the major lineages of Gracillariidae plus outgroups. After finding that many deep divergences remained weakly supported, we sequenced eleven additional genes (6,375 bp) for a 27-taxon subset. We then compared results from different data sets to assess whether one sampling design can be favored over another. The concatenated data set comprising all genes and all taxa and three other data sets of different taxon and gene sub-sampling design were analyzed with maximum likelihood. Each data set was subject to five different models and partitioning schemes of non-synonymous and synonymous changes. Statistical significance of non-monophyly was examined with the Approximately Unbiased (AU) test.

RESULTS

Partial augmentation of genes led to high support for deep divergences, especially when non-synonymous changes were analyzed alone. Increasing the number of taxa without an increase in number of characters led to lower bootstrap support; increasing the number of characters without increasing the number of taxa generally increased bootstrap support. More than three-quarters of nodes were supported with bootstrap values greater than 80% when all taxa and genes were combined. Gracillariidae, Lithocolletinae + Leucanthiza, and Acrocercops and Parectopa groups were strongly supported in nearly every analysis. Gracillaria group was well supported in some analyses, but less so in others. We find strong evidence for the exclusion of Douglasiidae from Gracillarioidea sensu Davis and Robinson (1998). Our results strongly support the monophyly of a G.B.R.Y. clade, a group comprised of Gracillariidae + Bucculatricidae + Roeslerstammiidae + Yponomeutidae, when analyzed with non-synonymous changes only, but this group was frequently split when synonymous and non-synonymous substitutions were analyzed together.

CONCLUSIONS

1) Partially or fully augmenting a data set with more characters increased bootstrap support for particular deep nodes, and this increase was dramatic when non-synonymous changes were analyzed alone. Thus, the addition of sites that have low levels of saturation and compositional heterogeneity can greatly improve results. 2) Gracillarioidea, as defined by Davis and Robinson (1998), clearly do not include Douglasiidae, and changes to current classification will be required. 3) Gracillariidae were monophyletic in all analyses conducted, and nearly all species can be placed into one of six strongly supported clades though relationships among these remain unclear. 4) The difficulty in determining the phylogenetic placement of Bucculatricidae is probably attributable to compositional heterogeneity at the third codon position. From our tests for compositional heterogeneity and strong bootstrap values obtained when synonymous changes are excluded, we tentatively conclude that Bucculatricidae is closely related to Gracillariidae + Roeslerstammiidae + Yponomeutidae.

Molecular approaches to identify cryptic species and polymorphic species within a complex community of fig wasps

Cryptic and polymorphic species can complicate traditional taxonomic research and both of these concerns are common in fig wasp communities. Species identification is very difficult, despite great effort and the ecological importance of fig wasps. Herein, we try to identify all chalcidoid wasp species hosted by one species of fig, using both morphological and molecular methods. We compare the efficiency of four different DNA regions and find that ITS2 is highly effective for species identification, while mitochondrial COI and Cytb regions appear less reliable, possibly due to the interference signals from either nuclear copies of mtDNA, i.e. NUMTs, or the effects of Wolbachia infections. The analyses suggest that combining multiple markers is the best choice for inferring species identifications as any one marker may be unsuitable in a given case.

Assessing the value of DNA barcodes and other priority gene regions for molecular phylogenetics of Lepidoptera

BACKGROUND

Despite apparently abundant amounts of observable variation and species diversity, the order Lepidoptera exhibits a morphological homogeneity that has provided only a limited number of taxonomic characters and led to widespread use of nucleotides for inferring relationships. This study aims to characterize and develop methods to quantify the value of priority gene regions designated for Lepidoptera molecular systematics. In particular, I assess how the DNA barcode segment of the mitochondrial COI gene performs across a broad temporal range given its number one position of priority, most sequenced status, and the conflicting opinions on its phylogenetic performance.

METHODOLOGY/PRINCIPAL FINDINGS

Gene regions commonly sequenced for lepidoptera phylogenetics were scored using multiple measures across three categories: practicality, which includes universality of primers and sequence quality; phylogenetic utility; and phylogenetic signal. I found that alternative measures within a category often appeared correlated, but high scores in one category did not necessarily translate into high scores in another. The DNA barcode was easier to sequence than other genes, and had high scores for utility but low signal above the genus level.

CONCLUSIONS/SIGNIFICANCE

Given limited financial resources and time constraints, careful selection of gene regions for molecular phylogenetics is crucial to avoid wasted effort producing partially informative data. This study introduces an approach to assessing the value of gene regions prior to the initiation of new studies and presents empirical results to help guide future selections.

Phylogeny and biogeography of hawkmoths (Lepidoptera: Sphingidae): evidence from five nuclear genes

BACKGROUND

The 1400 species of hawkmoths (Lepidoptera: Sphingidae) comprise one of most conspicuous and well-studied groups of insects, and provide model systems for diverse biological disciplines. However, a robust phylogenetic framework for the family is currently lacking. Morphology is unable to confidently determine relationships among most groups. As a major step toward understanding relationships of this model group, we have undertaken the first large-scale molecular phylogenetic analysis of hawkmoths representing all subfamilies, tribes and subtribes.

METHODOLOGY/PRINCIPAL FINDINGS

The data set consisted of 131 sphingid species and 6793 bp of sequence from five protein-coding nuclear genes. Maximum likelihood and parsimony analyses provided strong support for more than two-thirds of all nodes, including strong signal for or against nearly all of the fifteen current subfamily, tribal and sub-tribal groupings. Monophyly was strongly supported for some of these, including Macroglossinae, Sphinginae, Acherontiini, Ambulycini, Philampelini, Choerocampina, and Hemarina. Other groupings proved para- or polyphyletic, and will need significant redefinition; these include Smerinthinae, Smerinthini, Sphingini, Sphingulini, Dilophonotini, Dilophonotina, Macroglossini, and Macroglossina. The basal divergence, strongly supported, is between Macroglossinae and Smerinthinae+Sphinginae. All genes contribute significantly to the signal from the combined data set, and there is little conflict between genes. Ancestral state reconstruction reveals multiple separate origins of New World and Old World radiations.

CONCLUSIONS/SIGNIFICANCE

Our study provides the first comprehensive phylogeny of one of the most conspicuous and well-studied insects. The molecular phylogeny challenges current concepts of Sphingidae based on morphology, and provides a foundation for a new classification. While there are multiple independent origins of New World and Old World radiations, we conclude that broad-scale geographic distribution in hawkmoths is more phylogenetically conserved than previously postulated.

New Guinea highland origin of a widespread arthropod supertramp

The biologically and geologically extremely diverse archipelagos of Wallacea, Australasia and Oceania have long stimulated ecologists and evolutionary biologists. Yet, few molecular phylogenetic analyses of the terrestrial fauna have been carried out to understand the evolutionary patterns. We use dense taxon and character sampling of more than 7000 bp DNA sequence data for a group of diving beetles ranging from the Holarctic throughout Asia to as far east as French Polynesia. We here show that an ecologically diverse, common and widespread (Portugal to New Zealand) arthropod supertramp species originated in the highlands of New Guinea, ca 6.0-2.7 Myr ago. The approximately 25 closely related species are narrow endemics in Australasia/Oceania. The ancestor of this clade colonized that region from Eurasia ca 9-7 Myr ago. Our finding contradicts the widely held view of local endemism as an evolutionary dead end, as we find multiple peripatric speciation events within the Pleistocene and complex colonization patterns between the Oriental and Australian zoogeographic regions, including the recolonization of Eurasia, jumping across Wallace's line and colonization of continental Australia out of New Guinea. Our study strongly highlights the importance of dispersal over water gaps in shaping biogeographic patterns.

Oscheius tipulae as an example of eEF1A gene diversity in nematodes

We characterized four eEF1A genes in the alternative rhabditid nematode model organism Oscheius tipulae. This is twice the copy number of eEF1A genes in C. elegans, C. briggsae, and, probably, many other free-living and parasitic nematodes. The introns show features remarkably different from those of other metazoan eEF1A genes. Most of the introns in the eEF1A genes are specific to O. tipulae and are not shared with any of the other genes described in metazoans. Most of the introns are phase 0 (inserted between two codons), and few are inserted in protosplice sites (introns inserted between the nucleotide sequence A/CAG and G/A). Two of these phase 0 introns are conserved in sequence in two or more of the four eEF1A gene copies, and are inserted in the same position in the genes. Neither of these characteristics has been detected in any of the nematode eEF1A genes characterized to date. The coding sequences were also compared with other eEF1A cDNAs from 11 different nematodes to determine the variability of these genes within the phylum Nematoda. Parsimony and distance trees yielded similar topologies, which were similar to those created using other molecular markers. The presence of more than one copy of the eEF1A gene with nearly identical coding regions makes it difficult to define the orthologous cDNAs. As shown by our data on O. tipulae, careful and extensive examination of intron positions in the eEF1A gene across the phylum is necessary to define their potential for use as valid phylogenetic markers.

Mitochondrial DNA analysis of field populations of Helicoverpa armigera (Lepidoptera: Noctuidae) and of its relationship to H. zea

BACKGROUND

Helicoverpa armigera and H. zea are amongst the most significant polyphagous pest lepidopteran species in the Old and New Worlds respectively. Separation of H. armigera and H. zea is difficult and is usually only achieved through morphological differences in the genitalia. They are capable of interbreeding to produce fertile offspring. The single species status of H. armigera has been doubted, due to its wide distribution and plant host range across the Old World. This study explores the global genetic diversity of H. armigera and its evolutionary relationship to H zea.

RESULTS

We obtained partial (511 bp) mitochondrial DNA (mtDNA) Cytochrome Oxidase-I (COI) sequences for 249 individuals of H. armigera sampled from Australia, Burkina Faso, Uganda, China, India and Pakistan which were associated with various host plants. Single nucleotide polymorphisms (SNPs) within the partial COI gene differentiated H. armigera populations into 33 mtDNA haplotypes. Shared haplotypes between continents, low F-statistic values and low nucleotide diversity between countries (0.0017-0.0038) suggests high mobility in this pest. Phylogenetic analysis of four major Helicoverpa pest species indicates that H. punctigera is basal to H. assulta, which is in turn basal to H. armigera and H. zea. Samples from North and South America suggest that H. zea is also a single species across its distribution. Our data reveal short genetic distances between H. armigera and H. zea which seem to have been established via a founder event from H. armigera stock at around 1.5 million years ago.

CONCLUSION

Our mitochondrial DNA sequence data supports the single species status of H. armigera across Africa, Asia and Australia. The evidence for inter-continental gene flow observed in this study is consistent with published evidence of the capacity of this species to migrate over long distances. The finding of high genetic similarity between Old World H. armigera and New World H. zea emphasises the need to consider work on both pests when building pest management strategies for either.

Characterization of the Helicoverpa assulta nucleopolyhedrovirus genome and sequence analysis of the polyhedrin gene region

A local strain of Helicoverpa assulta nucleopolyhedrovirus (HasNPV) was isolated from infected H.assulta larvae in Korea. Restriction endonuclease fragment analysis, using 4 restriction enzymes, estimated that the total genome size of HasNPV is about 138 kb. A degenerate polymerase chain reaction (PCR) primer set for the polyhedrin gene successfully amplified the partial polyhedrin gene of HasNPV. The sequencing results showed that the about 430 bp PCR product was a fragment of the corresponding polyhedrin gene. Using HasNPV partial predicted polyhedrin to probe the Southern blots, we identified the location of the polyhedrin gene within the 6 kb Eco RI, 15 kb Nco I, 20 kb Xho I, 17 kb Bgl II and 3 kb Cla I fragments, respectively. The 3 kb Cla I fragment was cloned and the nucleotide sequences of the polyhedrin coding region and its flaking regions were determined. Nucleotide sequence analysis indicated the presence of an open reading frame of 735 nucleotides which could encode 245 amino acids with a predicted molecular mass of 29 kDa. The nucleotide sequences within the coding region of HasNPV polyhedrin shared 73.7% identity with the polyhedrin gene from Autographa californica NPV but were most closely related to Helicoverpa and Heliothis species NPVs with over 99% sequence identity.

Dissecting the ancient rapid radiation of microgastrine wasp genera using additional nuclear genes

Previous estimates of a generic level phylogeny for the ubiquitous parasitoid wasp subfamily Microgastrinae (Hymenoptera) have been problematic due to short internal branches deep in the phylogeny. These short branches might be attributed to a rapid radiation among the taxa, the use of genes that are unsuitable for the levels of divergence being examined, or insufficient quantity of data. We added over 1200 nucleotides from four nuclear genes to a dataset derived from three genes to produce a dataset of over 3000 nucleotides per taxon. While the number of well-supported short branches in the phylogeny increased, we still did not obtain strong bootstrap support for every node. Parametric and nonparametric bootstrap simulations projected that an enormous, and likely unobtainable, amount of data would be required to get bootstrap support greater than 50% for every node. However, a marked increase in the number of well-supported nodes was seen when we conducted a Bayesian analysis of a combined dataset generated from morphological characters added to the seven gene dataset. Our results suggest that, in some cases, combining morphological and genetic characters may be the most practical way to increase support for short branches deep in a phylogeny.

Bilaterian phylogeny based on analyses of a region of the sodium-potassium ATPase beta-subunit gene

Molecular investigations of deep-level relationships within and among the animal phyla have been hampered by a lack of slowly evolving genes that are amenable to study by molecular systematists. To provide new data for use in deep-level metazoan phylogenetic studies, primers were developed to amplify a 1.3-kb region of the alpha subunit of the nuclear-encoded sodium-potassium ATPase gene from 31 bilaterians representing several phyla. Maximum parsimony, maximum likelihood, and Bayesian analyses of these sequences (combined with ATPase sequences for 23 taxa downloaded from GenBank) yield congruent trees that corroborate recent findings based on analyses of other data sets (e.g., the 18S ribosomal RNA gene). The ATPase-based trees support monophyly for several clades (including Lophotrochozoa, a form of Ecdysozoa, Vertebrata, Mollusca, Bivalvia, Gastropoda, Arachnida, Hexapoda, Coleoptera, and Diptera) but do not support monophyly for Deuterostomia, Arthropoda, or Nemertea. Parametric bootstrapping tests reject monophyly for Arthropoda and Nemertea but are unable to reject deuterostome monophyly. Overall, the sodium-potassium ATPase alpha-subunit gene appears to be useful for deep-level studies of metazoan phylogeny.

The phylogeny of the social wasp subfamily Polistinae: evidence from microsatellite flanking sequences, mitochondrial COI sequence, and morphological characters

BACKGROUND

Social wasps in the subfamily Polistinae (Hymenoptera: Vespidae) have been important in studies of the evolution of sociality, kin selection, and within colony conflicts of interest. These studies have generally been conducted within species, because a resolved phylogeny among species is lacking. We used nuclear DNA microsatellite flanking sequences, mitochondrial COI sequence, and morphological characters to generate a phylogeny for the Polistinae (Hymenoptera) using 69 species.

RESULTS

Our phylogeny is largely concordant with previous phylogenies at higher levels, and is more resolved at the species level. Our results support the monophyly of the New World subgenera of Polistini, while the Old World subgenera are a paraphyletic group. All genera for which we had more than one exemplar were supported as monophyletic except Polybia which is not resolved, and may be paraphyletic.

CONCLUSION

The combination of DNA sequences from flanks of microsatellite repeats with mtCOI sequences and morphological characters proved to be useful characters establishing relationships among the different subgenera and species of the Polistini. This is the first detailed hypothesis for the species of this important group.

Defining a synthetic pheromone blend attractive to male Heliothis subflexa under wind tunnel conditions

Heliothis subflexa males were flown in a wind tunnel to a variety of combinations of synthetic pheromone components admixed (on a filter paper disk. Blends containing (Z)-l1-hexadecenal (Z11-16:Ald, 1000 ng), (Z)-9-hexadecenal (Z9-16:Ald, 500 ng) and (Z)-11-hexadecenol (Z11-16:OH, 10-500 ng) elicited upwind flight and source contact in 52-69% of males. All these compounds have previously been isolated and identified from female H. subflexa gland extracts and volatile pheromone emissions. Males were not attracted by blends in which Z9-16:Ald was omitted (0% source contact). Similarly, blends lacking Z11-16:OH were unattractive to male H. subflexa (39% or less source contact). Males were extremely sensitive to the presence of Z11-16:OH; ever, responding in high numbers (57-69% source contact) to blends containing a dosage of 1% (10 ng) or greater Z11-16:OH. Males were unresponsive to blends in which Z9-16:Ald was replaced with a variety of dosages of (Z)-9tetradecenal, a secondary component of a closely-related congeneric species, Heliothis virescens. Another compound present in the blend emitted by con-specific females, (Z)-11-hexadecenyl acetate (Z11-16:Ac), did not inhibit H. subflexa males (69% source contact) when added to the three-component mixture (1:0.5:0.1) at a ratio of 0.1 (100 ng) with respect to Z11-16:Ald. These results indicate that Z9-16:Ald and Z11-16:OH are required in addition to Z11-16:Ald to elicit significant levels of upwind flight in H. subflexa males. The effects of Z11-16:Ac are more subtle, but at the dosage tested in these experiments, this compound does not have an antagonistic effect on upwind flight and source location by H. subflexa males.

Multiple evolutionary origins of the fungus causing Panama disease of banana: concordant evidence from nuclear and mitochondrial gene genealogies

Panama disease of banana, caused by the fungus Fusarium oxysporum f. sp. cubense, is a serious constraint both to the commercial production of banana and cultivation for subsistence agriculture. Previous work has indicated that F. oxysporum f. sp. cubense consists of several clonal lineages that may be genetically distant. In this study we tested whether lineages of the Panama disease pathogen have a monophyletic origin by comparing DNA sequences of nuclear and mitochondrial genes. DNA sequences were obtained for translation elongation factor 1alpha and the mitochondrial small subunit ribosomal RNA genes for F. oxysporum strains from banana, pathogenic strains from other hosts and putatively nonpathogenic isolates of F. oxysporum. Cladograms for the two genes were highly concordant and a partition-homogeneity test indicated the two datasets could be combined. The tree inferred from the combined dataset resolved five lineages corresponding to "F. oxysporum f. sp. cubense" with a large dichotomy between two taxa represented by strains most commonly isolated from bananas with Panama disease. The results also demonstrate that the latter two taxa have significantly different chromosome numbers. F. oxysporum isolates collected as nonpathogenic or pathogenic to other hosts that have very similar or identical elongation factor 1alpha and mitochondrial small subunit genotypes as banana pathogens were shown to cause little or no disease on banana. Taken together, these results indicate Panama disease of banana is caused by fungi with independent evolutionary origins.