Transverse Alpine Speciation Driven by Glaciation

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Comparative mitogenomic analysis provides evolutionary insights into Formica (Hymenoptera: Formicidae)

Formica is a large genus in the family Formicidae with high diversity in its distribution, morphology, and physiology. To better understand evolutionary characteristics of Formica, the complete mitochondrial genomes (mitogenomes) of two Formica species were determined and a comparative mitogenomic analysis for this genus was performed. The two newly sequenced Formica mitogenomes each included 37 typical mitochondrial genes and a large non-coding region (putative control region), as observed in other Formica mitogenomes. Base composition, gene order, codon usage, and tRNA secondary structure were well conserved among Formica species, whereas diversity in sequence size and structural characteristics was observed in control regions. We also observed several conserved motifs in the intergenic spacer regions. These conserved genomic features may be related to mitochondrial function and their highly conserved physiological constraints, while the diversity of the control regions may be associated with adaptive evolution among heterogenous habitats. A negative AT-skew value on the majority chain was presented in each of Formica mitogenomes, indicating a reversal of strand asymmetry in base composition. Strong codon usage bias was observed in Formica mitogenomes, which was predominantly determined by nucleotide composition. All 13 mitochondrial protein-coding genes of Formica species exhibited molecular signatures of purifying selection, as indicated by the ratio of non-synonymous substitutions to synonymous substitutions being less than 1 for each protein-coding gene. Phylogenetic analyses based on mitogenomic data obtained fairly consistent phylogenetic relationships, except for two Formica species that had unstable phylogenetic positions, indicating mitogenomic data are useful for constructing phylogenies of ants. Beyond characterizing two additional Formica mitogenomes, this study also provided some key evolutionary insights into Formica.

Across the highest mountain on earth: discordant phylogeographic patterns and recent dispersal of Tibetan stone loaches (Triplophysa) in the Himalayas

Phylogeographic congruence among co-distributed taxa is regarded as an inherent inference to vicariance events. Nonetheless, incongruent patterns of contemporary lineage divergence among taxa indicated that species differ in their response to common past events. To investigate the role of past events, ecological traits and lineage diversification time in shaping the contemporary phylogeographic patterns, comparative analyses were conducted for Tibetan stone loaches in the Himalayas using three gene markers and two ecological traits (depth of caudal peduncle in their length and presence/absence of posterior chamber of the air bladder). By a thorough sampling in two flanks of the Himalayas, the authors detected that phylogenetic breaks were spatially discordant and divergences of populations were also temporally asynchronous in co-distributed loaches. Estimated divergence time using fossil-calibrated node dating indicated that the Tibetan stone loaches colonised into the south flank of the Himalayas until the Pleistocene. The demographic expansions were also disconcerted between populations in north and south flanks, or east and west Himalayas. Ongoing gene flows between populations in north and south sides implied that the Himalayas do not strictly impede dispersal of cold-adapted species. The results highlight that the quaternary climatic oscillation, in conjunction with ecological traits and lineage diversification time, shaped contemporary phylogenetic patterns of stone loaches in the Himalayas and provide new insights into the biodiversity and composition of species in the Himalayas and surrounding region.

Punctuated equilibrium at 50: Anything there for evolutionary anthropology? Yes; definitely

The theory of punctuated equilibrium (PE) was developed a little over 50 years ago to explain long-term, large-scale appearance and disappearance of species in the fossil record. A theory designed specifically for that purpose cannot be expected, out of the box, to be directly applicable to biocultural evolution, but in revised form, PE offers a promising approach to incorporating not only a wealth of recent empirical research on genetic, linguistic, and technological evolution but also large databases that document human biological and cultural diversity across time and space. Here we isolate the fundamental components of PE and propose which pieces, when reassembled or renamed, can be highly useful in evolutionary anthropology, especially as humanity faces abrupt ecological challenges on an increasingly larger scale.

Genomic insights into the Montseny brook newt (Calotriton arnoldi), a Critically Endangered glacial relict

The Montseny brook newt (Calotriton arnoldi), considered the most endangered amphibian in Europe, is a relict salamandrid species endemic to a small massif located in northeastern Spain. Although conservation efforts should always be guided by genomic studies, those are yet scarce among urodeles, hampered by the extreme sizes of their genomes. Here, we present the third available genome assembly for the order Caudata, and the first genomic study of the species and its sister taxon, the Pyrenean brook newt (Calotriton asper), combining whole-genome and ddRADseq data. Our results reveal significant demographic oscillations which accurately mirrored Europe's climatic history. Although severe bottlenecks have led to depauperate genomic diversity and long runs of homozygosity along a gigantic genome, inbreeding might have been avoided by assortative mating strategies. Other life history traits, however, seem to have been less advantageous, and the lack of land dispersal has driven to exceptional levels of population fragmentation.

Are Appearances Deceiving? Morpho-Genetic Complexity of the Eumerus tricolor Group (Diptera: Syrphidae) in Europe, with a Focus on the Iberian Peninsula

Eumerus Meigen, 1822 is one of the largest Syrphidae genera in the Palaearctic Region, with the highest levels of taxonomic diversity found in the Eumerus tricolor species group. Despite its high diversity, the interspecific levels of morphological variability can be low. Additionally, some species may show certain levels of intraspecific variability. Hence, species delimitation may become challenging. In this work, we assessed the diversity of the E. tricolor group in the Iberian Peninsula through an integrative analysis of nomenclature, morphology and the 5' (COI-5') and 3' (COI-3') end regions of the Cytochrome c oxidase subunit I gene. Two new species, Eumerus ancylostylus Aguado-Aranda & Ricarte sp. n. and Eumerus petrarum Aguado-Aranda, Nedeljković & Ricarte sp. n., were described, and their intra- and interspecific variations discussed. In addition, the first barcodes of Iberian members of the E. tricolor group were obtained, and the distribution ranges of all species were mapped within the study area. The systematic position of the new species is discussed based on the resulting COI-based trees. The male genitalia of Eumerus hispanicus van der Goot, 1966 and Eumerus bayardi Séguy, 1961 were studied and illustrated. A lectotype was designated for Eumerus lateralis (Zetterstedt, 1819). An updated dichotomous key for all known European species of the E. tricolor group is provided. The egg of E. petrarum sp. n. is also described.

Integrating kelp genomic analyses and geological data to reveal ancient earthquake impacts

Detached buoyant kelp can disperse thousands of kilometres at sea and can colonize newly available shores in the wake of disturbances that wipe out competitors. Localized earthquake uplift can cause extirpation of intertidal kelp populations followed by recolonization. Sources of recolonizing kelp can be detectable in genomic structure of contemporary populations. Our field observations combined with LiDAR mapping identified a previously unrecognized zone of uplifted rocky coastline in a region that is slowly subsiding. Intertidal kelp (Durvillaea antarctica) on the uplifted section of coast is genetically distinctive from nearby populations, with genomic signatures most similar to that of kelp 300 km to the south. Genetic divergence between these locations suggests reproductive isolation for thousands of years. Combined geological and genetic data suggest that this uplift event occurred during one of four major earthquakes between 6000 and 2000 years ago, with one of the younger events most likely. Extirpation of the pre-existing kelp required sudden uplift of approximately 2 metres, precluding several small incremental uplift events. Our results show the power of integrating biological (genomic) analyses with geological data to understand ancient geological processes and their ecological impacts.

A clue to the evolutionary history of modern East Asian flora: insights from phylogeography and diterpenoid alkaloid distribution pattern of the Spiraea japonica complex

Each subkingdom of East Asian flora (EAF) has a unique evolutionary history, but which has rarely been described based on phylogeographic studies of EAF species. The Spiraea japonica L. complex, which is widespread in East Asia (EA), has received considerable attention because of the presence of diterpenoid alkaloids (DAs). It provides a proxy for understanding the genetic diversity and DA distribution patterns of species under various environmental conditions associated with the geological background in EA. In the present study, the plastome and chloroplast/nuclear DNA of 71 populations belonging to the S. japonica complex and its congeners were sequenced, combined with DA identification, environmental analyses, and ecological niche modelling, to investigate their phylogenetic relationships, genetic and DAs distribution patterns, biogeography, and demographic dynamics. An "ampliative" S. japonica complex was put forward, comprising all species of Sect. Calospira Ser. Japonicae, of which three evolutionary units carrying their respective unique types of DAs were identified and associated with the regionalization of EAF (referring to the Hengduan Mountains, central China, and east China). Moreover, a transition belt in central China with its biogeographic significance was revealed by genetic and DA distribution patterns from the perspective of ecological adaptation. The origin and onset differentiation of the "ampliative" S. japonica complex was estimated in the early Miocene (22.01/19.44 Ma). The formation of Japanese populations (6.75 Ma) was facilitated by the land bridge, which subsequently had a fairly stable demographic history. The populations in east China have undergone a founder effect after the Last Glacial Maximum, which may have been promoted by the expansion potential of polyploidization. Overall, the in-situ origin and diversification of the "ampliative" S. japonica complex since the early Miocene is a vertical section of the formation and development of modern EAF and was shaped by the geological history of each subkingdom.

Mitochondrial genomes reveal mid-Pleistocene population divergence, and post-glacial expansion, in Australasian snapper (Chrysophrys auratus)

Glacial cycles play important roles in determining the phylogeographic structure of terrestrial species, however, relatively little is known about their impacts on the distribution of marine biota. This study utilised modern (n = 350) and ancient (n = 26) mitochondrial genomes from Australasian snapper (Chrysophrys auratus) sampled in New Zealand to assess their demographic and phylogeographic history. We also tested for changes in genetic diversity using the up to 750-year-old mitochondrial genomes from pre-European archaeological sites to assess the potential impacts of human exploitation. Nucleotide diversity and haplotype diversity was high (π = 0.005, h = 0.972). There was no significant change in nucleotide diversity over the last 750 years (p = 0.343), with no detectable loss of diversity as a result of indigenous and industrial-scale fishing activity. While there was no evidence for contemporary population structure (AMOVA, p = 0.764), phylogeographic analyses identified two distinct mitochondrial clades that diverged approximately 650,000 years ago during the mid-Pleistocene, suggesting the species experienced barriers to gene flow when sea levels dropped over 120 m during previous glacial maxima. An exponential population increase was also observed around 8000 years ago consistent with a post-glacial expansion, which was likely facilitated by increased ocean temperatures and rising sea levels. This study demonstrates that glacial cycles likely played an important role in the demographic history of C. auratus and adds to our growing understanding of how dynamic climatic changes have influenced the evolution of coastal marine species.

New insights on patterns of genetic admixture and phylogeographic history in Iberian high mountain populations of midwife toads

Multiple Quaternary glacial refugia in the Iberian Peninsula, commonly known as "refugia within refugia", allowed diverging populations to come into contact and admix, potentially boosting substantial mito-nuclear discordances. In this study, we employ a comprehensive set of mitochondrial and nuclear markers to shed light onto the drivers of geographical differentiation in Iberian high mountain populations of the midwife toads Alytes obstetricans and A. almogavarii from the Pyrenees, Picos de Europa and Guadarrama Mountains. In the three analysed mountain regions, we detected evidence of extensive mito-nuclear discordances and/or admixture between taxa. Clustering analyses identified three major divergent lineages in the Pyrenees (corresponding to the eastern, central and central-western Pyrenees), which possibly recurrently expanded and admixed during the succession of glacial-interglacial periods that characterised the Late Pleistocene, and that currently follow a ring-shaped diversification pattern. On the other hand, populations from the Picos de Europa mountains (NW Iberian Peninsula) showed a mitochondrial affinity to central-western Pyrenean populations and a nuclear affinity to populations from the central Iberian Peninsula, suggesting a likely admixed origin for Picos de Europa populations. Finally, populations from the Guadarrama Mountain Range (central Iberian Peninsula) were depleted of genetic diversity, possibly as a consequence of a recent epidemic of chytridiomycosis. This work highlights the complex evolutionary history that shaped the current genetic composition of high mountain populations, and underscores the importance of using a multilocus approach to better infer the dynamics of population divergence.

A low-latitude species pump: Peripheral isolation, parapatric speciation and mating-system evolution converge in a marine radiation

Geologically recent radiations can shed light on speciation processes, but incomplete lineage sorting and introgressive gene flow render accurate evolutionary reconstruction and interpretation challenging. Independently evolving metapopulations of low dispersal taxa may provide an additional level of phylogeographic information, given sufficiently broad sampling and genome-wide sequencing. Evolution in the marine brown algal genus Fucus in the south-eastern North Atlantic was shaped by Quaternary climate-driven range shifts. Over this timescale, divergence and speciation occurred against a background of expansion-contraction cycles from multiple refugia, together with mating-system shifts from outcrossing (dioecy) to selfing hermaphroditism. We tested the hypothesis that peripheral isolation of range edge (dioecious) F. vesiculosus led to parapatric speciation and radiation of hermaphrodite lineages. Species tree methods using 876 single-copy nuclear genes and extensive geographic coverage produced conflicting topologies with respect to geographic clades of F. vesiculosus. All methods, however, revealed a new and early diverging hermaphrodite species, Fucus macroguiryi sp. nov. Both the multispecies coalescent and polymorphism-aware models (in contrast to concatenation) support sequential paraphyly in F. vesiculosus resulting from distinct evolutionary processes. Our results support (1) peripheral isolation of the southern F. vesiculosus clade prior to parapatric speciation and radiation of hermaphrodite lineages-a "low-latitude species pump". (2) Directional introgressive gene flow into F. vesiculosus around the present-day secondary contact zone (sympatric-allopatric boundary) between dioecious/hermaphrodite lineages as hermaphrodites expanded northwards, supported by concordance analysis and statistical tests of introgression. (3) Species boundaries in the extensive sympatric range are probably maintained by reproductive system (selfing in hermaphrodites) and reinforcement.

Adaptive responses drive the success of polyploid yellowcresses (Rorippa, Brassicaceae) in the Hengduan Mountains, a temperate biodiversity hotspot

Polyploids contribute substantially to plant evolution and biodiversity; however, the mechanisms by which they succeed are still unclear. According to the polyploid adaptation hypothesis, successful polyploids spread by repeated adaptive responses to new environments. Here, we tested this hypothesis using two tetraploid yellowcresses (Rorippa), the endemic Rorippa elata and the widespread Rorippa palustris, in the temperate biodiversity hotspot of the Hengduan Mountains. Speciation modes were resolved by phylogenetic modeling using 12 low-copy nuclear loci. Phylogeographical patterns were then examined using haplotypes phased from four plastid and ITS markers, coupled with historical niche reconstruction by ecological niche modeling. We inferred the time of hybrid origins for both species as the mid-Pleistocene, with shared glacial refugia within the southern Hengduan Mountains. Phylogeographic and ecological niche reconstruction indicated recurrent northward colonization by both species after speciation, possibly tracking denuded habitats created by glacial retreat during interglacial periods. Common garden experiment involving perennial R. elata conducted over two years revealed significant changes in fitness-related traits across source latitudes or altitudes, including latitudinal increases in survival rate and compactness of plant architecture, suggesting gradual adaptation during range expansion. These findings support the polyploid adaptation hypothesis and suggest that the spread of polyploids was aided by adaptive responses to environmental changes during the Pleistocene. Our results thus provide insight into the evolutionary success of polyploids in high-altitude environments.

Parallel recolonizations generate distinct genomic sectors in kelp following high-magnitude earthquake disturbance

Large-scale disturbance events have the potential to drastically reshape biodiversity patterns. Notably, newly vacant habitat space cleared by disturbance can be colonized by multiple lineages, which can lead to the evolution of distinct spatial "sectors" of genetic diversity within a species. We test for disturbance-driven sectoring of genetic diversity in intertidal southern bull kelp, Durvillaea antarctica (Chamisso) Hariot, following the high-magnitude 1855 Wairarapa earthquake in New Zealand. Specifically, we use genotyping-by-sequencing (GBS) to analyse fine-scale population structure across the uplift zone and apply machine learning to assess the fit of alternative recolonizaton models. Our analysis reveals that specimens from the uplift zone carry distinctive genomic signatures potentially linked to post-earthquake recolonization processes. Specifically, our analysis identifies two parapatric spatial-genomic sectors of D. antarctica at Turakirae Head, which experienced the most dramatic uplift. Based on phylogeographical modelling, we infer that bull kelp in the Wellington region was probably a source for recolonization of the heavily uplifted Turakirae Head coastline, via two parallel, eastward recolonization events. By identifying multiple parapatric genotypic sectors within a recently recolonized coastal region, the current study provides support for the hypothesis that competing lineage expansions can generate striking spatial structuring of genetic diversity, even in highly dispersive taxa.

Mitogenomes resolve the phylogeography and divergence times within the endemic New Zealand Callaeidae (Aves: Passerida)

The biogeographical origins of the endemic birds of New Zealand (Aotearoa) are of great interest, particularly Palaeogene lineages such as Callaeidae, a passerine family characterized by brightly coloured wattles behind the beak and, in some cases, extreme sexual dimorphism in bill size and shape. Ancestral representatives of Callaeidae are thought to have split from their closest relatives outside New Zealand in the Oligocene, but little is known about the timing of divergences within the family. We present a fully dated molecular phylogeny of Callaeidae mitogenomes and discuss the biogeographical implications. Our results suggest that formation of Pliocene marine seaways, such as the Manawatu Strait, are likely to have played a significant role in the differentiation of North Island and South Island kōkako (Callaeas spp.) and saddlebacks/tīeke (Philesturnus spp.).

The role of multiple Pleistocene refugia in promoting diversification in the Pacific Northwest

Pleistocene glacial cycles drastically changed the distributions of taxa endemic to temperate rainforests in the Pacific Northwest, with many experiencing reduced habitat suitability during glacial periods. In this study, we investigate whether glacial cycles promoted intraspecific divergence and whether subsequent range changes led to secondary contact and gene flow. For seven invertebrate species endemic to the PNW, we estimated Species Distribution Models (SDMs) and projected them onto current and historical climate conditions to assess how habitat suitability changed during glacial cycles. Using single nucleotide polymorphism (SNP) data from these species, we assessed population genetic structure and used a machine-learning approach to compare models with and without gene flow between populations upon secondary contact after the Last Glacial Maximum (LGM). Finally, we estimated divergence times and rates of gene flow between populations. SDMs suggest that there was less suitable habitat in the North Cascades and Northern Rocky Mountains during glacial compared to interglacial periods, resulting in reduced habitat suitability and habitat fragmentation during the LGM. Our genomic data identify population structure in all taxa, and support gene flow upon secondary contact in five of the seven taxa. Parameter estimates suggest that population divergences date to the later Pleistocene for most populations. Our results support a role of refugial dynamics in driving intraspecific divergence in the Cascades Range. In these invertebrates, population structure often does not correspond to current biogeographic or environmental barriers. Rather, population structure may reflect refugial lineages that have since expanded their ranges, often leading to secondary contact between once isolated lineages.

Phylogeography and Cryptic Species Structure of a Locally Adapted Parasite in New Zealand

Phylogeographic patterns of many taxa on New Zealand's South Island are characterised by disjunct distributions that have been attributed to Pleistocene climatic cycles and the formation of the Southern Alps. Pleistocene glaciation has also been implicated in shaping the contemporary genetic differentiation between populations of the aquatic snail Potamopyrgus antipodarum. We investigated whether similar phylogeographic patterns exist for the snail's locally adapted trematode parasite, Atriophallophorus winterbourni. We found evidence for a barrier to gene-flow in sympatry between cryptic, but ecologically divergent species. When focusing on the most common of these species, disjunct geographic distributions are found for mitochondrial lineages that diverged during the Pleistocene. The boundary between these distributions is found in the central part of the South Island and reinforced by low cross-alpine migration. Further support for a vicariant origin of the phylogeographic pattern was found when assessing nuclear multilocus SNP data. Nuclear and mitochondrial population differentiation was concordant in pattern, except for populations in a potential secondary contact zone. Additionally, we found larger than expected differentiation between nuclear- and mitochondrial-based empirical Bayes F_ST estimates (global F_ST : 0.02 versus 0.39 for nuclear and mitochondrial data, respectively). Population subdivision is theoretically expected to be stronger for mitochondrial genomes due to a smaller effective population size, but the strong difference here, together with mito-nuclear discordance in a putative contact zone, is potentially indicative of divergent gene flow of nuclear and mitochondrial genomes.

Diversity patterns of cushion plants on the Qinghai-Tibet Plateau: A basic study for future conservation efforts on alpine ecosystems

The Qinghai-Tibet Plateau (QTP) is an important cushion plant hotspot. However, the distribution of cushion plants on the QTP is unknown, as are the factors that drive cushion plant distribution, limiting our understanding of the evolution of cushion species in the region. In this study, we assessed spatial patterns of total cushion plant diversity (including taxonomic and phylogenetic) over the entire QTP and compared patterns of diversity of cushion plants with different typologies (i.e., compact vs. loose). We also examined how these patterns were related to climatic features. Our results indicate that the southern QTP hosts the highest total cushion plant richness, especially in the south-central Hengduan Mountains subregion. The total number of cushion species declines from south to north and from southeast to northwest. Compact cushion plants exhibit similar patterns as the total cushion plant richness, whereas loose cushion plants show random distribution. Cushion plant phylogenetic diversity showed a similar pattern as that of the total cushion plant richness. In addition, cushion plant phylogenetic community structure was clustered in the eastern and southwestern QTP, whereas random or overdispersed in other areas. Climatic features represented by annual energy and water trends, seasonality and extreme environmental factors, had significant effects on cushion plant diversity patterns but limited effects on the phylogenetic community structure, suggesting that climatic features indeed promote the formation of cushion plants. Because cushion plants play vital roles in alpine ecosystems, our findings not only promote our understanding of the evolution and formation of alpine cushion plant diversity but also provide an indispensable foundation for future studies on cushion plant functions and thus alpine ecosystem sustainability in the entire QTP region.

Fine-scale genome-wide signature of Pleistocene glaciation in Thitarodes moths (Lepidoptera: Hepialidae), host of Ophiocordyceps fungus in the Hengduan Mountains

The Hengduan Mountains region is a biodiversity hotspot known for its topologically complex, deep valleys and high mountains. While landscape and glacial refugia have been evoked to explain patterns of inter-species divergence, the accumulation of intra-species (i.e. population level) genetic divergence across the mountain-valley landscape in this region has received less attention. We used genome-wide restriction site-associated DNA sequencing (RADseq) to reveal signatures of Pleistocene glaciation in populations of Thitarodes shambalaensis (Lepidoptera: Hepialidae), the host moth of parasitic Ophiocordyceps sinensis (Hypocreales: Ophiocordycipitaceae) or "caterpillar fungus" endemic to the glacier of eastern Mt. Gongga. We used moraine history along the glacier valleys to model the distribution and environmental barriers to gene flow across populations of T. shambalaensis. We found that moth populations separated by less than 10 km exhibited valley-based population genetic clustering and isolation-by-distance (IBD), while gene flow among populations was best explained by models using information about their distributions at the local last glacial maximum (LGM_L , 58 kya), not their contemporary distribution. Maximum likelihood lineage history among populations, and among subpopulations as little as 500 meters apart, recapitulated glaciation history across the landscape. We also found signals of isolated population expansion following the retreat of LGM_L glaciers. These results reveal the fine-scale, long-term historical influence of landscape and glaciation on the genetic structuring of populations of an endangered and economically important insect species. Similar mechanisms, given enough time and continued isolation, could explain the contribution of glacier refugia to the generation of species diversity among the Hengduan Mountains.

Major population splits coincide with episodes of rapid climate change in a forest-dependent bird

Climate change influences population demography by altering patterns of gene flow and reproductive isolation. Direct mutation rates offer the possibility for accurate dating on the within-species level but are currently only available for a handful of vertebrate species. Here, we use the first directly estimated mutation rate in birds to study the evolutionary history of pied flycatchers (Ficedula hypoleuca). Using a combination of demographic inference and species distribution modelling, we show that all major population splits in this forest-dependent system occurred during periods of increased climate instability and rapid global temperature change. We show that the divergent Spanish subspecies originated during the Eemian-Weichselian transition 115-104 thousand years ago (kya), and not during the last glacial maximum (26.5-19 kya), as previously suggested. The magnitude and rates of climate change during the glacial-interglacial transitions that preceded population splits in pied flycatchers were similar to, or exceeded, those predicted to occur in the course of the current, human-induced climate crisis. As such, our results provide a timely reminder of the strong impact that episodes of climate instability and rapid temperature changes can have on species' evolutionary trajectories, with important implications for the natural world in the Anthropocene.

Geographical isolation versus dispersal: Relictual alpine grasshoppers support a model of interglacial diversification with limited hybridization

Alpine biotas are paradigmatic of the countervailing roles of geographical isolation and dispersal during diversification. In temperate regions, repeated distributional shifts driven by Pleistocene climatic oscillations produced both recurrent pulses of population fragmentation and opportunities for gene flow during range expansions. Here, we test whether a model of divergence in isolation vs. with gene flow is more likely in the diversification of flightless alpine grasshoppers of the genus Podisma from the Iberian Peninsula. The answer to this question can also provide key insights about the pace of evolution. Specifically, if the data fit a divergence in isolation model, this suggests rapid evolution of reproductive isolation. Genomic data confirm a Pleistocene origin of the species complex, and multiple analytical approaches revealed limited asymmetric historical hybridization between two taxa. Genomic-based demographic reconstructions, spatial patterns of genetic structure and range shifts inferred from palaeodistribution modelling suggest severe range contraction accompanied by declines in effective population sizes during interglacials (i.e., contemporary populations confined to sky islands are relicts) and expansions during the coldest stages of the Pleistocene in each taxon. Although limited hybridization during secondary contact leads to phylogenetic uncertainty if gene flow is not accommodated when estimating evolutionary relationships, all species exhibit strong genetic cohesiveness. Our study lends support to the notion that the accumulation of incipient differences during periods of isolation were sufficient to lead to lineage persistence, but also that the demographic changes, dispersal constraints and spatial distribution of the sky islands themselves mediated species diversification in temperate alpine biotas.

Phylogeography of the endemic red-tailed cicadas of New Zealand (Hemiptera: Cicadidae: Rhodopsalta), and molecular, morphological and bioacoustical confirmation of the existence of Hudson’s Rhodopsalta microdora

Why do some genera radiate, whereas others do not? The genetic structure of present-day populations can provide clues for developing hypotheses. In New Zealand, three Cicadidae genera are depauperate [Amphipsalta (three species), Notopsalta (one species) and Rhodopsalta (three species)], whereas two have speciated extensively [Kikihia (~30 species/subspecies) and Maoricicada (~20 species/subspecies). Here, we examine the evolution of Rhodopsalta, the last New Zealand genus to be studied phylogenetically and phylogeographically. We use Bayesian and maximum-likelihood analyses of mitochondrial cox1 and nuclear EF1α gene sequences. Concatenated and single-gene phylogenies for 70 specimens (58 localities) support its monophyly and three described species: Rhodopsalta cruentata, Rhodopsalta leptomera and Rhodopsalta microdora, the last taxon previously regarded as uncertain. We provide distribution maps, biological notes and the first descriptions of diagnostic songs. We show that both R. cruentata and R. microdora exhibit northern and southern genetic subclades. Subclades of the dry-adapted R. microdora clade show geographical structure, whereas those of the mesic R. cruentata and sand-dune specialist R. leptomera have few discernible patterns. Genetic, bioacoustical and detailed distributional evidence for R. microdora add to the known biodiversity of New Zealand. We designate a lectotype for Tettigonia cruentata  Fabricius, 1775, the type species of Rhodopsalta.

Inter-glacial isolation caused divergence of cold-adapted species: the case of the snow partridge

Deciphering the role of climatic oscillations in species divergence helps us understand the mechanisms that shape global biodiversity. The cold-adapted species may have expanded their distribution with the development of glaciers during glacial period. With the retreat of glaciers, these species were discontinuously distributed in the high-altitude mountains and isolated by geographical barriers. However, the study that focuses on the speciation process of cold-adapted species is scant. To fill this gap, we combined population genetic data and ecological niche models (ENMs) to explore divergence process of snow partridge (Lerwa lerwa). Lerwa lerwa is a cold-adapted bird that is distributed from 4,000 to 5,500 m. We found 2 genetic populations within L. lerwa, and they diverged from each other at about 0.40–0.44 million years ago (inter-glacial period after Zhongliangan glaciation). The ENMs suggested that L. lerwa expanded to the low elevations of the Himalayas and Hengduan mountains during glacial period, whereas it contracted to the high elevations, southern of Himalayas, and Hengduan mountains during inter-glacial periods. Effective population size trajectory also suggested that L. lerwa expanded its population size during the glacial period. Consistent with our expectation, the results support that inter-glacial isolation contributed to the divergence of cold-adapted L. lerwa on Qinghai-Tibetan Plateau. This study deepens our understanding of how climatic oscillations have driven divergence process of cold-adapted Phasianidae species distributed on mountains.

Signature of mid-Pleistocene lineages in the European silver fir (Abies alba Mill.) at its geographic distribution margin

In a conservation and sustainable management perspective, we identify the ecological, climatic, and demographic factors responsible for the genetic diversity patterns of the European silver fir (Abies alba Mill.) at its southwestern range margin (Pyrenees Mountains, France, Europe). We sampled 45 populations throughout the French Pyrenees and eight neighboring reference populations in the Massif Central, Alps, and Corsica. We genotyped 1,620 individuals at three chloroplast and ten nuclear microsatellite loci. We analyzed within- and among-population genetic diversity using phylogeographic reconstructions, tests of isolation-by-distance, Bayesian population structure inference, modeling of demographic scenarios, and regression analyses of genetic variables with current and past environmental variables. Genetic diversity decreased from east to west suggesting isolation-by-distance from the Alps to the Pyrenees and from the Eastern to the Western Pyrenees. We identified two Pyrenean lineages that diverged from a third Alpine-Corsica-Massif Central lineage 0.8 to 1.1 M years ago and subsequently formed a secondary contact zone in the Central Pyrenees. Population sizes underwent contrasted changes, with a contraction in the west and an expansion in the east. Glacial climate affected the genetic composition of the populations, with the western genetic cluster only observed in locations corresponding to the coldest past climate and highest elevations. The eastern cluster was observed over a larger range of temperatures and elevations. All demographic events shaping the current spatial structure of genetic diversity took place during the Mid-Pleistocene Transition, long before the onset of the Holocene. The Western Pyrenees lineage may require additional conservation efforts, whereas the eastern lineage is well protected in in situ gene conservation units. Due to past climate oscillations and the likely emergence of independent refugia, east-west oriented mountain ranges may be important reservoir of genetic diversity in a context of past and ongoing climate change in Europe.

Spatiotemporal maintenance of flora in the Himalaya biodiversity hotspot: Current knowledge and future perspectives

Mountain ecosystems support a significant one-third of all terrestrial biodiversity, but our understanding of the spatiotemporal maintenance of this high biodiversity remains poor, or at best controversial. The Himalaya hosts a complex mountain ecosystem with high topographic and climatic heterogeneity and harbors one of the world's richest floras. The high species endemism, together with increasing anthropogenic threats, has qualified the Himalaya as one of the most significant global biodiversity hotspots. The topographic and climatic complexity of the Himalaya makes it an ideal natural laboratory for studying the mechanisms of floral exchange, diversification, and spatiotemporal distributions. Here, we review literature pertaining to the Himalaya in order to generate a concise synthesis of the origin, distribution, and climate change responses of the Himalayan flora. We found that the Himalaya supports a rich biodiversity and that the Hengduan Mountains supplied the majority of the Himalayan floral elements, which subsequently diversified from the late Miocene onward, to create today's relatively high endemicity in the Himalaya. Further, we uncover links between this Miocene diversification and the joint effect of geological and climatic upheavals in the Himalaya. There is marked variance regarding species dispersal, elevational gradients, and impact of climate change among plant species in the Himalaya, and our review highlights some of the general trends and recent advances on these aspects. Finally, we provide some recommendations for conservation planning and future research. Our work could be useful in guiding future research in this important ecosystem and will also provide new insights into the maintenance mechanisms underpinning other mountain systems.

Phylogeography of the stonefly Kamimuria tibialis: multiple glacial refugia and sympatric occurrence of different lineages in the southern islands of Japan

To elucidate the effect of Pleistocene climatic fluctuations on the historical distribution and geographical genetic structure of temperate Japanese species, we performed phylogeographical and demographic analyses using mitochondrial gene sequences obtained from the stonefly species Kamimuria tibialis, sampled from four main islands of the Japanese Archipelago (i.e. Hokkaido, Honshu, Shikoku and Kyushu) and Tsushima Island. We detected three main clades with distinct geographical distributions, including the Tsushima, Kyushu and Hokkaido–Honshu–Shikoku phylogroups. These groups were estimated to have diverged from one another 0.54–2.02 Mya, suggesting they have undergone several glacial cycles in different refugia. Our results showed that during the glacial epochs and with a fall in sea-level, gene flow was limited among Tsushima and Kyushu, and among Hokkaido and Honshu, probably because the straits between these islands are deep. The population in Kyushu and Shikoku, the southernmost islands, exhibited high genetic diversity, with two distinct haplotype lineages occurring sympatrically. These results suggest that the population division into multiple refugia and the existence of stable southern refugia have contributed to the high genetic diversity of the species in this region.

Global changes explain the long-term demographic trend of the Eurasian common lizard (Squamata: Lacertidae)

The demographic trend of a species depends on the dynamics of its local populations, which can be compromised by local or by global phenomena. However, the relevance of local and global phenomena has rarely been investigated simultaneously. Here, we tested whether local phenomena compromised a species’ demographic trend using the Eurasian common lizard Zootoca vivipara, the terrestrial reptile exhibiting the widest geographic distribution, as a model species. We analyzed the species’ ancient demographic trend using genetic data from its 6 allopatric genetic clades and tested whether its demographic trend mainly depended on single clades or on global phenomena. Zootoca vivipara’s effective population size increased since 2.3 million years ago and started to increase steeply and continuously from 0.531 million years ago. Population growth rate exhibited 2 maxima, both occurring during global climatic changes and important vegetation changes on the northern hemisphere. Effective population size and growth rate were negatively correlated with global surface temperatures, in line with global parameters driving long-term demographic trends. Zootoca vivipara’s ancient demography was neither driven by a single clade, nor by the 2 clades that colonized huge geographic areas after the last glaciation. The low importance of local phenomena, suggests that the experimentally demonstrated high sensitivity of this species to short-term ecological changes is a response in order to cope with short-term and local changes. This suggests that what affected its long-term demographic trend the most, were not these local changes/responses, but rather the important and prolonged global climatic changes and important vegetation changes on the northern hemisphere, including the opening up of the forest by humans.

Spatial phylogenetics of two topographic extremes of the Hengduan Mountains in southwestern China and its implications for biodiversity conservation

Previous attempts to elucidate the drivers of speciation mechanisms and spatial distribution patterns of biodiversity in mountain regions have treated different floras within a single geological region as one flora, ignoring the potential contributions of high habitat/ecosystem heterogeneity. Furthermore, current conservation strategies largely focus on forest ecosystems and/or specific flagship species, ignoring marginal ecosystems, leaving species in these ecosystems at risk. Here, we compared the spatial patterns of biodiversity and the potential drivers of these patterns in the river valley and subnival ecosystems of the Hengduan Mountains region (HDM) in southwestern China. Specifically, we compared spatial patterns of diversity, endemism, and threatened species in these ecosystems based on both traditional measurements and recent phylogenetic approaches. We then examined how those patterns were related to environmental factors and human activity in these same regions. We found that the middle-southern HDM supports the highest diversity and endemism for the river valley and subnival ecosystems; however, the distribution patterns of neo- and paleo-endemism in these two ecosystems differ. Regression models indicate that habitat diversity and paleo-climatic fluctuation are important drivers of diversity and endemism for these two ecosystems. Temperature and precipitation, however, showed different influences on the spatial patterns in different ecosystems. Categorical analysis of neo- and paleo-endemism (CANAPE) indicated that most endemism centers are not covered by current nature reserves. Moreover, the intensity of human activity is highest in the southern and southeastern HDM, which coincides with the distribution patterns of diversity, mixed-endemism and high-priority (and threatened) species. These findings suggest that different floras within a single geographic/floristic region respond differently to environmental factors and show different spatial phylogenetic patterns. We, therefore, recommend that future research into the drivers of biodiversity consider the contributions of various ecosystem types within a single geological region. This study also provides a theoretical basis for protecting habitat diversity. Our work confirms that current conservation efforts are insufficient to protect ecosystem diversity in the river valley and subnival ecosystems of the Hengduan Mountains. Therefore, we recommend the establishment of nature reserves in the regions identified in this study; furthermore, we strongly recommend improving current and establishing new management policies for biodiversity conservation in this region.

Alpine speciation and morphological innovations: revelations from a species-rich genus in the northern hemisphere

A large number of studies have attempted to determine the mechanisms driving plant diversity and distribution on a global scale, but the diverse and endemic alpine herbs found in harsh environments, showing adaptive evolution, require more studies. Here, we selected 466 species from the genus Saussurea, one of the northern hemisphere's highest-altitude plant genera with high species richness and striking morphological traits, to explore the mechanisms driving speciation and adaptative evolution. We conducted phylogenetic signals analysis and ancestral character estimation to explore the phylogenetic significance of ecological factors. Moreover, we used spatial simultaneous autoregressive (SAR) error models, modified t-tests and partial regression models to quantify the relative effects of ecological factors and morphological diversity upon diversity and endemism of Saussurea. Phylogenetic analyses reveal that geological influences and climate stability exhibit significant phylogenetic signals and that Saussurea originated at a relatively high elevation. Regression models indicate that geological influences and climatic stability significantly affect the diversity and endemism patterns of Saussurea and its morphological innovations. Moreover, morphological innovations in an area show significant contributions to the local diversity and endemism of Saussurea. We conclude that geological influences (mean altitude and topographic heterogeneity), glacial-interglacial climate stability and phylogenetic conservatism have together promoted the speciation and adaptive evolution of the genus Saussurea. In addition, adaptively morphological innovations of alpine species also promote diversification in local regions. Our findings improve the understanding of the distribution pattern of diversity/endemism and adaptive evolution of alpine specie in the whole northern hemisphere.

Genomic Signature of Shifts in Selection in a Subalpine Ant and Its Physiological Adaptations

Understanding how organisms adapt to extreme environments is fundamental and can provide insightful case studies for both evolutionary biology and climate-change biology. Here, we take advantage of the vast diversity of lifestyles in ants to identify genomic signatures of adaptation to extreme habitats such as high altitude. We hypothesized two parallel patterns would occur in a genome adapting to an extreme habitat: 1) strong positive selection on genes related to adaptation and 2) a relaxation of previous purifying selection. We tested this hypothesis by sequencing the high-elevation specialist Tetramorium alpestre and four other phylogenetically related species. In support of our hypothesis, we recorded a strong shift of selective forces in T. alpestre, in particular a stronger magnitude of diversifying and relaxed selection when compared with all other ants. We further disentangled candidate molecular adaptations in both gene expression and protein-coding sequence that were identified by our genome-wide analyses. In particular, we demonstrate that T. alpestre has 1) a higher level of expression for stv and other heat-shock proteins in chill-shock tests and 2) enzymatic enhancement of Hex-T1, a rate-limiting regulatory enzyme that controls the entry of glucose into the glycolytic pathway. Together, our analyses highlight the adaptive molecular changes that support colonization of high-altitude environments.

Global topographic uplift has elevated speciation in mammals and birds over the last 3 million years

Topographic change shapes the evolution of biodiversity by influencing both habitat connectivity and habitat diversity as well as abiotic factors like climate. However, its role in creating global biodiversity gradients remains poorly characterized because geology, climate and evolutionary data have rarely been integrated across concordant timescales. Here we show that topographic uplift over the last 3 million years explains more spatial variation in the speciation of all mammals and birds than do the direct effects of palaeoclimate change and both present-day elevation and present-day temperature. By contrast, the effects of topographic changes are much smaller than those of present-day temperatures in eroded areas. Together, our results stress that historical geological processes rather than traditionally studied macroecological gradients may ultimately generate much of the world's biodiversity. More broadly, as the Earth's surface continues to rise and fall, topography will remain an important driver of evolutionary change and novelty.

Drainage basins serve as multiple glacial refugia for alpine habitats in the Sierra Nevada Mountains, California

The evolutionary histories of alpine species are often directly associated with responses to glaciation. Deep divergence among populations and complex patterns of genetic variation have been inferred as consequences of persistence within glacier boundaries (i.e., on nunataks), while shallow divergence and limited genetic variation are assumed to result from expansion from large refugia at the edge of ice shields (i.e., massifs de refuge). However, for some species, dependence on specific microhabitats could profoundly influence their spatial and demographic response to glaciation, and such a simple dichotomy may obscure the localization of actual refugia. In this study, we use the Nebria ingens complex (Coleoptera: Carabidae), a water-affiliated ground beetle lineage, to test how drainage basins are linked to their observed population structure. By analysing mitochondrial COI gene sequences and genome-wide single nucleotide polymorphisms, we find that the major drainage systems of the Sierra Nevada Mountains in California best explain the population structure of the N. ingens complex. In addition, we find that an intermediate morphotype within the N. ingens complex is the product of historical hybridization of N. riversi and N. ingens in the San Joaquin basin during glaciation. This study highlights the importance of considering ecological preferences in how species respond to climate fluctuations and provides an explanation for discordances that are often observed in comparative phylogeographical studies.

Dispersal barriers and opportunities drive multiple levels of phylogeographic concordance in the Southern Alps of New Zealand

Phylogeographic concordance, or the sharing of phylogeographic patterns among codistributed species, suggests similar responses to topography or climatic history. While the orientation and timing of breaks between lineages are routinely compared, spatial dynamics within regions occupied by individual lineages provide a second opportunity for comparing responses to past events. In environments with complex topography and glacial history, such as New Zealand's South Island, geographically nested comparisons can identify the processes leading to phylogeographic concordance between and within regional genomic clusters. Here, we used single nucleotide polymorphisms (obtained via ddRADseq) for two codistributed forest beetle species, Agyrtodes labralis (Leiodidae) and Brachynopus scutellaris (Staphylinidae), to evaluate the role of climate change and topography in shaping phylogeographic concordance at two, nested spatial scales: do species diverge over the same geographic barriers, with similar divergence times? And within regions delimited by these breaks, do species share similar spatial dynamics of directional expansion or isolation-by-distance? We found greater congruence of phylogeographic breaks between regions divided by the strongest dispersal barriers (i.e., the Southern Alps). However, these shared breaks were not indicative of shared spatial dynamics within the regions they delimit, and the most similar spatial dynamics between species occurred within regions with the strongest gradients in historical climatic stability. Our results indicate that lack of concordance as traditionally detected by lineage turnover does not rule out the possibility of shared histories, and variation in the presence and type of concordance may provide insights into the different processes shaping phylogeographic patterns across geologically dynamic regions.

Multiple glacial refugia and contemporary dispersal shape the genetic structure of an endemic amphibian from the Pyrenees

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species' present-day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present-day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex-biased dispersal. This work highlights the importance of integrating past evolutionary processes and present-day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats.

The genomic footprint of coastal earthquake uplift

Theory suggests that catastrophic earth-history events can drive rapid biological evolution, but empirical evidence for such processes is scarce. Destructive geological events such as earthquakes can represent large-scale natural experiments for inferring such evolutionary processes. We capitalized on a major prehistoric (800 yr BP) geological uplift event affecting a southern New Zealand coastline to test for the lasting genomic impacts of disturbance. Genome-wide analyses of three co-distributed keystone kelp taxa revealed that post-earthquake recolonization drove the evolution of novel, large-scale intertidal spatial genetic 'sectors' which are tightly linked to geological fault boundaries. Demographic simulations confirmed that, following widespread extirpation, parallel expansions into newly vacant habitats rapidly restructured genome-wide diversity. Interspecific differences in recolonization mode and tempo reflect differing ecological constraints relating to habitat choice and dispersal capacity among taxa. This study highlights the rapid and enduring evolutionary effects of catastrophic ecosystem disturbance and reveals the key role of range expansion in reshaping spatial genetic patterns.

An exploration of the complex biogeographical history of the Neotropical banner-wing damselflies (Odonata: Polythoridae)

BACKGROUND

The New World Tropics has experienced a dynamic landscape across evolutionary history and harbors a high diversity of flora and fauna. While there are some studies addressing diversification in Neotropical vertebrates and plants, there is still a lack of knowledge in arthropods. Here we examine temporal and spatial diversification patterns in the damselfly family Polythoridae, which comprises seven genera with a total of 58 species distributed across much of Central and South America.

RESULTS

Our time-calibrated phylogeny for 48 species suggests that this family radiated during the late Eocene (~ 33 Ma), diversifying during the Miocene. As with other neotropical groups, the Most Recent Common Ancestor (MRCA) of most of the Polythoridae genera has a primary origin in the Northern Andes though the MRCA of at least one genus may have appeared in the Amazon Basin. Our molecular clock suggests correlations with some major geographical events, and our biogeographical modeling (with BioGeoBEARS and RASP) found a significant influence of the formation of the Pebas and Acre systems on the early diversification of these damselflies, though evidence for the influence of the rise of the different Andean ranges was mixed. Diversification rates have been uniform in all genera except one-Polythore-where a significant increase in the late Pliocene (~ 3 mya) may have been influenced by recent Andean uplift.

CONCLUSION

The biogeographical models implemented here suggest that the Pebas and Acre Systems were significant geological events associated with the diversification of this damselfly family; while diversification in the tree shows some correlation with mountain building events, it is possible that other abiotic and biotic changes during our study period have influenced diversification as well. The high diversification rate observed in Polythore could be explained by the late uplift of the Northern Andes. However, it is possible that other intrinsic factors like sexual and natural selection acting on color patterns could be involved in the diversification of this genus.

River Capture and Freshwater Biological Evolution: A Review of Galaxiid Fish Vicariance

Geological processes can strongly affect the distribution and diversification of freshwater-limited species. In particular, a combination of geological and biological data has suggested that Earth history processes can drive vicariant isolation and speciation in non-migratory freshwater fishes. Here, we synthesise recently published geological and freshwater phylogeographic data to illustrate that changes in river drainage geometry are important drivers of galaxiid diversification, both in New Zealand and elsewhere. Major river capture events have led to the isolation and divergence of unique and geographically-restricted lineages, including taxa that are now of prime conservation concern. The parallel phylogeographic effects of drainage shifts have been verified by observations of concordant patterns in co-distributed species. Broadly, this study highlights the dynamic interplay between physical and biological processes in geologically active settings.

Using demographic model selection to untangle allopatric divergence and diversification mechanisms in the Rheum palmatum complex in the Eastern Asiatic Region

Allopatric divergence is often initiated by geological uplift and restriction to sky-islands, climate oscillations, or river capture. However, it can be difficult to establish which mechanism was the most likely to generate the current phylogeographical structure of a species. Recently, genomic data in conjunction with a model testing framework have been applied to address this issue in animals. To test whether such an approach is also likely to be successful in plants, we used population genomic data of the Rheum palmatum complex from the Eastern Asiatic Region, in conjunction with biogeographical reconstruction and demographic model selection, to identify the potential mechanism(s) which have led to the current level of divergence. Our results indicate that the R. palmatum complex originated in the central Hengduan Mts and possibly in regions further east, and then dispersed westward and eastward resulting in genetically distinct lineages. Populations are likely to have diverged in refugia during climate oscillations followed by subsequent expansion and secondary contact. However, model simulations within the western lineage of the R. palmatum complex cannot reject a restriction to sky-islands as a possible mechanism of diversification due to the genetically ambiguous position of one population. This highlights that genetically mixed populations might introduce ambiguity regarding the best diversification model in some cases. Although it might be possible to resolve this ambiguity using other data, sometimes this could prove to be difficult in complex biogeographical areas.

Rapid speciation and ecological divergence into North American alpine habitats: the Nippononebria (Coleoptera: Carabidae) species complex

The climate-driven species pump hypothesis has been supported in a number of phylogeographic studies of alpine species. Climate-driven shifts in distribution, coupled with rapid demographic change, have led to strong genetic drift and lineage diversification. Although the species pump has been linked to rapid speciation in a number of studies, few studies have demonstrated that ecological divergence accompanies rapid speciation. Here we examine genetic, morphological and physiological variation in members of the ground beetle taxon Nippononebria, to test three competing hypotheses of evolutionary diversification: isolation and incomplete lineage sorting (no speciation), recent speciation without ecological divergence, or recent speciation with ecological divergence into alpine habitats. Genetic data are consistent with recent divergence, with major lineages forming in the last million years. A species tree analysis, in conjunction with morphological divergence in male reproductive traits, support the formation of three recognized Nippononebria taxa. Furthermore, both morphological and physiological traits demonstrate ecological divergence in alpine lineages, with convergent shifts in body shape and thermal tolerance breadth. This provides strong evidence that the climate-driven species pump can generate ecological novelty, though it is argued that spatial scale may be a key determinant of broader patterns of macroevolution in alpine communities.

Speciation in protists: Spatial and ecological divergence processes cause rapid species diversification in a freshwater chrysophyte

Although eukaryotic microorganisms are extremely numerous, diverse and essential to global ecosystem functioning, they are largely understudied by evolutionary biologists compared to multicellular macroscopic organisms. In particular, very little is known about the speciation mechanisms which may give rise to the diversity of microscopic eukaryotes. It was postulated that the enormous population sizes and ubiquitous distribution of these organisms could lead to a lack of population differentiation and therefore very low speciation rates. However, such assumptions have traditionally been based on morphospecies, which may not accurately reflect the true diversity, missing cryptic taxa. In this study, we aim to articulate the major diversification mechanisms leading to the contemporary molecular diversity by using a colonial freshwater flagellate, Synura sphagnicola, as an example. Phylogenetic analysis of five sequenced loci showed that S. sphagnicola differentiated into two morphologically distinct lineages approximately 15.4 million years ago, which further diverged into several evolutionarily recent haplotypes during the late Pleistocene. The most recent haplotypes are ecologically and biogeographically much more differentiated than the old lineages, presumably because of their persistent differentiation after the allopatric speciation events. Our study shows that in microbial eukaryotes, species diversification via the colonization of new geographical regions or ecological resources occurs much more readily than was previously thought. Consequently, divergence times of microorganisms in some lineages may be equivalent to the estimated times of speciation in plants and animals.

Anthropocene refugia: integrating history and predictive modelling to assess the space available for biodiversity in a human-dominated world

During periods of strong environmental change, some areas may serve as refugia, where components of biodiversity can find protection, persist and potentially expand from should conditions again become favourable. The refugia concept has previously been used in the context of climatic change, to describe climatically stable areas in which taxa survived past Quaternary glacial-interglacial oscillations, or where they might persist in the future under anthropogenic climate change. However, with the recognition that Earth has entered the Anthropocene, an era in which human activities are the dominant driving force on ecosystems, it is critical to also consider human pressures on the environment as factors limiting species distributions. Here, we present a novel concept, Anthropocene refugia, to refer to areas that provide spatial and temporal protection from human activities and that will remain suitable for a given taxonomic unit in the long-term. It integrates a deep-time perspective on species biogeography that provides information on the natural rather than current-day relictual distribution of species, with spatial information on modern and future anthropogenic threats. We define the concept and propose a methodology to effectively identify and map realized and potential current and future refugia, using examples for two megafaunal species as a proof of concept. We argue that identifying Anthropocene refugia will improve biodiversity conservation and restoration by allowing better prediction of key areas for conservation and potential for re-expansions today and in the future. More generally, it forms a new conceptual framework to assess and manage the impact of anthropogenic activities on past, current and future patterns of species distributions. This article is part of a discussion meeting issue 'The past is a foreign country: how much can the fossil record actually inform conservation?'

Off-target capture data, endosymbiont genes and morphology reveal a relict lineage that is sister to all other singing cicadas

Phylogenetic asymmetry is common throughout the tree of life and results from contrasting patterns of speciation and extinction in the paired descendant lineages of ancestral nodes. On the depauperate side of a node, we find extant ‘relict’ taxa that sit atop long, unbranched lineages. Here, we show that a tiny, pale green, inconspicuous and poorly known cicada in the genus Derotettix, endemic to degraded salt-plain habitats in arid regions of central Argentina, is a relict lineage that is sister to all other modern cicadas. Nuclear and mitochondrial phylogenies of cicadas inferred from probe-based genomic hybrid capture data of both target and non-target loci and a morphological cladogram support this hypothesis. We strengthen this conclusion with genomic data from one of the cicada nutritional bacterial endosymbionts, Sulcia, an ancient and obligate endosymbiont of the larger plant-sucking bugs (Auchenorrhyncha) and an important source of maternally inherited phylogenetic data. We establish Derotettiginae subfam. nov. as a new, monogeneric, fifth cicada subfamily, and compile existing and new data on the distribution, ecology and diet of Derotettix. Our consideration of the palaeoenvironmental literature and host-plant phylogenetics allows us to predict what might have led to the relict status of Derotettix over 100 Myr of habitat change in South America.