Oceanography and life history predict contrasting genetic population structure in two Antarctic fish species

Combining population genomics and biophysical modelling to assess connectivity patterns in an Antarctic fish

Connectivity is a fundamental process of population dynamics in marine ecosystems. In the last decade, with the emergence of new methods, combining different approaches to understand the patterns of connectivity among populations and their regulation has become increasingly feasible. The Western Antarctic Peninsula (WAP) is characterized by complex oceanographic dynamics, where local conditions could act as barriers to population connectivity. Here, the notothenioid fish Harpagifer antarcticus, a demersal species with a complex life cycle (adults with poor swim capabilities and pelagic larvae), was used to assess connectivity along the WAP by combining biophysical modelling and population genomics methods. Both approaches showed congruent patterns. Areas of larvae retention and low potential connectivity, observed in the biophysical model output, coincide with four genetic groups within the WAP: (1) South Shetland Islands, (2) Bransfield Strait, (3) the central and (4) the southern area of WAP (Marguerite Bay). These genetic groups exhibited limited gene flow between them, consistent with local oceanographic conditions, which would represent barriers to larval dispersal. The joint effect of geographic distance and larval dispersal by ocean currents had a greater influence on the observed population structure than each variable evaluated separately. The combined effect of geographic distance and a complex oceanographic dynamic would be generating limited levels of population connectivity in the fish H. antarcticus along the WAP. Based on this, population connectivity estimations and priority areas for conservation were discussed, considering the marine protected area proposed for this threatened region of the Southern Ocean.

Genomics of Secondarily Temperate Adaptation in the Only Non-Antarctic Icefish

White-blooded Antarctic icefishes, a family within the adaptive radiation of Antarctic notothenioid fishes, are an example of extreme biological specialization to both the chronic cold of the Southern Ocean and life without hemoglobin. As a result, icefishes display derived physiology that limits them to the cold and highly oxygenated Antarctic waters. Against these constraints, remarkably one species, the pike icefish Champsocephalus esox, successfully colonized temperate South American waters. To study the genetic mechanisms underlying secondarily temperate adaptation in icefishes, we generated chromosome-level genome assemblies of both C. esox and its Antarctic sister species, Champsocephalus gunnari. The C. esox genome is similar in structure and organization to that of its Antarctic congener; however, we observe evidence of chromosomal rearrangements coinciding with regions of elevated genetic divergence in pike icefish populations. We also find several key biological pathways under selection, including genes related to mitochondria and vision, highlighting candidates behind temperate adaptation in C. esox. Substantial antifreeze glycoprotein (AFGP) pseudogenization has occurred in the pike icefish, likely due to relaxed selection following ancestral escape from Antarctica. The canonical AFGP locus organization is conserved in C. esox and C. gunnari, but both show a translocation of two AFGP copies to a separate locus, previously unobserved in cryonotothenioids. Altogether, the study of this secondarily temperate species provides an insight into the mechanisms underlying adaptation to ecologically disparate environments in this otherwise highly specialized group.

Species' attributes predict the relative magnitude of ecological and genetic recovery following mass mortality

Theoretically, species' characteristics should allow estimation of dispersal potential and, in turn, explain levels of population genetic differentiation. However, a mismatch between traits and genetic patterns is often reported for marine species, and interpreted as evidence that life-history traits do not influence dispersal. Here, we couple ecological and genomic methods to test the hypothesis that species with attributes favouring greater dispersal potential-e.g., longer pelagic duration, higher fecundity and larger population size-have greater realized dispersal overall. We used a natural experiment created by a large-scale and multispecies mortality event which created a "clean slate" on which to study recruitment dynamics, thus simplifying a usually complex problem. We surveyed four species of differing dispersal potential to quantify the abundance and distribution of recruits and to genetically assign these recruits to probable parental sources. Species with higher dispersal potential recolonized a broader extent of the impacted range, did so more quickly and recovered more genetic diversity than species with lower dispersal potential. Moreover, populations of taxa with higher dispersal potential exhibited more immigration (71%-92% of recruits) than taxa with lower dispersal potential (17%-44% of recruits). By linking ecological with genomic perspectives, we demonstrate that a suite of interacting life-history and demographic attributes do influence species' realized dispersal and genetic neighbourhoods. To better understand species' resilience and recovery in this time of global change, integrative eco-evolutionary approaches are needed to more rigorously evaluate the effect of dispersal-linked attributes on realized dispersal and population genetic differentiation.

The Genetic Population Structure of Lake Tanganyika's Lates Species Flock, an Endemic Radiation of Pelagic Top Predators

Understanding genetic connectivity plays a crucial role in species conservation decisions, and genetic connectivity is an important component of modern fisheries management. In this study, we investigated the population genetics of four endemic Lates species of Lake Tanganyika (Lates stappersii, L. microlepis, L. mariae, and L. angustifrons) using reduced-representation genomic sequencing methods. We find the four species to be strongly differentiated from one another (mean interspecific FST = 0.665), with no evidence for contemporary admixture. We also find evidence for strong genetic structure within L. mariae, with the majority of individuals from the most southern sampling site forming a genetic group that is distinct from the individuals at other sampling sites. We find evidence for much weaker structure within the other three species (L. stappersii, L. microlepis, and L. angustifrons). Our ability to detect this weak structure despite small and unbalanced sample sizes and imprecise geographic sampling locations suggests the possibility for further structure undetected in our study. We call for further research into the origins of the genetic differentiation in these four species-particularly that of L. mariae-which may be important for conservation and management of this culturally and economically important clade of fishes.

Productivity and Change in Fish and Squid in the Southern Ocean

Southern Ocean ecosystems are globally important and vulnerable to global drivers of change, yet they remain challenging to study. Fish and squid make up a significant portion of the biomass within the Southern Ocean, filling key roles in food webs from forage to mid-trophic species and top predators. They comprise a diverse array of species uniquely adapted to the extreme habitats of the region. Adaptations such as antifreeze glycoproteins, lipid-retention, extended larval phases, delayed senescence, and energy-conserving life strategies equip Antarctic fish and squid to withstand the dark winters and yearlong subzero temperatures experienced in much of the Southern Ocean. In addition to krill exploitation, the comparatively high commercial value of Antarctic fish, particularly the lucrative toothfish, drives fisheries interests, which has included illegal fishing. Uncertainty about the population dynamics of target species and ecosystem structure and function more broadly has necessitated a precautionary, ecosystem approach to managing these stocks and enabling the recovery of depleted species. Fisheries currently remain the major local driver of change in Southern Ocean fish productivity, but global climate change presents an even greater challenge to assessing future changes. Parts of the Southern Ocean are experiencing ocean-warming, such as the West Antarctic Peninsula, while other areas, such as the Ross Sea shelf, have undergone cooling in recent years. These trends are expected to result in a redistribution of species based on their tolerances to different temperature regimes. Climate variability may impair the migratory response of these species to environmental change, while imposing increased pressures on recruitment. Fisheries and climate change, coupled with related local and global drivers such as pollution and sea ice change, have the potential to produce synergistic impacts that compound the risks to Antarctic fish and squid species. The uncertainty surrounding how different species will respond to these challenges, given their varying life histories, environmental dependencies, and resiliencies, necessitates regular assessment to inform conservation and management decisions. Urgent attention is needed to determine whether the current management strategies are suitably precautionary to achieve conservation objectives in light of the impending changes to the ecosystem.

Climatic change drives dynamic source-sink relationships in marine species with high dispersal potential

While there is now strong evidence that many factors can shape dispersal, the mechanisms influencing connectivity patterns are species-specific and remain largely unknown for many species with a high dispersal potential. The rock lobsters Jasus tristani and Jasus paulensis have a long pelagic larval duration (up to 20 months) and inhabit seamounts and islands in the southern Atlantic and Indian Oceans, respectively. We used a multidisciplinary approach to assess the genetic relationships between J. tristani and J. paulensis, investigate historic and contemporary gene flow, and inform fisheries management. Using 17,256 neutral single nucleotide polymorphisms we found low but significant genetic differentiation. We show that patterns of connectivity changed over time in accordance with climatic fluctuations. Historic migration estimates showed stronger connectivity from the Indian to the Atlantic Ocean (influenced by the Agulhas Leakage). In contrast, the individual-based model coupled with contemporary migration estimates inferred from genetic data showed stronger inter-ocean connectivity in the opposite direction from the Atlantic to the Indian Ocean driven by the Subtropical Front. We suggest that the J. tristani and J. paulensis historical distribution might have extended further north (when water temperatures were lower) resulting in larval dispersal between the ocean basis being more influenced by the Agulhas Leakage than the Subtropical Front. As water temperatures in the region increase in accordance with anthropogenic climate change, a southern shift in the distribution range of J. tristani and J. paulensis could further reduce larval transport from the Indian to the Atlantic Ocean, adding complexity to fisheries management.

Connectivity and population structure of albacore tuna across southeast Atlantic and southwest Indian Oceans inferred from multidisciplinary methodology

Albacore tuna (Thunnus alalunga) is an important target of tuna fisheries in the Atlantic and Indian Oceans. The commercial catch of albacore is the highest globally among all temperate tuna species, contributing around 6% in weight to global tuna catches over the last decade. The accurate assessment and management of this heavily exploited resource requires a robust understanding of the species’ biology and of the pattern of connectivity among oceanic regions, yet Indian Ocean albacore population dynamics remain poorly understood and its level of connectivity with the Atlantic Ocean population is uncertain. We analysed morphometrics and genetics of albacore (n = 1,874) in the southwest Indian (SWIO) and southeast Atlantic (SEAO) Oceans to investigate the connectivity and population structure. Furthermore, we examined the species’ dispersal potential by modelling particle drift through major oceanographic features. Males appear larger than females, except in South African waters, yet the length–weight relationship only showed significant male–female difference in one region (east of Madagascar and Reunion waters). The present study produced a genetic differentiation between the southeast Atlantic and southwest Indian Oceans, supporting their demographic independence. The particle drift models suggested dispersal potential of early life stages from SWIO to SEAO and adult or sub-adult migration from SEAO to SWIO.

Fine-scale seascape genomics of an exploited marine species, the common cockle Cerastoderma edule, using a multimodelling approach

Population dynamics of marine species that are sessile as adults are driven by oceanographic dispersal of larvae from spawning to nursery grounds. This is mediated by life-history traits such as the timing and frequency of spawning, larval behaviour and duration, and settlement success. Here, we use 1725 single nucleotide polymorphisms (SNPs) to study the fine-scale spatial genetic structure in the commercially important cockle species Cerastoderma edule and compare it to environmental variables and current-mediated larval dispersal within a modelling framework. Hydrodynamic modelling employing the NEMO Atlantic Margin Model (AMM15) was used to simulate larval transport and estimate connectivity between populations during spawning months (April-September), factoring in larval duration and interannual variability of ocean currents. Results at neutral loci reveal the existence of three separate genetic clusters (mean F ST = 0.021) within a relatively fine spatial scale in the north-west Atlantic. Environmental association analysis indicates that oceanographic currents and geographic proximity explain over 20% of the variance observed at neutral loci, while genetic variance (71%) at outlier loci was explained by sea surface temperature extremes. These results fill an important knowledge gap in the management of a commercially important and overexploited species, bringing us closer to understanding the role of larval dispersal in connecting populations at a fine geographic scale.

Unraveling hierarchical genetic structure in a marine metapopulation: A comparison of three high-throughput genotyping approaches

Marine metapopulations often exhibit subtle population structure that can be difficult to detect. Given recent advances in high-throughput sequencing, an emerging question is whether various genotyping approaches, in concert with improved sampling designs, will substantially improve our understanding of genetic structure in the sea. To address this question, we explored hierarchical patterns of structure in the coral reef fish Elacatinus lori using a high-resolution approach with respect to both genetic and geographic sampling. Previously, we identified three putative E. lori populations within Belize using traditional genetic markers and sparse geographic sampling: barrier reef and Turneffe Atoll; Glover's Atoll; and Lighthouse Atoll. Here, we systematically sampled individuals at ~10 km intervals throughout these reefs (1,129 individuals from 35 sites) and sequenced all individuals at three sets of markers: 2,418 SNPs; 89 microsatellites; and 57 nonrepetitive nuclear loci. At broad spatial scales, the markers were consistent with each other and with previous findings. At finer spatial scales, there was new evidence of genetic substructure, but our three marker sets differed slightly in their ability to detect these patterns. Specifically, we found subtle structure between the barrier reef and Turneffe Atoll, with SNPs resolving this pattern most effectively. We also documented isolation by distance within the barrier reef. Sensitivity analyses revealed that the number of loci (and alleles) had a strong effect on the detection of structure for all three marker sets, particularly at small spatial scales. Taken together, these results illustrate empirically that high-throughput genotyping data can elucidate subtle genetic structure at previously-undetected scales in a dispersive marine fish.

Genetic diversity and spatial population structure of a deepwater snapper, Pristipomoides filamentosus in the south-west Indian Ocean

The crimson jobfish, Pristipomoides filamentosus Valenciennes, 1830 is an economically important demersal species largely distributed in the Indo-Pacific region. Pristipomoides filamentosus constitutes a significant portion of catch landed in demersal fisheries throughout the species' distribution range. Despite the species' economic importance, there is insufficient data to guide the species' conservation management, especially within the south-western (SW) Indian Ocean region. The aims of the present study were to conduct a population genetic analysis to determine the spatial genetic structure of the species and, whether different management units could be established in the region, using an analysis of both mitochondrial DNA fragment (mtDNA), and nuclear microsatellite loci. A total of 193 fin clips were collected from Seychelles, Kenya, Tanzania, Comoros, Madagascar, Mauritius and South Africa, with each having an established fishery of the species. Both haplotype diversity (h) and expected heterozygosity (H_E) for mtDNA and microsatellite loci respectively were generally high for all localities, except for Seychelles where both diversity indices were at the lowest (i.e. h = 0.429 ± 0.134; H_E = 0.647 ± 0.059). Even though mtDNA failed to detect population differentiation, the hypervariable microsatellite loci consistently indicated presence of four genetic clusters irrespective of the clustering approach applied. Based on present results, we propose recognising the four clusters as distinct fisheries management units of the species in the SW Indian Ocean region.

Along-shelf connectivity and circumpolar gene flow in Antarctic silverfish (Pleuragramma antarctica)

The Antarctic silverfish (Pleuragramma antarctica) is a critically important forage species with a circumpolar distribution and is unique among other notothenioid species for its wholly pelagic life cycle. Previous studies have provided mixed evidence of population structure over regional and circumpolar scales. The aim of the present study was to test the recent population hypothesis for Antarctic silverfish, which emphasizes the interplay between life history and hydrography in shaping connectivity. A total of 1067 individuals were collected over 25 years from different locations on a circumpolar scale. Samples were genotyped at fifteen microsatellites to assess population differentiation and genetic structuring using clustering methods, F-statistics, and hierarchical analysis of variance. A lack of differentiation was found between locations connected by the Antarctic Slope Front Current (ASF), indicative of high levels of gene flow. However, gene flow was significantly reduced at the South Orkney Islands and the western Antarctic Peninsula where the ASF is absent. This pattern of gene flow emphasized the relevance of large-scale circulation as a mechanism for circumpolar connectivity. Chaotic genetic patchiness characterized population structure over time, with varying patterns of differentiation observed between years, accompanied by heterogeneous standard length distributions. The present study supports a more nuanced version of the genetic panmixia hypothesis that reflects physical-biological interactions over the life history.

Challenging Dogma Concerning Biogeographic Patterns of Antarctica and the Southern Ocean

Antarctica is enormous, cold, remote, and particularly sensitive to climate change. Most biological research below 60°S has focused on the isolated nature of the biota and how organisms have adapted to the cold and ice. However, biogeographic patterns in Antarctica and the Southern Ocean, and the processes explaining how those patterns came about, still await adequate explanation. Both terrestrial and marine organisms have been influenced by climatic change (e.g., glaciation), physical phenomena (e.g., oceanic currents), and/or potential barriers to gene flow (e.g., steep thermal gradients). Whereas the Antarctic region contains diverse and complex marine communities, terrestrial systems tend to be comparatively simple with limited diversity. Here, we challenge the current dogma used to explain the diversity and biogeographic patterns present in the Antarctic. We assert that relatively modern processes within the last few million years, rather than geo-logical events that occurred in the Eocene and Miocene, account for present patterns of biodiversity in the region. Additionally, reproductive life history stages appear to have little influence in structuring genetic patterns in the Antarctic, as currents and glacial patterns are noted to be more important drivers of organismal patterns of distribution. Finally, we highlight the need for additional sampling, high-throughput genomic approaches, and broad, multinational cooperation for addressing outstanding questions of Antarctic biogeography and biodiversity.

Geographic isolation and larval dispersal shape seascape genetic patterns differently according to spatial scale

Genetic variation, as a basis of evolutionary change, allows species to adapt and persist in different climates and environments. Yet, a comprehensive assessment of the drivers of genetic variation at different spatial scales is still missing in marine ecosystems. Here, we investigated the influence of environment, geographic isolation, and larval dispersal on the variation in allele frequencies, using an extensive spatial sampling (47 locations) of the striped red mullet (Mullus surmuletus) in the Mediterranean Sea. Univariate multiple regressions were used to test the influence of environment (salinity and temperature), geographic isolation, and larval dispersal on single nucleotide polymorphism (SNP) allele frequencies. We used Moran's eigenvector maps (db-MEMs) and asymmetric eigenvector maps (AEMs) to decompose geographic and dispersal distances in predictors representing different spatial scales. We found that salinity and temperature had only a weak effect on the variation in allele frequencies. Our results revealed the predominance of geographic isolation to explain variation in allele frequencies at large spatial scale (>1,000 km), while larval dispersal was the major predictor at smaller spatial scale (<1,000 km). Our findings stress the importance of including spatial scales to understand the drivers of spatial genetic variation. We suggest that larval dispersal allows to maintain gene flows at small to intermediate scale, while at broad scale, genetic variation may be mostly shaped by adult mobility, demographic history, or multigenerational stepping-stone dispersal. These findings bring out important spatial scale considerations to account for in the design of a protected area network that would efficiently enhance protection and persistence capacity of marine species.

Stepping stones to isolation: Impacts of a changing climate on the connectivity of fragmented fish populations

In the marine environment, understanding the biophysical mechanisms that drive variability in larval dispersal and population connectivity is essential for estimating the potential impacts of climate change on the resilience and genetic structure of populations. Species whose populations are small, isolated and discontinuous in distribution will differ fundamentally in their response and resilience to environmental stress, compared with species that are broadly distributed, abundant and frequently exchange conspecifics. Here, we use an individual-based modelling approach, combined with a population genetics projection model, to consider the impacts of a warming climate on the population connectivity of two contrasting Antarctic fish species, Notothenia rossii and Champsocephalus gunnari. Focussing on the Scotia Sea region, sea surface temperatures are predicted to increase significantly by the end of the 21st century, resulting in reduced planktonic duration and increased egg and larval mortality. With shorter planktonic durations, the results of our study predict reduced dispersal of both species across the Scotia Sea, from Antarctic Peninsula sites to islands in the north and east, and increased dispersal among neighbouring sites, such as around the Antarctic Peninsula. Increased mortality modified the magnitude of population connectivity but had little effect on the overall patterns. Whilst the predicted changes in connectivity had little impact on the projected regional population genetic structure of N. rossii, which remained broadly genetically homogeneous within distances of ~1,500 km, the genetic isolation of C. gunnari populations in the northern Scotia Sea was predicted to increase with rising sea temperatures. Our study highlights the potential for increased isolation of island populations in a warming world, with implications for the resilience of populations and their ability to adapt to ongoing environmental change, a matter of high relevance to fisheries and ecosystem-level management.

Independent estimates of marine population connectivity are more concordant when accounting for uncertainties in larval origins

Marine larval dispersal is a complex biophysical process that depends on the effects of species biology and oceanography, leading to logistical difficulties in estimating connectivity among populations of marine animals with biphasic life cycles. To address this challenge, the application of multiple methodological approaches has been advocated, in order to increase confidence in estimates of population connectivity. However, studies seldom account for sources of uncertainty associated with each method, which undermines a direct comparative approach. In the present study we explicitly account for the statistical uncertainty in observed connectivity matrices derived from elemental chemistry of larval mussel shells, and compare these to predictions from a biophysical model of dispersal. To do this we manipulate the observed connectivity matrix by applying different confidence levels to the assignment of recruits to source populations, while concurrently modelling the intrinsic misclassification rate of larvae to known sources. We demonstrate that the correlation between the observed and modelled matrices increases as the number of observed recruits classified as unknowns approximates the observed larval misclassification rate. Using this approach, we show that unprecedented levels of concordance in connectivity estimates (r = 0.96) can be achieved, and at spatial scales (20-40 km) that are ecologically relevant.

Closely related octopus species show different spatial genetic structures in response to the Antarctic seascape

Determining whether comparable processes drive genetic divergence among marine species is relevant to molecular ecologists and managers alike. Sympatric species with similar life histories might be expected to show comparable patterns of genetic differentiation and a consistent influence of environmental factors in shaping divergence. We used microsatellite loci to quantify genetic differentiation across the Scotia Arc in three species of closely related benthic octopods, Pareledone turqueti, P. charcoti, and Adelieledone polymorpha. The relative importance of environmental factors (latitude, longitude, depth, and temperature) in shaping genetic structure was investigated when significant spatial genetic structure was uncovered. Isolated populations of P. turqueti and A. polymorpha at these species' range margins were genetically different to samples close to mainland Antarctica; however, these species showed different genetic structures at a regional scale. Samples of P. turqueti from the Antarctic Peninsula, Elephant Island, and Signy Island were genetically different, and this divergence was associated primarily with sample collection depth. By contrast, weak or nonsignificant spatial genetic structure was evident across the Antarctic Peninsula, Elephant Island, and Signy Island region for A. polymorpha, and slight associations between population divergence and temperature or depth (and/or longitude) were detected. Pareledone charcoti has a limited geographic range, but exhibited no genetic differentiation between samples from a small region of the Scotia Arc (Elephant Island and the Antarctic Peninsula). Thus, closely related species with similar life history strategies can display contrasting patterns of genetic differentiation depending on spatial scale; moreover, depth may drive genetic divergence in Southern Ocean benthos.

Influence of Benthic Macrofauna as a Spatial Structuring Agent for Juvenile Haddock (Melanogrammus aeglefinus) on the Eastern Scotian Shelf, Atlantic Canada

We examined the habitat of juvenile haddock on the eastern Scotian Shelf (off Nova Scotia, Canada) in relation to grab-sampled benthic macrofaunal invertebrate species assemblages in order to determine whether there were significant differences in benthic macrofauna between areas of historically persistent high and low juvenile haddock abundance. Our analyses were conducted over two spatial scales in each of two years: among banks (Emerald, Western and Sable Island), approximately 60 km distant from each other, and between areas of high and low juvenile haddock abundance at distances of 10 to 30 km-all in an area that had not experienced groundfishing in the decade prior to sampling. We also examined fine-scale (10s of metres) within-site variability in the macrofauna and used surficial sediment characteristics, along with hydrographic variables, to identify environmental correlates. PERMANOVA identified statistically significant differences in biomass, density and composition of the benthos associated with juvenile haddock abundance; however it was difficult to determine whether the results had biological relevance. Post hoc tests showed that these differences occurred only on Sable Island Bank where both fish and benthos may have been independently responding to sediment type which was most different there (100% sand in the area of low haddock abundance vs. 22% gravel in the area of high haddock abundance). In total, 383 benthic taxa representing 13 phyla were identified. Annelida was the most specious phylum (36.29% of taxa, representing 33 families), followed by Arthropoda (with Crustaceans, mostly Amphipoda, accounting for 25.07% of the total number of taxa). The strongest pattern in the macrofauna was expressed at the largest scale, between banks, accounting for approximately 25% of the variation in the data. Emerald Bank, deeper, warmer and saltier than the Western and Sable Island Banks, had a distinctive fauna.

Non-Antarctic notothenioids: Past phylogenetic history and contemporary phylogeographic implications in the face of environmental changes

The non-Antarctic Notothenioidei families, Bovichtidae, Pseudaphritidae and Eleginopsidae, diverged early from the main notothenioid lineage. They are important in clarifying the early evolutionary processes that triggered notothenioid evolution in the Antarctic. The early-diverged group represents 8% of all notothenioid species and never established themselves on the Antarctic shelf. Most attention has been paid to the Antarctic notothenioids and their limited physiological tolerance to climate change and increased temperatures. In this review, we discuss key life history traits that are characteristic of the non-Antarctic early-diverged notothenioid taxa as well as the genetic resources and population differentiation information available for this group. We emphasise the population fitness and dynamics of these species and indicate how resource management and conservation of the group can be strengthened through an integrative approach. Both Antarctic waters and the non-Antarctic regions face rapid temperature rises combined with strong anthropogenic exploitation. While it is expected that early-diverged notothenioid species may have physiological advantages over high Antarctic species, it is difficult to predict how climate changes might alter the geographic range, behaviour, phenology and ultimately genetic variability of these species. It is possible, however, that their high degree of endemism and dependence on local environmental specificities to complete their life cycles might enhance their vulnerability.