skip to main content

Title: Diversity, divergence and density: How habitat and hybrid zone dynamics maintain a genomic cline in an intertidal barnacle
Abstract Aim

As within‐species genomic data have been shown useful in interpreting broader biogeographic trends, we analysed the mode of population genomic isolation involved in a well‐studied intertidal genomic cline to better understand the mechanisms maintaining it. These results were interpreted in the context of spatial variation in habitat use and availability as well as likely fitness consequences for hybridization between the two lineages.


Pacific coast of North America.


Arthropods (Class Maxillopoda, Order Sessilia, Family Balanidae;Balanus glandula).


Genotype‐by‐sequencing approaches were used to generate single‐nucleotide polymorphism markers across sites sampled between southern Alaska and Southern California. Inference using standard population genomic methods, including analysis of population structure, inbreeding and linkage disequilibrium, was used to identify the steepest transitions across the largest number of loci examined. These data were put in the context of observed population density and habitat availability.


We show that the majority of markers analysed show strong clinal transitions in a very narrow portion of the California coast. Patterns of linkage disequilibrium among markers, along with prior evidence of variation in reproductive potential by latitude and by mitochondrial lineage, suggest some reproductive isolation among the northern and southern lineages ofB. glandulathat are concordant with the drop in population density and habitat availability in central California.

Main Conclusions

A significant clinal transition in genomic diversity is stronger and more localized than previously recognized and exhibits statistical patterns suggesting that the lineages are reproductively and phenotypically distinct in ways that may be ecologically important. As this species has been used to infer process in coastal biogeography, further study of concordant patterns will be important for advancing our understanding of this region.

more » « less
Award ID(s):
1924599 1759797
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Date Published:
Journal Name:
Journal of Biogeography
Medium: X Size: p. 2174-2185
["p. 2174-2185"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Context

    Processes that shape genomic and ecological divergence can reveal important evolutionary dynamics to inform the conservation of threatened species.Fontaineais a genus of rainforest shrubs and small trees including critically endangered and threatened species restricted to narrow, but complex geographic and ecological regions. Several species ofFontaineaare subject to spatially explicit conditions and experience limited intra-specific gene flow, likely generating genetic differentiation and local adaptation.


    Here, we explored the genetic and ecological mechanisms underlying patterns of diversification in two, closely related threatenedFontaineaspecies. Our aim was to compare spatial patterns of genetic variation between the vulnerableFontainea australis(Southern Fontainea) and critically endangeredF. oraria(Coastal Fontainea), endemic to the heterogeneous subtropical region of central, eastern Australia, where large-scale clearing has severely reduced rainforest habitat to a fraction (< 1%) of its pre-European settlement extent.


    We used a set of 10,000 reduced-representation markers to infer genetic relationships and the drivers of spatial genetic variation across the two species. In addition, we employed a combination of univariate and multivariate genome-environment association analysis using a set of topo-climatic variables to explore potential patterns of local adaptation as a factor impacting genomic divergence.


    Our study revealed that Coastal Fontainea have a close genetic relationship with Southern Fontainea. We showed that isolation by distance has played a key role in their genetic variation, indicating that vicariance can explain the spatial genetic distribution of the two species. Genotype-environment analyses showed a strong association with temperature and topographic features, suggesting adaptation to localised thermal environments. We used a multivariate redundancy analysis to identify a range of putatively adapted loci associated with local environmental conditions.


    Divergent selection at the local-habitat scale as a result of dispersal limitations and environmental heterogeneity (including physical barriers) are likely contributors to adaptive divergence between the twoFontaineaspecies. Our findings have presented evidence to indicate that Southern and Coastal Fontainea were comprised of distinct genetic groups and ecotypes, that together may form a single species continuum, with further phenotype research suggested to confirm the current species boundaries. Proactive conservation actions, including assisted migration to enhance the resilience of populations lacking stress-tolerant single nucleotide polymorphisms (SNPs) may be required to secure the long-term future of both taxa. This is especially vital for the critically endangered Coastal Fontainea given projections of habitat decline for the species under future climate scenarios.

    more » « less
  2. Abstract

    Pairs of species that exhibit broadly overlapping distributions, and multiple geographically isolated contact zones, provide opportunities to investigate the mechanisms of reproductive isolation. Such naturally replicated systems have demonstrated that hybridization rates can vary substantially among populations, raising important questions about the genetic basis of reproductive isolation. The topminnows,Fundulus notatusandF. olivaceus, are reciprocally monophyletic, and co‐occur in drainages throughout much of the central and southern United States. Hybridization rates vary substantially among populations in isolated drainage systems. We employed genome‐wide sampling to investigate geographic variation in hybridization, and to assess the possible importance of chromosome fusions to reproductive isolation among nine separate contact zones. The species differ by chromosomal rearrangements resulting from Robertsonian (Rb) fusions, so we hypothesized that Rb fusion chromosomes would serve as reproductive barriers, exhibiting steeper genomic clines than the rest of the genome. We observed variation in hybridization dynamics among drainages that ranged from nearly random mating to complete absence of hybridization. Contrary to predictions, our use of genomic cline analyses on mapped species‐diagnostic SNP markers did not indicate consistent patterns of variable introgression across linkage groups, or an association between Rb fusions and genomic clines that would be indicative of reproductive isolation. We did observe a relationship between hybridization rates and population phylogeography, with the lowest rates of hybridization tending to be found in populations inferred to have had the longest histories of drainage sympatry. Our results, combined with previous studies of contact zones between the species, support population history as an important factor in explaining variation in hybridization rates.

    more » « less
  3. 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
  4. Abstract Aim

    To investigate the structure and rate of gene flow among populations of habitat‐specialized species to understand the ecological and evolutionary processes underpinning their population dynamics and historical demography, including speciation and extinction.


    Peruvian and Argentine Andes.


    Two subspecies of torrent duck (Merganetta armata).


    We sampled 156 individuals in Peru (M. a. leucogenis;Chillón River,n = 57 and Pachachaca River,n = 49) and Argentina (M. a. armata;Arroyo Grande River,n = 33 and Malargüe River,n = 17), and sequenced the mitochondrial DNA (mtDNA) control region to conduct coarse and fine‐scale demographic analyses of population structure. Additionally, to test for differences between subspecies, and across genetic markers with distinct inheritance patterns, a subset of individuals (Peru,n = 10 and Argentina,n = 9) was subjected to partial genome resequencing, obtaining 4,027 autosomal and 189 Z‐linked double‐digest restriction‐associated DNA sequences.


    Haplotype and nucleotide diversities were higher in Peru than Argentina across all markers. Peruvian and Argentine subspecies showed concordant species‐level differences (ΦSTmtDNA = 0.82; ΦSTautosomal = 0.30; ΦSTZ chromosome = 0.45), including no shared mtDNA haplotypes. Demographic parameters estimated for mtDNA using IM and IMa2 analyses, and for autosomal markers using∂a∂i(isolation‐with‐migration model), supported an old divergence (mtDNA = 600,000 years before present (ybp), 95% HPD range = 1.2 Mya to 200,000 ybp; and autosomal∂a∂i = 782,490 ybp), between the two subspecies, characteristic of deeply diverged lineages. The populations were well‐differentiated in Argentina but moderately differentiated in Peru, with low unidirectional gene flow in each country.

    Main conclusions

    We suggest that the South American Arid Diagonal was preexisting and remains a current phylogeographic barrier between the ranges of the two torrent duck subspecies, and the adult territoriality and breeding site fidelity to the rivers define their population structure.

    more » « less
  5. Abstract

    Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between thealbaandpersonatasubspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437SNPloci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.

    more » « less