skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: Phylogeography and population genetics of pine butterflies: Sky islands increase genetic divergence

The sky islands of southeastern Arizona (AZ) mark a major transition zone between tropical and temperate biota and are considered a neglected biodiversity hotspot. Dispersal ability and host plant specificity are thought to impact the history and diversity of insect populations across the sky islands. We aimed to investigate the population structure and phylogeography of two pine‐feeding pierid butterflies, the pine white (Neophasia menapia) and the Mexican pine white (Neophasia terlooii), restricted to these “islands” at this transition zone. Given their dependence on pines as the larval hosts, we hypothesized that habitat connectivity affects population structure and is at least in part responsible for their allopatry. We sampled DNA from freshly collected butterflies from 17 sites in the sky islands and adjacent high‐elevation habitats and sequenced these samples using ddRADSeq. Up to 15,399 SNPs were discovered and analyzed in population genetic and phylogenetic contexts with Stacks and pyRAD pipelines. Low genetic differentiation inN. menapiasuggests that it is panmictic. Conversely, there is strong evidence for population structure withinN. terlooii. Each sky island likely contains a population ofN. terlooii, and clustering is hierarchical, with populations on proximal mountains being more related to each other. TheN. menapiahabitat, which is largely contiguous, facilitates panmixia, while theN. terlooiihabitat, restricted to the higher elevations on each sky island, creates distinct population structure. Phylogenetic results corroborate those from population genetic analyses. The historical climate‐driven fluxes in forest habitat connectivity have implications for understanding the biodiversity of fragmented habitats.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Ecology and Evolution
Page Range / eLocation ID:
p. 13389-13401
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Anthropogenic alterations to landscape structure and composition can have significant impacts on biodiversity, potentially leading to species extinctions. Population‐level impacts of landscape change are mediated by animal behaviors, in particular dispersal behavior. Little is known about the dispersal habits of rails (Rallidae) due to their cryptic behavior and tendency to occupy densely vegetated habitats. The effects of landscape structure on the movement behavior of waterbirds in general are poorly studied due to their reputation for having high dispersal abilities. We used a landscape genetic approach to test hypotheses of landscape effects on dispersal behavior of the Hawaiian gallinule (Gallinula galeata sandvicensis), an endangered subspecies endemic to the Hawaiian Islands. We created a suite of alternative resistance surfaces representing biologically plausible a priori hypotheses of how gallinules might navigate the landscape matrix and ranked these surfaces by their ability to explain observed patterns in genetic distance among 12 populations on the island of O`ahu. We modeled effective distance among wetland locations on all surfaces using both cumulative least‐cost‐path and resistance‐distance approaches and evaluated relative model performance using Mantel tests, a causal modeling approach, and the mixed‐model maximum‐likelihood population‐effects framework. Across all genetic markers, simulation methods, and model comparison metrics, surfaces that treated linear water features like streams, ditches, and canals as corridors for gallinule movement outperformed all other models. This is the first landscape genetic study on the movement behavior of any waterbird species to our knowledge. Our results indicate that lotic water features, including drainage infrastructure previously thought to be of minimal habitat value, contribute to habitat connectivity in this listed subspecies.

    more » « less
  2. Abstract Aim

    Periodic lowering of sea levels and formation of land bridges can reshape phylogeographic patterns of insular biotas. Using archipelago‐wide sampling, we aimed to test if phylogeography of an old‐endemic bat lineage reflected Pleistocene land bridges.


    Solomon Islands and Papua New Guinea.


    MelonycterisandNesonycterisbats (Pteropodidae).


    We sequenced genome‐wide RADseq data for 49 specimens from 15 islands. We assessed phylogenetic relationships using maximum likelihood inRAxMLand quartet‐based methods inSVDquartets, population structure usingStructure, and admixture using maximum likelihood methods inTreeMix. We tested for genetic and geographic distance correlations using distance‐based redundancy analyses (dbRDA), identifying best‐fit models using stepwise model selection.


    Phylogenetic analyses identified fiveNesonycterisclades corresponding to Greater Bukida, Guadalcanal, Makira, Malaita and New Georgia group. Makira samples were sister to remainingNesonycteris.Structureidentified four populations: New IrelandMelonycteris melanops; andNesonycterisfrom Greater Bukida (including Guadalcanal); Malaita and Makira; and New Georgia group. Genetic backgrounds of Mono, Ngella and Guadalcanal separated from remaining Greater Bukida islands. Makira and Malaita separated into two populations. New Georgia group lacked structure, and genetic and geographic distances were not correlated. The best‐fit geographic distance models forNesonycterisand a Greater Bukida subset were least shore‐to‐shore distance; and Euclidean and least‐cost distances respectively.

    Main Conclusions

    Influences of modern and Pleistocene island isolation and connectivity were evident in the overall Phylogeography ofNesonycteris. The lack of structure or geographic distance correlations within the New Georgia group indicated all islands were interconnected during the Last Glacial Maximum or contemporary oceanic divides are ineffective barriers. Conversely, genetic divergence across Greater Bukida islands reflected land‐bridge constrained dispersal. A Makira clade sister to allNesonycterispossibly indicates an origin on Makira. Alternately it reflects Makira's long‐isolated geographic status, as similar results exist for a range of taxa.

    more » « less
  3. null (Ed.)
    Climate refugia, where local populations of species can persist through periods of unfavorable regional climate, play a key role in the maintenance of regional biodiversity during times of environmental change. However, the ability of refugia to buffer biodiversity change may be mediated by the landscape context of refugial habitats. Here, we examined how plant communities restricted to refugial sky islands of alpine tundra in the Colorado Rockies are changing in response to rapid climate change in the region (increased temperature, declining snowpack, and earlier snow melt-out) and if these biodiversity changes are mediated by the area or geographic isolation of the sky island. We resampled plant communities in 153 plots at seven sky islands distributed across the Colorado Rockies at two time points separated by 12 years (2007/2008–2019/2020) and found changes in taxonomic, phylogenetic, and functional diversity over time. Specifically, we found an increase in species richness, a trend toward increased phylogenetic diversity, a shift toward leaf traits associated with the stress-tolerant end of leaf economics spectrum (e.g., lower specific leaf area, higher leaf dry matter content), and a decrease in the functional dispersion of specific leaf area. Importantly, these changes were partially mediated by refugial area but not by geographic isolation, suggesting that dispersal from nearby areas of tundra does not play a strong role in mediating these changes, while site characteristics associated with a larger area (e.g., environmental heterogeneity, larger community size) may be relatively more important. Taken together, these results suggest that considering the landscape context (area and geographic isolation) of refugia may be critical for prioritizing the conservation of specific refugial sites that provide the most conservation value. 
    more » « less
  4. Abstract

    Combating the current biodiversity crisis requires the accurate documentation of population responses to human‐induced ecological change. However, our ability to pinpoint population responses to human activities is often limited to the analysis of populations studied well after the fact. Museum collections preserve a record of population responses to anthropogenic change that can provide critical baseline data on patterns of genetic diversity, connectivity, and population structure prior to the onset of human perturbation. Here, we leverage a spatially replicated time series of specimens to document population genomic responses to the destruction of nearly 90% of coastal habitats occupied by the Savannah sparrow (Passerculus sandwichensis) in California. We sequenced 219 sparrows collected from 1889 to 2017 across the state of California using an exome capture approach. Spatial–temporal analyses of genetic diversity found that the amount of habitat lost was not predictive of genetic diversity loss. Sparrow populations from southern California historically exhibited lower levels of genetic diversity and experienced the most significant temporal declines in genetic diversity. Despite experiencing the greatest levels of habitat loss, we found that genetic diversity in the San Francisco Bay area remained relatively high. This was potentially related to an observed increase in gene flow into the Bay Area from other populations. While gene flow may have minimized genetic diversity declines, we also found that immigration from inland freshwater‐adapted populations into tidal marsh populations led to the erosion of divergence at loci associated with tidal marsh adaptation. Shifting patterns of gene flow through time in response to habitat loss may thus contribute to negative fitness consequences and outbreeding depression. Together, our results underscore the importance of tracing the genomic trajectories of multiple populations over time to address issues of fundamental conservation concern.

    more » « less
  5. Abstract Aim

    To investigate how putative barriers, forest refugia, and ecological gradients across the lower Guineo‐Congolian rain forest shape genetic and phenotypic divergence in the leaf‐folding frogAfrixalus paradorsalis, and examine the role of adjacent land bridge and sky‐islands in diversification.


    The Lower Guineo‐Congolian Forest, the Cameroonian Volcanic Line (CVL), and Bioko Island, Central Africa.


    Afrixalus paradorsalis(Family: Hyperoliidae), an African leaf‐folding frog.


    We used molecular and phenotypic data to investigate diversity and divergence among theA. paradorsalisspecies complex distributed across lowland rain forests, a land bridge island, and mountains in Central Africa. We examined the coincidence of population boundaries, landscape features, divergence times, and spatial patterns of connectivity and diversity, and subsequently performed demographic modelling using genome‐wide SNP variation to distinguish among divergence mechanisms in mainland (riverine barriers, forest refugia, ecological gradients) and land bridge island populations (vicariance, overwater dispersal).


    We detected four genetically distinct allopatric populations corresponding to Bioko Island, the CVL, and two lowland rain forest populations split by the Sanaga River. Although lowland populations are phenotypically indistinguishable, pronounced body size evolution occurs at high elevation, and the timing of the formation of the high elevation population coincides with mountain uplift in the CVL. Spatial analyses and demographic modelling revealed population divergence across mainland Lower Guinea is best explained by forest refugia rather than riverine barriers or ecological gradients, and that the Bioko Island population divergence is best explained by vicariance (marine incursion) rather than overseas dispersal.

    Main conclusions

    We provide growing support for the important role of forest refugia in driving intraspecific divergences in the Guineo‐Congolian rain forest. InA. paradorsalis, sky‐islands in the CVL have resulted in greater genetic and phenotypic divergences than marine incursions of the land bridge Bioko Island, highlighting important differences in patterns of island‐driven diversification in Lower Guinea.

    more » « less