Maternal genomic variability of the wild boar (Sus scrofa) reveals the uniqueness of East‐Caucasian and Central Italian populations

The phylogeography of the European wild boar was mainly determined by postglacial recolonization patterns from Mediterranean refugia after the last ice age. Here we present the first analysis of SNP polymorphism within the complete mtDNA genome of West Russian (n = 8), European (n = 64), and North African (n = 5) wild boar. Our analyses provided evidence of unique lineages in the East‐Caucasian (Dagestan) region and in Central Italy. A phylogenetic analysis revealed that these lineages are basal to the other European mtDNA sequences. We also show close connection between the Western Siberian and Eastern European populations. Also, the North African samples were clustered with the Iberian population. Phylogenetic trees and migration modeling revealed a high proximity of Dagestan sequences to those of Central Italy and suggested possible gene flow between Western Asia and Southern Europe which was not directly related to Northern and Central European lineages. Our results support the presence of old maternal lineages in two Southern glacial refugia (i.e., Caucasus and the Italian peninsula), as a legacy of an ancient wave of colonization of Southern Europe from an Eastern origin.     

Signatures of de‐domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution

Autochthonous pig breeds are usually reared in extensive or semi‐extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo‐Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south‐eastern European countries (Kr?kopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro‐geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild‐type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of ‘de‐domestication’ process, and wild resources are challenged by a ‘domestication’ drift. Both need to be further investigated and managed.     

Introgression between invasive and native blue mussels (genus Mytilus) in the central California hybrid zone

The ecological and genetic factors determining the extent of introgression between species in secondary contact zones remain poorly understood. Here, we investigate the relative importance of isolating barriers and the demographic expansion of invasive Mytilus galloprovincialis on the magnitude and the direction of introgression with the native Mytilus trossulus in a hybrid zone in central California. We use double‐digest restriction‐site‐associated DNA sequencing (ddRADseq) to genotype 1337 randomly selected single nucleotide polymorphisms and accurately distinguish early and advanced generation hybrids for the first time in the central California Mytilus spp. hybrid zone. Weak levels of introgression were observed in both directions but were slightly more prevalent from the native M. trossulus into the invasive M. galloprovincialis. Few early and advanced backcrossed individuals were observed across the hybrid zone confirming the presence of strong barriers to interbreeding. Heterogeneous patterns of admixture across the zone of contact were consistent with the colonization history of M. galloprovincialis with more extensive introgression in northern localities furthest away from the putative site of introduction in southern California. These observations reinforce the importance of dynamic spatial and demographic expansions in determining patterns of introgression between close congeners, even in those with high dispersal potential and well‐developed reproductive barriers. Our results suggest that the threat posed by invasive M. galloprovincialis is more ecological than genetic as it has displaced, and continues to displace the native M. trossulus from much of central and southern California.     

Testing models of speciation from genome sequences: divergence and asymmetric admixture in Island South‐East Asian Sus species during the Plio‐Pleistocene climatic fluctuations

In many temperate regions, ice ages promoted range contractions into refugia resulting in divergence (and potentially speciation), while warmer periods led to range expansions and hybridization. However, the impact these climatic oscillations had in many parts of the tropics remains elusive. Here, we investigate this issue using genome sequences of three pig (Sus) species, two of which are found on islands of the Sunda‐shelf shallow seas in Island South‐East Asia (ISEA). A previous study revealed signatures of interspecific admixture between these Sus species (Genome biology, 14 , 2013, R107). However, the timing, directionality and extent of this admixture remain unknown. Here, we use a likelihood‐based model comparison to more finely resolve this admixture history and test whether it was mediated by humans or occurred naturally. Our analyses suggest that interspecific admixture between Sunda‐shelf species was most likely asymmetric and occurred long before the arrival of humans in the region. More precisely, we show that these species diverged during the late Pliocene but around 23% of their genomes have been affected by admixture during the later Pleistocene climatic transition. In addition, we show that our method provides a significant improvement over D‐statistics which are uninformative about the direction of admixture.     

Genetic diversity and structuring across the range of a widely distributed ladybird: focus on rear‐edge populations phenotypically divergent

Population genetics and phenotypic structures are often predicted to vary along the geographic range of a species. This phenomenon would be accentuated for species with large range areas, with discontinuities and marginal populations. We herein compare the genetic patterns of central populations of Coccinella septempunctata L. with those of two phenotypically differentiated populations considered as rear‐edge populations and subspecies based on phenotype (Algeria and Japan). According to the central‐marginal model and expected characteristics of rear‐edge populations, we hypothesize that these rear‐edge populations have (1) a reduced genetic diversity, resulting from their relative isolation over long periods of time, (2) a higher population genetic differentiation, explained by low contemporary gene flow levels, and (3) a relationship between genetic diversity characteristics and phenotypes, due to historical isolation and/or local adaptation. Based on genotyping of 28 populations for 18 microsatellite markers, several levels of regional genetic diversity and differentiation are observed between and within populations, according to their localization: low within‐population genetic diversity and higher genetic differentiation of rear‐edge populations. The genetic structuring clearly dissociates the Algerian and Eastern Asia populations from the others. Geographical patterns of genetic diversity and differentiation support the hypothesis of the central‐marginal model. The pattern observed is in agreement with the phenotypic structure across species range. A clear genetic break between populations of Algeria, the Eastern Asia, and the remaining populations is a dominant feature of the data. Differential local adaptations, absence of gene flow between marginal and central populations, and/or incapacity to mate after colonization, have contributed to their distinct genotypic and phenotypic characteristics.     

Evolutionary history of Purple cone spruce (Picea purpurea) in the Qinghai–Tibet Plateau: homoploid hybrid origin and Pleistocene expansion

Hybridization and introgression can play an important role in speciation. Here, we examine their roles in the origin and evolution of Picea purpurea, a diploid spruce species occurring on the Qinghai–Tibet Plateau (QTP). Phylogenetic relationships and ecological differences between this species and its relatives, P. schrenkiana, P. likiangensis and P. wilsonii, are unclear. To clarify them, we surveyed sequence variation within and between them for 11 nuclear loci, three chloroplast (cp) and two mitochondrial (mt) DNA fragments, and examined their ecological requirements using ecological niche modelling. Initial analyses based on 11 nuclear loci rejected a close relationship between P. schrenkiana and P. purpurea. BP&P tests and ecological niche modelling indicated substantial divergence between the remaining three species and supported the species status of P. purpurea, which contained many private alleles as expected for a well‐established species. Sequence variation for cpDNA and mtDNA suggested a close relationship between P. purpurea and P. wilsonii, while variation at the nuclear se1364 gene suggested P. purpurea was more closely related to P. likiangensis. Analyses of genetic divergence, Bayesian clustering and model comparison using approximate Bayesian computation (ABC) of nuclear (nr) DNA variation all supported the hypothesis that P. purpurea originated by homoploid hybrid speciation from P. wilsonii and P. likiangensis. The ABC analysis dated the origin of P. purpurea at the Pleistocene, and the estimated hybrid parameter indicated that 69% of its nuclear composition was contributed by P. likiangensis and 31% by P. wilsonii. Our results further suggested that during or immediately following its formation, P. purpurea was subject to organelle DNA introgression from P. wilsonii such that it came to possess both mtDNA and cpDNA of P. wilsonii. The estimated parameters indicated that following its origin, P. purpurea underwent an expansion during/after the largest Pleistocene glaciation recorded for the QTP.     

Environmental (in)dependence of a hybrid zone: Insights from molecular markers and ecological niche modeling in a hybrid zone of Origanum (Lamiaceae) on the island of Crete

The role of environment and the relative significance of endogenous versus exogenous selection in shaping hybrid zones have been crucial issues in the studies of hybridization. Recent advances in ecological niche modeling (ENM) offer new methodological tools, especially in combination with the genotyping of individuals in the hybrid zone. Here, we study the hybrid zone between the widely known spices Origanum onites and Origanum vulgare ssp. hirtum in Crete. We analyze the genetic structure of both parental taxa and their hybrid Origanum × intercendens using AFLP markers on 15 sympatric and 12 allopatric populations and employ ecological niche modeling and niche similarity tests to study their niche patterns. We complement these analyses with seed viability measurements. Our study revealed that the hybridizing taxa O. onites and O. vulgare ssp. hirtum and the resulting genotypic classes showed geographical and environmental niche similarities based on the predictions of ENMs and the subsequent similarity tests. The occurrence of the hybrid zone is not directly dependent on environmental factors which favor the fitness of the hybrid compared to the parental taxa, but rather on aspects such as historical factors and management practices, which may contribute to the localization and maintenance of the contact zone between parental species. Our results suggest that if a minimum required niche differentiation between genotypic classes is not achieved, environmental dependence might not have a prominent role on the outcome of the hybridization.     

Optimizing the trade‐off between spatial and genetic sampling efforts in patchy populations: towards a better assessment of functional connectivity using an individual‐based sampling scheme

Genetic data are increasingly used in landscape ecology for the indirect assessment of functional connectivity, that is, the permeability of landscape to movements of organisms. Among available tools, matrix correlation analyses (e.g. Mantel tests or mixed models) are commonly used to test for the relationship between pairwise genetic distances and movement costs incurred by dispersing individuals. When organisms are spatially clustered, a population‐based sampling scheme (PSS) is usually performed, so that a large number of genotypes can be used to compute pairwise genetic distances on the basis of allelic frequencies. Because of financial constraints, this kind of sampling scheme implies a drastic reduction in the number of sampled aggregates, thereby reducing sampling coverage at the landscape level. We used matrix correlation analyses on simulated and empirical genetic data sets to investigate the efficiency of an individual‐based sampling scheme (ISS) in detecting isolation‐by‐distance and isolation‐by‐barrier patterns. Provided that pseudo‐replication issues are taken into account (e.g. through restricted permutations in Mantel tests), we showed that the use of interindividual measures of genotypic dissimilarity may efficiently replace interpopulation measures of genetic differentiation: the sampling of only three or four individuals per aggregate may be sufficient to efficiently detect specific genetic patterns in most situations. The ISS proved to be a promising methodological alternative to the more conventional PSS, offering much flexibility in the spatial design of sampling schemes and ensuring an optimal representativeness of landscape heterogeneity in data, with few aggregates left unsampled. Each strategy offering specific advantages, a combined use of both sampling schemes is discussed.     

Phylogeny and biogeography of the American live oaks (Quercus subsection Virentes): a genomic and population genetics approach

The nature and timing of evolution of niche differentiation among closely related species remains an important question in ecology and evolution. The American live oak clade, Virentes, which spans the unglaciated temperate and tropical regions of North America and Mesoamerica, provides an instructive system in which to examine speciation and niche evolution. We generated a fossil‐calibrated phylogeny of Virentes using RADseq data to estimate divergence times and used nuclear microsatellites, chloroplast sequences and an intron region of nitrate reductase (NIA‐i3) to examine genetic diversity within species, rates of gene flow among species and ancestral population size of disjunct sister species. Transitions in functional and morphological traits associated with ecological and climatic niche axes were examined across the phylogeny. We found the Virentes to be monophyletic with three subclades, including a southwest clade, a southeastern US clade and a Central American/Cuban clade. Despite high leaf morphological variation within species and transpecific chloroplast haplotypes, RADseq and nuclear SSR data showed genetic coherence of species. We estimated a crown date for Virentes of 11 Ma and implicated the formation of the Sea of Cortés in a speciation event ~5 Ma. Tree height at maturity, associated with fire tolerance, differs among the sympatric species, while freezing tolerance appears to have diverged repeatedly across the tropical–temperate divide. Sympatric species thus show evidence of ecological niche differentiation but share climatic niches, while allopatric and parapatric species conserve ecological niches, but diverge in climatic niches. The mode of speciation and/or degree of co‐occurrence may thus influence which niche axis plants diverge along.     

Pollen transfer in fragmented plant populations: insight from the pollen loads of pollinators and stigmas in a mass‐flowering species

Pollinator and/or mate scarcity affects pollen transfer, with important ecological and evolutionary consequences for plant reproduction. However, the way in which the pollen loads transported by pollinators and deposited on stigmas are affected by pollination context has been little studied. We investigated the impacts of plant mate and visiting insect availabilities on pollen transport and receipt in a mass‐flowering and facultative autogamous shrub (Rhododendron ferrugineum). First, we recorded insect visits to R. ferrugineum in plant patches of diverse densities and sizes. Second, we analyzed the pollen loads transported by R. ferrugineum pollinators and deposited on stigmas of emasculated and intact flowers, in the same patches. Overall, pollinators (bumblebees) transported much larger pollen loads than the ones found on stigmas, and the pollen deposited on stigmas included a high proportion of conspecific pollen. However, comparing pollen loads of emasculated and intact flowers indicated that pollinators contributed only half the conspecific pollen present on the stigma. At low plant density, we found the highest visitation rate and the lowest proportion of conspecific pollen transported and deposited by pollinators. By contrast, at higher plant density and lower visitation rate, pollinators deposited larger proportion of conspecific pollen, although still far from sufficient to ensure that all the ovules were fertilized. Finally, self‐pollen completely buffered the detrimental effects on pollination of patch fragmentation and pollinator failure. Our results indicate that pollen loads from pollinators and emasculated flowers should be quantified for an accurate understanding of the relative impacts of pollinator and mate limitation on pollen transfer in facultative autogamous species.     

Past and future demographic dynamics of alpine species: limited genetic consequences despite dramatic range contraction in a plant from the Spanish Sierra Nevada

Anthropogenic global climate change is expected to cause severe range contractions among alpine plants. Alpine areas in the Mediterranean region are of special concern because of the high abundance of endemic species with narrow ranges. This study combined species distribution models, population structure analyses and Bayesian skyline plots to trace the past and future distribution and diversity of Linaria glacialis, an endangered narrow endemic species that inhabits summits of Sierra Nevada (Spain). The results showed that: (i) the habitat of this alpine‐Mediterranean species in Sierra Nevada suffered little changes during glacial and interglacial stages of late Quaternary; (ii) climatic oscillations in the last millennium (Medieval Warm Period and Little Ice Age) moderately affected the demographic trends of L. glacialis; (iii) future warming conditions will cause severe range contractions; and (iv) genetic diversity will not diminish at the same pace as the distribution range. As a consequence of the low population structure of this species, genetic impoverishment in the alpine zones of Sierra Nevada should be limited during range contraction. We conclude that maintenance of large effective population sizes via high mutation rates and high levels of gene flow may promote the resilience of alpine plant species when confronted with global warming.