Phylogenetics and evolution of the eastern Asian-eastern North American disjunct aphid tribe, Hormaphidini (Hemiptera: Aphididae)

A chromosome-level genome assembly of the Rhus gall aphid Schlechtendalia chinensis provides insight into the endogenization of Parvovirus-like DNA sequences

The Rhus gall aphid, Schlechtendalia chinensis, feeds on its primary host plant Rhus chinensis to induce galls, which have economic importance in medicines and the food industry. Rhus gall aphids have a unique life cycle and are economically beneficial but there is huge gap in genomic information about this group of aphids. Schlechtendalia chinensis induces rich-tannin galls on its host plant and is emerging as a model organism for both commercial applications and applied research in the context of gall production by insects. Here, we generated a high-quality chromosome-level assembly for the S. chinensis genome, enabling the comparison between S. chinensis and non-galling aphids. The final genome assembly is 344.59 Mb with 91.71% of the assembled sequences anchored into 13 chromosomes. We predicted 15,013 genes, of which 14,582 (97.13%) coding genes were annotated, and 99% of the predicted genes were anchored to the 13 chromosomes. This assembly reveals the endogenization of parvovirus-related DNA sequences (PRDs) in the S. chinensis genome, which could play a role in environmental adaptations. We demonstrated the characterization and classification of cytochrome P450s in the genome assembly, which are functionally crucial for sap-feeding insects and have roles in detoxification and insecticide resistance. This genome assembly also revealed the whole genome duplication events in S. chinensis, which can be considered in comparative evolutionary analysis. Our work represents a reference genome for gall-forming aphids that could be used for comparative genomic studies between galling and non-galling aphids and provides the first insight into the endogenization of PRDs in the genome of galling aphids. It also provides novel genetic information for future research on gall-formation and insect-plant interactions.

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 novel family of secreted insect proteins linked to plant gall development

In an elaborate form of inter-species exploitation, many insects hijack plant development to induce novel plant organs called galls that provide the insect with a source of nutrition and a temporary home. Galls result from dramatic reprogramming of plant cell biology driven by insect molecules, but the roles of specific insect molecules in gall development have not yet been determined. Here, we study the aphid Hormaphis cornu, which makes distinctive "cone" galls on leaves of witch hazel Hamamelis virginiana. We found that derived genetic variants in the aphid gene determinant of gall color (dgc) are associated with strong downregulation of dgc transcription in aphid salivary glands, upregulation in galls of seven genes involved in anthocyanin synthesis, and deposition of two red anthocyanins in galls. We hypothesize that aphids inject DGC protein into galls and that this results in differential expression of a small number of plant genes. dgc is a member of a large, diverse family of novel predicted secreted proteins characterized by a pair of widely spaced cysteine-tyrosine-cysteine (CYC) residues, which we named BICYCLE proteins. bicycle genes are most strongly expressed in the salivary glands specifically of galling aphid generations, suggesting that they may regulate many aspects of gall development. bicycle genes have experienced unusually frequent diversifying selection, consistent with their potential role controlling gall development in a molecular arms race between aphids and their host plants.

A novel family of secreted insect proteins linked to plant gall development

In an elaborate form of inter-species exploitation, many insects hijack plant development to induce novel plant organs called galls that provide the insect with a source of nutrition and a temporary home. Galls result from dramatic reprogramming of plant cell biology driven by insect molecules, but the roles of specific insect molecules in gall development have not yet been determined. Here we study the aphidHormaphis cornu, which makes distinctive “cone” galls on leaves of witch hazelHamamelis virginiana. We found that derived genetic variants in the aphid genedeterminant of gall color(dgc) are associated with strong downregulation ofdgctranscription in aphid salivary glands, upregulation in galls of seven genes involved in anthocyanin synthesis, and deposition of two red anthocyanins in galls. We hypothesize that aphids inject DGC protein into galls, and that this results in differential expression of a small number of plant genes.Dgcis a member of a large, diverse family of novel predicted secreted proteins characterized by a pair of widely spaced cysteine-tyrosine-cysteine (CYC) residues, which we named BICYCLE proteins.Bicyclegenes are most strongly expressed in the salivary glands specifically of galling aphid generations, suggesting that they may regulate many aspects of gall development.Bicyclegenes have experienced unusually frequent diversifying selection, consistent with their potential role controlling gall development in a molecular arms race between aphids and their host plants.

One Sentence Summary

Aphidbicyclegenes, which encode diverse secreted proteins, contribute to plant gall development.

Morphometric and molecular discrimination of the sugarcane aphid, Melanaphis sacchari, (Zehntner, 1897) and the sorghum aphid Melanaphis sorghi (Theobald, 1904)

Melanaphis sacchari (Zehntner, 1897) and Melanaphis sorghi (Theobald, 1904) are major worldwide crop pests causing direct feeding damage on sorghum and transmitting viruses to sugarcane. It is common in the scientific literature to consider these two species as synonyms, referred to as the 'sugarcane aphid', although no formal study has validated this synonymy. In this study, based on the comparison of samples collected from their whole distribution area, we use both morphometric and molecular data to better characterize the discrimination between M. sacchari and M. sorghi. An unsupervised multivariate analysis of morphometric data clearly confirmed the separation of the two species. The best discriminating characters separating these species were length of the antenna processus terminalis relative to length of hind tibia, siphunculus or cauda. However, those criteria sometimes do not allow an unambiguous identification. Bayesian clustering based on microsatellite data delimited two clusters, which corresponded to the morphological species separation. The DNA sequencing of three nuclear and three mitochondrial regions revealed slight divergence between species. In particular, the COI barcode region proved to be uninformative for species separation because one haplotype is shared by both species. In contrast, one SNP located on the nuclear EF1-α gene was diagnostic for species separation. Based on morphological and molecular evidence, the invasive genotype damaging to sorghum in the US, Mexico and the Caribbean since 2013 is found to be M. sorghi.

Differentiation in the eastern Asian Periphyllus koelreuteriae (Hemiptera: Aphididae) species complex driven by climate and host plant

Divergent adaptation to different ecological conditions is regarded as important for speciation. For phytophagous insects, there is limited empirical evidence on species differentiation driven by climate and host plant. The recent application of molecular data and integrative taxonomic practice may improve our understanding of population divergence and speciation. Periphyllus koelreuteriae aphids feed exclusively on Koelreuteria (Sapindaceae) in temperate and subtropical regions of eastern Asia, and show morphological and phenological variations in different regions. In this study, phylogenetic and haplotype network analyses based on four genes revealed that P. koelreuteriae populations comprised three distinct genetic clades corresponding to climate and host plants, with the populations from subtropical highland regions and on Koelreuteria bipinnata host plants representing the most basal clade. These genetic lineages also showed distinct characteristics in terms of morphology and life cycle. The results indicate that P. koelreuteriae is a species complex with previously unrevealed lineages, whose differentiation may have been driven by climatic difference and host plant.

Diversity of bacteria associated with Hormaphidinae aphids (Hemiptera: Aphididae)

Bacteria are ubiquitous inhabitants of animals. Hormaphidinae is a particular aphid group exhibiting very diverse life history traits. However, the microbiota in this group is poorly known. In the present study, using high-throughput sequencing of bacterial 16S ribosomal RNA gene amplicons, we surveyed the bacterial flora in hormaphidine aphids and explored whether the aphid tribe, host plant and geographical distribution are associated with the distribution of secondary symbionts. The most dominant bacteria detected in hormaphidine species are heritable symbionts. As expected, the primary endosymbiont Buchnera aphidicola is the most abundant symbiont across all species and has cospeciated with its host aphids. Six secondary symbionts were detected in Hormaphidinae. Arsenophonus is widespread in Hormaphidinae species, suggesting the possibility of ancient acquisition of this symbiont. Ordination analyses and statistical tests show that the symbiont composition does not seem to relate to any of the aphid tribes, host plants or geographical distributions, which indicate that horizontal transfers might occur for these symbionts in Hormaphidinae. Correlation analysis exhibits negative interference between Buchnera and coexisting secondary symbionts, while the interactions between different secondary symbionts are complicated. These findings display a comprehensive picture of the microbiota in Hormaphidinae and may be helpful in understanding the symbiont diversity within a group of aphids.

The evolution of climate tolerance in conifer-feeding aphids in relation to their host's climatic niche

Climate adaptation has major consequences in the evolution and ecology of all living organisms. Though phytophagous insects are an important component of Earth's biodiversity, there are few studies investigating the evolution of their climatic preferences. This lack of research is probably because their evolutionary ecology is thought to be primarily driven by their interactions with their host plants. Here, we use a robust phylogenetic framework and species-level distribution data for the conifer-feeding aphid genus Cinara to investigate the role of climatic adaptation in the diversity and distribution patterns of these host-specialized insects. Insect climate niches were reconstructed at a macroevolutionary scale, highlighting that climate niche tolerance is evolutionarily labile, with closely related species exhibiting strong climatic disparities. This result may suggest repeated climate niche differentiation during the evolutionary diversification of Cinara. Alternatively, it may merely reflect the use of host plants that occur in disparate climatic zones, and thus, in reality the aphid species' fundamental climate niches may actually be similar but broad. Comparisons of the aphids' current climate niches with those of their hosts show that most Cinara species occupy the full range of the climatic tolerance exhibited by their set of host plants, corroborating the hypothesis that the observed disparity in Cinara species' climate niches can simply mirror that of their hosts. However, 29% of the studied species only occupy a subset of their hosts' climatic zone, suggesting that some aphid species do indeed have their own climatic limitations. Our results suggest that in host-specialized phytophagous insects, host associations cannot always adequately describe insect niches and abiotic factors must be taken into account.

Molecular and morphological evidence for the identity of two nominal species of Astegopteryx (Hemiptera, Aphididae, Hormaphidinae)

The morphology of many insect species is usually influenced by environmental factors and therefore high phenotypic variation exists even within a species. This causes difficulty and uncertainty in species taxonomy, which can be remedied by using molecular data and integrative taxonomy. Astegopteryxbambusae and A.bambucifoliae are currently regarded as two closely related aphid species with similar bamboo hosts and overlapping distributions in the oriental region. However, in practice it is hard to distinguish between them. By incorporating molecular data from four mitochondrial and nuclear genes as well as morphological information from an extensive collection of live specimens, the present study indicates that A.bambucifoliae is a junior synonym of A.bambusae. The data also indicate that large-scale geographic patterns of population differentiation may exist within this species.

Congruent phylogenetic relationships of Melaphidina aphids (Aphididae: Eriosomatinae: Fordini) according to nuclear and mitochondrial DNA data with taxonomic implications on generic limits

Melaphidina aphids (Rhus-gall aphids; Eriosomatinae: Fordini) comprise five genera from eastern Asia and one monotypic genus from eastern North America. Melaphidina are unique in feeding on plant species of Rhus subgenus Rhus (Anacardiaceae), on which they form galls during the summer. The phylogenetic relationships among some species of Melaphidina aphids remain controversial. In this study, we sought to resolve the backbone phylogeny of Melaphidina aphids by sampling 15 accessions representing all six genera, all species, and all subspecies except Meitanaphis microgallis using 20 gene regions: five nuclear genes as well as 13 protein-coding genes and two rRNA genes of the mitochondrial genome. Phylogenetic analyses included Bayesian and maximum likelihood methods. Independent analyses of nuclear and mitochondrial genes returned congruent topologies, and analyses of all gene regions combined showed well-supported relationships among Melaphidina species. In particular, these were: (1) Nurudea (excluding N. ibofushi) is sister to a clade composed of the five remaining genera; (2) the monotypic North American genus Melaphis is sister to a clade comprising the four remaining genera; and (3) (Schlechtendalia + N. ibofushi) is sister to the clade (Floraphis (Meitanaphis + Kaburagia). Our results support the transfer of Meitanaphis flavogallis to Kaburagia as an additional subspecies or species, and the recognition of Floraphis as a distinct genus. This study provides important molecular resources for subsequent evolutionary studies using more nuclear genes on the Melaphidina aphids and their close allies.

Islands conserve high species richness and areas of endemism of Hormaphidinae aphids

Patterns of biodiversity and endemism are important and form the foundation for biogeography and conservation studies. Hormaphidinae is an aphid group mainly distributed in the Sino-Japanese, Oriental, and Nearctic zoogeographic realms. To infer the areas of endemism of Hormaphidinae aphids in the Sino-Japanese and Oriental realms, we employed a geographical distribution dataset covering all 225 species in subfamily Hormaphidinae. In total, 1,245 distribution occurrence records for all species were analyzed in addition to the number of species in a certain grid cell to calculate species richness. Two methods (parsimony analysis of endemicity (PAE) and the use of an optimization criterion-NDM/VNDM) using a total of 6 grid sizes were applied to detect possible areas of endemism and to assess the probable effects of the 2 different methods and 6 grid sizes on the results of the patterns of the areas of endemism. The results revealed that islands presented most of the areas with high species richness and endemic species, particularly Japan, Taiwan Island, Java Island, the Malaysian Peninsula, southeast Himalaya, and the Hengduan Mountains. Most of these areas of endemism were located on islands or a peninsula. Islands were therefore shown to play an essential role in facilitating the formation of high species richness and endemism. Different grid sizes directly determined the regions of the areas of endemism, with small grid sizes tending to detect small and discrete areas of endemism and large grid sizes identifying continuous areas of endemism. In terms of the two methods, NDM/VNDM was found to identify more areas of endemism than PAE.

Genetic Structure of the Aphis craccivora (Hemiptera: Aphididae) From Thailand Inferred From Mitochondrial COI Gene Sequence

The cowpea aphid, Aphis craccivora Koch (Hemiptera: Aphididae), is one of the most destructive insect pests of legume plants worldwide. Although outbreaks of this pest occur annually in Thailand causing heavy damage, its genetic structure and demographic history are poorly understood. In order to determine genetic structure and genetic relationship of the geographic populations of this species, we examined sequences of mitochondrial cytochrome c oxidase subunit I (COI) gene of 51 individuals collected from 32 localities throughout Thailand. Within the sequences of these geographic populations, 32 polymorphic sites defined 17 haplotypes, ranging in sequence divergence from 0.2% (1 nucleotide) to 2.7% (16 nucleotides). A relatively high haplotype diversity but low nucleotide diversity was detected in the populations of A. craccivora, a finding that is typical for migratory species. Phylogenetic analysis revealed a weak phylogeographic structuring among the geographic populations and among the haplotypes, indicating their close relationship. Considering the distance between the sampling sites, the occurrence of identical haplotypes over wide areas is noteworthy. Moreover, the low genetic distance (FST ranging from -0.0460 to 0.3263) and high rate of per-generation female migration (Nm ranging from 1.0323 to 20.3333) suggested population exchange and gene flow between the A. craccivora populations in Thailand.

Another look at the phylogenetic relationships and intercontinental biogeography of eastern Asian - North American Rhus gall aphids (Hemiptera: Aphididae: Eriosomatinae): Evidence from mitogenome sequences via genome skimming

The Rhus gall aphids are sometimes referred to as subtribe Melaphidina (Aphididae: Eriosomatinae: Fordini) and comprise a unique group that forms galls on the primary host plants, Rhus. We examined the evolutionary relationships within the Melaphidina aphids using sequences of the complete mitochondrial genome and with samples of 11 of the 12 recognized species representing all six genera. Bayesian, maximum likelihood and parsimony analyses of the mitochondrial genome data support five well-supported clades within Melaphidina: (1) Nurudea (except N. ibofushi), (2) Schlechtendalia-Nurudea ibofushi, (3) Meitanaphis-Kaburagia, (4) Floraphis, and (5) Melaphis. Nurudea shiraii and N. yanoniella are sister to each other, but N. ibofushi is nested within Schlechtendalia. The Nurudea shiraii-N. yanoniella clade is sister to the large clade of the remaining taxa of Melaphidina aphids. The Bayesian and maximum likelihood analyses support the North American Melaphis rhois as sister to the clade of Floraphis-Kaburagia-Meitanaphis-Schlechtendalia from eastern Asia, whereas the parsimony analysis suggests Melaphis sister to Floraphis with low support (bootstrap support 38%), and the amino acid data weakly place it sister to Schlechtendalia-Nurudea ibofushi. The Melaphis position needs to be further tested with nuclear data. Meitanaphis flavogallis is sister to Kaburagia species instead of grouping with Meitanaphis elongallis. Using the Bayesian method, the North American Melaphis was estimated to have diverged from its closest Asian relatives around 64.6 (95% HPD 59.4-69.8) Ma, which is in the early Paleocene near the Cretaceous and Paleogene boundary (K/Pg boundary). At the K/Pg boundary, mass extinctions caused many types of insect-plant associations to disappear, and these extinctions may explain some of the difficulties in the phylogenetic placement of Melaphis within the analyses.

The relationships within the Chaitophorinae and Drepanosiphinae (Hemiptera, Aphididae) inferred from molecular-based phylogeny and comprehensive morphological data

The Chaitophorinae is a bionomically diverse Holarctic subfamily of Aphididae. The current classification includes two tribes: the Chaitophorini associated with deciduous trees and shrubs, and Siphini that feed on monocotyledonous plants. We present the first phylogenetic hypothesis for the subfamily, based on molecular and morphological datasets. Molecular analyses were based on the mitochondrial gene cytochrome oxidase subunit I (COI) and the nuclear gene elongation factor-1α (EF-1α). Phylogenetic inferences were obtained individually on each of genes and joined alignments using Bayesian inference (BI) and Maximum likelihood (ML). In phylogenetic trees reconstructed on the basis of nuclear and mitochondrial genes as well as a morphological dataset, the monophyly of Siphini and the genus Chaitophorus was supported. Periphyllus forms independent lineages from Chaitophorus and Siphini. Within this genus two clades comprising European and Asiatic species, respectively, were indicated. Concerning relationships within the subfamily, EF-1α and joined COI and EF-1α genes analysis strongly supports the hypothesis that Chaitophorini do not form a monophyletic clade. Periphyllus is a sister group to a clade containing Chaitophorus and Siphini. The Asiatic unit of Periphyllus also includes Trichaitophorus koyaensis. The analysis of morphological dataset under equally weighted parsimony also supports the view that Chaitophorini is an artificial taxon, as Lambersaphis pruinosae and Pseudopterocomma hughi, both traditionally included in the Chaitophorini, formed independent lineages. COI analyses support consistent groups within the subfamily, but relationships between groups are poorly resolved. These analyses were extended to include the species of closely related and phylogenetically unstudied subfamily Drepanosiphinae, which produced congruent results. Genera Drepanosiphum and Depanaphis are monophyletic and sister. The position of Yamatocallis tokyoensis differs in the molecular and morphological analyses, i.e. it is either an independent lineage (EF-1α, COI, joined COI and EF-1α genes) or is nested inside this unit (morphology). Our data also support separation of Chaitophorinae from Drepanosiphinae.

Changing Names with Changed Address: Integrated Taxonomy and Species Delimitation in the Holarctic Colymbetes paykulli Group (Coleoptera: Dytiscidae)

Species delimitation of geographically isolated forms is a long-standing problem in less studied insect groups. Often taxonomic decisions are based directly on morphologic variation, and lack a discussion regarding sample size and the efficiency of migration barriers or dispersal/migration capacity of the studied species. These problems are here exemplified in a water beetle complex from the Bering Sea region that separates North America from Eurasia. Only a few sampled specimens occur from this particular area and they are mostly found in museum and private collections. Here we utilize the theory of integrated taxonomy to discuss the speciation of the Holarctic Colymbetes paykulli water beetle complex, which historically has included up to five species of which today only two are recognized. Three delimitation methods are used; landmark based morphometry of body shape, variation in reticulation patterns of the pronotum exo-skeleton and sequence variation of the partial mitochondrial gene Cyt b. Our conclusion is that the Palearctic and Nearctic populations of C. paykulli are given the status of separate species, based on the fact that all methods showed significant separation between populations. As a consequence the name of the Palearctic species is C. paykulli Erichson and the Nearctic species should be known as C. longulus LeConte. There is no clear support for delineation between Palearctic and Nearctic populations of C. dahuricus based on mtDNA. However, significant difference in size and reticulation patterns from the two regions is shown. The combined conclusion is that the C. dahuricus complex needs a more thorough investigation to fully disentangle its taxonomic status. Therefore it is here still regarded as a Holarctic species. This study highlights the importance to study several diagnosable characters that has the potential to discriminate evolutionary lineage during speciation.

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.

Diversity of bacterial symbionts in populations of Sitobion miscanthi (Hemiptera: Aphididae) in China

Aphids are a group of insects frequently associated with bacterial symbionts. Although Chinese aphids harbor a high level of species diversity, the associations between Chinese aphids and bacterial symbionts are less known. In this study, we uncovered the diversity of bacterial symbionts in a Chinese widespread aphid, Sitobion miscanthi (Takahashi). In this study, we detected the aphid obligate symbiont Buchnera aphidicola, and two secondary symbionts, Hamiltonella defensa and Regiella insecticola, with the diagnostic polymerase chain reaction method in S. miscanthi samples. In addition, symbiotic species of Acinetobacter, Aeromonas, Enterobacter, Pantoea, and Pseudomonas, and the family Enterobacteriaceae were also found. Geographically, sporadic occurrences were detected for H. defensa and R. insecticola. Moreover, the infection rates of them vary widely among the infected populations: H. defensa (5.26-95.2%) and R. insecticola (5.26-46.7%). Phylogenetic analyses indicated that the strain of B. aphidicola mirrored the history and divergence of S. miscanthi; however, the H. defensa and R. insecticola strains were probably experienced horizontal transmission among S. miscanthi and its distantly related species.

Evolutionary relationships of Pemphigus and allied genera (Hemiptera: Aphididae: Eriosomatinae) and their primary endosymbiont, Buchnera aphidicola

Aphids harbor primary endosymbionts, Buchnera aphidicola, in specialized cells within their body cavities. Aphids and Buchnera have strict mutualistic relationships in nutrition exchange. This ancient association has received much attention from researchers who are interested in endosymbiotic evolution. Previous studies have found parallel phylogenetic relationships between non-galling aphids and Buchnera at lower taxonomic levels (genus, species). To understand whether relatively isolated habitats such as galls have effect on the parallel relationships between aphids and Buchnera, the present paper investigated the phylogenetic relationships of gall aphids from Pemphigus and allied genera, which induce pseudo-galls or galls on Populus spp. (poplar) and Buchnera. The molecular phylogenies inferred from three aphid genes (COI, COII and EF-1α) and two Buchnera genes (gnd, 16S rRNA gene) indicated significant congruence between aphids and Buchnera at generic as well as interspecific levels. Interestingly, both aphid and Buchnera phylogenies supported three main clades corresponding to the galling locations of aphids, namely leaf, the joint of leaf blade and petiole, and branch of the host plant. The results suggest phylogenetic conservatism of gall characters, which indicates gall characters are more strongly affected by aphid phylogeny, rather than host plants.

A total-evidence phylogenetic analysis of Hormaphidinae (Hemiptera: Aphididae), with comments on the evolution of galls

A phylogenetic analysis of Hormaphidinae is presented based on a total-evidence approach. Four genes (two mitochondrial, COI and CytB, and two nuclear, EF-1α and LWO) are combined with 65 morphological and seven biological characters. Sixty-three hormaphidine species representing three tribes and 36 genera as well as nine outgroups are included. Parsimony and model-based approaches are used, and several support values and implied weighting schemes are explored to assess clade stability. The monophyly of Hormaphidinae and Nipponaphidini is supported, but Cerataphidini and Hormaphidini are not recovered as monophyletic. Based on the parsimony hypothesis from the total-evidence analysis, the phylogenetic relationships within Hormaphidinae are discussed. Cerataphidini is re-delimited to exclude Doraphis and Tsugaphis, and Hormaphidini is redefined to include Doraphis. Ceratocallis Qiao & Zhang is established as a junior synonym of Ceratoglyphina van der Goot, syn. nov. Lithoaphis quercisucta Qiao, Guo & Zhang is transferred to the genus Neohormaphis Noordam as Neohormaphis quercisucta (Qiao, Guo & Zhang) comb. nov. Galls have evolved independently within three tribes of Hormaphidinae. In Cerataphidini, pseudogalls are ancestral, both single-cavity and multiple-cavity galls have evolved once, and galls appear to have evolved towards greater complexity. Galling on secondary hosts has evolved twice in hormaphidines.

Molecular phylogenetic evidence for paraphyly of Ceratovacuna and Pseudoregma (Hemiptera, Hormaphidinae) reveals late Tertiary radiation

Ceratovacuna and Pseudoregma are important groups in Cerataphidini (Hemiptera, Hormaphidinae) that not only produce soldier aphids in galls on the primary hosts but also produce horned soldiers on the herbaceous secondary hosts. However, due to sampling bias in previous studies, the phylogenetic relationships of these two genera remain inconclusive. In this study, based on more extensive sampling and examination of both mitochondrial (cytochrome c oxidase subunit I (COI); cytochrome b (Cytb)) and nuclear (elongation factor-1α (EF-1α); long-wavelength opsin (LWO)) genes, we reconstructed the phylogenetic relationships of Ceratovacuna and Pseudoregma. Phylogenetic analyses, along with morphological evidence, suggested that these two genera belong to the paraphyletic groups with species clustered into three main groups. The monophyly of Ceratovacuna and Pseudoregma as a whole was generally supported by all analyses. Monophyly of Pseudoregma was also supported. The estimated divergence times demonstrated that diversification of Ceratovacuna and Pseudoregma occurred approximately at 10 mya. The relatively low resolution for the basal relationships of the three main clades may indicate that these two genera have experienced a rapid radiation along with speciation burst of their secondary hosts during the late Tertiary.

Relationships within aphids Cinara (Cupressobium) (Hemiptera) based on mitochondrial and nuclear DNA sequences

The relationships between Cinara (Cupressobium) aphids inhabiting woody parts and leaves of conifers belonging to Cupressaceae have been studied using a mitochondrial gene (COI) and a nuclear gene (EF1-α). Based on the COI sequences, genetic distances between species ranged from 5.6 % between Cinara (C.) tujafilina (del Guercio) and Cinara (C.) juniperi (De Geer) to 10.5 % between C. (C.) tujafilina and Cinara (C.) mordvilkoi (Pašek). Genetic distances among EF1-α sequences were lower and showed from 0.1 % between C. cupressi and C. juniperi to 2.3 % between C. tujafilina and C. mordvilkoi. Molecular phylogenetic trees were constructed using the Bayesian inference (BI) phylogenetic analysis and maximum parsimony (MP) criterion. Phylogenetic trees obtained based on COI and EF1-α marker genes created two sister clades. Our results indicate that Cinara (Cupressobium) are a monophyletic group of aphids. Phylogenetic relationships amongst Cupressobium aphids do not result from the association with the host plant, but from the feeding site on the host plant or an ability to change the microhabitat on the plant. As closely related species inhabit similar microhabitats on different host plants, it suggests that the host switching is the main mode of speciation in this subgenus.

Is ecological speciation a major trend in aphids? Insights from a molecular phylogeny of the conifer-feeding genus Cinara

Introduction

In the past decade ecological speciation has been recognized as having an important role in the diversification of plant-feeding insects. Aphids are host-specialised phytophagous insects that mate on their host plants and, as such, they are prone to experience reproductive isolation linked with host plant association that could ultimately lead to species formation. The generality of such a scenario remains to be tested through macroevolutionary studies. To explore the prevalence of host-driven speciation in the diversification of the aphid genus Cinara and to investigate alternative modes of speciation, we reconstructed a phylogeny of this genus based on mitochondrial, nuclear and Buchnera aphidicola DNA sequence fragments and applied a DNA-based method of species delimitation. Using a recent software (PhyloType), we explored evolutionary transitions in host-plant genera, feeding sites and geographic distributions in the diversification of Cinara and investigated how transitions in these characters have accompanied speciation events.

Results

The diversification of Cinara has been constrained by host fidelity to conifer genera sometimes followed by sequential colonization onto different host species and by feeding-site specialisation. Nevertheless, our analyses suggest that, at the most, only half of the speciation events were accompanied by ecological niche shifts. The contribution of geographical isolation in the speciation process is clearly apparent in the occurrence of species from two continents in the same clades in relatively terminal positions in our phylogeny. Furthermore, in agreement with predictions from scenarios in which geographic isolation accounts for speciation events, geographic overlap between species increased significantly with time elapsed since their separation.

Conclusions

The history of Cinara offers a different perspective on the mode of speciation of aphids than that provided by classic models such as the pea aphid. In this genus of aphids, the role of climate and landscape history has probably been as important as host-plant specialisation in having shaped present-day diversity.

Historical biogeography of Eastern Asian-Eastern North American disjunct Melaphidina aphids (Hemiptera: Aphididae: Eriosomatinae) on Rhus hosts (Anacardiaceae)

Intercontinental biotic disjunctions have been documented and analyzed in numerous Holarctic taxa. Patterns previously synthesized for animals compared to plants suggest that the timing of animal disjunctions are mostly Early Tertiary and were generated by migration and vicariance events occurring in the North Atlantic, while plant disjunctions are mostly Mid-Late Tertiary and imply migration and vicariance over Beringia. Melaphidina aphids (Hemiptera: Aphididae: Fordini) exhibit host-alternating life cycles comprising an obligate seasonal shift between Rhus subgenus Rhus species (Anacardiaceae) and mosses (Bryophyta). Similar to their Rhus hosts, melaphidines are distributed disjunctly between Eastern Asia and Eastern North America. We examined evolutionary relationships within Melaphidina to determine the position of the North American lineage, date its divergence from Asian relatives, and compare these results to a previous historical biogeographic study of Rhus. We sampled nine species and three subspecies representing all six genera of Melaphidina. Data included sequences of mitochondrial cytochrome c oxidase subunits I and II+leucine tRNA, cytochrome b, and nuclear elongation factor 1α genes. Phylogenetic analyses (Bayesian, maximum-likelihood, parsimony) of the combined data (3282 bp) supported the monophyly of all genera except Nurudea and Schlechtendalia, due to the position of N. ibofushi. While the exact position of the North American Melaphis was not well resolved, there was high support for a derived position within Asian taxa. The divergence of Melaphis from Asian relatives centered on the Eocene-Oligocene boundary (~33-35Ma), which coincides with closure of Beringian Land Bridge I. This also corresponded to the Asian-North American disjunction previously estimated for subgenus Rhus spp. We suggest the late-Eocene Bering Land Bridge as the most likely migration route for Melaphis ancestors, as was also hypothesized for North American Rhus ancestors. Results for the Melaphidina disjunction depart from the modal pattern in animal lineages, and present a case where insect and host-plant taxa apparently responded similarly to Tertiary climate change.

Host phylogeny and specialisation in parasitoids

The host range of insect parasitoids and herbivores is influenced by both preference-related traits which mediate host choice behaviour, and performance-related traits which mediate the physiological suitability of the consumer-resource interaction. In a previous study, we characterised the influence of preference- and performance-related traits on the host range of the aphid parasitoid Binodoxys communis (Hymenoptera: Braconidae) and herein we build upon those data sets by mapping a series of these traits onto the phylogeny of the (aphid) host species. We found a strong effect of host phylogeny on overall parasitoid reproduction on the 20 host species tested, but no effect of the phylogeny of host plants of the aphids. We found an effect of aphid phylogeny on host acceptance and sting rates (related to preference) from behavioural observations and for pupal survivorship (related to performance), showing that both classes of traits show phylogenetic conservatism with respect to host species.

Host race evolution in Schizaphis graminum (Hemiptera: Aphididae): nuclear DNA sequences

The greenbug aphid, Schizaphis graminum (Rondani) was introduced into the United States in the late 1880s, and quickly was established as a pest of wheat, oat, and barley. Sorghum was also a host, but it was not until 1968 that greenbug became a serious pest of it as well. The most effective control method is the planting of resistant varieties; however, the occurrence of greenbug biotypes has hampered the development and use of plant resistance as a management technique. Until the 1990s, the evolutionary status of greenbug biotypes was obscure. Four mtDNA cytochrome oxidase subunit I (COI) haplotypes were previously identified, suggesting that S. graminum sensu lato was comprised of host-adapted races. To elucidate the current evolutionary and taxonomic status of the greenbug and its biotypes, two nuclear genes and introns were sequenced; cytochrome c (CytC) and elongation factor 1-α (EF1-α). Phylogenetic analysis of CytC sequences were in complete agreement with COI sequences and demonstrated three distinct evolutionary lineages in S. graminum. EF1-α DNA sequences were in partial agreement with COI and CytC sequences, and demonstrated two distinct evolutionary lineages. Host-adapted races in greenbug are sympatric and appear reproductively isolated. Agricultural biotypes in S. graminum likely arose by genetic recombination via meiosis during sexual reproduction within host-races. The 1968 greenbug outbreak on sorghum was the result of the introduction of a host race adapted to sorghum, and not selection by host resistance genes in crops.

Macroevolutionary patterns in the Aphidini aphids (Hemiptera: Aphididae): diversification, host association, and biogeographic origins

Background

Due to its biogeographic origins and rapid diversification, understanding the tribe Aphidini is key to understanding aphid evolution. Major questions about aphid evolution include origins of host alternation as well as age and patterns of diversification in relation to host plants. To address these questions, we reconstructed the phylogeny of the Aphidini which contains Aphis, the most diverse genus in the family. We used a combined dataset of one nuclear and four mitochondrial DNA regions. A molecular dating approach, calibrated with fossil records, was used to estimate divergence times of these taxa.

Principal Findings

Most generic divergences in Aphidini occurred in the Middle Tertiary, and species-level divergences occurred between the Middle and Late Tertiary. The ancestral state of host use for Aphidini was equivocal with respect to three states: monoecy on trees, heteroecy, and monoecy on grasses. The ancestral state of Rhopalosiphina likely included both heteroecy and monoecy, whereas that of Aphidina was most likely monoecy. The divergence times of aphid lineages at the generic or subgeneric levels are close to those of their primary hosts. The species-level divergences in aphids are consistent with the diversification of the secondary hosts, as a few examples suggest. The biogeographic origin of Aphidini as a whole was equivocal, but the major lineages within Aphidina likely separated into Nearctic, Western Palearctic, and Eastern Palearctic regions.

Conclusions

Most generic divergences in Aphidini occurred in the Middle Tertiary when primary hosts, mainly in the Rosaceae, were diverging, whereas species-level divergences were contemporaneous with diversification of the secondary hosts such as Poaceae in the Middle to Late Tertiary. Our results suggest that evolution of host alternation within Aphidini may have occurred during the Middle Tertiary (Oligocene) when the secondary hosts emerged.

Taxonomic position of Hormaphis similibetulae Qiao & Zhang, 2004 (Hemiptera, Aphididae): molecular and biological evidences

The taxonomic position of Hormaphis similibetulae Qiao & Zhang, 2004 has been reexamined. The phylogenetic position of Hormaphis similibetulae was inferred by maximum parsimony, maximum likelihood and Bayesian analyses on the basis of partial nuclear elongation factor-1α and mitochondrial tRNA leucine/cytochrome oxidase II sequences. The results showed that this species fell into the clade of Hamamelistes species, occupying a basal position, and was clearly distinct from other Hormaphis species. A closer relationship between Hormaphis similibetulae and Hamamelistes species was also revealed by life cycle analysis. Therefore, we conclude that Hormaphis similibetulae should be transferred to the genus Hamamelistes as Hamamelistes similibetulae (Qiao & Zhang), comb. n.

Use of a mitochondrial COI sequence to identify species of the subtribe Aphidina (Hemiptera, Aphididae)

Aphids of the subtribe Aphidina are found mainly in the North Temperate Zone. The relative lack of diagnostic morphological characteristics has hindered the identification of species in this group. However, DNA-based taxonomic methods can clarify species relationships within this group. Sequence variation in a partial segment of the mitochondrial COI gene was highly effective for identifying species within Aphidina. Thirty-six species of Aphidina were identified in a neighbor-joining tree. Mean intraspecific sequence divergence in Aphidina was 0.52%, with a range of 0.00% to 2.95%, and the divergences of most species were less than 1%. Mean interspecific divergence within previously recognized genera or morphologically similar species groups was 6.80%, with a range of 0.68% to 11.40%, with variation mainly in the range of 3.50% to 8.00%. Possible reasons for anomalous levels of mean nucleotide divergence within or between some taxa are discussed.

A new aphid genus Neoaulacorthum (Hemiptera: Aphididae: Macrosiphini), determined by molecular and morphometric analyses

We performed molecular and morphological analyses to determine the generic limit of the genus Aulacorthum, including several species with controversial taxonomic histories. The sequences of four mitochondrial genes, COI, COII, srRNA and lrRNA, and one nuclear gene, EF1a, implied that Aulacorthum is not monophyletic, with Aulacorthum magnoliae and Aulacorthum nipponicum forming a clade that is not sister to other currently recognized Aulacorthum species. Morphometric analysis based on 20 morphological characters also showed that A. magnoliae and A. nipponicum exhibited morphological characteristics distinct from congeneric species. Based on these results, we propose a new genus, Neoaulacorthum ge. n., for A. magnoliae and A. nipponicum.

A phylogeny and classification of the Muhlenbergiinae (Poaceae: Chloridoideae: Cynodonteae) based on plastid and nuclear DNA sequences

PREMISE OF THE STUDY

To understand the origins of C(4) grasslands, we must have a better interpretation of plant traits via phylogenetic reconstruction. Muhlenbergiinae, the largest subtribe of C(4) grasses in Mexico and the southwestern United States (with 176 species), is taxonomically poorly understood. •

METHODS

We conducted a phylogenetic analysis of 47 genera and 174 species using six plastid regions (ndhA intron, ndhF, rps16-trnK, rps16 intron, rps3, and rpl32-trnL) and the nuclear ITS 1 and 2 (ribosomal internal transcribed spacer) regions to infer evolutionary relationships and revise the classification. •

KEY RESULTS

In our analyses, Muhlenbergia (ca. 153 species) is paraphyletic, with nine genera (Aegopogon, Bealia, Blepharoneuron, Chaboissaea, Lycurus, Muhlenbergia, Pereilema, Redfieldia, Schaffnerella, and Schedonnardus) found nested within. We recognized the following five well-supported monophyletic lineages within Muhlenbergia: subg. Muhlenbergia, with species that have phosphoenolpyruvate carboxykinase-like leaf anatomy and long, scaly rhizomes; subg. Trichochloa with long-lived species that are relatively tall (up to 3 m); subg. Clomena with 3-nerved upper glumes; sect. Pseudosporobolus species with narrow panicles and plumbeous spikelets; and sect. Bealia species with lemmas with hairy margins and midveins. •

CONCLUSIONS

We propose expanding the circumscription of Muhlenbergia to include the other nine genera in this subtribe and make the following new combinations: Muhlenbergia subg. Bealia, M. diandra, M. geminiflora, M. paniculata, M. phleoides, M. subg. Pseudosporobolus (also lectotipified), M. solisii, M. tricholepis. We also propose several new names: M. ammophila, M. columbi, M. plumosa. Our phylograms suggest that Muhlenbergia originated in North America because the sister (Sohnsia filifolia and Scleropogoninae) is composed of predominantly North American species.

Gall-induction in insects: evolutionary dead-end or speciation driver?

Background

The tree of life is significantly asymmetrical - a result of differential speciation and extinction - but general causes of such asymmetry are unclear. Differences in niche partitioning are thought to be one possible general explanation. Ecological specialization might lead to increases in diversification rate or, alternatively, specialization might limit the evolutionary potential of specialist lineages and increase their extinction risk. Here we compare the diversification rates of gall-inducing and non-galling insect lineages. Compared with other insect herbivores feeding on the same host plant, gall-inducing insects feed on plant tissue that is more nutritious and less defended, and they do so in a favorable microhabitat that may also provide some protection from natural enemies. We use sister-taxon comparisons to test whether gall-inducing lineages are more host-specific than non-galling lineages, and more or less diverse than non-gallers. We evaluate the significance of diversity bipartitions under Equal Rates Markov models, and use maximum likelihood model-fitting to test for shifts in diversification rates.

Results

We find that, although gall-inducing insect groups are more host-specific than their non-galling relatives, there is no general significant increase in diversification rate in gallers. However, gallers are found at both extremes - two gall-inducing lineages are exceptionally diverse (Euurina sawflies on Salicaceae and Apiomorpha scale insects on Eucalytpus), and one gall-inducing lineage is exceptionally species-poor (Maskellia armored scales on Eucalyptus).

Conclusions

The effect of ecological specialization on diversification rates is complex in the case of gall-inducing insects, but host range may be an important factor. When a gall-inducing lineage has a host range approximate to that of its non-galling sister, the gallers are more diverse. When the non-galler clade has a much wider host range than the galler, the non-galler is also much more diverse. There are also lineage-specific effects, with gallers on the same host group exhibiting very different diversities. No single general model explains the observed pattern.

Molecular phylogeny of Pemphiginae (Hemiptera: Aphididae) inferred from nuclear gene EF-1alpha sequences

Three traditional tribes of Fordini, Pemphigini and Eriosomatini comprise Pemphiginae, and there are two subtribes in Fordini and Pemphigini, respectively. Most of the species in this subfamily live heteroecious holocyclic lives with distinct primary host specificity. The three tribes of Pemphigini (except Prociphilina), Eriosomatini and Fordini use three families of plants, Salicaceae (Populus), Ulmaceae (Ulums) and Anacardiaceae (Pistacia and Rhus), as primary hosts, respectively, and form galls on them. Therefore, the Pemphigids are well known as gall makers, and their galls can be divided into true galls and pseudo-galls in type. We performed the first molecular phylogenetic study of Pemphiginae based on molecular data (EF-1alpha sequences). Results show that Pemphiginae is probably not a monophylum, but the monophyly of Fordini is supported robustly. The monophyly of Pemphigini is not supported, and two subtribes in it, Pemphigina and Prociphilina, are suggested to be raised to tribal level, equal with Fordini and Eriosomatini. The molecular phylogenetic analysis does not show definite relationships among the four tribes of Pemphiginae, as in the previous phylogenetic study based on morphology. It seems that the four tribes radiated at nearly the same time and then evolved independently. Based on this, we can speculate that galls originated independently four times in the four tribes, and there is no evidence to support that true galls are preceded by pseudo-galls, as in the case of thrips and willow sawflies.

Molecular phylogeny of Fordini (Hemiptera: Aphididae: Pemphiginae) inferred from nuclear gene EF-1 alpha and mitochondrial gene COI

The tribe Fordini is a fascinating group because of its complicated life history, primary host specificity and gall-forming characteristic. Different species produce galls with different morphology on different parts of the host plants. The EF-1alpha-based, COI-based and combined sequences-based phylogenetic trees with three algorithms MP, ML and Bayes all strongly suggest that Fordini is a monophyletic group with two clades corresponding to two subtribes, Fordina and Melaphidina, each also monophyletic. Some important morphological characters and primary host plants of aphids were mapped onto the phylogenetic tree to analyse the division of subtribes and to uncover at which level the aphids correspond to their primary hosts, Pistacia and Rhus. Results suggest that the division of subtribes in Fordini is closely related to host selection of aphids. The evolution of gall morphology and the probable driving force behind it in this tribe were also discussed. The Fordini aphids seem to have evolved towards a better ability to manipulate their host plant, induce strong sinks and gain high reproductive success. Galls in this tribe evolved mainly along two directions to attain this goal: (i) by enlarging the gall from small bag to spherical, even big cauliflower-like, and changing the galls' location or forming two galls in their life cycle (Fordina); (ii) by moving the gall position from midrib, petiole of the leaflet, and eventually to the common petiole of the compound leaf (Melaphidina).

Assembly of the eastern North American herpetofauna: new evidence from lizards and frogs

Darwin first recognized the importance of episodic intercontinental dispersal in the establishment of worldwide biotic diversity. Faunal exchange across the Bering Land Bridge is a major example of such dispersal. Here, we demonstrate with mitochondrial DNA evidence that three independent dispersal events from Asia to North America are the source for almost all lizard taxa found in continental eastern North America. Two other dispersal events across Beringia account for observed diversity among North American ranid frogs, one of the most species-rich groups of frogs in eastern North America. The contribution of faunal elements from Asia via dispersal across Beringia is a dominant theme in the historical assembly of the eastern North American herpetofauna.

Evolution of host specialization in the Adelgidae (Insecta: Hemiptera) inferred from molecular phylogenetics

The Adelgidae form a small group of insects in the Aphidoidea. They are cyclically parthenogenetic with host alternating, multiple-generation complex life cycles and are restricted to certain host genera in the Pinaceae. Species that host alternate always have Picea as the primary host where sexual reproduction and gall formation occur, and another genus in the Pinaceae as the secondary host where a series of parthenogenetic generations are produced. Other species that do not host alternate complete their entire life cycle on one host and only reproduce parthenogenetically. We studied relationships within Adelgidae using DNA sequences from the mitochondrial COI, COII, and cytb genes, and the nuclear EF1alpha gene. Analysis of the combined data resulted in a well-resolved phylogeny in which the major adelgid clades correspond neatly to their association with secondary host genera. Specialization on each secondary host genus occurred only once and was followed by diversification on the host genus. Molecular dating of divergence times in the Adelgidae suggest that diversification among host genera occurred in the Late Cretaceous and Early Tertiary when the Pinaceae genera were diverging. It is not clear, however, whether the Adelgidae and Pinaceae co-diversified because the relationships among the Pinaceae genera are not fully resolved. We discuss implications for adelgid taxonomy, life cycle evolution, and evolution of the interaction between adelgids and their host plants.

Phylogeny of the genus Aphis Linnaeus, 1758 (Homoptera: Aphididae) inferred from mitochondrial DNA sequences

Aphis is the largest aphid genus in the world and contains several of the most injurious aphid pests. It is also the most reluctant aphid genus to any comprehensive taxonomic treatment: while most species are easily classified into "species groups" that form well defined entities, numerous species within these groups are difficult to tell apart morphologically and identification keys remain ambiguous and mostly rely on host plant affiliation. In this paper, we used partial sequences of COI/COII and CytB genes to reconstruct the first phylogeny of Aphis and discuss the present systematics. The monophyly of the subgenus Bursaphis and of the tree major species groups, Black aphid, Black backed aphid and frangulae-like species was recovered by all phylogenetic analyses. However our data suggested that the nominal subgenus was not monophyletic. Relationships between major species groups were often ambiguous but "Black" and "Black backed" species groups appeared as sister clades. The most striking result of this study was that our molecular data met the same limits as the morphological characters used in classifications: mitochondrial DNA did not allow the differentiation of species that are difficult to identify. Further, interspecies relationships within groups of species for which taxonomic treatment is difficult stayed unresolved. This suggests that species delineation in the genus Aphis is often ambiguous and that diversification might have been a rapid process.

Molecular systematics and phylogeography of the cryptic species complex Baetis rhodani (Ephemeroptera, Baetidae)

Genetic studies have highlighted cryptic diversity in many well-known taxa including aquatic insects, with the general implication that there are more species than are currently recognised. Baetis rhodani Pictet are among the most widespread, abundant and ecologically important of all European mayflies (Ephemeroptera), and used widely as biological indicators of stream quality. Traditional taxonomy and systematics have never fully resolved differences among suspected cryptic species in the B. rhodani complex because morphological characters alone do not allow reliable distinction. This is particularly true among larvae, the life-stage used most widely in monitoring. Here, we assess the molecular diversity of this complex in one of the largest such studies of cryptic species in the order Ephemeroptera to date. Phylogenies were constructed using data from the mitochondrial cytochrome oxidase subunit I (COI) gene. Two monophyletic groups were recovered consisting of one major haplogroup and a second clade of 6 smaller but distinct haplogroups. Haplogroup divergence ranged from 0.2-3% (within) to 8-19% (among) with the latter values surpassing maxima typically reported for other insects, and provides strong evidence for cryptic species in the B. rhodani complex. The taxonomic status of these seven haplogroups remains undefined. Their distributions across Western Europe reveal no obvious geographic pattern, suggesting widespread diffusion of genetic lineages since the last glacial maximum. The implications of these findings are far-reaching given the ecological and bioindicator significance of what now appears to be several taxa.