Using ddRAD-seq phylogeography to test for genetic effects of headwater river capture in suckermouth armored catfish (Loricariidae:Hypostomus) from the central Brazilian Shield

Abstract

AbstractRiver capture is a geological process of potentially great importance in shaping the genetic diversity, distributions, and community composition of freshwater taxa. Using phylogeographic analyses of ddRAD-seq data from suckermouth armored catfish (Hypostomussp. 2) populations, we tested for predicted genetic effects of headwater river capture events in central Brazil, previously supported by geological and community ecological data. We analyzed 227 ddRAD tags (3829 SNP loci) across 42 samples. Molecular results strongly supported sixHypostomusgenetic clusters/lineages, with the deepest divergence ∼1.25 million years ago in the early Pleistocene between a clade from the Upper Paraná and Upper São Francisco river basins versus all other lineages. Consistent with the ‘Paraná Capture Hypothesis’, several lines of evidence supported mid-Pleistocene colonization and vicariant isolation ofHypostomuspopulations from an ancestral Upper Paraná population, including: (1) significant phylogeographic structure, with predicted phylogenetic patterns, (2) higher Paraná lineage diversity, (3) ancestral geographic locations reconstructed in the Paraná basin, and (4) non-random interdrainage dispersal and vicariance events, indicating river captures primarily into the Tocantins and Upper São Francisco basinsc.∼220,000–145,500 years ago. Phylogeographic inference was complicated by lack of lineage monophyly across loci and lineages distributed in multiple basins, the latter of which lent support to the non-mutually exclusive ‘Frequent Interdrainage Dispersal Hypothesis’. However, species tree and demographic modeling results suggested these were artefacts of incomplete sorting of alleles in large ancestral populations over a geologically recent timeframe of divergence. Qualitative and quantitative sensitivity analyses demonstrated that our downstream genetic results were robust to effects of varying ddRAD-seq assembly parameters, which heavily influenced the number of output loci. We predict that codistributed freshwater taxa in Central Brazil may not exhibit phylogeographic patterns similar toHypostomussp. 2 due to complex patterns of superimposed river capture events, or if smaller ancestral population sizes have allowed more complete lineage sorting in other taxa.