In the light of evolution X: Comparative phylogeography

A phylogeographic approach of three species of Clinostomum Leidy, 1856, (Trematoda: Clinostomidae) from the Neotropical region of Mexico, with the description of a new species from Ardea herodias (Ardeidae)

Clinostomum Leidy, 1856 is a globally distributed group of endoparasites that reside attached to the mucosa of the mouth cavity and oesophagus of fish-eating birds (mainly herons) in the adult stage, whereas larvae known as metacercariae ("yellow grubs") damage freshwater fish and amphibians that serve as secondary intermediate hosts. In the Nearctic and Neotropical regions of Mexico, species of the genus Clinostomum have been intensively studied, revealing seven species that share a common ancestor. In the present study, sequences of the ITS from nuclear ribosomal DNA and cox 1 from mitochondrial DNA were generated for C. tataxumui, C. cichlidorum, and an unidentified species identified as lineage 3 sensu Pérez-Ponce de León et al. 2026, which is distributed sympatrically in Mexico and Middle America. The objectives of the present study were to 1) describe a new species of Clinostomum, named herein as Clinostomum chaacci n. sp., which combines morphological and molecular characters; 2) analyse the population genetic structure of two previously described species in addition to C. chaacci n. sp.; and 3) infer phylogeographic patterns using cox 1 sequences. The inferred phylogeographic results indicated that the populations of C. tataxumui, and C. cichlidorum lack a geographic structure with high haplotype diversity, low nucleotide diversity, and low Fst values among the populations with negative values on the neutrality test, suggesting that the populations are expanding. In contrast, the populations of C. chaacci n. sp., are less connected and subdivided into smaller populations, leading to slight genetic differentiation, suggesting that the populations have undergone significant demographic changes, such as recent expansion or selection effects.

Artificial intelligence enables unified analysis of historical and landscape influences on genetic diversity

While genetic variation in any species is potentially shaped by a range of processes, phylogeography and landscape genetics are largely concerned with inferring how environmental conditions and landscape features impact neutral intraspecific diversity. However, even as both disciplines have come to utilize SNP data over the last decades, analytical approaches have remained for the most part focused on either broad-scale inferences of historical processes (phylogeography) or on more localized inferences about environmental and/or landscape features (landscape genetics). Here we demonstrate that an artificial intelligence model-based analytical framework can consider both deeper historical factors and landscape-level processes in an integrated analysis. We implement this framework using data collected from two Brazilian anurans, the Brazilian sibilator frog (Leptodactylus troglodytes) and granular toad (Rhinella granulosa). Our results indicate that historical demographic processes shape most the genetic variation in the sibulator frog, while landscape processes primarily influence variation in the granular toad. The machine learning framework used here allows both historical and landscape processes to be considered equally, rather than requiring researchers to make an a priori decision about which factors are important.

Phylogenetic assessment within a species complex of a subterranean rodent (Geomys bursarius) with conservation implications for isolated subspecies

Range contraction and expansion from glaciation have led to genetic divergence that may be particularly pronounced in fossorial species with low dispersal. The plains pocket gopher (Geomys bursarius) is a fossorial species that ranges widely across North America but has a poorly understood phylogeny. We used mitogenomes (14,996 base pairs) from 56 individuals across seven subspecies, plus two outgroup species, to assess genetic divergence from minimum spanning trees, measure genetic distances, and infer phylogenetic trees using BEAST. We found G. b. wisconsinensis was monophyletic with recent divergence. Further assessment is needed for G. b. major because it was paraphyletic and exhibited inconsistent groupings with other clades. Importantly, we identified G. b. illinoensis as being genetically distinct and monophyletic likely due to a unique colonization event eastward across the Mississippi River. Because G. b. illinoensis faces continued pressures from niche reduction and habitat loss, we recommend that G. b. illinoensis be considered an evolutionary significant unit warranting conservation actions to promote connectivity and restore suitable habitat. Such conservation efforts should benefit other grassland species including those originating from clades west of the Mississippi River that may also be evolutionary significant units.

Phylogeography of the desert scorpion illuminates a route out of Central Asia

A comprehensive understanding of phylogeography requires the integration of knowledge across different organisms, ecosystems, and geographic regions. However, a critical knowledge gap exists in the arid biota of the vast Asian drylands. To narrow this gap, here we test an "out-of-Central Asia" hypothesis for the desert scorpion Mesobuthus mongolicus by combining Bayesian phylogeographic reconstruction and ecological niche modeling. Phylogenetic analyses of one mitochondrial and three nuclear loci and molecular dating revealed that M. mongolicus represents a coherent lineage that diverged from its most closely related lineage in Central Asia about 1.36 Ma and underwent radiation ever since. Bayesian phylogeographic reconstruction indicated that the ancestral population dispersed from Central Asia gradually eastward to the Gobi region via the Junggar Basin, suggesting that the Junggar Basin has severed as a corridor for Quaternary faunal exchange between Central Asia and East Asia. Two major dispersal events occurred probably during interglacial periods (around 0.8 and 0.4 Ma, respectively) when climatic conditions were analogous to present-day status, under which the scorpion achieved its maximum distributional range. M. mongolicus underwent demographic expansion during the Last Glacial Maximum, although the predicted distributional areas were smaller than those at present and during the Last Interglacial. Development of desert ecosystems in northwest China incurred by intensified aridification might have opened up empty habitats that sustained population expansion. Our results extend the spatiotemporal dimensions of trans-Eurasia faunal exchange and suggest that species' adaptation is an important determinant of their phylogeographic and demographic responses to climate changes.

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large‐scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD‐sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species‐specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

Plant migration under long-lasting hyperaridity - phylogenomics unravels recent biogeographic history in one of the oldest deserts on Earth

The post-Miocene climatic histories of arid environments have been identified as key drivers of dispersal and diversification. Here, we investigate how climatic history correlates with the historical biogeography of the Atacama Desert genus Cristaria (Malvaceae). We analyze phylogenetic relationships and historical biogeography by using next-generation sequencing (NGS), molecular clock dating, Dispersal Extinction Cladogenesis and Bayesian sampling approaches. We employ a novel way to identify biogeographically meaningful regions as well as a rarely utilized program permitting the use of dozens of ancestral areas. Partial incongruence between the established taxonomy and our phylogenetic data argue for a complex historical biogeography with repeated introgression and incomplete lineage sorting. Cristaria originated in the central southern part of the Atacama Desert, from there the genus colonized other areas from the late Miocene onwards. The more recently diverged lineages appear to have colonized different habitats in the Atacama Desert during pluvial phases of the Pliocene and early Pleistocene. We show that NGS combined with near-comprehensive sampling can provide an unprecedented degree of phylogenetic resolution and help to correlate the historical biogeography of plant communities with cycles of arid and pluvial phases.

Climate Warming Since the Holocene Accelerates West-East Communication for the Eurasian Temperate Water Strider Species Aquarius paludum

Holocene climate warming has dramatically altered biological diversity and distributions. Recent human-induced emissions of greenhouse gases will exacerbate global warming and thus induce threats to cold-adapted taxa. However, the impacts of this major climate change on transcontinental temperate species are still poorly understood. Here, we generated extensive genomic datasets for a water strider, Aquarius paludum, which was sampled across its entire distribution in Eurasia and used these datasets in combination with ecological niche modeling (ENM) to elucidate the influence of the Holocene and future climate warming on its population structure and demographic history. We found that A. paludum consisted of two phylogeographic lineages that diverged in the middle Pleistocene, which resulted in a “west–east component” genetic pattern that was probably triggered by Central Asia-Mongoxin aridification and Pleistocene glaciations. The diverged western and eastern lineages had a second contact in the Holocene, which shaped a temporary hybrid zone located at the boundary of the arid–semiarid regions of China. Future predictions detected a potentially novel northern corridor to connect the western and eastern populations, indicating west–east gene flow would possibly continue to intensify under future warming climate conditions. Further integrating phylogeographic and ENM analyses of multiple Eurasian temperate taxa based on published studies reinforced our findings on the “west–east component” genetic pattern and the predicted future northern corridor for A. paludum. Our study provided a detailed paradigm from a phylogeographic perspective of how transcontinental temperate species differ from cold-adapted taxa in their response to climate warming.

Population genomics for symbiotic anthozoans: can reduced representation approaches be used for taxa without reference genomes?

Population genetic studies of symbiotic anthozoans have been historically challenging because their endosymbioses with dinoflagellates have impeded marker development. Genomic approaches like reduced representation sequencing alleviate marker development issues but produce anonymous loci, and without a reference genome, it is unknown which organism is contributing to the observed patterns. Alternative methods such as bait-capture sequencing targeting Ultra-Conserved Elements are now possible but costly. Thus, RADseq remains attractive, but how useful are these methods for symbiotic anthozoan taxa without a reference genome to separate anthozoan from algal sequences? We explore this through a case-study using a double-digest RADseq dataset for the sea anemone Bartholomea annulata. We assembled a holobiont dataset (3854 loci) for 101 individuals, then used a reference genome to create an aposymbiotic dataset (1402 loci). For both datasets, we investigated population structure and used coalescent simulations to estimate demography and population parameters. We demonstrate complete overlap in the spatial patterns of genetic diversity, demographic histories, and population parameter estimates for holobiont and aposymbiotic datasets. We hypothesize that the unique combination of anthozoan biology, diversity of the endosymbionts, and the manner in which assembly programs identify orthologous loci alleviates the need for reference genomes in some circumstances. We explore this hypothesis by assembling an additional 21 datasets using the assembly programs pyRAD and Stacks. We conclude that RADseq methods are more tractable for symbiotic anthozoans without reference genomes than previously realized.

Variation in intraspecific demography drives localised concordance but species-wide discordance in response to past climatic change

BACKGROUND

Understanding how species biology may facilitate resilience to climate change remains a critical factor in detecting and protecting species at risk of extinction. Many studies have focused on the role of particular ecological traits in driving species responses, but less so on demographic history and levels of standing genetic variation. Additionally, spatial variation in the interaction of demographic and adaptive factors may further complicate prediction of species responses to environmental change. We used environmental and genomic datasets to reconstruct the phylogeographic histories of two ecologically similar and largely co-distributed freshwater fishes, the southern (Nannoperca australis) and Yarra (N. obscura) pygmy perches, to assess the degree of concordance in their responses to Plio-Pleistocene climatic changes. We described contemporary genetic diversity, phylogenetic histories, demographic histories, and historical species distributions across both species, and statistically evaluated the degree of concordance in co-occurring populations.

RESULTS

Marked differences in contemporary genetic diversity, historical distribution changes and historical migration were observed across the species, with a distinct lack of genetic diversity and historical range expansion suggested for N. obscura. Although several co-occurring populations within a shared climatic refugium demonstrated concordant demographic histories, idiosyncratic population size changes were found at the range edges of the more spatially restricted species. Discordant responses between species were associated with low standing genetic variation in peripheral populations. This might have hindered adaptive potential, as documented in recent demographic declines and population extinctions for the two species.

CONCLUSION

Our results highlight both the role of spatial scale in the degree of concordance in species responses to climate change, and the importance of standing genetic variation in facilitating range shifts. Even when ecological traits are similar between species, long-term genetic diversity and historical population demography may lead to discordant responses to ongoing and future climate change.

The usefulness of maternally inherited genetic markers for phylogeographic studies in village chicken

Phylogeography plays a major role in understanding micro and macroevolutionary processes dealing with evolutionary interpretations of geographical distribution. This field integrates information from molecular genetics, population genetics, demography, and phylogeny for the interpretation of the geographical distribution of lineages. The full mtDNA sequence and W chromosome polymorphisms were exploited to assess the usefulness of two maternally-inherited genetic markers for phylogeographic studies of village chickens. We studied 243 full mtDNA sequences from three countries (Iraq, n = 27; Ethiopia, n = 211; and Saudi Arabia, n = 5) and a 13-kb fragment of the W chromosome from 20 Iraqi and 137 Ethiopian female chickens. The results show a high level of genetic diversity for the mtDNA within and among countries as well as within populations. On the other hand, sequence analysis of the W chromosome shows low genetic diversity both within and among populations. Six full mtDNA haplogroups (A, B, C1, C2, D1, and E1) were observed and 25 distinct W haplotypes. The results support the effectiveness of full mtDNA sequences but not the W chromosome in tracing the maternal historical genome background with, however, weak within a country phylogeographic signal.

Phylogeography of the white-crowned parrot (Pionus senilis)

Abstract The white-crowned parrot Pionus senilis (von Spix, 1824) is distributed throughout Middle America, inhabiting the Gulf of Mexico coastal area from Tamaulipas (Mexico) to northern Panama. We used mitochondrial data (COI, ND2 and ND4) from 55 specimens to infer phylogenetic relationships, and analyzed the phylogeographic structure, genetic diversity, divergence periods, and historical demography to explore phylogeographic patterns. We found three divergent lineages: two geographically separated by the Isthmus of Tehuantepec, and the third, in Costa Rica by the Nicaragua Depression. The analysis of molecular variance and statistical analyses were consistent with genetically distinct populations. The Central American lineage diverged 1.33 million years ago, whereas the other two lines branched off 1.19 million years ago. This phylogenetic pattern has been reported in other species of Middle American birds.

The evolution of comparative phylogeography: putting the geography (and more) into comparative population genomics

Comparative population genomics is an ascendant field using genomic comparisons between species to draw inferences about forces regulating genetic variation. Comparative phylogeography, by contrast, focuses on the shared lineage histories of species codistributed geographically and is decidedly organismal in perspective. Comparative phylogeography is approximately 35 years old, and, by some metrics, is showing signs of reduced growth. Here, we contrast the goals and methods of comparative population genomics and comparative phylogeography and argue that comparative phylogeography offers an important perspective on evolutionary history that succeeds in integrating genomics with landscape evolution in ways that complement the suprageographic perspective of comparative population genomics. Focusing primarily on terrestrial vertebrates, we review the history of comparative phylogeography, its milestones and ongoing conceptual innovations, its increasingly global focus, and its status as a bridge between landscape genomics and the process of speciation. We also argue that, as a science with a strong “sense of place,” comparative phylogeography offers abundant “place-based” educational opportunities with its focus on geography and natural history, as well as opportunities for collaboration with local communities and indigenous peoples. Although comparative phylogeography does not yet require whole-genome sequencing for many of its goals, we conclude that it nonetheless plays an important role in grounding our interpretation of genetic variation in the fundamentals of geography and Earth history.

Haplotype network branch diversity, a new metric combining genetic and topological diversity to compare the complexity of haplotype networks

A common way of illustrating phylogeographic results is through the use of haplotype networks. While these networks help to visualize relationships between individuals, populations, and species, evolutionary studies often only quantitatively analyze genetic diversity among haplotypes and ignore other network properties. Here, we present a new metric, haplotype network branch diversity (HBd), as an easy way to quantifiably compare haplotype network complexity. Our metric builds off the logic of combining genetic and topological diversity to estimate complexity previously used by the published metric haplotype network diversity (HNd). However, unlike HNd which uses a combination of network features to produce complexity values that cannot be defined in probabilistic terms, thereby obscuring the values' implication for a sampled population, HBd uses frequencies of haplotype classes to incorporate topological information of networks, keeping the focus on the population and providing easy-to-interpret probabilistic values for randomly sampled individuals. The goal of this study is to introduce this more intuitive metric and provide an R script that allows researchers to calculate diversity and complexity indices from haplotype networks. A group of datasets, generated manually (model dataset) and based on published data (empirical dataset), were used to illustrate the behavior of HBd and both of its terms, haplotype diversity, and a new index called branch diversity. Results followed a predicted trend in both model and empirical datasets, from low metric values in simple networks to high values in complex networks. In short, the new combined metric joins genetic and topological diversity of haplotype networks, into a single complexity value. Based on our analysis, we recommend the use of HBd, as it makes direct comparisons of network complexity straightforward and provides probabilistic values that can readily discriminate situations that are difficult to resolve with available metrics.

A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present-day populations. In Europe, where both Mediterranean and extra-Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra-Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range-wide colonization. The results provide clear evidence that both extra-Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid-latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%-47%) and the Carpathian (53%-72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end-glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end-glacial colonization history of this well-studied organism and underscores the importance of genomic data in phylogeographic interpretation.

Common patterns in the molecular phylogeography of western palearctic birds: a comprehensive review

A plethora of studies have offered crucial insights in the phylogeographic status of Western Palearctic bird species. However, an overview integrating all this information and analyzing the combined results is still missing. In this study, we compiled all published peer-reviewed and grey literature available on the phylogeography of Western Palearctic bird species. Our literature review indicates a total number of 198 studies, with the overwhelming majority published as journal articles (n = 186). In total, these literature items offer information on 145 bird species. 85 of these species are characterized by low genetic differentiation, 46 species indicate genetic variation but no geographic structuring i.e. panmixia, while 14 species show geographically distinct lineages and haplotypes. Majority of bird species inhabiting the Western Palearctic display genetic admixture. The glaciation cycles in the past few million years were pivotal factors in shaping this situation: during warm periods many species expanded their distribution range to the north over wide areas of Eurasia; whereas, during ice ages most areas were no longer suitable and species retreated to refugia, where lineages mixed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01893-x.

Life history and ecology might explain incongruent population structure in two co-distributed montane bird species of the Atlantic Forest

Comparative phylogeography is a powerful approach to investigate the role of historical and environmental processes in the evolution of biodiversity within a region. In this regard, comparative studies of species with similar habitat preferences are valuable to reduce the confounding influence of habitat association when interpreting phylogeographic patterns. In the Atlantic Forest of South America, phylogeographic studies of highland and lowland species have shown distinct population structure patterns so far, suggesting that such species have responded differently to Pleistocene glacial cycles. Herein, we performed a comparative analysis using molecular data and paleodistribution models of two Montane Atlantic Forest (MAF) co-distributed passerine birds with similar habitat requirements but with distinct life-history traits and ecologies: the frugivore lek-breeding Blue Manakin (Chiroxiphia caudata) and the insectivore and socially monogamous Drab-Breasted Bamboo Tyrant (Hemitriccus diops). We aimed to shed light on the role of contrasting life histories and ecologies onto the demography and population structure of MAF species. We sampled both species throughout most of their distribution range, sequenced a mitochondrial and a nuclear molecular marker, and used standard phylogeographic methods to investigate population structure and ecological niche modeling (ENM) to infer the species' paleodistributions. Our analyses recovered a phylogeographic break in H. diops in the region of the Doce River, but no genetic structure in C. caudata. We also found higher differentiation among subpopulations within each lineage of H. diops than among subpopulations of C. caudata. We suggest that these discrepancies in population structure might be due to distinct life-history traits and their impact on gene flow and generation time. For example, while H. diops is an insectivore species, C. caudata is a frugivore and the latter ecological aspect likely selects for a higher dispersion distance. Additionally, because C. caudata is a lek-breeding species, it has a longer generation time than H. diops. These traits could hinder genetic differentiation when populations become geographically isolated. Nonetheless, both species showed some common biological features, such as signatures of synchronous population expansion and larger distribution ranges during the Last Glacial Maximum, possibly due to similar cold tolerance.

Circumpolar diversification of the Ixodes uriae tick virome

Ticks (order: Ixodida) are a highly diverse and ecologically important group of ectoparasitic blood-feeding organisms. One such species, the seabird tick (Ixodes uriae), is widely distributed around the circumpolar regions of the northern and southern hemispheres. It has been suggested that Ix. uriae spread from the southern to the northern circumpolar region millions of years ago and has remained isolated in these regions ever since. Such a profound biographic subdivision provides a unique opportunity to determine whether viruses associated with ticks exhibit the same evolutionary patterns as their hosts. To test this, we collected Ix. uriae specimens near a Gentoo penguin (Pygoscelis papua) colony at Neko harbour, Antarctica, and from migratory birds—the Razorbill (Alca torda) and the Common murre (Uria aalge)—on Bonden island, northern Sweden. Through meta-transcriptomic next-generation sequencing we identified 16 RNA viruses, seven of which were novel. Notably, we detected the same species, Ronne virus, and two closely related species, Bonden virus and Piguzov virus, in both hemispheres indicating that there have been at least two cross-circumpolar dispersal events. Similarly, we identified viruses discovered previously in other locations several decades ago, including Gadgets Gully virus, Taggert virus and Okhotskiy virus. By identifying the same or closely related viruses in geographically disjunct sampling locations we provide evidence for virus dispersal within and between the circumpolar regions. In marked contrast, our phylogenetic analysis revealed no movement of the Ix. uriae tick hosts between the same locations. Combined, these data suggest that migratory birds are responsible for the movement of viruses at both local and global scales.

As host populations diverge, so may those microorganisms, including viruses, that are dependent on those hosts. To examine this key issue in host-microbe evolution we compared the co-phylogenies of the seabird tick, Ixodes uriae, and their RNA viruses sampled from the far northern and southern hemispheres. Despite the huge geographic distance between them, phylogeographic analysis reveals that the same and closely related viruses were found both within and between the northern and southern circumpolar regions, most likely reflecting transfer by virus-infected migratory birds. In contrast, genomic data suggested that the Ix. uriae populations were phylogenetically distinct between the northern and southern hemispheres. This work emphasises the importance of migratory birds and ticks as vectors and sources of virus dispersal and introduction at both the local and global scales.

The Quaternary evolutionary history of Bristol rock cress (Arabis scabra, Brassicaceae), a Mediterranean element with an outpost in the north-western Atlantic region

BACKGROUND AND AIMS

Bristol rock cress is among the few plant species in the British Isles considered to have a Mediterranean-montane element. Spatiotemporal patterns of colonization of the British Isles since the last interglacial and after the Last Glacial Maximum (LGM) from mainland Europe are underexplored and have not yet included such floristic elements. Here we shed light on the evolutionary history of a relic and outpost metapopulation of Bristol rock cress in the south-western UK.

METHODS

Amplified fragment length polymorphisms (AFLPs) were used to identify distinct gene pools. Plastome assembly and respective phylogenetic analysis revealed the temporal context. Herbarium material was largely used to exemplify the value of collections to obtain a representative sampling covering the entire distribution range.

KEY RESULTS

The AFLPs recognized two distinct gene pools, with the Iberian Peninsula as the primary centre of genetic diversity and the origin of lineages expanding before and after the LGM towards mountain areas in France and Switzerland. No present-day lineages are older than 51 ky, which is in sharp contrast to the species stem group age of nearly 2 My, indicating severe extinction and bottlenecks throughout the Pleistocene. The British Isles were colonized after the LGM and feature high genetic diversity.

CONCLUSIONS

The short-lived perennial herb Arabis scabra, which is restricted to limestone, has expanded its distribution range after the LGM, following corridors within an open landscape, and may have reached the British Isles via the desiccated Celtic Sea at about 16 kya. This study may shed light on the origin of other rare and peculiar species co-occurring in limestone regions in the south-western British Isles.

Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome-wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans-Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter-riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.

Phylogeography of Ticks (Acari: Ixodida)

Improved understanding of tick phylogeny has allowed testing of some biogeographical patterns. On the basis of both literature data and a meta-analysis of available sequence data, there is strong support for a Gondwanan origin of Ixodidae, and probably Ixodida. A particularly strong pattern is observed for the genus Amblyomma, which appears to have originated in Antarctica/southern South America, with subsequent dispersal to Australia. The endemic Australian lineages of Ixodidae (no other continent has such a pattern) appear to result from separate dispersal events, probably from Antarctica. Minimum ages for a number of divergences are determined as part of an updated temporal framework for tick evolution. Alternative hypotheses for tick evolution, such as a very old Pangean group, a Northern hemisphere origin, or an Australian origin, fit less well with observed phylogeographic patterns.

Genotyping by sequencing reveals contrasting patterns of population structure, ecologically mediated divergence, and long-distance dispersal in North American palms

Comparative studies can provide powerful insights into processes that affect population divergence and thereby help to elucidate the mechanisms by which contemporary populations may respond to environmental change. Furthermore, approaches such as genotyping by sequencing (GBS) provide unprecedented power for resolving genetic differences among species and populations. We therefore used GBS to provide a genomewide perspective on the comparative population structure of two palm genera, Washingtonia and Brahea, on the Baja California peninsula, a region of high landscape and ecological complexity. First, we used phylogenetic analysis to address taxonomic uncertainties among five currently recognized species. We resolved three main clades, the first corresponding to W. robusta and W. filifera, the second to B. brandegeei and B. armata, and the third to B. edulis from Guadalupe Island. Focusing on the first two clades, we then delved deeper by investigating the underlying population structure. Striking differences were found, with GBS uncovering four distinct Washingtonia populations and identifying a suite of loci associated with temperature, consistent with ecologically mediated divergence. By contrast, individual mountain ranges could be resolved in Brahea and few loci were associated with environmental variables, implying a more prominent role of neutral divergence. Finally, evidence was found for long-distance dispersal events in Washingtonia but not Brahea, in line with knowledge of the dispersal mechanisms of these palms including the possibility of human-mediated dispersal. Overall, our study demonstrates the power of GBS together with a comparative approach to elucidate markedly different patterns of genomewide divergence mediated by multiple effectors.

Climate Is Not All: Evidence From Phylogeography of Rhodiola fastigiata (Crassulaceae) and Comparison to Its Closest Relatives

How geological events and climate oscillations in the Pleistocene glaciation shaped the geographic distribution of genetic variation of species on the Qinghai-Tibetan Plateau (QTP) and its adjacent areas has been extensively studied. However, little studies have investigated whether closely related species in the same genus with similar physiological and life history traits responded similarly to the glacial climatic oscillations. If this is not the case, we would expect that the population histories of studied species were not driven by extrinsic environmental changes alone. Here we conducted a phylogeographic study of a succulent alpine plant Rhodiola fastigiata, using sequences from chloroplast genome and nrITS region, as well as ecological niche modeling. The results of R. fastigiata were compared to other congeneric species that have been studied, especially to R. alsia and R. crenulata. We found that for both markers, two geographic groups could be revealed, corresponding to the QTP plateau and the Hengduan Mountains, respectively, indicating isolated refugia in those two areas. The two groups diverged 1.23 Mya during the Pleistocene. We detected no significant population expansion by mismatch distribution analysis and Bayesian Skyline Plot. We found that even these similar species with similar physiological and life history traits have had different demographic histories in the Quaternary glacial periods. Our comparative phylogeographic study sheds new lights into phylogeographic research that extrinsic environmental changes are not the only factor that can drive population demography, and other factors, such as coevolved interactions between plants and their specialized pathogens, that probably played a role need to be examined with more case studies.

Demographically idiosyncratic responses to climate change and rapid Pleistocene diversification of the walnut genus Juglans (Juglandaceae) revealed by whole-genome sequences

Whether species demography and diversification are driven primarily by extrinsic environmental changes such as climatic oscillations in the Quaternary or by intrinsic biological interactions like coevolution between antagonists is a matter of active debate. In fact, their relative importance can be assessed by tracking past population fluctuations over considerable time periods. We applied the pairwise sequentially Markovian coalescent approach on the genomes of 11 temperate Juglans species to estimate trajectories of changes in effective population size (N_e ) and used a Bayesian-coalescent based approach that simultaneously considers multiple genomes (G-PhoCS) to estimate divergence times between lineages. N_e curves of all study species converged 1.0 million yr ago, probably reflecting the time when the walnut genus last shared a common ancestor. This estimate was confirmed by the G-PhoCS estimates of divergence times. But all species did not react similarly to the dramatic climatic oscillations following early Pleistocene cooling, so the timing and amplitude of changes in N_e differed among species and even among conspecific lineages. The population histories of temperate walnut species were not driven by extrinsic environmental changes alone, and a key role was probably played by species-specific factors such as coevolutionary interactions with specialized pathogens.