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

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini)

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

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.

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.

Genomic evidence of paternal genome elimination in the globular springtail Allacma fusca

Paternal genome elimination (PGE) - a type of reproduction in which males inherit but fail to pass on their father's genome - evolved independently in six to eight arthropod clades. Thousands of species, including several important for agriculture, reproduce via this mode of reproduction. While PGE is well established in some of the clades, the evidence in globular springtails (Symphypleona) remains elusive, even though they represent the oldest and most species rich clade putatively reproducing via PGE. We sequenced genomic DNA from whole bodies of Allacma fusca males with high fractions (>27.5%) of sperm to conclusively confirm that all the sperm carry one parental haplotype only. Although it is suggestive that the single haplotype present in sperm is maternally inherited, definitive genetic proof of the parent of origin is still needed. The genomic approach we developed allows for detection of genotypic differences between germline and soma in all species with sufficiently high fraction of germline in their bodies. This opens new opportunities for scans of reproductive modes in small organisms.

Rampant nuclear-mitochondrial-plastid phylogenomic discordance in globally distributed calcifying microalgae

Incongruent phylogenies have been widely observed between nuclear and plastid or mitochondrial genomes in terrestrial plants and animals. However, few studies have examined these patterns in microalgae or the discordance between the two organelles. Here we investigated the nuclear-mitochondrial-plastid phylogenomic incongruence in Emiliania-Gephyrocapsa, a group of cosmopolitan calcifying phytoplankton with enormous populations and recent speciations. We assembled mitochondrial and plastid genomes of 27 strains from across global oceans and temperature regimes, and analyzed the phylogenomic histories of the three compartments using concatenation and coalescence methods. Six major clades with varying morphology and distribution are well recognized in the nuclear phylogeny, but such relationships are absent in the mitochondrial and plastid phylogenies, which also differ substantially from each other. The rampant phylogenomic discordance is due to a combination of organellar capture (introgression), organellar genome recombination, and incomplete lineage sorting of ancient polymorphic organellar genomes. Hybridization can lead to replacements of whole organellar genomes without introgression of nuclear genes and the two organelles are not inherited as a single cytoplasmic unit. This study illustrates the convoluted evolution and inheritance of organellar genomes in isogamous haplodiplontic microalgae and provides a window into the phylogenomic complexity of marine unicellular eukaryotes.

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.

Avian introgression patterns are consistent with Haldane's Rule

According to Haldane’s Rule, the heterogametic sex will show the greatest fitness reduction in a hybrid cross. In birds, where sex is determined by a ZW system, female hybrids are expected to experience lower fitness compared to male hybrids. This pattern has indeed been observed in several bird groups, but it is unknown whether the generality of Haldane’s Rule also extends to the molecular level. First, given the lower fitness of female hybrids, we can expect maternally inherited loci (i.e., mitochondrial and W-linked loci) to show lower introgression rates than biparentally inherited loci (i.e., autosomal loci) in females. Second, the faster evolution of Z-linked loci compared to autosomal loci and the hemizygosity of the Z-chromosome in females might speed up the accumulation of incompatible alleles on this sex chromosome, resulting in lower introgression rates for Z-linked loci than for autosomal loci. I tested these expectations by conducting a literature review which focused on studies that directly quantified introgression rates for autosomal, sex-linked, and mitochondrial loci. Although most studies reported introgression rates in line with Haldane’s Rule, it remains important to validate these genetic patterns with estimates of hybrid fitness and supporting field observations to rule out alternative explanations. Genomic data provide exciting opportunities to obtain a more fine-grained picture of introgression rates across the genome, which can consequently be linked to ecological and behavioral observations, potentially leading to novel insights into the genetic mechanisms underpinning Haldane’s Rule.

Spatial and temporal patterns of genetic diversity in Bombus terrestris populations of the Iberian Peninsula and their conservation implications

The bumblebee Bombus terrestris is used worldwide for crop pollination. Despite its positive impact on crop yield, it has become a widespread threat to biodiversity due to its interactions with local bumblebee populations. Commercial subspecies introduced to the Iberian Peninsula since the 1990s without any regulation have colonized the environment, with evidence of naturalization and introgression with the endemic subspecies Bombus terrestris lusitanicus. We have used mitochondrial and nuclear genetic data to describe the current genetic diversity of the Iberian population and to estimate the expansion of commercial bumblebees. Samples from the natural distribution range of the commercial subspecies, the natural intergradation area between the two subspecies and from a period prior to the use of commercial colonies (i.e., before the 1990s) have been used for comparison. Our results show that the mitochondrial haplotype of the commercial breeds has spread throughout the territory, which, together with subtle changes observed in the nuclear genetic diversity of the populations, indicates that hybridization and consequent introgression are occurring in most of the peninsula. It is, therefore, necessary to improve the existing legislation concerning the management and exportation of commercial bumblebees to conserve locally adapted populations.

Limited dispersal and local adaptation promote allopatric speciation in a biodiversity hotspot

Recently diverged or diverging populations can offer unobstructed insights into early barriers to gene flow during the initial stages of speciation. The current study utilised a novel insect system (order Mantophasmatodea) to shed light on the early drivers of speciation. The members of this group have limited dispersal abilities, small allopatric distributions and strong habitat associations in the Cape Floristic Region biodiversity hotspot in South Africa. Sister taxa from the diverse family Austrophasmatidae were chosen as focal species (Karoophasma biedouwense, K. botterkloofense). Population genetics and Generalized Dissimilarity Modelling (GDM) were used to characterise spatial patterns of genetic variation and evaluate the contribution of environmental factors to population divergence and speciation. Extensive sampling confirmed the suspected allopatry of these taxa. However, hybrids were identified in a narrow region occurring between the species' distributions. Strong population structure was found over short geographic distances; particularly in K. biedouwense in which geographic distance accounted for 32% of genetic variation over a scale of 50 km (r = .56, p < .001). GDM explained 42%-78% of the deviance in observed genetic dissimilarities. Geographic distance was consistently indicated to be important for between species and within population differentiation, suggesting that limited dispersal ability may be an important neutral driver of divergence. Temperature, altitude, precipitation and vegetation were also indicated as important factors, suggesting the possible role of adaptation to local environmental conditions for species divergence. The discovery of the hybrid-zone, and the multiple allopatric species pairs in Austrophasmatidae support the idea that this could be a promising group to further our understanding of speciation modes.

Impact of invasive ant species on native fauna across similar habitats under global environmental changes

Biotic invasions can predominantly alter the dynamics, composition, functions, and structure of natural ecosystems. Social insects, particularly ants, are among the most damaging invasive alien species. Invasive ant species are among the supreme threats to ecosystems. There are about 23 species of invasive ants recorded worldwide, according to the ant invasive databases. The ecological impacts of invasive ants comprise predation, hybridization, and competition with native species that changes the ecosystem processes with the biodiversity loss and upsurge of pests. The effects of invasion on native fauna in the same habitats might be catastrophic for the native community through various ecological mechanisms, e.g., habitat disturbance, resource competition, limiting the foraging activity of native species, and various other indirect mechanisms of invasive species. Invasive species may have harmful impacts on habitats and devastating effects on natural flora and fauna, and stopping these new species from being introduced is the most effective way to deter future invasions and maintain biodiversity. This paper reviews the literature to evaluate the effects of invasive ant species on the native species, including vertebrates, invertebrates, and plants sharing the same habitats as the non-native species under global environmental changes. We also highlighted the various management strategies that could be adopted in minimizing the adverse effects of these invasive ant species on the natural ecosystem. To this end, strategies that could regulate the mode and rate of invasion by these alien ant species are the most effective ways to deter future invasions and maintain biodiversity.

Patterns of reproductive isolation in a haplodiploid mite, Amphitetranychus viennensis: prezygotic isolation, hybrid inviability and hybrid sterility

BACKGROUND

Evolution of reproductive isolation is an important process, generating biodiversity and driving speciation. To better understand this process, it is necessary to investigate factors underlying reproductive isolation through various approaches but also in various taxa. Previous studies, mainly focusing on diploid animals, supported the prevalent view that reproductive barriers evolve gradually as a by-product of genetic changes accumulated by natural selection by showing a positive relationship between the degree of reproductive isolation and genetic distance. Haplodiploid animals are expected to generate additional insight into speciation, but few studies investigated the prevalent view in haplodiploid animals. In this study, we investigate whether the relationship also holds in a haplodiploid spider mite, Amphitetranychus viennensis (Zacher).

RESULTS

We sampled seven populations of the mite in the Palaearctic region, measured their genetic distance (mtDNA) and carried out cross experiments with all combinations. We analyzed how lack of fertilization rate (as measure of prezygotic isolation) as well as hybrid inviability and hybrid sterility (as measures of postzygotic isolation) varies with genetic distance. We found that the degree of reproductive isolation varies among cross combinations, and that all three measures of reproductive isolation have a positive relationship with genetic distance. Based on the mtDNA marker, lack of fertilization rate, hybrid female inviability and hybrid female sterility were estimated to be nearly complete (99.0-99.9% barrier) at genetic distances of 0.475-0.657, 0.150-0.209 and 0.145-0.210, respectively. Besides, we found asymmetries in reproductive isolation.

CONCLUSIONS

The prevalent view on the evolution of reproductive barriers is supported in the haplodiploid spider mite we studied here. According to the estimated minimum genetic distance for total reproductive isolation in parent population crosses in this study and previous work, a genetic distance of 0.15-0.21 in mtDNA (COI) appears required for speciation in spider mites. Variations and asymmetries in the degree of reproductive isolation highlight the importance of reinforcement of prezygotic reproductive isolation through incompatibility and the importance of cytonuclear interactions for reproductive isolation in haplodiploid spider mites.

Heterogeneous taxonomic resolution of cytochrome b gene identification of bats from Argentina: Implications for field studies

Bats are among the most diverse, widespread, and abundant mammals. In Argentina, 67 species of bats have been recorded, belonging to 5 families and 29 genera. These high levels of biodiversity are likely to complicate identification at fieldwork, especially between closely related species, where external morphology-based approaches are the only immediate means for a priori species assignment. The use of molecular markers can enhance species identification, and acquires particular relevance in capture-release studies. In this study, we discuss the extent of the use of the mitochondrial cytochrome b gene for species identification, comparing external morphology identification with a molecular phylogenetic classification based on this marker, under the light of current bat systematics. We analyzed 33 samples collected in an eco-epidemiological survey in the province of Santa Fe (Argentina). We further sequenced 27 museum vouchers to test the accuracy of cytochrome b -based phylogenies in taxonomic identification of bats occurring in the Pampean/Chacoan regions of Argentina. The cytochrome b gene was successfully amplified in all Molossid and Vespertilionid species except for Eptesicus, for which we designed a new reverse primer. The resulting Bayesian phylogeny was congruent with current systematics. Cytochrome b proved useful for species-level delimitation in non-conflicting genera (Eumops, Dasypterus, Molossops) and has infrageneric resolution in more complex lineages (Eptesicus, Myotis, Molossus). We discuss four sources of incongruence that may act separately or in combination: 1) molecular processes, 2) biology, 3) limitations in identification, and 4) errors in the current taxonomy. The present study confirms the general applicability of cytochrome b -based phylogenies in eco-epidemiological studies, but its resolution and reliability depend mainly, but not solely, on the level of genetic differentiation within each bat genus.

First microsatellite markers for the pine catkin sawfly Xyela concava (Hymenoptera, Xyelidae) and their application in phylogeography and population genetics

Microsatellites are widely used as powerful markers in population genetics because of their ability to access recent genetic variation and to resolve subtle population genetic structures. However, their development, especially for non-model organisms with no available genome-wide sequence data has been difficult and time-consuming. Here, a commercial high-throughput sequencing approach (HTS) was used for the very first identification of microsatellite motifs in the genome of Xyela concava and the design of primer pairs flanking those motifs. Sixteen of those primer pairs were selected and implemented successfully to answer questions on the phylogeography and population genetics of X. concava. The markers were characterized in three geographically distinct populations of X. concava and tested for cross-species amplification in two additional Xyela and one Pleroneura species (Xyelidae). All markers showed substantial polymorphism as well as revealing subtle genetic structures among the three genotyped populations. We also analyzed a fragment of the nuclear gene region of sodium/potassium-transporting ATPase subunit alpha (NaK) and a partial mitochondrial gene region coding for cytochrome oxidase subunit I (COI) to demonstrate different genetic resolutions and sex-biased patterns of these markers, and their potential for combined use in future studies on the phylogeography and population genetics of X. concava. Although a limited number of populations was analyzed, we nevertheless obtained new insights on the latter two topics. The microsatellites revealed a generally high gene flow between the populations, but also suggested a deep historical segregation into two genetic lineages. This deep genetic segregation was confirmed by NaK. While the high gene flow was unexpected, because of assumed restricted dispersal ability of X. concava and the discontinuous distribution of the host trees between the populations, the segregation of two lineages is comprehensible and could be explained by different refuge areas of the hosts during glacial times. The COI results showed a discordant strong genetic structure between all populations, which might be explained by the smaller effective population size of the mitochondrial genome. However, given the frequent evidence of a similar nature in recent studies on sawflies, we also consider and discuss mitochondrial introgression on population level as an alternative explanation.

Cross-contamination and strong mitonuclear discordance in Empria sawflies (Hymenoptera, Tenthredinidae) in the light of phylogenomic data

In several sawfly taxa strong mitonuclear discordance has been observed, with nuclear genes supporting species assignments based on morphology, whereas the barcode region of the mitochondrial COI gene suggests different relationships. As previous studies were based on only a few nuclear genes, the causes and the degree of mitonuclear discordance remain ambiguous. Here, we obtained genomic-scale ddRAD data together with Sanger sequences of mitochondrial COI and two to three nuclear protein coding genes to investigate species limits and mitonuclear discordance in two closely related species groups of the sawfly genus Empria. As found previously based on nuclear ITS and mitochondrial COI sequences, species are in most cases supported as monophyletic based on new nuclear data reported here, but not based on mitochondrial COI. This mitonuclear discordance can be explained by occasional mitochondrial introgression with little or no nuclear gene flow, a pattern that might be common in haplodiploid taxa with slowly evolving mitochondrial genomes. Some species in the E. immersa group are not recovered as monophyletic according to either mitochondrial or nuclear data, but this could partly be because of unresolved taxonomy. Preliminary analyses of ddRAD data did not recover monophyly of E. japonica within the E. longicornis group (three Sanger sequenced nuclear genes strongly supported monophyly), but closer examination of the data and additional Sanger sequencing suggested that both specimens were substantially (possibly 10-20% of recovered loci) cross-contaminated. A reason could be specimen identification tag jumps during sequencing library preparation that in previous studies have been shown to affect up to 2.5% of the sequenced reads. We provide an R script to examine patterns of identical loci among the specimens and estimate that the cross-contamination rate is not unusually high for our ddRAD dataset as a whole (based on counting of identical sequences in the immersa and longicornis groups, which are well separated from each other and probably do not hybridise). The high rate of cross-contamination for both E. japonica specimens might be explained by the small number of recovered loci (~1000) compared to most other specimens (>10 000 in some cases) because of poor sequencing results. We caution against drawing unexpected biological conclusions when closely related specimens are pooled before sequencing and tagged only at one end of the molecule or at both ends using a unique combination of limited number of tags (less than the number of specimens).

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.

Holobiont suture zones: Parasite evidence across the European house mouse hybrid zone

Parasite hybrid zones resulting from host secondary contact have never been described in nature although parasite hybridization is well known and secondary contact should affect them similarly to free-living organisms. When host populations are isolated, diverge and recontact, intimate parasites (host specific, direct life cycle) carried during isolation will also meet and so may form parasite hybrid zones. If so, we hypothesize these should be narrower than the host's hybrid zone as shorter parasite generation time allows potentially higher divergence. We investigate multilocus genetics of two parasites across the European house mouse hybrid zone. We find each host taxon harbours its own parasite taxa. These also hybridize: Parasite hybrid zones are significantly narrower than the host's. Here, we show a host hybrid zone is a suture zone for a subset of its parasite community and highlight the potential of such systems as windows on the evolutionary processes of host-parasite interactions and recombinant pathogen emergence.

Conflict between heterozygote advantage and hybrid incompatibility in haplodiploids (and sex chromosomes)

In many diploid species, the sex chromosomes play a special role in mediating reproductive isolation. In haplodiploids, where females are diploid and males haploid, the whole genome behaves similarly to the X/Z chromosomes of diploids. Therefore, haplodiploid systems can serve as a model for the role of sex chromosomes in speciation and hybridization. A previously described population of Finnish Formica wood ants displays genome-wide signs of ploidally and sexually antagonistic selection resulting from hybridization. Here, hybrid females have increased survivorship but hybrid males are inviable. To understand how the unusual hybrid population may be maintained, we developed a mathematical model with hybrid incompatibility, female heterozygote advantage, recombination and assortative mating. The rugged fitness landscape resulting from the co-occurrence of heterozygote advantage and hybrid incompatibility results in a sexual conflict in haplodiploids, which is caused by the ploidy difference. Thus, whereas heterozygote advantage always promotes long-term polymorphism in diploids, we find various outcomes in haplodiploids in which the population stabilizes either in favour of males, females or via maximizing the number of introgressed individuals. We discuss these outcomes with respect to the potential long-term fate of the Finnish wood ant population and provide approximations for the extension of the model to multiple incompatibilities. Moreover, we highlight the general implications of our results for speciation and hybridization in haplodiploids versus diploids and how the described fitness relationships could contribute to the outstanding role of sex chromosomes as hotspots of sexual antagonism and genes involved in speciation.

DNA barcoding a nightmare taxon: assessing barcode index numbers and barcode gaps for sweat bees

There is an ongoing campaign to DNA barcode the world’s >20 000 bee species. Recent revisions of Lasioglossum (Dialictus) (Hymenoptera: Halictidae) for Canada and the eastern United States were completed using integrative taxonomy. DNA barcode data from 110 species of L. (Dialictus) are examined for their value in identification and discovering additional taxonomic diversity. Specimen identification success was estimated using the best close match method. Error rates were 20% relative to current taxonomic understanding. Barcode Index Numbers (BINs) assigned using Refined Single Linkage Analysis (RESL) and barcode gaps using the Automatic Barcode Gap Discovery (ABGD) method were also assessed. RESL was incongruent for 44.5% of species, although some cryptic diversity may exist. Forty-three of 110 species were part of merged BINs with multiple species. The barcode gap is non-existent for the data set as a whole and ABGD showed levels of discordance similar to the RESL. The viridatum species-group is particularly problematic, so that DNA barcodes alone would be misleading for species delimitation and specimen identification. Character-based methods using fixed nucleotide substitutions could improve specimen identification success in some cases. The use of DNA barcoding for species discovery for standard taxonomic practice in the absence of a well-defined barcode gap is discussed.

Widespread hybridization within mound-building wood ants in Southern Finland results in cytonuclear mismatches and potential for sex-specific hybrid breakdown

Hybridization and gene flow between diverging lineages are increasingly recognized as common evolutionary processes, and their consequences can vary from hybrid breakdown to adaptive introgression. We have previously found a population of wood ant hybrids between Formica aquilonia and F. polyctena that shows antagonistic effects of hybridization: females with introgressed alleles show hybrid vigour, whereas males with the same alleles show hybrid breakdown. Here, we investigate whether hybridization is a general phenomenon in this species pair and analyse 647 worker samples from 16 localities in Finland using microsatellite markers and a 1200-bp mitochondrial sequence. Our results show that 27 sampled nests contained parental-like gene pools (six putative F. polyctena and 21 putative F. aquilonia) and all remaining nests (69), from nine localities, contained hybrids of varying degrees. Patterns of genetic variation suggest these hybrids arise from several hybridization events or, instead, have backcrossed to the parental gene pools to varying extents. In contrast to expectations, the mitochondrial haplotypes of the parental species were not randomly distributed among the hybrids. Instead, nests that were closer to parental-like F. aquilonia for nuclear markers preferentially had F. polyctena's mitochondria and vice versa. This systematic pattern suggests there may be underlying selection favouring cytonuclear mismatch and hybridization. We also found a new hybrid locality with strong genetic differences between the sexes similar to those predicted under antagonistic selection on male and female hybrids. Further studies are needed to determine the selective forces that act on male and female genomes in these newly discovered hybrids.

The on-again, off-again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co-evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.

Sex Determination, Sex Chromosomes, and Karyotype Evolution in Insects

Insects harbor a tremendous diversity of sex determining mechanisms both within and between groups. For example, in some orders such as Hymenoptera, all members are haplodiploid, whereas Diptera contain species with homomorphic as well as male and female heterogametic sex chromosome systems or paternal genome elimination. We have established a large database on karyotypes and sex chromosomes in insects, containing information on over 13000 species covering 29 orders of insects. This database constitutes a unique starting point to report phylogenetic patterns on the distribution of sex determination mechanisms, sex chromosomes, and karyotypes among insects and allows us to test general theories on the evolutionary dynamics of karyotypes, sex chromosomes, and sex determination systems in a comparative framework. Phylogenetic analysis reveals that male heterogamety is the ancestral mode of sex determination in insects, and transitions to female heterogamety are extremely rare. Many insect orders harbor species with complex sex chromosomes, and gains and losses of the sex-limited chromosome are frequent in some groups. Haplodiploidy originated several times within insects, and parthenogenesis is rare but evolves frequently. Providing a single source to electronically access data previously distributed among more than 500 articles and books will not only accelerate analyses of the assembled data, but also provide a unique resource to guide research on which taxa are likely to be informative to address specific questions, for example, for genome sequencing projects or large-scale comparative studies.

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.

Species-Level Para- and Polyphyly in DNA Barcode Gene Trees: Strong Operational Bias in European Lepidoptera

The proliferation of DNA data is revolutionizing all fields of systematic research. DNA barcode sequences, now available for millions of specimens and several hundred thousand species, are increasingly used in algorithmic species delimitations. This is complicated by occasional incongruences between species and gene genealogies, as indicated by situations where conspecific individuals do not form a monophyletic cluster in a gene tree. In two previous reviews, non-monophyly has been reported as being common in mitochondrial DNA gene trees. We developed a novel web service “Monophylizer” to detect non-monophyly in phylogenetic trees and used it to ascertain the incidence of species non-monophyly in COI (a.k.a. cox1) barcode sequence data from 4977 species and 41,583 specimens of European Lepidoptera, the largest data set of DNA barcodes analyzed from this regard. Particular attention was paid to accurate species identification to ensure data integrity. We investigated the effects of tree-building method, sampling effort, and other methodological issues, all of which can influence estimates of non-monophyly. We found a 12% incidence of non-monophyly, a value significantly lower than that observed in previous studies. Neighbor joining (NJ) and maximum likelihood (ML) methods yielded almost equal numbers of non-monophyletic species, but 24.1% of these cases of non-monophyly were only found by one of these methods. Non-monophyletic species tend to show either low genetic distances to their nearest neighbors or exceptionally high levels of intraspecific variability. Cases of polyphyly in COI trees arising as a result of deep intraspecific divergence are negligible, as the detected cases reflected misidentifications or methodological errors. Taking into consideration variation in sampling effort, we estimate that the true incidence of non-monophyly is ∼23%, but with operational factors still being included. Within the operational factors, we separately assessed the frequency of taxonomic limitations (presence of overlooked cryptic and oversplit species) and identification uncertainties. We observed that operational factors are potentially present in more than half (58.6%) of the detected cases of non-monophyly. Furthermore, we observed that in about 20% of non-monophyletic species and entangled species, the lineages involved are either allopatric or parapatric—conditions where species delimitation is inherently subjective and particularly dependent on the species concept that has been adopted. These observations suggest that species-level non-monophyly in COI gene trees is less common than previously supposed, with many cases reflecting misidentifications, the subjectivity of species delimitation or other operational factors.

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.