Lineage range estimation method reveals fine-scale endemism linked to Pleistocene stability in Australian rainforest herpetofauna

Disentangling the drivers of continental mammalian endemism

Endemic species and species with small ranges are ecologically and evolutionarily distinct and are vulnerable to extinction. Determining which abiotic and biotic factors structure patterns of endemism on continents can advance our understanding of global biogeographic processes, but spatial patterns of mammalian endemism have not yet been effectively predicted and reconstructed. Using novel null model techniques, we reconstruct trends in mammalian endemism and describe the isolated and combined effects of physiographic, ecological, and evolutionary factors on endemism. We calculated weighted endemism for global continental ecoregions and compared the spatial distribution of endemism to niche-based, geographic null models of endemism. These null models distribute species randomly across continents, simulating their range sizes from their degree of climatic specialization. They isolate the effects of physiography (topography and climate) and species richness on endemism. We then ran linear and structural models to determine how topography and historical climate stability influence endemism. The highest rates of mammalian endemism were found in topographically rough, climatically stable ecoregions with many species. The null model that isolated physiography did not closely approximate the observed distribution of endemism (r²  = .09), whereas the null model that incorporated both physiography and species richness did (r²  = .59). The linear models demonstrate that topography and climatic stability both influenced endemism values, but that average climatic niche breadth was not highly correlated with endemism. Climate stability and topography both influence weighted endemism in mammals, but the spatial distribution of mammalian endemism is driven by a combination of physiography and species richness. Despite its relationship to individual range size, average climate niche breadth has only a weak influence on endemism. The results highlight the importance of historical biogeographic processes (e.g. centers of speciation) and geography in driving endemism patterns, and disentangle the mechanisms structuring species ranges worldwide.

Molecular ecology meets systematic conservation planning

Integrative and proactive conservation approaches are critical to the long-term persistence of biodiversity. Molecular data can provide important information on evolutionary processes necessary for conserving multiple levels of biodiversity (genes, populations, species, and ecosystems). However, molecular data are rarely used to guide spatial conservation decision-making. Here, we bridge the fields of molecular ecology (ME) and systematic conservation planning (SCP) (the 'why') to build a foundation for the inclusion of molecular data into spatial conservation planning tools (the 'how'), and provide a practical guide for implementing this integrative approach for both conservation planners and molecular ecologists. The proposed framework enhances interdisciplinary capacity, which is crucial to achieving the ambitious global conservation goals envisioned for the next decade.

Environmental and topographic drivers of amphibian phylogenetic diversity and endemism in the Iberian Peninsula

Understanding the ecological and evolutionary processes driving biodiversity patterns and allowing their persistence is of utmost importance. Many hypotheses have been proposed to explain spatial diversity patterns, including water-energy availability, habitat heterogeneity, and historical climatic refugia. The main goal of this study is to identify if general spatial drivers of species diversity patterns of phylogenetic diversity (PD) and phylogenetic endemism (PE) at the global scale are also predictive of PD and PE at regional scales, using Iberian amphibians as a case study. Our main hypothesis assumes that topography along with contemporary and historical climate are drivers of phylogenetic diversity and endemism, but that the strength of these predictors may be weaker at the regional scale than it tends to be at the global scale. We mapped spatial patterns of Iberian amphibians' phylogenetic diversity and endemism, using previously published phylogenetic and distribution data. Furthermore, we compiled spatial data on topographic and climatic variables related to the water-energy availability, topography, and historical climatic instability hypotheses. To test our hypotheses, we used Spatial Autoregressive Models and selected the best model to explain diversity patterns based on Akaike Information Criterion. Our results show that, out of the variables tested in our study, water-energy availability and historical climate instability are the most important drivers of amphibian diversity in Iberia. However, as predicted, the strength of these predictors in our case study is weaker than it tends to be at global scales. Thus, additional drivers should also be investigated and we suggest caution when interpreting these predictors as surrogates for different components of diversity.

Congruent evolutionary responses of European steppe biota to late Quaternary climate change

Quaternary climatic oscillations had a large impact on European biogeography. Alternation of cold and warm stages caused recurrent glaciations, massive vegetation shifts, and large-scale range alterations in many species. The Eurasian steppe biome and its grasslands are a noteworthy example; they underwent climate-driven, large-scale contractions during warm stages and expansions during cold stages. Here, we evaluate the impact of these range alterations on the late Quaternary demography of several phylogenetically distant plant and insect species, typical of the Eurasian steppes. We compare three explicit demographic hypotheses by applying an approach combining convolutional neural networks with approximate Bayesian computation. We identified congruent demographic responses of cold stage expansion and warm stage contraction across all species, but also species-specific effects. The demographic history of the Eurasian steppe biota reflects major paleoecological turning points in the late Quaternary and emphasizes the role of climate as a driving force underlying patterns of genetic variance on the biome level.

Ecological niche modelling and climate change in two species groups of huntsman spider genus Eusparassus in the Western Palearctic

The huntsman spiders' genus Eusparassus are apex arthropod predators in desert ecosystems of the Afrotropical and Palearctic ecoregions. The Eusparassus dufouri and E. walckenaeri clades are two distinct taxonomic, phylogenetic, and geographic units concerning morphology, molecular phylogeny, and spatial data; but little is known about their ecological niche. We applied the maximum-entropy approach and modelled ecologic niches of these two phylogenetically closely related clades. Ecological niches of the two clades were compared using identity and background tests and two different metrics, the Schooner's D and Warren's I. We also predicted the impacts of climate change on the distribution of the two clades. The results of the identity test showed that the ecological niches of the two clades were different in geographic space but were similar in environmental space. While results of the background test revealed that the ecological niches of the two clades were similar in geographic and environmental space. This indicated that "niche conservatism" had an important role over the evolutionary time of allopatric diversification. However, the normalized difference vegetation index vs. topographic heterogeneity had influenced the niches of the dufouri and walckenaeri clades, respectively. The analyses recovered that the two clades' climatically suitable habitats will increase under future climate (the year 2070). However, since the two clades are characterized by the narrow range of environmental optimum and the accordingly high limits of tolerance, they are vulnerable to climate change.

Lineage diversity within a widespread endemic Australian skink to better inform conservation in response to regional-scale disturbance

Much attention is paid in conservation planning to the concept of a species, to ensure comparability across studies and regions when classifying taxa against criteria of endangerment and setting priorities for action. However, various jurisdictions now allow taxonomic ranks below the level of species and nontaxonomic intraspecific divisions to be factored into conservation planning—subspecies, key populations, evolutionarily significant units, or designatable units. Understanding patterns of genetic diversity and its distribution across the landscape is a key component in the identification of species boundaries and determination of substantial geographic structure within species. A total of 12,532 reliable polymorphic SNP loci were generated from 63 populations (286 individuals) covering the distribution of the Australian eastern three‐lined skink, Bassiana duperreyi, to assess genetic population structure in the form of diagnosable lineages and their distribution across the landscape, with particular reference to the recent catastrophic bushfires of eastern Australia. Five well‐supported diagnosable operational taxonomic units (OTUs) existed within B. duperreyi. Low levels of divergence of B. duperreyi between mainland Australia and Tasmania (no fixed allelic differences) support the notion of episodic exchange of alleles across Bass Strait (ca 60 m, 25 Kya) during periods of low sea level during the Upper Pleistocene rather than the much longer period of isolation (1.7 My) indicated by earlier studies using mitochondrial sequence variation. Our study provides foundational work for the detailed taxonomic re‐evaluation of this species complex and the need for biodiversity assessment to include an examination of cryptic species and/or cryptic diversity below the level of species. Such information on lineage diversity within species and its distribution in the context of disturbance at a regional scale can be factored into conservation planning regardless of whether a decision is made to formally diagnose new species taxonomically and nomenclaturally.

Addressing the global snakebite crisis with geo-spatial analyses - Recent advances and future direction

Venomous snakebite is a neglected tropical disease that annually leads to hundreds of thousands of deaths or long-term physical and mental ailments across the developing world. Insufficient data on spatial variation in snakebite risk, incidence, human vulnerability, and accessibility of medical treatment contribute substantially to ineffective on-ground management. There is an urgent need to collect data, fill knowledge gaps and address on-ground management problems. The use of novel, and transdisciplinary approaches that take advantage of recent advances in spatio-temporal models, 'big data', high performance computing, and fine-scale spatial information can add value to snakebite management by strategically improving our understanding and mitigation capacity of snakebite. We review the background and recent advances on the topic of snakebite related geospatial analyses and suggest avenues for priority research that will have practical on-ground applications for snakebite management and mitigation. These include streamlined, targeted data collection on snake distributions, snakebites, envenomings, venom composition, health infrastructure, and antivenom accessibility along with fine-scale models of spatio-temporal variation in snakebite risk and incidence, intraspecific venom variation, and environmental change modifying human exposure. These measures could improve and 'future-proof' antivenom production methods, antivenom distribution and stockpiling systems, and human-wildlife conflict management practices, while simultaneously feeding into research on venom evolution, snake taxonomy, ecology, biogeography, and conservation.

Quantitative estimates of glacial refugia for chimpanzees (Pan troglodytes) since the Last Interglacial (120,000 BP)

Paleoclimate reconstructions have enhanced our understanding of how past climates have shaped present-day biodiversity. We hypothesize that the geographic extent of Pleistocene forest refugia and suitable habitat fluctuated significantly in time during the late Quaternary for chimpanzees (Pan troglodytes). Using bioclimatic variables representing monthly temperature and precipitation estimates, past human population density data, and an extensive database of georeferenced presence points, we built a model of changing habitat suitability for chimpanzees at fine spatio-temporal scales dating back to the Last Interglacial (120,000 BP). Our models cover a spatial resolution of 0.0467° (approximately 5.19 km² grid cells) and a temporal resolution of between 1000 and 4000 years. Using our model, we mapped habitat stability over time using three approaches, comparing our modeled stability estimates to existing knowledge of Afrotropical refugia, as well as contemporary patterns of major keystone tropical food resources used by chimpanzees, figs (Moraceae), and palms (Arecacae). Results show habitat stability congruent with known glacial refugia across Africa, suggesting their extents may have been underestimated for chimpanzees, with potentially up to approximately 60,000 km² of previously unrecognized glacial refugia. The refugia we highlight coincide with higher species richness for figs and palms. Our results provide spatio-temporally explicit insights into the role of refugia across the chimpanzee range, forming the empirical foundation for developing and testing hypotheses about behavioral, ecological, and genetic diversity with additional data. This methodology can be applied to other species and geographic areas when sufficient data are available.

The roles of aridification and sea level changes in the diversification and persistence of freshwater fish lineages

While the influence of Pleistocene climatic changes on divergence and speciation has been well-documented across the globe, complex spatial interactions between hydrology and eustatics over longer timeframes may also determine species evolutionary trajectories. Within the Australian continent, glacial cycles were not associated with changes in ice cover and instead largely resulted in fluctuations from moist to arid conditions across the landscape. We investigated the role of hydrological and coastal topographic changes brought about by Plio-Pleistocene climatic changes on the biogeographic history of a small Australian freshwater fish, the southern pygmy perch Nannoperca australis. Using 7958 ddRAD-seq (double digest restriction-site associated DNA) loci and 45,104 filtered SNPs, we combined phylogenetic, coalescent and species distribution analyses to assess the various roles of aridification, sea level and tectonics and associated biogeographic changes across southeast Australia. Sea-level changes since the Pliocene and reduction or disappearance of large waterbodies throughout the Pleistocene were determining factors in strong divergence across the clade, including the initial formation and maintenance of a cryptic species, N. 'flindersi'. Isolated climatic refugia and fragmentation due to lack of connected waterways maintained the identity and divergence of inter- and intraspecific lineages. Our historical findings suggest that predicted increases in aridification and sea level due to anthropogenic climate change might result in markedly different demographic impacts, both spatially and across different landscape types.

Benchmarking Taxonomic and Genetic Diversity After the Fact: Lessons Learned From the Catastrophic 2019–2020 Australian Bushfires

Environmental catastrophes are increasing in frequency and severity under climate change, and they substantially impact biodiversity. Recovery actions after catastrophes depend on prior benchmarking of biodiversity and that in turn minimally requires critical assessment of taxonomy and species-level diversity. Long-term recovery of species also requires an understanding of within-species diversity. Australia’s 2019–2020 bushfires were unprecedented in their extent and severity and impacted large portions of habitats that are not adapted to fire. Assessments of the fires’ impacts on vertebrates identified 114 species that were a high priority for management. In response, we compiled explicit information on taxonomic diversity and genetic diversity within fire-impacted vertebrates to provide to government agencies undertaking rapid conservation assessments. Here we discuss what we learned from our effort to benchmark pre-fire taxonomic and genetic diversity after the event. We identified a significant number of candidate species (genetic units that may be undescribed species), particularly in frogs and mammals. Reptiles and mammals also had high levels of intraspecific genetic structure relevant to conservation management. The first challenge was making published genetic data fit for purpose because original publications often focussed on a different question and did not provide raw sequence read data. Gaining access to analytical files and compiling appropriate individual metadata was also time-consuming. For many species, significant unpublished data was held by researchers. Identifying which data existed was challenging. For both published and unpublished data, substantial sampling gaps prevented areas of a species’ distribution being assigned to a conservation unit. Summarising sampling gaps across species revealed that many areas were poorly sampled across taxonomic groups. To resolve these issues and prepare responses to future catastrophes, we recommend that researchers embrace open data principles including providing detailed metadata. Governments need to invest in a skilled taxonomic workforce to document and describe biodiversity before an event and to assess its impacts afterward. Natural history collections should also target increasing their DNA collections based on sampling gaps and revise their collection strategies to increasingly take population-scale DNA samples in order to document within-species genetic diversity.

Intraspecific variation and taxonomy of northern pika Ochotona hyperborea (Mammalia, Lagomorpha)

Rock-dwelling pikas of the genus Ochotona have fragmented montane geographical distribution presumably associated with geological and climatic events of Late Quaternary. Within this genus, the northern pika, O. hyperborea has the widest distribution. Taxonomic composition and subspecific structure of this species has been controversial during the past century. Sixteen available names were combined in different manners into a variable number of weakly overlapping subspecies by previous authors. A comparison of data sets from the entire distribution range of northern pikas is presented here for the first time, including: craniometric, bioacoustic, and multilocus genetic data. We examined 301 skulls, and the mitochondrial COI gene and two nuclear introns in 79 specimens, as well as the structure of alarm calls from 107 individuals. Our results show that the six subspecies within O. hyperborea correspond to the six genetic parapatric lineages. Five of these are hypothesized to be involved in hybridization at the edges of their distribution. Three acoustic races completely correspond with the six genetic lineages: each acoustic race consists of two lineages. Morphometric data do not display any phylogenetic signal in our study. The subspecies from the Mountains of Khabarovsk Territory, north of Amur River, is described de novo as O. h. fedoseevi ssp. n. The type locality of O. hyperborea is defined more accurately. A neotype is designated for O. h. cinereofusca to stabilize use of names of pikas from the Amur region.

First phylogenetic assessment and taxonomic synopsis of the open-holed trapdoor spider genus Namea (Mygalomorphae: Anamidae): a highly diverse mygalomorph lineage from Australia’s tropical eastern rainforests

The tropical and subtropical rainforests of Australia’s eastern mesic zone have given rise to a complex and highly diverse biota. Numerous old endemic, niche-conserved groups persist in the montane rainforests south of Cooktown, where concepts of serial allopatric speciation resulting from the formation of xeric interzones have largely driven our biogeographic understanding of the region. Among invertebrate taxa, studies on less vagile arachnid lineages now complement extensive research on vertebrate taxa, and phylogenetic studies on mygalomorph spiders in particular are revealing significant insights about the biogeographic history of the Australian continent since the Eocene. One mygalomorph lineage entirely endemic to Australia’s tropical and subtropical eastern rainforests is the open-holed trapdoor spider genus Namea Raven, 1984 (family Anamidae). We explore, for the first time, the phylogenetic diversity and systematics of this group of spiders, with the aims of understanding patterns of rainforest diversity in Namea, of exploring the relative roles of lineage overlap versus in situ speciation in driving predicted high levels of congeneric sympatry, and of broadly reconciling morphology with evolutionary history. Original and legacy sequences were obtained for three mtDNA and four nuDNA markers from 151 specimens, including 82 specimens of Namea. We recovered a monophyletic genus Namea sister to the genus Teyl Main, 1975, and monophyletic species clades corresponding to 30 morphospecies OTUs, including 22 OTUs nested within three main species-complex lineages. Remarkable levels of sympatry for a single genus of mygalomorph spiders were revealed in rainforest habitats, with upland subtropical rainforests in south-eastern Queensland often home to multiple (up to six) congeners of usually disparate phylogenetic affinity living in direct sympatry or close parapatry, likely the result of simultaneous allopatric speciation in already co-occurring lineages, and more recent dispersal in a minority of taxa. In situ speciation, in contrast, appears to have played a relatively minor role in generating sympatric diversity within rainforest ‘islands’. At the population level, changes in the shape and spination of the male first leg relative to evolutionary history reveal subtle but consistent interspecific morphological shifts in the context of otherwise intraspecific variation, and understanding this morphological variance provides a useful framework for future taxonomic monography. Based on the phylogenetic results, we further provide a detailed taxonomic synopsis of the genus Namea, formally diagnosing three main species-complexes (the brisbanensis-complex, the dahmsi-complex and the jimna-complex), re-illustrating males of all 15 described species, and providing images of live spiders and burrows where available. In doing so, we reveal a huge undescribed diversity of Namea species from tropical and subtropical rainforest habitats, and an old endemic fauna that is beginning to shed light on more complex patterns of rainforest biogeography.

Insights into Platypus Population Structure and History from Whole-Genome Sequencing

The platypus is an egg-laying mammal which, alongside the echidna, occupies a unique place in the mammalian phylogenetic tree. Despite widespread interest in its unusual biology, little is known about its population structure or recent evolutionary history. To provide new insights into the dispersal and demographic history of this iconic species, we sequenced the genomes of 57 platypuses from across the whole species range in eastern mainland Australia and Tasmania. Using a highly improved reference genome, we called over 6.7 M SNPs, providing an informative genetic data set for population analyses. Our results show very strong population structure in the platypus, with our sampling locations corresponding to discrete groupings between which there is no evidence for recent gene flow. Genome-wide data allowed us to establish that 28 of the 57 sampled individuals had at least a third-degree relative among other samples from the same river, often taken at different times. Taking advantage of a sampled family quartet, we estimated the de novo mutation rate in the platypus at 7.0 × 10-9/bp/generation (95% CI 4.1 × 10-9-1.2 × 10-8/bp/generation). We estimated effective population sizes of ancestral populations and haplotype sharing between current groupings, and found evidence for bottlenecks and long-term population decline in multiple regions, and early divergence between populations in different regions. This study demonstrates the power of whole-genome sequencing for studying natural populations of an evolutionarily important species.

Assessing global range expansion in a cryptic species complex: insights from the red seaweed genus Asparagopsis (Florideophyceae)

The mitochondrial genetic diversity, distribution and invasive potential of multiple cryptic operational taxonomic units (OTUs) of the red invasive seaweed Asparagopsis were assessed by studying introduced Mediterranean and Hawaiian populations. Invasive behavior of each Asparagopsis OTU was inferred from phylogeographic reconstructions, past historical demographic dynamics, recent range expansion assessments and future distributional predictions obtained from demographic models. Genealogical networks resolved Asparagopsis gametophytes and tetrasporophytes into four A. taxiformis and one A. armata cryptic OTUs. Falkenbergia isolates of A. taxiformis L3 were recovered for the first time in the western Mediterranean Sea and represent a new introduction for this area. Neutrality statistics supported past range expansion for A. taxiformis L1 and L2 in Hawaii. On the other hand, extreme geographic expansion and an increase in effective population size were found only for A. taxiformis L2 in the western Mediterranean Sea. Distribution models predicted shifts of the climatically suitable areas and population expansion for A. armata L1 and A. taxiformis L1 and L2. Our integrated study confirms a high invasive risk for A. taxiformis L1 and L2 in temperate and tropical areas. Despite the differences in predictions among modelling approaches, a number of regions were identified as zones with high invasion risk for A. taxiformis L2. Since range shifts are likely climate-driven phenomena, future invasive behavior cannot be excluded for the rest of the lineages.

Chasing the phantom: biogeography and conservation of Vipera latastei-monticola in the Maghreb (North Africa)

The Maghreb region (North Africa) constitutes a major component of the Mediterranean Basin biodiversity hotspot. During the last centuries, a consistent human population growth has led to an unprecedented rate of habitat transformation and loss in the region and thus, threatening its biodiversity. The Western Mediterranean viper Vipera latastei-monticola inhabits humid and subhumid areas in the main mountain ranges of the Maghreb, facing such threatening factors; however, its elusive character and rarity hindered data collection for distinct biological purposes. Here, we study the biogeographical patterns and conservation status of the Maghrebian V. latastei-monticola resulting from recent sampling campaigns in Morocco and Tunisia. We update species distribution, and integrate phylogeographic and ecological niche modelling analyses at both species and lineage level to identify suitable areas, and to evaluate the impact of anthropogenic transformation and level of protection of their suitable space. We identified four highly divergent mitochondrial lineages, including a new lineage endemic to the Western High Atlas, with allopatric distributions and restricted to mountain ranges, supporting the role of mountains as past climatic refugia. Despite the remoteness of suitable areas, we report widespread habitat degradation and identify the low effectiveness of the current protected areas system in preserving the species and lineages range. Our study shows the urgent need to apply management actions for the long-term conservation of this vulnerable species and suggests a revaluation of the specific status of V. monticola, as these populations likely represent an ecotype of V. latastei.

Spatial conservation prioritization of biodiversity spanning the evolutionary continuum

Accounting for evolutionary relationships between and within species is important for biodiversity conservation planning, but is rarely considered in practice. Here we introduce a novel framework to identify priority conservation areas accounting for phylogenetic and intraspecific diversity, integrating concepts from phylogeny, phylogeography, spatial statistics and spatial conservation prioritization. The framework allows planners to incorporate and combine different levels of evolutionary diversity and can be applied to any taxonomic group and to any region in the world. We illustrate our approach using amphibian and reptile species occurring in a biodiversity hotspot region, the Iberian Peninsula. We found that explicitly incorporating phylogenetic and intraspecific diversity in systematic conservation planning provides advantages in terms of maximizing overall biodiversity representation while enhancing its persistence and evolutionary potential. Our results emphasize the need to account for the evolutionary continuum in order to efficiently implement biodiversity conservation planning decisions.

Correction: Lineage Range Estimation Method Reveals Fine-Scale Endemism Linked to Pleistocene Stability in Australian Rainforest Herpetofauna

[This corrects the article DOI: 10.1371/journal.pone.0126274.].