skip to main content

Title: A legacy of geo-climatic complexity and genetic divergence along the lower Colorado River: Insights from the geological record and 33 desert-adapted animals
Aim: To review the histories of the Colorado River and North American monsoon system to ascertain their effects on the genetic divergence of desert‐adapted animals. Location: Lower Colorado River region, including Mojave and Sonoran deserts, United States. Methods: We synthesized recent geological literature to summarize initiation phases of lower Colorado River evolution, their discrepancies, and potential for post‐vicariance dispersal of animals across the river. We simulated data under geological models and performed a meta‐analysis of published and unpublished genetic data including population diversity metrics, relatedness and historical migration rates to assess alternative divergence hypotheses. Results: The two models for arrival of the Colorado River into the Gulf of California impose east‐west divergence ages of 5.3 and 4.8 Ma, respectively. We found quantifiable river‐associated differentiation in the lower Colorado River region in reptiles, arachnids and mammals relative to flying insects. However, topological statistics, historical migration rates and cross‐river extralimital populations suggest that the river should be considered a leaky barrier that filters, rather than prevents, gene flow. Most markers violated neutrality tests. Differential adaptation to monsoon‐based precipitation differences may contribute to divergence between Mojave and Sonoran populations and should be tested. Main Conclusions: Rivers are dynamic features that can both limit and facilitate gene flow through time, the impacts of which are mitigated by species‐specific life history and dispersal traits. The Southwest more » is a geo‐climatically complex region with the potential to produce pseudocongruent patterns of genetic divergence, offering a good setting to evaluate intermediate levels of geological‐biological (geobiological) complexity. « less
Authors:
; ;
Award ID(s):
1754587
Publication Date:
NSF-PAR ID:
10157737
Journal Name:
Journal of biogeography
Volume:
46
Issue:
11
Page Range or eLocation-ID:
2479–2505
ISSN:
1365-2699
Sponsoring Org:
National Science Foundation
More Like this
  1. Identifying the evolutionary and ecological mechanisms that drive lineage diversification in the species-rich tropics is of broad interest to evolutionary biologists. Here, we use phylogeographic and demographic analyses of genomic scale RADseq data to assess the impact of a large geographic feature, the Amazon River, on lineage formation in a venomous pitviper, Bothrops atrox. We compared genetic differentiation in samples from four sites near Santarem, Brazil that spanned the Amazon and represented major habitat types. A species delimitation analysis identified each population as a distinct evolutionary lineage while a species tree analysis with populations as taxa revealed a phylogenetic tree consistent with dispersal across the Amazon from north to south. Phylogenetic analyses of mtDNA variation confirmed this pattern and suggest that all lineages originated during the mid- to late-Pleistocene. Historical demographic analyses support a population model of lineage formation through isolation between lineages with low ongoing migration between large populations and reject a model of differentiation through isolation by distance alone. Our results provide a rare example of a phylogeographic pattern demonstrating dispersal over evolutionary time scales across a large tropical river and suggest a role for the Amazon River as a driver of in-situ divergence by both impeding (butmore »not preventing) gene flow and through parapatric differentiation along an ecological gradient.« less
  2. Abstract

    Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these, we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulfmore »of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro’s current range reflects postglacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects. [Phylogenomics; phylogeography; rooting; Sonoran Desert.]

    « less
  3. Understanding the timescales on which different geologic processes influence genetic divergence is crucial to defining and testing geogenomic hypotheses and characterizing Earth- life evolution. To see if we can recover a genetic signal produced by a hypothetical physical barrier to gene flow, we used a geographically explicit simulation approach. We used the CDMetaPop software to simulate heritable genetic, nonadaptive, data for 20 geographically distinct populations distributed throughout the Baja California peninsula of Mexico, a landscape where a transpeninsular seaway barrier has been proposed to have isolated the southern peninsula and caused the observed latitudinal genetic divergence in over 80 terrestrial species. We simulated 10,000 generations of isolation by a barrier under two dispersal scenarios (1 km and 100 km of max. dispersal from population of origin per generation) and three DNA substitution rates (10-7, 10-8 and 10-9 nucleotide substitutions per site per generation). Our simulations indicate that a physical barrier can produce strong genetic divergence within 10,000 generations, comparable to the continuum of values observed in nature for different taxonomic groups and geological settings. We found that the generation time of the organism was by far the most important factor dictating the rate of divergence. Evaluating different generation times (0.02,more »0.2, 2 and 20 years), showed that species with longer generation times require longer periods of isolation to accumulate genetic divergence over 10k generations (~1 My). Simulating 10,000 generations of gene flow following removal of the barrier showed that the divergence signal eroded quickly, in less than 1,000 generations in every scenario, a pattern supported by theory from population genetics. These results are particularly relevant to geogenomic studies because they show that ephemeral gene flow barriers produce different magnitudes of genetic signals depending on attributes of the organism, particularly generation time, and that if reproductive isolation is not achieved during isolation, then the evolutionary signal of an ephemeral barrier may not develop. This work helps guide the limits of detectability when integrating genomic data with geological and climatic processes.« less
  4. Many freshwater gastropod species face extinction, including 79% of species in the family Pleuroceridae. The Oblong Rocksnail, Leptoxis compacta , is a narrow range endemic pleurocerid from the Cahaba River basin in central Alabama that has seen rapid range contraction in the last 100 years. Such a decline is expected to negatively affect genetic diversity in the species. However, precise patterns of genetic variation and gene flow across the restricted range of L. compacta are unknown. This lack of information limits our understanding of human impacts on the Cahaba River system and Pleuroceridae. Here, we show that L. compacta has likely seen a species-wide decline in genetic diversity, but remaining populations have relatively high genetic diversity. We also report a contemporary range extension compared to the last published survey. Our findings indicate that historical range contraction has resulted in the absence of common genetic patterns seen in many riverine taxa like isolation by distance as the small distribution of L. compacta allows for relatively unrestricted gene flow across its remaining range despite limited dispersal abilities. Two collection sites had higher genetic diversity than others, and broodstock sites for future captive propagation and reintroduction efforts should utilize sites identified here asmore »having the highest genetic diversity. Broadly, our results support the hypothesis that range contraction will result in the reduction of species-wide genetic diversity, and common riverscape genetic patterns cannot be assumed to be present in species facing extinction risk.« less
  5. Reproductive isolation is necessary for population divergence to lead to the formation of separate species. This can occur due to physical isolation of populations, which drives allopatric speciation, or other methods of isolation, such as sympatric speciation where the diverging species are physically in the same range, but structural genomic changes or mutations cause the population to diverge into two different species. Parapatric speciation occurs when populations that are geographically adjacent to each other diverge, which can be driven by adaptations to environmental differences, even with ongoing gene flow. Two desert- adapted brittlebush species, Encelia farinosa and Encelia californica, diverged less than 1 million years ago (Singhal et al., 2020) and have a parapatric distribution, residing in different environments in the Mojave and Sonoran deserts. Encelia farinosa (Brittlebush) has unique silvery leaves that are covered in tiny hairs (leaf pubescence) to better control leaf temperature in the hot and arid conditions of the Sonoran Desert. Encelia californica (California Brittlebush) does not display leaf pubescence and is found in a smaller region of the Mediterranean-like environment of the west coast of North America. Encelia californica is exposed to more precipitation than most other Encelia species. Even with their different morphologies, thesemore »two species are still able to hybridize and create fertile offspring (Clark, 1998). Using PacBio sequencing and Hi-C scaffolding, we assembled and annotated reference genomes for both species to investigate the genomic basis of reproductive isolation in these two species. The scaffold N50/L50 are 10 scaffolds and 76.3 Mbp, and 12 scaffolds and 64.5 Mbp for E. farinosa and E. californica, respectively. Using comparative genomic analyses such as tests for differential adaptation and chromosomal translocations will help reveal whether the drivers of speciation in the Encelia radiation were external (e.g., geologic/climatic) or internal (e.g., genomic rearrangement). These analyses will also help answer how accumulated genomic differences can cause speciation in populations that are not geographically isolated. Analyses such as these are new, exciting sources of information for testing geogenomic and other Earth- life hypotheses.« less