The Tangled Evolutionary Legacies of Range Expansion and Hybridization

The ecological importance of hybridization

Hybridization as an evolutionary process has been studied in depth over the past few decades. Research has focused on its role in shaping reproductive barriers, its adaptive value, and its genomic consequences. In contrast, our knowledge of ecological dimensions of hybridization is still in its infancy, despite hybridization being an inherently ecological interaction. Using examples from various organisms, we show that hybridization can affect and be affected by non-reproductive interactions, including predation, competition, parasitism, mutualism, commensalism, and organism-environment interactions, with significant implications for community structure and ecosystem functioning. However, since these dimensions of hybridization have mostly been revealed from studies designed to decipher other evolutionary processes, we argue that much of the eco-evolutionary importance of hybridization is yet to be discovered.

Patterns of Genetic And Epigenetic Diversity Across A Range Expansion in The White-Footed Mouse (Peromyscus Leucopus)

Populations at the leading front of a range expansion must rapidly adapt to novel conditions. Increased epigenetic diversity has been hypothesized to facilitate adaptation and population persistence via non-genetic phenotypic variation, especially if there is reduced genetic diversity when populations expand (i.e., epigenetic diversity compensates for low genetic diversity). In this study, we use the spatial distribution of genetic and epigenetic diversity to test this hypothesis in populations of the white-footed mouse (Peromyscus leucopus) sampled across a purported recent range expansion gradient. We found mixed support for the epigenetic compensation hypothesis and a lack of support for expectations for expansion populations of mice at the range edge, which likely reflects a complex history of expansion in white-footed mice in the Upper Peninsula of Michigan. Specifically, epigenetic diversity was not increased in the population at the purported edge of the range expansion in comparison to the other expansion populations. However, input from an additional ancestral source populations may have increased genetic diversity at this range edge population, counteracting the expected genetic consequences of expansion, as well as reducing the benefit of increased epigenetic diversity at the range edge. Future work will expand the focal populations to include expansion areas with a single founding lineage to test for the robustness of a general trend that supports the hypothesized compensation of reduced genetic diversity by epigenetic variation observed in the expansion population that was founded from a single historical source.

Evolution of personality and locomotory performance traits during a late Pleistocene island colonization in a tree frog

Recent empirical and theoretical studies suggest that personality and locomotory performance traits linked to dispersal abilities are crucial components of the dispersal syndromes, and that they can evolve during range expansions and colonization processes. Island colonization is one of the best characterized processes in dispersal biogeography, and its implication in the evolution of phenotypic traits has been investigated over a wide range of temporal scales. However, the effect of island colonization on personality and performance traits of natural populations, and how these traits could drive island colonization, has been little explored. Noteworthy, no studies have addressed these processes in the context of late Pleistocene range expansions. Here, we investigated the contribution of island colonization triggered by postglacial range expansions to intraspecific variation in personality and locomotory performance traits. We compared boldness, exploration, jumping performance, and stickiness abilities in populations from 3 equidistant areas of the Tyrrhenian tree frog Hyla sarda, 2 from the main island (Corsica Island), and 1 from the recently colonized island of Elba. Individuals from Elba were significantly bolder than individuals from Corsica, as they emerged sooner from a shelter (P = 0.028), while individuals from Corsica showed markedly higher jumping and stickiness performance (both P < 0.001), resulting as more performing than those of Elba. We discuss these results in the context of the major microevolutionary processes at play during range expansion, including selection, spatial sorting, founder effects, and their possible interaction with local adaptation processes.

Different patterns of introgression in a three species hybrid zone among European cave salamanders

Hybrid zones occur where genetically distinct populations meet, mate and produce offspring with mixed ancestry. In Plethodontid salamanders, introgressive hybridization is a common phenomenon, where hybrids backcross with parental populations leading to the spread of new alleles into the parental genomes. Whereas many hybrid zones have been reported in American Plethodontid salamanders, only a single hybrid zone has been documented in European plethodontids so far, which is located at the Apuan Alps in the Italian Peninsula. Here, we describe a previously unreported hybrid zone in the Northern Apennines involving all the three Plethodontid salamander species inhabiting the Italian Peninsula. We found 21 new Speleomantes sites of occurrence, from a hitherto unexplored area located at the boundaries between three Speleomantes species ranges. Using mitochondrial (Cytb and ND2 genes) and nuclear markers (two diagnostic SNPs at the NCX1 gene), we revealed a three-way contact zone where all the three mainland species hybridize: S. strinatii, S. ambrosii and S. italicus. We observed a strong mitonuclear discordance, with mitochondrial markers showing a conspicuous geographic pattern, while diagnostic nuclear SNPs coexisted in both the same populations and individuals, providing evidence of hybridization in many possible combinations. The introgression is asymmetric, with S. italicus mitogenome usually associated with S. a. ambrosii and, to a lesser extent, to S. strinatii nuclear alleles. This finding confirms that Plethodontid are a group of choice to investigate hybridization mechanisms and suggests that behavioural, genetic and ecological components may concur in determining the direction and extent of introgression.

Surrounded by Kindred: Spermophilus major Hybridization with Other Spermophilus Species in Space and Time

Among the numerous described cases of hybridization in mammals, the most intriguing are (a) cases of introgressive hybridization deeply affecting the evolutionary history of species, and (b) models involving not a pair of species but a multi-species complex. Therefore, the hybridization history of the russet ground squirrel Spermophilus major, whose range has repeatedly changed due to climatic fluctuations and now borders the ranges of four related species, is of great interest. The main aims of this study were to determine the direction and intensity of gene introgression, the spatial depth of the infiltration of extraneous genes into the S. major range, and to refine the hypothesis of the hybridogenic replacement of mitochondrial genomes in the studied group. Using phylogenetic analysis of the variability of mitochondrial (CR, cytb) and nuclear (SmcY, BGN, PRKCI, c-myc, i6p53) markers, we determined the contribution of neighboring species to the S. major genome. We showed that 36% of S. major individuals had extraneous alleles. All peripheral species that were in contact with S. major contributed towards its genetic variability. We also proposed a hypothesis for the sequence and localization of serial hybridization events. Our assessment of the S. major genome implications of introgression highlights the importance of implementing conservation measures to protect this species.

De novo transcriptome assembly and annotation for gene discovery in Salamandra salamandra at the larval stage

Dispersal is a key process in ecology and evolutionary biology, as it shapes biodiversity patterns over space and time. Attitude to disperse is unevenly distributed among individuals within populations, and that individual personality can have pivotal roles in the shaping of this attitude. Here, we assembled and annotated the first de novo transcriptome of the head tissues of Salamandra salamandra from individuals, representative of distinct behavioral profiles. We obtained 1,153,432,918 reads, which were successfully assembled and annotated. The high-quality of the assembly was confirmed by three assembly validators. The alignment of contigs against the de novo transcriptome led to a mapping percentage higher than 94%. The homology annotation with DIAMOND led to 153,048 (blastx) and 95,942 (blastp) shared contigs, annotated on NR, Swiss-Prot and TrEMBL. The domain and site protein prediction led to 9850 GO-annotated contigs. This de novo transcriptome represents reliable reference for comparative gene expression studies between alternative behavioral types, for comparative gene expression studies within Salamandra, and for whole transcriptome and proteome studies in amphibians.

Molecular insights into the dynamics of species invasion by hybridisation in Tasmanian eucalypts

In plants where seed dispersal is limited compared with pollen dispersal, hybridisation may enhance gene exchange and species dispersal. We provide genetic evidence of hybridisation contributing to the expansion of the rare Eucalyptus risdonii into the range of the widespread Eucalyptus amygdalina. These closely related tree species are morphologically distinct, and observations suggest that natural hybrids occur along their distribution boundaries and as isolated trees or in small patches within the range of E. amygdalina. Hybrid phenotypes occur outside the range of normal dispersal for E. risdonii seed, yet in some hybrid patches small individuals resembling E. risdonii occur and are hypothesised to be a result of backcrossing. Using 3362 genome-wide SNPs assessed from 97 individuals of E. risdonii and E. amygdalina and 171 hybrid trees, we show that (i) isolated hybrids match the genotypes expected of F₁ /F₂ hybrids, (ii) there is a continuum in the genetic composition among the isolated hybrid patches from patches dominated by F₁ /F₂ -like genotypes to those dominated by E. risdonii-backcross genotypes, and (iii) the E. risdonii-like phenotypes in the isolated hybrid patches are most-closely related to proximal larger hybrids. These results suggest that the E. risdonii phenotype has been resurrected in isolated hybrid patches established from pollen dispersal, providing the first steps in its invasion of suitable habitat by long-distance pollen dispersal and complete introgressive displacement of E. amygdalina. Such expansion accords with the population demographics, common garden performance data, and climate modelling which favours E. risdonii and highlights a role of interspecific hybridisation in climate change adaptation and species expansion.

Secondary contact rather than coexistence-Erebia butterflies in the Alps

Secondary contact zones are ideal systems to study the processes that govern the evolution of reproductive barriers, especially at advanced stages of the speciation process. An increase in reproductive isolation resulting from selection against maladaptive hybrids is thought to contribute to reproductive barrier buildup in secondary contact zones. Although such processes have been invoked for many systems, it remains unclear to which extent they influence contact zone dynamics in nature. Here, we study a very narrow contact zone between the butterfly species Erebia cassioides and Erebia tyndarus in the Swiss Alps. We quantified phenotypic traits related to wing shape and reproduction as well as ecology to compare the degree of intra- and interspecific differentiation. Even though only very few first-generation hybrids occur, we find no strong indications for current reinforcing selection, suggesting that if reinforcement occurred in our system, it likely operated in the past. Additionally, we show that both species differ less in their ecological niche at the contact zone than elsewhere, which could explain why coexistence between these butterflies may currently not be possible.

Breaking the reproductive barrier of divergent species to explore the genomic landscape

Background: Climate change will have significant consequences for species. Species range shifts induce the emergence of new hybrid zones or the spatial displacement of pre-existing ones. These hybrid zones may become more porous as alleles are passed from one species to another. Currently, hybridization between highly divergent species living in sympatry seems extremely limited. Indeed, this phenomenon involves breaking two barriers. The first is the pre-mating barrier, related to the reproductive phenology of the two species. The second is the post-zygotic barrier, related to the genetic divergence between these species. Here, we were interested in identifying new hybridization patterns and potential implications, especially in the context of environmental modifications.

Methods: We sampled Telestes souffia and Parachondrostoma toxostoma wild specimens from different locations across France and genotyped them for SNP markers. We identified discriminant loci using F1-hybrid specimens and parental species and performed principal component analysis and Bayesian model-based clustering to analyze phylogenetic information. Furthermore, we assessed deviation in allele frequency from F1 to F2 and for Hardy–Weinberg equilibrium for F2 and assessed gene function associated with two F2 cohorts.

Results: We demonstrate that by breaking the ecological barrier, massive introgressive hybridization is possible between two endemic lineages of Cyprinidae belonging to two distinct genera. For both cohorts studied (=2 cm and >2 cm), a large majority of loci (>88%) presented no deviation in allele frequency and no departure from the Hardy–Weinberg equilibrium. For individuals beyond the 2 cm stage, two phenomena were observed. The first was an allelic imbalance in favor of P. toxostoma, for some genomic regions, with genes involved in developmental regulatory processes, cytoskeletal organization, and chromosome organization. The second was an excess of heterozygous loci coupled with an equilibrium of allelic frequencies for genes involved in immune response and kidney/liver development. Moreover, the 2 cm-sized specimens with high mortality yielded a particular genomic signature.

Conclusion: Our study displayed important results for understanding the early stages of hybridization between divergent lineages and predicting the emergence of future hybrid zones in the wild. Moreover, this hybridization generates a wide spectrum of hybrids that are a potential source of important evolutionary novelties.

Adaptive introgression facilitate adaptation to high latitudes in European aspen (Populus tremula L.)

Understanding local adaptation has become a key research area given the ongoing climate challenge and the concomitant requirement to conserve genetic resources. Perennial plants, such as forest trees, are good models to study local adaptation given their wide geographic distribution, largely outcrossing mating systems, and demographic histories. We evaluated signatures of local adaptation in European aspen (Populus tremula) across Europe by means of whole-genome resequencing of a collection of 411 individual trees. We dissected admixture patterns between aspen lineages and observed a strong genomic mosaicism in Scandinavian trees, evidencing different colonization trajectories into the peninsula from Russia, Central and Western Europe. As a consequence of the secondary contacts between populations after the last glacial maximum, we detected an adaptive introgression event in a genome region of ∼500 kb in chromosome 10, harboring a large-effect locus that has previously been shown to contribute to adaptation to the short growing seasons characteristic of Northern Scandinavia. Demographic simulations and ancestry inference suggest an Eastern origin—probably Russian—of the adaptive Nordic allele which nowadays is present in a homozygous state at the north of Scandinavia. The strength of introgression and positive selection signatures in this region is a unique feature in the genome. Furthermore, we detected signals of balancing selection, shared across regional populations, that highlight the importance of standing variation as a primary source of alleles that facilitate local adaptation. Our results, therefore, emphasize the importance of migration–selection balance underlying the genetic architecture of key adaptive quantitative traits.

Positive selection on two mitochondrial coding genes and adaptation signals in hares (genus Lepus) from China

BACKGROUND

Animal mitochondria play a central role in energy production in the cells through the oxidative phosphorylation (OXPHOS) pathway. Recent studies of selection on different mitochondrial OXPHOS genes have revealed the adaptive implications of amino acid changes in these subunits. In hares, climatic variation and/or introgression were suggested to be at the origin of such adaptation. Here we looked for evidence of positive selection in three mitochondrial OXPHOS genes, using tests of selection, protein structure modelling and effects of amino acid substitutions on the protein function and stability. We also used statistical models to test for climate and introgression effects on sites under positive selection.

RESULTS

Our results revealed seven sites under positive selection in ND4 and three sites in Cytb. However, no sites under positive selection were observed in the COX1 gene. All three subunits presented a high number of codons under negative selection. Sites under positive selection were mapped on the tridimensional structure of the predicted models for the respective mitochondrial subunit. Of the ten amino acid replacements inferred to have evolved under positive selection for both subunits, six were located in the transmembrane domain. On the other hand, three codons were identified as sites lining proton translocation channels. Furthermore, four codons were identified as destabilizing with a significant variation of Δ vibrational entropy energy between wild and mutant type. Moreover, our PROVEAN analysis suggested that among all positively selected sites two fixed amino acid replacements altered the protein functioning. Our statistical models indicated significant effects of climate on the presence of ND4 and Cytb protein variants, but no effect by trans-specific mitochondrial DNA introgression, which is not uncommon in a number of hare species.

CONCLUSIONS

Positive selection was observed in several codons in two OXPHOS genes. We found that substitutions in the positively selected codons have structural and functional impacts on the encoded proteins. Our results are concordantly suggesting that adaptations have strongly affected the evolution of mtDNA of hare species with potential effects on the protein function. Environmental/climatic changes appear to be a major trigger of this adaptation, whereas trans-specific introgressive hybridization seems to play no major role for the occurrence of protein variants.

Range and niche expansion through multiple interspecific hybridization: a genotyping by sequencing analysis of Cherleria (Caryophyllaceae)

BACKGROUND

Cherleria (Caryophyllaceae) is a circumboreal genus that also occurs in the high mountains of the northern hemisphere. In this study, we focus on a clade that diversified in the European High Mountains, which was identified using nuclear ribosomal (nrDNA) sequence data in a previous study. With the nrDNA data, all but one species was monophyletic, with little sequence variation within most species. Here, we use genotyping by sequencing (GBS) data to determine whether the nrDNA data showed the full picture of the evolution in the genomes of these species.

RESULTS

The overall relationships found with the GBS data were congruent with those from the nrDNA study. Most of the species were still monophyletic and many of the same subclades were recovered, including a clade of three narrow endemic species from Greece and a clade of largely calcifuge species. The GBS data provided additional resolution within the two species with the best sampling, C. langii and C. laricifolia, with structure that was congruent with geography. In addition, the GBS data showed significant hybridization between several species, including species whose ranges did not currently overlap.

CONCLUSIONS

The hybridization led us to hypothesize that lineages came in contact on the Balkan Peninsula after they diverged, even when those lineages are no longer present on the Balkan Peninsula. Hybridization may also have helped lineages expand their niches to colonize new substrates and different areas. Not only do genome-wide data provide increased phylogenetic resolution of difficult nodes, they also give evidence for a more complex evolutionary history than what can be depicted by a simple, branching phylogeny.

Land cover, individual's age and spatial sorting shape landscape resistance in the invasive frog Xenopus laevis

The description of functional connectivity is based on the quantification of landscape resistance, which represents species-specific movement costs across landscape features. Connectivity models use these costs to identify movement corridors at both individual and population levels and provide management recommendations for populations of conservation interest. Typically, resistance costs assigned to specific land cover types are assumed to be valid for all individuals of the population. Little attention has been paid to intraspecific variation in resistance costs due to age or dispersal syndrome, which may significantly affect model predictions. We quantified resistance costs in an expanding invasive population of the African clawed frog Xenopus laevis in Western France. In this principally aquatic amphibian, juveniles, sub-adults and adults disperse overland. The enhancement of dispersal traits via spatial sorting has been also observed at the range periphery of the population. Resistance costs, and thus connectivity, might vary as a function of life stage and position within the invaded range. We assessed multiple dimensions of functional connectivity. On various land cover types, we measured locomotion, as crossing speed, in different post-metamorphic age classes, and dehydration, sensitivity of locomotion to dehydration and substrate preference in juveniles. We also tested the effect of the position in the invaded range (core vs. periphery) on individual performances. In juveniles, general trends towards higher resistance costs on grass and lower resistance costs on bare soil and asphalt were observed, although not all experiments provided the same cost configurations. Resistance to locomotion varied between age classes, with adults and sub-adults facing lower costs than juveniles, particularly when crossing structurally complex land cover types such as grass and leaf litter. The position in the range had a minor effect on landscape resistance, and only in the dehydration experiment, where water loss in juveniles was lower at the range periphery. Depicting functional connectivity requires (a) assessing multiple dimensions of behavioural and physiological challenges faced by animals during movement; (b) considering factors, such as age and dispersal syndrome, that may affect movement at both individual and population levels. Ignoring this complexity might generate unreliable connectivity models and provide unsupported management recommendations for conservation.

Genetic Rescue and the Plight of Ponui Hybrids

Long-term sustainable and resilient populations is a key goal of conservation. How to best achieve this is controversial. There are, for instance, polarized views concerning the fitness and conservation value of hybrid populations founded through multi-origin translocations. A classic example concerns Apteryx (kiwi) in New Zealand. The A. mantelli of Ponui Island constitute a hybrid population where the birds are highly successful in their island habitat. A key dilemma for managers is understanding the reason for this success. Are the hybrid birds of Ponui Island of “no future conservation value” as recently asserted, or do they represent an outstanding example of genetic rescue and an important resource for future translocations? There has been a paradigm shift in scientific thinking concerning hybrids, but the ecological significance of admixed genomes remains difficult to assess. This limits what we can currently predict in conservation science. New understanding from genome science challenges the sufficiency of population genetic models to inform decision making and suggests instead that the contrasting outcomes of hybridization, “outbreeding depression” and “heterosis,” require understanding additional factors that modulate gene and protein expression and how these factors are influenced by the environment. We discuss these findings and the investigations that might help us to better understand the birds of Ponui, inform conservation management of kiwi and provide insight relevant for the future survival of Apteryx.

Biogeographic history moulds population differentiation in ageing of oxidative status in an amphibian

Regulation of oxidative status plays a substantial role in physiological ageing. However, we know little about age-related changes of oxidative status in wild animals, and even less about the role of population history in moulding ageing rates. We addressed these questions by means of a common garden experiment, using the Tyrrhenian tree frog Hyla sarda as the study species. This species underwent a range expansion from northern Sardinia (source) up to Corsica (newly founded) during the Late Pleistocene, and then the two populations became geographically isolated. We found that, at the beginning of the experiment, Sardinian and Corsican frogs had similar concentrations of all oxidative status markers analysed. One year later, Corsican frogs had higher oxidative stress and suffered higher mortality than Sardinian frogs. Our results suggest the intriguing scenario that population differentiation in rates of physiological ageing owing to oxidative stress might be an overlooked legacy of past biogeographic processes.

Secondary contact zones of closely-related Erebia butterflies overlap with narrow phenotypic and parasitic clines

Zones of secondary contact between closely related taxa are a common legacy of the Quaternary ice ages. Despite their abundance, the factors that keep species apart and prevent hybridization are often unknown. Here, we study a very narrow contact zone between three closely related butterfly species of the Erebia tyndarus species complex. Using genomic data, we first determined whether gene flow occurs and then assessed whether it might be hampered by differences in chromosome number between some species. We found interspecific gene flow between sibling species that differ in karyotype by one chromosome. Conversely, only F1 hybrids occurred between two species that have the same karyotype, forming a steep genomic cline. In a second step, we fitted clines to phenotypic, ecological and parasitic data to identify the factors associated with the genetic cline. We found clines for phenotypic data and the prevalence of the endosymbiont parasite Wolbachia to overlap with the genetic cline, suggesting that they might be drivers for separating the two species. Overall, our results highlight that some gene flow is possible between closely related species despite different chromosome numbers, but that other barriers restrict such gene flow.

You're Just My Type: Mate Choice and Behavioral Types

Consistent individual differences in behavior [i.e., behavioral types (BTs)], are common across the animal kingdom. Consistency can make behavior an adaptive trait for mate choice decisions. Here, we present a conceptual framework to explain how and why females might evaluate a male's BT before mating. Because BTs are consistent across time or context, a male's BT can be a reliable indicator of his potential to provide direct benefits. Heritable BTs can enable informed mate choice via indirect benefits. Many key issues regarding patterns of mate choice, including sensory biases, context dependence, and assortative mating apply to BT-dependent mate choice. Understanding the relationship between BTs and mate choice may offer insights into patterns of variation and consistency common in behavioral traits.

Getting Back to Nature: Feralization in Animals and Plants

Formerly domesticated organisms and artificially selected genes often escape controlled cultivation, but their subsequent evolution is not well studied. In this review, we examine plant and animal feralization through an evolutionary lens, including how natural selection, artificial selection, and gene flow shape feral genomes, traits, and fitness. Available evidence shows that feralization is not a mere reversal of domestication. Instead, it is shaped by the varied and complex histories of feral populations, and by novel selection pressures. To stimulate further insight we outline several future directions. These include testing how 'domestication genes' act in wild settings, studying the brains and behaviors of feral animals, and comparative analyses of feral populations and taxa. This work offers feasible and exciting research opportunities with both theoretical and practical applications.

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.

Assessing Behavioral Associations in a Hybrid Zone through Social Network Analysis: Complex Assortative Behaviors Structure Associations in a Hybrid Quail Population

Behavior can strongly influence rates and patterns of hybridization between animal populations and species. Yet few studies have examined reproductive behaviors in natural hybrid zones within the fine-scale social context in which they naturally occur. We use radio-frequency identification tags with social network analyses to test whether phenotypic similarity in plumage and mass correlate with social behavior throughout a breeding season in a California and Gambel's quail hybrid zone. We use a novel approach to partition phenotypic variation in a way that does not confound differences between sexes and species, and we illustrate the complex ways that phenotype and behavior structure the social environment, mating opportunities, and male-male associations. Associations within the admixed population were random with respect to species-specific plumage but showed strong patterns of assortment based on sexually dimorphic plumage, monomorphic plumage, and mass. Weak behavioral reproductive isolation in this admixed population may be the result of complex patterns of phenotypic assortment based on multiple traits rather than a lack of phenotypic discrimination. More generally, our results support the utility of social network analyses for analyzing behavioral factors affecting genetic exchange between populations and species.

Asymmetric hybridization in Cordulegaster (Odonata: Cordulegastridae): Secondary postglacial contact and the possible role of mechanical constraints

Two Cordulegaster dragonflies present in Italy, the Palaearctic and northern distributed Cordulegaster boltonii and the endemic to the south of the peninsula Cordulegaster trinacriae, meet in central Italy and give rise to individuals of intermediate morphology. By means of mitochondrial and nuclear markers and of Geometric Morphometrics applied to sexual appendages, we defined i) the geographical boundaries between the two species in Italy and ii) we determined the presence, the extent, and the genetic characteristics of the hybridization. Genetic data evidenced asymmetric hybridization with the males of C. trinacriae able to mate both interspecifically and intraspecifically. The results contrast with expectations under neutral gene introgression and sexual selection. This data, along with the morphological evidence of significant differences in size and shape of sexual appendages between the males of the two species, seem indicative of the role of mechanical constraints in intraspecific matings. The origin of the two species is dated about to 1.32 Mya and the hybridization resulted related to range expansion of the two species after Last Glacial Maximum and this led to the secondary contact between the two taxa in central Italy. At last, our results indicate that the range of C. trinacriae, a threatened and protected species, has been moving northward probably driven by climate changes. As a result, the latter species is currently intruding into the range of C. boltonii. The hybrid area is quite extended and the hybrids seem well adapted to the environment. From a conservation point of view, even if C. trinacriae has a strong genetic identity, the discovery of hybridization between the two species should be considered in a future species management.

Hybridization and extensive mitochondrial introgression among fire salamanders in peninsular Italy

Discordance between mitochondrial and nuclear patterns of population genetic structure is providing key insights into the eco-evolutionary dynamics between and within species, and their assessment is highly relevant to biodiversity monitoring practices based on DNA barcoding approaches. Here, we investigate the population genetic structure of the fire salamander Salamandra salamandra in peninsular Italy. Both mitochondrial and nuclear markers clearly identified two main population groups. However, nuclear and mitochondrial zones of geographic transition between groups were located 600 km from one another. Recent population declines in central Italy partially erased the genetic imprints of past hybridization dynamics. However, the overall pattern of genetic variation, together with morphological and fossil data, suggest that a rampant mitochondrial introgression triggered the observed mitonuclear discordance, following a post-glacial secondary contact between lineages. Our results clearly show the major role played by reticulate evolution in shaping the structure of Salamandra salamandra populations and, together with similar findings in other regions of the species' range, contribute to identify the fire salamander as a particularly intriguing case to investigate the complexity of mechanisms triggering patterns of mitonuclear discordance in animals.

Experimental evidence that thermal selection shapes mitochondrial genome evolution

Mitochondria are essential organelles, found within eukaryotic cells, which contain their own DNA. Mitochondrial DNA (mtDNA) has traditionally been used in population genetic and biogeographic studies as a maternally-inherited and evolutionary-neutral genetic marker. However, it is now clear that polymorphisms within the mtDNA sequence are routinely non-neutral, and furthermore several studies have suggested that such mtDNA polymorphisms are also sensitive to thermal selection. These observations led to the formulation of the “mitochondrial climatic adaptation” hypothesis, for which all published evidence to date is correlational. Here, we use laboratory-based experimental evolution in the fruit fly, Drosophila melanogaster, to test whether thermal selection can shift population frequencies of two mtDNA haplogroups whose natural frequencies exhibit clinal associations with latitude along the Australian east-coast. We present experimental evidence that the thermal regime in which the laboratory populations were maintained drove changes in haplogroup frequencies across generations. Our results strengthen the emerging view that intra-specific mtDNA variants are sensitive to selection, and suggest spatial distributions of mtDNA variants in natural populations of metazoans might reflect adaptation to climatic environments rather than within-population coalescence and diffusion of selectively-neutral haplotypes across populations.

Molecular and ecological signatures of an expanding hybrid zone

Many species are currently changing their distributions and subsequently form sympatric zones with hybridization between formerly allopatric species as one possible consequence. The damselfly Ischnura elegans has recently expanded south into the range of its ecologically and morphologically similar sister species Ischnura graellsii. Molecular work shows ongoing introgression between these species, but the extent to which this species mixing is modulated by ecological niche use is not known. Here, we (1) conduct a detailed population genetic analysis based on molecular markers and (2) model the ecological niche use of both species in allopatric and sympatric regions. Population genetic analyses showed chronic introgression between I. elegans and I. graellsii across a wide part of Spain, and admixture analysis corroborated this, showing that the majority of I. elegans from the sympatric zone could not be assigned to either the I. elegans or I. graellsii species cluster. Niche modeling demonstrated that I. elegans has modified its environmental niche following hybridization and genetic introgression with I. graellsii, making niche space of introgressed I. elegans populations more similar to I. graellsii. Taken together, this corroborates the view that adaptive introgression has moved genes from I. graellsii into I. elegans and that this process is enabling Spanish I. elegans to occupy a novel niche, further facilitating its expansion. Our results add to the growing evidence that hybridization can play an important and creative role in the adaptive evolution of animals.

Interspecific competition, hybridization, and reproductive isolation in secondary contact: missing perspectives on males and females

Research on sexual selection and hybridization has focused on female mate choice and male-male competition. While the evolutionary outcomes of interspecific female preference have been well explored, we are now gaining a better understanding of the processes by which male-male competition between species in secondary contact promotes reproductive isolation versus hybridization. What is relatively unexplored is the interaction between female choice and male competition, as they can oppose one another or align with similar outcomes for reproductive isolation. The role of female-female competition in hybridization is also not well understood, but could operate similarly to male-male competition in polyandrous and other systems where costs to heterospecific mating are low for females. Reproductive competition between either sex of sympatric species can cause the divergence and/or convergence of sexual signals and recognition, which in turn influences the likelihood for interspecific mating. Future work on species interactions in secondary contact should test the relative influences of both mate choice and competition for mates on hybridization outcomes, and should not ignore the possibilities that females can compete over mating resources, and males can exercise mate choice.

Cryptic diversity in an Atlantic Forest malaria vector from the mountains of South-East Brazil

BACKGROUND

Anopheles (Kerteszia) cruzii is the primary vector of human and simian malarias in Brazilian regions covered by the Atlantic Rainforest. Previous studies found that An. cruzii presents high levels of behavioural, chromosomal and molecular polymorphisms, which led to the hypothesis that it may be a complex of cryptic species. Here, An. cruzii specimens were collected in five sites in South-East Brazil located at different altitudes on the inner and coastal slopes of two mountain ranges covered by Atlantic Rainforest, known as Serra do Mar and Serra da Mantiqueria. Partial sequences for two genes (Clock and cpr) were generated and compared with previously published sequences from Florianópolis (southern Brazil). Genetic diversity was analysed with estimates of population structure (F _ST ) and haplotype phylogenetic trees in order to understand how many species of the complex may occur in this biome and how populations across the species distribution are related.

RESULTS

The sequences from specimens collected at sites located on the lower coastal slopes of Serra do Mar (Guapimirim, Tinguá and Sana) clustered together in the phylogenetic analysis, while the major haplotypes from sites located on higher altitude and at the continental side of the same mountains (Bocaina) clustered with those from Serra da Mantiqueira (Itatiaia), an inner mountain range. These two An. cruzii lineages showed statistically significant genetic differentiation and fixed characters, and have high F _ST values typical of between species comparisons. Finally, in Bocaina, where the two lineages occur in sympatry, we found deviations from Hardy-Weinberg equilibrium due to a deficit of heterozygotes, indicating partial reproductive isolation. These results strongly suggest that at least two distinct lineages of An. cruzii (provisorily named "Group 1" and "Group 2") occur in the mountains of South-East Brazil.

CONCLUSIONS

At least two genetically distinct An. cruzii lineages occur in the Atlantic Forest covered mountains of South-East Brazil. The co-occurrence of distinct lineages of An. cruzii (possibly incipient species) in those mountains is an interesting biological phenomenon and may have important implications for malaria prevalence, Plasmodium transmission dynamics and control.

Molecular Epidemiology of Anisakis and Anisakiasis: An Ecological and Evolutionary Road Map

This review addresses the biodiversity, biology, distribution, ecology, epidemiology, and consumer health significance of the so far known species of Anisakis, both in their natural hosts and in human accidental host populations, worldwide. These key aspects of the Anisakis species' biology are highlighted, since we consider them as main driving forces behind which most of the research in this field has been carried out over the past decade. From a public health perspective, the human disease caused by Anisakis species (anisakiasis) appears to be considerably underreported and underestimated in many countries or regions around the globe. Indeed, when considering the importance of marine fish species as part of the everyday diet in many coastal communities around the globe, there still exist significant knowledge gaps as to local epidemiological and ecological drivers of the transmission of Anisakis spp. to humans. We further identify some key knowledge gaps related to Anisakis species epidemiology in both natural and accidental hosts, to be filled in light of new 'omic' technologies yet to be fully developed. Moreover, we suggest that future Anisakis research takes a 'holistic' approach by integrating genetic, ecological, immunobiological, and environmental factors, thus allowing proper assessment of the epidemiology of Anisakis spp. in their natural hosts, in human populations, and in the marine ecosystem, in both space and time.

Admixture on the northern front: population genomics of range expansion in the white-footed mouse (Peromyscus leucopus) and secondary contact with the deer mouse (Peromyscus maniculatus)

Range expansion has genetic consequences expected to result in differentiated wave-front populations with low genetic variation and potentially introgression from a local species. The northern expansion of Peromyscus leucopus in southern Quebec provides an opportunity to test these predictions using population genomic tools. Our results show evidence of recent and post-glacial expansion. Genome-wide variation in P. leucopus indicates two post-glacial lineages are separated by the St. Lawrence River, with a more recent divergence of populations isolated by the Richelieu River. In two of three transects we documented northern populations with low diversity in at least one genetic measure, although most relationships were not significant. Consistent with bottlenecks and allele surfing during northward expansion, we document a northern-most population with low nucleotide diversity, divergent allele frequencies and the most private alleles, and observed heterozygosity indicates outcrossing. Ancestry proportions revealed putative hybrids of P. leucopus and P. maniculatus. A formal test for gene flow confirmed secondary contact, showing that a reticulate population phylogeny between P. maniculatus and P. leucopus was a better fit to the data than a bifurcating model without gene flow. Thus, we provide the first genomic evidence of gene flow between this pair of species in natural populations. Understanding the evolutionary consequences of secondary contact is an important conservation concern as climate-induced range expansions are expected to result in new hybrid zones between closely related species.

Climate change promotes hybridisation between deeply divergent species

Rare hybridisations between deeply divergent animal species have been reported for decades in a wide range of taxa, but have often remained unexplained, mainly considered chance events and reported as anecdotal. Here, we combine field observations with long-term data concerning natural hybridisations, climate, land-use, and field-validated species distribution models for two deeply divergent and naturally sympatric toad species in Europe (Bufo bufo and Bufotes viridis species groups). We show that climate warming and seasonal extreme temperatures are conspiring to set the scene for these maladaptive hybridisations, by differentially affecting life-history traits of both species. Our results identify and provide evidence of an ultimate cause for such events, and reveal that the potential influence of climate change on interspecific hybridisations goes far beyond closely related species. Furthermore, climate projections suggest that the chances for these events will steadily increase in the near future.

Digging up the roots of an insular hotspot of genetic diversity: decoupled mito-nuclear histories in the evolution of the Corsican-Sardinian endemic lizard Podarcis tiliguerta

BACKGROUND

Mediterranean islands host a disproportionately high level of biodiversity and endemisms. Growing phylogeographic evidence on island endemics has unveiled unexpectedly complex patterns of intra-island diversification, which originated at diverse spatial and temporal scales. We investigated multilocus genetic variation of the Corsican-Sardinian endemic lizard Podarcis tiliguerta with the aim of shedding more light on the evolutionary processes underlying the origin of Mediterranean island biodiversity.

RESULTS

We analysed DNA sequences of mitochondrial (12S and nd4) and nuclear (acm4 and mc1r) gene fragments in 174 individuals of P. tiliguerta from 81 localities across the full range of the species in a geographic and genealogical framework. We found surprisingly high genetic diversity both at mitochondrial and nuclear loci. Seventeen reciprocally monophyletic allopatric mitochondrial haplogroups were sharply divided into four main mitochondrial lineages (two in Corsica and two in Sardinia) of Miocene origin. In contrast, shallow divergence and shared diversity within and between islands was observed at the nuclear loci. We evaluated alternative biogeographic and evolutionary scenarios to explain such profound discordance in spatial and phylogenetic patterning between mitochondrial and nuclear genomes. While neutral models provided unparsimonious explanations for the observed pattern, the hypothesis of environmental selection driving mitochondrial divergence in the presence of nuclear gene flow is favoured.

CONCLUSIONS

Our study on the genetic variation of P. tiliguerta reveals surprising levels of diversity underlining a complex phylogeographic pattern with a striking example of mito-nuclear discordance. These findings have profound implications, not only for the taxonomy and conservation of P. tiliguerta. Growing evidence on deep mitochondrial breaks in absence of geographic barriers and of climatic factors associated to genetic variation of Corsican-Sardinian endemics warrants additional investigation on the potential role of environmental selection driving the evolution of diversity hotspots within Mediterranean islands.

Range expansion underlies historical introgressive hybridization in the Iberian hare

Introgressive hybridization is an important and widespread evolutionary process, but the relative roles of neutral demography and natural selection in promoting massive introgression are difficult to assess and an important matter of debate. Hares from the Iberian Peninsula provide an appropriate system to study this question. In its northern range, the Iberian hare, Lepus granatensis, shows a northwards gradient of increasing mitochondrial DNA (mtDNA) introgression from the arctic/boreal L. timidus, which it presumably replaced after the last glacial maximum. Here, we asked whether a south-north expansion wave of L. granatensis into L. timidus territory could underlie mtDNA introgression, and whether nuclear genes interacting with mitochondria (“mitonuc” genes) were affected. We extended previous RNA-sequencing and produced a comprehensive annotated transcriptome assembly for L. granatensis. We then genotyped 100 discovered nuclear SNPs in 317 specimens spanning the species range. The distribution of allele frequencies across populations suggests a northwards range expansion, particularly in the region of mtDNA introgression. We found no correlation between variants at 39 mitonuc genes and mtDNA introgression frequency. Whether the nuclear and mitochondrial genomes coevolved will need a thorough investigation of the hundreds of mitonuc genes, but range expansion and species replacement likely promoted massive mtDNA introgression.