Genes meet geology: fish phylogeographic pattern reflects ancient, rather than modern, drainage connections

Gorges partition diversity within New Zealand flathead Galaxias populations

Understanding the landscape factors governing population connectivity in riverine ecosystems represents an ongoing challenge for freshwater biologists. We used DNA sequence analysis to test the hypothesis that major geomorphological features underpin freshwater-limited fish diversity in a tectonically dynamic region of New Zealand. Phylogeographic analysis of 101 Galaxias depressiceps cytochrome b sequences, incorporating 55 localities from southern New Zealand, revealed 26 haplotypes, with only one shared among rivers. We detect strong hierarchical genetic differentiation both among and within river systems. Genetic structuring is particularly pronounced across the Taieri River system (63 individuals from 35 sites, 18 haplotypes), with 92% of variation partitioned among locations. Distinctive within-river genetic clusters are invariably associated with major subcatchment units, typically isolated by substantial gorges. The anomalous distribution of a single lineage across a major drainage divide is consistent with local, tectonically driven headwater capture. We conclude that major landscape features such as gorges can strongly partition riverine fish diversity and constrain freshwater biodiversity.

The roles of aridification and sea level changes in the diversification and persistence of freshwater fish lineages

While the influence of Pleistocene climatic changes on divergence and speciation has been well-documented across the globe, complex spatial interactions between hydrology and eustatics over longer timeframes may also determine species evolutionary trajectories. Within the Australian continent, glacial cycles were not associated with changes in ice cover and instead largely resulted in fluctuations from moist to arid conditions across the landscape. We investigated the role of hydrological and coastal topographic changes brought about by Plio-Pleistocene climatic changes on the biogeographic history of a small Australian freshwater fish, the southern pygmy perch Nannoperca australis. Using 7958 ddRAD-seq (double digest restriction-site associated DNA) loci and 45,104 filtered SNPs, we combined phylogenetic, coalescent and species distribution analyses to assess the various roles of aridification, sea level and tectonics and associated biogeographic changes across southeast Australia. Sea-level changes since the Pliocene and reduction or disappearance of large waterbodies throughout the Pleistocene were determining factors in strong divergence across the clade, including the initial formation and maintenance of a cryptic species, N. 'flindersi'. Isolated climatic refugia and fragmentation due to lack of connected waterways maintained the identity and divergence of inter- and intraspecific lineages. Our historical findings suggest that predicted increases in aridification and sea level due to anthropogenic climate change might result in markedly different demographic impacts, both spatially and across different landscape types.

Paraphyly and evolutionary independent lineages in Gymnotus pantherinus (Gymnotiformes: Gymnotidae) in the Brazilian Atlantic Forest Coastal Streams

In the Brazilian Atlantic Forest (AF), many terrestrial species with broad geographical distributions show high diversity and endemism of intraspecific lineages, as revealed by molecular genetic data. This pattern, however, is less explored in freshwater fishes. Gymnotus pantherinus is an electric fish endemic to the Brazilian coastal drainages that shows a wide distribution, ranging from the states of Bahia to Santa Catarina, an unusual pattern for AF fishes. It has been hypothesized that G. pantherinus is a species complex because distinct morphotypes were described for the species based on morphometric and meristic features. We used mitochondrial and nuclear data to test this hypothesis. Based on phylogenetic inference and multi-locus, multispecies coalescent methods, we identified six independent lineages, flagging them as candidate species. One such lineage is the recently described species G. refugio that is nested within G. pantherinus and renders it paraphyletic, showing it is a species complex. We named G. pantherinus stricto sensu the lineage that includes samples from the type locality (Santos, SP). Our results show that genetic lineages correspond only partially and far exceed the number of previously reported morphotypes. Genetic breaks in the group correspond to landscape features associated with the Serra do Mar mountain range and with riverine dynamics caused by sea level changes during the last glacial maximum. Moreover, we found evidence of river capture events affecting phylogeographic structure in the group. We uncovered an important dimension of diversity in the group and encourage further integration of genetic and phenotypic data. Such integration is a fruitful approach not only to reduce the gap between taxonomy and evolutionary history in Gymnotidae, but also to uncover the real AF biodiversity.

The distribution of Cteniloricaria (Siluriformes: Loricariidae): known and new records in Brazil suggest headwater captures as drivers of disjoint distribution

Abstract Cteniloricaria is a genus of Neotropical armored catfishes belonging to the Loricariinae, currently including two valid species: C. platystoma and C. napova. Cteniloricaria platystoma is presently recorded across the main coastal drainages of the Guiana Shield, from the Sinnamary River, French Guiana, to the Essequibo River basin, Guyana, and is considered to be restricted to the region. Cteniloricaria napova is only known from its type locality at the headwaters of the Paru de Oeste River, Amazon basin, Sipaliwini Savannah, Trio Amerindian territory in Suriname, close to the Brazilian border. Based on a specimen of C. napova, captured in the Cuminapanema River, a tributary to the Curuá River, within Brazilian territory, the geographic distribution of the species and the genus is extended, representing the first record of Cteniloricaria in Brazil. The genus shows a disjoint distribution, and divergence between populations in the north-flowing coastal rivers of the Guianas and the south-flowing Amazon tributaries, and more recent headwater capture between south-flowing Amazon tributaries, may have played a key role in shaping its current distribution. Illustrations, diagnostic characters, morphometrics, description of the habitat where the new specimen was captured, extinction risk assessment, and a discussion of the distribution of the genus are provided.

The role of the south-western Alps as a unidirectional corridor for Mediterranean brown trout (Salmo truttacomplex) lineages

The role of the south-western Alps as a corridor for Mediterranean trout (Salmo trutta complex Linnaeus, 1758) was evaluated in order to understand the influence of the last glacial events in shaping the spatial distribution of the genetic diversity of this salmonid. For this, the allochthonous hypothesis of a man-mediated French origin (19th century) of the Mediterranean trout inhabiting the Po tributaries in the Italian side of the south-western Alps was tested. A total of 412 individuals were analysed at the mitochondrial control region. The phylogenetic classification was carried out by using a Median-Joining Network analysis. Mismatch pair-wise analysis, molecular dating and Kernel density distribution analysis of the main mitochondrial lineages were evaluated to compare past demographic dynamics with the current spatial distribution of genetic diversity. The main outcomes resulted strongly in agreement with a biogeographic scenario where the south-western Alps acted as a unidirectional corridor that permitted the colonization of the upper Durance (Rhône River basin) by trout from the Po River basin. Therefore, the Mediterranean trout should be considered as native also along the Italian side of the south-western Alps and the allochthonous hypothesis should be rejected.

River capture or ancestral polymorphism: an empirical genetic test in a freshwater fish using approximate Bayesian computation

A headwater or river capture is a phenomenon commonly invoked to explain the absence of reciprocal monophyly of genetic lineages among isolated hydrographic basins in freshwater fish. Under the assumption of river capture, a secondary contact between populations previously isolated in different basins explains the observed genetic pattern. However, the absence of reciprocal monophyly could also arise under population isolation through the retention of ancestral of polymorphisms. Here, we applied an approximate Bayesian computation (ABC) framework for estimating the relative probability of scenarios with and without secondary contact. We used Cnesterodon decemmaculatus as a study model because of the multiple possible cases of river capture and the demographic parameters estimated in a previous mitochondrial DNA study that are useful for simulating scenarios to test both hypotheses using the ABC framework. Our results showed that, in general, mitochondrial DNA is useful for distinguishing between these alternative demographic scenarios with reasonable confidence, but in extreme cases (e.g. recent divergence or large population size) there is no power to discriminate between scenarios. Testing hypotheses of drainage rearrangement under a statistically rigorous framework is fundamental for understanding the evolution of freshwater fish fauna as a complement to, or in the absence of, geological evidence.

River Capture and Freshwater Biological Evolution: A Review of Galaxiid Fish Vicariance

Geological processes can strongly affect the distribution and diversification of freshwater-limited species. In particular, a combination of geological and biological data has suggested that Earth history processes can drive vicariant isolation and speciation in non-migratory freshwater fishes. Here, we synthesise recently published geological and freshwater phylogeographic data to illustrate that changes in river drainage geometry are important drivers of galaxiid diversification, both in New Zealand and elsewhere. Major river capture events have led to the isolation and divergence of unique and geographically-restricted lineages, including taxa that are now of prime conservation concern. The parallel phylogeographic effects of drainage shifts have been verified by observations of concordant patterns in co-distributed species. Broadly, this study highlights the dynamic interplay between physical and biological processes in geologically active settings.

Fish diversification at the pace of geomorphological changes: evolutionary history of western Iberian Leuciscinae (Teleostei: Leuciscidae) inferred from multilocus sequence data

The evolutionary history of western Iberian Leuciscinae, obligatory freshwater fish, is directly linked to the evolution of the hydrographic network of the Iberian Peninsula after its isolation from the rest of Europe, which involved dramatic rearrangements such as the transition from endorheic lakes to open basins draining to the Atlantic. Previous phylogenetic research on western Iberian leuciscines, using mainly mitochondrial DNA and more recently one or two nuclear genes, has found contradictory results and there remain many unresolved issues regarding species relationships, taxonomy, and evolutionary history. Moreover, there is a lack of integration between phylogenetic and divergence time estimates and information on the timing of geomorphological changes and paleobasin rearrangements in the Iberian Peninsula. This study presents the first comprehensive fossil-calibrated multilocus coalescent species tree of western Iberian Leuciscinae (including 14 species of Achondrostoma, Iberochondrostoma, Pseudochondrostoma and Squalius endemic to the Iberian Peninsula, seven of which endemic to Portugal) based on seven nuclear genes, and from which we infer their biogeographic history by comparing divergence time estimates to known dated geological events. The phylogenetic pattern suggests slow-paced evolution of leuciscines during the Early-Middle Miocene endorheic phase of the main Iberian river basins, with the shift to exorheism in the late Neogene-Quaternary allowing westward dispersals that resulted in many cladogenetic events and a high rate of endemism in western Iberia. The results of this study also: (i) confirm the paraphyly of S. pyrenaicus with respect to S. carolitertii, and thus the possible presence of a new taxon in the Portuguese Tagus currently assigned to S. pyrenaicus; (ii) support the taxonomic separation of the Guadiana and Sado populations of S. pyrenaicus; (iii) show the need for further population sampling and taxonomic research to clarify the phylogenetic status of A. arcasii from the Minho basin and of the I. lusitanicum populations in the Sado and Tagus basins; and (iv) indicate that A. occidentale, I. olisiponensis and P. duriensis are the most ancient lineages within their respective genera.

Phylogenomic history of enigmatic pygmy perches: implications for biogeography, taxonomy and conservation

Pygmy perches (Percichthyidae) are a group of poorly dispersing freshwater fishes that have a puzzling biogeographic disjunction across southern Australia. Current understanding of pygmy perch phylogenetic relationships suggests past east-west migrations across a vast expanse of now arid habitat in central southern Australia, a region lacking contemporary rivers. Pygmy perches also represent a threatened group with confusing taxonomy and potentially cryptic species diversity. Here, we present the first study of the evolutionary history of pygmy perches based on genome-wide information. Data from 13 991 ddRAD loci and a concatenated sequence of 1 075 734 bp were generated for all currently described and potentially cryptic species. Phylogenetic relationships, biogeographic history and cryptic diversification were inferred using a framework that combines phylogenomics, species delimitation and estimation of divergence times. The genome-wide phylogeny clarified the biogeographic history of pygmy perches, demonstrating multiple east-west events of divergence within the group across the Australian continent. These results also resolved discordance between nuclear and mitochondrial data from a previous study. In addition, we propose three cryptic species within a southwestern species complex. The finding of potentially new species demonstrates that pygmy perches may be even more susceptible to ecological and demographic threats than previously thought. Our results have substantial implications for improving conservation legislation of pygmy perch lineages, especially in southwestern Western Australia.

Barrier Displacement on a Neutral Landscape: Toward a Theory of Continental Biogeography

Macroevolutionary theory posits three processes leading to lineage diversification and the formation of regional biotas: dispersal (species geographic range expansion), speciation (species lineage splitting), and extinction (species lineage termination). The Theory of Island Biogeography (TIB) predicts species richness values using just two of these processes; dispersal and extinction. Yet most species on Earth live on continents or continental shelves, and the dynamics of evolutionary diversification at regional and continental scales are qualitatively different from those that govern the formation of species richness on biogeographic islands. Certain geomorphological processes operating perennially on continental platforms displace barriers to gene flow and organismal dispersal, and affect all three terms of macroevolutionary diversification. For example, uplift of a dissected landscape and river capture both merge and separate portions of adjacent areas, allowing dispersal and larger geographic ranges, vicariant speciation and smaller geographic ranges, and extinction when range sizes are subdivided below a minimum persistence threshold. The TIB also does not predict many biogeographic and phylogenetic patterns widely observed in continentally distributed taxa, including: (i) power function-like species-area relationships; (ii) log-normal distribution of species geographic range sizes, in which most species have restricted ranges (are endemic) and few species have broad ranges (are cosmopolitan); (iii) mid-domain effects with more species toward the geographic center, and more early-branching, species-poor clades toward the geographic periphery; (iv) exponential rates of net diversification with log-linear accumulation of lineages through geological time; and (v) power function-like relationships between species-richness and clade diversity, in which most clades are species-poor and few clades are species-rich. Current theory does not provide a robust mechanistic framework to connect these seemingly disparate patterns. Here we present SEAMLESS (Spatially Explicit Area Model of Landscape Evolution by SimulationS) that generates clade diversification by moving geographic barriers on a continuous, neutral landscape. SEAMLESS is a neutral Landscape Evolution Model (LEM) that treats species and barriers as functionally equivalent with respect to model parameters. SEAMLESS differs from other model-based biogeographic methods (e.g., Lagrange, GeoSSE, BayArea, and BioGeoBEARS) by modeling properties of dispersal barriers rather than areas, and by modeling the evolution of species lineages on a continuous landscape, rather than the evolution of geographic ranges along branches of a phylogeny. SEAMLESS shows how dispersal is required to maintain species richness and avoid clade-wide extinction, demonstrates that ancestral range size does not predict species richness, and provides a unified explanation for the suite of commonly observed biogeographic and phylogenetic patterns listed above. SEAMLESS explains how a simple barrier-displacement mechanism affects lineage diversification under neutral conditions, and is advanced here toward the formulation of a general theory of continental biogeography. [Diversification, extinction, geodispersal, macroevolution, river capture, vicariance.].

Multilocus approach reveals cryptic lineages in the goby Rhinogobius duospilus in Hong Kong streams: Role of paleodrainage systems in shaping marked population differentiation in a city

Drainage history is a well-demonstrated factor that influences the population structure of freshwater inhabitants over a broad geographic scale. However, there has been little research undertaken on such a relationship with freshwater fish on a small geographical scale, especially in Asia. In this study, we investigated the role of local, small drainage systems in affecting the population genetic structure of a freshwater goby, Rhinogobius duospilus, in Hong Kong streams using a multilocus approach. Analyses on nine genetic markers (2 mitochondrial and 7 nuclear markers, including 5 microsatellite markers) reveal prominent and intensive genetic structuring (2.1-5.4% mtDNA sequence divergence) in R. duospilus in Hong Kong. The lineages and clusters recovered from mtDNA data and assignment analysis of nuclear markers coincide with the paleodrainage networks. Furthermore, marked population subdivision between streams located on different side branches (<20km apart) within the same paleodrainage area is observed and gene flow occurs only between closely situated streams that share common paleodrainage confluences. In an extreme case, gene flow is limited between streams that are less than 5km apart. Apparently, such an intensive population structure is attributed to the regional paleodrainage pattern, together with the highly sedentary life style of R. duospilus, which reduces contemporary gene flow and dispersal between populations in neighbouring streams.

The role of river drainages in shaping the genetic structure of capybara populations

The capybara, Hydrochoerus hydrochaeris, is an herbivorous rodent widely distributed throughout most of South American wetlands that lives closely associated with aquatic environments. In this work, we studied the genetic structure of the capybara throughout part of its geographic range in Argentina using a DNA fragment of the mitochondrial control region. Haplotypes obtained were compared with those available for populations from Paraguay and Venezuela. We found 22 haplotypes in 303 individuals. Hierarchical AMOVAs were performed to evaluate the role of river drainages in shaping the genetic structure of capybara populations at the regional and basin scales. In addition, two landscape genetic models, isolation by distance and isolation by resistance, were used to test whether genetic distance was associated with Euclidean distance (i.e. isolation by distance) or river corridor distance (i.e. isolation by resistance) at the basin scale. At the regional scale, the results of the AMOVA grouping populations by mayor river basins showed significant differences between them. At the basin scale, we also found significant differences between sub-basins in Paraguay, together with a significant correlation between genetic and river corridor distance. For Argentina and Venezuela, results were not significant. These results suggest that in Paraguay, the current genetic structure of capybaras is associated with the lack of dispersion corridors through permanent rivers. In contrast, limited structuring in Argentina and Venezuela is likely the result of periodic flooding facilitating dispersion.

Comparative riverscape genetics reveals reservoirs of genetic diversity for conservation and restoration of Great Plains fishes

We used comparative landscape genetics to examine the relative roles of historical events, intrinsic traits and landscape factors in determining the distribution of genetic diversity of river fishes across the North American Great Plains. Spatial patterns of diversity were overlaid on a patch-based graphical model and then compared within and among three species that co-occurred across five Great Plains watersheds. Species differing in reproductive strategy (benthic vs. pelagic-spawning) were hypothesized to have different patterns of genetic diversity, but the overriding factor shaping contemporary patterns of diversity was the signature of past climates and geological history. Allelic diversity was significantly higher at southern latitudes for Cyprinella lutrensis and Hybognathus placitus, consistent with northward expansion from southern Pleistocene refugia. Within the historical context, all species exhibited lowered occupancy and abundance in heavily fragmented and drier upstream reaches, particularly H. placitus; a pelagic-spawning species, suggesting rates of extirpation have outpaced losses of genetic diversity in this species. Within most tributary basins, genetically diverse populations of each species persisted. Hence, reconnecting genetically diverse populations with those characterized by reduced diversity (regardless of their position within the riverine network) would provide populations with greater genetic and demographic resilience. We discuss cases where cross-basin transfer may be appropriate to enhance genetic diversity and mitigate negative effects of climate change. Overall, striking similarities in genetic patterns and in response to fragmentation and dewatering suggest a common strategy for genetic resource management in this unique riverine fish assemblage.

Dynamic reorganization of river basins

River networks evolve as migrating drainage divides reshape river basins and change network topology by capture of river channels. We demonstrate that a characteristic metric of river network geometry gauges the horizontal motion of drainage divides. Assessing this metric throughout a landscape maps the dynamic states of entire river networks, revealing diverse conditions: Drainage divides in the Loess Plateau of China appear stationary; the young topography of Taiwan has migrating divides driving adjustment of major basins; and rivers draining the ancient landscape of the southeastern United States are reorganizing in response to escarpment retreat and coastal advance. The ability to measure the dynamic reorganization of river basins presents opportunities to examine landscape-scale interactions among tectonics, erosion, and ecology.

Influence of drainage divides versus arid corridors on genetic structure and demography of a widespread freshwater turtle, Emydura macquarii krefftii, from Australia

The influence of Pleistocene climatic cycles on Southern Hemisphere biotas is not yet well understood. Australia's eastern coastal margin provides an ideal setting for examining the relative influence of landscape development, sea level fluctuation, and cyclic climatic aridity on the evolution of freshwater biodiversity. We examined the impact of climatic oscillations and physical biogeographic barriers on the evolutionary history of the wide-ranging Krefft's river turtle (Emydura macquarii krefftii), using range-wide sampling (649 individuals representing 18 locations across 11 drainages) and analysis of mitochondrial sequences (∼1.3-kb control region and ND4) and nuclear microsatellites (12 polymorphic loci). A range of phylogeographic (haplotype networks, molecular dating), demographic (neutrality tests, mismatch distributions), and population genetic analyses (pairwise F ST, analysis of molecular variance, Bayesian clustering analysis) were implemented to differentiate between competing demographic (local persistence vs. range expansion) and biogeographic (arid corridor vs. drainage divide) scenarios. Genetic data reveal population genetic structure in Krefft's river turtles primarily reflects isolation across drainage divides. Striking north-south regional divergence (2.2% ND4 p-distance; c. 4.73 Ma, 95% higher posterior density (HPD) 2.08-8.16 Ma) was consistent with long-term isolation across a major drainage divide, not an adjacent arid corridor. Ancient divergence among regional lineages implies persistence of northern Krefft's populations despite the recurrent phases of severe local aridity, but with very low contemporary genetic diversity. Stable demography and high levels of genetic diversity are inferred for southern populations, where aridity was less extreme. Range-wide genetic structure in Krefft's river turtles reflects contemporary and historical drainage architecture, although regional differences in the extent of Plio-Pleistocene climatic aridity may be reflected in current levels of genetic diversity.

Phylogeography of the endangered Otago skink, Oligosoma otagense: population structure, hybridisation and genetic diversity in captive populations

Climatic cooling and substantial tectonic activity since the late Miocene have had a pronounced influence on the evolutionary history of the fauna of New Zealand's South Island. However, many species have recently experienced dramatic range reductions due to habitat fragmentation and the introduction of mammalian predators and competitors. These anthropogenic impacts have been particularly severe in the tussock grasslands of the Otago region. The Otago skink (Oligosoma otagense), endemic to the region, is one of the most critically endangered vertebrates in New Zealand. We use mitochondrial DNA sequence data to investigate the evolutionary history of the Otago skink, examine its population genetic structure, and assess the level of genetic diversity in the individuals in the captive breeding program. Our data indicate that the Otago skink diverged from its closest relatives in the Miocene, consistent with the commencement of tectonic uplift of the Southern Alps. However, there is evidence for past introgression with the scree skink (O. waimatense) in the northern Otago-southern Canterbury region. The remnant populations in eastern Otago and western Otago are estimated to have diverged in the mid-Pliocene, with no haplotypes shared between these two regions. This divergence accounts for 95% of the genetic diversity in the species. Within both regions there is strong genetic structure among populations, although shared haplotypes are generally evident between adjacent localities. Although substantial genetic diversity is present in the captive population, all individuals originate from the eastern region and the majority had haplotypes that were not evident in the intensively managed populations at Macraes Flat. Our data indicate that eastern and western populations should continue to be regarded as separate management units. Knowledge of the genetic diversity of the breeding stock will act to inform the captive management of the Otago skink and contribute to a key recovery action for the species.

Watershed boundaries and geographic isolation: patterns of diversification in cutthroat trout from western North America

Background

For wide-ranging species, intraspecific variation can occur as a result of reproductive isolation from local adaptive differences or from physical barriers to movement. Cutthroat trout (Oncorhynchus clarkii), a widely distributed fish species from North America, has been divided into numerous putative subspecies largely based on its isolation in different watersheds. In this study, we examined mtDNA sequence variation of cutthroat trout to determine the major phylogenetic lineages of this polytypic species. We use these data as a means of testing whether geographic isolation by watershed boundaries can be a primary factor organizing intraspecific diversification.

Results

We collected cutthroat trout from locations spanning almost the entire geographic range of this species and included samples from all major subspecies of cutthroat trout. Based on our analyses, we reveal eight major lineages of cutthroat trout, six of which correspond to subspecific taxonomy commonly used to describe intraspecific variation in this species. The Bonneville cutthroat trout (O. c. utah) and Yellowstone cutthroat trout (O. c. bouvieri) did not form separate monophyletic lineages, but instead formed an intermixed clade. We also document the geographic distribution of a Great Basin lineage of cutthroat trout; a group typically defined as Bonneville cutthroat trout, but it appears more closely related to the Colorado River lineage of cutthroat trout.

Conclusion

Our study indicates that watershed boundaries can be an organizing factor isolating genetic diversity in fishes; however, historical connections between watersheds can also influence the template of isolation. Widely distributed species, like cutthroat trout, offer an opportunity to assess where historic watershed connections may have existed, and help explain the current distribution of biological diversity across a landscape.

Watershed boundaries and geographic isolation: patterns of diversification in cutthroat trout from western North America

BACKGROUND

For wide-ranging species, intraspecific variation can occur as a result of reproductive isolation from local adaptive differences or from physical barriers to movement. Cutthroat trout (Oncorhynchus clarkii), a widely distributed fish species from North America, has been divided into numerous putative subspecies largely based on its isolation in different watersheds. In this study, we examined mtDNA sequence variation of cutthroat trout to determine the major phylogenetic lineages of this polytypic species. We use these data as a means of testing whether geographic isolation by watershed boundaries can be a primary factor organizing intraspecific diversification.

RESULTS

We collected cutthroat trout from locations spanning almost the entire geographic range of this species and included samples from all major subspecies of cutthroat trout. Based on our analyses, we reveal eight major lineages of cutthroat trout, six of which correspond to subspecific taxonomy commonly used to describe intraspecific variation in this species. The Bonneville cutthroat trout (O. c. utah) and Yellowstone cutthroat trout (O. c. bouvieri) did not form separate monophyletic lineages, but instead formed an intermixed clade. We also document the geographic distribution of a Great Basin lineage of cutthroat trout; a group typically defined as Bonneville cutthroat trout, but it appears more closely related to the Colorado River lineage of cutthroat trout.

CONCLUSION

Our study indicates that watershed boundaries can be an organizing factor isolating genetic diversity in fishes; however, historical connections between watersheds can also influence the template of isolation. Widely distributed species, like cutthroat trout, offer an opportunity to assess where historic watershed connections may have existed, and help explain the current distribution of biological diversity across a landscape.

Speciation with gene flow and the genetics of habitat transitions

Whether speciation can advance to completion in the face of initially high levels of gene flow is a very controversial topic in evolutionary biology. Extensive gene exchange is generally considered to homogenize populations and counteract divergence. Moreover, the role of introgressive hybridization in evolution remains largely unexplored in animals, particularly in freshwater zooplankton in which allopatric speciation is considered to be the norm. Our work investigates the genetic structure of two young ecological species: the pond species, Daphnia pulex and the lake species, Daphnia pulicaria. Phylogenetic and population genetics analyses were conducted on mitochondrial NADH dehydrogenase 5 (ND5) gene, the nuclear Lactate dehydrogenase (Ldh) gene and 21 nuclear microsatellite markers in 416 individuals from habitats with various degrees of permanence. The strong and consistent phylogenetic discordance between nuclear and mitochondrial markers suggests a complex evolutionary history of multiple independent habitat transition events that involved hybridization and introgression between lake and pond Daphnia. On the other hand, the low level of contemporary gene flow between adjacent populations indicates the presence of effective habitat isolating barriers. The Daphnia system provides strong evidence for a divergence-with-gene flow speciation model that involves multiple habitat transition events.

Swimming against the current: genetic structure, host mobility and the drift paradox in trematode parasites

Life-cycle characteristics and habitat processes can potentially interact to determine gene flow and genetic structuring of parasitic species. In this comparative study, we analysed the genetic structure of two freshwater trematode species with different life histories using cytochrome c oxidase I gene (COI) sequences and examined the effect of a unidirectional river current on their genetic diversity at 10 sites along the river. We found moderate genetic structure consistent with an isolation-by-distance pattern among subpopulations of Coitocaecum parvum but not in Stegodexamene anguillae. These contrasting parasite population structures were consistent with the relative dispersal abilities of their most mobile hosts (i.e. their definitive hosts). Genetic diversity decreased, as a likely consequence of unidirectional river flow, with increasing distance upstream in C. parvum, which utilizes a definitive host with only restricted mobility. The absence of such a pattern in S. anguillae suggests that unidirectional river flow affects parasite species differently depending on the dispersal abilities of their most mobile host. In conclusion, genetic structure, genetic diversity loss and drift are stronger in parasites whose most mobile hosts have low dispersal abilities and small home ranges. An additional prediction can be made for parasites under unidirectional drift: those parasites that stay longer in their benthic intermediate host or have more than one benthic intermediate hosts would have relatively high local recruitment and hence increased retention of upstream genetic diversity.

A new species of Austrolebias Costa (Cyprinodontiformes: Rivulidae) from northeastern Uruguay, with comments on distribution patterns

A new species of Austrolebias is described based on individuals from the middle and upper río Negro (río Uruguay basin) and río Yaguarón (Patos-Merín system). The new species can be differentiated from all other species of the genus by the unique presence in males of uniform bluish gray pigmentation on flanks (without vertical bands) and unpaired fins. The new species is also distinguished by the combination of characters associated with a reduction of the squamation of the abdominal, preopercular, and opercular regions. The new species presents some morphological characteristics similar to A. gymnoventris and A. luteoflammulatus. The distribution of the new species is concordant with three other species of Austrolebias and may represent a case of drainage rearrangement of the río Negro upstream tributaries (río Uruguay basin) and tributaries of laguna Merín system.

A likelihood framework for estimating phylogeographic history on a continuous landscape

Due to lack of an adequate statistical framework, biologists studying phylogeography are abandoning traditional methods of estimating phylogeographic history in favor of statistical methods designed to test a priori hypotheses. These new methods may, however, have limited descriptive utility. Here, we develop a new statistical framework that can be used to both test a priori hypotheses and estimate phylogeographic history of a gene (and the statistical confidence in that history) in the absence of such hypotheses. The statistical approach concentrates on estimation of geographic locations of the ancestors of a set of sampled organisms. Now we use (2) to derive the likelihood of the ancestral geographic coordinates and the value of the scaled dispersal parameter, given the observed geographic coordinates (assuming known topology and branch lengths). Using a maximum likelihood approach, which is implemented in the new program PhyloMapper, we apply this statistical framework to a 246-taxon mitochondrial genealogy of North American chorus frogs, focusing in detail on one of these species. We demonstrate three lines of evidence for recent northward expansion of the mitochondrion of the coastal clade of Pseudacris feriarum: higher per-generation dispersal distance in the recently colonized region, a noncentral ancestral location, and directional migration. After illustrating one method of accommodating phylogenetic uncertainty, we conclude by discussing how extensions of this framework could function to incorporate a priori ecological and geological information into phylogeographic analyses.

Surviving historical Patagonian landscapes and climate: molecular insights from Galaxias maculatus

BACKGROUND

The dynamic geological and climatic histories of temperate South America have played important roles in shaping the contemporary distributions and genetic diversity of endemic freshwater species. We use mitochondria and nuclear sequence variation to investigate the consequences of mountain barriers and Quaternary glacial cycles for patterns of genetic diversity in the diadromous fish Galaxias maculatus in Patagonia (approximately 300 individuals from 36 locations).

RESULTS

Contemporary populations of G. maculatus, east and west of the Andes in Patagonia, represent a single monophyletic lineage comprising several well supported groups. Mantel tests using control region data revealed a strong positive relationship when geographic distance was modeled according to a scenario of marine dispersal. (r = 0.69, P = 0.055). By contrast, direct distance between regions was poorly correlated with genetic distance (r = -0.05, P = 0.463). Hierarchical AMOVAs using mtDNA revealed that pooling samples according to historical (pre-LGM) oceanic drainage (Pacific vs. Atlantic) explained approximately four times more variance than pooling them into present-day drainage (15.6% vs. 3.7%). Further post-hoc AMOVA tests revealed additional genetic structure between populations east and west of the Chilean Coastal Cordillera (coastal vs. interior). Overall female effective population size appears to have remained relatively constant until roughly 0.5 Ma when population size rapidly increased several orders of magnitude [100x (60x-190x)] to reach contemporary levels. Maximum likelihood analysis of nuclear alleles revealed a poorly supported gene tree which was paraphyletic with respect to mitochondrial-defined haplogroups.

CONCLUSIONS

First diversifying in the central/north-west region of Patagonia, G. maculatus extended its range into Argentina via the southern coastal regions that join the Pacific and Atlantic oceans. More recent gene flow between northern populations involved the most ancient and most derived lineages, and was likely facilitated by drainage reversal(s) during one or more cooling events of the late Pleistocene. Overall female effective population size represents the end result of a widespread and several hundred-fold increase over approximately 0.5 Ma, spanning several climatic fluctuations of the Pleistocene. The minor influence of glacial cycles on the genetic structure and diversity of G. maculatus likely reflects the access to marine refugia during repeated bouts of global cooling. Evidence of genetic structure that was detected on a finer scale between lakes/rivers is most likely the result of both biological attributes (i.e., resident non-migratory behavior and/or landlocking and natal homing in diadromous populations), and the Coastal Cordillera as a dispersal barrier.

Surviving historical Patagonian landscapes and climate: molecular insights from Galaxias maculatus

Background

The dynamic geological and climatic histories of temperate South America have played important roles in shaping the contemporary distributions and genetic diversity of endemic freshwater species. We use mitochondria and nuclear sequence variation to investigate the consequences of mountain barriers and Quaternary glacial cycles for patterns of genetic diversity in the diadromous fish Galaxias maculatus in Patagonia (~300 individuals from 36 locations).

Results

Contemporary populations of G. maculatus, east and west of the Andes in Patagonia, represent a single monophyletic lineage comprising several well supported groups. Mantel tests using control region data revealed a strong positive relationship when geographic distance was modeled according to a scenario of marine dispersal. (r = 0.69, P = 0.055). By contrast, direct distance between regions was poorly correlated with genetic distance (r = -0.05, P = 0.463). Hierarchical AMOVAs using mtDNA revealed that pooling samples according to historical (pre-LGM) oceanic drainage (Pacific vs. Atlantic) explained approximately four times more variance than pooling them into present-day drainage (15.6% vs. 3.7%). Further post-hoc AMOVA tests revealed additional genetic structure between populations east and west of the Chilean Coastal Cordillera (coastal vs. interior). Overall female effective population size appears to have remained relatively constant until roughly 0.5 Ma when population size rapidly increased several orders of magnitude [100× (60×-190×)] to reach contemporary levels. Maximum likelihood analysis of nuclear alleles revealed a poorly supported gene tree which was paraphyletic with respect to mitochondrial-defined haplogroups.

Conclusions

First diversifying in the central/north-west region of Patagonia, G. maculatus extended its range into Argentina via the southern coastal regions that join the Pacific and Atlantic oceans. More recent gene flow between northern populations involved the most ancient and most derived lineages, and was likely facilitated by drainage reversal(s) during one or more cooling events of the late Pleistocene. Overall female effective population size represents the end result of a widespread and several hundred-fold increase over approximately 0.5 Ma, spanning several climatic fluctuations of the Pleistocene. The minor influence of glacial cycles on the genetic structure and diversity of G. maculatus likely reflects the access to marine refugia during repeated bouts of global cooling. Evidence of genetic structure that was detected on a finer scale between lakes/rivers is most likely the result of both biological attributes (i.e., resident non-migratory behavior and/or landlocking and natal homing in diadromous populations), and the Coastal Cordillera as a dispersal barrier.

Population subdivision in Siamese mud carp Henicorhynchus siamensis in the Mekong River basin: implications for management

A molecular approach was employed to investigate stock structure in Siamese mud carp Henicorhynchus siamensis populations collected from 14 sites across mainland south-east Asia, with the major focus being the lower Mekong River basin. Spatial analysis of a mitochondrial DNA fragment (ATPase 6 and 8) identified four stocks in the Mekong River basin that were all significantly differentiated from a population in the nearby Khlong River, Thailand. In the Mekong River basin, populations in northern Lao People's Democratic Republic and northern Thailand represent two independent stocks, and samples from Thai tributaries group with those from adjacent Mekong sites above the Khone Falls to form a third stock. All sites below the Khone Falls constituted a single vast stock that includes Cambodia and the Mekong Delta in Vietnam. While H. siamensis is considered currently to undertake extensive annual migrations across the Mekong River basin, the data presented here suggest that natural gene flow may occur over much more restricted geographical scales within the basin, and hence populations may need to be managed at finer spatial scales than at the whole-of-drainage-basin level.

Evolutionary dynamics of transferrin in Notropis

This study reports on the development of comparative data for the transferrin (TF) gene in cyprinid fishes, focusing on the genus Notropis. While previous studies have suggested varied patterns of adaptation influencing the diversity at this gene locus both within and among species, sequence data for (TF) in Notropis exhibit limited evidence for selection. However, there are significant effects of sequence variation associated with Notropis lutipinnis, suggesting some form of diversifying selection acting among populations of this species. Overall, the gene performs well as a second locus for phylogenetic and biogeographic inference and may help improve description of the pattern and the process of diversification in Notropis.