The biogeography of mitochondrial and nuclear discordance in animals

Geography and past climate changes have shaped the evolution of a widespread lizard in arid Central Asia

The complex orogenic history and structure of Central Asia, coupled with Pleistocene glacial cycles have generated its stepwise aridification. Such events would have significantly influenced the evolution of many mid-latitude species in arid Central Asia (ACA). In this study, we employed two mitochondrial genes (CO1 and ND2) and genome-wide SNPs, coupled with ecological niche modeling, to investigate the lineage diversification and historical demography within a widespread lizard Phrynocepahlus helioscopus, and their associations with geography and past climate change. We obtained themtDNA dataset for 300 individuals from 96 localities within the known range of the lizard, among which 51 individuals from 27 localities were selected for generating the SNP dataset via genotyping-by-sequencing approach. Phylogenetic analyses of the concatenated mtDNA dataset revealed eight geographically correlated lineages that diverged by 4.21-10.41% for the CO1 gene, which were estimated to have coalesced ∼4.47 million years ago. However, we observed mito-nuclear discordance pattern regarding the population of Clade V (P. helioscopus sergeevi) from Tajikistan. Ancestral area estimations suggested that P. helioscopus originated from the Fergana Valley and then dispersed into the adjacent areas in ACA along with a history of multiple allopatric divergence processes, suggesting that Fergana may have been the cradle of diversification of P. helioscopus. The intensification of aridification across Central Asia during the Late Pliocene may have facilitated the rapid radiation of this arid-adapted lizard throughout this vast territory. Subsequently, the geological events (e.g., uplift of the Hissar-Alay, transgressions of the Caspian Sea) and geographic barriers (e.g., Amu Darya River, Zerarshan River) during the Pleistocene triggered the progressive diversification of P. helioscopus. Interestingly, Clade VIII (P. helioscopus varius) experienced rapid population growth coupled with range expansion while Clade IV (P. helioscopus cameranoi) underwent drastic population expansion associated with range contraction during the Last Glacial Maximum. In Clade IV, but not in Clade VIII, environmental turnover contributes more to mitochondrial genetic distinctiveness than geographic distance. Overall, the SNP dataset demonstrates that geographic distance plays a greater role than environmental distance. Both the mtDNA dataset and the SNP dataset suggest local-scale genetic differentiation in Clade IV and Clade VIII, revealing potential geographic barriers in the Ili River Valley and the Junggar Basin, respectively. Twenty-seven outlier SNPs associated with environmental factors (precipitation and temperature) were identified, which supports the signature of local adaptation to the arid desert environment. Finally, our finding suggests taxonomic implications, such as support for full species status for P. saidalievi (Clade II) and P. meridionalis (Clade I). Future analyses based on further evidence and increased taxon and geographic sampling should be carried out to corroborate our findings.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Morphology, nuclear SNPs and mate selection reveal that COI barcoding overestimates species diversity in a Mediterranean freshwater amphipod by an order of magnitude

DNA sequence information has revealed many morphologically cryptic species worldwide. For animals, DNA-based assessments of species diversity usually rely on the mitochondrial cytochrome c oxidase subunit I (COI) gene. However, a growing amount of evidence indicate that mitochondrial markers alone can lead to misleading species diversity estimates due to mito-nuclear discordance. Therefore, reports of putative species based solely on mitochondrial DNA should be verified by other methods, especially in cases where COI sequences are identical for different morphospecies or where divergence within the same morphospecies is high. Freshwater amphipods are particularly interesting in this context because numerous putative cryptic species have been reported. Here, we investigated the species status of the numerous mitochondrial molecular operational taxonomic units (MOTUs) found within Echinogammarus sicilianus. We used an integrative approach combining DNA barcoding with mate selection observations, detailed morphometrics and genome-wide double digest restriction site-associated DNA sequencing (ddRAD-seq). Within a relatively small sampling area, we detected twelve COI MOTUs (divergence = 1.8-20.3%), co-occurring in syntopy at two-thirds of the investigated sites. We found that pair formation was random and there was extensive nuclear gene flow among the ten MOTUs co-occurring within the same river stretch. The four most common MOTUs were also indistinguishable with respect to functional morphology. Therefore, the evidence best fits the hypothesis of a single, yet genetically diverse, species within the main river system. The only two MOTUs sampled outside the focal area were genetically distinct at the nuclear level and may represent distinct species. Our study reveals that COI-based species delimitation can significantly overestimate species diversity, highlighting the importance of integrative taxonomy for species validation, especially in hyperdiverse complexes with syntopically occurring mitochondrial MOTUs.

Multiple introgression events during the diversification history of the edible mexican grasshopper genus sphenarium (orthoptera: pyrgomorphidae)

Speciation with gene flow often leads to ambiguous phylogenetic reconstructions, reticulate patterns of relatedness and conflicting nuclear versus mitochondrial (mt) lineages. Here we employed a fragment of the COI mtDNA gene and nuclear genome-wide data (3RAD) to assess the diversification history of Sphenarium, an orthopteran genus of great economic importance in Mexico that is presumed to have experienced hybridisation events in some of its species. We carried out separate phylogenetic analyses to evaluate the existence of mito-nuclear discordance in the species relationships, and also assessed the genomic diversity and population genomic structure and investigated the existence of interspecific introgression and species limits of the taxa involved based on the nuclear dataset. The species delineation analyses discriminated all the currently recognised species, but also supported the existence of four undescribed species. The mt and nuclear topologies had four discordant species relationships that can be explained by mt introgression, where the mt haplotypes of S. purpurascens appear to have replaced those of S. purpurascens A and B, S. variabile and S. zapotecum. Moreover, our analyses supported the existence of nuclear introgression events between four species pairs that are distributed in the Sierra Madre del Sur province in southeast Mexico, with three of them occurring in the Tehuantepec Isthmus region. Our study highlights the relevance of genomic data to address the relative importance of allopatric isolation versus gene flow in speciation.

Temporal stability of the hybrid zone between Calocitta magpie-jays revealed through comparison of museum specimens and iNaturalist photos

Hybrid zones are natural experiments for the study of avian evolution. Hybrid zones can be dynamic, moving as species adjust to new climates and habitats, with unknown implications for species and speciation. There are relatively few studies that have comparable modern and historic sampling to assess change in hybrid zone location and width over time, and those studies have generally found mixed results, with many hybrid zones showing change over time, but others showing stability. The white-throated magpie-jay (Calocitta formosa) and black-throated magpie-jay (Calocitta colliei) occur along the western coast of Mexico and Central America. The two species differ markedly in throat color and tail length, and prior observation suggests a narrow hybrid zone in southern Jalisco where individuals have mixed throat color. This study aims to assess the existence and temporal stability of this putative hybrid zone by comparing throat color between georeferenced historical museum specimens and modern photos from iNaturalist with precise locality information. Our results confirm the existence of a narrow hybrid zone in Jalisco, with modern throat scores gradually increasing from the parental ends of the cline toward the cline center in a sigmoidal curve characteristic of hybrid zones. Our temporal comparison suggests that the hybrid zone has not shifted its position between historical (pre-1973) and modern (post-2005) time periods-a surprising result given the grand scale of habitat change to the western Mexican lowlands during this time. An anomalous pocket of white-throated individuals in the northern range of the black-throated magpie-jay hints at the possibility of prehistorical long-distance introduction. Future genomic data will help disentangle the evolutionary history of these lineages and better characterize how secondary contact is affecting both the DNA and the phenotype of these species.

Multiple introgressions shape mitochondrial evolutionary history in Drosophila paulistorum and the Drosophila willistoni group

Hybridization and the consequent introgression of genomic elements is an important source of genetic diversity for biological lineages. This is particularly evident in young clades in which hybrid incompatibilities are still incomplete and mixing between species is more likely to occur. Drosophila paulistorum, a representative of the Neotropical Drosophila willistoni subgroup, is a classic model of incipient speciation. The species is divided into six semispecies that show varying degrees of pre- and post-mating incompatibility with each other. In the present study, we investigate the mitochondrial evolutionary history of D. paulistorum and the willistoni subgroup. For that, we perform phylogenetic and comparative analyses of the complete mitochondrial genomes and draft nuclear assemblies of 25 Drosophila lines of the willistoni and saltans species groups. Our results show that the mitochondria of D. paulistorum are polyphyletic and form two non-sister clades that we name α and β. Identification and analyses of nuclear mitochondrial insertions further reveal that the willistoni subgroup has an α-like mitochondrial ancestor and strongly suggest that both the α and β mitochondria of D. paulistorum were acquired through introgression from unknown fly lineages of the willistoni subgroup. We also uncover multiple mitochondrial introgressions across D. paulistorum semispecies and generate novel insight into the evolution of the species.

New Insights into the Taxonomy of Malacopsylloidea Superfamily (Siphonaptera) Based on Morphological, Molecular and Phylogenetic Characterization of Phthiropsylla agenoris (Malacopsyllidae) and Polygenis (Polygenis) rimatus (Rhopalopsyllidae)

From a phylogenetic point of view, the Malacopsyllidae family and the Rhopalopsillidae family (comprising Parapsyllinae and Rhopalopsyllinae subfamilies) have been traditionally classified within the Malacopsylloidea superfamily, mostly restricted to South America. The phylogenetic relationships and taxonomic status of Malacopsyllidae and Rhopalopsillidae have never been assessed since no molecular loci of Malacopsyllidae have been sequenced by any authors, and the phylogeny provided so far was not based on any sort of formal quantitative analysis of flea morphology. Based on these precedents, the objective of this study was to carry out a comparative phylogenetic, molecular and morphological study of two different species belonging to each family, Phthiropsylla agenoris (Malacopsylla) and Polygenis (Polygenis) rimatus (Rhopalopsyllidae, Rhopalopsyllinae). In this study, we demonstrated the usefulness of several morphological features as diagnostic characters to differentiate between P. (P.) rimatus and P. agenoris. Using molecular and phylogenetic data, we easily discriminated between the two taxa (P. agenoris and P. (P.) rimatus) by comparing both nuclear and mitochondrial markers. This fact proves the usefulness of ITS2, EF1−α, cox1, cytb and cox2 as molecular diagnostic markers to characterize and identify different Siphonaptera taxa. Additionally, the phylogenetic results confirm, for the first time, the monophyly of the Malacopsyllidae family and suggest a clear paraphyletic position of the Paraspsyllinae subfamily and, consequently, the Rhopalopsyllidae family.

Mitochondriomics of Clarias Fishes (Siluriformes: Clariidae) with a New Assembly of Clarias camerunensis: Insights into the Genetic Characterization and Diversification

The mitogenome of an endemic catfish Clarias camerunensis was determined from the Cameroon water. This circular mitogenome was 16,511 bp in length and comprised 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a single AT-rich control region. The heavy strand accommodates 28 genes, whereas the light strand is constituted by ND6 and eight transfer RNA (tRNA) genes. The C. camerunensis mitochondrial genome is AT biased (56.89%), as showcased in other Clarias species. The comparative analyses revealed that most of the Clarias species have 6 overlapping and 11 intergenic spacer regions. Most of the PCGs were initiated and terminated with the ATG start codon and TAA stop codon, respectively. The tRNAs of C. camerunensis folded into the distinctive cloverleaf secondary structure, except trnS1. The placement of the conserved domains in the control region was similar in all the Clarias species with highly variable nucleotides in CSB-I. Both maximum likelihood and Bayesian-based matrilineal phylogenies distinctly separated all Clarias species into five clades on the basis of their known distributions (South China, Sundaland, Indochina, India, and Africa). The TimeTree analysis revealed that the two major clades (Indo-Africa and Asia) of Clarias species might have diverged during the Paleogene (≈28.66 MYA). Our findings revealed the separation of Indian species (C. dussumieri) and African species (C. camerunensis and Clarias gariepinus) took place during the Paleogene, as well as the South Chinese species (Clarias fuscus) and Sundaland species (Clarias batrachus) splits from the Indochinese species (Clarias macrocephalus) during the Neogene through independent colonization. This pattern of biotic relationships highlights the influence of topography and geological events in determining the evolutionary history of Clarias species. The enrichment of mitogenomic data and multiple nuclear loci from their native range or type locality will confirm the true diversification of Clarias species in African and Asian countries.

Range-wide whole-genome resequencing of the brown bear reveals drivers of intraspecies divergence

Population-genomic studies can shed new light on the effect of past demographic processes on contemporary population structure. We reassessed phylogeographical patterns of a classic model species of postglacial recolonisation, the brown bear (Ursus arctos), using a range-wide resequencing dataset of 128 nuclear genomes. In sharp contrast to the erratic geographical distribution of mtDNA and Y-chromosomal haplotypes, autosomal and X-chromosomal multi-locus datasets indicate that brown bear population structure is largely explained by recent population connectivity. Multispecies coalescent based analyses reveal cases where mtDNA haplotype sharing between distant populations, such as between Iberian and southern Scandinavian bears, likely results from incomplete lineage sorting, not from ancestral population structure (i.e., postglacial recolonisation). However, we also argue, using forward-in-time simulations, that gene flow and recombination can rapidly erase genomic evidence of former population structure (such as an ancestral population in Beringia), while this signal is retained by Y-chromosomal and mtDNA, albeit likely distorted. We further suggest that if gene flow is male-mediated, the information loss proceeds faster in autosomes than in X chromosomes. Our findings emphasise that contemporary autosomal genetic structure may reflect recent population dynamics rather than postglacial recolonisation routes, which could contribute to mtDNA and Y-chromosomal discordances.

Species delimitation using genomic data to resolve taxonomic uncertainties in a speciation continuum of pelagic seabirds

Speciation is a continuous and complex process shaped by the interaction of numerous evolutionary forces. Despite the continuous nature of the speciation process, the implementation of conservation policies relies on the delimitation of species and evolutionary significant units (ESUs). Puffinus shearwaters are globally distributed and threatened pelagic seabirds. Due to remarkable morphological status the group has been under intense taxonomic debate for the past three decades. Here, we use double digest Restriction-Site Associated DNA sequencing (ddRAD-Seq) to genotype species and subspecies of North Atlantic and Mediterranean Puffinus shearwaters across their entire geographical range. We assess the phylogenetic relationships and population structure among and within the group, evaluate species boundaries, and characterise the genomic landscape of divergence. We find that current taxonomies are not supported by genomic data and propose a more accurate taxonomy by integrating genomic information with other sources of evidence. Our results show that several taxon pairs are at different stages of a speciation continuum. Our study emphasises the potential of genomic data to resolve taxonomic uncertainties, which can help to focus management actions on relevant taxa, even if they do not necessarily coincide with the taxonomic rank of species.

Complex cycles of divergence and migration shape lineage structure in the common kingsnake species complex

Aim

The Nearctic is a complex patchwork of habitats and geologic features that form barriers to gene flow resulting in phylogeographic structure and speciation in many lineages. Habitats are rarely stable over geologic time, and the Nearctic has undergone major climatic changes in the past few million years. We use the common kingsnake species complex to study how climate, geography, and history influence lineage formation over a large, complex landscape.

Location

Nearctic/North America.

Taxon

Common kingsnake, Lampropeltis getula, species complex.

Methods

We analyzed genome-wide sequence data from 51 snakes spanning the majority of the species complex's range. We used population clustering, generalized dissimilarity modeling, and coalescent methods to identify the number of genetic clusters within the L. getula complex, infer the environmental correlates of genetic differentiation, and estimate models of divergence and gene flow among lineages.

Results

We identified three major lineages within the L. getula complex and further continuous spatial structure within lineages. The most important ecological correlates of genetic distance in the complex are related to aridity and precipitation, consistent with lineage breaks at the Great Plains/Desert ecotone and the Cochise Filter Barrier. Lineages are estimated to have undergone multiple rounds of isolation and secondary contact, with highly asymmetric migration occurring at present.

Main conclusions

Changing climates combined with a large and geologically complex landscape have resulted in a mosaic of discrete and spatially continuous genetic structure. Multiple rounds of isolation and secondary contact as climate fluctuated over the past ~4.4 My have likely driven the evolution of discrete lineages that maintain high levels of gene flow. Continuous structure is strongly shaped by aridity and precipitation, suggesting roles for major precipitation gradients in helping to maintain lineage identity in the face of gene flow when lineages are in geographic contact.

Molecular taxonomy of phlebotomine sand flies (Diptera, Psychodidae) with emphasis on DNA barcoding: A review

The taxonomy and systematics of sand flies (Diptera, Psychodidae, Phlebotominae) are one of the pillars of research aimed to identifying vector populations and the agents transmitted by these insects. Traditionally, the use of morphological traits has been the main line of evidence for the definition of species, but the use of DNA sequences is useful as an integrative approach for their delimitation. Here, we discuss the current status of the molecular taxonomy of sand flies, including their most sequenced molecular markers and the main results. Only about 37% of all sand fly species have been processed for any molecular marker and are publicly available in the NCBI GenBank or BOLD Systems databases. The genera Phlebotomus, Nyssomyia, Psathyromyia and Psychodopygus are well-sampled, accounting for more than 56% of their sequenced species. However, less than 34% of the species of Sergentomyia, Lutzomyia, Trichopygomyia and Trichophoromyia have been sampled, representing a major gap in the knowledge of these groups. The most sequenced molecular markers are those within mtDNA, especially the DNA barcoding fragment of the cytochrome c oxidase subunit I (coi) gene, which has shown promising results in detecting cryptic diversity within species. Few sequences of conserved genes have been generated, which hampers higher-level phylogenetic inferences. We argue that sand fly species should be sequenced for at least the coi DNA barcoding marker, but multiple markers with different mutation rates should be assessed, whenever possible, to generate multilocus analysis.

Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe

Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.

Genetic evidence for widespread population size expansion in North American boreal birds prior to the Last Glacial Maximum

Pleistocene climate cycles are well documented to have shaped contemporary species distributions and genetic diversity. Northward range expansions in response to deglaciation following the Last Glacial Maximum (LGM; approximately 21 000 years ago) are surmised to have led to population size expansions in terrestrial taxa and changes in seasonal migratory behaviour. Recent findings, however, suggest that some northern temperate populations may have been more stable than expected through the LGM. We modelled the demographic history of 19 co-distributed boreal-breeding North American bird species from full mitochondrial gene sets and species-specific molecular rates. We used these demographic reconstructions to test how species with different migratory strategies were affected by glacial cycles. Our results suggest that effective population sizes increased in response to Pleistocene deglaciation earlier than the LGM, whereas genetic diversity was maintained throughout the LGM despite shifts in geographical range. We conclude that glacial cycles prior to the LGM have most strongly shaped contemporary genetic diversity in these species. We did not find a relationship between historic population dynamics and migratory strategy, contributing to growing evidence that major switches in migratory strategy during the LGM are unnecessary to explain contemporary migratory patterns.

A wide hybrid zone mediated by precipitation contributed to confused geographical structure of Scutiger boulengeri

Confused geographical structure of a population and mitonuclear discordance are shaped by a combination of rapid changes in population demographics and shifts in ecology. In this study, we generated a time-calibrated phylogeny of Scutiger boulengeri, an endemic Xizang alpine toad occurring in mountain streams on the Qinghai-Xizang (Tibet) Plateau (QTP). Based on three mitochondrial DNA (mtDNA) genes, eight clades were assigned to three deeply divergent lineages. Analysis of nuclear DNA (nuDNA) genes revealed three distinct clusters without geographic structure, indicating significantly high rates of gene flow. Coalescent theory framework analysis (approximate Bayesian computation model DIYABC and Migrate-N) suggested that divergence of the main intraspecific clusters was the result of hybridization after secondary contact in the Holocene around 0.59 million years ago (Ma). The ratio of mtDNA F _ST (fixation index) to nuDNA F _ST was 2.3, thus failing to show male-biased dispersal. Geographic cline analysis showed that a wide hybrid zone was initially established in southwestern China, without significant reproductive isolation but with strong introgression in S. boulengeri, suggesting high hybrid fitness. Furthermore, mtDNA genes exhibited isolation by distance (IBD) while nuDNA genes exhibited significant isolation by environment (IBE). Results suggested that mitonuclear discordance may have initially been caused by geographic isolation, followed by precipitation-mediated hybridization, producing a wide hybrid zone and geographic structure confusion of nuDNA genes in S. boulengeri. This study indicated that complicated historical processes may have led to specific genetic patterns, with a specific climate factor facilitating gene flow in the system.

Contrasting Phylogeographic Patterns of Mitochondrial and Genome-Wide Variation in the Groundwater Amphipod Crangonyx islandicus That Survived the Ice Age in Iceland

The analysis of phylogeographic patterns has often been based on mitochondrial DNA variation, but recent analyses dealing with nuclear DNA have in some instances revealed mito-nuclear discordances and complex evolutionary histories. These enigmatic scenarios, which may involve stochastic lineage sorting, ancestral hybridization, past dispersal and secondary contacts, are increasingly scrutinized with a new generation of genomic tools such as RADseq, which also poses additional analytical challenges. Here, we revisited the previously inconclusive phylogeographic history, showing the mito-nuclear discordance of an endemic groundwater amphipod from Iceland, Crangonyx islandicus, which is the only metazoan known to have survived the Pleistocene beneath the glaciers. Previous studies based on three DNA markers documented a mitochondrial scenario with the main divergence occurring between populations in northern Iceland and an ITS scenario with the main divergence between the south and north. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to clarify this mito-nuclear discordance by applying several statistical methods while estimating the sensitivity to different analytical approaches (data-type, differentiation indices and base call uncertainty). A majority of nuclear markers and methods support the ITS divergence. Nevertheless, a more complex scenario emerges, possibly involving introgression led by male-biased dispersal among northern locations or mitochondrial capture, which may have been further strengthened by natural selection.

Genomic evidence for three distinct species in the Erebia manto complex in Central Europe (Lepidoptera, Nymphalidae)

A problem to implement conservation strategies is that in many cases recognized taxa are in fact complexes of several cryptic species. Failure to properly delineate species may lead to misplaced priorities or to inadequate conservation measures. One such species complex is the yellow-spotted ringlet Erebia manto, which comprises several phenotypically distinct lineages, whose degree of genomic isolation has so far not been assessed. Some of these lineages are geographically restricted and thus possibly represent distinct units with conservation priorities. Using several thousand nuclear genomic markers, we evaluated to which degree the bubastis lineage from the Alps and the vogesiaca lineage from the Vosges, are genetically isolated from the widespread manto lineage. Our results suggest that both lineages are genetically as strongly differentiated from manto as other taxonomically well separated sibling species in this genus from each other, supporting a delineation of bubastis and vogesiaca as independent species. Given the restricted and isolated range of vogesiaca as well as the disjunct distribution of bubastis, our findings have significant implication for future conservation efforts on these formerly cryptic species and highlight the need to investigate the genomic identity within species complexes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-023-01501-w.

A Phylogenomic Assessment of Processes Underpinning Convergent Evolution in Open-Habitat Chats

Insights into the processes underpinning convergent evolution advance our understanding of the contributions of ancestral, introgressed, and novel genetic variation to phenotypic evolution. Phylogenomic analyses characterizing genome-wide gene tree heterogeneity can provide first clues about the extent of ILS and of introgression and thereby into the potential of these processes or (in their absence) the need to invoke novel mutations to underpin convergent evolution. Here, we were interested in understanding the processes involved in convergent evolution in open-habitat chats (wheatears of the genus Oenanthe and their relatives). To this end, based on whole-genome resequencing data from 50 taxa of 44 species, we established the species tree, characterized gene tree heterogeneity, and investigated the footprints of ILS and introgression within the latter. The species tree corroborates the pattern of abundant convergent evolution, especially in wheatears. The high levels of gene tree heterogeneity in wheatears are explained by ILS alone only for 30% of internal branches. For multiple branches with high gene tree heterogeneity, D-statistics and phylogenetic networks identified footprints of introgression. Finally, long branches without extensive ILS between clades sporting similar phenotypes provide suggestive evidence for the role of novel mutations in the evolution of these phenotypes. Together, our results suggest that convergent evolution in open-habitat chats involved diverse processes and highlight that phenotypic diversification is often complex and best depicted as a network of interacting lineages.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Evaluating a species phylogeny using ddRAD SNPs: Cyto-nuclear discordance and introgression in the salmonid genus Thymallus (Salmonidae)

Hybridization and introgression are very common among freshwater fishes due to the dynamic nature of hydrological landscapes. Cyclic patterns of allopatry and secondary contact provide numerous opportunities for interspecific gene flow, which can lead to discordant paths of evolution for mitochondrial and nuclear genomes. Here, we used double digest restriction-site associated DNA sequencing (ddRADseq) to obtain a genome-wide single nucleotide polymorphism (SNP) dataset comprehensive for allThymallus (Salmonidae)species to infer phylogenetic relationships and evaluate potential recent and historical gene flow among species. The newly obtained nuclear phylogeny was largely concordant with a previously published mitogenome-based topology but revealed a few cyto-nuclear discordances. These incongruencies primarily involved the placement of internal nodes rather than the resolution of species, except for one European species where anthropogenic stock transfers are thought to be responsible for the observed pattern. The analysis of four contact zones where multiple species are found revealed a few cases of mitochondrial capture and limited signals of nuclear introgression. Interestingly, the mechanisms restricting interspecific gene flow might be distinct; while in zones of secondary contact, small-scale physical habitat separation appeared as a limiting factor, biologically based reinforcement mechanisms are presumed to be operative in areas where species presumably evolved in sympatry. Signals of historical introgression were largely congruent with the routes of species dispersal previously inferred from mitogenome data. Overall, the ddRADseq dataset provided a robust phylogenetic reconstruction of the genus Thymallus including new insights into historical hybridization and introgression, opening up new questions concerning their evolutionary history.

The potential and shortcomings of mitochondrial DNA analysis for cheetah conservation management

There are only about 7,100 adolescent and adult cheetahs (Acinonyx jubatus) remaining in the wild. With the majority occurring outside protected areas, their numbers are rapidly declining. Evidence-based conservation measures are essential for the survival of this species. Genetic data is routinely used to inform conservation strategies, e.g., by establishing conservation units (CU). A commonly used marker in conservation genetics is mitochondrial DNA (mtDNA). Here, we investigated the cheetah's phylogeography using a large-scale mtDNA data set to refine subspecies distributions and better assign individuals to CUs. Our dataset mostly consisted of historic samples to cover the cheetah's whole range as the species has been extinct in most of its former distribution. While our genetic data largely agree with geography-based subspecies assignments, several geographic regions show conflicting mtDNA signals. Our analyses support previous findings that evolutionary forces such as incomplete lineage sorting or mitochondrial capture likely confound the mitochondrial phylogeography of this species, especially in East and, to some extent, in Northeast Africa. We caution that subspecies assignments solely based on mtDNA should be treated carefully and argue for an additional standardized nuclear single nucleotide polymorphism (SNP) marker set for subspecies identification and monitoring. However, the detection of the A. j. soemmeringii specific haplogroup by a newly designed Amplification-Refractory Mutation System (ARMS) can already provide support for conservation measures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-022-01483-1.

High heterogeneity in genomic differentiation between phenotypically divergent songbirds: a test of mitonuclear co-introgression

Comparisons of genomic variation among closely related species often show more differentiation in mitochondrial DNA (mtDNA) and sex chromosomes than in autosomes, a pattern expected due to the differing effective population sizes and evolutionary dynamics of these genomic components. Yet, introgression can cause species pairs to deviate dramatically from general differentiation trends. The yellowhammer (Emberiza citrinella) and pine bunting (E. leucocephalos) are hybridizing avian sister species that differ greatly in appearance and moderately in nuclear DNA, but that show no mtDNA differentiation. This discordance is best explained by adaptive mtDNA introgression-a process that can select for co-introgression at nuclear genes with mitochondrial functions (mitonuclear genes). To better understand these discordant differentiation patterns and characterize nuclear differentiation in this system, we investigated genome-wide differentiation between allopatric yellowhammers and pine buntings and compared it to what was seen previously in mtDNA. We found significant nuclear differentiation that was highly heterogeneous across the genome, with a particularly wide differentiation peak on the sex chromosome Z. We further investigated mitonuclear gene co-introgression between yellowhammers and pine buntings and found support for this process in the direction of pine buntings into yellowhammers. Genomic signals indicative of co-introgression were common in mitonuclear genes coding for subunits of the mitoribosome and electron transport chain complexes. Such introgression of mitochondrial DNA and mitonuclear genes provides a possible explanation for the patterns of high genomic heterogeneity in genomic differentiation seen among some species groups.

Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group)

Discordance between different genomic regions, often identified through multilocus sequencing of selected markers, presents particular difficulties in identifying historical processes which drive species diversity and boundaries. Mechanisms causing discordance, such as incomplete lineage sorting or introgression due to interspecific hybridization, are better identified based on population-level genomic datasets. In the toads of the Rhinella granulosa species group, patterns of mito-nuclear discordance and potential hybridization have been reported by several studies. However, these patterns were proposed based on few loci, such that alternative mechanisms behind gene-tree heterogeneity cannot be ruled out. Using genome-wide ddRADseq loci from a subset of species within this clade, we found only partial concordance between currently recognized species-level taxon boundaries and patterns of genetic structure. While most taxa within the R. granulosa group correspond to clades, genetic clustering analyses sometimes grouped distinct taxonomic units into a single cluster. Moreover, levels of admixture between inferred clusters were limited and restricted to a single taxon pair which is best explained by incomplete lineage sorting as opposed to introgressive hybridization, according to D-statistics results. These findings contradict previous assertions of widespread cryptic diversity and gene flow within the R. granulosa clade. Lastly, our analyses suggest that diversification events within the Rhinella granulosa group mostly dated back to the early Pliocene, being generally younger than species divergences in other closely related clades that present high levels of cross-species gene flow. This finding uniquely contradicts common assertions that this young clade of toads exhibits interspecific hybridization.

Evolution of mitochondrial and nuclear genomes in Pennatulacea

We examine the phylogeny of sea pens using sequences of whole mitochondrial genomes and the nuclear ribosomal cluster generated through low coverage Illumina sequencing. Taxon sampling includes 30 species in 19 genera representing 13 families. Ancestral state reconstruction shows that most sea pen mitochondrial genomes have the ancestral gene order, and that Pennatulacea with diverse gene orders are found in a single clade. The monophyly of Pennatulidae and Protoptilidae are rejected by both the mitochondrial and nuclear dataset, while the mitochondrial dataset further rejects monophyly of Virgulariidae, and the nuclear dataset rejects monophyly of Kophobelemnidae. We show discordance between nuclear ribosomal gene cluster phylogenies and whole mitochondrial genome phylogenies and highlight key Pennatulacea taxa that could be included in cnidarian genome-wide studies to better resolve the sea pen tree of life. We further illustrate how well frequently sequenced markers capture the overall diversity of the mitochondrial genome and the nuclear ribosomal genes in sea pens.

Whole-genome Analysis Reveals Contrasting Relationships Among Nuclear and Mitochondrial Genomes Between Three Sympatric Bat Species

Understanding mechanisms involved in speciation can be challenging, especially when hybridization or introgression blurs species boundaries. In bats, resolving relationships of some closely related groups has proved difficult due subtle interspecific variation both in morphometrics and molecular data sets. The endemic South American Histiotus bats, currently considered a subgenus of Eptesicus, harbor unresolved phylogenetic relationships and of those is a trio consisting of two closely related species: Eptesicus (Histiotus) macrotus and Eptesicus (Histiotus) montanus, and their relationship with a third, Eptesicus (Histiotus) magellanicus. The three sympatric species bear marked resemblance to each other, but can be differentiated morphologically. Furthermore, previous studies have been unable to differentiate the species from each other at a molecular level. In order to disentangle the phylogenetic relationships of these species, we examined the differentiation patterns and evolutionary history of the three Eptesicus (H.) species at the whole-genome level. The nuclear DNA statistics between the species suggest strong gene flow and recent hybridization between E. (H.) montanus and E. (H.) macrotus, whereas E. (H.) magellanicus shows a higher degree of isolation. In contrast, mitochondrial DNA shows a closer relationship between E. (H.) magellanicus and E. (H.) montanus. Opposing patterns in mtDNA and nuclear markers are often due to differences in dispersal, and here it could be both as a result of isolation in refugia during the last glacial maximum and female philopatry and male-biased dispersal. In conclusion, this study shows the importance of both the nuclear and mitochondrial DNA in resolving phylogenetic relationships and species histories.

New insights on patterns of genetic admixture and phylogeographic history in Iberian high mountain populations of midwife toads

Multiple Quaternary glacial refugia in the Iberian Peninsula, commonly known as "refugia within refugia", allowed diverging populations to come into contact and admix, potentially boosting substantial mito-nuclear discordances. In this study, we employ a comprehensive set of mitochondrial and nuclear markers to shed light onto the drivers of geographical differentiation in Iberian high mountain populations of the midwife toads Alytes obstetricans and A. almogavarii from the Pyrenees, Picos de Europa and Guadarrama Mountains. In the three analysed mountain regions, we detected evidence of extensive mito-nuclear discordances and/or admixture between taxa. Clustering analyses identified three major divergent lineages in the Pyrenees (corresponding to the eastern, central and central-western Pyrenees), which possibly recurrently expanded and admixed during the succession of glacial-interglacial periods that characterised the Late Pleistocene, and that currently follow a ring-shaped diversification pattern. On the other hand, populations from the Picos de Europa mountains (NW Iberian Peninsula) showed a mitochondrial affinity to central-western Pyrenean populations and a nuclear affinity to populations from the central Iberian Peninsula, suggesting a likely admixed origin for Picos de Europa populations. Finally, populations from the Guadarrama Mountain Range (central Iberian Peninsula) were depleted of genetic diversity, possibly as a consequence of a recent epidemic of chytridiomycosis. This work highlights the complex evolutionary history that shaped the current genetic composition of high mountain populations, and underscores the importance of using a multilocus approach to better infer the dynamics of population divergence.

Genomic Analyses Implicate the Amazon-Orinoco Plume as the Driver of Cryptic Speciation in a Swimming Crab

The Amazon-Orinoco plume (AOP) is the world's largest freshwater and sediment discharge into the ocean. Previous studies limited to mtDNA suggest that the swimming crab Callinectes ornatus Ordway, 1863 exists as two distinct genetic clusters separated by the AOP. However, questions concerning migration, diversification time, and species delimitation are unresolved. Densely sampling markers across the genome (SNPs) could elucidate the evolutionary processes within this species. Here, we combined mtDNA data and ddRAD-seq to explore the diversification patterns and processes within the swimming crab C. ornatus. We show great genetic differentiation between groups on the north and south sides of the plume but also signs of hybridization. Demographic modeling indicates the divergence between groups starting around 8 Mya following the AOP's formation. After a period of isolation, we detect two incidences of secondary contact with stronger migration in concordance with the North Brazil Current flow. Our results suggest speciation with gene flow explained by the interplay among the AOP, oceanographic currents, and long larval dispersal. This work represents the first investigation employing ddRAD-seq in a marine invertebrate species with distribution encompassing the north and south Atlantic and sheds light on the role of the AOP in the diversification of a marine species.

Multiple invasions, Wolbachia and human-aided transport drive the genetic variability of Aedes albopictus in the Iberian Peninsula

The Asian tiger mosquito, Aedes albopictus, is one of the most invasive species in the world. Native to the tropical forests of Southeast Asia, over the past 30 years it has rapidly spread throughout tropical and temperate regions of the world. Its dramatic expansion has resulted in public health concerns as a consequence of its vector competence for at least 16 viruses. Previous studies showed that Ae. albopictus spread has been facilitated by human-mediated transportation, but much remains unknown about how this has affected its genetic attributes. Here we examined the factors that contributed to shaping the current genetic constitution of Ae. albopictus in the Iberian Peninsula, where the species was first found in 2004, by combining population genetics and Bayesian modelling. We found that both mitochondrial and nuclear DNA markers showed a lack of genetic structure and the presence of worldwide dominant haplotypes, suggesting regular introductions from abroad. Mitochondrial DNA showed little genetic diversity compared to nuclear DNA, likely explained by infection with maternally transmitted bacteria of the genus Wolbachia. Multilevel models revealed that greater mosquito fluxes (estimated from commuting patterns and tiger mosquito population distribution) and spatial proximity between sampling sites were associated with lower nuclear genetic distance, suggesting that rapid short- and medium-distance dispersal is facilitated by humans through vehicular traffic. This study highlights the significant role of human transportation in shaping the genetic attributes of Ae. albopictus and promoting regional gene flow, and underscores the need for a territorially integrated surveillance across scales of this disease-carrying mosquito.

Characterising Mitochondrial Capture in an Iberian Shrew

Mitochondrial introgression raises questions of biogeography and of the extent of reproductive isolation and natural selection. Previous phylogenetic work on the Sorex araneus complex revealed apparent mitonuclear discordance in Iberian shrews, indicating past hybridisation of Sorex granarius and the Carlit chromosomal race of S. araneus, enabling introgression of the S. araneus mitochondrial genome into S. granarius. To further study this, we genetically typed 61 Sorex araneus/coronatus/granarius from localities in Portugal, Spain, France, and Andorra at mitochondrial, autosomal, and sex-linked loci and combined our data with the previously published sequences. Our data are consistent with earlier data indicating that S. coronatus and S. granarius are the most closely related of the three species, confirming that S. granarius from the Central System mountain range in Spain captured the mitochondrial genome from a population of S. araneus. This mitochondrial capture event can be explained by invoking a biogeographical scenario whereby S. araneus was in contact with S. granarius during the Younger Dryas in central Iberia, despite the two species currently having disjunct distributions. We discuss whether selection favoured S. granarius with an introgressed mitochondrial genome. Our data also suggest recent hybridisation and introgression between S. coronatus and S. granarius, as well as between S. araneus and S. coronatus.

Complete mitochondrial genomes and updated divergence time of the two freshwater clupeids endemic to Lake Tanganyika (Africa) suggest intralacustrine speciation

Background

The hydrogeological history of Lake Tanganyika paints a complex image of several colonization and adaptive radiation events. The initial basin was formed around 9–12 million years ago (MYA) from the predecessor of the Malagarasi–Congo River and only 5–6 MYA, its sub-basins fused to produce the clear, deep waters of today. Next to the well-known radiations of cichlid fishes, the lake also harbours a modest clade of only two clupeid species, Stolothrissatanganicae and Limnothrissamiodon. They are members of Pellonulini, a tribe of clupeid fishes that mostly occur in freshwater and that colonized West and Central-Africa during a period of high sea levels during the Cenozoic. There is no consensus on the phylogenetic relationships between members of Pellonulini and the timing of the colonization of Lake Tanganyika by clupeids.

Results

We use short-read next generation sequencing of 10X Chromium libraries to sequence and assemble the full mitochondrial genomes of S.tanganicae and L.miodon. We then use Maximum likelihood and Bayesian inference to place them into the phylogeny of Pellonulini and other clupeiforms, taking advantage of all available full mitochondrial clupeiform genomes. We identify Potamothrissaobtusirostris as the closest living relative of the Tanganyika sardines and confirm paraphyly for Microthrissa. We estimate the divergence of the Tanganyika sardines around 3.64 MYA [95% CI: 0.99, 6.29], and from P.obtusirostris around 10.92 MYA [95% CI: 6.37–15.48].

Conclusions

These estimates imply that the ancestor of the Tanganyika sardines diverged from a riverine ancestor and entered the proto-lake Tanganyika around the time of its formation from the Malagarasi–Congo River, and diverged into the two extant species at the onset of deep clearwater conditions. Our results prompt a more thorough examination of the relationships within Pellonulini, and the new mitochondrial genomes provide an important resource for the future study of this tribe, e.g. as a reference for species identification, genetic diversity, and macroevolutionary studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02085-8.

Genomic Signatures of Mitonuclear Coevolution in Mammals

Mitochondrial (mt) and nuclear-encoded proteins are integrated in aerobic respiration, requiring co-functionality among gene products from fundamentally different genomes. Different evolutionary rates, inheritance mechanisms, and selection pressures set the stage for incompatibilities between interacting products of the two genomes. The mitonuclear coevolution hypothesis posits that incompatibilities may be avoided if evolution in one genome selects for complementary changes in interacting genes encoded by the other genome. Nuclear compensation, in which deleterious mtDNA changes are offset by compensatory nuclear changes, is often invoked as the primary mechanism for mitonuclear coevolution. Yet, direct evidence supporting nuclear compensation is rare. Here, we used data from 58 mammalian species representing eight orders to show strong correlations between evolutionary rates of mt and nuclear-encoded mt-targeted (N-mt) proteins, but not between mt and non-mt-targeted nuclear proteins, providing strong support for mitonuclear coevolution across mammals. N-mt genes with direct mt interactions also showed the strongest correlations. Although most N-mt genes had elevated dN/dS ratios compared to mt genes (as predicted under nuclear compensation), N-mt sites in close contact with mt proteins were not overrepresented for signs of positive selection compared to noncontact N-mt sites (contrary to predictions of nuclear compensation). Furthermore, temporal patterns of N-mt and mt amino acid substitutions did not support predictions of nuclear compensation, even in positively selected, functionally important residues with direct mitonuclear contacts. Overall, our results strongly support mitonuclear coevolution across ∼170 million years of mammalian evolution but fail to support nuclear compensation as the major mode of mitonuclear coevolution.

Small islands and large biogeographic barriers have driven contrasting speciation patterns in Indo-Pacific sunbirds (Aves: Nectariniidae)

Birds of the Indo-Pacific have provided biologists with many foundationalinsights. This study presents evidence for strong phylogeographic structure in two sunbird species from the heart of this region, the olive-backed sunbird, Cinnyris jugularis, and the black sunbird, Leptocoma aspasia. We assessed population divergence using morphological, plumage, bioacoustic and molecular data (mitochondrial ND2/ND3). Our findings indicate that the olive-backed sunbird should be recognized as multiple species, because birds from Sulawesi and the Sahul Shelf are closely related to each other, but widely separated from those in other regions. In addition, we provide evidence for an endemic species on the Wakatobi Islands, an archipelago of deep-sea islands off south-east Sulawesi. That a small bird could exhibit a range all the way from Sulawesi to Australia, while diverging on a small archipelago within this range, illustrates the complex interplay between dispersal and speciation. Our black sunbird genetic data also suggest unrecognized population structure, despite relatively weak plumage divergence. Black sunbirds in Sulawesi are likely to be a separate species from those in New Guinea, with a mean genetic distance of 9.1%. Current taxonomy suggests these sunbird species transcend classic biogeographic barriers, but our results suggest that these barriers are not easily bypassed.

Mitochondrial phylogenomics of Acanthocephala: nucleotide alignments produce long-branch attraction artefacts

BACKGROUND

Classification of the Acanthocephala, a clade of obligate endoparasites, remains unresolved because of insufficiently strong resolution of morphological characters and scarcity of molecular data with a sufficient resolution. Mitochondrial genomes may be a suitable candidate, but they are available for a small number of species and their suitability for the task has not been tested thoroughly.

METHODS

Herein, we sequenced the first mitogenome for the large family Rhadinorhynchidae: Micracanthorhynchina dakusuiensis. These are also the first molecular data generated for this entire genus. We conducted a series of phylogenetic analyses using concatenated nucleotides (NUC) and amino acids (AAs) of all 12 protein-coding genes, three different algorithms, and the entire available acanthocephalan mitogenomic dataset.

RESULTS

We found evidence for strong compositional heterogeneity in the dataset, and Micracanthorhynchina dakusuiensis exhibited a disproportionately long branch in all analyses. This caused a long-branch attraction artefact (LBA) of M. dakusuiensis resolved at the base of the Echinorhynchida clade when the NUC dataset was used in combination with standard phylogenetic algorithms, maximum likelihood (ML) and Bayesian inference (BI). Both the use of the AA dataset (BI-AAs and ML-AAs) and the CAT-GTR model designed for suppression of LBA (CAT-GTR-AAs and CAT-GTR-NUC) at least partially attenuated this LBA artefact. The results support Illiosentidae as the basal radiation of Echinorhynchida and Rhadinorhynchidae forming a clade with Echinorhynchidae and Pomporhynchidae. The questions of the monophyly of Rhadinorhynchidae and its sister lineage remain unresolved. The order Echinorhynchida was paraphyletic in all of our analyses.

CONCLUSIONS

Future studies should take care to attenuate compositional heterogeneity-driven LBA artefacts when applying mitogenomic data to resolve the phylogeny of Acanthocephala.

The emergence of a cryptic lineage and cytonuclear discordance through past hybridization in the Japanese fire-bellied newt, Cynops pyrrhogaster (Amphibia: Urodela)

Discrepancies in geographic variation patterns between nuclear DNA and mitochondrial DNA (mtDNA) are the result of the complicated differentiation processes in organisms and are the key to understanding their true evolutionary processes. The genetic differentiation of the northern and Southern-Izu lineages of the Japanese newt, Cynops pyrrhogaster, was investigated through their single nucleotide polymorphism variations obtained via multiplexed ISSR genotyping by sequencing (MIG-seq). We found three genetic groups (Tohoku, N-Kanto and S-Kanto), that were not detected by mtDNA variations, in the northern lineage. N-Kanto has intermediate genetic characteristics between Tohoku and S-Kanto. The genetic groups are now moderately isolated from each other and have unique genetic characteristics. An estimation of the evolutionary history using the approximate Bayesian computation (ABC) approach suggested that Tohoku diverged from the common ancestor of S-Kanto and S-Izu. Then, S-Kanto and S-Izu split, and the recent hybridization between Tohoku and S-Kanto gave rise to N-Kanto. The origin of N-Kanto through the hybridization is relatively young and seems to be related to changes in the distributions of Tohoku and S-Kanto as a result of climatic oscillation in the Pleistocene. We conclude that the mitochondrial genome of S-Kanto was captured in Tohoku and that the original mitochondrial genome of Tohoku was entirely removed through hybridization.

Subspecies at crossroads: the evolutionary significance of genomic and phenotypic variation in a wide-ranging Australian lizard (Ctenotus pantherinus)

Many subspecies were described to capture phenotypic variation in wide-ranging taxa, with some later being found to correspond to divergent genetic lineages. We investigate whether currently recognized subspecies correspond to distinctive and coherent evolutionary lineages in the widespread Australian lizard Ctenotus pantherinus based on morphological, mitochondrial and genome-wide nuclear variation. We find weak and inconsistent correspondence between morphological patterns and the presumed subspecies ranges, with character polymorphism within regions and broad morphological overlap across regions. Phylogenetic analyses suggest paraphyly of populations assignable to each subspecies, mitonuclear discordance and little congruence between subspecies ranges and the distribution of inferred clades. Genotypic clustering supports admixture across regions. These results undermine the presumed phenotypic and genotypic coherence and distinctiveness of C. pantherinus subspecies. Based on our findings, we comment on the operational and conceptual shortcomings of morphologically defined subspecies and discuss practical challenges in applying the general notion of subspecies as incompletely separated population lineages. We conclude by highlighting a historical asymmetry that has implications for ecology, evolution and conservation: subspecies proposed in the past are difficult to falsify even in the face of new data that challenge their coherence and distinctiveness, whereas modern researchers appear hesitant to propose new subspecies.

Speciation of pelagic zooplankton: Invisible boundaries can drive isolation of oceanic ctenophores

The study of evolution and speciation in non-model systems provides us with an opportunity to expand our understanding of biodiversity in nature. Connectivity studies generally focus on species with obvious boundaries to gene flow, but in open-ocean environments, such boundaries are difficult to identify. Due to the lack of obvious boundaries, speciation and population subdivision in the pelagic environment remain largely unexplained. Comb jellies (Phylum Ctenophora) are mostly planktonic gelatinous invertebrates, many of which are considered to have freely interbreeding distributions worldwide. It is thought that the lobate ctenophore Bolinopsis infundibulum is distributed throughout cooler northern latitudes and B. vitrea warmer. Here, we examined the global population structure for species of Bolinopsis with genetic and morphological data. We found distinct evolutionary patterns within the genus, where B. infundibulum had a broad distribution from northern Pacific to Atlantic waters despite many physical barriers, while other species were geographically segregated despite few barriers. Divergent patterns of speciation within the genus suggest that oceanic currents, sea-level, and geological changes over time can act as either barriers or aids to dispersal in the pelagic environment. Further, we used population genomic data to examine evolution in the open ocean of a distinct lineage of Bolinopsis ctenophores from the North Eastern Pacific. Genetic information and morphological observations validated this as a separate species, Bolinopsis microptera, which was previously described but has recently been called B. infundibulum. We found that populations of B. microptera from California were in cytonuclear discordance, which indicates a secondary contact zone for previously isolated populations. Discordance at this scale is rare, especially in a continuous setting.

Phylogenetic placement of Turkish populations of Ixodes ricinus and Ixodes inopinatus

Studies on phylogeography and population structure of Ixodes ricinus have been carried out in Europe for decades, but the number of specimens from the Middle East included in these analyses is relatively small, despite the wide distribution of the species in this area. This study aimed to clarify the phylogenetic positions of I. ricinus from Turkey as well as to investigate the presence of Ixodes inopinatus in Anatolia. For this purpose, one mitochondrial (mt 16S rDNA) and one nuclear gene (defensin) were used to generate molecular data from I. ricinus samples, which were collected from 17 locations across the species' distributional range in Turkey. Bayesian inference was used to investigate phylogenetic relationships. Globally, the mt 16S rDNA lineages correspond to the lineages revealed by defensin; I. ricinus and I. inopinatus sequences clustered separately. However, a discordant genetic pattern was observed between the phylogenetic position of turkish I. ricinus revealed by nuclear versus mitochondrial genes. All Turkish haplotypes of mt 16SrDNA clustered with I. ricinus samples from Europe, which might be the result of extensive gene flow between populations of Europe and the Middle East. On the other hand, a sample from Thrace Region grouped within I. inopinatus clade. Thus, the occurrence of I. inopinatus in Turkey was demonstrated for the first time using molecular data. Moreover, four individuals were found to be heterozygous for the defensin. The potential evolutionary processes that underlie this observed discrepancy between the phylogenetic trees of two genes have been discussed.

Phylogenomics reveals extensive introgression and a case of mito-nuclear discordance in the killifish genus Kryptolebias

Introgression is a widespread evolutionary process leading to phylogenetic inconsistencies among distinct parts of the genomes, particularly between mitochondrial and nuclear-based phylogenetic reconstructions (e.g., mito-nuclear discordances). Here, we used mtDNA and genome-wide nuclear sites to provide the first phylogenomic-based hypothesis on the evolutionary relationships within the killifish genus Kryptolebias. In addition, we tested for evidence of past introgression in the genus given the multiple reports of undergoing hybridization between its members. Our mtDNA phylogeny generally agreed with the relationships previously proposed for the genus. However, our reconstruction based on nuclear DNA revealed an unknown lineage - Kryptolebias sp. 'ESP' - as the sister group of the self-fertilizing mangrove killifishes, K. marmoratus and K. hermaphroditus. All individuals sequenced of Kryptolebias sp. 'ESP' had the same mtDNA haplotype commonly observed in K. hermaphroditus, demonstrating a clear case of mito-nuclear discordance. Our analysis further confirmed extensive history of introgression between Kryptolebias sp. 'ESP' and K. hermaphroditus. Population genomics analyses indicate no current gene flow between the two lineages, despite their current sympatry and history of introgression. We also confirmed introgression between other species pairs in the genus that have been recently reported to form hybrid zones. Overall, our study provides a phylogenomic reconstruction covering most of the Kryptolebias species, reveals a new lineage hidden in a case of mito-nuclear discordance, and provides evidence of multiple events of ancestral introgression in the genus. These findings underscore the importance of investigating different genomic information in a phylogenetic framework, particularly in taxa where introgression is common as in the sexually diverse mangrove killifishes.

Widespread genomic signatures of reproductive isolation and sex-specific selection in the Eastern Yellow Robin, Eopsaltria australis

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane’s Rule.

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.

Introgression at the emerging secondary contact zone of magpie Pica pica subspecies (Aves: Corvidae): integrating data on nuclear and mitochondrial markers, vocalizations, and field observations

Zones of secondary contact provide a good opportunity to investigate the origin and dynamics of reproductive isolation between related populations. We analyzed genetic and phenotypic patterns and gene flow between two subspecies of the Eurasian magpie Pica pica s.l. which recently came into contact after presumably long periods of isolation. We describe the distribution of subspecies in a young contact zone at Argun’ river basin in southern Siberia where populations occur in parapatry and an older hybrid population in eastern Mongolia. Based on genome-wide SNP data, we analyzed patterns and strength of gene flow between the subspecies. Our results indicate occasional hybridization with backcrossing and asymmetric introgression along a wide range in Transbaikalia and locally in eastern Mongolia. Males of P. p. jankowskii apparently exhibit higher dispersal ability towards the west compared to P. p. leucoptera (towards the east). The former occasionally migrates to eastern Mongolia and Transbaikalia where introgression of nuclear, but not mitochondrial DNA was evident. Bioacoustic investigations showed differences between the subspecies in speed and structure of vocalization. We discovered intermediate calls of hybrid magpies and bilingual birds alternating calls that are typical for the two taxa. Furthermore, we found dramatically decreased reproductive success in hybridogeneous populations. By complementing our results with established phylogeographic patterns of P. pica s.l. based on a mitochondrial marker sequence, and considering indications of sterility of hybrids in the contact zone, we propose to elevate the two corresponding subspecies to species level: P. pica for the western form and P. serica for the eastern form.