How reticulated are species?

Geographical isolation versus dispersal: Relictual alpine grasshoppers support a model of interglacial diversification with limited hybridization

Alpine biotas are paradigmatic of the countervailing roles of geographical isolation and dispersal during diversification. In temperate regions, repeated distributional shifts driven by Pleistocene climatic oscillations produced both recurrent pulses of population fragmentation and opportunities for gene flow during range expansions. Here, we test whether a model of divergence in isolation vs. with gene flow is more likely in the diversification of flightless alpine grasshoppers of the genus Podisma from the Iberian Peninsula. The answer to this question can also provide key insights about the pace of evolution. Specifically, if the data fit a divergence in isolation model, this suggests rapid evolution of reproductive isolation. Genomic data confirm a Pleistocene origin of the species complex, and multiple analytical approaches revealed limited asymmetric historical hybridization between two taxa. Genomic-based demographic reconstructions, spatial patterns of genetic structure and range shifts inferred from palaeodistribution modelling suggest severe range contraction accompanied by declines in effective population sizes during interglacials (i.e., contemporary populations confined to sky islands are relicts) and expansions during the coldest stages of the Pleistocene in each taxon. Although limited hybridization during secondary contact leads to phylogenetic uncertainty if gene flow is not accommodated when estimating evolutionary relationships, all species exhibit strong genetic cohesiveness. Our study lends support to the notion that the accumulation of incipient differences during periods of isolation were sufficient to lead to lineage persistence, but also that the demographic changes, dispersal constraints and spatial distribution of the sky islands themselves mediated species diversification in temperate alpine biotas.

Plotting for change: an analytical framework to aid decisions on which lineages are candidate species in phylogenomic species discovery

A recent study argued that coalescent-based models of species delimitation mostly delineate population structure, not species, and called for the validation of candidate species using biological information additional to the genetic information, such as phenotypic or ecological data. Here, we introduce a framework to interrogate genomic datasets and coalescent-based species trees for the presence of candidate species in situations where additional biological data are unavailable, unobtainable or uninformative. For de novo genomic studies of species boundaries, we propose six steps: (1) visualize genetic affinities among individuals to identify both discrete and admixed genetic groups from first principles and to hold aside individuals involved in contemporary admixture for independent consideration; (2) apply phylogenetic techniques to identify lineages; (3) assess diagnosability of those lineages as potential candidate species; (4) interpret the diagnosable lineages in a geographical context (sympatry, parapatry, allopatry); (5) assess significance of difference or trends in the context of sampling intensity; and (6) adopt a holistic approach to available evidence to inform decisions on species status in the difficult cases of allopatry. We adopt this approach to distinguish candidate species from within-species lineages for a widespread species complex of Australian freshwater fishes (Retropinna spp.). Our framework addresses two cornerstone issues in systematics that are often not discussed explicitly in genomic species discovery: diagnosability and how to determine it, and what criteria should be used to decide whether diagnosable lineages are conspecific or represent different species.

Genome-wide data reveal extensive gene flow during the diversification of the western rattlesnakes (Viperidae: Crotalinae: Crotalus)

Hybridization and introgression are important, but often overlooked processes when inferring phylogenies. When these processes are not accounted for and a strictly diverging phylogenetic model is applied to groups with a history of hybridization, phylogenetic inference and parameter estimation can be inaccurate. Recent developments in phylogenetic network approaches coupled with the increasing availability of genomic data allow inferences of reticulate evolutionary histories across the tree of life. The western rattlesnake species group (C. viridis species complex, C. mitchellii species complex, C. scutulutas, and C. tigris) is an iconic snake lineage that is widespread across western North America. This group is composed of several species complexes with unclear species limits, likely the result of ongoing gene flow among nascent lineages. Here I generate reduced representation genomic data and test for a history of reticulation within this group. I demonstrate that all species have undergone hybridization with at least one other lineage, suggesting introgression is widespread in this group. Topologies differ between phylogenies estimated under the multispecies coalescent and multispecies network coalescent methods, indicating that gene flow has obscured phylogenetic relationships within this group. These past introgression events are predominantly restricted to species that co-occur geographically. However, within species that have a history of introgression, this signature is detected regardless of specimen sampling across geography. Overall, my results suggest the accumulation of reproductive isolating barriers occurs slowly in rattlesnakes which likely leads to the difficulty in delimiting species, furthermore, the results of this study have implications for trait evolution in this group.

Most Genomic Loci Misrepresent the Phylogeny of an Avian Radiation Because of Ancient Gene Flow

Phylogenetic trees based on genome-wide sequence data may not always represent the true evolutionary history for a variety of reasons. One process that can lead to incorrect reconstruction of species phylogenies is gene flow, especially if interspecific gene flow has affected large parts of the genome. We investigated phylogenetic relationships within a clade comprising eight species of passerine birds (Phylloscopidae, Phylloscopus, leaf warblers) using one de novo genome assembly and 78 resequenced genomes. On the basis of hypothesis-exclusion trials based on D-statistics, phylogenetic network analysis, and demographic inference analysis, we identified ancient gene flow affecting large parts of the genome between one species and the ancestral lineage of a sister species pair. This ancient gene flow consistently caused erroneous reconstruction of the phylogeny when using large amounts of genome-wide sequence data. In contrast, the true relationships were captured when smaller parts of the genome were analyzed, showing that the "winner-takes-all democratic majority tree" is not necessarily the true species tree. Under this condition, smaller amounts of data may sometimes avoid the effects of gene flow due to stochastic sampling, as hidden reticulation histories are more likely to emerge from the use of larger data sets, especially whole-genome data sets. In addition, we also found that genomic regions affected by ancient gene flow generally exhibited higher genomic differentiation but a lower recombination rate and nucleotide diversity. Our study highlights the importance of considering reticulation in phylogenetic reconstructions in the genomic era.[Bifurcation; introgression; recombination; reticulation; Phylloscopus.].

Morphological ghosts of introgression in Darwin's finch populations

Many species of plants, animals, and microorganisms exchange genes well after the point of evolutionary divergence at which taxonomists recognize them as species. Genomes contain signatures of past gene exchange and, in some cases, they reveal a legacy of lineages that no longer exist. But genomic data are not available for many organisms, and particularly problematic for reconstructing and interpreting evolutionary history are communities that have been depleted by extinctions. For these, morphology may substitute for genes, as exemplified by the history of Darwin's finches on the Galápagos islands of Floreana and San Cristóbal. Darwin and companions collected seven specimens of a uniquely large form of Geospiza magnirostris in 1835. The populations became extinct in the next few decades, partly due to destruction of Opuntia cactus by introduced goats, whereas Geospiza fortis has persisted to the present. We used measurements of large samples of G. fortis collected for museums in the period 1891 to 1906 to test for unusually large variances and skewed distributions of beak and body size resulting from introgression. We found strong evidence of hybridization on Floreana but not on San Cristóbal. The skew is in the direction of the absent G. magnirostris We estimate introgression influenced 6% of the frequency distribution that was eroded by selection after G. magnirostris became extinct on these islands. The genetic residuum of an extinct species in an extant one has implications for its future evolution, as well as for a conservation program of reintroductions in extinction-depleted communities.

Phylogenomic assessment of the role of hybridization and introgression in trait evolution

Trait evolution among a set of species-a central theme in evolutionary biology-has long been understood and analyzed with respect to a species tree. However, the field of phylogenomics, which has been propelled by advances in sequencing technologies, has ushered in the era of species/gene tree incongruence and, consequently, a more nuanced understanding of trait evolution. For a trait whose states are incongruent with the branching patterns in the species tree, the same state could have arisen independently in different species (homoplasy) or followed the branching patterns of gene trees, incongruent with the species tree (hemiplasy). Another evolutionary process whose extent and significance are better revealed by phylogenomic studies is gene flow between different species. In this work, we present a phylogenomic method for assessing the role of hybridization and introgression in the evolution of polymorphic or monomorphic binary traits. We apply the method to simulated evolutionary scenarios to demonstrate the interplay between the parameters of the evolutionary history and the role of introgression in a binary trait's evolution (which we call xenoplasy). Very importantly, we demonstrate, including on a biological data set, that inferring a species tree and using it for trait evolution analysis in the presence of gene flow could lead to misleading hypotheses about trait evolution.

Deep ancestral introgression shapes evolutionary history of dragonflies and damselflies

Introgression is an important biological process affecting at least 10% of the extant species in the animal kingdom. Introgression significantly impacts inference of phylogenetic species relationships where a strictly binary tree model cannot adequately explain reticulate net-like species relationships. Here we use phylogenomic approaches to understand patterns of introgression along the evolutionary history of a unique, non-model insect system: dragonflies and damselflies (Odonata). We demonstrate that introgression is a pervasive evolutionary force across various taxonomic levels within Odonata. In particular, we show that the morphologically "intermediate" species of Anisozygoptera (one of the three primary suborders within Odonata besides Zygoptera and Anisoptera), which retain phenotypic characteristics of the other two suborders, experienced high levels of introgression likely coming from zygopteran genomes. Additionally, we find evidence for multiple cases of deep inter-superfamilial ancestral introgression.

Rapid radiation and rampant reticulation: Phylogenomics of South American Liolaemus lizards

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction-site associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.

Gene flow in phylogenomics: Sequence capture resolves species limits and biogeography of Afromontane forest endemic frogs from the Cameroon Highlands

Puddle frogs of the Phrynobatrachus steindachneri species complex are a useful group for investigating speciation and phylogeography in Afromontane forests of the Cameroon Volcanic Line, western Central Africa. The species complex is represented by six morphologically relatively cryptic mitochondrial DNA lineages, only two of which are distinguished at the species level - southern P. jimzimkusi and Lake Oku endemic P. njiomock, leaving the remaining four lineages identified as 'P. steindachneri'. In this study, the six mtDNA lineages are subjected to genomic sequence capture analyses and morphological examination to delimit species and to study biogeography. The nuclear DNA data (387 loci; 571,936 aligned base pairs) distinguished all six mtDNA lineages, but the topological pattern and divergence depths supported only four main clades: P. jimzimkusi, P. njiomock, and only two divergent evolutionary lineages within the four 'P. steindachneri' mtDNA lineages. One of the two lineages is herein described as a new species, P. amieti sp. nov. Reticulate evolution (hybridization) was detected within the species complex with morphologically intermediate hybrid individuals placed between the parental species in phylogenomic analyses, forming a ladder-like phylogenetic pattern. The presence of hybrids is undesirable in standard phylogenetic analyses but is essential and beneficial in the network multispecies coalescent. This latter approach provided insight into the reticulate evolutionary history of these endemic frogs. Introgressions likely occurred during the Middle and Late Pleistocene climatic oscillations, due to the cyclic connections (likely dominating during cold glacials) and separations (during warm interglacials) of montane forests. The genomic phylogeographic pattern supports the separation of the southern (Mt. Manengouba to Mt. Oku) and northern mountains at the onset of the Pleistocene. Further subdivisions occurred in the Early Pleistocene, separating populations from the northernmost (Tchabal Mbabo, Gotel Mts.) and middle mountains (Mt. Mbam, Mt. Oku, Mambilla Plateau), as well as the microendemic lineage restricted to Lake Oku (Mt. Oku). This unique model system is highly threatened as all the species within the complex have exhibited severe population declines in the past decade, placing them on the brink of extinction. In addition, Mount Oku is identified to be of particular conservation importance because it harbors three species of this complex. We, therefore, urge for conservation actions in the Cameroon Highlands to preserve their diversity before it is too late.

HybPhaser: A workflow for the detection and phasing of hybrids in target capture data sets

PREMISE

Hybrids contain divergent alleles that can confound phylogenetic analyses but can provide insights into reticulated evolution when identified and phased. We developed a workflow to detect hybrids in target capture data sets and phase reads into parental lineages using a similarity and phylogenetic framework.

METHODS

We used Angiosperms353 target capture data for Nepenthes, including known hybrids to test the novel workflow. Reference mapping was used to assess heterozygous sites across the data set and to detect hybrid accessions and paralogous genes. Hybrid samples were phased by mapping reads to multiple references and sorting reads according to similarity. Phased accessions were included in the phylogenetic framework.

RESULTS

All known Nepenthes hybrids and nine additional samples had high levels of heterozygous sites, had reads associated with multiple divergent clades, and were phased into accessions resembling divergent haplotypes. Phylogenetic analysis including phased accessions increased clade support and confirmed parental lineages of hybrids.

DISCUSSION

HybPhaser provides a novel approach to detect and phase hybrids in target capture data sets, which can provide insights into reticulations by revealing origins of hybrids and reduce conflicting signal, leading to more robust phylogenetic analyses.

Hybridization Outcomes Have Strong Genomic and Environmental Contingencies

AbstractExtreme F₂ phenotypes known as transgressive segregants can cause increased or decreased fitness in hybrids beyond the ranges seen in parental populations. Despite the usefulness of transgression for plant and animal breeding and its potential role in hybrid speciation, the genetic mechanisms and predictors of transgressive segregation remain largely untested. We generated seven hybrid crosses between five widely divergent Saccharomyces yeast species and measured the fitness of the parents and their viable F₁ and F₂ hybrids in seven stressful environments. We found that on average 16.6% of all replicate F₂ hybrids had higher fitness than both parents. Against our predictions, transgression frequency was not a function of parental genetic and phenotypic distances across test environments. Within environments, some relationships were significant, but not in the predicted direction; for example, genetic distance was negatively related to transgression in ethanol and hydrogen peroxide. Significant effects of hybrid cross, test environment, and cross × environment interactions suggest that the amount of transgression produced in a hybrid cross is highly context specific and that outcomes of hybridization differ even among crosses made from the same two parents. If the goal is to reliably predict hybrid fitness and forecast the evolutionary potential of admixed populations, we need more efforts to identify patterns beyond the idiosyncrasies caused by specific genomic or environmental contexts.

Hybrid Sterility, Genetic Conflict and Complex Speciation: Lessons From the Drosophila simulans Clade Species

The three fruitfly species of the Drosophila simulans clade- D. simulans, D. mauritiana, and D. sechellia- have served as important models in speciation genetics for over 40 years. These species are reproductively isolated by geography, ecology, sexual signals, postmating-prezygotic interactions, and postzygotic genetic incompatibilities. All pairwise crosses between these species conform to Haldane's rule, producing fertile F1 hybrid females and sterile F1 hybrid males. The close phylogenetic proximity of the D. simulans clade species to the model organism, D. melanogaster, has empowered genetic analyses of their species differences, including reproductive incompatibilities. But perhaps no phenotype has been subject to more continuous and intensive genetic scrutiny than hybrid male sterility. Here we review the history, progress, and current state of our understanding of hybrid male sterility among the D. simulans clade species. Our aim is to integrate the available information from experimental and population genetics analyses bearing on the causes and consequences of hybrid male sterility. We highlight numerous conclusions that have emerged as well as issues that remain unresolved. We focus on the special role of sex chromosomes, the fine-scale genetic architecture of hybrid male sterility, and the history of gene flow between species. The biggest surprises to emerge from this work are that (i) genetic conflicts may be an important general force in the evolution of hybrid incompatibility, (ii) hybrid male sterility is polygenic with contributions of complex epistasis, and (iii) speciation, even among these geographically allopatric taxa, has involved the interplay of gene flow, negative selection, and positive selection. These three conclusions are marked departures from the classical views of speciation that emerged from the modern evolutionary synthesis.

Genomic analyses overturn two long-standing homoploid hybrid speciation hypotheses

The importance of hybridization in generating biological diversity has been historically controversial. Previously, inference about hybridization was limited by dependence on morphological data; with the advent of the next-generation sequencing tools for nonmodel organisms, the evolutionary significance of hybridization is more evident. Here, we test classic hypotheses of hybrid origins of two species in the Phlox pilosa complex. Morphological intermediacy motivated the hypotheses that Phlox amoena lighthipei and Phlox pilosa deamii were independent homoploid hybrid lineages derived from P. amoena amoena and P. pilosa pilosa. We use double-digest restriction site-associated DNA sequencing of individuals from throughout the range of these taxa to conduct the most thorough analysis of evolutionary history in this system to date. Surprisingly, we find no support for the hybrid origin of P. pilosa deamii or P. amoena lighthipei. Our data do identify a history of admixture in individuals collected at a contemporary hybrid zone between the putative parent lineages. We show that three very different evolutionary histories, only one of which involves hybrid origin, have produced intermediate or recombinant morphological traits between P. amoena amoena and P. pilosa pilosa. Although morphological data are still an efficient means of generating hypotheses about past gene flow, genomic data are now the standard of evidence for elucidating evolutionary history.

Pure species discriminate against hybrids in the Drosophila melanogaster species subgroup

Introgression, the exchange of alleles between species, is a common event in nature. This transfer of alleles between species must happen through fertile hybrids. Characterizing the traits that cause defects in hybrids illuminates how and when gene flow is expected to occur. Inviability and sterility are extreme examples of fitness reductions but are not the only type of defects in hybrids. Some traits specific to hybrids are more subtle but are important to determine their fitness. In this report, we study whether F1 hybrids between two species pairs of Drosophila are as attractive as the parental species. We find that in both species pairs, the sexual attractiveness of the F1 hybrids is reduced and that pure species discriminate strongly against them. We also find that the cuticular hydrocarbon (CHC) profile of the female hybrids is intermediate between the parental species. Perfuming experiments show that modifying the CHC profile of the female hybrids to resemble pure species improves their chances of mating. Our results show that behavioral discrimination against hybrids might be an important component of the persistence of species that can hybridize.

Signatures of Introgression across the Allele Frequency Spectrum

The detection of introgression from genomic data is transforming our view of species and the origins of adaptive variation. Among the most widely used approaches to detect introgression is the so-called ABBA-BABA test or D-statistic, which identifies excess allele sharing between nonsister taxa. Part of the appeal of D is its simplicity, but this also limits its informativeness, particularly about the timing and direction of introgression. Here we present a simple extension, D frequency spectrum or DFS, in which D is partitioned according to the frequencies of derived alleles. We use simulations over a large parameter space to show how DFS carries information about various factors. In particular, recent introgression reliably leads to a peak in DFS among low-frequency derived alleles, whereas violation of model assumptions can lead to a lack of signal at low frequencies. We also reanalyze published empirical data from six different animal and plant taxa, and interpret the results in the light of our simulations, showing how DFS provides novel insights. We currently see DFS as a descriptive tool that will augment both simple and sophisticated tests for introgression, but in the future it may be usefully incorporated into probabilistic inference frameworks.

A Study of Hybridization Between Marmota baibacina and M. sibirica in Their Secondary Contact Zone in Mongolian Altai

The role of hybridization as one of the factors of speciation in mammals has been underestimated for a long time, but now there is a lot of data on its impact in mammalian evolution. Hybridization of species often occurs in their secondary contact zones, which is a natural model for testing factors that ensure species integrity. Studies of hybrid zones are increasingly revealing the essential role of ecological and behavioral features both in initiating crossbreeding and in maintaining interspecific barriers. We studied the hybridization of two species of marmots Marmota baibacina and M. sibirica in the zone of sympatry in Mongolian Altai Mountains. We used a bioacoustic approach to determine the localization of individuals of different species and their cohabitation sites. Genetic typing with two diploid nuclear markers and one marker each of paternal and maternal lines was used to identify hybrids. Habitat preferences of marmots were studied to understand the conditions for the formation of heterospecific pairs. We found a high proportion of hybrid individuals in boulder screes where conditions for the formation of heterospecific pairs probably exist. Our data indicate the viability and fertility of F1 hybrids and their descendants. We hypothesize that the environmental preferences and behavioral features of both species of marmots are important factors that both create conditions for hybridization and limit hybrid dispersal.

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

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

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

Extensive introgression and mosaic genomes of Mediterranean endemic lizards

The Mediterranean basin is a hotspot of biodiversity, fuelled by climatic oscillation and geological change over the past 20 million years. Wall lizards of the genus Podarcis are among the most abundant, diverse, and conspicuous Mediterranean fauna. Here, we unravel the remarkably entangled evolutionary history of wall lizards by sequencing genomes of 34 major lineages covering 26 species. We demonstrate an early (>11 MYA) separation into two clades centred on the Iberian and Balkan Peninsulas, and two clades of Mediterranean island endemics. Diversification within these clades was pronounced between 6.5-4.0 MYA, a period spanning the Messinian Salinity Crisis, during which the Mediterranean Sea nearly dried up before rapidly refilling. However, genetic exchange between lineages has been a pervasive feature throughout the entire history of wall lizards. This has resulted in a highly reticulated pattern of evolution across the group, characterised by mosaic genomes with major contributions from two or more parental taxa. These hybrid lineages gave rise to several of the extant species that are endemic to Mediterranean islands. The mosaic genomes of island endemics may have promoted their extraordinary adaptability and striking diversity in body size, shape and colouration, which have puzzled biologists for centuries.

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (d_XY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.

Evolutionary impacts of introgressive hybridization in a rapidly evolving group of jumping spiders (F. Salticidae, Habronattus americanus group)

Introgressive hybridization can be a powerful force impacting patterns of evolution at multiple taxonomic levels. We aimed to understand how introgression has affected speciation and diversification within a species complex of jumping spiders. The Habronattus americanus subgroup is a recently radiating group of jumping spiders, with species now in contact after hypothesized periods of isolation during glaciation cycles of the Pleistocene. Effects of introgression on genomes and morphology were investigated using phylogenomic and clustering methods using RADseq, ultraconserved elements (UCEs), and morphological data. We characterized 14 unique species/morphs using non-metric multidimensional scaling of morphological data, a majority of which were not recovered as monophyletic in our phylogenomic analyses. Morphological clusters and genetic lineages are highly incongruent, such that geographic region was a greater predictor of phylogenetic relatedness and genomic similarity than species or morph identity. STRUCTURE analyses support this pattern, revealing clusters corresponding to larger geographic regions. A history of rapid radiation in combination with frequent introgression seems to have mostly homogenized the genomes of species in this system, while selective forces maintain distinct male morphologies. GEMMA analyses support this idea by identifying SNPs correlated with distinct male morphologies. Overall, we have uncovered a system at odds with a typical bifurcating evolutionary model, instead supporting one where closely related species evolve together connected through multiple introgression events, creating a reticulate evolutionary history.

Mosquitoes of the Maculipennis complex in Northern Italy

The correct identification of mosquito vectors is often hampered by the presence of morphologically indiscernible sibling species. The Maculipennis complex is one of these groups that include both malaria vectors of primary importance and species of low/negligible epidemiological relevance, of which distribution data in Italy are outdated. Our study was aimed at providing an updated distribution of Maculipennis complex in Northern Italy through the sampling and morphological/molecular identification of specimens from five regions. The most abundant species was Anopheles messeae (2032), followed by Anopheles maculipennis s.s. (418), Anopheles atroparvus (28) and Anopheles melanoon (13). Taking advantage of ITS2 barcoding, we were able to finely characterize tested mosquitoes, classifying all the Anopheles messeae specimens as Anopheles daciae, a taxon with debated rank to which we referred as species inquirenda (sp. inq.). The distribution of species was characterized by Ecological Niche Models (ENMs), fed by recorded points of presence. ENMs provided clues on the ecological preferences of the detected species, with An. daciae sp. inq. linked to stable breeding sites and An. maculipennis s.s. more associated to ephemeral breeding sites. We demonstrate that historical Anopheles malaria vectors are still present in Northern Italy.

Biased Hybridization and Its Impact on Adaptive Introgression

Gene exchange between species can influence ecological and evolutionary processes ranging from population rescue to adaptive radiation. Genomic tools have provided new insights into the prevalence and nature of gene exchange between species. However, much remains unknown of how ecological, behavioral, and evolutionary factors determine what genetic variation moves between species in the first place. In particular, more research is needed that evaluates whether such factors bias gene flow from one species to another, and whether any such biases affect how genetic variation from another species is ultimately retained in the genome of a given species. Addressing this issue is crucial in a changing world where hybridization and introgression might determine which species succeed and which become extinct.

A Bayesian Implementation of the Multispecies Coalescent Model with Introgression for Phylogenomic Analysis

Recent analyses suggest that cross-species gene flow or introgression is common in nature, especially during species divergences. Genomic sequence data can be used to infer introgression events and to estimate the timing and intensity of introgression, providing an important means to advance our understanding of the role of gene flow in speciation. Here, we implement the multispecies-coalescent-with-introgression model, an extension of the multispecies-coalescent model to incorporate introgression, in our Bayesian Markov chain Monte Carlo program Bpp. The multispecies-coalescent-with-introgression model accommodates deep coalescence (or incomplete lineage sorting) and introgression and provides a natural framework for inference using genomic sequence data. Computer simulation confirms the good statistical properties of the method, although hundreds or thousands of loci are typically needed to estimate introgression probabilities reliably. Reanalysis of data sets from the purple cone spruce confirms the hypothesis of homoploid hybrid speciation. We estimated the introgression probability using the genomic sequence data from six mosquito species in the Anopheles gambiae species complex, which varies considerably across the genome, likely driven by differential selection against introgressed alleles.

Genomic evidence supports the introgression between two sympatric stickleback species inhabiting the White Sea basin

Interspecies hybridization is driven by a complex interplay of factors where introgression plays an important role. In the present study, the transfer of genetic material, between two quite distant fish species from different genera, through spontaneous hybridization was documented with dedicated molecular and bioinformatics tools. We investigate the genomic landscape of putative stickleback-relative introgression by carefully analyzing the tractable transposable elements (TE) on the admixed genome of some individuals of two sympatric stickleback species inhabiting northwestern Russia, namely the three-spined (Gasterosteus aculeatus) and the nine-spined (Pungitius pungitius) sticklebacks. Our data revealed that unique TE amplification types exist, supporting our proposed hypothesis that infers on the interspecific introgression. By running a restriction site-associated DNA sequencing (RAD-Seq) with eight samples of G. aculeatus and P. pungitius and subjecting further the results to a contrasting analysis by variated bioinformatic tools, we identified the related introgression-linked markers. The admixture nature observed in a single sample of the nine-spined stickleback demonstrated the possible traces of remote introgression between these two species. Our work reveals the potential that introgression has on providing particular variants at a high-frequency speed while linking blocks of sequence with multiple functional mutations. However, even though our results are of significant interest, an increased number of samples displaying the introgression are required to further ascertain our conclusions.

The number of species of degus (genus Octodon) is currently underestimated: An appraisal of species limits and their phylogenetic relationships (Rodentia: Hystricomorpha: Octodontidae)

As currently understood, the genus Octodon contains five species degus, lunatus, bridgesii, pacificus, and ricardojeda. Previous phylogenetic studies suggest that genus specific diversity is underestimated. In order to evaluate the taxonomic diversity of Octodon, we implemented unilocus (cytochrome-b) and multilocus (cytochrome-b + 4 nuclear genes) species delimitation methods. Octodon degus was recovered as a sister of the other species of the genus. The unilocus bGMYC and mPTP methods, based on cytochrome-b sequences, delimits 11 and 7 candidate species respectively, and both methods fail to recognize O. pacificus from O. ricardojeda. Results of the multilocus analysis (BPP) vary as a function of the dataset used. When the five genes are used 11 species are delimited, while eight species are delimited when only the nuclear genes are used. Octodon bridgesii is shown as comprising at least two species (one on the Pacific coast and the typical form found on the Andean slopes), while O. ricardojeda may comprise two species (one on the Chilean side of the Andes and the other in Argentina). Likewise, both multilocus matrices recover O. pacificus as a distinct species. This shows that species diversity of Octodon is underestimated. Remarkably, many of the delimited species based on genetic data are morphologically differentiated in cranio-dental characteristics. However, a pair of species has not achieved morphological differentiation, being cryptic species. Finally, the incongruence between mitochondrial and nuclear phylogenies suggests that processes such as incomplete lineage sorting and/or introgression have been present during the radiation of the genus.

A Comprehensive Approach to Detect Hybridization Sheds Light on the Evolution of Earth's Largest Lizards

Hybridization between species occurs more frequently in vertebrates than traditionally thought but distinguishing ancient hybridization from other phenomena that generate similar evolutionary patterns remains challenging. Here, we used a comprehensive workflow to discover evidence of ancient hybridization between the Komodo dragon (Varanus komodoensis) from Indonesia and a common ancestor of an Australian group of monitor lizards known colloquially as sand monitors. Our data comprises >300 nuclear loci, mitochondrial genomes, phenotypic data, fossil and contemporary records, and past/present climatic data. We show that the four sand monitor species share more nuclear alleles with V. komodoensis than expected given a bifurcating phylogeny, likely as a result of hybridization between the latter species and a common ancestor of sand monitors. Sand monitors display phenotypes that are intermediate between their closest relatives and V. komodoensis. Biogeographic analyses suggest that V. komodoensis and ancestral sand monitors co-occurred in northern Australia. In agreement with the fossil record, this provides further evidence that the Komodo dragon once inhabited the Australian continent. Our study shows how different sources of evidence can be used to thoroughly characterize evolutionary histories that deviate from a treelike pattern, that hybridization can have long-lasting effects on phenotypes and that detecting hybridization can improve our understanding of evolutionary and biogeographic patterns.

Relevance of ddRADseq method for species and population delimitation of closely related and widely distributed wolf spiders (Araneae, Lycosidae)

Although species delimitation is often controversial, emerging DNA-based and classical morphology-based methods are rarely compared using large-scale samplings, even less in the case of widely distributed species that have distant, allopatric populations. In the current study, we examined species boundaries within two wolf spider species of the genus Pardosa (Araneae, Lycosidae), P. riparia and P. palustris. Wolf spiders constitute an excellent model for testing the relevance of traditional vs. modern methods in species and population delimitation because several closely related species are distributed over cross-continental geographic ranges. Allopatric populations of the two Pardosa species were sampled across Europe to Far East Russia (latitudinal range > 150°) and several dozen individuals were studied using morphological characters (morphometry of three measures for both sexes, plus five in males only and two in females only), DNA barcoding (COI sequencing) and double-digest restriction site associated DNA sequencing (ddRADseq). The results obtained allow for changing the taxonomic status of two Far East Russian populations to subspecies and ddRADseq proved to be a powerful tool for taxonomic research despite scarce sampling and inherent subjectivity of species delimitation in allopatry. Overall, this study pleads for both multi-criteria and more population-based studies in taxonomy.

An integrative approach to address species limits in the southernmost members of the Liolaemus kingii group (Squamata: Liolaemini)

Recent conceptual and methodological advances have enabled an increasing number of studies to address the problem of species delimitation in a comprehensive manner. This is of particular interest in cases of species whose divergence times are recent and/or effective population sizes are large, where the conclusions obtained from a single source of evidence may lead to erroneous estimations of true species numbers or incorrect assignment of individuals to species. Iguanian lizards of the Liolaemus kingii group (13 species) comprise an important component of the endemic fauna of Patagonia. The southernmost species of this group (namely L. baguali, L. escarchadosi, L. sarmientoi, and L. tari) show widely overlapping distributions across southern Patagonia, also, their phylogenetic relationships are ambiguous and species boundaries have not been explicitly tested. Here we use a comprehensive approach to assess species limits through the use of molecular and morphological information (mitochondrial cytb, nuclear sequences collected by ddRADseq, and linear, meristic and landmark-based morphometrics). We found support for the current taxonomy given that the different analyses recognized the nominal species (4 entities), also a candidate species was supported by mitochondrial and morphological data. In addition, we detected signs of admixture between some of the species. Our results indicate that the L. kingii group can serve as a model system in studies of diversification accompanied by hybridization in nature, which in turn might have been promoted by past climatic oscillations and generalist morphologies. We emphasize the importance of using multiple lines of evidence in order to solve evolutionary stories, and minimizing potential erroneous results that may arise when relying on a single source of information.

The Gene and Gene Expression (GAGE) Species Concept: An Universal Approach for All Eukaryotic Organisms

The Gene and Gene Expression (GAGE) species concept, a new version of the Pragmatic Species Concept of Seifert (2014), is proposed as a concept applicable to any described recent or fossil eukaryotic organism independent from its mode of reproduction or evolutionary history. In addition to presenting the concept as such, the article also provides practical recommendations for taxonomists when delimiting species and describing taxa. The wording of the new concept contains a heading core sentence plus five attached sentences addressing essential conditions for its translation into a sound taxonomic practice: "Species are separable clusters that have passed a threshold of evolutionary divergence and are exclusively defined by nuclear DNA sequences and/or their expression products. Nuclear DNA sequences and their expression products are different character systems but have a highly correlated indicative function. Character systems with the least risk of epigenetic or ontogenetic modification have superior indicative value when conflicts between character systems of integrative studies arise. All character systems have to be described by an adequate numerics allowing cluster formation and determination of thresholds. Thresholds for each character system should be fixed by consensus among the experts under the principle of avoiding oversplitting or lumping. Clusters must not be the expression of intraspecific polymorphism." Recognizing the distortions and conflicts caused to taxonomy through barcoding or through assessment on the basis of association with other organisms, the GAGE species concept strongly downgrades the use of cytoplasmic DNA of endosymbiotic origin (mtDNA, cpDNA) or DNA of closely associated microbes (e.g., Wolbachia bacteria) for final taxonomic decision-making. Recognizing the distortion of phylogenies by the high frequency of reticulate evolution, it is argued that delimiting and naming species has to be separated from constructing bifurcating phylogenetic trees. [Cytoplasmic DNA; lumping; nuclear DNA; numeric taxonomy; oversplitting; reticulate evolution.].

Determining the probability of hemiplasy in the presence of incomplete lineage sorting and introgression

The incongruence of character states with phylogenetic relationships is often interpreted as evidence of convergent evolution. However, trait evolution along discordant gene trees can also generate these incongruences - a phenomenon known as hemiplasy. Classic comparative methods do not account for discordance, resulting in incorrect inferences about the number, timing, and direction of trait transitions. Biological sources of discordance include incomplete lineage sorting (ILS) and introgression, but only ILS has received theoretical consideration in the context of hemiplasy. Here, we present a model that shows introgression makes hemiplasy more likely, such that methods that account for ILS alone will be conservative. We also present a method and software (HeIST) for making statistical inferences about the probability of hemiplasy and homoplasy in large datasets that contain both ILS and introgression. We apply our methods to two empirical datasets, finding that hemiplasy is likely to contribute to the observed trait incongruences in both.

Cytonuclear discordance, reticulation and cryptic diversity in one of North America's most common frogs

Complicated phylogenetic histories benefit from diverse sources of inference. Pseudacris crucifer (spring peeper) spans most of eastern North America and comprises six mtDNA lineages that form multiple contact zones. The putative Miocene or early Pliocene origins of the oldest lineages within Pseudacris crucifer imply sufficient time for species-level divergence. To understand why this species appears unified while congeners have radiated, we analyze and compare male advertisement calls, mitochondrial, and nuclear markers and speak to the complex processes that have potentially influenced its contemporary patterns. We find extensive geographic and topological mitonuclear discordance, with three nuclear lineages containing 6 more-structured mtDNA lineages, and nuclear introgression at some contact zones. Male advertisement call differentiation is incongruent with the genetic structure as only one lineage appears differentiated. Occupying the Interior Highlands of the central United States, this Western lineage also has the most concordant mitochondrial and nuclear geographic patterns. Based on our findings we suggest that the antiquity of common ancestors was not as important as the maintenance of allopatry in the divergence in P. crucifer genetic lineages. We use multiple lines of evidence to generate hypotheses of isolation, reticulation, and discordance within this species and to expand our understanding of the early stages of speciation.

Hybrid Speciation and Introgression Both Underlie the Genetic Structures and Evolutionary Relationships of Three Morphologically Distinct Species of Lilium (Liliaceae) Forming a Hybrid Zone Along an Elevational Gradient

We studied hybrid interactions of Lilium meleagrinum, Lilium gongshanense, and Lilium saluenense using an integrative approach combining population genetics, fieldwork, and phenological research. These three species occur along an elevational gradient, with L. meleagrinum occurring at lower elevations, L. saluenense at higher elevations, and L. gongshanense between them. The species show strong morphological differentiation despite there being no clear environmental barriers to gene flow among them. Lilium gongshanense is likely to have a hybrid origin based on our prior work, but its progenitors remain uncertain. We sought to determine whether gene flow occurs among these three parapatric species, and, if so, whether L. gongshanense is a hybrid of L. meleagrinum and/or L. saluenense. We analyzed data from multiple chloroplast genes and spacers, nuclear internal transcribed spacer (ITS), and 18 nuclear Expressed Sequence Tag-Simple Sequence Repeat (EST-SSR) microsatellites for accessions of the three species representing dense population-level sampling. We also inferred phenology by examining species in the field and using herbarium specimens. We found that there are only two types of chloroplast genomes shared among the three species and that L. gongshanense forms two distinct groups with closest links to other species of Lilium based on ITS. Taken together, L. gongshanense is unlikely to be a hybrid species resulting from a cross between L. meleagrinum and L. saluenense, but gene flow is occurring among the three species. The gene flow is likely to be rare according to evidence from all molecular datasets, and this is corroborated by detection of only one putative hybrid individual in the field and asynchronous phenology. We suspect that the rarity of hybridization events among the species facilitates their continued genetic separation.

Primate phylogenomics uncovers multiple rapid radiations and ancient interspecific introgression

Our understanding of the evolutionary history of primates is undergoing continual revision due to ongoing genome sequencing efforts. Bolstered by growing fossil evidence, these data have led to increased acceptance of once controversial hypotheses regarding phylogenetic relationships, hybridization and introgression, and the biogeographical history of primate groups. Among these findings is a pattern of recent introgression between species within all major primate groups examined to date, though little is known about introgression deeper in time. To address this and other phylogenetic questions, here, we present new reference genome assemblies for 3 Old World monkey (OWM) species: Colobus angolensis ssp. palliatus (the black and white colobus), Macaca nemestrina (southern pig-tailed macaque), and Mandrillus leucophaeus (the drill). We combine these data with 23 additional primate genomes to estimate both the species tree and individual gene trees using thousands of loci. While our species tree is largely consistent with previous phylogenetic hypotheses, the gene trees reveal high levels of genealogical discordance associated with multiple primate radiations. We use strongly asymmetric patterns of gene tree discordance around specific branches to identify multiple instances of introgression between ancestral primate lineages. In addition, we exploit recent fossil evidence to perform fossil-calibrated molecular dating analyses across the tree. Taken together, our genome-wide data help to resolve multiple contentious sets of relationships among primates, while also providing insight into the biological processes and technical artifacts that led to the disagreements in the first place.

Combining three newly sequenced primate genomes with other published genomes, this study adapts a little-known method for detecting ancient introgression to genome-scale data, revealing multiple previously unknown examples of hybridization between primate species.

Radiation with reticulation marks the origin of a major malaria vector

Advances in genomics have led to an appreciation that introgression is common, but its evolutionary consequences are poorly understood. In recent species radiations the sharing of genetic variation across porous species boundaries can facilitate adaptation to new environments and generate novel phenotypes, which may contribute to further diversification. Most Anopheles mosquito species that are of major importance as human malaria vectors have evolved within recent and rapid radiations of largely nonvector species. Here, we focus on one of the most medically important yet understudied anopheline radiations, the Afrotropical Anopheles funestus complex (AFC), to investigate the role of introgression in its diversification and the possible link between introgression and vector potential. The AFC comprises at least seven morphologically similar species, yet only An. funestus sensu stricto is a highly efficient malaria vector with a pan-African distribution. Based on de novo genome assemblies and additional whole-genome resequencing, we use phylogenomic and population genomic analyses to establish species relationships. We show that extensive interspecific gene flow involving multiple species pairs has shaped the evolutionary history of the AFC since its diversification. The most recent introgression event involved a massive and asymmetrical movement of genes from a distantly related AFC lineage into An. funestus, an event that predated and plausibly facilitated its subsequent dramatic geographic range expansion across most of tropical Africa. We propose that introgression may be a common mechanism facilitating adaptation to new environments and enhancing vectorial capacity in Anopheles mosquitoes.

Natural variation in meiotic recombination rate shapes introgression patterns in intraspecific hybrids between wild and domesticated barley

Meiotic recombination rates vary considerably between species, populations and individuals. The genetic exchange between homologous chromosomes plays a major role in evolution by breaking linkage between advantageous and deleterious alleles in the case of introgressions. Identifying recombination rate modifiers is thus of both fundamental and practical interest to understand and utilize variation in meiotic recombination rates. We investigated recombination rate variation in a large intraspecific hybrid population (named HEB-25) derived from a cross between domesticated barley and 25 wild barley accessions. We observed quantitative variation in total crossover number with a maximum of a 1.4-fold difference between subpopulations and increased recombination rates across pericentromeric regions. The meiosis-specific α-kleisin cohesin subunit REC8 was identified as a candidate gene influencing crossover number and patterning. Furthermore, we quantified wild barley introgression patterns and revealed how local and genome-wide recombination rate variation shapes patterns of introgression. The identification of allelic variation in REC8 in combination with the observed changes in crossover patterning suggest a difference in how chromatin loops are tethered to the chromosome axis, resulting in reduced crossover suppression across pericentromeric regions. Local and genome-wide recombination rate variation is shaping patterns of introgressions and thereby directly influences the consequences of linkage drag.

Allopatry, hybridization, and reproductive isolation in Arctostaphylos

PREMISE

Hybridization plays a key role in introgressive adaptation, speciation, and adaptive radiation as a source of evolutionary innovation. Hybridization is considered common in Arctostaphylos, yet species boundaries are retained in stands containing multiple species. Arctostaphylos contains diploids and tetraploids, and recent phylogenies indicate two clades; we hypothesize combinations of these traits limit or promote hybridization rates.

METHODS

We statistically analyzed co-occurrence patterns of species by clade membership and ploidy level from 87 random 0.1 ha plots. We sampled multiple sites to analyze for percent hybridization based on morphology. Finally, phenophases were analyzed by scoring herbarium sheets for a large number of taxa from both clades as well as tetraploids, and second, surveying three field sites over two years for divergence in phenological stages between co-occurring taxa.

RESULTS

Most taxa in Arctostaphylos are allopatric relative to other congenerics. When two taxa co-occur, the patterns are a diploid with a tetraploid, or two diploids from different clades. When three taxa co-occur, the pattern is two diploids from different clades and a tetraploid. Field and herbarium data both indicate flowering phenology is displaced between diploids from the two clades; one of the diploid clades and tetraploids overlap considerably.

CONCLUSIONS

The two deep clades in Arctostaphylos are genetically distant, with hybrids rare or non-existent when taxa co-occur. Reproductive isolation between clades is enhanced by displaced flowering phenology for co-occurring species. Within clades, taxa appear to have few reproductive barriers other than an allopatric distribution or different ploidy levels.

Genetic Variation and Hybridization in Evolutionary Radiations of Cichlid Fishes

Evolutionary radiations are responsible for much of the variation in biodiversity across taxa. Cichlid fishes are well known for spectacular evolutionary radiations, as they have repeatedly evolved into large and phenotypically diverse arrays of species. Cichlid genomes carry signatures of past events and, at the same time, are the substrate for ongoing evolution. We survey genome-wide data and the available literature covering 438 cichlid populations (412 species) across multiple radiations to synthesize information about patterns and sharing of genetic variation. Nucleotide diversity within species is low in cichlids, with 92% of surveyed populations having less diversity than the median value found in other vertebrates. Divergence within radiations is also low, and a large proportion of variation is shared among species due to incomplete lineage sorting and widespread hybridization. Population genetics therefore provides a suitable conceptual framework for evolutionary genomic studies of cichlid radiations. We focus in detail on the roles of hybridization in shaping the patterns of genetic variation and in promoting cichlid diversification.

Phylogenomics and species delimitation for effective conservation of manta and devil rays

Practical biodiversity conservation relies on delineation of biologically meaningful units. Manta and devil rays (Mobulidae) are threatened worldwide, yet morphological similarities and a succession of recent taxonomic changes impede the development of an effective conservation strategy. Here, we generate genome-wide single nucleotide polymorphism (SNP) data from a geographically and taxonomically representative set of manta and devil ray samples to reconstruct phylogenetic relationships and evaluate species boundaries under the general lineage concept. We show that nominal species units supported by alternative data sources constitute independently evolving lineages, and find robust evidence for a putative new species of manta ray in the Gulf of Mexico. Additionally, we uncover substantial incomplete lineage sorting indicating that rapid speciation together with standing variation in ancestral populations has driven phylogenetic uncertainty within Mobulidae. Finally, we detect cryptic diversity in geographically distinct populations, demonstrating that management below the species level may be warranted in certain species. Overall, our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.

Gene flow between diploid and tetraploid junipers - two contrasting evolutionary pathways in two Juniperus populations

BACKGROUND

Gene flow and polyploidy have been found to be important in Juniperus evolution. However, little evidence has been published elucidating the association of both phenomena in juniper taxa in the wild. Two main areas were studied in Spain (Eastern Iberian Range and Sierra de Baza) with both diploid and tetraploid taxa present in sympatry. Gene flow and ploidy level were assessed for these taxa and the resulted offspring.

RESULTS

Twenty-two allo-triploid hybrids between J. sabina var. sabina and J. thurifera were found in the Eastern Iberian Range population. However, in the Sierra de Baza population no triploids were found. Instead, 18 allo-tetraploid hybrids between two tetraploid taxa: J. sabina var. balkanensis and J. thurifera were discovered. High genetic diversity was exhibited among the tetraploid hybrids at Sierra de Baza, in contrast to the genetically identical triploid hybrids at the Eastern Iberian Range; this suggests meiotic difficulties within the triploid hybrids. In addition, unidirectional gene flow was observed in both studied areas.

CONCLUSION

Polyploidy and hybridization can be complementary partners in the evolution of Juniperus taxa in sympatric occurrences. Juniperus was shown to be an ideal coniferous model to study these two phenomena, independently or in concert.

Effect of Whole-Genome Duplication on the Evolutionary Rescue of Sterile Hybrid Monkeyflowers

Hybridization is a creative evolutionary force, increasing genomic diversity and facilitating adaptation and even speciation. Hybrids often face significant challenges to establishment, including reduced fertility that arises from genomic incompatibilities between their parents. Whole-genome duplication in hybrids (allopolyploidy) can restore fertility, cause immediate phenotypic changes, and generate reproductive isolation. Yet the survival of polyploid lineages is uncertain, and few studies have compared the performance of recently formed allopolyploids and their parents under field conditions. Here, we use natural and synthetically produced hybrid and polyploid monkeyflowers (Mimulus spp.) to study how polyploidy contributes to the fertility, reproductive isolation, phenotype, and performance of hybrids in the field. We find that polyploidization restores fertility and that allopolyploids are reproductively isolated from their parents. The phenotype of allopolyploids displays the classic gigas effect of whole-genome duplication, in which plants have larger organs and are slower to flower. Field experiments indicate that survival of synthetic hybrids before and after polyploidization is intermediate between that of the parents, whereas natural hybrids have higher survival than all other taxa. We conclude that hybridization and polyploidy can act as sources of genomic novelty, but adaptive evolution is key in mediating the establishment of young allopolyploid lineages.

Genomic evidence suggests further changes of butterfly names

Further genomic sequencing of butterflies by our research group expanding the coverage of species and specimens from different localities, coupled with genome-scale phylogenetic analysis and complemented by phenotypic considerations, suggests a number of changes to the names of butterflies, mostly those recorded from the United States and Canada. Here, we present evidence to support these changes. The changes are intended to make butterfly classification more internally consistent at the genus, subgenus and species levels. I.e., considering all available evidence, we attempt to assign similar taxonomic ranks to the clades of comparable genetic differentiation, which on average is correlated with the age of phylogenetic groups estimated from trees. For species, we use criteria devised by genomic analysis of the genetic differentiation across suture zones and comparison of sympatric populations of closely related species. As a result, we resurrect 4 genera and 1 subgenus from subgeneric status or synonymy, change the rank of 8 currently used genera to subgenus, synonymize 7 genus-group names, summarize evidence to support 19 taxa as species instead of subspecies and 1 taxon as subspecies instead of species, along with a number of additional changes. One new genus and one new subspecies are described. Namely, the following taxa are treated as genera Tharsalea Scudder, 1876, Helleia Verity, 1943, Apangea Zhdanko, 1995, and Boldenaria Zhdanko, 1995. Tetracharis Grote, 1898 is a valid subgenus (not a synonym of Anthocharis Boisduval, Rambur, [Duménil] & Graslin, [1833]) that consists of Anthocharis cethura C. Felder & R. Felder, 1865 (Müller, 1764), Anthocharis midea (Hübner, [1809]), and Anthocharis limonea (A. Butler, 1871). The following are subgenera: Speyeria Scudder, 1872 of Argynnis Fabricius, 1807; Aglais Dalman, 1816 and Polygonia Hübner, [1819] of Nymphalis Kluk, 1780; Palaeonympha Butler, 1871 of Megisto Hübner, [1819]; Hyponephele Muschamp, 1915 of Cercyonis Scudder, 1875; Pyronia Hübner, [1819] and Aphantopus Wallengren, 1853 of Maniola Schrank, 1801 and Pseudonymphidia Callaghan, 1985 of Pachythone. Lafron Grishin, gen. n. (type species Papilio orus Stoll, [1780], parent subfamily Lycaeninae [Leach], [1815]) is described. Dipsas japonica Murray, 1875 is fixed as the type species of Neozephyrus Sibatani & Ito, 1942. The following taxa are junior subjective synonyms: Falcapica Klots, 1930 of Tetracharis Grote, 1898; Habrodais Scudder, 1876, Favonius Sibatani & Ito, 1942, Neozephyrus Sibatani & Ito, 1942, Quercusia Verity, 1943, Chrysozephyrus Shirôzu & Yamamoto, 1956, and Sibataniozephyrus Inomata, 1986 of Hypaurotis Scudder, 1876; Plesioarida Trujano & García, 2018 of Roeberella Strand, 1932; Papilio temenes Godart, 1819 (lectotype designated herein) of Heraclides aristodemus (Esper, 1794), Speyeria hydaspe conquista dos Passos & Grey, 1945 of Argynnis hesperis tetonia (dos Passos & Grey, 1945), and Erycides imbreus Plötz, 1879 of Phocides polybius polybius (Fabricius, 1793). The following are revised genus-species combinations: Pachythone lencates (Hewitson, 1875) Pachythone flocculus (Brévignon & Gallard, 1993), Pachythone floccus (Brévignon, 2013), Pachythone heberti (P. Jauffret & J. Jauffret, 2007), Pachythone marajoara (P. Jauffret & J. Jauffret, 2007) and Cissia cleophes (Godman & Salvin, 1889). The following species are transferred between subgenera: Anthocharis lanceolata Lucas, 1852 belongs to Anthocharis Boisduval, Rambur, [Duménil] & Graslin, [1833] instead of Paramidea Kuznetsov, 1929 and Danaus eresimus (Cramer, 1777) belongs to Danaus Kluk, 1780, and not to Anosia Hübner, 1816. The following taxa are distinct species rather than subspecies (of species shown in parenthesis): Heraclides ponceana (Schaus, 1911) (not Heraclides aristodemus (Esper, 1794)), Colias elis Strecker, 1885 (not Colias meadii W. H. Edwards, 1871), Argynnis irene Boisduval, 1869 and Argynnis nausicaa W. H. Edwards, 1874 (not Argynnis hesperis W. H. Edwards, 1864), Coenonympha california Westwood, [1851] (not Coenonympha tullia (Müller, 1764)), Dione incarnata N. Riley, 1926 (not Dione vanillae (Linnaeus, 1758)), Chlosyne coronado (M. Smith & Brock, 1988) (not Chlosyne fulvia (W. H. Edwards, 1879)), Chlosyne chinatiensis (Tinkham, 1944) (not Chlosyne theona (Ménétriés, 1855)), Phocides lilea (Reakirt, [1867]) (not Phocides polybius (Fabricius, 1793)), Cecropterus nevada (Scudder, 1872) and Cecropterus dobra (Evans, 1952) (not Cecropterus mexicana (Herrich-Schäffer, 1869)), Telegonus anausis Godman & Salvin, 1896, (not Telegonus anaphus (Cramer, 1777)), Epargyreus huachuca Dixon, 1955 (not Epargyreus clarus (Cramer, 1775)), Nisoniades bromias (Godman & Salvin, 1894) (not Nisoniades rubescens (Möschler, 1877)), Pholisora crestar J. Scott & Davenport, 2017 (not Pholisora catullus (Fabricius, 1793)), Carterocephalus mandan (W. H. Edwards, 1863) and Carterocephalus skada (W. H. Edwards, 1870) (not Carterocephalus palaemon (Pallas, 1771)), Amblyscirtes arizonae H. Freeman, 1993 (not Amblyscirtes elissa Godman, 1900), and Megathymus violae D. Stallings & Turner, 1956 (not Megathymus ursus Poling, 1902). Resulting from these changes, the following are revised species-subspecies combinations: Heraclides ponceana bjorndalae (Clench, 1979), Heraclides ponceana majasi L. Miller, 1987, Argynnis irene dodgei Gunder, 1931, Argynnis irene cottlei J. A. Comstock, 1925, Argynnis irene hanseni (J. Emmel, T. Emmel & Mattoon, 1998), Argynnis nausicaa elko (Austin, 1984), Argynnis nausicaa greyi (Moeck, 1950), Argynnis nausicaa viola (dos Passos & Grey, 1945), Argynnis nausicaa tetonia (dos Passos & Grey, 1945), Argynnis nausicaa chitone W. H. Edwards, 1879, Argynnis nausicaa schellbachi (Garth, 1949), Argynnis nausicaa electa W. H. Edwards, 1878, Argynnis nausicaa dorothea (Moeck, 1947), and Argynnis nausicaa capitanensis (R. Holland, 1988), Argynnis zerene atossa W. H. Edwards, 1890, Dione incarnata nigrior (Michener, 1942), Chlosyne coronado pariaensis (M. Smith & Brock, 1988), Cecropterus nevada aemilea (Skinner, 1893), Cecropterus nevada blanca (J. Scott, 1981), Telegonus anausis annetta (Evans, 1952), Telegonus anausis anoma (Evans, 1952), Telegonus anausis aniza (Evans, 1952), Epargyreus huachuca profugus Austin, 1998, Carterocephalus mandan mesapano (Scudder, 1868) and Carterocephalus skada magnus Mattoon & Tilden, 1998. American Coenonympha subspecies placed under C. tullia other than Coenonympha tullia kodiak W. H. Edwards, 1869, Coenonympha tullia mixturata Alpheraky, 1897 and Coenonympha tullia yukonensis W. Holland, 1900 belong to C. california. Heraclides ponceana latefasciatus Grishin, ssp. n. is described from Cuba. Argynnis coronis carolae dos Passos & Grey, 1942 is considered a subspecies-level taxon. Unless stated otherwise, all subgenera, species, subspecies and synonyms of mentioned genera and species are transferred together with their parent taxa, and others remain as previously classified.

Model-Based Inference of Punctuated Molecular Evolution

In standard models of molecular evolution, DNA sequences evolve through asynchronous substitutions according to Poisson processes with a constant rate (called the molecular clock) or a rate that can vary (relaxed clock). However, DNA sequences can also undergo episodes of fast divergence that will appear as synchronous substitutions affecting several sites simultaneously at the macroevolutionary timescale. Here, we develop a model, which we call the Relaxed Clock with Spikes model, combining basal, clock-like molecular substitutions with episodes of fast divergence called spikes arising at speciation events. Given a multiple sequence alignment and its time-calibrated species phylogeny, our model is able to detect speciation events (including hidden ones) cooccurring with spike events and to estimate the probability and amplitude of these spikes on the phylogeny. We identify the conditions under which spikes can be distinguished from the natural variance of the clock-like component of molecular substitutions and from variations of the clock. We apply the method to genes underlying snake venom proteins and identify several spikes at gene-specific locations in the phylogeny. This work should pave the way for analyses relying on whole genomes to inform on modes of species diversification.

Mimicry and mitonuclear discordance in nudibranchs: New insights from exon capture phylogenomics

Phylogenetic inference and species delimitation can be challenging in taxonomic groups that have recently radiated and where introgression produces conflicting gene trees, especially when species delimitation has traditionally relied on mitochondrial data and color pattern. Chromodoris, a genus of colorful and toxic nudibranch in the Indo-Pacific, has been shown to have extraordinary cryptic diversity and mimicry, and has recently radiated, ultimately complicating species delimitation. In these cases, additional genome-wide data can help improve phylogenetic resolution and provide important insights about evolutionary history. Here, we employ a transcriptome-based exon capture approach to resolve Chromodoris phylogeny with data from 2,925 exons and 1,630 genes, derived from 15 nudibranch transcriptomes. We show that some previously identified mimics instead show mitonuclear discordance, likely deriving from introgression or mitochondrial capture, but we confirm one "pure" mimic in Western Australia. Sister-species relationships and species-level entities were recovered with high support in both concatenated maximum likelihood (ML) and summary coalescent phylogenies, but the ML topologies were highly variable while the coalescent topologies were consistent across datasets. Our work also demonstrates the broad phylogenetic utility of 149 genes that were previously identified from eupulmonate gastropods. This study is one of the first to (a) demonstrate the efficacy of exon capture for recovering relationships among recently radiated invertebrate taxa, (b) employ genome-wide nuclear markers to test mimicry hypotheses in nudibranchs and (c) provide evidence for introgression and mitochondrial capture in nudibranchs.

Fasciola Species Introgression: Just a Fluke or Something More?

The threats posed by a range of viral and bacterial zoonotic diseases inevitably receive renewed attention in the wake of global pandemic events due to their overt and devastating impacts on human health and the economy. Parasitic zoonoses, however, many of which affect millions of people each day, are frequently ignored. In the case of fasciolosis, caused by infection with Fasciola hepatica or Fasciola gigantica, this oversight has allowed for the expansion of areas of parasite sympatry and thus increased the incidence of hybridization and possible introgression between the two species. Here we highlight how an increased demand for animal-derived protein, combined with a lack of appropriate tools for detection of these events, is changing the status quo of these zoonotic parasites.

Increased demand for animal-derived protein from Fasciola hepatica-endemic countries has led to a growing number of reports of hybridization between F. hepatica and Fasciola gigantica in Southeast Asia.

Hybridization and eventual introgression have been reported in a range of protozoan, helminth, and arthropod parasites and act as important drivers of evolutionary change and adaptation.

Introgression between Fasciola spp. remains unproven but has potentially serious human and animal health consequences as seen in other parasites.

New tools for the characterization of hybridization and introgression events between Fasciola spp. are needed.

Genetic and morphological variation of Oxymycterus (Rodentia: Sigmodontinae) in the Brazilian Atlantic Forest

We present a new assessment of the genetic and morphological variation within Oxymycterus quaestor Thomas, 1903, which currently includes the junior synonyms O. judex Thomas, 1909 and O. misionalis Sanborn, 1931. We integrate distinct lines of evidence, including variation of mitochondrial (Cytochrome b [Cytb]) and nuclear (intron 7 of beta fibrinogen gene [Fgb]) sequences, and the assessment of skull quantitative traits based on geometric morphometrics, throughout the Atlantic Forest of Southeastern-Southern Brazil, Argentina, and Paraguay. Phylogenetic relationships based on Cytb indicate that O. quaestor is structured in four well-supported clades (lineages A–D), one of them (lineage C) including topotypes of a previously associated nominal form (O. judex). However, these Cytb lineages exhibit lower levels of differentiation based on the Fgb locus, and are not recovered in the genealogies of this nuclear marker, representing a case of mitonuclear discordance. The Cytb lineages also broadly overlapped in the morphospace both in skull shape and size, which sustain the current wider concept of O. quaestor as one single young species (0.947 Myr) that is recently expanding, and ultimately branching out, in the Atlantic Forest.

Phylotranscriptomic evidence for pervasive ancient hybridization among Old World salamanders

Hybridization can leave genealogical signatures in an organism's genome, originating from the parental lineages and persisting over time. This potentially confounds phylogenetic inference methods that aim to represent evolution as a strictly bifurcating tree. We apply a phylotranscriptomic approach to study the evolutionary history of, and test for inter-lineage introgression in the Salamandridae, a Holarctic salamanders group of interest in studies of toxicity and aposematism, courtship behavior, and molecular evolution. Although the relationships between the 21 currently recognized salamandrid genera have been the subject of numerous molecular phylogenetic studies, some branches have remained controversial and sometimes affected by discordances between mitochondrial vs. nuclear trees. To resolve the phylogeny of this family, and understand the source of mito-nuclear discordance, we generated new transcriptomic (RNAseq) data for 20 salamandrids and used these along with published data, including 28 mitochondrial genomes, to obtain a comprehensive nuclear and mitochondrial perspective on salamandrid evolution. Our final phylotranscriptomic data set included 5455 gene alignments for 40 species representing 17 of the 21 salamandrid genera. Using concatenation and species-tree phylogenetic methods, we find (1) Salamandrina sister to the clade of the "True Salamanders" (consisting of Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), (2) Ichthyosaura sister to the Near Eastern genera Neurergus and Ommatotriton, (3) Triturus sister to Lissotriton, and (4) Cynops paraphyletic with respect to Paramesotriton and Pachytriton. Combining introgression tests and phylogenetic networks, we find evidence for introgression among taxa within the clades of "Modern Asian Newts" and "Modern European Newts". However, we could not unambiguously identify the number, position, and direction of introgressive events. Combining evidence from nuclear gene analysis with the observed mito-nuclear phylogenetic discordances, we hypothesize a scenario with hybridization and mitochondrial capture among ancestral lineages of (1) Lissotriton into Ichthyosaura and (2) Triturus into Calotriton, plus introgression of nuclear genes from Triturus into Lissotriton. Furthermore, both mitochondrial capture and nuclear introgression may have occurred among lineages assigned to Cynops. More comprehensive genomic data will, in the future, allow testing this against alternative scenarios involving hybridization with other, extinct lineages of newts.

Speciation in North American Junonia from a genomic perspective

Delineating species boundaries in phylogenetic groups undergoing recent radiation is a daunting challenge akin to discretizing continuity. Here, we propose a general approach exemplified by American butterflies from the genus Junonia Hübner notorious for the variety of similar phenotypes, ease of hybridization, and the lack of consensus about their classification. We obtain whole-genome shotgun sequences of about 200 specimens. We reason that discreteness emerges from continuity by means of a small number of key players, and search for the proteins that diverged markedly between sympatric populations of different species, while keeping low polymorphism within these species. Being 0.25% of the total number, these three dozen 'speciation' proteins indeed partition pairs of Junonia populations into two clusters with a prominent break in between, while all proteins taken together fail to reveal this discontinuity. Populations with larger divergence from each other, comparable to that between two sympatric species, form the first cluster and correspond to different species. The other cluster is characterized by smaller divergence, similar to that between allopatric populations of the same species and comprise conspecific pairs. Using this method, we conclude that J. genoveva (Cramer), J. litoralis Brévignon, J. evarete (Cramer), and J. divaricata C. & R. Felder are restricted to South America. We find that six species of Junonia are present in the United States, one of which is new: Junonia stemosa Grishin, sp.n. (i), found in south Texas and phenotypically closest to J. nigrosuffusa W. Barnes & McDunnough (ii) in its dark appearance. In the pale nudum of the antennal club, these two species resemble J. zonalis C. & R. Felder (iii) from Florida and the Caribbean Islands. The pair of sister species, J. grisea Austin & J. Emmel (iv) and J. coenia Hübner (v), represent the classic west/east U.S.A. split. The mangrove feeder (as caterpillar), dark nudum J. neildi Brévignon (vi) enters south Texas as a new subspecies Junonia neildi varia Grishin ssp.n. characterized by more extensive hybridization with and introgression from J. coenia, and, as a consequence, more variable wing patterns compared with the nominal J. n. neildi in Florida. Furthermore, a new mangrove-feeding species from the Pacific Coast of Mexico is described as Junonia pacoma Grishin sp.n. Finally, genomic analysis suggests that J. nigrosuffusa may be a hybrid species formed by the ancestors of J. grisea and J. stemosa sp.n.

The effect of hybridization on transposable element accumulation in an undomesticated fungal species

Transposable elements (TEs) are mobile genetic elements that can profoundly impact the evolution of genomes and species. A long-standing hypothesis suggests that hybridization could deregulate TEs and trigger their accumulation, although it received mixed support from studies mostly in plants and animals. Here, we tested this hypothesis in fungi using incipient species of the undomesticated yeast Saccharomyces paradoxus. Population genomic data revealed no signature of higher transposition in natural hybrids. As we could not rule out the elimination of past transposition increase signatures by natural selection, we performed a laboratory evolution experiment on a panel of artificial hybrids to measure TE accumulation in the near absence of selection. Changes in TE copy numbers were not predicted by the level of evolutionary divergence between the parents of a hybrid genotype. Rather, they were highly dependent on the individual hybrid genotypes, showing that strong genotype-specific deterministic factors govern TE accumulation in yeast hybrids.