First record of the bacterial endosymbiont Wolbachia for phytophagous hoverflies from genus Merodon (Diptera: Syrphidae)

Wolbachia is a widespread bacterial endosymbiont among arthropod species. It influences the reproduction of the host species and also mitochondrial DNA diversity. Until now there were only a few studies that detected Wolbachia infections in hoverflies (Diptera: Syrphidae), and this is the first broader study with the aim of examining the incidence of Wolbachia in the hoverfly genus Merodon. The obtained results indicate an infection rate of 96% and the presence of both Wolbachia supergroup A and B, which are characteristic for most of the infected arthropod species. Additionally, the presence of multiple Wolbachia strains in the Merodon aureus group species was detected and the mitochondrial DNA COI‐based relationships of the group are discussed in the light of infection. Finally, we discuss plant‐mediated horizontal transmission of Wolbachia strains among the studied hoverfly species.     

Protein–protein interactions of the cytochrome c oxidase DNA barcoding region

Enzymatic amplification of homologous regions of DNA using ‘universal’ polymerase chain reaction primers has provided insight into insect systematics, phylogeography, molecular evolution and species identification. One of the more commonly amplified and sequenced regions is a short region of the cytochrome c oxidase subunit I gene (COI), commonly called the barcoding region. COI is one of three mitochondrial‐encoded subunits of complex IV (Cox) of the electron transport chain. In addition to the mitochondrial subunits there are nine nuclear‐encoded subunits of the complex in Drosophila. Whereas a number of phylogenetic biases associated with this region have been examined and the quaternary structure of Cox has been modelled, the influence of protein–protein interactions on the observed patterns of evolution in this barcoding region of insects has never been examined critically. Using a well‐resolved independently derived phylogeny of 38 Diptera species, we examined the homogeneity of the substitution processes within the barcoding region. We show that, within Diptera, amino acid residues interacting with nuclear‐encoded subunits of Cox are evolving at elevated rates across the phylogeny. Furthermore, we show that codon position two is biased by protein–protein interactions. In contrast, third codon positions provide a less biased estimate of genetic variation in the region. This study highlights the need to examine the potential for systematic bias in DNA barcoding regions as part of the critical assessment of evidence in systematics and in biodiversity assessments.     

Genetic variation amongst biotypes of Dactylopius tomentosus

The tomentose cochineal scale insect, Dactylopius tomentosus (Lamarck) (Hemiptera: Dactylopiidae), is an important biological control agent against invasive species of Cylindropuntia (Caryophyllales: Cactaceae). Recent studies have demonstrated that this scale is composed of host‐affiliated biotypes with differential host specificity and fitness on particular host species. We investigated genetic variation and phylogenetic relationships among D. tomentosus biotypes and provenances to examine the possibility that genetic diversity may be related to their host‐use pattern, and whether their phylogenetic relationships would give insights into taxonomic relatedness of their host plants. Nucleotide sequence comparison was accomplished using sequences of the mitochondrial cytochrome c oxidase I (COI) gene. Sequences of individuals from the same host plant within a region were identical and characterized by a unique haplotype. Individuals belonging to the same biotype but from different regions had similar haplotypes. However, haplotypes were not shared between different biotypes. Phylogenetic analysis grouped the monophyletic D. tomentosus into 3 well‐resolved clades of biotypes. The phylogenetic relationships and clustering of biotypes corresponded with known taxonomic relatedness of their hosts. Two biotypes, Fulgida and Mamillata, tested positive for Wolbachia (α‐Proteobacteria), a common endosymbiont of insects. The Wolbachia sequences were serendipitously detected by using insect‐specific COI DNA barcoding primers and are most similar to Wolbachia Supergroup F strains. This study is the first molecular characterization of cochineal biotypes that, together with Wolbachia sequences, contribute to the better identification of the biotypes of cochineal insects and to the biological control of cacti using host‐specific biotypes of the scale.     

DNA barcoding of crickets,katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria,Germany and Switzerland

We present a DNA barcoding study on the insect order Orthoptera that was generated in collaboration between four barcoding projects in three countries, viz. Barcoding Fauna Bavarica (Germany), German Barcode of Life, Austrian Barcode of Life and Swiss Barcode of Life. Our data set includes 748 COI sequences from 127 of the 162 taxa (78.4%) recorded in the three countries involved. Ninety‐three of these 122 species (76.2%, including all Ensifera) can be reliably identified using DNA barcodes. The remaining 26 caeliferan species (families Acrididae and Tetrigidae) form ten clusters that share barcodes among up to five species, in three cases even across different genera, and in six cases even sharing individual barcodes. We discuss incomplete lineage sorting and hybridization as most likely causes of this phenomenon, as the species concerned are phylogenetically young and hybridization has been previously observed. We also highlight the problem of nuclear mitochondrial pseudogenes (numts), a known problem in the barcoding of orthopteran species, and the possibility of Wolbachia infections. Finally, we discuss the possible taxonomic implications of our barcoding results and point out future research directions.     

The gnd gene of Buchnera as a new,effective DNA barcode for aphid identification

DNA barcoding uses a standard DNA sequence to facilitate species identification. Although the COI gene has been adopted as the standard, COI alone is imperfect due to several shortcomings. The primary endosymbiont of aphids, Buchnera, has higher evolutionary rates and interspecies divergence than its co‐diverging aphid hosts, making it a potential tool for resolving the ambiguities in aphid taxonomy. We compared the effectiveness of employing two different DNA regions, gnd and COI, for the discrimination of over 100 species of aphids. The mean interspecific divergence of the gnd region was significantly higher than the mean intraspecific variation; there were nearly nonoverlapping distributions between the intra‐ and interspecific samples. In contrast, COI showed a lower interspecific divergence, which led to difficulties in identifying closely related species. Our results show that gnd can identify species in the Aphididae, which suggests that the gnd region of Buchnera is a potentially effective barcode for aphid species identification. We also recommend the 2‐locus combination of gnd + COI as the aphid barcode. This will provide a universal framework for the routine use of DNA sequence data to identify specimens and contribute toward the discovery of overlooked species of aphids.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Morphology versus DNA barcoding: two sides of the same coin. A case study of Ceutorhynchus erysimi and C. contractus identification

Genotyping of 2 well‐known weevil species from the genus Ceutorhynchus (Coleoptera: Curculionidae) distributed in west Palearctic, C. erysimi and C. contractus, revealed phenotype versus genotype inconsistencies in a set of 56 specimens (25 C. erysimi and 31 C. contractus) collected from 25 locations in Serbia and Montenegro. An analysis of the mitochondrial cytochrome oxidase subunit I gene (COI), widely used as a barcoding region, and a nuclear gene, elongation factor‐1α (EF‐1α), revealed stable genetic divergence among these species. The average uncorrected pairwise distances for the COI and EF‐1α genes were 3.8%, and 1.3%, respectively, indicating 2 genetically well‐segregated species. However, the genetic data were not congruent with the phenotypic characteristics of the studied specimens. In the first place, C. erysimi genotypes were attached to specimens with phenotypic characteristics of C. contractus. Species‐specific PCR‐RFLP assays for the barcoding gene COI were applied for the molecular identification of 101 additional specimens of both morphospecies (33 C. erysimi and 68 C. contractus) and were found to confirm this incongruity. The discrepancy between the genetic and morphological data raises the question of the accuracy of using a barcoding approach, as it may result in misleading conclusions about the taxonomic position of the studied organism. Additionally, the typological species concept shows considerable weakness when genetic data are not supported with phenotypic characteristics as in case of asymmetric introgression, which may cause certain problems, especially in applied studies such as biological control programs in which the biological properties of the studied organisms are the main focus.     

Cold Code: the global initiative to DNA barcode amphibians and nonavian reptiles

DNA barcoding facilitates the identification of species and the estimation of biodiversity by using nucleotide sequences, usually from the mitochondrial genome. Most studies accomplish this task by using the gene encoding cytochrome oxidase subunit I (COI; Entrez COX1). Within this barcoding framework, many taxonomic initiatives exist, such as those specializing in fishes, birds, mammals, and fungi. Other efforts center on regions, such as the Arctic, or on other topics, such as health. DNA barcoding initiatives exist for all groups of vertebrates except for amphibians and nonavian reptiles. We announce the formation of Cold Code, the international initiative to DNA barcode all species of these ‘cold‐blooded’ vertebrates. The project has a Steering Committee, Coordinators, and a home page. To facilitate Cold Code, the Kunming Institute of Zoology, Chinese Academy of Sciences will sequence COI for the first 10 specimens of a species at no cost to the steward of the tissues.     

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.     

First report on Wolbachia endosymbiosis in freshwater Aphelocheirus aestivalis (Heteroptera: Aphelocheiridae) and its potential impact on genetic diversity of host

Alphaproteobacteria Wolbachia have been described as endosymbionts of approximately half of all aquatic insect species. These bacteria might affect not only reproduction but also the genetic diversity of its hosts. In the present study we identified Wolbachia endosymbiosis in freshwater true bug Aphelocheirus aestivalis F., 1794 (Heteroptera: Aphelocheiridae). Despite the fact that A. aestivalis is widely distributed in Europe, it occurs rather locally, often in isolated populations. Taking into account that Wolbachia, close relationships and past demographic phenomena could affect the genetic diversity of its host, we analyzed mitochondrial (COI and 16S) and nuclear (internal transcribed spacer 2) markers determined for A. aestivalis individuals collected from five populations. Moreover, we compared obtained COI sequences with those deposited in GenBank. Analyses revealed low genetic differentiation among samples tested, as well as low variation among determined COI sequences and those downloaded from the database. Although Wolbachia infection could correlate with decreasing mitochondrial diversity of its host, we suggest that low genetic variation observed in tested A. aestivalis samples (at both mitochondrial and nuclear levels) is a result of populations’ close relationships, past demographic phenomena or is characteristic for this species. Detailed analysis of the wsp gene fragment revealed two distinct strains of Wolbachia infecting A. aestivalis. Both of them belong to supergroup A, also found in other arthropods.     

Molecular phylogenetic analysis of New Zealand mosquito species

ABSTRACT

Taxonomic information for New Zealand mosquito species is predominantly morphological with very few molecular data available to date. In this study, the mitochondrial COI gene and nuclear ITS1 were amplified and sequenced from DNA templates representing 17 species; 15 previously known New Zealand species, a newly discovered undescribed Aedes species from the Chatham Islands and a recently eradicated invader, Ae. camptorhynchus. This paper provides DNA barcoding sequences for the entire known New Zealand mosquito fauna, the first for the majority of these species. Phylogenetic analysis of COI and ITS1 indicated that the endemic species are all genetically distinct from the exotic species examined including vector species of pathogens of human diseases. The genus Opifex is distant from the genus Aedes based on analysis of ITS1 sequences, and Ae. chathamicus is more closely related to species within Aedes than to the genus Opifex. Results show Culex asteliae to be a valid species but Cx. rotoruae is not necessarily so. The Aedes species appears to be closely related to Ae. subalbirostris. The introduced Cx. quinquefasciatus was clearly shown to be related to Cx. pipiens; however, the endemic Cx. pervigilans was not. No evidence of population variation based on geographic location was detected.     

High mitochondrial diversity in geographically widespread butterflies of Madagascar: A test of the DNA barcoding approach

The standardized use of mitochondrial cytochrome c oxidase subunit I (COI) gene sequences as DNA barcodes has been widely promoted as a high-throughput method for species identification and discovery. Species delimitation has been based on the following criteria: (1) monophyletic association and less frequently (2) a minimum 10× greater divergence between than within species. Divergence estimates, however, can be inflated if sister species pairs are not included and the geographic extent of variation within any given taxon is not sampled comprehensively. This paper addresses both potential biases in DNA divergence estimation by sampling range-wide variation in several morphologically distinct, endemic butterfly species in the genus Heteropsis, some of which are sister taxa. We also explored the extent to which mitochondrial DNA from the barcode region can be used to assess the effects of historical rainforest fragmentation by comparing genetic variation across Heteropsis populations with an unrelated forest-associated taxon Saribia tepahi. Unexpectedly, generalized primers led to the inadvertent amplification of the endosymbiont Wolbachia, undermining the use of universal primers and necessitating the design of genus-specific COI primers alongside a Wolbachia-specific PCR assay. Regardless of the high intra-specific genetic variation observed, most species satisfy DNA barcoding criteria and can be differentiated in the nuclear phylogeny. Nevertheless, two morphologically distinguishable candidate species fail to satisfy the barcoding 10× genetic distance criterion, underlining the difficulties of applying a standard distance threshold to species delimitation. Phylogeographic analysis of COI data suggests that forest fragmentation may have played an important role in the recent evolutionary diversification of these butterflies. Further work on other Malagasy taxa using both mitochondrial and nuclear data will provide better insight into the role of historical habitat fragmentation in species diversification and may potentially contribute to the identification of priority areas for conservation.     

Phylogeography of the Temminck’s Stint (Calidris temminckii): historical vicariance but little present genetic structure in a regionally endangered Palearctic wader

Aim We study the population differentiation and phylogeography of the Temminck’s Stint (Calidris temminckii). Specifically, we seek signs of past and present population size changes and dispersal events and evaluate management and conservation unit status of the populations. We also study the possibility of introgression as the origin of two mitochondrial DNA (mtDNA) lineages found and estimate the divergence time of the lineages. Location Northern Eurasia. Methods We analysed 583 bp of mtDNA control region domains I and II and 11 microsatellite loci from 13 localities throughout the breeding range. In addition, we used mitochondrial cytochrome c oxidase subunit I (COI), a barcoding gene, to search for signs of introgression. Results More population differentiation was found from microsatellites than from mtDNA, although differentiation was weak in both markers. Signs of past population growth were observed, in addition to more recent decline in some areas. Both control region and COI sequences revealed two maternal lineages coexisting in Fennoscandia and in north‐west Siberia. No signs of introgression were detected. Lineage divergence time was estimated to have occurred during the glacial periods of Pleistocene. Main conclusions Slight differences in mtDNA and microsatellite differentiation and diversity may reflect different features – such as the mutation rate and effective population size – of the markers used, or female‐biased dispersal pattern and high male site‐fidelity of the species. The coexistence of the two mitochondrial lineages is most likely a consequence of post‐glacial mixing of two refugial Pleistocene populations. Based on genetic information alone, global conservation concerns are not imminent. However, fast decline of a marginal Bothnian Bay population and the smallness and remoteness of a Central Yakutian population warrant conservation actions.     

Widespread mito‐nuclear discordance with evidence for introgressive hybridization and selective sweeps in Lycaeides

We investigated the extent and potential cause(s) of mitochondrial introgression within the polytypic North American Lycaeides species complex (Lepidoptera). By comparing population genetic structure based on mitochondrial DNA (COI and COII) and nuclear DNA (251 polymorphic amplified fragment length polymorphism markers), we detected substantial mito‐nuclear discordance, primarily involving a single mitochondrial haplotype (h01), which is likely due to mitochondrial introgression between differentiated Lycaeides populations and/or species. We detected reduced mitochondrial genetic diversity relative to nuclear genetic diversity in populations where mitochondrial haplotype h01 occurs, suggesting that the spread of this haplotype was facilitated by selection. We found no evidence that haplotype h01 is associated with increased fitness (in terms of survival to eclosion, fresh adult weight, and adult longevity) in a polymorphic Lycaeides melissa population. However, we did find a positive association between mitochondrial haplotype h01 and infection by the endoparasitic bacterium Wolbachia in one out of three lineages tested. Linkage disequilibrium between mitochondrial haplotype h01 and Wolbachia infection status may have resulted in indirect selection favouring the spread of haplotype h01 in at least one lineage of North American Lycaeides. These results illustrate the potential for introgressive hybridization to produce substantial mito‐nuclear discordance and demonstrate that an individual's mitochondrial and nuclear genome may have strikingly different evolutionary histories resulting from non‐neutral processes and intrinsic differences in the inheritance and biology of these genomes.     

Assessing the potential of candidate DNA barcodes for identifying non‐flowering seed plants

In plants, matK and rbcL have been selected as core barcodes by the Consortium for the Barcode of Life (CBOL) Plant Working Group (PWG), and ITS/ITS2 and psbA‐trnH were suggested as supplementary loci. Yet, research on DNA barcoding of non‐flowering seed plants has been less extensive, and the evaluation of DNA barcodes in this division has been limited thus far. Here, we evaluated seven markers (psbA‐trnH, matK, rbcL, rpoB, rpoC1, ITS and ITS2) from non‐flowering seed plants. The usefulness of each region was assessed using four criteria: the success rate of PCR amplification, the differential intra‐ and inter‐specific divergences, the DNA barcoding gap and the ability to discriminate species. Among the seven loci tested, ITS2 produced the best results in the barcoding of non‐flowering seed plants. In addition, we compared the abilities of the five most‐recommended markers (psbA‐trnH, matK, rbcL, ITS and ITS2) to identify additional species using a large database of gymnosperms from GenBank. ITS2 remained effective for species identification in a wide range of non‐flowering seed plants: for the 1531 samples from 608 species of 80 diverse genera, ITS2 correctly authenticated 66% of them at the species level. In conclusion, the ITS2 region can serve as a useful barcode to discriminate non‐flowering seed plants, and this study will contribute valuable information for the barcoding of plant species.     

Hidden diversity of Ctenophora revealed by new mitochondrial COI primers and sequences

The mitochondrial gene cytochrome-c-oxidase subunit 1 (COI) is useful in many taxa for phylogenetics, population genetics, metabarcoding, and rapid species identifications. However, the phylum Ctenophora (comb jellies) has historically been difficult to study due to divergent mitochondrial sequences and the corresponding inability to amplify COI with degenerate and standard COI “barcoding” primers. As a result, there are very few COI sequences available for ctenophores, despite over 200 described species in the phylum. Here, we designed new primers and amplified the COI fragment from members of all major groups of ctenophores, including many undescribed species. Phylogenetic analyses of the resulting COI sequences revealed high diversity within many groups that was not evident from more conserved 18S rDNA sequences, in particular among the Lobata (Ctenophora; Tentaculata; Lobata). The COI phylogenetic results also revealed unexpected community structure within the genus Bolinopsis, suggested new species within the genus Bathocyroe, and supported the ecological and morphological differences of some species such as Lampocteis cruentiventer and similar undescribed lobates (Lampocteis sp. “V” stratified by depth, and “A” differentiated by colour). The newly designed primers reported herein provide important tools to enable researchers to illuminate the diversity of ctenophores worldwide via quick molecular identifications, improve the ability to analyse environmental DNA by improving reference libraries and amplifications, and enable a new breadth of population genetic studies.     

DNA barcoding of microgastrine parasitoid wasps (Hymenoptera: Braconidae) using high‐throughput methods more than doubles the number of species known for Australia

The Microgastrinae are a hugely diverse subfamily of endoparasitoid wasps of lepidopteran caterpillars. They are important in agriculture as biological control agents and play a significant ecological role in the regulation of caterpillar populations. Whilst the group has been the focus of intensive rearing and DNA barcoding studies in the Northern Hemisphere, the Australian fauna has received little attention. In total, 99 species have been described from or have been introduced into Australia, but the real species diversity for the region is clearly much larger than this. In this study, museum ethanol samples and recent field collections were mined for hundreds of specimens of microgastrine wasps, which were then barcoded for the COI region, ITS2 ribosomal spacer and the wingless nuclear genes, using a pooled sequencing approach on an Illumina Miseq system. Full COI sequences were obtained for 525 individuals which, when combined with 162 publicly available sequences, represented 417 haplotypes, and a total of 236 species were delimited using a consensus approach. By more than doubling the number of known microgastrine wasp species in Australia, our study highlights the value of DNA barcoding in the context of employing high‐throughput sequencing methods of bulk ethanol museum collections for biodiversity assessment.     

Development of internal COI primers to improve and extend barcoding of fruit flies (Diptera: Tephritidae: Dacini)

Accurate species-level identifications underpin many aspects of basic and applied biology;however,identifications can be hampered by a lack of discriminating morphological characters,taxonomic expertise or time.Molecular approaches,such as DNA"barcoding"of the cytochrome c oxidase(COI)gene,are argued to overcome these issues.However,nuclear encoding of mitochondrial genes(numts)and poor amplification success of suboptimally preserved specimens can lead to erroneous identifications.One insect group for which these molecular and morphological problems are significant are the dacine fruit flies(Diptera:Tephritidae:Dacini).We addressed these issues associated with COI barcoding in the dacines by first assessing several"universal"COI primers against public mitochondrial genome and numt sequences for dacine taxa.We then modified a set of four primers that more closely matched true dacine COI sequence and amplified two overlapping portions of the COI barcode region.Our new primers were tested alongside universal primers on a selection of dacine species,including both fresh preserved and decades-old dry specimens.Additionally,Bactrocera tiyoni mitochondrial and nuclear genomes were compared to identify putative numts.Four numt clades were identified,three of which were amplified using existing universal primers.In contrast,our new primers preferentially amplified the"true"mitochondrial COI barcode in all dacine species tested.The new primers also successfully amplified partial barcodes from dry specimens for which full length barcodes were unobtainable.Thus we recommend these new primers be incorporated into the suites of primers used by diagnosticians and quarantine labs for the accurate identification of dacine species.     

Wolbachia have made it twice: Hybrid introgression between two sister species of Eurema butterflies

Wolbachia, cytoplasmically inherited endosymbionts of arthropods, are known to hijack their host reproduction in various ways to increase their own vertical transmission. This may lead to the selective sweep of associated mitochondria, which can have a large impact on the evolution of mitochondrial lineages. In Japan, two different Wolbacahia strains (wCI and wFem) are found in two sister species of pierid butterflies, Eurema mandarina and Eurema hecabe. In both species, females infected with wCI (C females) produce offspring with a nearly 1:1 sex ratio, while females infected with both wCI and wFem (CF females) produce all‐female offspring. Previous studies have suggested the historical occurrence of hybrid introgression in C individuals between the two species. Furthermore, hybrid introgression in CF individuals is suggested by the distinct mitochondrial lineages between C females and CF females of E. mandarina. In this study, we performed phylogenetic analyses based on nuclear DNA and mitochondrial DNA markers of E. hecabe with previously published data on E. mandarina. We found that the nuclear DNA of this species significantly diverged from that of E. mandarina. By contrast, mitochondrial DNA haplotypes comprised two clades, mostly reflecting Wolbachia infection status rather than the individual species. Collectively, our results support the previously suggested occurrence of two independent historical events wherein the cytoplasms of CF females and C females moved between E. hecabe and E. mandarina through hybrid introgression.