Environmental and geographic drivers of global bat phylogenetic diversity

This global study demonstrates that publicly available single-locus genetic data can effectively map bat phylogenetic diversity and reveals that historical and current temperature variables are the primary drivers of these patterns, with Southeast Asia and South America identified as key hotspots for conservation.

The Gist

Imagine the entire world's bat population as a massive, ancient library. For a long time, conservationists have tried to save this library by counting how many different titles (species) are on the shelves. But this paper argues that counting titles isn't enough. We also need to look at the history of the books.

Some books are brand new editions of the same story, while others are ancient, unique manuscripts that have no close relatives. If you lose a unique ancient manuscript, you lose a huge chunk of the library's history. If you lose a modern reprint, you lose less "history," even if the title count goes down by one.

This paper is like a team of detectives using a giant, global map and a time machine to figure out where the most unique, ancient bat history is hiding and why it's there.

Here is the breakdown of their investigation:

1. The Detective Work: Using "Genetic Fingerprints"

Instead of relying on old field guides (which often have wrong names or miss hidden species), the researchers went straight to the source: DNA.

• The Analogy: Imagine trying to sort a pile of mixed-up LEGO sets. Instead of looking at the box art (which might be wrong), they looked at the specific shapes of the bricks (the DNA).
• They gathered over 14,000 DNA samples from bats all over the world. They used a computer program to group these samples into "genetic families" (called OTUs). This helped them find "cryptic species"—bats that look identical but are actually genetically different, like twins who have been separated at birth.

2. The Map: Dividing the World into "Habitat Neighborhoods"

They didn't just look at random squares on a map. They divided the world into ecoregions—natural neighborhoods defined by their climate and plants (like "The Amazon Rainforest" or "The Himalayan Mountains").

• The Analogy: Think of the world as a giant patchwork quilt. Each patch is a different ecosystem. The researchers asked: "Which patches of the quilt hold the most unique, ancient bat stories?"

3. The Time Machine: What the Past Tells Us

The researchers didn't just look at today's weather. They used a "time machine" to look at what the climate was like during the Last Glacial Maximum (when ice sheets covered the Earth), the Last Interglacial (a warm period before that), and even the Pliocene (millions of years ago).

• The Analogy: To understand why a family lives in a specific town today, you have to know where their grandparents lived and what the weather was like when they were young. The bats' current locations are heavily influenced by where they could survive during ice ages and warm periods in the past.

4. The Big Discovery: Temperature is the Master Key

Using a powerful computer model (called Random Forest, which is like a super-smart guesser that learns from thousands of examples), they tried to predict where the most unique bat history would be found.

The Results:

• The Main Driver: Temperature is the boss. Specifically, how warm it is now, how warm it was during the ice ages, and how much the temperature has changed over time.
• The Hotspots: The places with the most unique bat history are Southeast Asia, South America, parts of Africa, and the Himalayas.
• The Surprise: While South America has a huge number of bat species, Southeast Asia and the Himalayas actually have a higher concentration of unique evolutionary history. It's like South America has a library with 1,000 copies of the same bestseller, while the Himalayas have 200 books, but every single one is a rare, one-of-a-kind original.

5. Why Does This Matter? (The "So What?")

The world is changing fast. Climate change is heating up our planet, and human development is shrinking habitats.

• The Problem: If we only protect areas with the most species, we might miss the areas with the most unique evolutionary history. If we lose a unique bat lineage in the Himalayas, we lose millions of years of evolution that can never be replaced.
• The Solution: This study gives conservationists a new tool. Instead of just saying "Save the place with the most bats," they can now say, "Save the place with the most unique bat history."
• The Takeaway: The study shows that publicly available DNA data is a powerful, cheap, and fast way to map out the world's evolutionary treasures. It proves that to save the future of biodiversity, we need to understand the past.

In a Nutshell

Think of bats as the living fossils of the sky. This paper used a giant genetic puzzle to show us that the most precious parts of the bat "family tree" are hanging out in the warm, stable, and historically complex corners of the world. By protecting these specific neighborhoods, we aren't just saving animals; we are saving the deep, unique history of life on Earth.

Technical Summary

1. Problem Statement

Understanding the drivers of biodiversity is critical for effective conservation, particularly in the face of anthropogenic habitat fragmentation and climate change. While species richness is a common metric, Phylogenetic Diversity (PD) offers a more robust measure by accounting for evolutionary history and genetic distinctness. However, global analyses of PD are often hindered by:

• Taxonomic instability: Rapidly changing species names and misidentifications in open-source databases.
• Data limitations: A lack of globally consistent, georeferenced genetic data.
• Scale dependency: Uncertainty regarding how environmental drivers (e.g., temperature, precipitation, human impact) influence PD at different spatial scales (local ecoregions vs. broad biomes).
• Historical context: A frequent omission of paleoclimatic data (e.g., Last Glacial Maximum) in predicting current biodiversity patterns.

The study aims to overcome these barriers by utilizing open-source single-locus genetic data to map global bat PD and identify the environmental and geographic predictors driving these patterns across multiple spatial scales.

2. Methodology

Data Curation and Processing

• Source: The authors curated a global dataset from GenBank and GBIF via the phylogatR platform, focusing on the mitochondrial Cytochrome oxidase I (COI) gene.
• Volume: The final dataset comprised 14,037 COI sequences representing 343 described species across 17 of 21 bat families.
• Species Delimitation: To bypass taxonomic inconsistencies, the authors did not rely on species names. Instead, they used Assemble Species by Automatic Partitioning (ASAP), a distance-based method, to define Operational Taxonomic Units (OTUs). This resulted in 372 genetic OTUs.
• Phylogeny Reconstruction: A Maximum Likelihood (ML) tree was reconstructed using RAxML-NG with the GTR+F+R8 model. Divergence times were estimated using TreePL (penalized likelihood) calibrated with fossil and molecular clock constraints from literature.

Spatial Units and Diversity Metrics

• Spatial Scales: Two datasets were analyzed:
  1. 109-eco set: Individual RESOLVE Ecoregions (smaller scale).
  2. 25-eco set: Aggregated continuous ecoregions grouped by biome (larger scale).
• Metrics: Three PD metrics were calculated using the R package picante:
  • Faith's Index: Total branch length (highly correlated with OTU richness, thus excluded from predictive modeling).
  • Mean Pairwise Distance (MPD): Average evolutionary distance between all pairs.
  • Mean Nearest Taxon Distance (MNTD): Average distance to the closest relative.
  • Note: Standardized Effect Sizes (SES) were used to control for missing taxa.

Predictive Modeling

• Predictors: 310 environmental and geographic variables were curated, including:
  • Current & Historical Climate: BIOCLIM layers for Present, Mid-Holocene, Last Glacial Maximum (LGM), Last Interglacial (LIG), and Pliocene.
  • Climate Change: Differences in median values between time periods (e.g., LGM to LIG).
  • Geography: Latitude, longitude, elevation, population density, flood risk, and habitat heterogeneity.
• Algorithm: Random Forest (RF) regression was used to predict MPD and MNTD.
  • Variables were grouped into 7 categories (Current, Geography, Holocene, LGM, LIG, Pliocene, Climate Change).
  • Correlated variables ($r \geq 0.75$) were removed to reduce multicollinearity.
  • Models were run 100 times to assess stability.
  • Generalized Least Squares (GLS) regressions were subsequently applied to the top predictors to test statistical significance while accounting for spatial autocorrelation.

3. Key Contributions

• Methodological Innovation: Demonstrated that single-locus (COI) open-source data, when processed through automated species delimitation (ASAP), can effectively reconstruct global phylogenetic patterns and identify diversity drivers, even with incomplete taxonomic sampling.
• Scale-Dependent Drivers: Highlighted that the drivers of PD vary significantly by spatial scale. Global patterns are driven by historical climate, while local patterns are influenced by contemporary human factors.
• Historical Climate Integration: Successfully integrated deep-time paleoclimatic data (Pliocene to LIG) to show that historical temperature stability and fluctuations are primary drivers of current bat PD.

4. Key Results

Phylogenetic Patterns

• Hotspots: The highest levels of PD were found in Southeast Asia, South America, parts of Africa, and the Himalayas.
• Lineage Through Time: The LTT plot confirmed a slowdown in speciation rates over time for most bat families, with the exception of the Phyllostomidae (New World leaf-nosed bats), which showed a more recent rapid diversification.
• OTU vs. Species: The ASAP method identified cryptic diversity, splitting 81 nominal species and merging 161 nominal species into 69 groups, revealing hidden biodiversity.

Predictive Modeling Results

• Model Performance: RF models using MPD as the response variable had high predictive power (explaining 33–68% of variation), whereas MNTD models had low predictive power.
• 25-Ecoregion Set (Global Scale):
  • Top Drivers: Temperature variables with a historical component were dominant.
  • Significant Predictors: Current median temperature of the wettest quarter, LGM median temperature of the wettest quarter, and the temperature change between the LIG and LGM.
  • Interpretation: Areas with historically stable, warm climates harbor higher evolutionary distinctness.
• 109-Ecoregion Set (Local Scale):
  • Top Drivers: While temperature remained important, latitude and human population density emerged as top predictors (though only temperature remained significant after spatial autocorrelation correction).
  • Interpretation: At finer scales, human impact and latitudinal gradients begin to influence PD patterns.
• Reduced Dataset (41 ecoregions with significant PD estimates):
  • Top predictors included Köppen-Geiger climate classification heterogeneity, Pliocene temperature variability, and LIG coldest quarter temperatures.

5. Significance and Implications

• Conservation Strategy: The study identifies specific ecoregions (e.g., Southeast Asia, Himalayas, parts of Africa) that are critical for preserving evolutionary history, not just species counts. These areas should be prioritized for protection.
• Climate Resilience: The strong association between PD and historical temperature stability suggests that regions with stable climates have acted as evolutionary cradles. Protecting these areas is vital for maintaining biodiversity resilience against future climate change.
• Efficiency of Single-Locus Data: The research validates that computationally efficient, single-locus genetic data can serve as a proxy for complex biodiversity assessments, making global conservation planning more accessible and scalable.
• Human Impact: The identification of population density as a predictor at smaller scales underscores the immediate threat anthropogenic pressure poses to phylogenetic diversity, necessitating land-use management strategies that integrate evolutionary history.

In conclusion, this paper provides a robust framework for using open-source genetic data to map and predict global biodiversity patterns, emphasizing that historical temperature stability is the primary engine of bat phylogenetic diversity, while human activity increasingly shapes these patterns at local scales.