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

Long-term climatic stability drives accumulation and maintenance of divergent freshwater fish lineages in a temperate biodiversity hotspot

Anthropogenic climate change is forecast to drive regional climate disruption and instability across the globe. These impacts are likely to be exacerbated within biodiversity hotspots, both due to the greater potential for species loss but also to the possibility that endemic lineages might not have experienced significant climatic variation in the past, limiting their evolutionary potential to respond to rapid climate change. We assessed the role of climatic stability on the accumulation and persistence of lineages in an obligate freshwater fish group endemic to the southwest Western Australia (SWWA) biodiversity hotspot. Using 19,426 genomic (ddRAD-seq) markers and species distribution modelling, we explored the phylogeographic history of western (Nannoperca vittata) and little (Nannoperca pygmaea) pygmy perches, assessing population divergence and phylogenetic relationships, delimiting species and estimating changes in species distributions from the Pliocene to 2100. We identified two deep phylogroups comprising three divergent clusters, which showed no historical connectivity since the Pliocene. We conservatively suggest these represent three isolated species with additional intraspecific structure within one widespread species. All lineages showed long-term patterns of isolation and persistence owing to climatic stability but with significant range contractions likely under future climate change. Our results highlighted the role of climatic stability in allowing the persistence of isolated lineages in the SWWA. This biodiversity hotspot is under compounding threat from ongoing climate change and habitat modification, which may further threaten previously undetected cryptic diversity across the region.

Biogeography of Greater Antillean freshwater fishes, with a review of competing hypotheses

In biogeography, vicariance and long-distance dispersal are often characterised as competing scenarios. However, they are related concepts, both relying on collective geological, ecological, and phylogenetic evidence. This is illustrated by freshwater fishes, which may immigrate to islands either when freshwater connections are temporarily present and later severed (vicariance), or by unusual means when ocean gaps are crossed (long-distance dispersal). Marine barriers have a strong filtering effect on freshwater fishes, limiting immigrants to those most capable of oceanic dispersal. The roles of vicariance and dispersal are debated for freshwater fishes of the Greater Antilles. We review three active hypotheses [Cretaceous vicariance, Greater Antilles-Aves Ridge (GAARlandia), long-distance dispersal] and propose long-distance dispersal to be an appropriate model due to limited support for freshwater fish use of landspans. Greater Antillean freshwater fishes have six potential source bioregions (defined from faunal similarity): Northern Gulf of México, Western Gulf of México, Maya Terrane, Chortís Block, Eastern Panamá, and Northern South America. Faunas of the Greater Antilles are composed of taxa immigrating from many of these bioregions, but there is strong compositional disharmony between island and mainland fish faunas (>90% of Antillean species are cyprinodontiforms, compared to <10% in Northern Gulf of México and Northern South America, and ≤50% elsewhere), consistent with a hypothesis of long-distance dispersal. Ancestral-area reconstruction analysis indicates there were 16 or 17 immigration events over the last 51 million years, 14 or 15 of these by cyprinodontiforms. Published divergence estimates and evidence available for each immigration event suggests they occurred at different times and by different pathways, possibly with rafts of vegetation discharged from rivers or washed to sea during storms. If so, ocean currents likely provide critical pathways for immigration when flowing from one landmass to another. On the other hand, currents create dispersal barriers when flowing perpendicularly between landmasses. In addition to high salinity tolerance, cyprinodontiforms collectively display a variety of adaptations that could enhance their ability to live with rafts (small body size, viviparity, low metabolism, amphibiousness, diapause, self-fertilisation). These adaptations likely also helped immigrants establish island populations after arrival and to persist long term thereafter. Cichlids may have used a pseudo bridge (Nicaragua Rise) to reach the Greater Antilles. Gars (Lepisosteidae) may have crossed the Straits of Florida to Cuba, a relatively short crossing that is not a barrier to gene flow for several cyprinodontiform immigrants. Indeed, widespread distributions of Quaternary migrants (Cyprinodon, Gambusia, Kryptolebias), within the Greater Antilles and among neighbouring bioregions, imply that long-distance dispersal is not necessarily inhibitory for well-adapted species, even though it appears to be virtually impossible for all other freshwater fishes.

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.

Freshwater fish diversity in the western Amazon basin shaped by Andean uplift since the Late Cretaceous

South America is home to the highest freshwater fish biodiversity on Earth, and the hotspot of species richness is located in the western Amazon basin. The location of this hotspot is enigmatic, as it is inconsistent with the pattern observed in river systems across the world of increasing species richness towards a river's mouth. Here we investigate the role of river capture events caused by Andean mountain building and repeated episodes of flooding in western Amazonia in shaping the modern-day richness pattern of freshwater fishes in South America, and in Amazonia in particular. To this end, we combine a reconstruction of river networks since 80 Ma with a mechanistic model simulating dispersal, allopatric speciation and extinction over the dynamic landscape of rivers and lakes. We show that Andean mountain building and consequent numerous small river capture events in western Amazonia caused freshwater habitats to be highly dynamic, leading to high diversification rates and exceptional richness. The history of marine incursions and lakes, including the Miocene Pebas mega-wetland system in western Amazonia, played a secondary role.

Phylogeography of a widespread Australian freshwater fish, western carp gudgeon (Eleotridae: Hypseleotris klunzingeri): Cryptic species, hybrid zones, and strong intra-specific divergences

Despite belonging to the most abundant and widespread genus of freshwater fishes in the region, the carp gudgeons of eastern Australia (genus Hypseleotris) have proved taxonomically and ecologically problematic to science since the 19th century. Several molecular studies and a recent taxonomic revision have now shed light on the complex biology and evolutionary history that underlies this group. These studies have demonstrated that carp gudgeons include a sexual/unisexual complex (five sexual species plus an assortment of hemiclonal lineages), many members of which also co-occur with an independent sexual relative, the western carp gudgeon (H. klunzingeri). Here, we fill yet another knowledge gap for this important group by presenting a detailed molecular phylogeographic assessment of the western carp gudgeon across its entire and extensive geographic range. We use a suite of nuclear genetic markers (SNPs and allozymes) plus a matrilineal genealogy (cytb) to demonstrate that H. klunzingeri s.l. also displays considerable taxonomic and phylogeographic complexity. All molecular datasets concur in recognizing the presence of multiple candidate species, two instances of historic between-species admixture, and the existence of a natural hybrid zone between two of the three candidate species found in the Murray-Darling Basin. We also discuss the major phylogeographic patterns evident within each taxon. Together, these analyses provide a robust molecular, taxonomic, and distributional framework to underpin future morphological and ecological investigations on this prominent member of regional freshwater ecosystems in eastern Australia.

Morphology, nuclear SNPs and mate selection reveal that COI barcoding overestimates species diversity in a Mediterranean freshwater amphipod by an order of magnitude

DNA sequence information has revealed many morphologically cryptic species worldwide. For animals, DNA-based assessments of species diversity usually rely on the mitochondrial cytochrome c oxidase subunit I (COI) gene. However, a growing amount of evidence indicate that mitochondrial markers alone can lead to misleading species diversity estimates due to mito-nuclear discordance. Therefore, reports of putative species based solely on mitochondrial DNA should be verified by other methods, especially in cases where COI sequences are identical for different morphospecies or where divergence within the same morphospecies is high. Freshwater amphipods are particularly interesting in this context because numerous putative cryptic species have been reported. Here, we investigated the species status of the numerous mitochondrial molecular operational taxonomic units (MOTUs) found within Echinogammarus sicilianus. We used an integrative approach combining DNA barcoding with mate selection observations, detailed morphometrics and genome-wide double digest restriction site-associated DNA sequencing (ddRAD-seq). Within a relatively small sampling area, we detected twelve COI MOTUs (divergence = 1.8-20.3%), co-occurring in syntopy at two-thirds of the investigated sites. We found that pair formation was random and there was extensive nuclear gene flow among the ten MOTUs co-occurring within the same river stretch. The four most common MOTUs were also indistinguishable with respect to functional morphology. Therefore, the evidence best fits the hypothesis of a single, yet genetically diverse, species within the main river system. The only two MOTUs sampled outside the focal area were genetically distinct at the nuclear level and may represent distinct species. Our study reveals that COI-based species delimitation can significantly overestimate species diversity, highlighting the importance of integrative taxonomy for species validation, especially in hyperdiverse complexes with syntopically occurring mitochondrial MOTUs.

The role of landscape evolution in the genetic diversification of a stream fish Sarcocheilichthys parvus from Southern China

Sarcocheilichthys parvus (Cypriniformes: Gobionidae) is a stream fish which is endemic to sub-tropical coastal drainages in southern China, thus offering a valuable model for understanding how genetic divergence arises in stream-adapting freshwater fishes in this region. Using the mitochondrial Cyt b gene, integrative analyses of phylogeny, population demography, and ancestral area and paleo-drainage reconstructions are carried out to explicitly explore the role of landscape evolution in genetic diversification of S. parvus. The time-calibrated phylogeny of S. parvus indicates the splitting of two major lineages (A and B) at ∼3.66 Ma. Lineage A inhabits the Poyang Lake sub-drainage of the middle Yangtze River, Han River and Pearl River, and can be split into two sub-lineages (A-I and A-II), where sub-lineage A-II can be further sub-divided into three infra-sub-lineages (A-IIa, A-IIb and A-IIc). Except for the infra-sub-lineage A-IIc, which is restricted to the Han River and Pearl River, the other sub-lineages and infra-sub-lineages live exclusively in the Poyang Lake sub-drainage. Lineage B lives in the lower Yangtze River, Qiantang River, Jiaojiang River and Ou River, displaying close genetic relationships among the drainages. Rapid population expansion has occurred since the Late Pleistocene. Our findings indicate that the splitting of lineages A and B could be attributed to geographic isolation due to the Zhe-Min Uplift, acting as a biogeographic barrier before the late Early Pleistocene. Furthermore, the strong genetic divergence within Lineage A could be explained by the isolation role of the Nanling Mountains and Poyang Lake acting as an ecological barrier; while the lack of phylogenetic structure within Lineage B may have been the result of paleo-drainage connections or episodic freshwater connections during the eustatic low stand of sea level in the late Middle-Late Pleistocene.

Extremely divergent COI sequences within an amphipod species complex: A possible role for endosymbionts?

Some heritable endosymbionts can affect host mtDNA evolution in various ways. Amphipods host diverse endosymbionts, but whether their mtDNA has been influenced by these endosymbionts has yet to be considered. Here, we investigated the role of endosymbionts (microsporidians and Rickettsia) in explaining highly divergent COI sequences in Paracalliope fluviatilis species complex, the most common freshwater amphipods in New Zealand. We first contrasted phylogeographic patterns using COI, ITS, and 28S sequences. While molecular species delimitation methods based on 28S sequences supported 3-4 potential species (N, C, SA, and SB) among freshwater lineages, COI sequences supported 17-27 putative species reflecting high inter-population divergence. The deep divergence between NC and S lineages (~20%; 28S) and the substitution saturation on the 3rd codon position of COI detected even within one lineage (SA) indicate a very high level of morphological stasis. Interestingly, individuals infected and uninfected by Rickettsia comprised divergent COI lineages in one of four populations tested, suggesting a potential influence of endosymbionts in mtDNA patterns. We propose several plausible explanations for divergent COI lineages, although they would need further testing with multiple lines of evidence. Lastly, due to common morphological stasis and the presence of endosymbionts, phylogeographic patterns of amphipods based on mtDNA should be interpreted with caution.

Phylogeography reveals the complex impact of the Last Glacial Maximum on New Zealand's terrestrial biota

We review the major phylogeographic patterns in Aotearoa New Zealand's terrestrial flora and fauna that have been associated with the Ōtira Glaciation of the Pleistocene, the end of which coincides with the global Last Glacial Maximum (LGM). We focus on (1) the complexity of biogeographic histories of New Zealand species, with LGM-driven phenomena overlaying the impacts of mountain-building and other drivers of phylogeographic structure; (2) the locations of glacial refugia and sets of taxa which may have shared refugia; and (3) the role of glaciation in driving diversification. We end with a brief focus on the next directions, including what can we learn about New Zealand's glacial history by expanding our phylogeographic toolbox to include genomic methods and hypothesis-driven inference methods. We provide follow-up questions which take advantage of the wealth of phylogeographic data for New Zealand.

Contrasting Patterns of Genetic Diversity and Divergence Between Landlocked and Migratory Populations of Fish Galaxias maculatus, Evaluated Through Mitochondrial DNA Sequencing and Nuclear DNA Microsatellites

Galaxias species are interesting biogeographic models due to their distribution and different types of life cycles, with migratory and landlocked populations. To obtain a better understanding of the genetic consequences of the Quaternary glacial cycles in Galaxias maculatus, in this work we compared landlocked and migratory populations collected in areas that were differentially affected by ice advances and retreats. We included nine populations of G. maculatus, four collected from lakes (landlocked) and five from their associated estuaries/rivers (migratory) in three estuary-lake systems across southern Chile. Genetic analyses were performed using the mitochondrial control region and nine microsatellite loci. Genetic diversity measured with both markers was significantly higher in migratory than in landlocked populations across the study area. The levels of genetic differentiation showed higher differentiation among lakes than estuaries. Genetic diversity was higher in migratory populations located in areas that were less impacted by ice during Quaternary glacial processes. These results may be the consequence of recent recolonization of small freshwater bodies following the Last Glacial Maximum (LGM). Finally, the greatest differentiation was observed in populations that were exposed to continental ice advances and retreats during the LGM. Thus, in the present work we corroborate a pattern of differentiation between lakes and estuaries, using mtDNA sequences and microsatellite nuclear markers. This pattern may be due to a combination of biological factors, i.e., resident non-migratory behaviour or landlocking and natal homing-in, as well as geological factors, i.e., Expansion-Contraction Quaternary glacial biogeographic processes.

Introgression and Species Delimitation in the Longear Sunfish Lepomis megalotis (Teleostei: Percomorpha: Centrarchidae)

Introgression and hybridization are major impediments to genomic-based species delimitation because many implementations of the multispecies coalescent framework assume no gene flow among species. The sunfish genus Lepomis, one of the world's most popular groups of freshwater sport fish, has a complicated taxonomic history. The results of ddRAD phylogenomic analyses do not provide support for the current taxonomy that recognizes two species, L. megalotis and L. peltastes, in the L. megalotis complex. Instead, evidence from phylogenomics and phenotype warrants recognizing six relatively ancient evolutionary lineages in the complex. The introgressed and hybridizing populations in the L. megalotis complex are localized and appear to be the result of secondary contact or rare hybridization events between non-sister species. Segregating admixed populations from our multispecies coalescent analyses identifies six species with moderate to high genealogical divergence, whereas including admixed populations drives all but one lineage below the species threshold of genealogical divergence. Segregation of admixed individuals also helps reveal phenotypic distinctiveness among the six species in morphological traits used by ichthyologists to discover and delimit species over the last two centuries. Our protocols allow for the identification and accommodation of hybridization and introgression in species delimitation. Genomic-based species delimitation validated with multiple lines of evidence provides a path towards the discovery of new biodiversity and resolving long-standing taxonomic problems.

Genomic fingerprints of palaeogeographic history: The tempo and mode of rift tectonics across tropical Africa has shaped the diversification of the killifish genus Nothobranchius (Teleostei: Cyprinodontiformes)

This paper reports a phylogeny of the African killifishes (Genus Nothobranchius, Order Cyprinodontiformes) informed by five genetic markers (three nuclear, two mitochondrial) of 80 taxa (seven undescribed and 73 of the 92 recognized species). These short-lived annual fishes occupy seasonally wet habitats in central and eastern Africa, and their distribution coincides largely with the East African Rift System (EARS). The fossil dates of sister clades used to constrain a chronometric tree of all sampled Nothobranchius recovered the origin of the genus at ~13.27 Mya. It was followed by the radiations of six principal clades through the Neogene. An ancestral area estimation tested competing biogeographical hypotheses to constrain the ancestral origin of the genus to the Nilo-Sudan Ecoregion, which seeded a mid-Miocene dispersal event into the Coastal ecoregion, followed closely (~10 Mya) by dispersals southward across the Mozambique coastal plain into the Limpopo Ecoregion. Extending westwards across the Tanzanian plateau, a pulse of radiations through the Pliocene were associated with dispersals and fragmentation of wetlands across the Kalahari and Uganda Ecoregions. We interpret this congruence of drainage rearrangements with dispersals and cladogenic events of Nothobranchius to reflect congruent responses to recurrent uplift and rifting. The coevolution of these freshwater fishes and wetlands is attributed to ultimate control by tectonics, as the EARS extended southwards during the Neogene. Geobiological consilience of the combined evidence supports a tectonic hypothesis for the evolution of Nothobranchius.