Whole genomes reveal how Andean climate history shapes genetic diversity and modern conservation risk in South American pumas

By integrating whole-genome sequencing with paleoclimate history, this study reveals how Andean climatic oscillations shaped distinct demographic trajectories and genetic risks in Ecuadorian pumas, providing a framework for targeted conservation strategies that balance corridor restoration with genetic rescue.

The Gist

Imagine the Andes Mountains as a giant, ancient stage where the climate has been acting like a dramatic director for millions of years. Sometimes the stage is warm and open, allowing actors to move freely. Other times, it's cold and icy, forcing them into separate, isolated dressing rooms.

This paper is the story of the Puma (the mountain lion), the star actor of this stage in Ecuador. Scientists used a high-tech "time machine" called whole-genome sequencing to read the pumas' DNA. Think of DNA not just as a code, but as a family photo album that has been passed down for thousands of years, recording every move, every separation, and every family reunion.

Here is the story of the Ecuadorian pumas, broken down into simple chapters:

1. The Three Neighborhoods

The researchers looked at pumas living in three very different "neighborhoods" in Ecuador:

• The Misty Northwest (Mashpi): A lush, wet cloud forest on the western side of the mountains.
• The High Northern Andes (Cayambe-Coca): A vast, rugged area high up in the mountains.
• The Dry Southern Coast (Arenillas): A hot, dry region on the Pacific coast.

2. The Climate Rollercoaster

For millions of years, the Earth went through "ice ages" (cold) and "warm periods" (hot).

• During the Cold: The mountains acted like giant walls. The pumas in the north and the south got stuck in their own corners, unable to meet.
• During the Warm: The ice melted, and the "corridors" opened up. It was like a giant highway opening between the north and south, allowing pumas to travel, meet, and have babies with each other.

3. The Family Secrets Revealed

By reading the DNA "photo albums," the scientists found three distinct stories:

Story A: The Hermit of the Mist (Mashpi)
The pumas in the northwest cloud forest are like a reclusive family that has lived in the same house for 3.5 million years. They never really left. Even though they are only a short drive away from the northern pumas, their DNA shows they have been isolated for a very long time. They are a unique, distinct lineage, perfectly adapted to their wet, misty home.

Story B: The Travelers of the North and South (Cayambe-Coca & Arenillas)
The pumas in the north and the south are like long-lost cousins. Even though they are currently 470 kilometers apart (a long way!), their DNA shows they used to mix freely when the climate was warm. They share a common history of meeting up during the "warm intermissions" of the ice ages. However, recently, human activity (roads, farms, cities) has built a wall between them, cutting off their ancient connection.

Story C: The Struggling Southern Group (Arenillas)
The southern pumas are in trouble. Imagine a small town that has been shrinking for thousands of years. Because they have always lived in a small, dry area, their population has been small for a long time.

• The Problem: In small populations, bad genetic "typos" (mutations) can pile up because there aren't enough new people to dilute them.
• The Evidence: These pumas have a lot of "homozygosity," which is a fancy way of saying their DNA is very repetitive, like a song on repeat. This means they have a higher burden of harmful genetic errors that could make them sick or less able to survive.

4. The "Kinked Tail" Mystery

One puma in the north, named Julio, has a very specific problem: his tail is bent or "kinked."

• The Clue: In the animal world, a kinked tail is often a sign of inbreeding (when close relatives have babies).
• The Twist: Julio's DNA doesn't show that his entire family is inbred. It looks like he might be the result of a very recent, local accident—perhaps his parents were related, but the rest of the population is fine. It's like finding one cracked egg in a basket of fresh ones. It suggests that while the northern population is generally healthy, there might be small pockets where pumas are getting stuck and mating with relatives.

5. What Should We Do? (The Conservation Plan)

The scientists are saying: "We need to fix the broken bridges."

• For the South (Arenillas): They are like a ship taking on water. They need help from the outside. The scientists suggest genetic rescue—bringing in healthy pumas from the north to mix with the southern group. This would introduce fresh DNA, wash out the bad "typos," and give them a fighting chance.
• For the North (Mashpi): They are unique and special. We shouldn't mix them with the others too much, or we might lose their special "flavor" (local adaptation). Instead, we need to protect their forest and maybe build a bridge to other similar wet forests so they can stay safe without losing their identity.
• For Everyone: We need to restore corridors. Imagine building wildlife overpasses or protecting strips of forest so pumas can walk from the north to the south again. This would let the "travelers" meet up naturally, keeping the gene pool healthy and diverse.

The Bottom Line

The Andes pumas are living history books. Their DNA tells us that climate change shaped them for millions of years, but human change is now shaping their future. If we don't fix the broken connections between their neighborhoods, the southern pumas might fade away, and the northern ones might start making mistakes. But if we restore the paths, these magnificent cats can continue their ancient dance across the mountains.

Technical Summary

1. Problem Statement

While the puma (Puma concolor) is widely distributed across the Americas, the long-term genomic consequences of Andean climatic oscillations on its population structure and contemporary conservation status remain unclear. Previous studies relied heavily on mitochondrial DNA and microsatellites, which lack the resolution to reconstruct detailed demographic histories or assess functional genetic load (deleterious mutations). Furthermore, it is unknown how historical climate-driven isolation versus connectivity has shaped current genetic diversity, inbreeding risks, and the accumulation of harmful variants in Ecuadorian puma populations.

2. Methodology

The study employed a comprehensive whole-genome sequencing (WGS) approach to analyze six pumas from three ecologically distinct regions in Ecuador:

• Sampling:
  • Mashpi: Northwestern Andean foothills (Cloud forest/Chocó region).
  • Cayambe-Coca: Northern Andes (Heterogeneous highland landscape).
  • Arenillas: Southern Pacific coast (Tumbesian dry forest).
• Sequencing & Assembly:
  • Generated 150 bp paired-end reads on an Illumina NovaSeq X 10B platform.
  • Mapped reads to a high-quality, chromosome-level Puma concolor reference genome (UCLA assembly, ~2.5 Gb, 46× coverage).
  • Used GATK Best Practices for variant calling, filtering for high-confidence SNPs (min depth 5, GQ ≥ 20).
• Analytical Framework:
  • Population Structure: Principal Component Analysis (PCA), kinship analysis (PLINK), and Bayesian clustering (fastStructure).
  • Phylogeography & Gene Flow: TreeMix to infer population splits and historical migration edges.
  • Demographic History: Multiple Sequentially Markovian Coalescent (MSMC) to estimate effective population size ($N_e$) trajectories over the last 500,000 years.
  • Inbreeding & Genetic Load:
    • Runs of Homozygosity (ROH): Categorized into short (<1 Mb), medium (1–10 Mb), and long (>10 Mb) to distinguish ancient bottlenecks from recent inbreeding.
    • Deleterious Variation: Functional annotation (SnpEff) to quantify Loss-of-Function (LOF) and missense mutations.
  • Paleoclimatic Reconstruction: Integration of PaleoClim data (Mid-Pliocene, Last Interglacial, Heinrich Stadial 1) with PCA to model historical climatic distances between populations.

3. Key Results

A. Population Structure and Historical Connectivity

• Distinct Lineages: PCA and TreeMix revealed three distinct genetic clusters. The Mashpi population is a deeply divergent, non-admixed lineage, isolated for millions of years.
• Historical Admixture: Despite current geographic separation (~470 km), the Arenillas (south) and Cayambe-Coca (north) populations show signals of historical gene flow. TreeMix detected a migration edge between them, suggesting they were connected during warm interglacial periods.
• Recent Fragmentation: Genetic clustering (fastStructure) at $K=5$ suggests recent subdivision within the northern populations, likely driven by contemporary habitat fragmentation.

B. Demographic History and Climate Correlation

• Mashpi (Northwest): Remained in a stable, humid refugium (Chocó biogeographic region) with persistent climatic isolation from other populations throughout the Pleistocene.
• Arenillas vs. Cayambe-Coca: Paleoclimatic modeling shows these regions experienced climatic convergence during warm periods (Mid-Pliocene ~3.2 Ma and Last Interglacial ~130 kya), facilitating dispersal. Conversely, during cold phases (e.g., Heinrich Stadial 1, 17–14.7 kya), they experienced divergent conditions (cold/high-elevation vs. warm/dry), leading to isolation.
• Population Size Trajectories:
  • Arenillas: Exhibited a chronically small effective population size ($N_e \approx 80,000$ historically, dropping to ~10,000 recently), consistent with long-term demographic constraint.
  • Northern Populations: Maintained larger $N_e$ (~100,000) until a sharp decline in the last 20,000 years, dropping to ~20,000.

C. Genetic Load and Inbreeding

• Arenillas (Southern): Showed the highest burden of homozygous Loss-of-Function (LOF) variants (~25% of alleles). This is attributed to long-term small $N_e$ reducing the efficiency of purifying selection, allowing deleterious alleles to drift to fixation.
• Cayambe-Coca (Northern): One individual, "Julio," displayed:
  • The highest burden of long ROH (~47 Mb, >2% of genome).
  • A kinked tail phenotype, a known marker of inbreeding in pumas.
  • Elevated missense mutations.
  • Interpretation: While the rest of the population shows healthy diversity, Julio's profile suggests recent, localized inbreeding or a specific locus effect, indicating early signs of inbreeding depression in this region.
• Mashpi: Showed high heterozygosity and low ROH, consistent with a stable, isolated population that has not suffered recent severe bottlenecks.

4. Key Contributions

1. First WGS Analysis of Ecuadorian Pumas: Moves beyond mtDNA/microsatellites to provide a high-resolution view of genomic diversity, demographic history, and functional load.
2. Climate-Genome Linkage: Directly correlates paleoclimatic oscillations (warm vs. cold periods) with genomic signatures of connectivity and isolation, demonstrating that Andean puma structure is a product of deep-time climate dynamics rather than just recent human activity.
3. Differentiation of Inbreeding Types: Distinguishes between ancient demographic constraint (Arenillas: small $N_e$, high LOF, medium ROH) and recent inbreeding (Julio: high long ROH, phenotypic defects), offering nuanced conservation diagnostics.
4. Identification of a Cryptic Lineage: Confirms the Mashpi population as a distinct, evolutionarily significant unit (ESU) that has been isolated for ~3.5 million years.

5. Significance and Conservation Implications

The study provides a framework for region-specific conservation strategies that balance connectivity with the preservation of local adaptation:

• Arenillas (South): Due to long-term isolation and high genetic load, this population is vulnerable. Genetic rescue (assisted gene flow) from genetically compatible northern populations may be necessary to reduce deleterious load, provided it does not swamp local adaptations.
• Mashpi (Northwest): As a distinct, locally adapted lineage, this population should be protected from admixture with other populations to preserve its unique evolutionary trajectory. Connectivity should be restored only with similar humid forest habitats (e.g., Chocó region).
• Cayambe-Coca (North): While generally robust, the detection of inbreeding in individual "Julio" signals that habitat fragmentation is disrupting historical connectivity. Restoring corridors in the western foothills is critical to prevent localized inbreeding depression.
• General: The findings highlight that "Least Concern" status on the IUCN Red List may mask severe, localized genomic erosion driven by both deep-time climate history and modern anthropogenic fragmentation.

Conclusion: The genomic health of Andean pumas is a legacy of Pleistocene climate cycles. Effective conservation requires maintaining historical connectivity corridors where appropriate while protecting unique, isolated lineages from genetic swamping.