Search

Genomics has revolutionized our understanding of hybridization and introgression, but most of the early evidence for these processes came from studies of mitochondrial introgression. To expand these evolutionary insights from mitochondrial patterns, we evaluate phylogenetic discordance across the nuclear genomes of a hybridizing system, the Papilio machaon group of swallowtail butterflies. This species group contains three hybrid lineages (P. brevicauda, P. joanae, and P. m. kahli) that are geographically disjunct across North America and have complete fixation of a mitochondrial lineage that is otherwise primarily found in P. m. hudsonianus, a boreal subspecies of the Holarctic P. machaon. Genome-wide nuclear markers place the three hybrid lineages as a monophyletic group that is sister to P. polyxenes/P. zelicaon rather than P. machaon, although ancient hybridization between a subspecies of P. machaon and the ancestor of these three lineages is also shown by their greater nuclear affinity to P. m. hudsonianus than to other subspecies of P. machaon. Individuals from contemporary hybrid swarms in Alberta, where mitochondrial DNA fixation has not occurred, were more intermediate between their respective parent species, demonstrating diversity in mito-nuclear discordance following hybrid interactions. Our new phylogenetic findings for the P. machaon species group also include: subspecific paraphyly within P. machaon itself across its Holarctic distribution; paraphyly of P. zelicaon relative to P. polyxenes; and more divergent placement of a Mediterranean species, P. hospiton. These results provide the first comprehensive genomic evaluation of relationships within this species group and provide insight into the evolutionary dynamics of hybridization and mitochondrial introgression.

International audience; Mitogenomes are useful markers for phylogenetic studies across a range of taxonomic levels. Here, we focus on mitogenome variation across the tortricid moth genus Choristoneura and particularly the spruce budworm (Choristoneura fumiferana) species complex, a notorious pest group of North American conifer forests. Phylogenetic relationships of Tortricidae, representing two subfamilies, four tribes and nine genera, were analyzed using 21 mitogenomes. These included six newly-sequenced mitogenomes for species in the spruce budworm complex plus three additional Choristoneura species and 12 previously published mitogenomes from other tortricids and one from the Cossidae. We evaluated the phylogenetic informativeness of the mitogenomes and reconstructed a time-calibrated tree with fossil and secondary calibrations. We found that tortricid mitogenomes had conserved protein and ribosomal regions, and analysis of all protein-coding plus ribosomal genes together provided an efficient marker at any taxonomic rank. The time-calibrated phylogeny showed evolutionary convergence of conifer feeding within Choristoneura, with two independent lineages, the Nearctic spruce budworm complex and the Palearctics pecies Choristoneura murinana, both shifting onto conifers about 11 million years ago from angiosperms. These two host-plant shifts both occurred after the formation of boreal forest in the late Miocene. Haplotype diversification within the spruce budworm complex occurred in the last 4 million years, and is probably linked to the initial cooling cycles of the Northern Hemisphere in the Pliocene.

High throughput sequencing technologies have revolutionized the potential to reconcile incongruence between gene and species trees, and numerous approaches have been developed to take advantage of these advances. Genotyping-by-sequencing is becoming a regular tool for gathering phylogenetic data, yet comprehensive evaluations of phylogenetic methods using these data are sparse. Here we use multiple phylogenetic and population genetic methods for genotyping-by-sequencing data to assess species relationships in a group of forest insect pests, the spruce budworm (Choristoneura fumiferana) species complex. With few exceptions, all methods agree on the same relationships, most notably placing C. pinus as basal to the remainder of the group, rather than C. fumiferana as previously suggested. We found strong support for the monophyly of C. pinus, C. fumiferana, and C. retiniana, but more ambiguous relationships and signatures of introgression in a clade of western lineages, including C. carnana, C. lambertiana, C. occidentalis occidentalis, C. occidentalis biennis, and C. orae. This represents the most taxonomically comprehensive genomic treatment of the spruce budworm species group, which is further supported by the broad agreement among multiple methodologies.

Pine cone beetles (Conophthorus spp.) feed and kill immature cones of Pinus species, thereby reducing seed production and seriously impairing reforestation of forest ecosystems. Population variation of Conophthorus reproductive behavior has hampered the development of semiochemical control of these pests. This difficulty is compounded by a lack of taxonomic knowledge and species diagnostic characters. Researchers and managers rely, in part, on host associations and geographic locality for species identifications and these have arguable taxonomic utility. However, host use and/or geographic separation may influence Conophthorus lineage diversification. To improve Conophthorus taxonomy and understand the association of host and geography with lineage diversification, a phylogeny of 43 individuals, including all valid species and a robust sample of C. ponderosae from different hosts, is reconstructed using 785 nucleotides of the 3'-end of the mitochondrial cytochrome oxidase I gene. Thirty trees were recovered in a parsimony analysis and the strict consensus was well resolved and supported by branch support measures. Conophthorus was monophyletic but mitochondrial polyphyly was uncovered for several species. The data also suggested an underestimation of species diversity. Phylogenetically related Conophthorus lineages were significantly associated with geographic proximity but not with host, as indicated by comparisons of character optimized geographic distributions and host associations against randomized distributions of these attributes on the parsimony tree. These results suggest that geographic separation better explains the mode of Conophthorus lineage diversification than does host specialization. Based on these results, researchers and managers of Conophthorus should consider populations as potentially different evolutionary entities until species boundaries are delineated via a robust phylogenetic revision of Conophthorus.

DNA taxonomy has been proposed as a method to quickly assess diversity and species limits in highly diverse, understudied taxa. Here we use five methods for species delimitation and two genetic markers (COI and ITS2) to assess species diversity within the parasitoid genus, Ophion. We searched for compensatory base changes (CBC’s) in ITS2, and determined that they are too rare to be of practical use in delimiting species in this genus. The other four methods used both COI and ITS2, and included distance-based (threshold analysis and ABGD) and tree-based (GMYC and PTP) models. We compared the results of these analyses to each other under various parameters and tested their performance with respect to 11 Nearctic species/morphospecies and 15 described Palearctic species. We also computed barcode accumulation curves of COI sequences to assess the completeness of sampling. The species count was highly variable depending on the method and parameters used, ranging from 47 to 168 species, with more conservative estimates of 89–121 species. Despite this range, many of the Nearctic test species were fairly robust with respect to method. We concluded that while there was often good congruence between methods, GMYC and PTP were less reliant on arbitrary parameters than the other two methods and more easily applied to genetic markers other than COI. However, PTP was less successful at delimiting test species than was GMYC. All methods, as well as the barcode accumulation curves, indicate that several Palearctic species remain undescribed and that we have scarcely begun to appreciate the Nearctic diversity within this genus.

Species identifications have been historically difficult in the Choristoneura fumiferana group, an important insect pest complex. We examined the utility of simple sequence repeats (SSRs, also referred to as microsatellites) and mitochondrial DNA (mtDNA) for delimiting and identifying eight currently recognized species sampled across North America. Four of these species formed discrete clusters using SSRs, while only two species were delimited with mtDNA. There was evidence for hybridization or incomplete lineage sorting between several species pairs. An integrative approach, using both phenotypic traits and molecular markers, allowed for the discrimination of more biologically relevant species units than did the use of molecular markers alone. As species are currently identified using putatively adaptive phenotypic traits, the differences observed between recognized species and neutral SSRs or mtDNA suggests that these species (or evolutionary significant units) have diverged via natural selection in spite of some gene flow.

DNA barcoding has focused increasing attention on the use of specific regions of mitochondrial cytochrome c oxidase I and II genes (COI–COII) to diagnose and delimit species. However, our understanding of patterns of molecular evolution within these genes is limited. Here we examine patterns of nucleotide divergence in COI–COII within species and between species pairs of Lepidoptera and Diptera using a sliding window analysis. We found that: (1) locations of maximum divergence within COI–COII were highly variable among taxa surveyed in this study; (2) there was major overlap in divergence within versus between species, including within individual COI–COII profiles; (3) graphical DNA saturation analysis showed variation in percent nucleotide transitions throughout COI–COII and only limited association with levels of DNA divergence. Ultimately, no single optimally informative 600 bp location was found within the 2.3kb of COI–COII, and the DNA barcoding region was no better than other regions downstream in COI. Consequently, we recommend that researchers should maximize sequence length to increase the probability of sampling regions of high phylogenetic informativeness, and to minimize stochastic variation in estimating total divergence.

We tested the taxonomic utility of morphology and seven mitochondrial or nuclear genes in a phylogenetic reconstruction of swallowtail butterflies in the subfamily Parnassiinae. Our data included 236 morphological characters and DNA sequences for seven genes that are commonly used to infer lepidopteran relationships (COI+COII, ND5, ND1, 16S, EF-1alpha, and wg; total 5775 bp). Nuclear genes performed best for inferring phylogenies, particularly at higher taxonomic levels, while there was substantial variation in performance among mitochondrial genes. Multiple analyses of molecular data (MP, ML and Bayesian) consistently produced a tree topology different from that obtained by morphology alone. Based on molecular evidence, sister-group relationships were confirmed between the genera Hypermnestra and Parnassius, as well as between Archon and Luehdorfia, while the monophyly of the subfamily was weakly supported. We recognize three tribes within Parnassiinae, with Archon and Luehdorfia forming the tribe Luehdorfiini Tutt, 1896 [stat. rev.]. Three fossil taxa were incorporated into a molecular clock analysis with biogeographic time constraints. Based on dispersal-vicariance (DIVA) analysis, the most recent common ancestor of Parnassiinae occurred in the Iranian Plateau and Central Asia to China. Early diversification of Parnassiinae took place at the same time that India collided into Eurasia, 65-42 million years ago.

We present a molecular phylogeny for the genus Hemileuca (Saturniidae), based on 624 bp of mitochondrial cytochrome oxidase I (COI) and 932 bp of the nuclear gene elongation factor 1 alpha (EF1alpha). Combined analysis of both gene sequences increased resolution and supported most of the phylogenetic relationships suggested by separate analysis of each gene. However, a maximum parsimony (MP) model for just COI sequence from one sample of most taxa produced a phylogeny incongruent with EF1alpha and combined dataset analyses under either MP or ML models. Time of year and time of day during which adult moths fly corresponded strongly with the phylogeny. Although most Hemileuca are diurnal, ancestral Hemileuca probably were nocturnal, fall-flying insects. The two-gene molecular phylogeny suggests that wing morphology is frequently homoplastic. There was no correlation between the primary larval hostplants and phylogenetic placement of taxa. No phylogenetic pattern of specialization was evident for single hostplant families across the genus. Our results suggest that phenological behavioral characters may be more conserved than the wing morphology characters that are more commonly used to infer phylogenetic relationships in Lepidoptera. Inclusion of a molecular component in the re-evaluation of systematic data is likely to alter prior assumptions of phylogenetic relationships in groups where such potentially homoplastic characters have been used.

We used 766 bp of DNA sequence data from the mitochondrial cytochrome oxidase I gene to reconstruct a phylogeny for 39 of 43 Ips species, many of which are economically important bark beetles. The phylogeny was reconstructed using equally weighted and weighted parsimony. In both analyses, peripheral clades were well supported while internal clades were poorly supported. Phylogenetic analysis of translated amino acids produced a poorly resolved tree that was discordant with trees reconstructed with nucleotide sequence data. Two main conclusions are drawn about the monophyly of Ips and traditional systematic groups within Ips. First, Ips is monophyletic only when I. mannsfeldi, I. nobilis, and the concinnus and latidens species groups are excluded. The latidens group, I. mannsfeldi, and I. nobilis form a monophyletic group with 3 Orthotomicus species, while the concinnus group has a more basal position. Second, the majority of the species groups in the current classification for Ips are not monophyletic. European Ips species do not form a monophyletic group, contrary to common usage, and are dispersed on the phylogenetic tree among North American species. These results indicate that a formal systematic revision of Ips is needed.

Butterflies of the genus Papilio have served as the basis for numerous studies in insect physiology, genetics, and ecology. However, phylogenetic work on relationships among major lineages in the genus has been limited and inconclusive. We have sequenced 2.3 kb of DNA from the mitochondrial cytochrome oxidase I and II genes (COI and COII) for 23 Papilio taxa and two outgroups, Pachliopta neptunus and Eurytides marcellus, in order to assess the potential of these genes for use in Papilio phylogenetics and to examine patterns of gene evolution across a broad taxonomic range. Nucleotide and amino acid variation is distributed heterogeneously, both within and between genes. Structural features of the proteins are not always reliable predictors of variation. In a combined analysis, these sequences support a nearly fully resolved topology within subgenera and species groups, though higher level relationships among species groups require additional study. The most noteworthy findings are that neither Papilio alexanor nor P. xuthus belongs in the machaon group and that the subgenus Pterourus is paraphyletic with respect to the subgenus Pyrrhosticta. We leave relationships among members of the phorcas species group as a trichotomy. These two protein coding genes, particularly COI, show excellent performance in resolving relationships at the level of species and species groups among Papilionidae. We strongly endorse a similar approach for future studies aimed at these levels.