Phylogeography and genetic structure of the red-legged partridge (Alectoris rufa): more evidence for refugia within the Iberian glacial refugium

The Pleistocene climatic oscillations promoted the diversification in avian species during the last glacial period. The red‐legged partridge (Alectoris rufa, Family Phasianidae) has a large natural distribution extending from the Mediterranean to humid temperate zones. However, the genetic structure for this species is unknown. The present study investigates the phylogeography, genetic structure and demographic history of A. rufa across its distribution, employing both mitochondrial DNA control region sequences and nuclear microsatellite loci. Our results propose that this species was greatly affected by Pleistocene glaciations. The mismatch analyses suggest that the current populations resulted from post‐glacial expansion and subsequent differentiation resulting in five diagnosable genetic clusters: Southwestern, Central‐eastern, Northwestern, Balearic and French and Italian. Further, we found evidence of three glacial refugia within the currently recognized Iberian glacial refugium. The intraspecific structure revealed by both maternal and biparental phylogeographic analyses was not resolved in the phylogenetic analyses. Based on all considerations, we recommended that five management units be recognized.     

Phylogeography of the Northern Alligator Lizard (Squamata,Anguidae): Hidden diversity in a western endemic

Western North America includes the California Floristic Province and the Pacific Northwest, biologically diverse regions highlighted by a complex topography, geology, climate and history. A number of animals span these regions and show distinctive patterns of dispersal, vicariance and lineage diversification. Examining phylogeographic patterns in the fauna of this area aids in our understanding of the forces that have contributed to the generation and maintenance of regional biodiversity. Here, we investigate the biogeography and population structure of the Northern Alligator Lizard (Elgaria coerulea), a wide‐ranging anguid endemic to western North America. We sequenced two mtDNA fragments (ND2 and ND4) for 181 individuals across the range of the species and analysed these data with phylogenetic approaches to infer population and biogeographic history, and date major divergences within the taxon. We further used Bayesian clustering methods to assess major patterns of population structure and performed ecological niche modelling (ENM) to aid in our interpretation of geographic structure and diversification of E. coerulea lineages. Our phylogeographic examination of E. coerulea uncovered surprising diversity and structure, recovering 10 major lineages, each with substantial geographic substructure. While some divergences within the species are relatively old (Pliocene, 5.3–2.6 mya), most intraspecific variation appears to be of more recent origin (Pleistocene, 2.6 mya‐11,700 ya). Current diversity appears to have arisen in the Sierra Nevada Mountains and spread west and north since the Pliocene. Finally, our ENMs suggest that much of the Coast Ranges in California provided ideal habitat during the Last Glacial Maxima (LGM) that has since contracted dramatically and shifted northwards, whereas significant portions of the Sierra Nevada were unsuitable during the LGM and have since become more suitable. Interestingly, E. coerulea shares a number of genetic boundaries with other sympatric taxa, suggesting common historical events and geomorphological features have shaped the biota of this region.     

Phylogeography of Schizopygopsis stoliczkai (Cyprinidae) in Northwest Tibetan Plateau area

Schizopygopsis stoliczkai (Cyprinidae, subfamily Schizothoracinae) is one of the major freshwater fishes endemic to the northwestern margin of the Tibetan Plateau. In the current study, we used mitochondrial DNA markers cytochrome b (Cyt b) and 16S rRNA (16S), as well as the nuclear marker, the second intron of the nuclear beta‐actin gene (Act2), to uncover the phylogeography of S. stoliczkai. In total, we obtained 74 haplotypes from 403 mitochondrial concatenated sequences. The mtDNA markers depict the phylogenetic structures of S. stoliczkai, which consist of clade North and clade South. The split time of the two clades is dated back to 4.27 Mya (95% HPD = 1.96–8.20 Mya). The estimated split time is earlier than the beginning of the ice age of Pleistocene (2.60 Mya), suggesting that the northwestern area of the Tibetan Plateau probably contain at least two glacial refugia for S. stoliczkai. SAMOVA supports the formation of four groups: (i) the Karakash River group; (ii) The Lake Pangong group; (iii) the Shiquan River group; (iv) the Southern Basin group. Clade North included Karakash River, Lake Pangong, and Shiquan River groups, while seven populations of clade South share the haplotypes. Genetic diversity, star‐like network, BSP analysis, as well as negative neutrality tests indicate recent expansions events of S. stoliczkai. Conclusively, our results illustrate the phylogeography of S. stoliczkai, implying the Shiquan River is presumably the main refuge for S. stoliczkai.     

Phylogeography and palaeodistribution modelling in the Patagonian steppe: the case of Mulinum spinosum (Apiaceae)

Aim An integrative study of the endemic, yet ubiquitous, Patagonian shrub Mulinum spinosum (Apiaceae) was performed: (1) to assess the historical processes that influenced its geographical pattern of genetic variation; (2) to test hypotheses of its survival in situ or in glacial refugia during glacial cycles; and (3) to model its extant and palaeoclimatic distributions to assess support for the phylogeographical patterns recovered. Location Chilean and Argentinian Andean region and Patagonian steppe. Methods Chloroplast DNA sequences, trnH–psbA, trnS–trnG and 3′trnV–ndhC, were obtained for 314 individuals of M. spinosum from 71 populations. The haplotype data matrix was analysed using nested clade analysis (NCA) to construct a network. Analysis of molecular variance (AMOVA), spatial analysis of molecular variance (SAMOVA) and neutrality tests were also used to test for genetic structure and range expansion in the species. The present potential geographical distribution of M. spinosum was modelled and projected onto a Last Glacial Maximum (LGM) model. Results Amongst the 29 haplotypes observed, one was widely distributed, but most were restricted to either northern or southern regions. The populations with highest haplotype diversity were found in southern Patagonia, the high Andean region, and northern Patagonia. AMOVA and SAMOVA showed latitudinal structure for Argentinian populations. NCA implied patterns of restricted gene flow or dispersal but with some long‐distance dispersal and also long‐distance colonization and/or past fragmentation. Neutrality tests did not support range expansions. The current distribution model was a fairly good representation of the extant geographical distribution of the species, and the distribution model for the LGM did not show important shifts of the extant range to lower latitudes, except for a shift towards the palaeoseashore. Main conclusions Based on agreement amongst phylogeographical patterns, distribution of genetic variability, equivocal evidence of putative refugia and palaeodistribution modelling, it is probable that glaciations did not greatly affect the distribution of Mulinum spinosum. Our results are consistent with the in situ survival hypothesis, and not with the latitudinal migration of plant communities to avoid adverse climate conditions during Pleistocene glaciations. It is possible that populations of northern Patagonia may have been isolated from the southern ones by the Chubut and Deseado basins.     

Toward ecologically realistic predictions of species distributions: A cross‐time example from tropical montane cloud forests

There is an urgent need for more ecologically realistic models for better predicting the effects of climate change on species’ potential geographic distributions. Here we build ecological niche models using MAXENT and test whether selecting predictor variables based on biological knowledge and selecting ecologically realistic response curves can improve cross‐time distributional predictions. We also evaluate how the method chosen for extrapolation into nonanalog conditions affects the prediction. We do so by estimating the potential distribution of a montane shrew (Mammalia, Soricidae, Cryptotis mexicanus) at present and the Last Glacial Maximum (LGM). Because it is tightly associated with cloud forests (with climatically determined upper and lower limits) whose distributional shifts are well characterized, this species provides clear expectations of plausible vs. implausible results. Response curves for the MAXENT model made using variables selected via biological justification were ecologically more realistic compared with those of the model made using many potential predictors. This strategy also led to much more plausible geographic predictions for upper and lower elevational limits of the species both for the present and during the LGM. By inspecting the modeled response curves, we also determined the most appropriate way to extrapolate into nonanalog environments, a previously overlooked factor in studies involving model transfer. This study provides intuitive context for recommendations that should promote more realistic ecological niche models for transfer across space and time.     

Phylogeography of a species complex of lowland Neotropical rain forest trees (Carapa,Meliaceae)

Aim Many tropical tree species have poorly delimited taxonomic boundaries and contain undescribed or cryptic species. We examined the genetic structure of a species complex in the tree genus Carapa in the Neotropics in order to evaluate age, geographic patterns of diversity and evolutionary relationships, and to quantify levels of introgression among currently recognized species. Location Lowland moist forests in the Guiana Shield, the Central and Western Amazon Basin, Chocó and Central America. Methods Genetic structure was analysed using seven nuclear simple sequence repeats (nuSSR), five chloroplast SSRs (cpSSR), and two chloroplast DNA (cpDNA) intergenic sequences (trnH–psbA and trnC–ycf6). Bayesian clustering analysis of the SSR data was used to infer population genetic structure and to assign 324 samples to their most likely genetic cluster. Bayesian coalescence analyses were performed on the two cpDNA markers to estimate evolutionary relationships and divergence times. Results Two genetic clusters (nu_guianensis and nu_surinamensis) were detected, which correspond to the Neotropical species C. guianensis (sensu latu) and C. surinamensis. Fourteen cpDNA haplotypes clustered into six haplogroups distributed between the two nuclear genetic clusters. Divergence between the haplogroups was initiated in the Miocene, with some haplotype structure evolving as recently as the Pleistocene. The absence of complete lineage sorting between the nuclear and chloroplast genomes and the presence of hybrid individuals suggest that interspecific reproductive barriers are incomplete. NuSSR diversity was highest in C. guianensis and, within C. guianensis, cpDNA diversity was highest in the Central and Western Amazon Basin. Regional genetic differentiation was strong but did not conform to an isolation‐by‐distance process or exhibit a phylogeographical signal. Main conclusions The biogeographical history of Neotropical Carapa appears to have been influenced by events that took place during the Neogene. Our results point to an Amazonian centre of origin and diversification of Neotropical Carapa, with subsequent migration to the Pacific coast of South America and Central America. Gene flow apparently occurs among species, and introgression events are supported by inconsistencies between chloroplast and nuclear lineage sorting. The absence of phylogeographical structure may be a result of the ineffectiveness of geographical barriers among populations and of reproductive isolation mechanisms among incipient and cryptic species in this species complex.     

Phylogeography and niche modelling of the relict plant Amborella trichopoda (Amborellaceae) reveal multiple Pleistocene refugia in New Caledonia

Amborella trichopoda Baill. (Amborellaceae, Amborellales), the sole living member of the sister group to all other extant angiosperms, is endemic to New Caledonia. We addressed the intraspecific phylogeography of Amborella by investigating whether its present population genetic structure could be related to its current and past habitats. We found moderate range‐wide genetic diversity based on nuclear microsatellite data and detected four well‐differentiated, geographically distinct genetic groups using Bayesian clustering analyses. We modelled the ecological niche of Amborella based on the current climatic and environmental conditions. The predictive ability of the model was very good throughout the Central East mainland zone, but Amborella was predicted in the northern part of the island where this plant has not been reported. Furthermore, no significant barrier was detected based on habitat suitability that could explain the genetic differentiation across the area. Conversely, we found that the main genetic clusters could be related to the distribution of the suitable habitat at the last glacial maximum (LGM, c. 21 000 years BP), when Amborella experienced a dramatic 96.5% reduction in suitable area. At least two lineages survived in distinct putative refugia located in the Massif des Lèvres and in the vicinity of Mount Aoupinié. Our findings finally confirmed the importance of LGM rainforest refugia in shaping the current intra‐ and interspecific diversity in New Caledonian plants and revealed the possibility of an as yet unreported refugium. The combination of niche modelling and population genetics thereby offered novel insight into the biogeographical history of an emblematic taxon.     

Phylogeography and species delineation of the genus Phoxinus Rafinesque, 1820 (Actinopterygii: Leuciscidae) in the Iberian Peninsula

The genus Phoxinus is comprised of at least 15 currently recognized species inhabiting Eurasia. Morphological traits have been traditionally used to delineate species in Phoxinus; however, the high level of phenotypic plasticity observed in the genus has confounded this process. Molecular genetic analyses have revealed a higher than expected genetic structure within Phoxinus. Here, we analyzed both nuclear and mitochondrial molecular genetic markers to infer the phylogeography and divergence times of Phoxinus in the Iberian Peninsula. Our results show that the Iberian lineages of Phoxinus were polyphyletic. They also support the co‐existence of three species in the Iberian Peninsula, two corresponding to two previously recognized species (Phoxinus bigerri and Phoxinus septimaniae) and a third undescribed species (Phoxinus sp.). Phoxinus bigerri is structured into western Cantabrian, eastern Cantabrian, and Artibai basins. We hypothesize that this structure is a consequence of glaciation–deglaciation cycles during the Pleistocene. While the presence of P. septimaniae in the Iberian Peninsula is possibly the result of human translocation, that of Phoxinus sp. in lower Ebro rivers may be attributed to past fluvial captures. Our study represents the first report to show a relationship among Phoxinus populations from central Pyrenean rivers of Spain and France. Furthermore, we found genetic hybridization between Phoxinus sp. and P. septimaniae in the shared localities, a likely consequence of anthropogenic activities. Overall, our findings provide insight into the genetic structure of Iberian Phoxinus populations, including the presence of an undescribed species and the putative introduction of some species that may have implications for conservation.