MITONUCLEAR DISCORDANCE IS CAUSED BY RAMPANT MITOCHONDRIAL INTROGRESSION IN NEODIPRION (HYMENOPTERA: DIPRIONIDAE) SAWFLIES

Whole Genomes Reveal Evolutionary Relationships and Mechanisms Underlying Gene-Tree Discordance in Neodiprion Sawflies

Rapidly evolving taxa are excellent models for understanding the mechanisms that give rise to biodiversity. However, developing an accurate historical framework for comparative analysis of such lineages remains a challenge due to ubiquitous incomplete lineage sorting (ILS) and introgression. Here, we use a whole-genome alignment, multiple locus-sampling strategies, and summary-tree and single nucleotide polymorphism-based species-tree methods to infer a species tree for eastern North American Neodiprion species, a clade of pine-feeding sawflies (Order: Hymenopteran; Family: Diprionidae). We recovered a well-supported species tree that-except for three uncertain relationships-was robust to different strategies for analyzing whole-genome data. Nevertheless, underlying gene-tree discordance was high. To understand this genealogical variation, we used multiple linear regression to model site concordance factors estimated in 50-kb windows as a function of several genomic predictor variables. We found that site concordance factors tended to be higher in regions of the genome with more parsimony-informative sites, fewer singletons, less missing data, lower GC content, more genes, lower recombination rates, and lower D-statistics (less introgression). Together, these results suggest that ILS, introgression, and genotyping error all shape the genomic landscape of gene-tree discordance in Neodiprion. More generally, our findings demonstrate how combining phylogenomic analysis with knowledge of local genomic features can reveal mechanisms that produce topological heterogeneity across genomes.

Sympatric diversity pattern driven by the secondary contact of two deeply divergent lineages of the soybean pod borer Leguminivora glycinivorella

The soybean pod borer, Leguminivora glycinivorella (Matsumura), is an important tortricid pest species widely distributed in most parts of China and its adjacent regions. Here, we analyzed the genetic diversity and population differentiation of L. glycinivorella using diverse genetic information including the standard cox1 barcode sequences, mitochondrial genomes (mitogenomes), and single-nucleotide polymorphisms (SNPs) from genotyping-by-sequencing. Based on a comprehensive sampling (including adults or larvae of L. glycinivorella newly collected at 22 of the total 30 localities examined) that covers most of the known distribution range of this pest, analyses of 543 cox1 barcode sequences and 60 mitogenomes revealed that the traditionally recognized and widely distributed L. glycinivorella contains two sympatric and widely distributed genetic lineages (A and B) that were estimated to have diverged ∼1.14 million years ago during the middle Pleistocene. Moreover, low but statistically significant correlations were recognized between genetic differentiation and geographic or environmental distances, indicating the existence of local adaptation to some extent. Based on SNPs, phylogenetic inference, principal component analysis, fixation index, and admixture analysis all confirm the two divergent sympatric lineages. Compared with the stable demographic history of Lineage B, the expansion of Lineage A had possibly made the secondary contact of the two lineages probable, and this process may be driven by the climate fluctuation during the late Pleistocene as revealed by ecological niche modeling.

Intraspecific diversity of Myrmecophilus acervorum (Orthoptera: Myrmecophilidae) indicating an ongoing cryptic speciation

Myrmecophilus acervorum, previously considered a parthenogenetic species widely-distributed in Europe, has been observed to have both sexes in populations inhabiting the central part of the distribution range. Specimens from those heterosexual populations have been found being infected with Wolbachia. New mitochondrial data (COI and 16S markers) revealed the well-supported differentiation of M. acervorum populations inhabiting western Polesie (Poland) and southern Europe. In turn, analyses of EF1α marker support the hypothesis on the unfinished lineage sorting at the nuclear DNA level. Interestingly, we found that parthenogenetic populations inhabiting western Polesie are infected with Wolbachia belonging to supergroup A, while endosymbionts occurring in sexual populations of M. acervorum observed in Romania belong to supergroup B. Furthermore, new and potentially diagnostic characteristics in the external structures of the eyes of M. acervorum were identified. The surface of ommatidia in specimens occurring in southern Europe was smooth. In contrast, the ommatidia surface of individuals collected in Poland was visibly sculptured. To sum up, the significant genetic variability found in the present case, and the differentiating morphological character, are almost certainly effects of cryptic species being present within M. acervorum. This is indicative of ongoing speciation within the populations of this insect, and of simultaneous unfinished lineage sorting at the nuclear DNA level.

Museum Skins Enable Identification of Introgression Associated with Cytonuclear Discordance

Increased sampling of genomes and populations across closely related species has revealed that levels of genetic exchange during and after speciation are higher than previously thought. One obvious manifestation of such exchange is strong cytonuclear discordance, where the divergence in mitochondrial DNA (mtDNA) differs from that for nuclear genes more (or less) than expected from differences between mtDNA and nuclear DNA (nDNA) in population size and mutation rate. Given genome-scale data sets and coalescent modeling, we can now confidently identify cases of strong discordance and test specifically for historical or recent introgression as the cause. Using population sampling, combining exon capture data from historical museum specimens and recently collected tissues we showcase how genomic tools can resolve complex evolutionary histories in the brachyotis group of rock-wallabies (Petrogale). In particular, applying population and phylogenomic approaches we can assess the role of demographic processes in driving complex evolutionary patterns and assess a role of ancient introgression and hybridization. We find that described species are well supported as monophyletic taxa for nDNA genes, but not for mtDNA, with cytonuclear discordance involving at least 4 operational taxonomic units across 4 species which diverged 183-278 kya. ABC modeling of nDNA gene trees supports introgression during or after speciation for some taxon pairs with cytonuclear discordance. Given substantial differences in body size between the species involved, this evidence for gene flow is surprising. Heterogenous patterns of introgression were identified but do not appear to be associated with chromosome differences between species. These and previous results suggest that dynamic past climates across the monsoonal tropics could have promoted reticulation among related species.

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Chromosome-Aware Phylogenomics of Assassin Bugs (Hemiptera: Reduvioidea) Elucidates Ancient Gene Conflict

Though the phylogenetic signal of loci on sex chromosomes can differ from those on autosomes, chromosomal-level genome assemblies for nonvertebrates are still relatively scarce and conservation of chromosomal gene content across deep phylogenetic scales has therefore remained largely unexplored. We here assemble a uniquely large and diverse set of samples (17 anchored hybrid enrichment, 24 RNA-seq, and 70 whole-genome sequencing samples of variable depth) for the medically important assassin bugs (Reduvioidea). We assess the performance of genes based on multiple features (e.g., nucleotide vs. amino acid, nuclear vs. mitochondrial, and autosomal vs. X chromosomal) and employ different methods (concatenation and coalescence analyses) to reconstruct the unresolved phylogeny of this diverse (∼7,000 spp.) and old (>180 Ma) group. Our results show that genes on the X chromosome are more likely to have discordant phylogenies than those on autosomes. We find that the X chromosome conflict is driven by high gene substitution rates that impact the accuracy of phylogenetic inference. However, gene tree clustering showed strong conflict even after discounting variable third codon positions. Alternative topologies were not particularly enriched for sex chromosome loci, but spread across the genome. We conclude that binning genes to autosomal or sex chromosomes may result in a more accurate picture of the complex evolutionary history of a clade.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

A Test of Haldane's Rule in Neodiprion Sawflies and Implications for the Evolution of Postzygotic Isolation in Haplodiploids

AbstractHaldane's rule-a pattern in which hybrid sterility or inviability is observed in the heterogametic sex of an interspecific cross-is one of the most widely obeyed rules in nature. Because inheritance patterns are similar for sex chromosomes and haplodiploid genomes, Haldane's rule may apply to haplodiploid taxa, predicting that haploid male hybrids will evolve sterility or inviability before diploid female hybrids. However, there are several genetic and evolutionary mechanisms that may reduce the tendency of haplodiploids to obey Haldane's rule. Currently, there are insufficient data from haplodiploids to determine how frequently they adhere to Haldane's rule. To help fill this gap, we crossed a pair of haplodiploid hymenopteran species (Neodiprion lecontei and Neodiprion pinetum) and evaluated the viability and fertility of female and male hybrids. Despite considerable divergence, we found no evidence of reduced fertility in hybrids of either sex, consistent with the hypothesis that hybrid sterility evolves slowly in haplodiploids. For viability, we found a pattern opposite to that of Haldane's rule: hybrid females, but not males, had reduced viability. This reduction was most pronounced in one direction of the cross, possibly due to a cytoplasmic-nuclear incompatibility. We also found evidence of extrinsic postzygotic isolation in hybrids of both sexes, raising the possibility that this form or reproductive isolation tends to emerge early in speciation in host-specialized insects. Our work emphasizes the need for more studies on reproductive isolation in haplodiploids, which are abundant in nature but underrepresented in the speciation literature.

Evaluating species boundaries using coalescent delimitation in pine-killing Monochamus (Coleoptera: Cerambycidae) sawyer beetles

Plant-feeding beetle species are diverse and often individually highly variable. Accurate classifications can be difficult to establish yet are essential for study of evolutionary patterns and processes. Molecular data are key to further characterizing morphologically difficult groups and defining genus and species boundaries. Monochamus Dejean species are ecologically and economically significant, and in coniferous forests they vector the nematode that causes Pine Wilt Disease. This study uses nuclear and mitochondrial genes to test the monophyly and relationships of Monochamus and applies coalescent methods to further delimit the conifer-feeding species. Monochamus has also included approximately 120 Old World species associated with diverse angiosperm tree species. We sample from these additional morphologically diverse species to determine their placement in the Lamiini. Through supermatrix and coalescent methods, the higher-level relationships of Monochamus show that conifer-feeders are a monophyletic group that includes the type species and has split into Nearctic and Palearctic clades. Molecular dating indicates a single dispersal of conifer-feeders to North America over the second Bering Land Bridge circa 5.3 Ma. All other Monochamus sampled fall in different parts of the Lamiini tree. Small-bodied angiosperm-feeding Monochamus group with the monotypic genus Microgoes Casey. The African Monochamus subgenera sampled are distantly related to the conifer-feeding clade. The multispecies coalescent delimitation methods BPP and STACEY delimit 17 conifer-feeding Monochamus species for a total of 18 species, and supports the retention of all current species. An interrogation with nuclear gene allele phasing reveals that unphased data can be unreliable for accurate delimitations and divergence times. The delimited species are discussed with integrative evidence, highlighting real-world challenges in recognizing the completion of speciation trajectories.

Genomic evidence for contrasting patterns of host-associated genetic differentiation across shared host-plant species in leaf- and bud-galling sawflies

Resource specialization and ecological speciation arising through host-associated genetic differentiation (HAD) are frequently invoked as an explanation for the high diversity of plant-feeding insects and other organisms with a parasitic lifestyle. While genetic studies have demonstrated numerous examples of HAD in insect herbivores, the rarity of comparative studies means that we still lack an understanding of how deterministic HAD is, and whether patterns of host shifts can be predicted over evolutionary timescales. We applied genome-wide single nucleotide polymorphism and mitochondrial DNA sequence data obtained through genome resequencing to define species limits and to compare host-plant use in population samples of leaf- and bud-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae) collected from seven shared willow (Salicaceae: Salix) host species. To infer the repeatability of long-term cophylogenetic patterns, we also contrasted the phylogenies of the two galler groups with each other as well as with the phylogeny of their Salix hosts estimated based on RADseq data. We found clear evidence for host specialization and HAD in both of the focal galler groups, but also that leaf gallers are more specialized to single host species compared with most bud gallers. In contrast to bud gallers, leaf gallers also exhibited statistically significant cophylogenetic signal with their Salix hosts. The observed discordant patterns of resource specialization and host shifts in two related galler groups that have radiated in parallel across a shared resource base indicate a lack of evolutionary repeatability in the focal system, and suggest that short- and long-term host use and ecological diversification in plant-feeding insects are dominated by stochasticity and/or lineage-specific effects.

Ceriodaphnia (Cladocera: Daphniidae) in China: Lineage diversity, phylogeography and possible interspecific hybridization

The distribution and species/lineage diversity of freshwater invertebrate zooplankton remains understudied in China. Here, we explored the species/lineage diversity and phylogeography of Ceriodaphnia species across China. The taxonomy of this genus is under-explored. Seven morphospecies of Ceriodaphnia (C. cornuta, C. laticaudata, C. megops, C. pulchella, C. quadrangula, C. rotunda and C. spinata) were identified across 45 of 422 water bodies examined. Rather little morphological variation was observed within any single morphospecies regardless of country of origin. Nevertheless, we recognized that some or all of these morphospecies might represent species complexes. To investigate this, phylogenetic relationships within and among these morphospecies were investigated based on mitochondrial (partial cytochrome c oxidase subunit I gene) and nuclear (partial 28S rRNA gene) markers. The mitochondrial marker placed these populations in nine lineages corresponding to the morphospecies: C. laticaudata and C. pulchella were each represented by two lineages, suggesting that both are species complexes. The remaining five morphospecies were each represented by a single mtDNA lineage. Three of the nine mitochondrial lineages (belonging to C. pulchella, C. rotunda and C. megops) are newly reported and exhibited a restricted distribution within China. The nuclear-DNA phylogeny also recognized seven Ceriodaphnia taxa within China. We detected occasional mito-nuclear discordances in Ceriodaphnia taxa across China, suggesting interspecific introgression and hybridization. Our study contributes to an understanding of the species/lineage diversity of Ceriodaphnia, a genus with understudied taxonomy.

Climatic Suitability and Distribution Overlap of Sawflies (Hymenoptera: Diprionidae) and Threatened Populations of Pinaceae

Ecological Niche Models (ENM) are tools used to predict suitability, based on climatic variables selected and occurrence data of the target species, and characterize the environmental space. Sawflies (Hymenoptera: Diprionidae) are one of the main factors threatening forest health in Mexico, with cyclical population outbreaks and a wide range of hosts. In the present paper, we calculate the climatic niche in Mexico of three diprionids, Neodiprion abietis (Harris), N. omosus Smith, and Zadiprion rohweri (Middleton); the first and the latter with recent records in Mexico, and N. omosus with presence in Mexico and Guatemala; contrasting them against the distribution records of host species in the country. The climatic suitability of N. abietis was high in the Sierra Madre Occidental where its hosts, Pinus ponderosa, P. strobiformis, and P. menziesii are distributed. For N. omosus, the environmental suitability was projected towards the Southern Altiplano, where it coincides with a small presence of its hosts P. leiophylla and P. ayacahuite; however, it was possible to calculate its coincidence with more hosts in other biogeographic provinces. Pinaceae species considered under threat, Abies concolor, P. monophylla, and P. strobiformis, have populations within the environmental suitability of the sawflies.

Partial Ribosomal Nontranscribed Spacer Sequences Distinguish Rhagoletis zephyria (Diptera: Tephritidae) From the Apple Maggot, R. pomonella

The apple maggot, Rhagoletis pomonella (Walsh), was introduced into the apple-growing regions of the Pacific Northwest in the U.S.A. during the past 60-100 yr. Apple maggot (larvae, puparia, and adults) is difficult to distinguish from its morphologically similar sister species, Rhagoletis zephyria Snow, which is native and abundant in the Pacific Northwest. While morphological identifications are common practice, a simple, inexpensive assay based on genetic differences would be very useful when morphological traits are unclear. Here we report nucleotide substitution and insertion-deletion mutations in the nontranscribed spacer (NTS) of the ribosomal RNA gene cistron of R. pomonella and R. zephyria that appear to be diagnostic for these two fly species. Insertion-deletion variation is substantial and results in a 49 base-pair difference in PCR amplicon size between R. zephyria and R. pomonella that can be scored using agarose gel electrophoresis. PCR amplification and DNA sequencing of 766 bp of the NTS region from 38 R. pomonella individuals and 35 R. zephyria individuals from across their geographic ranges led to the expected PCR fragments of approx. 840 bp and 790 bp, respectively, as did amplification and sequencing of a smaller set of 26 R. pomonella and 16 R. zephyria flies from a sympatric site in Washington State. Conversely, 633 bp mitochondrial COI barcode sequences from this set of flies were polyphyletic with respect to R. pomonella and R. zephyria. Thus, differences in NTS PCR products on agarose gels potentially provide a simple way to distinguish between R. pomonella and R. zephyria.

Faster-haplodiploid evolution under divergence-with-gene-flow: simulations and empirical data from pine-feeding hymenopterans

Although haplodiploidy is widespread in nature, the evolutionary consequences of this mode of reproduction are not well characterized. Here, we examine how genome-wide hemizygosity and a lack of recombination in haploid males affects genomic differentiation in populations that diverge via natural selection while experiencing gene flow. First, we simulated diploid and haplodiploid "genomes" (500-kb loci) evolving under an isolation-with-migration model with mutation, drift, selection, migration, and recombination; and examined differentiation at neutral sites both tightly and loosely linked to a divergently selected site. So long as there is divergent selection and migration, sex-limited hemizygosity and recombination cause elevated differentiation (i.e., produce a "faster-haplodiploid effect") in haplodiploid populations relative to otherwise equivalent diploid populations, for both recessive and codominant mutations. Second, we used genome-wide SNP data to model divergence history and describe patterns of genomic differentiation between sympatric populations of Neodiprion lecontei and N. pinetum, a pair of pine sawfly species (order: Hymenoptera; family: Diprionidae) that are specialized on different pine hosts. These analyses support a history of continuous gene exchange throughout divergence and reveal a pattern of heterogeneous genomic differentiation that is consistent with divergent selection on many unlinked loci. Third, using simulations of haplodiploid and diploid populations evolving according to the estimated divergence history of N. lecontei and N. pinetum, we found that divergent selection would lead to higher differentiation in haplodiploids. Based on these results, we hypothesize that haplodiploids undergo divergence-with-gene-flow and sympatric speciation more readily than diploids.

Genome-wide sequence data show no evidence of hybridization and introgression among pollinator wasps associated with a community of Panamanian strangler figs

The specificity of pollinator host choice influences opportunities for reproductive isolation in their host plants. Similarly, host plants can influence opportunities for reproductive isolation in their pollinators. For example, in the fig and fig wasp mutualism, offspring of fig pollinator wasps mate inside the inflorescence that the mothers pollinate. Although often host specific, multiple fig pollinator species are sometimes associated with the same fig species, potentially enabling hybridization between wasp species. Here, we study the 19 pollinator species (Pegoscapus spp.) associated with an entire community of 16 Panamanian strangler fig species (Ficus subgenus Urostigma, section Americanae) to determine whether the previously documented history of pollinator host switching and current host sharing predicts genetic admixture among the pollinator species, as has been observed in their host figs. Specifically, we use genome‐wide ultraconserved element (UCE) loci to estimate phylogenetic relationships and test for hybridization and introgression among the pollinator species. In all cases, we recover well‐delimited pollinator species that contain high interspecific divergence. Even among pairs of pollinator species that currently reproduce within syconia of shared host fig species, we found no evidence of hybridization or introgression. This is in contrast to their host figs, where hybridization and introgression have been detected within this community, and more generally, within figs worldwide. Consistent with general patterns recovered among other obligate pollination mutualisms (e.g. yucca moths and yuccas), our results suggest that while hybridization and introgression are processes operating within the host plants, these processes are relatively unimportant within their associated insect pollinators.

Cytonuclear discordance in the crowned-sparrows, Zonotrichia atricapilla and Zonotrichia leucophrys

The golden-crowned (Zonotrichia atricapilla) and white-crowned (Z. leucophrys) sparrows have been presented as a compelling case for rapid speciation. They display divergence in song and plumage with overlap in their breeding ranges implying reproductive isolation, but have almost identical mitochondrial genomes. Previous research proposed hybridization and subsequent mitochondrial introgression as an alternate explanation, but lacked robust nuclear gene trees to distinguish between introgression and incomplete lineage sorting. We test for signatures of these processes between Z. atricapilla and Z. leucophrys, and investigate the relationships among Z. leucophrys subspecies, using mitochondrial sequencing and a reduced representation nuclear genomic dataset. Contrary to the paraphyly evident in mitochondrial gene trees, we confirmed the reciprocal monophyly of Z. atricapilla and Z. leucophrys using large panels of single nucleotide polymorphisms (SNPs). The pattern of cytonuclear discordance is consistent with limited, historical hybridization and mitochondrial introgression, rather than a recent origin and incomplete lineage sorting between recent sister species. We found evidence of nuclear phylogeographic structure within Z. leucophrys with two distinct clades. Altogether, our results indicate deeper divergences between Z. atricapilla and Z. leucophrys than inferred using mitochondrial markers. Our results demonstrate the limitations of relying solely on mitochondrial DNA for taxonomy, and raise questions about the possibility of selection on the mitochondrial genome during temperature oscillations (e.g. during the Pleistocene). Historical mitochondrial introgression facilitated by past environmental changes could cause erroneous dating of lineage splitting in other taxa when based on mitochondrial DNA alone.

How mitonuclear discordance and geographic variation have confounded species boundaries in a widely studied snake

As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.

Gene flow creates a mirage of cryptic species in a Southeast Asian spotted stream frog complex

Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.

Phylogenomic test of mitochondrial clues to archaic ancestors in a group of hybridizing swallowtail butterflies

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.

Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis)

Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.

Understanding Admixture: Haplodiploidy to the Rescue

Hybridization has broad evolutionary consequences, from fueling or counteracting speciation to facilitating adaptation to novel environments. Hybridization and subsequent introgression appear widespread along the tree of life. However, our understanding of how distinct evolutionary forces shape admixed genomes and the fate of introgressed genetic variants remains scarce. Most admixture research in animals has focused on diploid organisms. We propose that haplodiploid organisms can help resolve open questions about the genomic consequences of hybridization in natural populations. The ploidy difference between haploid males and diploid females, the availability of genome-wide male haplotypes, and ongoing cases of admixture make haplodiploid organisms promising models to improve our knowledge with regards to the evolution of hybrid genomes.

Host and geography together drive early adaptive radiation of Hawaiian planthoppers

The interactions between insects and their plant host have been implicated in driving diversification of both players. Early arguments highlighted the role of ecological opportunity, with the idea that insects "escape and radiate" on new hosts, with subsequent hypotheses focusing on the interplay between host shifting and host tracking, coupled with isolation and fusion, in generating diversity. Because it is rarely possible to capture the initial stages of diversification, it is particularly difficult to ascertain the relative roles of geographic isolation versus host shifts in initiating the process. The current study examines genetic diversity between populations and hosts within a single species of endemic Hawaiian planthopper, Nesosydne umbratica (Hemiptera, Delphacidae). Given that the species was known as a host generalist occupying unrelated hosts, Clermontia (Campanulaceae) and Pipturus (Urticaceae), we set out to determine the relative importance of geography and host in structuring populations in the early stages of differentiation on the youngest islands of the Hawaiian chain. Results from extensive exon capture data showed that N. umbratica is highly structured, both by geography, with discrete populations on each volcano, and by host plant, with parallel radiations on Clermontia and Pipturus leading to extensive co-occurrence. The marked genetic structure suggests that populations can readily become established on novel hosts provided opportunity; subsequent adaptation allows monopolization of the new host. The results support the role of geographic isolation in structuring populations and with host shifts occurring as discrete events that facilitate subsequent parallel geographic range expansion.

Reconciling the Mitonuclear Compatibility Species Concept with Rampant Mitochondrial Introgression

The mitonuclear compatibility species concept defines a species as a population that is genetically isolated from other populations by uniquely coadapted mitochondrial (mt) and nuclear genes. A key prediction of this hypothesis is that the mt genotype of each species will be functionally distinct and that introgression of mt genomes will be prevented by mitonuclear incompatibilities that arise when heterospecific mt and nuclear genes attempt to cofunction to enable aerobic respiration. It has been proposed, therefore, that the observation of rampant introgression of mt genotypes from one species to another constitutes a strong refutation of the mitonuclear speciation. The displacement of a mt genotype from a nuclear background with which it co-evolved to a foreign nuclear background will necessarily lead to fitness loss due to mitonuclear incompatibilities. Here I consider two potential benefits of mt introgression between species that may, in some cases, overcome fitness losses arising from mitonuclear incompatibilities. First, the introgressed mt genotype may be better adapted to the local environment than the native mt genotype such that higher fitness is achieved through improved adaptation via introgression. Second, if the mitochondria of the recipient taxa carry a high mutational load, then introgression of a foreign, less corrupt mt genome may enable the recipient taxa to escape its mutational load and gain a fitness advantage. Under both scenarios, fitness gains from novel mt genotypes could theoretically compensate for the fitness that is lost via mitonuclear incompatibility. I also consider the role of endosymbionts in non-adaptive rampant introgression of mt genomes. I conclude that rampant introgression is not necessarily evidence against the idea of tight mitonuclear coadaptation or the mitonuclear compatibility species concept. Rampant mt introgression will typically lead to erasure of species but in some cases could lead to hybrid speciation.

Mitogenomics of historical type specimens of Australasian turtles: clarification of taxonomic confusion and old mitochondrial introgression

Diagnosability is central to taxonomy as are type specimens which define taxa. New advances in technologies and the discovery of new informative traits must be matched with previous taxonomic decisions based on name-bearing type specimens. Consequently, the challenge of sequencing highly degraded DNA from historical types becomes an inevitability to resolve the very many taxonomic issues arising from, by modern standards, poor historical species descriptions leading to difficulties to assign names to genetic clusters identified from fresh material. Here we apply high-throughput parallel sequencing and sequence baiting to reconstruct the mitogenomes from 18 type specimens of Australasian side-necked turtles (Chelidae). We resolve a number of important issues that have confused the taxonomy of this family, and analyse the mitogenomes of the types and those of fresh material to improve our understanding of the phylogenetic relationships of this morphologically conservative group. Together with previously published nuclear genomic data, our study provides evidence for multiple old mitochondrial introgressions.

New lineages and old species: Lineage diversity and regional distribution of Moina (Crustacea: Cladocera) in China

The distribution and genetic diversity of freshwater zooplankton is understudied in the Eastern Palearctic. Here, we explored the lineage diversity and regional distribution of the genus Moina in China. Members of this genus are often keystone components of freshwater ecosystems and have been frequently subjected to toxicological and physiological studies. Four species of Moina were identified, based on morphology, in 50 of 113 Chinese water bodies examined, and their phylogenetic position was analyzed using both a mitochondrial (mitochondrial cytochrome c oxidase subunit I; COI) and a nuclear marker (the nuclear internal transcribed spacer; ITS-1). Both molecular markers identified four clades corresponding broadly to the morphological species. Mitochondrial DNA analysis showed the presence of four species complexes with eleven lineages across China, five of which were new. However, some lineages (and even individual haplotypes) were widespread in Eurasia, suggesting an ability to disperse over long distances. In contrast, a few lineages exhibited restricted distributions. The nuclear phylogeny also recognized four species of Moina within China and seven very distinct clades. Interestingly, one specimen possessing Moina cf. micrura mtDNA had ITS-1 alleles of the M. cf. brachiata clade. This discordance between mtDNA and nuclear ITS-1 phylogenies is indicative of interspecific introgression and hybridization. Additionally, our COI phylogeny showed apparent paraphyly in two Moina species groups, suggesting introgression of their mitochondrial genomes. Our data shows the regional distribution/diversity of the Moina species complex in a Eurasian context.

Cytonuclear diversity and shared mitochondrial haplotypes among Daphnia galeata populations separated by seven thousand kilometres

BACKGROUND

The zooplanktonic cladocerans Daphnia, present in a wide range of water bodies, are an important component of freshwater ecosystems. In contrast to their high dispersal capacity through diapausing eggs carried by waterfowl, Daphnia often exhibit strong population genetic differentiation. Here, to test for common patterns in the population genetic structure of a widespread Holarctic species, D. galeata, we genotyped two sets of populations collected from geographically distant areas: across 13 lakes in Eastern China and 14 lakes in Central Europe. The majority of these populations were genotyped at two types of markers: a mitochondrial gene (for 12S rRNA) and 15 nuclear microsatellite loci.

RESULTS

Mitochondrial DNA demonstrated relatively shallow divergence within D. galeata, with distinct haplotype compositions in the two study regions but one widely distributed haplotype shared between several of the Chinese as well as European populations. At microsatellite markers, clear separation was observed at both large (between China and Europe) and small (within Europe) geographical scales, as demonstrated by Factorial Correspondence Analyses, Bayesian assignment and a clustering method based on genetic distances. Genetic diversity was comparable between the sets of Chinese and European D. galeata populations for both types of markers. Interestingly, we observed a significant association between genetic distance and geographical distance for D. galeata populations in China but not in Europe.

CONCLUSIONS

Our results indicate relatively recent spread of D. galeata across wide expanses of the Palaearctic, with one mtDNA lineage of D. galeata successfully establishing over large distances. Despite a clear differentiation of Chinese and European D. galeata at a nuclear level, the pattern of genetic variation is nevertheless similar between both regions. Overall, our findings provide insights into the genetic population structure of a cladoceran species with extremely wide geographical range.

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build-up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double-digest restriction site-associated DNA sequencing (ddRADseq) to obtain and analyse genomic-scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.

The first mitochondrial genome of the model echinoid Lytechinus variegatus and insights into Odontophoran phylogenetics

Assembly of publically available next-generation sequence data facilitated the generation of three camarodont echinoid mitogenomes: two for the Green Urchin (Lytechinus variegatus) and one for the Red Urchin (Mesocentrotus franciscanus). The data generated are exploited in a phylogenomic analysis of the superfamily Odontophora, originally proposed for echinoids with tooth supports on the epiphyses of the jaw. The analysis highly supports this taxon and its current subdivision into three families: the Echinometridae, Toxopneustidae, and Strongylocentrotidae. The analysis furthermore implies that historical taxonomic issues between two members of the genus Strongylocentrotus (S. pallidus and S. droebachiensis) may have a genetic basis. The novel mitogenomes for the model species L. variegatus complements the draft genome available for this taxon, one of only three genome-enabled echinoid species. The assembly method applied herein, follows a divide-and-conquer approach that provides for reduced computational requirements and facilitates resolving assembly problems when processing ultra-high coverage next-generation sequence data.

Genetic Basis of Body Color and Spotting Pattern in Redheaded Pine Sawfly Larvae (Neodiprion lecontei)

Pigmentation has emerged as a premier model for understanding the genetic basis of phenotypic evolution, and a growing catalog of color loci is starting to reveal biases in the mutations, genes, and genetic architectures underlying color variation in the wild. However, existing studies have sampled a limited subset of taxa, color traits, and developmental stages. To expand the existing sample of color loci, we performed QTL mapping analyses on two types of larval pigmentation traits that vary among populations of the redheaded pine sawfly (Neodiprion lecontei): carotenoid-based yellow body color and melanin-based spotting pattern. For both traits, our QTL models explained a substantial proportion of phenotypic variation and suggested a genetic architecture that is neither monogenic nor highly polygenic. Additionally, we used our linkage map to anchor the current N. lecontei genome assembly. With these data, we identified promising candidate genes underlying (1) a loss of yellow pigmentation in populations in the mid-Atlantic/northeastern United States [C locus-associated membrane protein homologous to a mammalian HDL receptor-2 gene (Cameo2) and lipid transfer particle apolipoproteins II and I gene (apoLTP-II/I)], and (2) a pronounced reduction in black spotting in Great Lakes populations [members of the yellow gene family, tyrosine hydroxylase gene (pale), and dopamine N-acetyltransferase gene (Dat)]. Several of these genes also contribute to color variation in other wild and domesticated taxa. Overall, our findings are consistent with the hypothesis that predictable genes of large effect contribute to color evolution in nature.

Divergent clades or cryptic species? Mito-nuclear discordance in a Daphnia species complex

BACKGROUND

Genetically divergent cryptic species are frequently detected by molecular methods. These discoveries are often a byproduct of molecular barcoding studies in which fragments of a selected marker are used for species identification. Highly divergent mitochondrial lineages and putative cryptic species are even detected in intensively studied animal taxa, such as the crustacean genus Daphnia. Recently, eleven such lineages, exhibiting genetic distances comparable to levels observed among well-defined species, were recorded in the D. longispina species complex, a group that contains several key taxa of freshwater ecosystems. We tested if three of those lineages represent indeed distinct species, by analyzing patterns of variation of ten nuclear microsatellite markers in six populations.

RESULTS

We observed a discordant pattern between mitochondrial and nuclear DNA, as all individuals carrying one of the divergent mitochondrial lineages grouped at the nuclear level with widespread, well-recognized species coexisting at the same localities (Daphnia galeata, D. longispina, and D. cucullata).

CONCLUSIONS

A likely explanation for this pattern is the introgression of the mitochondrial genome of undescribed taxa into the common species, either in the distant past or after long-distance dispersal. The occurrence of highly divergent but rare mtDNA lineages in the gene pool of widespread species would suggest that hybridization and introgression in the D. longispina species complex is frequent even across large phylogenetic distances, and that discoveries of such distinct clades must be interpreted with caution. However, maintenance of ancient polymorphisms through selection is another plausible alternative that may cause the observed discordance and cannot be entirely excluded.

Gregariousness does not vary with geography, developmental stage, or group relatedness in feeding redheaded pine sawfly larvae

Aggregations are widespread across the animal kingdom, yet the underlying proximate and ultimate causes are still largely unknown. An ideal system to investigate this simple, social behavior is the pine sawfly genus Neodiprion, which is experimentally tractable and exhibits interspecific variation in larval gregariousness. To assess intraspecific variation in this trait, we characterized aggregative tendency within a single widespread species, the redheaded pine sawfly (N. lecontei). To do so, we developed a quantitative assay in which we measured interindividual distances over a 90‐min video. This assay revealed minimal behavioral differences: (1) between early‐feeding and late‐feeding larval instars, (2) among larvae derived from different latitudes, and (3) between groups composed of kin and those composed of nonkin. Together, these results suggest that, during the larval feeding period, the benefits individuals derive from aggregating outweigh the costs and that this cost‐to‐benefit ratio does not vary dramatically across space (geography) or ontogeny (developmental stage). In contrast to the feeding larvae, our assay revealed a striking reduction in gregariousness following the final larval molt in N. lecontei. We also found some intriguing interspecific variation: While N. lecontei and N. maurus feeding larvae exhibit significant aggregative tendencies, feeding N. compar larvae do not aggregate at all. These results set the stage for future work investigating the proximate and ultimate mechanisms underlying developmental and interspecific variation in larval gregariousness across Neodiprion.

Ancient and modern genome shuffling: Reticulate mito-nuclear phylogeny of four related allopatric species of Gyrodactylus von Nordmann, 1832 (Monogenea: Gyrodactylidae), ectoparasites on the Eurasian minnow Phoxinus phoxinus (L.) (Cyprinidae)

Phylogenetic analyses including four allopatric species of Gyrodactylus von Nordmann, 1832 on the Eurasian minnow Phoxinus phoxinus (L.) (Cyprinidae) revealed incongruence between the nuclear ITS1-5.8S-ITS2 and mitochondrial cox1 phylogenies due to ancient hybridisation. Gyrodactylus pannonicus Molnár, 1968 was sampled close to its type-locality, the upper reaches of River Tisza, tributary of Danube in the Black Sea Basin. Faunistic search detected three new related species with maximum composite likelihood distances in cox1 between 16.8-23.2% (tentatively 1.3 to 1.8 My of divergence). Gyrodactylus albolacustris n. sp. recorded in the White Sea Basin, eastern Baltic Basin and Mongolia was close to G. pannonicus in the nuclear ITS (divergence of 0.9%), but diverged in cox1 by 19.8%. The Mongolian isolate of G. albolacustris n. sp. diverged from the European isolates in cox1 by 8.9%, suggesting 0.7 My of isolation. The two other new species differed from G. pannonicus by >4% in ITS and some large indels in ITS1, and by >20% in cox1. Gyrodactylus danastriae n. sp. was found in River Strwiąż, a tributary of the River Dniester (Black Sea Basin) and was characterised by smaller size of anchors and by 29-41 bp dimorphic insertion in ITS1. Gyrodactylus botnicus n. sp. is considered endemic in the Baltic Basin, but was also found in the White Sea Basin as a postglacial immigrant, where it had hybridised with G. albolacustris n. sp. in spite of the high divergence in ITS (3.9%) and cox1 (22%). The discordant nuclear and mitochondrial phylogenies revealed an ancient mitochondrial introgression: G. albolacustris n. sp. was derived from a hybridisation combining proto-pannonicus ITS with proto-danastriae mitochondria, perhaps 1.3 My ago. The postglacial hybridisation of G. albolacustris n. sp. (as the donor of mtDNA^alb and ITS^alb) and G. botnicus n. sp. (donor of the ITS^bot) offered a model of shuffling of the genomic components: the process of the homogenisation and stabilisation of nuclear ITS (concerted evolution) and the lineage sorting has hardly begun.

Oviposition traits generate extrinsic postzygotic isolation between two pine sawfly species

BACKGROUND

Although empirical data indicate that ecological speciation is prevalent in nature, the relative importance of different forms of reproductive isolation and the traits generating reproductive isolation remain unclear. To address these questions, we examined a pair of ecologically divergent pine-sawfly species: while Neodiprion pinetum specializes on a thin-needled pine (Pinus strobus), N. lecontei utilizes thicker-needled pines. We hypothesized that extrinsic postzygotic isolation is generated by oviposition traits. To test this hypothesis, we assayed ovipositor morphology, oviposition behavior, and host-dependent oviposition success in both species and in F1 and backcross females.

RESULTS

Compared to N. lecontei, N. pinetum females preferred P. strobus more strongly, had smaller ovipositors, and laid fewer eggs per needle. Additionally, we observed host- and trait-dependent reductions in oviposition success in F1 and backcross females. Hybrid females that had pinetum-like host preference (P. strobus) and lecontei-like oviposition traits (morphology and egg pattern) fared especially poorly.

CONCLUSIONS

Together, these data indicate that maladaptive combinations of oviposition traits in hybrids contribute to extrinsic postzygotic isolation between N. lecontei and N. pinetum, suggesting that oviposition traits may be an important driver of divergence in phytophagous insects.

Ancient Mitochondrial Capture as Factor Promoting Mitonuclear Discordance in Freshwater Fishes: A Case Study in the Genus Squalius (Actinopterygii, Cyprinidae) in Greece

Hybridization and incomplete lineage sorting are common confounding factors in phylogeny and speciation resulting in mitonuclear disparity. Mitochondrial introgression, a particular case of hybridization, may, in extreme cases, lead to replacement of the mitochondrial genome of one species with that of another (mitochondrial capture). We investigated mitochondrial introgression involving two species of the cyprinid genus Squalius in the western Peloponnese region of Greece using molecular and morphological data. We found evidence of complete mitochondrial introgression of Squalius keadicus into two populations recognized as Squalius peloponensis from the Miras and Pamissos River basins and a divergence of mitochondrial genomes of S. keadicus from the Evrotas basin from that of the introgressed populations dating from the Pleistocene. Secondary contact among basins is a possible factor in connection of the species and the introgression event. Morphological analyses support the hypothesis of mitochondrial introgression, as S. keadicus was different from the other three populations recognized as S. peloponensis, although significant differences were found among the four populations. Isolation by geographical barriers arose during Pleistocene in the western Peloponnese were the source of the evolution of the two reciprocally monophyletic subclades found in the S. keadicus mitochondrial clade, and the morphological differences found among the four populations. Along with the lack of structure in the nuclear genome in the three populations ascribed to S. peloponensis, this suggests an incipient speciation process occurring in these Squalius species in the western Peloponnese.

Identification of sawflies and horntails (Hymenoptera, 'Symphyta') through DNA barcodes: successes and caveats

The 'Symphyta' is a paraphyletic assemblage at the base of the order Hymenoptera, comprising 14 families and about 8750 species. All have phytophagous larvae, except for the Orussidae, which are parasitoids. This study presents and evaluates the results of DNA barcoding of approximately 5360 specimens of 'Symphyta', mainly adults, and 4362 sequences covering 1037 species were deemed of suitable quality for inclusion in the analysis. All extant families are represented, except for the Anaxyelidae. The majority of species and specimens are from Europe, but approximately 38% of the species and 13% of the specimens are of non-European origin. The utility of barcoding for species identification and taxonomy of 'Symphyta' is discussed on the basis of examples from each of the included families. A significant level of cryptic species diversity was apparent in many groups. Other attractive applications include the identification of immature stages without the need to rear them, community analyses based on metabarcoding of bulk samples and association of the sexes of adults.

What haplodiploids can teach us about hybridization and speciation

Most evolutionary theory focuses on species that reproduce through sexual reproduction where both sexes have a diploid chromosome count. Yet a substantial proportion of multicellular species display complex life cycles, with both haploid and diploid life stages. A classic example is haplodiploidy, where females develop from fertilized eggs and are diploid, while males develop from unfertilized eggs and are haploid. Although haplodiploids make up about 15% of all animals (de la Filia et al. ), this type of reproduction is rarely considered in evolutionary theory. In this issue of Molecular Ecology, Patten et al. () develop a theoretical model to compare the rate of nuclear and mitochondrial introgression in haplodiploid and diploid species. They show that when two haplodiploid species hybridize, nuclear genes are much less likely to cross the species barrier than if both species were to be diploids. The reason for this is that only half of the offspring resulting from matings between haplodiploid species are true hybrids: sons from such mating only inherit their mother genes and therefore only contain genes of the maternal species. Truly, hybrid males can only occur through backcrossing of a hybrid female to a male of one of the parental species. While this twist of haplodiploid transmission genetics limits nuclear introgression, mitochondrial genes, which are maternally inherited, are unaffected by the scarcity of hybrid males. In other words, the rate of mitochondrial introgression is the same for haplodiploid and diploid species. As a result, haplodiploid species on average show a bias of mitochondrial compared to nuclear introgression.

Biased introgression of mitochondrial and nuclear genes: a comparison of diploid and haplodiploid systems

Hybridization between recently diverged species, even if infrequent, can lead to the introgression of genes from one species into another. The rates of mitochondrial and nuclear introgression often differ, with some taxa showing biases for mitochondrial introgression and others for nuclear introgression. Several hypotheses exist to explain such biases, including adaptive introgression, sex differences in dispersal rates, sex-specific prezygotic isolation and sex-specific fitness of hybrids (e.g. Haldane's rule). We derive a simple population genetic model that permits an analysis of sex-specific demographic and fitness parameters and measures the relative rates of mitochondrial and nuclear introgression between hybridizing pairs. We do this separately for diploid and haplodiploid species. For diploid taxa, we recover results consistent with previous hypotheses: an excess of one sex among the hybridizing migrants or sex-specific prezygotic isolation causes a bias for one type of marker or the other; when Haldane's rule is obeyed, we find a mitochondrial bias in XY systems and a nuclear bias in ZW systems. For haplodiploid taxa, the model reveals that owing to their unique transmission genetics, they are seemingly assured of strong mitochondrial biases in introgression rates, unlike diploid taxa, where the relative fitness of male and female hybrids can tip the bias in either direction. This heretofore overlooked aspect of hybridization in haplodiploids provides what is perhaps the most likely explanation for differential introgression of mitochondrial and nuclear markers and raises concerns about the use of mitochondrial DNA barcodes for species delimitation in these taxa.

Evolutionary consequences of climate-induced range shifts in insects

Range shifts can rapidly create new areas of geographic overlap between formerly allopatric taxa and evidence is accumulating that this can affect species persistence. We review the emerging literature on the short- and long-term consequences of these geographic range shifts. Specifically, we focus on the evolutionary consequences of novel species interactions in newly created sympatric areas by describing the potential (i) short-term processes acting on reproductive barriers between species and (ii) long-term consequences of range shifts on the stability of hybrid zones, introgression and ultimately speciation and extinction rates. Subsequently, we (iii) review the empirical literature on insects to evaluate which processes have been studied, and (iv) outline some areas that deserve increased attention in the future, namely the genomics of hybridisation and introgression, our ability to forecast range shifts and the impending threat from insect vectors and pests on biodiversity, human health and crop production. Our review shows that species interactions in de novo sympatric areas can be manifold, sometimes increasing and sometimes decreasing species diversity. A key issue that emerges is that climate-induced hybridisations in insects are much more widespread than anticipated and that rising temperatures and increased anthropogenic disturbances are accelerating the process of species mixing. The existing evidence only shows the tip of the iceberg and we are likely to see many more cases of species mixing following range shifts in the near future.

Congruence of microsatellite and mitochondrial DNA variation in acrobat ants (Crematogaster subgenus Decacrema, Formicidae: Myrmicinae) inhabiting Macaranga (Euphorbiaceae) myrmecophytes

A previously reported mitochondrial DNA (mtDNA) phylogeny of Crematogaster (subgenus Decacrema) ants inhabiting Macaranga myrmecophytes indicated that the partners diversified synchronously and their specific association has been maintained for 20 million years. However, the mtDNA clades did not exactly match morphological species, probably owing to introgressive hybridization among younger species. In this study, we determined the congruence between nuclear simple sequence repeat (SSR, also called microsatellite) genotyping and mtDNA phylogeny to confirm the suitability of the mtDNA phylogeny for inferring the evolutionary history of Decacrema ants. Analyses of ant samples from Lambir Hills National park, northeastern Borneo, showed overall congruence between the SSR and mtDNA groupings, indicating that mtDNA markers are useful for delimiting species, at least at the local level. We also found overall high host-plant specificity of the SSR genotypes of Decacrema ants, consistent with the specificity based on the mtDNA phylogeny. Further, we detected cryptic genetic assemblages exhibiting high specificity toward particular plant species within a single mtDNA clade. This finding, which may be evidence for rapid ecological and genetic differentiation following a host shift, is a new insight into the previously suggested long-term codiversification of Decacrema ants and Macaranga plants.

A comparative analysis of genetic differentiation across six shared willow host species in leaf- and bud-galling sawflies

Genetic divergence and speciation in plant-feeding insects could be driven by contrasting selection pressures imposed by different plant species and taxa. While numerous examples of host-associated differentiation (HAD) have been found, the overall importance of HAD in insect diversification remains unclear, as few studies have investigated its frequency in relation to all speciation events. One promising way to infer the prevalence and repeatability of HAD is to estimate genetic differentiation in multiple insect taxa that use the same set of hosts. To this end, we measured and compared variation in mitochondrial COI and nuclear ITS2 sequences in population samples of leaf-galling Pontania and bud-galling Euura sawflies (Hymenoptera: Tenthredinidae) collected from six Salix species in two replicate locations in northern Fennoscandia. We found evidence of frequent HAD in both species complexes, as individuals from the same willow species tended to cluster together on both mitochondrial and nuclear phylogenetic trees. Although few fixed differences among the putative species were found, hierarchical AMOVAs showed that most of the genetic variation in the samples was explained by host species rather than by sampling location. Nevertheless, the levels of HAD measured across specific pairs of host species were not correlated in the two focal galler groups. Hence, our results support the hypothesis of HAD as a central force in herbivore speciation, but also indicate that evolutionary trajectories are only weakly repeatable even in temporally overlapping radiations of related insect taxa.

Assortative mating between two sympatric closely-related specialists: inferred from molecular phylogenetic analysis and behavioral data

Host plant shifting of phytophagous insects can lead to the formation of host associated differentiation and ultimately speciation. In some cases, host plant specificity alone acts as a nearly complete pre-mating isolating barrier among insect populations. We here test whether effective pre-mating isolation and host-independent behavioral isolation have evolved under the condition of extreme host specilization using two sympatric flea beetles with incomplete post-mating isolation under laboratory conditions. Phylogenetic analysis and coalescent simulation results showed that there is a limited interspecific gene flow, indicating effctive isolation between these species. Three types of mating tests in the absence of host plant cues showed that strong host-independent behavioral isolation has evolved between them. We conclude that almost perfect assortative mating between these two extreme host specialists results from a combination of reduced encounter rates due to differential host preference and strong sexual isolation.