Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers

Phylogeographic patterns of Cyphocharax from trans-Andean rivers and northward expansion to lower Central America (Teleostei, Curimatidae)

Phylogenetic analyses of mitochondrial and nuclear data of 31 specimens of Cyphocharax from trans-Andean rivers support the presence of one lineage of Cyphocharax aspilos in Lago Maracaibo and three cryptic lineages of Cyphocharax magdalenae: (1) Cauca-Magdalena and Ranchería, (2) León and Atrato, and (3) Chucunaque-Tuira, Santa María, and Chiriquí basins of Central America. Results suggest that the Serranía del Perijá facilitated Late Miocene cladogenetic events, whereas post-Isthmian C. magdalenae expansion was enabled by gene flow across the lower Magdalena valley and Central American lowlands. Time-calibrated phylogenetics indicate that the C. magdalenae colonized lower Central America in the Pliocene (3.7 MYA; Ma), the divergence Atrato-Magdalena occurred in Late Pliocene (3.0 Ma) and the split Ranchería-Magdalena during the Middle Pleistocene (1.3 Ma). Updated geographic distribution data support the hypothesis that the Cordillera de Talamanca functions as a barrier to northward expansion of C. magdalenae in Central America.

New records of native and introduced fish species in a river basin of Western Ecuador, the Chocó-Darien Ecoregion, using DNA barcoding

DNA barcoding, based on mitochondrial markers, is widely applied in species identification and biodiversity studies. The aim of this study was to establish a barcoding reference database of fishes inhabiting the Cube River from Western Ecuador in the Chocó-Darien Global Ecoregion (CGE), a threatened ecoregion with high diversity and endemism, and evaluate the applicability of using barcoding for the identification of fish species. Barcode sequences were obtained from seven orders, 17 families, 23 genera and 26 species, which were validated through phylogenetic analysis, morphological measurements, and literature review. Our results showed that 43% of fish species in this region are endemic, confirmed the presence of known species in the area, and included the addition of three new records of native (Hoplias microlepis, Rhamdia guatemalensis and Sicydium salvini) and an introduced species (Xiphophorus maculatus) to Ecuador. In addition, eight species were barcoded for the first time. Species identification based on barcoding and morphology showed discrepancy with species lists from previous studies in the CGE, suggesting that the current baseline of western fishes of Ecuador is still incomplete. Because this study analyzed fishes from a relatively small basin (165 km2), more molecular-based studies focusing on fish are needed to achieve a robust sequence reference library of species inhabiting Western Ecuador. The new sequences of this study will be useful for future comparisons and biodiversity monitoring, supporting the application of barcoding tools for studying fish diversity in genetically unexplored regions and to develop well-informed conservation programs.

Phylogenomics and classification of Notropis and related shiners (Cypriniformes: Leuciscidae) and the utility of exon capture on lower taxonomic groups

North American minnows of the Shiner Clade, within the family Leuciscidae, represent one of the most taxonomically complex clades of the order Cypriniformes due to the large number of taxa coupled with conserved morphologies. Species within this clade were moved between genera and subgenera until the community decided to lump many of the unclassified taxa with similar morphologies into one genus, Notropis, which has held up to 325 species. Despite phylogentic studies that began to re-elevate some genera merged into Notropis, such as Cyprinella, Luxilus, Lythrurus, and Pteronotropis, the large genus Notropis remained as a taxonomic repository for many shiners of uncertain placement. Recent molecular advances in sequencing technologies have provided the opportunity to re-examine the Shiner Clade using phylogenomic markers. Using a fish probe kit, we sequenced 90 specimens in 87 species representing 16 genera included in the Shiner Clade, with a resulting dataset of 1,004 loci and 286,455 base pairs. Despite the large dataset, only 32,349 bp (11.29%) were phylogenetically informative. In our maximum likelihood tree, 78% of nodes are 100% bootstrap supported demonstrating the utility of the phylogenomic markers at lower taxonomic levels. Unsurprisingly, species within Notropis as well as Hudsonius, Luxilus, and Alburnops are not resolved as monophyletic groups. Cyprinella is monophyletic if Cyprinella callistia is excluded, and Pteronotropis is monophyletic if it includes Hudsonius cummingsae. Taxonomic changes we propose are: restriction of species included in Alburnops and Notropis, elevation of the subgenus Hydrophlox, expansion of species included in Miniellus, movement of Hudsonius cummingsae to Pteronotropis, and resurrection of the genera Coccotis and Paranotropis. We additionally had two specimens of three species, Notropis atherinoides, Ericymba amplamala, and Pimephales vigilax and found signficant differences between the localities (1,086, 1,424, and 845 nucleotides respectively).

Demographic inference provides insights into the extirpation and ecological dominance of eusocial snapping shrimps

Although eusocial animals often achieve ecological dominance in the ecosystems where they occur, many populations are unstable, resulting in local extinction. Both patterns may be linked to the characteristic demography of eusocial species-high reproductive skew and reproductive division of labor support stable effective population sizes that make eusocial groups more competitive in some species, but also lower effective population sizes that increase susceptibility to population collapse in others. Here, we examine the relationship between demography and social organization in Synalpheus snapping shrimps, a group in which eusociality has evolved recently and repeatedly. We show using coalescent demographic modelling that eusocial species have had lower but more stable effective population sizes across 100,000 generations. Our results are consistent with the idea that stable population sizes may enable competitive dominance in eusocial shrimps, but they also suggest that recent population declines are likely caused by eusocial shrimps' heightened sensitivity to environmental changes, perhaps as a result of their low effective population sizes and localized dispersal. Thus, although the unique life histories and demography of eusocial shrimps have likely contributed to their persistence and ecological dominance over evolutionary timescales, these social traits may also make them vulnerable to contemporary environmental change.

Integrating multi-taxon palaeogenomes and sedimentary ancient DNA to study past ecosystem dynamics

Ancient DNA (aDNA) has played a major role in our understanding of the past. Important advances in the sequencing and analysis of aDNA from a range of organisms have enabled a detailed understanding of processes such as past demography, introgression, domestication, adaptation and speciation. However, to date and with the notable exception of microbiomes and sediments, most aDNA studies have focused on single taxa or taxonomic groups, making the study of changes at the community level challenging. This is rather surprising because current sequencing and analytical approaches allow us to obtain and analyse aDNA from multiple source materials. When combined, these data can enable the simultaneous study of multiple taxa through space and time, and could thus provide a more comprehensive understanding of ecosystem-wide changes. It is therefore timely to develop an integrative approach to aDNA studies by combining data from multiple taxa and substrates. In this review, we discuss the various applications, associated challenges and future prospects of such an approach.

The unknown diversity of the genus Characidium (Characiformes: Crenuchidae) in the Chocó biogeographic region, Colombian Andes: Two new species supported by morphological and molecular data

Two new species of Characidium, C. tatama and C. dule, are described from the biogeographic Chocó region in western Colombia. Both new species are supported by both morphological and molecular data. C. tatama from the San Juan River in the Pacific and C. dule from the Atrato River in the Caribbean portion of Colombia are both distributed in the upper and lower portions of these basins. An extensive comparison with other trans- and cis-Andean species of Characidium was made, in addition to species delimitation, using COI sequences by distinct methods (GMYC, ABGD, bPTP).

Exon probe sets and bioinformatics pipelines for all levels of fish phylogenomics

Exon markers have a long history of use in phylogenetics of ray-finned fishes, the most diverse clade of vertebrates with more than 35,000 species. As the number of published genomes increases, it has become easier to test exons and other genetic markers for signals of ancient duplication events and filter out paralogues that can mislead phylogenetic analysis. We present seven new probe sets for current target-capture phylogenomic protocols that capture 1,104 exons explicitly filtered for paralogues using gene trees. These seven probe sets span the diversity of teleost fishes, including four sets that target five hyperdiverse percomorph clades which together comprise ca. 17,000 species (Carangaria, Ovalentaria, Eupercaria, and Syngnatharia + Pelagiaria combined). We additionally included probes to capture legacy nuclear exons and mitochondrial markers that have been commonly used in fish phylogenetics (despite some exons being flagged for paralogues) to facilitate integration of old and new molecular phylogenetic matrices. We tested these probes experimentally for 56 fish species (eight species per probe set) and merged new exon-capture sequence data into an existing data matrix of 1,104 exons and 300 ray-finned fish species. We provide an optimized bioinformatics pipeline to assemble exon capture data from raw reads to alignments for downstream analysis. We show that legacy loci with known paralogues are at risk of assembling duplicated sequences with target-capture, but we also assembled many useful orthologous sequences that can be integrated with many PCR-generated matrices. These probe sets are a valuable resource for advancing fish phylogenomics because targeted exons can easily be extracted from increasingly available whole genome and transcriptome data sets, and also may be integrated with existing PCR-based exon and mitochondrial data.

Divergence, gene flow, and speciation in eight lineages of trans-Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans-Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split-migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST ), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.

Genetic Diversity of the Cichlid Andinoacara latifrons (Steindachner, 1878) as a Conservation Strategy in Different Colombian Basins

Cichlids constitute a diverse monophyletic group that have developed adaptive strategies to thrive in diverse environments. Andinoacara represents an example of diversification on the South American Andean uplift, providing a key model for understanding the evolution of biogeographic patterns. In this study, we analyzed the species Andinoacara latifrons using two mitochondrial markers (COI, cytb) and one nuclear marker (RAG1) in a populational level. Sequences were obtained through tissue collection and from the GenBank database. Populational analysis showed significant structuration among populations, also corroborated with population pairwise F_st results. Fu’s Fs and Tajima’s D results showed populations that seems to be under populational expansion. We identified 22 haplotypes using cytb. The population associations in the Cauca haplotype are related to the Momposina depression and the mixture of the Cauca-Magdalena river basins in the lower Cauca-Magdalena region. We constructed a new phylogenetic tree, which grouped mainly two A. latifrons lineages: (1) an upper Magdalena and Catatumbo clade and (2) an upper Cauca and upper Magdalena clade. Thus, A. latifrons represents a diverse entity that contributes to our understanding of the evolutionary history of northern South America. Our findings provide insight into devising public policies in determining refuges for the preservation of biodiversity in the lower Cauca and Magdalena regions in Colombia.

Complete Mitogenomes of Three Carangidae (Perciformes) Fishes: Genome Description and Phylogenetic Considerations

Carangidae are ecologically and economically important marine fish. The complete mitogenomes of three Carangidae species (Alectis indicus, Decapterus tabl, and Alepes djedaba) were sequenced, characterized, and compared with 29 other species of the family Carangidae in this study. The length of the three mitogenomes ranged from 16,530 to 16,610 bp, and the structures included 2 rRNA genes (12S rRNA and 16S rRNA), 1 control region (a non-coding region), 13 protein-coding genes, and 22 tRNA genes. Among the 22 tRNA genes, only tRNA-Ser (GCT) was not folded into a typical cloverleaf secondary structure and had no recognizable DHU stem. The full-length sequences and protein-coding genes (PCGs) of the mitogenomes of the three species all had obvious AT biases. The majority of the AT-skew and GC-skew values of the PCGs among the three species were negative, demonstrating bases T and C were more plentiful than A and G. Analyses of Ka/Ks and overall p-genetic distance demonstrated that ATP8 showed the highest evolutionary rate and COXI/COXII were the most conserved genes in the three species. The phylogenetic tree based on PCGs sequences of mitogenomes using maximum likelihood and Bayesian inference analyses showed that three clades were divided corresponding to the subfamilies Caranginae, Naucratinae, and Trachinotinae. The monophyly of each superfamily was generally well supported. The divergence time analyses showed that Carangidae evolved during three geological periods, the Cretaceous, Paleogene, and Neogene. A. indicus began to differentiate from other species about 27.20 million years ago (Mya) in the early Miocene, while D. tabl (21.25 Mya) and A. djedaba (14.67 Mya) differentiated in the middle Oligocene.

Ancient genetic divergence in bumblebee catfish of the genus Pseudopimelodus (Pseudopimelodidae: Siluriformes) from northwestern South America

Pseudopimelodus is a Neotropical genus of bumblebee catfish, composed of four valid species occurring in both trans- and cis-Andean rivers of South America. The orogeny of the Andes has led to diversification in the genus Pseudopimelodus in Colombia. This study analyzed partial sequences of mitochondrial cox1 and nuclear rag2 genes to test the hypothesis that the species, nominally recognized as P. schultzi and P. bufonius in Colombia, correspond to more than two different evolutionary lineages. Results indicate high levels of genetic divergence among individuals of nominal P. schultzi and P. bufonius, from trans- and cis-Andean basins in Colombia. In addition, five divergent lineages of Pseudopimelodus were confidently delimited by using a single-locus species-discovery approach and confirmed by species tree analyses. Additionally, molecular-clock dating showed that most diversification processes in Pseudopimelodus took place during the Miocene, when Andean tectonic evolution was occurring in northwestern South America. The present study provides, for the first time, phylogeographic insight into this Neotropical genus.

Comparative genomics approach to evolutionary process connectivity

The influence of species life history traits and historical demography on contemporary connectivity is still poorly understood. However, these factors partly determine the evolutionary responses of species to anthropogenic landscape alterations. Genetic connectivity and its evolutionary outcomes depend on a variety of spatially dependent evolutionary processes, such as population structure, local adaptation, genetic admixture, and speciation. Over the last years, population genomic studies have been interrogating these processes with increasing resolution, revealing a large diversity of species responses to spatially structured landscapes. In parallel, multispecies meta-analyses usually based on low-genome coverage data have provided fundamental insights into the ecological determinants of genetic connectivity, such as the influence of key life history traits on population structure. However, comparative studies still lack a thorough integration of macro- and micro-evolutionary scales to fully realize their potential. Here, I present how a comparative genomics framework may provide a deeper understanding of evolutionary process connectivity. This framework relies on coupling the inference of long-term demographic and selective history with an assessment of the contemporary consequences of genetic connectivity. Standardizing this approach across several species occupying the same landscape should help understand how spatial environmental heterogeneity has shaped the diversity of historical and contemporary connectivity patterns in different taxa with contrasted life history traits. I will argue that a reasonable amount of genome sequence data can be sufficient to resolve and connect complex macro- and micro-evolutionary histories. Ultimately, implementing this framework in varied taxonomic groups is expected to improve scientific guidelines for conservation and management policies.

Fine-Scale Population Genetic Structure and Parapatric Cryptic Species of Kuruma Shrimp (Marsupenaeus japonicus), Along the Northwestern Pacific Coast of China

The kuruma shrimp (Marsupenaeus japonicus) includes two cryptic species, which are distributed mostly allopatrically but co-occur in the northern South China Sea (from Huilai to Beihai). To obtain a better understanding of the fine-scale genetic structure and parapatric diversification of these two varieties in the northwestern Pacific region, we used a genotyping-by-sequencing (GBS) and comparative transcriptomics approach to establish their phylogenetic relationships. Using the GBS technique, we genotyped 28891 SNPs in 160 individuals in the Northwest Pacific. The results supported two highly diverged evolutionary lineages of kuruma shrimp (var. I and II). The ND and XM populations showed complex genetic patterns, which might be affected by the complex environment of the Taiwan Strait. In addition, the migration rates and inbreeding coefficients of XM and BH were much lower than those of the other populations, which might be related to the land-sea changes and complex ocean currents in the Taiwan Strait and Qiongzhou Strait. Based on the synonymous substitution rates (ds) of 2,491 candidate orthologs, we estimated that the divergence time between the two varieties was 0.26~0.69 Mya. Choice and no-choice interbreeding experiments provided support for the biological species concept, by showing the existence of reproductive isolation or incompatibility. In view of these differences between the two Marsupenaeus species, we believe that it is essential and urgent to establish a genetic database for each and reevaluate their ecological suitable conditions in order to improve species-specific culturing techniques. Moreover, this research can serve as a case study for future research on speciation and hybridization.

Species delimitation reveals an underestimated diversity of Andean catfishes of the family Astroblepidae (Teleostei: Siluriformes)

ABSTRACT Catfishes of the family Astroblepidae form a group composed by 82 valid species of the genus Astroblepus inhabiting high-gradient streams and rivers throughout tropical portions of the Andean Cordillera. Little has been advanced in the systematics and biodiversity of astroblepids other than an unpublished thesis, a single regional multilocus study and isolated species descriptions. Here, we examined 208 specimens of Astroblepus that apparently belong to 16 valid species from several piedmont rivers from northern Colombia to southern Peru. Using three single-locus approaches for species delimitation in combination with a species tree analysis estimated from three mitochondrial genes, we identified a total of 25 well-delimited lineages including eight valid and 17 potential undescribed species distributed in two monophyletic groups: the Central Andes Clade, which contains 14 lineages from piedmont rivers of the Peruvian Amazon, and the Northern Andes Clade with 11 lineages from trans- and cis-Andean rivers of Colombia and Ecuador, including the Orinoco, Amazon, and Magdalena-Cauca basins and Pacific coastal drainages. Results of species delimitation methods highlight several taxonomical incongruences in recently described species denoting potential synonymies.

Tempo and mode of allopatric divergence in the weakly electric fish Sternopygus dariensis in the Isthmus of Panama

Spatial isolation is one of the main drivers of allopatric speciation, but the extent to which spatially-segregated populations accumulate genetic differences relevant to speciation is not always clear. We used data from ultraconserved elements (UCEs) and whole mitochondrial genomes (i.e., mitogenomes) to explore genetic variation among allopatric populations of the weakly electric fish Sternopygus dariensis across the Isthmus of Panama. We found strong genetic divergence between eastern and western populations of S. dariensis. Over 77% of the UCE loci examined were differentially fixed between populations, and these loci appear to be distributed across the species' genome. Population divergence occurred within the last 1.1 million years, perhaps due to global glaciation oscillations during the Pleistocene. Our results are consistent with a pattern of genetic differentiation under strict geographic isolation, and suggest the presence of incipient allopatric species within S. dariensis. Genetic divergence in S. dariensis likely occurred in situ, long after the closure of the Isthmus of Panama. Our study highlights the contribution of spatial isolation and vicariance to promoting rapid diversification in Neotropical freshwater fishes. The study of spatially-segregated populations within the Isthmus of Panama could reveal how genetic differences accumulate as allopatric speciation proceeds.