The impact of aquaculture on the genetics and distribution of the onuphid annelid Diopatra biscayensis

AimEvolutionary history of natural populations can be confounded by human intervention such as the case of decorator worm species Diopatra (Onuphidae), which have a history of being transported through anthropogenic activities. Because they build tubes and act as ecosystem engineers, they can have a large impact on the overall ecosystem in which they occur. One conspicuous member, Diopatra biscayensis, which was only described in 2012, has a fragmented distribution that includes the Bay of Biscay and the Normanno‐Breton Gulf in the English Channel. This study explores the origin of these worms in the Normanno‐Breton region, which has been debated to either be the result of a historic range contraction from a relic continuous population or a more recent introduction.LocationNortheastern Atlantic, the Bay of Biscay, and the Normanno‐Breton Gulf.MethodsWe utilized a RAD‐tag‐based SNP approach to create a reduced genomic data set to recover fine‐scale population structure and infer which hypothesis best describes the D. biscayensis biogeographic distribution. The reduced genomic data set was used to calculate standard genetic diversities and genetic differentiation statistics, and utilized various clustering analyses, including PCAs, DAPC, and admixture.ResultsClustering analyses were consistent with D. biscayensis as a single population spanning the Bay of Biscay to the Normanno‐Breton Gulf in the English Channel, although unexpected genetic substructure was recovered from Arcachon Bay, in the middle of its geographic range. Consistent with a hypothesized introduction, the isolated Sainte‐Anne locality in the Normanno‐Breton Gulf was recovered to be a subset of the diversity found in the rest of the Bay of Biscay.Main conclusionsThese results are congruent with previous simulations that did not support connectivity from the Bay of Biscay to the Normanno‐Breton Gulf by natural dispersal. These genomic findings, with support from previous climatic studies, further support the hypothesis that D. biscayensis phylogeographic connectivity is the result of introductions, likely through the regions’ rich shellfish aquaculture, and not of a historically held range contraction.     

Sex differentiation in seasonal distribution of the starry smooth-hound Mustelus asterias

This mark-recapture study of starry smooth-hound Mustelus asterias tagged during the summer months near the Dutch coast demonstrates a large-scale spatial sex differentiation in their circannual migration patterns and small-scale spatial sex differentiation during summer. Overwintering occurs in the North Sea, English Channel and Bay of Biscay, with significantly more males in the Northern North Sea and more females in the Bay of Biscay. During summer, sheltered sea arms off the Dutch coast were almost exclusively used by adult females. In subsequent summers post-release, both sexes were mostly confined to the Southern North Sea, suggesting philopatry.     

Interglacial genetic diversification of Moussonia deppeana (Gesneriaceae), a hummingbird‐pollinated,cloud forest shrub in northern Mesoamerica

Recent empirical work on cloud forest‐adapted species supports the role of both old divergences across major geographical areas and more recent divergences attributed to Pleistocene climate changes. The shrub Moussonia deppeana is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling throughout the species range and employing plastid and nuclear markers, we (i) test whether the fragmented distribution is correlated with main evolutionary lineages, (ii) reconstruct its phylogeographical history to infer the history of cloud forest in northern Mesoamerica and (iii) evaluate a set of refugia/vicariance scenarios for the region and demographic patterns of the populations whose ranges expanded and tracked cloud forest conditions during the Last Glacial Maximum. We found a deep evolutionary split in M. deppeana about 6–3 Ma, which could be consistent with a Pliocene divergence. Comparison of variation in plastid and nuclear markers revealed several lineages mostly congruent with their isolated geographical distribution and restricted gene flow among groups. Results of species distribution modelling and coalescent simulations fit a model of multiple refugia diverging during interglacial cycles. The demographic history of M. deppeana is not consistent with an expanding–contracting cloud forest archipelago model during the Last Glacial Maximum. Instead, our data suggest that populations persisted across the geographical range throughout the glacial cycles, and experienced isolation and divergence during interglacial periods.     

The not‐so‐Irish spurge: Euphorbia hyberna (Euphorbiaceae) and the Littletonian plant ‘steeplechase’

The disjunct distributions of the Lusitanian flora, which are found only in south‐west Ireland and northern Iberia, and are generally absent from intervening regions, have been of great interest to biogeographers. There has been much debate as to whether Irish populations represent relicts that survived the Last Glacial Maximum (LGM; approximately 21 kya), or whether they recolonized from southern refugia subsequent to the retreat of the ice and, if so, whether this occurred directly (i.e. the result of long distance dispersal) or successively (i.e. in the manner of a ‘steeplechase’, with the English Channel and Irish Sea representing successive ‘water‐jumps’ that have to be successfully crossed). In the present study, we used a combined palaeodistribution modelling and phylogeographical approach to determine the glacial history of the Irish spurge, Euphorbia hyberna, the sole member of the Lusitanian flora that is also considered to occur naturally in south‐western England. Our findings suggest that the species persisted through the LGM in several southern refugia, and that northern populations are the result of successive recolonization of Britain and Ireland during the postglacial Littletonian warm stage, akin to the ‘steeplechase’ hypothesis. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 249–259.     

Long‐distance dispersal or postglacial contraction? Insights into disjunction between Himalaya–Hengduan Mountains and Taiwan in a cold‐adapted herbaceous genus,Triplostegia

Current disjunct patterns can result from long‐distance dispersal or postglacial contraction. We herein investigate the evolutionary history of Triplostegia to elucidate the disjunction between the Himalaya–Hengduan Mountain region (HHM) and Taiwan (TW). Genetic structure of Triplostegia was investigated for 48 populations using sequences from five chloroplast loci and the ribosomal nuclear internal transcribed spacer. Divergence time estimation, ancestral area reconstruction, and species distribution modeling (SDM) were employed to examine the biogeographic history of Triplostegia. Substantial genetic differentiation among populations from southwestern China (SW), Central China (CC), and TW was detected. Triplostegia was inferred to have originated in SW, and diversification began during the late Miocene; CC was colonized in the mid‐Pliocene, and TW was finally colonized in the early Pleistocene. SDM suggested an expansion of climatically suitable areas during the Last Glacial Maximum and range contraction during the Last interglacial in Triplostegia. Disjunction between HHM and TW in Triplostegia is most likely the consequence of topographic isolation and postglacial contraction. The potential climatic suitability areas for Triplostegia by 2070s (2061–2080) are predicted to slightly shrink and move northward. With continued global warming and human‐induced deforestation, extinction risk may increase for the cold‐adapted species, and appropriate strategies should be employed for ecosystem conservation.     

Phylogeography of an estuarine mysid, Neomysis integer (Crustacea, Mysida), along the north-east Atlantic coasts

Aim The brackish water mysid, Neomysis integer, is one of the most common mysid species along the coasts of the north‐east Atlantic. In the present study, the phylogeographical patterns were examined throughout the distribution range of N. integer. In particular, the latitudinal trends in genetic diversity and the distribution of genetic variation were examined in order to elucidate the imprints of the Pleistocene glaciations. Location North‐east Atlantic coasts from the Baltic Sea to the south of Spain. Methods A total of 461 specimens from 11 populations were analysed by means of single‐stranded conformation polymorphism analysis combined with DNA sequencing of a fragment of the mitochondrial cytochrome c oxidase I gene. The genetic structure was examined by using a progression of phylogenetic, demographic and population genetic analyses to elucidate not only the geographical structure, but also the evolutionary history producing that structure. Results The levels of genetic diversity were relatively uniform throughout the distribution range, with the exception of a decline at the northern and southern edges of distribution. A high heterogeneity was observed between the populations analysed (global Φ_ST = 0.787). This is caused by the disparate distribution of the cytochrome oxidase I haplotypes, with several population‐specific haplotypes. A clear genetic break (2.4% sequence divergence) occurred between the southernmost Guadalquivir population and all other populations. Main conclusions The present study corroborates the expectations of the genetic patterns typically observed in an estuarine species. The within‐population variability was low, whereas a significant (moderate to high) divergence was observed between populations. Phylogeographical analysis revealed that northern populations within the English Channel, North Sea and Baltic Sea are characterized by several widespread haplotypes, while the Irish population and all sites south of the Bay of Biscay consist solely of unique haplotypes. This pattern, combined with the relative high levels of genetic diversity, could be indicative for the presence of a glacial refugium in the English Channel region. Under this scenario N. integer must have survived the Last Glacial Maximum in the palaeoriver system present in that region.     

Spatial variability in digenean metacercariae infection of 0-group common sole Solea solea among nurseries along the French Atlantic coast

We sampled 0-group sole juveniles (N = 174) in September and October 2003 in 9 major nurseries located along the French Atlantic coast (English Channel and Bay of Biscay). 0-group sole were infected with 3 genera of digenean metacercariae, Timoniella spp. (Acanthostomidae), Prosorhynchus crucibulum and Prosorhynchus sp. A (Bucephalidae), Podocotyle sp. (Opecoelidae), and an unidentified species of Digenea. Parasite infection levels in the English Channel nurseries were lower, and the community composition was different from Bay of Biscay nurseries. We hypothesize that the difference between geographic areas was due to differences in first intermediate host communities and, to a lesser extent, because of lower temperatures in the English Channel compared to the Bay of Biscay. For the Bay of Biscay nurseries, mean total parasite abundance was strongly negatively correlated with mean annual river flow. This relationship may be the result of the more upstream location of 0-group sole nurseries in estuaries so that fish were further from local points of parasite transmission than those in embayed nurseries. Digenean metacercariae load may be influenced by 3 major local factors, i.e. abundance and proximity of the first intermediate hosts and cercariae dispersal capacities.     

Morphological variability of the shape of striped red mullet Mullus surmuletus in relation to stock discrimination between the Bay of Biscay and the eastern English Channel

Truss analysis and length measurements were made on 168 striped red mullet Mullus surmuletus. Multivariate statistical analyses with principal component analysis and partial redundancy analysis (pRDA) were used on these measurements to evaluate the influence of maturity, sex and geographical area distribution on body shape. Truss measurements were important to quantify and discriminate changing body shape, presumably due to changing environmental conditions. Sexual dimorphism was not observed and juveniles could be distinguished from adults based on their body shape. More importantly, M. surmuletus occurring in different geographical areas could be differentiated using this method. Based on pRDA, a significant difference of head morphological dimensions was observed between populations occurring in the eastern English Channel and those occurring in the Bay of Biscay, suggesting that fish from these areas could represent two subpopulations.     

Seasonal migration of the starry smooth‐hound shark Mustelus asterias as revealed from tag‐recapture data of an angler‐led tagging programme

The primary aim of this long‐term angler‐led tagging programme was to gain information about seasonal changes in distribution of the starry smooth‐hound shark Mustelus asterias, along the Dutch coast for management and conservation purposes. Between 2011 and 2014, M. asterias comprised 92·6% (n = 2418) of the total elasmobranch catch (n = 2612) by the licenced group of taggers within the Dutch Delta of which 2244 M. asterias were fin‐tagged with plastic rototags. Sex and total length (L_T) composition inside the eastern tidal basin (Oosterschelde) were significantly different, i.e. more females and larger individuals, than outside indicating a pupping ground, which was confirmed by the capture of 30 newborn pups (≤32 cm). The distribution pattern of reported recaptured M. asterias (return‐rate 3·6%, n = 80) showed a circannual migration between summering in the southern North Sea and wintering in the English Channel and the Bay of Biscay, and suggests that M. asterias is philopatric. The Dutch angling season for M. asterias runs from approximately mid‐May to mid‐October when the water temperature is above 13° C. Recaptures of eight mature females, but no males in the Bay of Biscay, indicate partial spatial segregation by sex, where mature females migrate over larger distances than immature females and males. These observations, with the absence of recaptures in other known summering areas (i.e. the Irish Sea and Bristol Channel), suggest that the southern North Sea is used by a separate population. Implications for management and recommendations to improve and expand the study approach are discussed.     

Genetic diversity in Crambe maritima along the English Channel: the role of ocean currents in determining population structure

We utilised genetic markers (ISSR's) to investigate population structuring in Crambe maritima , a coastal perennial found in isolated populations on either side of the English Channel. Despite the expectation that smaller populations would contain lower levels of genetic diversity, we found no correlation with total population size or the number of flowering plants and genetic diversity. Although populations were genetically differentiated from each other, isolation by distance was not apparent in the whole data set or within the UK samples.
Despite being sampled from a wide geographic range, the majority of diversity (71%) was contained within populations, with no unique genotype, or even allele, identifying any population. When genetic distances were plotted using MDS techniques, samples from each population had a tendency to cluster, but the spread of points was wide indicating high levels of gene flow. High levels of gene flow, through the dispersal of seeds, were verified by assignment testing, which showed that 8.7% of individuals were assigned to a population other than the one they were collected from. The pattern of genetic variation can be explained by examining the direction of currents in the English Channel and Bay of Biscay. Gene flow, via seed dispersal, is possible between UK and French populations on either side of the English Channel. In contrast, French populations on the Bay of Biscay coast are effectively isolated by the direction of currents and consequently show a greater degree of genetic differentiation.     

Southern montane populations did not contribute to the recolonization of West Siberian Plain by Siberian larch (Larix sibirica): a range‐wide analysis of cytoplasmic markers

While many species were confined to southern latitudes during the last glaciations, there has lately been mounting evidence that some of the most cold‐tolerant species were actually able to survive close to the ice sheets. The contribution of these higher latitude outposts to the main recolonization thrust remains, however, untested. In the present study, we use the first range‐wide survey of genetic diversity at cytoplasmic markers in Siberian larch (Larix sibirica; four mitochondrial (mt) DNA loci and five chloroplast (cp) DNA SSR loci) to (i) assess the relative contributions of southern and central areas to the current L. sibirica distribution range; and (ii) date the last major population expansion in both L. sibirica and adjacent Larix species. The geographic distribution of cpDNA variation was uninformative, but that of mitotypes clearly indicates that the southernmost populations, located in Mongolia and the Tien‐Shan and Sayan Mountain ranges, had a very limited contribution to the current populations of the central and northern parts of the range. It also suggests that the contribution of the high latitude cryptic refugia was geographically limited and that most of the current West Siberian Plain larch populations likely originated in the foothills of the Sayan Mountains. Interestingly, the main population expansion detected through Approximate Bayesian Computation (ABC) in all four larch species investigated here pre‐dates the LGM, with a mode in a range of 220 000–1 340 000 years BP. Hence, L. sibirica, like other major conifer species of the boreal forest, was strongly affected by climatic events pre‐dating the Last Glacial Maximum.     

Two Robin Populations

Capsule This study is the first ever documented evidence of an interglacial refugium during the Last Interglacial for birds in Anatolia and suggests the need of a re-examination of the effects of the Last Interglacial on the geographic distribution and genetic structure of species.

Aims We tested whether, in accordance with the ‘refugia within refugia’ model, multiple refugia existed for Kruper's Nuthatch Sitta krueperi during the Last Glacial Maximum or the species survived along the coastal belt of Anatolia through the Late Quaternary glacial–interglacial cycles.

Methods An ecological niche model was developed to predict the geographic distribution of Kruper's Nuthatch under reconstructed past (the Last Interglacial and the Last Glacial Maximum), present, and projected future bioclimatic conditions. Also, robust coalescent-based analyses were used to assess demographic events over the history of Kruper's Nuthatch.

Results Kruper's Nuthatch survived the Last Glacial Maximum almost along the coastal belt of Anatolia, but not in multiple refugia, and surprisingly, contrary to expectations, it survived the Last Interglacial in southern Anatolia, but not along the coastal belt of Anatolia.

Conclusion A kind of the ‘refugia within refugia’ model (i.e. the ‘refugium within refugium’ model) was supported because range shifts took place within Anatolia (itself also a refugium) for Kruper's Nuthatch.     

A new Late-glacial site with Picea abies in the northern Apennine foothills: an exception to the model of glacial refugia of trees

We describe a new palaeobotanical site at Bubano quarry on the easternmost Po plain, northern Italy. Pollen and macrofossils from river and marsh sediments demonstrate the occurrence of Picea in a Pinus sylvestris forest growing in a radius of some tens of kilometres south of the sedimentation place, at the beginning of the Late-glacial interstadial. The Late-glacial and Holocene history of Picea in the northern Apennines is reconstructed on the basis of the palaeobotanical record. The sharp climatic continentality increase eastwards across the northern Apennines from the Tyrrhenian to the Adriatic coast is considered significant for the survival of Picea during the Late-glacial. The most critical phase of survival is related to the moisture changes and consequent Abies competition associated with the last glacial-interglacial transition and the early Holocene. The residual spruce populations expanded during the middle Holocene. The history of Picea in the northern Apennines is a case of ineffective interglacial spread of tree populations from pre-existing stands of LGM (Last Glacial Maximum) and Late-glacial age.     

Do plant populations on distinct inselbergs talk to each other? A case study of genetic connectivity of a bromeliad species in an Ocbil landscape

Here, we explore the historical and contemporaneous patterns of connectivity among Encholirium horridum populations located on granitic inselbergs in an Ocbil landscape within the Brazilian Atlantic Forest, using both nuclear and chloroplast microsatellite markers. Beyond to assess the E. horridum population genetic structure, we built species distribution models across four periods (current conditions, mid‐Holocene, Last Glacial Maximum [LGM], and Last Interglacial) and inferred putative dispersal corridors using a least‐cost path analysis to elucidate biogeographic patterns. Overall, high and significant genetic divergence was estimated among populations for both nuclear and plastid DNA (Φ_ST(n) = 0.463 and Φ_ST(plastid) = 0.961, respectively, p < .001). For nuclear genome, almost total absence of genetic admixture among populations and very low migration rates were evident, corroborating with the very low estimates of immigration and emigration rates observed among E. horridum populations. Based on the cpDNA results, putative dispersal routes in Sugar Loaf Land across cycles of climatic fluctuations in the Quaternary period revealed that the populations’ connectivity changed little during those events. Genetic analyses highlighted the low genetic connectivity and long‐term persistence of populations, and the founder effect and genetic drift seemed to have been very important processes that shaped the current diversity and genetic structure observed in both genomes. The genetic singularity of each population clearly shows the need for in situ conservation of all of them.     

Quaternary climate change and habitat preference shaped the genetic differentiation and phylogeography of Rhodiola sect. Prainia in the southern Qinghai–Tibetan Plateau

There are two long‐standing biogeographic hypotheses regarding the glacial survival of plant species in the Qinghai–Tibetan Plateau (QTP): the in situ survival hypothesis and the tabula rasa hypothesis. We tested these two hypotheses in a phylogeographic study of Rhodiola sect. Prainia, a monophyletic section with ecologically divergent lineages. Molecular data from the nuclear internal transcribed spacer, six plastid markers and 13 nuclear microsatellite loci were analyzed for 240 individuals from 19 populations of this section. Environmental data were used to analyze the niches of major phylogenetic lineages within this section and to model changes in their distributions since the Last Glacial Maximum (LGM). We found that Rhodiola sect. Prainia consists of three evolutionary lineages: all populations of R. stapfii, R. prainii populations at the southern edge of the QTP, and R. prainii populations in the interior part of the QTP. During the LGM, the survival of R. prainii in the interior part of the QTP corresponded with the in situ survival hypothesis, while R. stapfii most probably survived the LGM in a manner corresponding with the tabula rasa hypothesis. The evolutionary history of different lineages of this section was shaped by topography, climate change, and lineage‐specific habitat preferences.     

East Asian allopatry and north Eurasian sympatry in Long‐tailed Tit lineages despite similar population dynamics during the late Pleistocene

Eurasia is a large continent characterized by heterogeneous environments. Glacial cycles during the late Pleistocene have had variable impacts on the avifauna across Eurasia. Bird populations from South‐East Asia show stability through the Last Glacial Maximum (LGM), while populations from Europe exhibit evidence of post‐LGM expansion. We investigated the phylogeography of the Long‐tailed Tit (Aegithalos caudatus), which spans the longitudinal breadth of Eurasia to test how climatic history and regional topographical complexity affected populations and diversification within the species complex. Our results show that two lineages from central and southern China (lineages C and D) segregate geographically, while lineages across northern Eurasia (lineage A and B) show substantial sympatry. Bayesian estimates for the timing of diversification suggest that the four lineages diverged during the middle Pleistocene, splitting in parallel and undergoing concurrent demographic histories since divergence. A. caudatus lineages experienced similar and synchronous population size dynamics during glacial cycles before the LGM. We conclude that the difference in geo‐topologic complexity may be an important factor that led to the variation in secondary admixture between northern Eurasian and eastern Asian lineages.     

The recolonization of Europe by brown bears Ursus arctos Linnaeus, 1758 after the Last Glacial Maximum

1. At the end of the Last Glacial Maximum brown bears Ursus arctos recolonized the glacial landscape of Central and Northern Europe faster than all other carnivorous mammal species of the Holocene fauna. Ursus arctos was recorded in Northern Europe from the beginning of the Late-Glacial. The recolonization of northern Central Europe may have taken place directly after the maximum glaciation. The distribution of the brown bear was restricted to glacial refugia only during the Last Glacial Maximum, for probably no more than 10 000 years. 2. Genetic analyses have suggested three glacial refugia for the brown bear: the Iberian Peninsula, the Italian Peninsula and the Balkans. Subfossil records of Ursus arctos from north-western Moldova as well as reconstructed environmental conditions during the Last Glacial Maximum in this area suggest to us a fourth glacial refuge for the brown bear. Because of its connection to the Carpathians, we designate this as the ‘Carpathian refuge’. 3. Due to the low genetic distance between brown bears of northern Norway, Finland, Estonia, north-eastern Russia and the northern Carpathians (the so-called eastern lineage), the Carpathians were considered the geographical origin of the recolonization of these regions. During the recolonization of northern Europe the brown bear probably reached these areas rapidly from the putative Carpathian refuge.     

Climate‐driven shifts in the distribution of koala‐browse species from the Last Interglacial to the near future

The koala's Phascolarctos cinereus distribution is currently restricted to eastern and south‐eastern Australia. However, fossil records dating from 70 ± 4 ka (ka = 10³ yr) from south‐western Australia and the Nullarbor Plain are evidence of subpopulation extinctions in the southwest at least after the Last Interglacial (~128–116 ka). We hypothesize that koala sub‐population extinctions resulted from the eastward retraction of the koala's main browse species in response to unsuitable climatic conditions. We further posit a general reduction in the distribution of main koala‐browse trees in the near future in response climate change. We modelled 60 koala‐browse species and constructed a set of correlative species distribution models for five time periods: Last Interglacial (~128–116 ka), Last Glacial Maximum (~23–19 ka), Mid‐Holocene (~7–5 ka), present (interpolations of observed data, representative of 1960–1990), and 2070. We based our projections on five hindcasts and one forecast of climatic variables extracted from WorldClim based on two general circulation models (considering the most pessimistic scenario of high greenhouse‐gas emissions) and topsoil clay fraction. We used 17 dates of koala fossil specimens identified as reliable from 70 (± 4) to 535 (± 49) ka, with the last appearance of koalas at 70 ka in the southwest. The main simulated koala‐browse species were at their greatest modelled extent of suitability during the Last Glacial Maximum, with the greatest loss of koala habitat occurring between the Mid‐Holocene and the present. We predict a similar habitat loss between the present and 2070. The spatial patterns of habitat change support our hypothesis that koala extinctions in the southwest, Nullarbor Plain and central South Australia resulted from the eastward retraction of the dominant koala‐browse species in response to long‐term climate changes. Future climate patterns will likely increase the extinction risk of koalas in their remaining eastern ranges.     

Phylogeography of North African Atlas cedar (Cedrus atlantica, Pinaceae): Combined molecular and fossil data reveal a complex Quaternary history

Northwest Africa is a major hotspot of plant biodiversity, but very little is known about the Quaternary range dynamics of plant species in this region. Here we investigate the range-wide population structure and phylogeography of Atlas cedar (Cedrus atlantica), an emblematic forest tree endemic to Morocco and Algeria. We genotyped 261 individuals from 11 populations using AFLP markers. Data were analyzed using both conventional F(ST)-based techniques and Bayesian clustering. Overall population differentiation was high (F(ST) = 0.25). Two major groups of populations were identified, one distributed through the Rif and Middle Atlas mountains in Morocco and the other through the Algerian Tell Atlas and Aurès mountains as well as the Middle Atlas. Combined molecular and fossil data indicate that C. atlantica survived the Last Glacial Maximum in at least three disjunct refugia along the coast of the Mediterranean Sea, whereas the Middle Atlas, today the core of the species range, has been colonized relatively recently (<10000 yr BP). The colonization history of individual populations has left clear imprints in their present-day diversity, which may vary greatly even between nearby stands. Our study illustrates how integrating different data sources and analytical approaches can help elucidate complex range dynamics that would otherwise remain undeciphered.     

Microsatellite DNA analyses of a highly disjunct New Zealand tree reveal strong differentiation and imply a formerly more continuous distribution

Although New Zealand is a biodiversity hotspot, there has been little genetic investigation of why so many of its threatened and uncommon plants have naturally disjunct distributions. We investigated the small tree Pseudopanax ferox (Araliaceae), which has a widespread but highly disjunct lowland distribution within New Zealand. Genotyping of nuclear microsatellites and a chloroplast locus revealed pronounced genetic differentiation and four principal genetic clusters. Our results indicate that the disjunct distribution is a product of vicariance rather than long‐distance dispersal. This highlights the need to preserve multiple populations when disjunct distributions are the result of vicariance, rather than focusing conservation efforts on a core area, in order to retain as much as possible of a species’ evolutionary legacy and potential. Additionally, based on our genetic findings and the ecology of P. ferox, we hypothesize that it was more continuously distributed during the drier (but not maximally colder) interstadials of glacial periods and/or on the fertile soils available immediately postglacial. We further hypothesize that P. ferox belongs to a suite of species of drought‐prone and/or fertile habitats whose distributions are actually restricted during warmer and wetter interglacial periods, despite being principally of the lowlands. Our genetic data for P. ferox are also the first consistent with the survival during the Last Glacial Maxima of a lowland tree at high latitudes in the south‐eastern South Island.