skip to main content

Title: Puma genomes from North and South America provide insights into the genomic consequences of inbreeding

Pumas are the most widely distributed felid in the Western Hemisphere. Increasingly, however, human persecution and habitat loss are isolating puma populations. To explore the genomic consequences of this isolation, we assemble a draft puma genome and a geographically broad panel of resequenced individuals. We estimate that the lineage leading to present-day North American pumas diverged from South American lineages 300–100 thousand years ago. We find signatures of close inbreeding in geographically isolated North American populations, but also that tracts of homozygosity are rarely shared among these populations, suggesting that assisted gene flow would restore local genetic diversity. The genome of a Florida panther descended from translocated Central American individuals has long tracts of homozygosity despite recent outbreeding. This suggests that while translocations may introduce diversity, sustaining diversity in small and isolated populations will require either repeated translocations or restoration of landscape connectivity. Our approach provides a framework for genome-wide analyses that can be applied to the management of similarly small and isolated populations.

more » « less
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; « less
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Urbanization is decreasing wildlife habitat and connectivity worldwide, including for apex predators, such as the puma (Puma concolor). Puma populations along California's central and southern coastal habitats have experienced rapid fragmentation from development, leading to calls for demographic and genetic management. To address urgent conservation genomic concerns, we used double‐digest restriction‐site associated DNA (ddRAD) sequencing to analyze 16,285 genome‐wide single‐nucleotide polymorphisms (SNPs) from 401 pumas sampled broadly across the state. Our analyses indicated support for 4–10 geographically nested, broad‐ to fine‐scale genetic clusters. At the broadest scale, the four genetic clusters had high genetic diversity and exhibited low linkage disequilibrium, indicating that pumas have retained genomic diversity statewide. However, multiple lines of evidence indicated substructure, including 10 finer‐scale genetic clusters, some of which exhibited fixed alleles and linkage disequilibrium. Fragmented populations along the Southern Coast and Central Coast had particularly low genetic diversity and strong linkage disequilibrium, indicating genetic drift and close inbreeding. Our results demonstrate that genetically at risk populations are typically nested within a broader‐scale group of interconnected populations that collectively retain high genetic diversity and heterogenous fixations. Thus, extant variation at the broader scale has potential to restore diversity to local populations if management actions can enhance vital gene flow and recombine locally sequestered genetic diversity. These state‐ and genome‐wide results are critically important for science‐based conservation and management practices. Our nested population genomic analysis highlights the information that can be gained from population genomic studies aiming to provide guidance for the conservation of fragmented populations.

    more » « less
  2. Enard, David (Ed.)
    Abstract High-quality reference genomes are fundamental tools for understanding population history, and can provide estimates of genetic and demographic parameters relevant to the conservation of biodiversity. The federally endangered Pacific pocket mouse (PPM), which persists in three small, isolated populations in southern California, is a promising model for studying how demographic history shapes genetic diversity, and how diversity in turn may influence extinction risk. To facilitate these studies in PPM, we combined PacBio HiFi long reads with Omni-C and Hi-C data to generate a de novo genome assembly, and annotated the genome using RNAseq. The assembly comprised 28 chromosome-length scaffolds (N50 = 72.6 MB) and the complete mitochondrial genome, and included a long heterochromatic region on chromosome 18 not represented in the previously available short-read assembly. Heterozygosity was highly variable across the genome of the reference individual, with 18% of windows falling in runs of homozygosity (ROH) >1 MB, and nearly 9% in tracts spanning >5 MB. Yet outside of ROH, heterozygosity was relatively high (0.0027), and historical Ne estimates were large. These patterns of genetic variation suggest recent inbreeding in a formerly large population. Currently the most contiguous assembly for a heteromyid rodent, this reference genome provides insight into the past and recent demographic history of the population, and will be a critical tool for management and future studies of outbreeding depression, inbreeding depression, and genetic load. 
    more » « less
  3. Abstract Background

    Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada.


    Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006–0.0009; θW = 0.0005–0.0007) relative to populations in California (π = 0.0014–0.0019; θW = 0.0011–0.0017) and the Rocky Mountains (π = 0.0025–0.0027; θW = 0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’sDwas positive for all sites (D = 0.240–0.811), consistent with recent contraction in population sizes range-wide.


    Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.

    more » « less
  4. Abstract Aim

    Numerous glacial refugia have been hypothesized along North America's North Pacific Coast that may have increased divergence of refugial taxa, leading to elevated endemism and subsequently clustered hybrid zones following deglaciation. The locations and community composition of these ice‐free areas remains controversial, but whole‐genome sequences now enable detailed analysis of the demographic and evolutionary histories of refugial taxa. Here, we use genomic data to test spatial and temporal processes of diversification among martens with respect to the Coastal Refugium Hypothesis, to understand the role of climate cycling in shaping diversity across complex landscapes.


    North America and North Pacific Coast archipelagos.


    North American martens (Martes).


    Short‐read whole‐genome resequencing data were generated for 11 martens: fourM. americana, fourM. caurina, two hybrids, and one outgroup (Martes zibellina). Sampling was representative of known genetic clades within New World martens, including sampling within insular and continental hybrid zones and along the North Pacific Coast (five island populations).ADMIXTURE, F‐statistics, andD‐statistics (ABBA‐BABA) were used to identify introgression and infer directionality. Heterozygosity densities, estimated via PSMC, were used to characterize historical demography at and below the species level to infer refugial and colonization processes.


    Forest‐associated Pacific martens (M. caurina) are divided into distinct insular and continental clades consistent with the Coastal Refugium Hypothesis. There was no evidence of introgression on islands that received historical translocations of American pine martens (M. americana), but introgression was detected in two active zones of secondary contact: one insular and one continental. Only early‐generational hybrids were identified across multiple hybrid zones, a pattern consistent with potential genetic swamping ofM. caurinabyM. americana.

    Main conclusions

    Despite an incomplete fossil record, genomic evidence supports the persistence of forest‐associated martens, likely the insular Pacific marten lineage, along the western edges of the Alexander Archipelago during the Last Glacial Maximum. This discovery informs our understanding of refugial paleoenvironments, critical to interpreting refugial timing, duration, and community composition. Genomic reevaluations of other taxa along North America's North Pacific Coast may yield new and deeper perspectives on the history of refugial forest communities and the role of dynamic climate shifts in shaping high‐latitude diversity across complex insular landscapes.

    more » « less
  5. Barton, Nick H. (Ed.)

    In the formation of species, adaptation by natural selection generates distinct combinations of traits that function well together. The maintenance of adaptive trait combinations in the face of gene flow depends on the strength and nature of selection acting on the underlying genetic loci. Floral pollination syndromes exemplify the evolution of trait combinations adaptive for particular pollinators. The North American wildflower genusPenstemondisplays remarkable floral syndrome convergence, with at least 20 separate lineages that have evolved from ancestral bee pollination syndrome (wide blue-purple flowers that present a landing platform for bees and small amounts of nectar) to hummingbird pollination syndrome (bright red narrowly tubular flowers offering copious nectar). Related taxa that differ in floral syndrome offer an attractive opportunity to examine the genomic basis of complex trait divergence. In this study, we characterized genomic divergence among 229 individuals from aPenstemonspecies complex that includes both bee and hummingbird floral syndromes. Field plants are easily classified into species based on phenotypic differences and hybrids displaying intermediate floral syndromes are rare. Despite unambiguous phenotypic differences, genome-wide differentiation between species is minimal. Hummingbird-adapted populations are more genetically similar to nearby bee-adapted populations than to geographically distant hummingbird-adapted populations, in terms of genome-widedXY. However, a small number of genetic loci are strongly differentiated between species. These approximately 20 “species-diagnostic loci,” which appear to have nearly fixed differences between pollination syndromes, are sprinkled throughout the genome in high recombination regions. Several map closely to previously established floral trait quantitative trait loci (QTLs). The striking difference between the diagnostic loci and the genome as whole suggests strong selection to maintain distinct combinations of traits, but with sufficient gene flow to homogenize the genomic background. A surprisingly small number of alleles confer phenotypic differences that form the basis of species identity in this species complex.

    more » « less