Local niches explain coexistence in environmentally-distinct contact zones between Western Mediterranean vipers.

Species' ecological niches are frequently analysed to gain insights into how anthropogenic changes affect biodiversity. Coping with these changes often involves shifts in niche expression, which can disrupt local biotic interactions. Secondary contact zones, where competition and ecological segregation commonly occur, are ideal for studying the ecological factors influencing species' niches. In this study, we investigated the effect of climate and landscape factors on the ecological niches of two viper species, Vipera aspis and Vipera latastei, across three contact zones in northern Iberia, characterized by varying levels of landscape alteration. Using niche overlap tests, ecological niche models and spatial analyses we observed local variation in the expression of the species' niches across the three contact zones, resulting from the different abiotic and biotic conditions of each area. Rather than spatial niche segregation, we observed high niche overlap, suggesting niche convergence. A pattern of asymmetrical niche variation was identified in all contact zones, driven by species' climatic tolerances and the environmental conditions of each area. V. aspis generally exhibited a wider niche, except in the southernmost zone where the pure Mediterranean climate favored V. latastei. Human-induced landscape changes intensified niche asymmetry, by favoring the most generalist V. aspis over the specialist V. latastei, increasing habitat overlap, and likely competition. This study presents a comprehensive analysis of niche expression at range margins, anticipating a heightened impact of landscape changes in V. latastei. The methodological framework implemented here, and our findings, hold significant relevance for biodiversity management and conservation in human-impacted areas.

Environmental and topographic drivers of amphibian phylogenetic diversity and endemism in the Iberian Peninsula.

Understanding the ecological and evolutionary processes driving biodiversity patterns and allowing their persistence is of utmost importance. Many hypotheses have been proposed to explain spatial diversity patterns, including water-energy availability, habitat heterogeneity, and historical climatic refugia. The main goal of this study is to identify if general spatial drivers of species diversity patterns of phylogenetic diversity (PD) and phylogenetic endemism (PE) at the global scale are also predictive of PD and PE at regional scales, using Iberian amphibians as a case study. Our main hypothesis assumes that topography along with contemporary and historical climate are drivers of phylogenetic diversity and endemism, but that the strength of these predictors may be weaker at the regional scale than it tends to be at the global scale. We mapped spatial patterns of Iberian amphibians' phylogenetic diversity and endemism, using previously published phylogenetic and distribution data. Furthermore, we compiled spatial data on topographic and climatic variables related to the water-energy availability, topography, and historical climatic instability hypotheses. To test our hypotheses, we used Spatial Autoregressive Models and selected the best model to explain diversity patterns based on Akaike Information Criterion. Our results show that, out of the variables tested in our study, water-energy availability and historical climate instability are the most important drivers of amphibian diversity in Iberia. However, as predicted, the strength of these predictors in our case study is weaker than it tends to be at global scales. Thus, additional drivers should also be investigated and we suggest caution when interpreting these predictors as surrogates for different components of diversity.

Diversity, distribution and conservation of land mammals in Mauritania, North-West Africa.

Detailed knowledge about biodiversity distribution is critical for monitoring the biological effects of global change processes. Biodiversity knowledge gaps hamper the monitoring of conservation trends and they are especially evident in the desert biome. Mauritania constitutes a remarkable example on how remoteness and regional insecurity affect current knowledge gaps. Mammals remain one of the least studied groups in this country, without a concerted species checklist, the mapping of regions concentrating mammal diversity, or a national assessment of their conservation status. This work assessed the diversity, distribution, and conservation of land mammals in Mauritania. A total of 6,718 published and original observations were assembled in a spatial database and used to update the occurrence status, distribution area, and conservation status. The updated taxonomic list comprises 107 species, including 93 extant, 12 Regionally Extinct, and 2 Extinct in the Wild. Mapping of species distributions allowed locating concentrations of extant mammal species richness in coastal areas, along the Senegal River valley, and in mountain plateaus. Recent regional extinction of large-sized Artiodactyla and Carnivora has been very high (11% extinct species). From the extant mammals, 11% are threatened, including flagship species (e.g., Addax nasomaculatus and Panthera pardus). Species richness is poorly represented by the current protected areas. Despite the strong advances made, 23% of species categorise as Data Deficient. Persisting systematics and distribution uncertainties require further research. Field surveys in currently unexplored areas (northern and south-eastern regions) are urgently needed to increase knowledge about threatened mammals. The long-term conservation of land mammals in Mauritania is embedded in a complex web of socioeconomic and environmental factors that call for collaborative action and investment in sustainable human development. The current work sets the baseline for the future development of detailed research studies and to address the general challenges faced by mammals and biodiversity in the country.

Landscape resistance constrains hybridization across contact zones in a reproductively and morphologically polymorphic salamander.

Explicitly accounting for phenotypic differentiation together with environmental heterogeneity is crucial to understand the evolutionary dynamics in hybrid zones. Species showing intra-specific variation in phenotypic traits that meet across environmentally heterogeneous regions constitute excellent natural settings to study the role of phenotypic differentiation and environmental factors in shaping the spatial extent and patterns of admixture in hybrid zones. We studied three environmentally distinct contact zones where morphologically and reproductively divergent subspecies of Salamandra salamandra co-occur: the pueriparous S. s. bernardezi that is mostly parapatric to its three larviparous subspecies neighbours. We used a landscape genetics framework to: (i) characterise the spatial location and extent of each contact zone; (ii) assess patterns of introgression and hybridization between subspecies pairs; and (iii) examine the role of environmental heterogeneity in the evolutionary dynamics of hybrid zones. We found high levels of introgression between parity modes, and between distinct phenotypes, thus demonstrating the evolution to pueriparity alone or morphological differentiation do not lead to reproductive isolation between these highly divergent S. salamandra morphotypes. However, we detected substantial variation in patterns of hybridization across contact zones, being lower in the contact zone located on a topographically complex area. We highlight the importance of accounting for spatial environmental heterogeneity when studying evolutionary dynamics of hybrid zones.

The demise of a wonder: Evolutionary history and conservation assessments of the Wonder Gecko Teratoscincus keyserlingii (Gekkota, Sphaerodactylidae) in Arabia.

Effective biodiversity conservation planning starts with genetic characterization within and among focal populations, in order to understand the likely impact of threats for ensuring the long-term viability of a species. The Wonder Gecko, Teratoscincus keyserlingii, is one of nine members of the genus. This species is distributed in Iran, Afghanistan, and Pakistan, with a small isolated population in the United Arab Emirates (UAE), where it is classified nationally as Critically Endangered. Within its Arabian range, anthropogenic activity is directly linked to the species' decline, with highly localised and severely fragmented populations. Here we describe the evolutionary history of Teratoscincus, by reconstructing its phylogenetic relationships and estimating its divergence times and ancestral biogeography. For conservation implications of T. keyserlingii we evaluate the genetic structure of the Arabian population using genomic data. This study supports the monophyly of most species and reveals considerable intraspecific variability in T. microlepis and T. keyserlingii, which necessitate broad systematic revisions. The UAE population of T. keyserlingii likely arrived from southern Iran during the Pleistocene and no internal structure was recovered within, implying a single population status. Regional conservation of T. keyserlingii requires improved land management and natural habitat restoration in the species' present distribution, and expansion of current protected areas, or establishment of new areas with suitable habitat for the species, mostly in northern Abu Dhabi Emirate.

Longitudinal sampling of external mucosae in farmed European seabass reveals the impact of water temperature on bacterial dynamics.

Fish microbiota are intrinsically linked to health and fitness, but they are highly variable and influenced by both biotic and abiotic factors. Water temperature particularly limits bacterial adhesion and growth, impacting microbial diversity and bacterial infections on the skin and gills. Aquaculture is heavily affected by infectious diseases, especially in warmer months, and industry practices often promote stress and microbial dysbiosis, leading to an increased abundance of potentially pathogenic bacteria. In this regard, fish mucosa health is extremely important because it provides a primary barrier against pathogens. We used 16 rRNA V4 metataxonomics to characterize the skin and gill microbiota of the European seabass, Dicentrarchus labrax, and the surrounding water over 12 months, assessing the impact of water temperature on microbial diversity and function. We show that the microbiota of external mucosae are highly dynamic with consistent longitudinal trends in taxon diversity. Several potentially pathogenic genera (Aliivibrio, Photobacterium, Pseudomonas, and Vibrio) were highly abundant, showing complex interactions with other bacterial genera, some of which with recognized probiotic activity, and were also significantly impacted by changes in temperature. The surrounding water temperature influenced fish microbial composition, structure and function over time (days and months). Additionally, dysbiosis was more frequent in warmer months and during transitions between cold/warm months. We also detected a strong seasonal effect in the fish microbiota, which is likely to result from the compound action of several unmeasured environmental factors (e.g., pH, nutrient availability) beyond temperature. Our results highlight the importance of performing longitudinal studies to assess the impact of environmental factors on fish microbiotas.

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt.

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south-western and south-eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south-western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south-eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150-km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2-27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population-level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.

Genes on the edge: A framework to detect genetic diversity imperiled by climate change.

Ongoing global warming is disrupting several ecological and evolutionary processes, spanning different levels of biological organization. Species are expected to shift their ranges as a response to climate change, with relevant implications to peripheral populations at the trailing and leading edges. Several studies have analyzed the exposure of species to climate change but few have explored exposure at the intraspecific level. We introduce a framework to forecast exposure to climate change at the intraspecific level. We build on existing methods by combining correlative species distribution models, a model of species range dynamics, and a model of phylogeographic interpolation. We demonstrate the framework by applying it to 20 Iberian amphibian and reptile species. Our aims were to: (a) identify which species and intraspecific lineages will be most exposed to future climate change; (b) test if nucleotide diversity at the edges of species ranges are significantly higher or lower than on the overall range; and (c) analyze if areas of higher species gain, loss, and turnover coincide with those predicted for lineages richness and nucleotide diversity. We found that about 80% of the studied species are predicted to contract their range. Within each species, some lineages were predicted to contract their range, while others were predicted to maintain or expand it. Therefore, estimating the impacts of climate change at the species level only can underestimate losses at the intraspecific level. Some species had significant high amount of nucleotide at the trailing or leading edge, or both, but we did not find a consistent pattern across species. Spatial patterns of species richness, gain, loss, and turnover were fairly concurrent with lineages richness and nucleotide diversity. Our results support the need for increased attention to intraspecific diversity regarding monitoring and conservation strategies under climate change.

An integrative assessment of the diversity, phylogeny, distribution, and conservation of the terrestrial reptiles (Sauropsida, Squamata) of the United Arab Emirates.

In the present study we use an unprecedented database of 5,535 distributional records to infer the diversity, ecological preferences and spatial distribution of the 60 species of terrestrial reptiles of the United Arab Emirates (UAE), and use the 57 native species to test the effectiveness of the protected areas' network in conserving this unique vertebrate fauna. We infer a time-calibrated phylogeny with 146 species of squamates and 15 genes including all UAE terrestrial reptile species to determine the phylogenetic diversity (PD) and evolutionary distinctiveness (ED) of the native species and to compare it with the distribution of the hotspots of native species richness. The results of this study indicate that the sampling effort is remarkable, covering 75% of the country's territory representing nearly the entire climatic space of the UAE defined by the mean annual temperature and the total annual precipitation, as well as the multivariate climatic space defined by a principal component analysis (PCA). Species richness is highest in the northeast of the country, in a transitional area from sandy desert to the mountainous terrain of the Hajar Mountains. The highest PD of a single square cell of 10 arc-minutes grid is of 2,430 million years (my) of accumulated evolutionary history and the strong correlation between PD and species richness suggests that the raw number of species is a good surrogate to quantify the evolutionary history (i.e., PD). The species with the highest values of ED are those in families represented by only one species in the UAE. Finally, the assessment of the UAE protected areas shows that, despite their relevance in protecting the terrestrial reptiles, they do not offer adequate protection for some threatened species. Therefore, a reassessment of some of the protected areas or the creation of species specific conservation action plans are recommended in order to ensure the preservation of the unique diversity of UAE terrestrial reptiles.

Phylin 2.0: Extending the phylogeographical interpolation method to include uncertainty and user-defined distance metrics.

Estimating geographical ranges of intra-specific evolutionary lineages is crucial to the fields of biogeography, evolution, and biodiversity conservation. Models of isolation mechanisms often consider multiple distances in order to explain genetic divergence. Yet, the available methods to estimate the geographical ranges of lineages are based on direct geographical distances, neglecting other distance metrics that can better explain the spatial genetic structure. We extended the phylogeographical interpolation method (phylin) in order to accommodate user-defined distance metrics and to incorporate the uncertainty associated with genetic distance calculation. These new features were tested with simulated and empirical data sets. Multiple distance matrices were generated including geographical, resistance, and environmental distances to derive maps of lineage occurrence. The new additions to this method improved the ability to predict lineage occurrence, even with low sample size. We used a regression framework to quantify the relationship between the genetic divergence and competing distance matrices representing potential isolation processes that are subsequently used in the interpolation process. Including uncertainty in tree topology and the different distance matrices improved the robustness of the variograms, allowing a better fit of the theoretical model of spatial dependence. The improvements to the method increase its potential application in other fields. Accurately mapping genetic divergence can help to locate potential contact zones between lineages as well as barriers to gene flow, which has a broad interest in biogeographical and evolutionary studies. Additionally, conservation efforts could benefit from the integration of genetic variation and landscape features in a spatially explicit framework.

Combining molecular and landscape tools for targeting evolutionary processes in reserve design: An approach for islands.

The importance of targeting ecological and evolutionary processes in reserve design has been widely acknowledged in the literature but rarely implemented on islands. Using Socotran reptiles as models, we aim to relate richness of widespread and restricted-range species directly with landscape variables and to compare the impact of setting conservation targets for lineages versus species. Socotra Island is a UNESCO Natural World Heritage Site, containing high levels of endemism in relation to its area, especially of reptiles, the vertebrates with the most comprehensive available genetic data. We predicted the occurrences of reptile species using distribution models and used a novel approach to interpolate maps of spatial phylogenetic patterns. Patterns of intra and interspecifc diversity and differences between spatial outputs of lineage and species richness were related to eco-geographic variables. We evaluated differences in target achievement for each conservation unit within protected areas (PAs) under the current Zoning Plan (ZP) using gap and reserve design analyses. Although intraspecific richness was strongly correlated with interspecific richness, differences in their spatial distribution reached ~30% in some areas. Differences were more pronounced for wide-ranging than restricted-range taxa. Gap analysis indicates that most conservation units are under-represented in sanctuaries and that intra and interspecific richness were significantly higher outside PAs. This work will guide local-scale conservation planning as the ZP is due to be re-evaluated. This is one of the few studies on islands using genetic data from an entire class of vertebrates to incorporate lineage diversity in reserve design. This study provides an alternative methodological framework for supporting the use of landscape and genetic tools in reserve design, circumventing the use of phylogenetic distances and deterministic spatial interpolation of lineage diversity that can be widely applied to other systems.

Diversity, distribution and conservation of the terrestrial reptiles of Oman (Sauropsida, Squamata).

In the present work, we use an exceptional database including 5,359 records of 101 species of Oman's terrestrial reptiles together with spatial tools to infer the spatial patterns of species richness and endemicity, to infer the habitat preference of each species and to better define conservation priorities, with especial focus on the effectiveness of the protected areas in preserving this unique arid fauna. Our results indicate that the sampling effort is not only remarkable from a taxonomic point of view, with multiple observations for most species, but also for the spatial coverage achieved. The observations are distributed almost continuously across the two-dimensional climatic space of Oman defined by the mean annual temperature and the total annual precipitation and across the Principal Component Analysis (PCA) of the multivariate climatic space and are well represented within 17 out of the 20 climatic clusters grouping 10% of the explained climatic variance defined by PC1 and PC2. Species richness is highest in the Hajar and Dhofar Mountains, two of the most biodiverse areas of the Arabian Peninsula, and endemic species richness is greatest in the Jebel Akhdar, the highest part of the Hajar Mountains. Oman's 22 protected areas cover only 3.91% of the country, including within their limits 63.37% of terrestrial reptiles and 50% of all endemics. Our analyses show that large areas of the climatic space of Oman lie outside protected areas and that seven of the 20 climatic clusters are not protected at all. The results of the gap analysis indicate that most of the species are below the conservation target of 17% or even the less restrictive 12% of their total area within a protected area in order to be considered adequately protected. Therefore, an evaluation of the coverage of the current network of protected areas and the identification of priority protected areas for reptiles using reserve design algorithms are urgently needed. Our study also shows that more than half of the species are still pending of a definitive evaluation by the International Union for Conservation of Nature (IUCN).

Spatial conservation prioritization of biodiversity spanning the evolutionary continuum.

Accounting for evolutionary relationships between and within species is important for biodiversity conservation planning, but is rarely considered in practice. Here we introduce a novel framework to identify priority conservation areas accounting for phylogenetic and intraspecific diversity, integrating concepts from phylogeny, phylogeography, spatial statistics and spatial conservation prioritization. The framework allows planners to incorporate and combine different levels of evolutionary diversity and can be applied to any taxonomic group and to any region in the world. We illustrate our approach using amphibian and reptile species occurring in a biodiversity hotspot region, the Iberian Peninsula. We found that explicitly incorporating phylogenetic and intraspecific diversity in systematic conservation planning provides advantages in terms of maximizing overall biodiversity representation while enhancing its persistence and evolutionary potential. Our results emphasize the need to account for the evolutionary continuum in order to efficiently implement biodiversity conservation planning decisions.

The complete mitochondrial genome of Pristurus rupestris rupestris.

Pristurus rupestris rupestris is a gecko of the family Sphaerodactylidae adapted to the arid habitat found in the Hajar Mountains of southeastern Arabia. The complete mitochondrial genome was obtained with Illumina sequencing. The sequenced mitogenome has 13 protein-coding genes, 22 tRNAs, two rRNA genes and two non-coding regions, totalling 16,993 bp. The AT content of the obtained sequence is 52.1% (A:28.7%, T:23.4%, G:14.7%, C:33.2%). The control region has an AT content of 54.3% and a length of 1558 bp.

Ecological and evolutionary influences on body size and shape in the Galápagos marine iguana (Amblyrhynchus cristatus).

Oceanic islands are often inhabited by endemic species that have undergone substantial morphological evolutionary change due to processes of multiple colonizations from various source populations, dispersal, and local adaptation. Galápagos marine iguanas are an example of an island endemic exhibiting high morphological diversity, including substantial body size variation among populations and sexes, but the causes and magnitude of this variation are not well understood. We obtained morphological measurements from marine iguanas throughout their distribution range. These data were combined with genetic and local environmental data from each population to investigate the effects of evolutionary history and environmental conditions on body size and shape variation and sexual dimorphism. Our results indicate that body size and shape are highly variable among populations. Sea surface temperature and island perimeter, but not evolutionary history as depicted by phylogeographic patterns in this species, explain variation in body size among populations. Conversely, evolutionary history, but not environmental parameters or island size, was found to influence variation in body shape among populations. Finally, in all populations except one, we found strong sexual dimorphism in body size and shape in which males are larger, with higher heads than females, while females have longer heads than males. Differences among populations suggest that plasticity and/or genetic adaptation may shape body size and shape variation in marine iguanas. This study will help target future investigations to address the contribution of plasticity versus genetic adaptation on size and shape variation in marine iguanas.

phylin: an r package for phylogeographic interpolation.

phylin is a package for the r programming environment which offers different methods to spatially interpolate genetic information from phylogeographic data. These interpolations can be used to predict the spatial occurrence of different lineages within a phylogeny using a modified method of kriging, which allows the usage of a genetic distance matrix to derive a model of spatial dependence. phylin improves the available methods to generate interpolated surfaces from a phylogenetic trees by assessing the autocorrelation structure of the genetic information, interpolating the genetic data based on a statistical model, estimating the uncertainty of the predictions and identifying lineage occurrence and contact zones probability without projection of pairwise genetic distances into mid-points between sample locations. The package also includes methods to plot interpolation surfaces and provide summary tables from the generated data and models. We provide an example of the usefulness of this tool by inferring the spatial occurrence of distinct historical evolutionary lineages of the Lataste's viper (Vipera latastei Boscá, 1878) in the Iberian Peninsula and identifying potential contact areas. The maps of phylogenetic patterns obtained with these methods provide a spatial context to test hypotheses related to processes underlying the geographic distribution of genetic diversity and to inform conservation planning.

Hybridization at an ecotone: ecological and genetic barriers between three Iberian vipers.

The formation of stable genetic boundaries between emerging species is often diagnosed by reduced hybrid fitness relative to parental taxa. This reduced fitness can arise from endogenous and/or exogenous barriers to gene flow. Although detecting exogenous barriers in nature is difficult, we can estimate the role of ecological divergence in driving species boundaries by integrating molecular and ecological niche modelling tools. Here, we focus on a three-way secondary contact zone between three viper species (Vipera aspis, V. latastei and V. seoanei) to test for the contribution of ecological divergence to the development of reproductive barriers at several species traits (morphology, nuclear DNA and mitochondrial DNA). Both the nuclear and mitochondrial data show that all taxa are genetically distinct and that the sister species V. aspis and V. latastei hybridize frequently and backcross over several generations. We find that the three taxa have diverged ecologically and meet at a hybrid zone coincident with a steep ecotone between the Atlantic and Mediterranean biogeographical provinces. Integrating landscape and genetic approaches, we show that hybridization is spatially restricted to habitats that are suboptimal for parental taxa. Together, these results suggest that niche separation and adaptation to an ecological gradient confer an important barrier to gene flow among taxa that have not achieved complete reproductive isolation.

Unravelling biodiversity, evolution and threats to conservation in the Sahara-Sahel.

Deserts and arid regions are generally perceived as bare and rather homogeneous areas of low diversity. The Sahara is the largest warm desert in the world and together with the arid Sahel displays high topographical and climatic heterogeneity, and has experienced recent and strong climatic oscillations that have greatly shifted biodiversity distribution and community composition. The large size, remoteness and long-term political instability of the Sahara-Sahel, have limited knowledge on its biodiversity. However, over the last decade, there have been an increasing number of published scientific studies based on modern geomatic and molecular tools, and broad sampling of taxa of these regions. This review tracks trends in knowledge about biodiversity patterns, processes and threats across the Sahara-Sahel, and anticipates needs for biodiversity research and conservation. Recent studies are changing completely the perception of regional biodiversity patterns. Instead of relatively low species diversity with distribution covering most of the region, studies now suggest a high rate of endemism and larger number of species, with much narrower and fragmented ranges, frequently limited to micro-hotspots of biodiversity. Molecular-based studies are also unravelling cryptic diversity associated with mountains, which together with recent distribution atlases, allows identifying integrative biogeographic patterns in biodiversity distribution. Mapping of multivariate environmental variation (at 1 km × 1 km resolution) of the region illustrates main biogeographical features of the Sahara-Sahel and supports recently hypothesised dispersal corridors and refugia. Micro-scale water-features present mostly in mountains have been associated with local biodiversity hotspots. However, the distribution of available data on vertebrates highlights current knowledge gaps that still apply to a large proportion of the Sahara-Sahel. Current research is providing insights into key evolutionary and ecological processes, including causes and timing of radiation and divergence for multiple taxa, and associating the onset of the Sahara with diversification processes for low-mobility vertebrates. Examples of phylogeographic patterns are showing the importance of allopatric speciation in the Sahara-Sahel, and this review presents a synthetic overview of the most commonly hypothesised diversification mechanisms. Studies are also stressing that biodiversity is threatened by increasing human activities in the region, including overhunting and natural resources prospection, and in the future by predicted global warming. A representation of areas of conflict, landmines, and natural resources extraction illustrates how human activities and regional insecurity are hampering biodiversity research and conservation. Although there are still numerous knowledge gaps for the optimised conservation of biodiversity in the region, a set of research priorities is provided to identify the framework data needed to support regional conservation planning.

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit.

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X-chromosome, Y-chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95-4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co-localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X- and Y-chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.