Dominance effects strengthen premating hybridization barriers between sympatric species of grasshoppers (Acrididae,Orthoptera)

Sexual selection can lead to the rapid evolution of premating hybridization barriers and allows accelerated diversification and speciation within an evolutionary lineage. Especially during early stages of divergence, hybridization may impede further divergence, which strongly depends on the reproductive success of hybrids. Behavioural sterility of hybrids can limit or even prevent homogenizing gene flow. In this study, we investigated the attractiveness of male courtship songs for females of the grasshopper species Chorthippus biguttulus and C. brunneus and their interspecific F1 and F2 hybrids. Song preferences of females of both species are highly species specific and differ in three parameters: shape of the preference function, preference for syllable pattern and phrase duration. F1 hybrid females of both reciprocal crosses as well as F2 hybrid females resembled closely pure C. biguttulus females in respect of shape of the preference function and preference for syllable pattern, while preference for phrase duration showed an intermediate expression. This resulted in song preferences of hybrid females that closely resembled those of one parental species, that is C. biguttulus females. Such strong dominance effects were rarely reported so far. They represent an effective barrier limiting gene flow between the two species, since hybrid females will backcross to only one parental species and discriminate against hybrid males, which are behaviourally sterile. Such taxon‐specific modes of inheritance may have facilitated the rapid divergence of acoustically communicating grasshoppers of the species group of Chorthippus biguttulus. Our findings have novel implications on the expression of neuronal filters and the evolution of complex courtship signals.     

Large‐scale introduction of the Indo‐Pacific damselfish Abudefduf vaigiensis into Hawai'i promotes genetic swamping of the endemic congener A. abdominalis

Hybridization in the ocean was once considered rare, a process prohibited by the rapid evolution of intrinsic reproductive barriers in a high‐dispersal medium. However, recent genetic surveys have prompted a reappraisal of marine hybridization as an important demographic and evolutionary process. The Hawaiian Archipelago offers an unusual case history in this arena, due to the recent arrival of the widely distributed Indo‐Pacific sergeant (Abudefduf vaigiensis), which is hybridizing with the endemic congener, A. abdominalis. Surveys of mtDNA and three nuclear loci across Hawai'i (N = 396, Abudefduf abdominalis and N = 314, A. vaigiensis) reveal that hybridization is significantly higher in the human‐perturbed southeast archipelago (19.8%), tapering off to 5.9% in the pristine northwest archipelago. While densities of the two species varied throughout Hawai'i, hybridization was highest in regions with similar species densities, contradicting the generalization that the rarity of one species promotes interspecific mating. Our finding of later generation hybrids throughout the archipelago invokes the possibility of genetic swamping of the endemic species. Exaptation, an adaptation with unintended consequences, may explain these findings: the endemic species has transient yellow coloration during reproduction, whereas the introduced species has yellow coloration continuously as adults, in effect a permanent signal of reproductive receptivity. Haplotype diversity is higher in Hawaiian A. vaigiensis than in our samples from the native range, indicating large‐scale colonization almost certainly facilitated by the historically recent surge of marine debris. In this chain of events, marine debris promotes colonization, exaptation promotes hybridization, and introgression invokes the possible collapse of an endemic species.     

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross‐species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next‐generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross‐genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter‐ and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.     

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species‐specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri‐specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure ), bidirectional estimates of gene flow (migrate ‐n and IMa2) and phylogeny‐based methods (*beast , buck y) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species‐specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri‐like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.     

Does genetic distance between parental species influence outcomes of hybridization among coral reef butterflyfishes?

Christmas Island is located at the overlap of the Indian and Pacific Ocean marine provinces and is a hot spot for marine hybridization. Here, we evaluate the ecological framework and genetic consequences of hybridization between butterflyfishes Chaetodon guttatissimus and Chaetodon punctatofasciatus. Further, we compare our current findings to those from a previous study of hybridization between Chaetodon trifasciatus and Chaetodon lunulatus. For both species groups, habitat and dietary overlap between parental species facilitate frequent heterospecific encounters. Low abundance of potential mates promotes heterospecific pair formation and the breakdown of assortative mating. Despite similarities in ecological frameworks, the population genetic signatures of hybridization differ between the species groups. Mitochondrial and nuclear data from C. guttatissimus × C. punctatofasciatus (1% divergence at cyt b) show bidirectional maternal contributions and relatively high levels of introgression, both inside and outside the Christmas Island hybrid zone. In contrast, C. trifasciatus × C. lunulatus (5% cyt b divergence) exhibit unidirectional mitochondrial inheritance and almost no introgression. Back‐crossing of hybrid C. guttatissimus × C. punctatofasciatus and parental genotypes may eventually confound species‐specific signals within the hybrid zone. In contrast, hybrids of C. trifasciatus and C. lunulatus may coexist with and remain genetically distinct from the parents. Our results, and comparisons with hybridization studies in other reef fish families, indicate that genetic distance between hybridizing species may be a factor influencing outcomes of hybridization in reef fish, which is consistent with predictions from terrestrially derived hybridization theory.     

Females discriminate against heterospecific sperm in a natural hybrid zone

When hybridization is maladaptive, species‐specific mate preferences are selectively favored, but low mate availability may constrain species‐assortative pairing. Females paired to heterospecifics may then benefit by copulating with multiple males and subsequently favoring sperm of conspecifics. Whether such mechanisms for biasing paternity toward conspecifics act as important reproductive barriers in socially monogamous vertebrate species remains to be determined. We use a combination of long‐term breeding records from a natural hybrid zone between collared and pied flycatchers (Ficedula albicollis and F. hypoleuca), and an in vitro experiment comparing conspecific and heterospecific sperm performance in female reproductive tract fluid, to evaluate the potential significance of female cryptic choice. We show that the females most at risk of hybridizing (pied flycatchers) frequently copulate with multiple males and are able to inhibit heterospecific sperm performance. The negative effect on heterospecific sperm performance was strongest in pied flycatcher females that were most likely to have been previously exposed to collared flycatcher sperm. We thus demonstrate that a reproductive barrier acts after copulation but before fertilization in a socially monogamous vertebrate. While the evolutionary history of this barrier is unknown, our results imply that there is opportunity for it to be accentuated via a reinforcement‐like process.     

Evolutionary neutrality of mtDNA introgression: evidence from complete mitogenome analysis in roe deer

Introgressive hybridization offers a unique platform for studying the molecular basis of natural selection acting on mitogenomes. Most of the mtDNA protein‐coding genes are extremely conserved; however, some of the observed variations have potentially adaptive significance. Here, we evaluated whether the evolution of mtDNA in closely related roe deer species affected by widespread mtDNA introgression is neutral or adaptive. We characterized and compared 16 complete mitogenomes of European (Capreolus capreolus) and Siberian (C. pygargus) roe deer, including four of Siberian origin introgressed into European species. The average sequence divergence of species‐specific lineages was estimated at 2.8% and varied across gene classes. Only 21 of 315 fixed differences identified in protein‐coding genes represented nonsynonymous changes. Only three of them were determined to have arisen in the C. pygargus lineage since the time to the most recent common ancestor (TMRCA) of both Capreolus species, reflecting a decelerated evolutionary ratio. The almost four‐fold higher d_N/d_S ratio described for the European roe deer lineage is constrained by overall purifying selection, especially pronounced in the ND4 and ND5 genes. We suggest that the highly divergent C. capreolus lineage could have maintained a capability for genomic incorporation of the well‐preserved and almost ancestral type of mtDNA present in C. pygargus. Our analyses did not indicate any signs of positive selection for Siberian roe deer mtDNA, suggesting that the present widespread introgression is evolutionarily neutral.     

Hierarchical behaviour,habitat use and species size differences shape evolutionary outcomes of hybridization in a coral reef fish

Hybridization is an important evolutionary process, with ecological and behavioural factors influencing gene exchange between hybrids and parent species. Patterns of hybridization in anemonefishes may result from living in highly specialized habitats and breeding status regulated by size‐based hierarchal social groups. Here, morphological, ecological and genetic analyses in Kimbe Bay, Papua New Guinea, examine the hybrid status of Amphiprion leucokranos, a nominal species and presumed hybrid between Amphiprion sandaracinos and Amphiprion chrysopterus. We test the hypothesis that habitat use and relative size differences of the parent species and hybrids determine the patterns of gene exchange. There is strong evidence that A. leucokranos is a hybrid of smaller A. sandaracinos and larger A. chrysopterus, where A. chrysopterus is exclusively the mother to each hybrid, based on mtDNA cytochrome b and multiple nDNA microsatellite loci. Overlap in habitat, depth and host anemone use was found, with hybrids intermediate to parents and cohabitation in over 25% of anemones sampled. Hybrids, intermediate in body size, colour and pattern, were classified 55% of the time as morphologically first‐generation hybrids relative to parents, whereas 45% of hybrids were more A. sandaracinos‐like, suggesting backcrossing. Unidirectional introgression of A. chrysopterus mtDNA into A. sandaracinos via hybrid backcrosses was found, with larger female hybrids and small male A. sandaracinos mating. Potential nDNA introgression was also evident through distinct intermediate hybrid genotypes penetrating both parent species. Findings support the hypothesis that anemonefish hierarchical behaviour, habitat use and species‐specific size differences determine how hybrids form and the evolutionary consequences of hybridization.     

Experimental evolution of competing bean beetle species reveals long‐term reversals of short‐term evolution,but no consistent character displacement

Interspecific competition for shared resources should select for evolutionary divergence in resource use between competing species, termed character displacement. Many purported examples of character displacement exist, but few completely rule out alternative explanations. We reared genetically diverse populations of two species of bean beetles, Callosobruchus maculatus and Callosobruchus chinensis, in allopatry and sympatry on a mixture of adzuki beans and lentils, and assayed oviposition preference and other phenotypic traits after four, eight, and twelve generations of (co)evolution. C. maculatus specializes on adzuki beans; the generalist C. chinensis uses both beans. C. chinensis growing in allopatry emerged equally from both bean species. In sympatry, the two species competing strongly and coexisted via strong realized resource partitioning, with C. chinensis emerging almost exclusively from lentils and C. maculatus emerging almost exclusively from adzuki beans. However, oviposition preferences, larval survival traits, and larval development rates in both beetle species did not vary consistently between allopatric versus sympatric treatments. Rather, traits evolved in treatment‐independent fashion, with several traits exhibiting reversals in their evolutionary trajectories. For example, C. chinensis initially evolved a slower egg‐to‐adult development rate on adzuki beans in both allopatry and sympatry, then subsequently evolved back toward the faster ancestral development rate. Lack of character displacement is consistent with a previous similar experiment in bean beetles and may reflect lack of evolutionary trade‐offs in resource use. However, evolutionary reversals were unexpected and remain unexplained. Together with other empirical and theoretical work, our results illustrate the stringency of the conditions for character displacement.     

The role of functional constraints in nonrandom mating patterns for a dance fly with female ornaments

Most hypotheses to explain nonrandom mating patterns invoke mate choice, particularly in species that display elaborate ornaments. However, conflicting selection pressures on traits can result in functional constraints that can also cause nonrandom mating patterns. We tested for functional load‐lifting constraints during aerial copulation in Rhamphomyia longicauda, a species of dance fly that displays multiple extravagant female‐specific ornaments that are unusual among sexual traits because they are under stabilizing selection. R. longicauda males provide females with a nuptial gift before engaging in aerial mating, and the male bears the entire weight of the female and nuptial gift for the duration of copulation. In theory, a male's ability to carry females and nuptial gifts could constrain pairing opportunities for the heaviest females, as reported for nonornamented dance flies. In concert with directional preferences for large females with mature eggs, such a load‐lifting constraint could produce the stabilizing selection on female size previously observed in this species. We therefore tested whether wild‐caught male R. longicauda collected during copulation were experiencing load‐lift limitations by comparing the mass carried by males during copulation with the male's wing loading traits. We also performed permutation tests to determine whether the loads carried by males during copulation were lighter than expected. We found that heavier males are more often found mating with heavier females suggesting that whereas R. longicauda males do not experience a load‐lift constraint, there is a strong relationship of assortative mating by mass. We suggest that active male mate choice for intermediately adorned females is more likely to be causing the nonrandom mating patterns observed in R. longicauda.     

Molecular identification of natural hybridization between Melastoma malabathricum and Melastoama beccarianum in Sarawak,Malaysia

Hybridization is very common in flowering plants and it plays a significant role in plant evolution and adaptation. Melastoma L. (Melastomataceae) comprises about 80–90 species in tropical Asia and Oceania, among which 41 species occur in Borneo. Natural hybridization is frequently reported in Melastoma in China, but so far there have been no confirmed cases of hybridization in Southeast Asia (including Borneo), where most species occur. Here, we identified a case of natural hybridization between Melastoma malabathricum L. and Melastoma beccarianum Cogn. in Sarawak, Malaysia, by using sequence data of three nuclear genes and one chloroplast intergenic spacer. Melastoma malabathricum is the most widespread species of this genus, occurring in almost the whole range of this genus, while M. beccarianum is a local species endemic to northern Borneo. Our results showed that natural hybridization and introgression occur between M. malabathricum and M. beccarianum, and the introgression was asymmetrical, mainly from M. malabathricum to M. beccarianum. As adaptive traits can be transferred by introgression, our study suggests that natural hybridization should be a significant mechanism for the evolution and adaptation of Melastoma in Southeast Asia. However, introgression from the common species M. malabathricum to the relatively rare species M. beccarianum may cause the decline of M. beccarianum, incurring conservation concern. With a large number of species of Melastoma and almost year‐around flowering in Southeast Asia, more cases of natural hybridization are expected to be found and identified in near future.     

Detection,Identification and Phylogenetic Analysis of Indian Ralstonia solanacearum Isolates by Comparison of Partial hrpB Gene Sequences

A species‐specific Polymerase Chain Reaction (sPCR) method was developed to identify and detect isolates of Ralstonia solanacearum, the cause of bacterial wilt disease in chilli. PCR primers for R. solanacearum were identified by alignment of hrpB gene sequences and selection of sequences specific for R. solanacearum at their 3′ ends. The primers were shown to be specific for R. solanacearum, as no PCR product was obtained when genomic DNA from other bacterial species including closely related Ralstonia species, were used as test species. Lone pair of primers (RshrpBF and RshrpBR) was designed using hrpB gene sequence, unique to R. solanacearum which amplified a predicted PCR product of 810 bp from 20 different isolates. Phylogenetic analysis was also attempted to understand the evolutionary divergence of Indian R. solanacearum isolates. Based on phylogenetic analysis, Indian isolates showed homology with the standard reference isolates from other countries but, interestingly, one new isolate showed complete evolutionary divergence by forming an out‐group.     

Phylogeography and lineage‐specific patterns of genetic diversity and molecular evolution in a group of North American skinks

Geography influences the evolutionary trajectory of species by mediating opportunities for hybridization, gene flow, demographic shifts and adaptation. We sought to understand how geography and introgression can generate species‐specific patterns of genetic diversity by examining phylogeographical relationships in the North American skink species Plestiodon multivirgatus and P. tetragrammus (Squamata: Scincidae). Using a multilocus dataset (three mitochondrial genes, four nuclear genes; a total of 3455 bp) we discovered mito‐nuclear discordance, consistent with mtDNA introgression. We further tested for evidence of species‐wide mtDNA introgression by using comparisons of genetic diversity, selection tests and extended Bayesian skyline analyses. Our findings suggest that P. multivirgatus acquired its mitochondrial genome from P. tetragrammus after their initial divergence. This putative species‐wide mitochondrial capture was further evidenced by statistically indistinguishable substitution rates between mtDNA and nDNA in P. multivirgatus. This rate discrepancy was observed in P. multivirgatus but not P. tetragrammus, which has important implications for studies that combine mtDNA and nDNA sequences when inferring time since divergence between taxa. Our findings suggest that by facilitating opportunities for interspecific introgression, geography can alter the course of molecular evolution between recently diverged lineages.     

Rapid evolution and the genomic consequences of selection against interspecific mating

While few species introduced into a new environment become invasive, those that do provide critical information on ecological mechanisms that determine invasions success and the evolutionary responses that follow invasion. Aedes albopictus (the Asian tiger mosquito) was introduced into the naturalized range of Aedes aegypti (the yellow fever mosquito) in the United States in the mid‐1980s, resulting in the displacement of A. aegypti in much of the south‐eastern United States. The rapid displacement was likely due to the superior competitive ability of A. albopictus as larvae and asymmetric mating interference competition, in which male A. albopictus mate with and sterilize A. aegypti females, a process called “satyrization.” The goal of this study was to examine the genomic responses of a resident species to an invasive species in which the mechanism of character displacement is understood. We used double‐digest restriction enzyme DNA sequencing (ddRADseq) to analyse outlier loci between selected and control lines of laboratory‐reared A. aegypti females from two populations (Tucson, AZ and Key West, Florida, USA), and individual females classified as either “resisted” or “mated with” A. albopictus males via mating trials of wild‐derived females from four populations in Florida. We found significant outlier loci in comparing selected and control lines and between mated and nonmated A. aegypti females in the laboratory and wild‐derived populations, respectively. We found overlap in specific outlier loci between different source populations that support consistent genomic signatures of selection within A. aegypti. Our results point to regions of the A. aegypti genome and potential candidate genes that may be involved in mating behaviour, and specifically in avoiding interspecific mating choices.     

Genetic variation in Plethodon cinereus and Plethodon hubrichti from in and around a contact zone

Climate change poses several challenges to biological communities including changes in the frequency of encounters between closely related congeners as a result of range shifts. When climate change leads to increased hybridization, hybrid dysfunction or genetic swamping may increase extinction risk—particularly in range‐restricted species with low vagility. The Peaks of Otter Salamander, Plethodon hubrichti, is a fully terrestrial woodland salamander that is restricted to ~18 km of ridgeline in the mountains of southwestern Virginia, and its range is surrounded by the abundant and widespread Eastern Red‐backed Salamander, Plethodon cinereus. In order to determine whether these two species are hybridizing and how their range limits may be shifting, we assessed variation at eight microsatellite loci and a 1,008 bp region of Cytochrome B in both species at allopatric reference sites and within a contact zone. Our results show that hybridization between P. hubrichti and P. cinereus either does not occur or is very rare. However, we find that diversity and differentiation are substantially higher in the mountaintop endemic P. hubrichti than in the widespread P. cinereus, despite similar movement ability for the two species as assessed by a homing experiment. Furthermore, estimation of divergence times between reference and contact zone populations via approximate Bayesian computation is consistent with the idea that P. cinereus has expanded into the range of P. hubrichti. Given the apparent recent colonization of the contact zone by P. cinereus, future monitoring of P. cinereus range limits should be a priority for the management of P. hubrichti populations.     

Genetic structure of Micromeria (Lamiaceae) in Tenerife,the imprint of geological history and hybridization on within‐island diversification

Geological history of oceanic islands can have a profound effect on the evolutionary history of insular flora, especially in complex islands such as Tenerife in the Canary Islands. Tenerife results from the secondary connection of three paleo‐islands by a central volcano, and other geological events that further shaped it. This geological history has been shown to influence the phylogenetic history of several taxa, including genus Micromeria (Lamiaceae). Screening 15 microsatellite markers in 289 individuals representing the eight species of Micromeria present in Tenerife, this study aims to assess the genetic diversity and structure of these species and its relation with the geological events on the island. In addition, we evaluate the extent of hybridization among species and discuss its influence on the speciation process. We found that the species restricted to the paleo‐islands present lower levels of genetic diversity but the highest levels of genetic differentiation suggesting that their ranges might have contracted over time. The two most widespread species in the island, M. hyssopifolia and M. varia, present the highest genetic diversity levels and a genetic structure that seems correlated with the geological composition of the island. Samples from M. hyssopifolia from the oldest paleo‐island, Adeje, appear as distinct while samples from M. varia segregate into two main clusters corresponding to the paleo‐islands of Anaga and Teno. Evidence of hybridization and intraspecific migration between species was found. We argue that species boundaries would be retained despite hybridization in response to the habitat's specific conditions causing postzygotic isolation and preserving morphological differentiation.     

Does selection on horn length of males and females differ in protected and hunted populations of a weakly dimorphic ungulate?

Weaponry in ungulates may be costly to grow and maintain, and different selective pressures in males and females may lead to sex‐biased natural survival. Sexual differences in the relationship between weapon growth and survival may increase under anthropogenic selection through culling, for example because of trophy hunting. Selection on weaponry growth under different scenarios has been largely investigated in males of highly dimorphic ungulates, for which survival costs (either natural or hunting related) are thought to be greatest. Little is known, however, about the survival costs of weaponry in males and females of weakly dimorphic species. We collected information on horn length and age at death/shooting of 407 chamois Rupicapra rupicapra in a protected population and in two hunted populations with different hunting regimes, to explore sexual differences in the selection on early horn growth under contrasting selective pressures. We also investigated the variation of horn growth and body mass in yearling males (n = 688) and females (n = 539) culled in one of the hunted populations over 14 years. The relationship between horn growth and survival showed remarkable sexual differences under different evolutionary scenarios. Within the protected population, under natural selection, we found no significant trade‐off in either males or females. Under anthropogenic pressure, selection on early horn growth of culled individuals showed diametrically opposed sex‐biased patterns, depending on the culling regime and hunters’ preferences. Despite the selective bias between males and females in one of the hunted populations, we did not detect significant sex‐specific differences in the long‐term pattern of early growth. The relationship between early horn growth and natural survival in either sex might suggest stabilizing selection on horn size in chamois. Selection through culling can be strongly sex‐biased also in weakly dimorphic species, depending on hunters’ preferences and hunting regulations, and long‐term data are needed to reveal potential undesirable evolutionary consequences.     

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.