Genomewide SNP data reveal cryptic phylogeographic structure and microallopatric divergence in a rapids-adapted clade of cichlids from the Congo River

Impacts of Quaternary Climatic Changes on the Diversification of Riverine Cichlids in the Lower Congo River

Climatic and geomorphological changes during the Quaternary period impacted global patterns of speciation and diversification across a wide range of taxa, but few studies have examined these effects on African riverine fish. The lower Congo River is an excellent natural laboratory for understanding complex speciation and population diversification processes, as it is hydrologically extremely dynamic and recognized as a continental hotspot of diversity harboring many narrowly endemic species. A previous study using genome-wide SNP data highlighted the importance of dynamic hydrological regimes to the diversification and speciation in lower Congo River cichlids. However, historical climate and hydrological changes (e.g., reduced river discharge during extended dry periods) have likely also influenced ichthyofaunal diversification processes in this system. The lower Congo River offers a unique opportunity to study climate-driven changes in river discharge, given the massive volume of water from the entire Congo basin flowing through this short stretch of the river. Here, we, for the first time, investigate the impacts of paleoclimatic factors on ichthyofaunal diversification in this system by inferring divergence times and modeling patterns of gene flow in four endemic lamprologine cichlids, including the blind cichlid, Lamprologus lethops. Our results suggest that Quaternary climate changes associated with river discharge fluctuations may have impacted the diversification of species along the system and the emergence of cryptophthalmic phenotype in some endemic species. Our study, using reduced representation sequencing (2RADseq), indicates that the lower Congo River lamprologines emerged during the Early-Middle Pleistocene transition, characterized as one of the earth's major climatic transformation periods. Modeling results suggest that gene flow across populations and between species was not constant but occurred in temporally constrained pulses. We show that these results correlate with glacial-interglacial fluctuations. The current hyper-diverse fish assemblages of the lower Congo River riverscape likely reflect the synergistic effects of multiple drivers fueling complex evolutionary processes through time.

Ecological genetics of isolated loach populations indicate compromised adaptive potential

Many endangered species live in fragmented and isolated populations with low genetic variability, signs of inbreeding, and small effective population sizes - all features elevating their extinction risk. The flat-headed loach (Oreonectes platycephalus), a small noemacheilid fish, is widely across southern China, but only in the headwaters of hillstreams; as a result, they are spatially isolated from conspecific populations. We surveyed single nucleotide polymorphisms in 16 Hong Kong populations of O. platycephalus to determine whether loach populations from different streams were genetically isolated from each other, showed low levels of genetic diversity, signs of inbreeding, and had small contemporary effective population sizes. Estimates of average observed heterozygosity (HO = 0.0473), average weighted nucleotide diversity (πw = 0.0546) and contemporary effective population sizes (Ne = 10.2 ~ 129.8) were very low, and several populations showed clear signs of inbreeding as judged from relatedness estimates. The degree of genetic differentiation among populations was very high (average FST = 0.668), even over short geographic distances (<1.5 km), with clear patterns of isolation by distance. These results suggest that Hong Kong populations of O. platycephalus have experienced strong genetic drift and loss of genetic variability because sea-level rise after the last glaciation reduced connectedness among paleodrainages, isolating populations in headwaters. All this, together with the fact that the levels of genetic diversity and contemporary effective population sizes within O. platycephalus populations are lower than most other freshwater fishes, suggests that they face high local extinction risk and have limited capacity for future adaptation.

Genetic structure and origin of emu populations in Japanese farms inferred from large-scale SNP genotyping based on double-digest RAD-seq

The emu is a novel poultry species in Japan. However, Japanese farmed emu populations have reduced genetic diversity owing to inbreeding. We have previously suggested that there are genetic resources in the Tohoku Safari Park (TSP) and Fuji/Kakegawa Kachoen Garden Park (FGP/KGP) to extend the genetic diversity of commercial emu farms based on microsatellite (SSR) and mitochondrial DNA. However, those markers provide relatively poor information. Thus, we investigated the genetic structure of farmed Japanese populations based on a large-scale genotyping system using RAD-seq and verified the usefulness of TSP and FGP/KGP as genetic resources for expanding genetic diversity. Admixture, phylogenetic, and principal component analyses based on 28,676 SNPs showed that TSP individuals were ancestors in the Okhotsk Emu Farm (OEF). FGP/KGP individuals showed a unique genetic component that differed from that of the others. We have previously reported that the mitochondrial haplotypes of FGP/KGP were shared with an isolated wild population in eastern Australia. These results suggest that FGP/KGP individuals originated from an eastern Australia isolated population different from other populations including ancestral of OEF/TSP. Our results would provide information for the development of Japanese emu farms and industry and for the conservation of genetic resources in the Australian wild emu.

Phenotypic rate and state are decoupled in response to river-to-lake transitions in cichlid fishes

Geographic access to isolated ecosystems is an important catalyst of adaptive radiation. Cichlid fishes repeatedly colonized rift, crater, and volcanic lakes from surrounding rivers. We test the "lake effect" on the phenotypic rate and state across 253 cichlid species. The rate of evolution was consistently higher (~10-fold) in lakes, and consistent across different dimensions of the phenotype. Rate shifts tended to occur coincident with or immediately following river-to-lake transitions, generally resulting in 2- to 5-fold faster rates than in the founding riverine lineage. By contrast, river- and lake-dwelling cichlids exhibit considerable overlap in phenotypes, generally with less disparity in lakes, but often different evolutionary optima. Taken together, these results suggest that lake radiations rapidly expand into niches largely already represented by ancestral riverine lineages, albeit in different frequencies. Lakes may provide ecological opportunity via ecological release (e.g., from predators/competitors) but need not be coupled with access to novel ecological niches.

Mitogenomic Characterization of Cameroonian Endemic Coptodon camerunensis (Cichliformes: Cichlidae) and Matrilineal Phylogeny of Old-World Cichlids

The mitogenomic evolution of old-world cichlids is still largely incomplete in Western Africa. In this present study, the complete mitogenome of the Cameroon endemic cichlid, Coptodon camerunensis, was determined by next-generation sequencing. The mitogenome was 16,557 bp long and encoded with 37 genes (13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a control region). The C. camerunensis mitogenome is AT-biased (52.63%), as exhibited in its congener, Coptodon zillii (52.76% and 53.04%). The majority of PCGs start with an ATG initiation codon, except COI, which starts with a GTG codon and five PCGs and ends with the TAA termination codon and except seven PCGs with an incomplete termination codon. In C. camerunensis mitogenome, most tRNAs showed classical cloverleaf secondary structures, except tRNA-serine with a lack of DHU stem. Comparative analyses of the conserved blocks of two Coptodonini species control regions revealed that the CSB-II block was longer than other blocks and contained highly variable sites. Using 13 concatenated PCGs, the mitogenome-based Bayesian phylogeny easily distinguished all the examined old-world cichlids. Except for Oreochromini and Coptodinini tribe members, the majority of the taxa exhibited monophyletic clustering within their respective lineages. C. camerunensis clustered closely with Heterotilapia buttikoferi (tribe Heterotilapiini) and had paraphyletic clustering with its congener, C. zillii. The Oreochromini species also displayed paraphyletic grouping, and the genus Oreochromis showed a close relationship with Coptodinini and Heterotilapiini species. In addition, illustrating the known distribution patterns of old-world cichlids, the present study is congruent with the previous hypothesis and proclaims that prehistoric geological evolution plays a key role in the hydroclimate of the African continent during Mesozoic, which simultaneously disperses and/or colonizes cichlids in different ichthyological provinces and Rift Lake systems in Africa. The present study suggests that further mitogenomes of cichlid species are required, especially from western Africa, to understand their unique evolution and adaptation.

Exceptional levels of species discovery ameliorate inferences of the biogeography and diversification of an Afrotropical catfish family

Endeavours in species discovery, particularly the characterisation of cryptic species, have been greatly aided by the application of DNA molecular sequence data to phylogenetic reconstruction and inference of evolutionary and biogeographic processes. However, the extent of cryptic and undescribed diversity remains unclear in tropical freshwaters, where biodiversity is declining at alarming rates. To investigate how data on previously undiscovered biodiversity impacts inferences of biogeography and diversification dynamics, we generated a densely sampled species-level family tree of Afrotropical Mochokidae catfishes (220 valid species) that was ca. 70 % complete. This was achieved through extensive continental sampling specifically targeting the genus Chiloglanis a specialist of the relatively unexplored fast-flowing lotic habitat. Applying multiple species-delimitation methods, we report exceptional levels of species discovery for a vertebrate genus, conservatively delimiting a staggering ca. 50 putative new Chiloglanis species, resulting in a near 80 % increase in species richness for the genus. Biogeographic reconstructions of the family identified the Congo Basin as a critical region in the generation of mochokid diversity, and further revealed complex scenarios for the build-up of continental assemblages of the two most species rich mochokid genera, Synodontis and Chiloglanis. While Syndontis showed most divergence events within freshwater ecoregions consistent with largely in situ diversification, Chiloglanis showed much less aggregation of freshwater ecoregions, suggesting dispersal as a key diversification process in this older group. Despite the significant increase in mochokid diversity identified here, diversification rates were best supported by a constant rate model consistent with patterns in many other tropical continental radiations. While our findings highlight fast-flowing lotic freshwaters as potential hotspots for undescribed and cryptic species diversity, a third of all freshwater fishes are currently threatened with extinction, signifying an urgent need to increase exploration of tropical freshwaters to better characterise and conserve its biodiversity.

Strong species structure but weak geographical structure in demersal Lake Victoria cichlids

Studying phenotypic and genetic differentiation between very young species can be very informative with regard to learning about processes of speciation. Identifying and characterizing genetic species structure and distinguishing it from spatial genetic structure within a species is a prerequisite for this and is often not given sufficient attention. Young radiations of cichlid fish are classical speciation study systems. However, it is only during the past decade that population genomics based on next-generation sequencing has begun to provide the power to resolve species and distinguish speciation from spatial population structure for the youngest of these radiations. The Lake Victoria haplochromine cichlids constitute the youngest large cichlid fish radiation, probably <20,000 years old. Earlier work showed that communities of rocky reef cichlids are composed of many reciprocally monophyletic species despite their very recent origins. Here, we build on this work by studying assemblages of offshore demersal cichlids, adding analyses of within-species spatial structure to the sympatric species structure. We sampled seven multispecies communities along a 6-km-long transect from one side of the Mwanza Gulf to the other side. We investigated whether phenotypically diagnosed putative species are reciprocally monophyletic and whether such monophyly is stable across species geographic ranges. We show that all species are genetically strongly differentiated in sympatry, that they are reciprocally monophyletic, and that monophyly is stable across distribution ranges. We found significant differentiation between geographically distinct populations in two species, but no or weak isolation by distance. We further found subtle but significant morphological differences between all species and a linear relationship between genomic and morphological distance which suggests that differences in morphology begin to accumulate after speciation has already affected genome-wide restrictions of gene flow.

Evolutionary convergence in body shape obscures taxonomic diversity in species of the African Labeo forskalii group: Case study of L. parvus Boulenger 1902 and L. ogunensis Boulenger 1910

Labeo is the third most diverse genus of African cyprinids and is widely distributed across the continent. Labeo parvus, a small species originally described from the Congo basin, has been considered the only species of the L. forskalii group distributed across five African ichthyofaunal provinces (Nilo-Sudan, Congo, Cuanza, and Upper and Lower Guinea). However, morphological similarity between L. parvus and numerous congeners remains a central cause of taxonomic confusion within the genus. Here we employed a phylogenetic comparative approach to assess phenotypic convergence among species of the L. forskalii group, investigate the taxonomic status of L. parvus sensu lato (sl) in west Africa, and reevaluate the composition and distribution of L. parvus sensu stricto (ss). Our phylogenetic analysis provides no support for a sister relationship between L. parvus ss and any of the west African Labeo parvus-like species. Geometric morphometric and molecular phylogenetic data indicate that L. parvus ss is a Congo basin endemic, and seemingly ecologically equivalent species found in west Africa are L. ogunensis, L. obscurus and other undescribed or previously synonymized species. We discuss our findings in terms of convergent evolution using phylomorphospace and tests for phylogenetic signal.

Exon-based phylogenomics and the relationships of African cichlid fishes: tackling the challenges of reconstructing phylogenies with repeated rapid radiations

African cichlids (subfamily: Pseudocrenilabrinae) are among the most diverse vertebrates, and their propensity for repeated rapid radiation has made them a celebrated model system in evolutionary research. Nonetheless, despite numerous studies, phylogenetic uncertainty persists, and riverine lineages remain comparatively underrepresented in higher-level phylogenetic studies. Heterogeneous gene histories resulting from incomplete lineage sorting (ILS) and hybridization are likely sources of uncertainty, especially during episodes of rapid speciation. We investigate relationships of Pseudocrenilabrinae and its close relatives while accounting for multiple sources of genetic discordance using species tree and hybrid network analyses with hundreds of single-copy exons. We improve sequence recovery for distant relatives, thereby extending the taxonomic reach of our probes, with a hybrid reference guided/de novo assembly approach. Our analyses provide robust hypotheses for most higher-level relationships and reveal widespread gene heterogeneity, including in riverine taxa. ILS and past hybridization are identified as sources of genetic discordance in different lineages. Sampling of various Blenniiformes (formerly Ovalentaria) adds strong phylogenomic support for convict blennies (Pholidichthyidae) as sister to Cichlidae, and points to other potentially useful protein-coding markers across the order. A reliable phylogeny with representatives from diverse environments will support ongoing taxonomic and comparative evolutionary research in the cichlid model system.

ddRAD-Seq reveals evolutionary insights into population differentiation and the cryptic phylogeography of Hyporhamphus intermedius in Mainland China

Species differentiation and local adaptation in heterogeneous environments have attracted much attention, although little is known about the mechanisms involved. Hyporhamphus intermedius is an anadromous, brackish-water halfbeak that is widely distributed in coastal areas and hyperdiverse freshwater systems in China, making it an interesting model for research on phylogeography and local adaptation. Here, 156 individuals were sampled at eight sites from heterogeneous aquatic habitats to examine environmental and genetic contributions to phenotypic divergence. Using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in the specimens from the different watersheds, 5498 single nucleotide polymorphisms (SNPs) were found among populations, with obvious population differentiation. We find that present-day Mainland China populations are structured into distinct genetic clusters stretching from southern and northern ancestries, mirroring geography. Following a transplant event in Plateau Lakes, there were virtually no variations of genetic diversity occurred in two populations, despite the fact two main splits were unveiled in the demographic history. Additionally, dorsal, and anal fin traits varied widely between the southern group and the others, which highlighted previously unrecognized lineages. We then explore genotype-phenotype-environment associations and predict candidate loci. Subgroup ranges appeared to correspond to geographic regions with heterogeneous hydrological factors, indicating that these features are likely important drivers of diversification. Accordingly, we conclude that genetic and phenotypic polymorphism and a moderate amount of genetic differentiation occurred, which might be ascribed to population subdivision, and the impact of abiotic factors.

Riverscape genomics of cichlid fishes in the lower Congo: Uncovering mechanisms of diversification in an extreme hydrological regime

Freshwater fishes are notably diverse, given that freshwater habitat represents a tiny fraction of the earth's surface, but the mechanisms generating this diversity remain poorly understood. Rivers provide excellent models to understand how freshwater diversity is generated and maintained across heterogeneous habitats. In particular, the lower Congo River (LCR) consists of a dynamic hydroscape exhibiting extraordinary aquatic biodiversity, endemicity, morphological and ecological specialization. Previous studies have suggested that the numerous high-energy rapids throughout the LCR form physical barriers to gene flow, thus facilitating diversification and speciation, generating ichthyofaunal diversity. However, this hypothesis has not been fully explored using genome-wide SNPs for fish species distributed across the LCR. Here, we examined four lamprologine cichlids endemic to the LCR that are distributed along the river without range overlap. Using genome-wide SNP data, we tested the hypotheses that high-energy rapids serve as physical barriers to gene flow that generate genetic divergence at inter- and intraspecific levels, and that gene flow occurs primarily in a downstream direction. Our results are consistent with the prediction that powerful rapids sometimes act as a barrier to gene flow but also suggest that, at certain temporal and spatial scales, they may provide multidirectional dispersal opportunities for riverine rheophilic cichlid fishes. These results highlight the complexity of diversification processes in rivers and the importance of assessing such processes across different riverscapes.

Rapid Parallel Morphological and Mechanical Diversification of South American Pike Cichlids (Crenicichla)

Explosive bouts of diversification are one of the most conspicuous features of the tree of life. When such bursts are repeated in similar environments it suggests some degree of predictability in the evolutionary process. We assess parallel adaptive radiation of South American pike cichlids (Crenicichla) using phylogenomics and phylogenetic comparative methods. We find that species flocks in the Uruguay and Iguazú River basins rapidly diversified into the same set of ecomorphs that reflect feeding ecology. Both adaptive radiations involve expansion of functional morphology, resulting in unique jaw phenotypes. Yet, form and function were decoupled such that most ecomorphs share similar mechanical properties of the jaws (i.e., jaw motion during a feeding strike). Prey mobility explained six to nine-fold differences in the rate of morphological evolution, but had no effect on the rate of mechanical evolution. We find no evidence of gene flow between species flocks or with surrounding coastal lineages that may explain their rapid diversification. When compared to cichlids of the East African Great Lakes and other prominent adaptive radiations, pike cichlids share many themes, including rapid expansion of phenotypic diversity, specialization along the benthic-to-pelagic habitat and soft-to-hard prey axes, and the evolution of conspicuous functional innovations. Yet, decoupled evolution of form and function and the absence of hybridization as a catalyzing force are departures from patterns observed in other adaptive radiations. Many-to-one mapping of morphology to mechanical properties is a mechanism by which pike cichlids exhibit a diversity of feeding ecologies while avoiding exacerbating underlying mechanical trade-offs.

Phylogenomics of bonytongue fishes (Osteoglossomorpha) shed light on the craniofacial evolution and biogeography of the weakly electric clade Mormyridae

Bonytongues (Osteoglossomorpha) constitute an ancient clade of teleost fishes distributed in freshwater habitats throughout the world. The group includes well-known species such as arowanas, featherbacks, pirarucus, and the weakly electric fishes in the family Mormyridae. Their disjunct distribution, extreme morphologies, and electrolocating capabilities (Gymnarchidae and Mormyridae) have attracted much scientific interest, but a comprehensive phylogenetic framework for comparative analysis is missing, especially for the species-rich family Mormyridae. Of particular interest are disparate craniofacial morphologies among mormyrids which might constitute an exceptional model system to study convergent evolution. We present a phylogenomic analysis based on 546 exons of 179 species (out of 260), 28 out of 29 genera, and all six families of extant bonytongues. Based on a recent reassessment of the fossil record of osteoglossomorphs, we inferred dates of divergence among trans-continental clades and the major groups. The estimated ages of divergence among extant taxa (e.g., Osteoglossomorpha, Osteoglossiformes, Mormyroidea) are older than previous reports, but most of the divergence dates obtained for clades on separate continents are too young to be explained by simple vicariance hypotheses. Biogeographic analysis of mormyrids indicates that their high species diversity in the Congo Basin is a consequence of range reductions of previously widespread ancestors and that the highest diversity of craniofacial morphologies among mormyrids originated in this basin. Special emphasis on a taxon-rich representation for mormyrids revealed pervasive misalignment between our phylogenomic results and mormyrid taxonomy due to repeated instances of convergence for extreme craniofacial morphologies. Estimation of ancestral phenotypes revealed contingent evolution of snout elongation and unique projections from the lower jaw to form the distinctive Schnauzenorgan. Synthesis of comparative analyses suggests that the remarkable craniofacial morphologies of mormyrids evolved convergently due to niche partitioning, likely enabled by interactions between their exclusive morphological and electrosensory adaptations.

The evolution of comparative phylogeography: putting the geography (and more) into comparative population genomics

Comparative population genomics is an ascendant field using genomic comparisons between species to draw inferences about forces regulating genetic variation. Comparative phylogeography, by contrast, focuses on the shared lineage histories of species codistributed geographically and is decidedly organismal in perspective. Comparative phylogeography is approximately 35 years old, and, by some metrics, is showing signs of reduced growth. Here, we contrast the goals and methods of comparative population genomics and comparative phylogeography and argue that comparative phylogeography offers an important perspective on evolutionary history that succeeds in integrating genomics with landscape evolution in ways that complement the suprageographic perspective of comparative population genomics. Focusing primarily on terrestrial vertebrates, we review the history of comparative phylogeography, its milestones and ongoing conceptual innovations, its increasingly global focus, and its status as a bridge between landscape genomics and the process of speciation. We also argue that, as a science with a strong “sense of place,” comparative phylogeography offers abundant “place-based” educational opportunities with its focus on geography and natural history, as well as opportunities for collaboration with local communities and indigenous peoples. Although comparative phylogeography does not yet require whole-genome sequencing for many of its goals, we conclude that it nonetheless plays an important role in grounding our interpretation of genetic variation in the fundamentals of geography and Earth history.

Comparing diversification rates in lakes, rivers, and the sea

The diversity of species inhabiting freshwater relative to marine habitats is striking, given that freshwater habitats encompass <1% of Earth's water. The most commonly proposed explanation for this pattern is that freshwater habitats are more fragmented than marine habitats, allowing more opportunities for allopatric speciation and thus increased diversification rates in freshwater. However, speciation may be generally faster in sympatry than in allopatry, as illustrated by lacustrine radiations such as African cichlids. Such differences between rivers and lakes may be important to consider when comparing diversification broadly among freshwater and marine groups. Here I compared diversification rates of teleost fishes in marine, riverine and lacustrine habitats. I found that lakes had faster speciation and net diversification rates than other aquatic habitats. However, most freshwater diversity arose in rivers. Surprisingly, riverine and marine habitats had similar rates of net diversification on average. Biogeographic models suggest that lacustrine habitats are evolutionarily unstable, explaining the dearth of lacustrine species in spite of their rapid diversification. Collectively, these results suggest that strong diversification rate differences are unlikely to explain the freshwater paradox. Instead, this pattern may be attributable to the comparable amount of time spent in riverine and marine habitats over the 200-million-year history of teleosts.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

Impact of Pleistocene Eustatic Fluctuations on Evolutionary Dynamics in Southeast Asian Biodiversity Hotspots

Pleistocene Climatic Fluctuations (PCF) are frequently highlighted as important evolutionary engines that triggered cycles of biome expansion and contraction. While there is ample evidence of the impact of PCF on biodiversity of continental biomes, the consequences in insular systems depend on the geology of the islands and the ecology of the taxa inhabiting them. The idiosyncratic aspects of insular systems are exemplified by the islands of the Sunda Shelf in Southeast Asia (Sundaland), where PCF-induced eustatic fluctuations had complex interactions with the geology of the region, resulting in high species diversity and endemism. Emergent land in Southeast Asia varied drastically with sea level fluctuations during the Pleistocene. Climate-induced fluctuations in sea level caused temporary connections between insular and continental biodiversity hotspots in Southeast Asia. These exposed lands likely had freshwater drainage systems that extended between modern islands: the Paleoriver Hypothesis. Built upon the assumption that aquatic organisms are among the most suitable models to trace ancient river boundaries and fluctuations of landmass coverage, the present study aims to examine the evolutionary consequences of PCF on the dispersal of freshwater biodiversity in Southeast Asia. Time-calibrated phylogenies of DNA-delimited species were inferred for six species-rich freshwater fish genera in Southeast Asia (Clarias, Channa, Glyptothorax, Hemirhamphodon, Dermogenys, Nomorhamphus). The results highlight rampant cryptic diversity and the temporal localization of most speciation events during the Pleistocene, with 88% of speciation events occurring during this period. Diversification analyses indicate that sea level-dependent diversification models poorly account for species proliferation patterns for all clades excepting Channa. Ancestral area estimations point to Borneo as the most likely origin for most lineages, with two waves of dispersal to Sumatra and Java during the last 5 Myr. Speciation events are more frequently associated with boundaries of the paleoriver watersheds, with 60%, than islands boundaries, with 40%. In total, one-third of speciation events are inferred to have occured within paleorivers on a single island, suggesting that habitat heterogeneity and factors other than allopatry between islands substantially affected diversification of Sundaland fishes. Our results suggest that species proliferation in Sundaland is not wholly reliant on Pleistocene sea-level fluctuations isolating populations on different islands.

Following the Mangroves: diversification in the banded lampeye Aplocheilichthys spilauchen (Duméril, 1861) (Cyprinodontiformes: Procatopodidae) along the Atlantic coast of Africa

Available ecological information, an extensive distributional range, conflicting osteological data, and a proposed early Miocene origin provide the impetus for the present study which investigates genetic structuring, biogeographic, and phylogenetic relationships within the Aplocheilichthys spilauchen lineage. Through the analysis of the mitochondrial gene COI, species delimitation methods (ABGD, GB, GMYC, bPTP) were applied, recognizing 6-7 OTUs with absolute pairwise genetic distances ranging between 8 and 22%. The onset of diversification is estimated to be within the middle Miocene and both dispersal and vicariance-shaped A. spilauchen diversity and distribution, as suggested by time-calibrated and ancestral range reconstruction (S-DIVA) analyses. We report for the first time, a pattern of diversification within a lineage of brackish water fish that is concordant with the historical distribution of coastal mangroves forests, shaped by a series of historical events that likely affected forest cover since the middle Miocene (e.g. major climate shifts and sea-level fluctuations, onset of the modern Congo River outlet, increased volcanism in the Cameroon Volcanic Line).

SUPPLEMENTARY INFORMATION

The online version of this article (doi:10.1007/s10750-020-04497-3) contains supplementary material, which is available to authorized users.

Genetic Variation and Hybridization in Evolutionary Radiations of Cichlid Fishes

Evolutionary radiations are responsible for much of the variation in biodiversity across taxa. Cichlid fishes are well known for spectacular evolutionary radiations, as they have repeatedly evolved into large and phenotypically diverse arrays of species. Cichlid genomes carry signatures of past events and, at the same time, are the substrate for ongoing evolution. We survey genome-wide data and the available literature covering 438 cichlid populations (412 species) across multiple radiations to synthesize information about patterns and sharing of genetic variation. Nucleotide diversity within species is low in cichlids, with 92% of surveyed populations having less diversity than the median value found in other vertebrates. Divergence within radiations is also low, and a large proportion of variation is shared among species due to incomplete lineage sorting and widespread hybridization. Population genetics therefore provides a suitable conceptual framework for evolutionary genomic studies of cichlid radiations. We focus in detail on the roles of hybridization in shaping the patterns of genetic variation and in promoting cichlid diversification.

Population structure of Apodemus flavicollis and comparison to Apodemus sylvaticus in northern Poland based on RAD-seq

Background

Mice of the genus Apodemus are one the most common mammals in the Palaearctic region. Despite their broad range and long history of ecological observations, there are no whole-genome data available for Apodemus, hindering our ability to further exploit the genus in evolutionary and ecological genomics context.

Results

Here we present results from the double-digest restriction site-associated DNA sequencing (ddRAD-seq) on 72 individuals of A. flavicollis and 10 A. sylvaticus from four populations, sampled across 500 km distance in northern Poland. Our data present clear genetic divergence of the two species, with average p-distance, based on 21377 common loci, of 1.51% and a mutation rate of 0.0011 - 0.0019 substitutions per site per million years. We provide a catalogue of 117 highly divergent loci that enable genetic differentiation of the two species in Poland and to a large degree of 20 unrelated samples from several European countries and Tunisia. We also show evidence of admixture between the three A. flavicollis populations but demonstrate that they have negligible average population structure, with largest pairwise F_ST<0.086.

Conclusion

Our study demonstrates the feasibility of ddRAD-seq in Apodemus and provides the first insights into the population genomics of the species.

Local adaptation fuels cryptic speciation in terrestrial annelids

Uncovering the genetic and evolutionary basis of cryptic speciation is a major focus of evolutionary biology. Next Generation Sequencing (NGS) allows the identification of genome-wide local adaptation signatures, but has rarely been applied to cryptic complexes - particularly in the soil milieu - as it is the case with integrative taxonomy. The earthworm genus Carpetania, comprising six previously suggested putative cryptic lineages, is a promising model to study the evolutionary phenomena shaping cryptic speciation in soil-dwelling lineages. Genotyping-By-Sequencing (GBS) was used to provide genome-wide information about genetic variability between 17 populations, and geometric morphometrics analyses of genital chaetae were performed to investigate unexplored cryptic morphological evolution. Genomic analyses revealed the existence of three cryptic species, with half of the previously-identified potential cryptic lineages clustering within them. Local adaptation was detected in more than 800 genes putatively involved in a plethora of biological functions (most notably reproduction, metabolism, immunological response and morphogenesis). Several genes with selection signatures showed shared mutations for each of the cryptic species, and genes under selection were enriched in functions related to regulation of transcription, including SNPs located in UTR regions. Finally, geometric morphometrics approaches partially confirmed the phylogenetic signal of relevant morphological characters such as genital chaetae. Our study therefore unveils that local adaptation and regulatory divergence are key evolutionary forces orchestrating genome evolution in soil fauna.

Genomic and mitochondrial evidence of ancient isolations and extreme introgression in the four-lined snake

Comparing mitochondrial and genomic phylogenies is an essential tool for investigating speciation processes, because each genome carries different inheritance properties and evolutionary characteristics. Furthermore, mitonuclear discordance may arise from ecological adaptation, historic isolation, population size changes, and sex-biased dispersal. Closely related taxa are expected to experience gene flow; however, this may not be true for insular populations or populations isolated in refugia. The four-lined snake Elaphe quatuorlineata has a fragmented distribution, separating populations of the Italian and Balkan Peninsulas, whereas several insular Aegean populations of significantly smaller body size (Cyclades island group and Skyros Island, Greece) are currently considered distinct subspecies. We constructed the species-tree phylogeny of this species utilizing genome-wide single nucleotide polymorphisms and a gene-tree based on complete cytochrome b sequences, aiming to detect convergence and discrepancies between biparentally and maternally inherited genomes. Population structuring, phylogenetic patterns and migration events among geographically defined lineages supported our hypothesis of isolation in multiple sub-refugia. Where biogeographical barriers did not restrict subsequent dispersal, extensive genetic exchange occurred between mainland Balkan populations. This process has led to the mitochondrial sweep of an ancestral mitolineage that survived only in peripheral (East Greece) and insular populations (North Cyclades and Skyros). The Central Cyclades represent an ancient lineage for both molecular markers that emerged almost 3.3 Mya. Considering their distinct morphology, insular E. quatuorlineata populations should be the future focus of an extensive sampling, especially since the mitonuclear discordance observed in this species could be related to ecological adaptations, such as the island-dwarfism phenomenon.

Phylogenomics of pike cichlids (Cichlidae: Crenicichla): the rapid ecological speciation of an incipient species flock

The rapid rise of phenotypic and ecological diversity in independent lake-dwelling groups of cichlids is emblematic of the East African Great Lakes. In this study, we show that similar ecologically based diversification has occurred in pike cichlids (Crenicichla) throughout the Uruguay River drainage of South America. We collected genomic data from nearly 500 ultraconserved element (UCEs) loci and >260 000 base pairs across 33 species, to obtain a phylogenetic hypothesis for the major species groups and to evaluate the relationships and genetic structure among five closely related, endemic, co-occurring species (the Uruguay River species flock; URSF). Additionally, we evaluated ecological divergence of the URSF based on body and lower pharyngeal jaw (LPJ) shape and gut contents. Across the genus, we recovered novel relationships among the species groups. We found strong support for the monophyly of the URSF; however, relationships among these species remain problematic, likely because of the rapid and recent evolution of this clade. Clustered co-ancestry analysis recovered most species as well delimited genetic groups. The URSF species exhibit species-specific body and LPJ shapes associated with specialized trophic roles. Collectively, our results suggest that the URSF consists of incipient species that arose via ecological speciation associated with the exploration of novel trophic roles.

Shift from slow- to fast-water habitats accelerates lineage and phenotype evolution in a clade of Neotropical suckermouth catfishes (Loricariidae: Hypoptopomatinae)

Identifying habitat characteristics that accelerate organismal evolution is essential to understanding both the origins of life on Earth and the ecosystem properties that are most critical to maintaining life into the future. Searching for these characteristics on a large scale has only recently become possible via advances in phylogenetic reconstruction, time-calibration, and comparative analyses. In this study, we combine these tools with habitat and phenotype data for 105 species in a clade of Neotropical suckermouth catfishes commonly known as cascudinhos. Our goal was to determine whether riverine mesohabitats defined by different flow rates (i.e., pools vs. rapids) and substrates (plants vs. rocks) have affected rates of cascudinho cladogenesis and morphological diversification. In contrast to predictions based on general theory related to life in fast-flowing, rocky riverine habitats, Neoplecostomini lineages associated with these habitats exhibited increased body size, head shape diversity, and lineage and phenotype diversification rates. These findings are consistent with a growing understanding of river rapids as incubators of biological diversification and specialization. They also highlight the urgent need to conserve rapids habitats throughout the major rivers of the world.