Rivers as barriers to gene flow: A comparative population genomic analysis using causal modeling
Abstract
Earth’s physiographic features shape the genetic evolution of organisms, but understanding how such features act as barriers to gene flow requires quantifying characteristics of both the barrier and the organism. Many barrier characteristics, however, are interdependent and not fully captured by traditional multivariate statistics. Here, we evaluate the use of directed acyclic (causal) graphs and structural equation modeling (SEM) to test the Riverine Barrier Hypothesis using 27 river-spanning population genomic datasets of terrestrial plants and animals associated with 24 rivers across the contiguous United States. These data were paired with seasonality, river width, and river discharge data. SEM analysis revealed patterns not captured by standard approaches. River width had a strong effect on population differentiation, with distinct direct and indirect effects for high and low dispersers. Results suggest a negative width-Fst relationship for low dispersers, which we interpret to be due to topographic context of higher elevation or bedrock-incised rivers. In contrast, high dispersers had a positive relationship, indicating wider rivers present a greater barrier to dispersal. The total effect of river discharge was negligible because its direct effects on population differentiation were canceled out by indirect effects on other river features. Overall, the best-fitting SEM explained 52% of population differentiation for low dispersers and 13% for high dispersers, consistent with the idea that high-dispersing species are less impacted by river barrier effects. This proof of concept shows the utility of causal graphs and SEM at modeling complex relationships between Earth’s physiographic features and the organisms that evolve with them.
Related
- Comparative population genomics reveal stark differences in the genomic structure and drivers of isolation for two sky island bumble bees.
- Machine learning-assisted genomic characterization of population differentiation between Chinese Kirghiz and Kyrgyzstan local sheep
- Genomic diversity, population structure, and relatedness of critically endangered Philippine eagles under conservation management
- A comprehensive comparative genomic analysis of the endogenous Florendovirus-1 in Prunus
- Isolation, identification, and comparative genomic analysis of Vibrio harveyi as a causative agent for a skin ulcer disease of endangered fish Chinese bahaba (Bahaba taipingensis)
BibTeX
@article{maag2026rivers,
title = {{Rivers as barriers to gene flow: A comparative population genomic analysis using causal modeling}},
author = {Garett L Maag and Austin R Biddy and Maya F Stokes and Greer A Dolby},
journal = {Journal of Heredity},
year = {2026},
month = jul,
doi = {10.1093/jhered/esag056},
url = {https://doi.org/10.1093/jhered/esag056},
note = {CC BY-NC 4.0},
}
RIS
TY - JOUR
TI - Rivers as barriers to gene flow: A comparative population genomic analysis using causal modeling
AU - Garett L Maag
AU - Austin R Biddy
AU - Maya F Stokes
AU - Greer A Dolby
PY - 2026
DA - 2026/07/14
JO - Journal of Heredity
DO - 10.1093/jhered/esag056
UR - https://doi.org/10.1093/jhered/esag056
AB - Earth’s physiographic features shape the genetic evolution of organisms, but understanding how such features act as barriers to gene flow requires quantifying characteristics of both the barrier and the organism. Many barrier characteristics, however, are interdependent and not fully captured by traditional multivariate statistics. Here, we evaluate the use of directed acyclic (causal) graphs and structural equation modeling (SEM) to test the Riverine Barrier Hypothesis using 27 river-spanning population genomic datasets of terrestrial plants and animals associated with 24 rivers across the contiguous United States. These data were paired with seasonality, river width, and river discharge data. SEM analysis revealed patterns not captured by standard approaches. River width had a strong effect on population differentiation, with distinct direct and indirect effects for high and low dispersers. Results suggest a negative width-Fst relationship for low dispersers, which we interpret to be due to topographic context of higher elevation or bedrock-incised rivers. In contrast, high dispersers had a positive relationship, indicating wider rivers present a greater barrier to dispersal. The total effect of river discharge was negligible because its direct effects on population differentiation were canceled out by indirect effects on other river features. Overall, the best-fitting SEM explained 52% of population differentiation for low dispersers and 13% for high dispersers, consistent with the idea that high-dispersing species are less impacted by river barrier effects. This proof of concept shows the utility of causal graphs and SEM at modeling complex relationships between Earth’s physiographic features and the organisms that evolve with them.
SN - 0022-1503
SN - 1465-7333
ER -