Phylogenomic incongruence, hypothesis testing, and taxonomic sampling: The monophyly of characiform fishes

Phylogenomics of Neogastropoda: The Backbone Hidden in the Bush

The molluskan order Neogastropoda encompasses over 15,000 almost exclusively marine species playing important roles in benthic communities and in the economies of coastal countries. Neogastropoda underwent intensive cladogenesis in the early stages of diversification, generating a "bush" at the base of their evolutionary tree, which has been hard to resolve even with high throughput molecular data. In the present study to resolve the bush, we use a variety of phylogenetic inference methods and a comprehensive exon capture dataset of 1817 loci (79.6% data occupancy) comprising 112 taxa of 48 out of 60 Neogastropoda families. Our results show consistent topologies and high support in all analyses at (super)family level, supporting monophyly of Muricoidea, Mitroidea, Conoidea, and, with some reservations, Olivoidea and Buccinoidea. Volutoidea and Turbinelloidea as currently circumscribed are clearly paraphyletic. Despite our analyses consistently resolving most backbone nodes, 3 prove problematic: First, the uncertain placement of Cancellariidae, as the sister group to either a Ficoidea-Tonnoidea clade or to the rest of Neogastropoda, leaves monophyly of Neogastropoda unresolved. Second, relationships are contradictory at the base of the major "core Neogastropoda" grouping. Third, coalescence-based analyses reject monophyly of the Buccinoidea in relation to Vasidae. We analyzed phylogenetic signal of targeted loci in relation to potential biases, and we propose the most probable resolutions in the latter 2 recalcitrant nodes. The uncertain placement of Cancellariidae may be explained by orthology violations due to differential paralog loss shortly after the whole genome duplication, which should be resolved with a curated set of longer loci.

Genome-scale approach to reconstructing the phylogenetic tree of psyllids (superfamily Psylloidea) with account of systematic bias

Psyllids (class Insecta: order Hemiptera: superfamily Psylloidea) are a taxonomically and phylogenetically challenging clade. Recent studies have largely advanced the phylogeny of this group, yet the family-level relationships among Aphalaridae, Carsidaridae, and others remain unresolved. Genome-scale phylogenetic analysis is known to provide a finer resolution for problems like that. However, such phylogenomics also introduces new problems: incorrect trees with high confidence yielded due to systematic error (bias). Here we addressed these issues using hundreds of single-copy orthologous (SCO) genes in psyllid transcriptomes and genomes. Our analyses revealed conflicts between the nucleotide-based and amino-acid-based phylogenetic trees. While the nucleotide-based phylogeny strongly supported the (Aphalaridae + Carsidaridae) + Others relationship, the amino-acid-based one recovered Aphalaridae + (Carsidaridae + Others) with 100% support. Further inspection revealed significant compositional heterogeneity in nucleotide sequences for 67% of SCO genes, but not in the corresponding translated amino acid sequences. We then used different strategies to combat this compositional bias, and found that using the RY-coding strategy (coding the standard nucleotides as purines and pyrimidines) the nucleotide-based phylogeny became consistent with the amino-acid-based one. We further applied RY-coding to a published concatenated nucleotide dataset and recovered the Aphalaridae monophyly (which is refuted by the original literature on non-recoded sequences) at the base of psyllid tree. Moreover, it was found that variations in evolutionary rate could lead to errors in nucleotide-based phylogeny. The fast-evolving Heteropsylla cubana (Psyllidae: Ciriacreminae) was incorrectly placed within the subfamily Psyllinae. This bias can be avoided by using data removal or RY-coding strategies. Together, our results strongly support the family relationship of Aphalaridae + (Carsidaridae + Others), and show that the amino-acid-based concatenation analysis is more robust than nucleotide-based one. Future phylogenomic analysis of psyllid nucleotide sequences should take into account methods such as the RY-coding scheme to address potential systematic biases arising from composition and rate heterogeneities.

Unravelling the genomic features, phylogeny and genetic basis of tooth ontogenesis in Characiformes through analysis of four genomes

Characiformes is a diverse and evolutionarily significant order of freshwater fish encompassing over 2,300 species. Despite its diversity, our understanding of Characiformes' evolutionary relationships and adaptive mechanisms is limited due to insufficient genome sequences. In this study, we sequenced and assembled the genomes of four Characiformes species, three of which were chromosome-level assemblies. Our analyses revealed dynamic changes in gene family evolution, repeat sequences and variations in chromosomal collinearity within these genomes. With the assembled genomes, we were not only able to elucidate the evolutionary relationship of the four main orders in Otophysi but also indicated Characiformes as the paraphyletic group. Comparative genomic analysis with other available fish genomes shed light on the evolution of genes related to tooth development in Characiformes. Notably, variations in the copy number of secretory calcium-binding phosphoproteins (SCPP) genes were observed among different orders of Otophysi, indicating their potential contribution to the diversity of tooth types. Our study offers invaluable genome sequences and novel insights into Characiformes' evolution, paving the way for further genomic and evolutionary research in fish.

Major Revisions in Pancrustacean Phylogeny and Evidence of Sensitivity to Taxon Sampling

The clade Pancrustacea, comprising crustaceans and hexapods, is the most diverse group of animals on earth, containing over 80% of animal species and half of animal biomass. It has been the subject of several recent phylogenomic analyses, yet relationships within Pancrustacea show a notable lack of stability. Here, the phylogeny is estimated with expanded taxon sampling, particularly of malacostracans. We show small changes in taxon sampling have large impacts on phylogenetic estimation. By analyzing identical orthologs between two slightly different taxon sets, we show that the differences in the resulting topologies are due primarily to the effects of taxon sampling on the phylogenetic reconstruction method. We compare trees resulting from our phylogenomic analyses with those from the literature to explore the large tree space of pancrustacean phylogenetic hypotheses and find that statistical topology tests reject the previously published trees in favor of the maximum likelihood trees produced here. Our results reject several clades including Caridoida, Eucarida, Multicrustacea, Vericrustacea, and Syncarida. Notably, we find Copepoda nested within Allotriocarida with high support and recover a novel relationship between decapods, euphausiids, and syncarids that we refer to as the Syneucarida. With denser taxon sampling, we find Stomatopoda sister to this latter clade, which we collectively name Stomatocarida, dividing Malacostraca into three clades: Leptostraca, Peracarida, and Stomatocarida. A new Bayesian divergence time estimation is conducted using 13 vetted fossils. We review our results in the context of other pancrustacean phylogenetic hypotheses and highlight 15 key taxa to sample in future studies.

Phylogenomics, evolution of trophic traits and divergence times of hemiodontid fishes (Ostariophysi: Characiformes)

The South American characiform family Hemiodontidae comprises five genera and 34 species. The family lacks comprehensive phylogenetic hypotheses resolving its species relationships. The studies that addressed these questions exhibited a narrow taxon sampling or used single-locus markers. Herein we surveyed hundreds of ultraconserved elements (UCEs) loci to provide the first molecular phylogenetic hypothesis and divergence time estimates for hemiodontids encompassing all its genera and most species (27 of the 34 valid species). We also tracked the history of the protractile upper jaw in the genera Argonectes and Bivibranchia across the recovered phylogenies through ancestral state reconstruction. Our results corroborate the monophyly of Hemiodontidae and the genera Argonectes and Bivibranchia in all phylogenetic methods with maximum clade support. The genera Anodus and Hemiodus were not monophyletic because Anodus elongatus was sister to the monotypic Micromischodus instead of A. orinocensis, and H. immaculatus did not form a clade with its other congeners, but instead was sister to the clade including Anodus and Micromischodus. All remaining species of Hemiodus were placed together into a monophyletic group, where they were arranged into four major subclades. The relationship in the family is summarised as: (Bivibranchia, (Argonectes, ((H. immaculatus, (Anodus, Micromischodus)), Hemiodus clade))), in discordance with the morphological phylogeny that placed all genera monophyletic and resolved the family as: ((Anodus, Micromischodus), (Hemiodus, (Argonectes, Bivibranchia))). The origin of Hemiodontidae was estimated from the Late Cretaceous to the Middle Paleogene, with the mean age in the Paleocene, while the origin of most hemiodontid genera except Bivibranchia occurred in the Miocene. Unordered parsimony and likelihood reconstruction indicates that Argonectes and Bivibranchia developed their protractile upper jaw independently.

Statistical evaluation of character support reveals the instability of higher-level dinosaur phylogeny

The interrelationships of the three major dinosaur clades (Theropoda, Sauropodomorpha, and Ornithischia) have come under increased scrutiny following the recovery of conflicting phylogenies by a large new character matrix and its extensively modified revision. Here, we use tools derived from recent phylogenomic studies to investigate the strength and causes of this conflict. Using maximum likelihood as an overarching framework, we examine the global support for alternative hypotheses as well as the distribution of phylogenetic signal among individual characters in both the original and rescored dataset. We find the three possible ways of resolving the relationships among the main dinosaur lineages (Saurischia, Ornithischiformes, and Ornithoscelida) to be statistically indistinguishable and supported by nearly equal numbers of characters in both matrices. While the changes made to the revised matrix increased the mean phylogenetic signal of individual characters, this amplified rather than reduced their conflict, resulting in greater sensitivity to character removal or coding changes and little overall improvement in the ability to discriminate between alternative topologies. We conclude that early dinosaur relationships are unlikely to be resolved without fundamental changes to both the quality of available datasets and the techniques used to analyze them.

Multiple Outgroups Can Cause Random Rooting in Phylogenomics

Outgroup selection has been a major challenge since the rise of phylogenetics, and it has remained so in the phylogenomic era. Our goal here is to use large phylogenomic animal datasets to examine the impact of outgroup selection on the final topology. The results of our analyses further solidify the fact that distant outgroups can cause random rooting, and that this holds for concatenated and coalescent-based methods. The results also indicate that the standard practice of using multiple outgroups often causes random rooting. Most researchers go out of their way to get multiple outgroups, as this has been standard practice for decades. Based on our findings, this practice should stop. Instead, our results suggest that a single (most closely) related relative should be selected as the outgroup, unless all outgroups are roughly equally closely related to the ingroup.

Lip service: Histological phenotypes correlate with diet and feeding ecology in herbivorous pacus

Complex prey processing requires the repositioning of food between the teeth, as modulated by a soft tissue appendage like a tongue or lips. In this study, we trace the evolution of lips and ligaments, which are used during prey capture and prey processing in an herbivorous group of fishes. Pacus (Serrasalmidae) are Neotropical freshwater fishes that feed on leaves, fruits, and seeds. These prey are hard or tough, require high forces to fracture, contain abrasive or caustic elements, or deform considerably before failure. Pacus are gape-limited and do not have the pharyngeal jaws many bony fishes use to dismantle and/or transport prey. Despite their gape limitation, pacus feed on prey larger than their mouths, relying on robust teeth and a hypertrophied lower lip for manipulation and breakdown of food. We used histology to compare the lip morphology across 14 species of pacus and piranhas to better understand this soft tissue. We found that frugivorous pacus have larger, more complex lips which are innervated and folded at their surface, while grazing species have callused, mucus-covered lips. Unlike mammalian lips or tongues, pacu lips lack any intrinsic skeletal or smooth muscle. This implies that pacu lips lack dexterity; however, we found a novel connection to the primordial ligament which suggests that the lips are actuated by the jaw adductors. We propose that pacus combine hydraulic repositioning of prey inside the buccal cavity with direct oral manipulation, the latter using a combination of a morphologically heterodont dentition and compliant lips for reorienting food.

The Phylogenetic Relationships of the Family Sciaenidae Based on Genome-Wide Data Analysis

Uncertainty and controversy exist in the phylogenetic status of the Sciaenidae family because of the limited genetic data availability. In this study, a data set of 69,098 bp, covering 309 shared orthologous genes, was extracted from 18 genomes and 5 transcriptomes of 12 species belonging to the Sciaenidae family and used for phylogenetic analysis. The maximum likelihood (ML) and Bayesian approach (BA) methods were used to reconstruct the phylogenetic trees. The resolved ML and BA trees showed similar topology, thus revealing two major evolutionary lineages within the Sciaenidae family, namely, Western Atlantic (WA) and Eastern Atlantic−Indo−West Pacific (EIP). The WA group included four species belonging to four genera: Cynoscion nebulosus, Equetus punctatus, Sciaenops ocellatus, and Micropogonias undulatus. Meanwhile, the EIP group formed one monophyletic clade, harboring eight species (Argyrosomus regius, A. japonicus, Pennahia anea, Nibea albiflora, Miichthys miiuy, Collichthys lucidus, Larimichthys polyactis, and L. crocea) from six genera. Our results indicated that the Western Atlantic (WA) group was more ancient in the studied species, while the Eastern Atlantic−Indo−West Pacific (EIP) group was a younger group. Within the studied species, the genera Collichthys and Larmichthys were the youngest lineages, and we do not suggest that Collichthys and Larmichthys should be considered as one genus. However, the origin of the Sciaenidae family and problems concerning the basal genus were not resolved because of the lack of genomes. Therefore, further sampling and sequencing efforts are needed.

First cytogenetic data on Afrotropical lutefishes (Citharinidae) in the light of karyotype evolution in Characiformes

The Afrotropical lutefish family Citharinidae (Citharinoidei, Characiformes) comprises three genera with eight species in total. Although Citharinidae have been studied in terms of taxonomy and systematics, no cytogenetic information was available for any representative of the family. Furthermore, only one species out of 116 in Citharinoidei (Distichodusaffinis Günther, 1873) has been studied cytogenetically. Here, we report the karyotypes of Citharinuscitharus (Geoffroy St. Hilaire, 1809) from West Africa and Citharinuslatus Müller et Troschel, 1844 from Northeast Africa. The former has the diploid chromosome number 2n = 40 and the fundamental number FN = 80, while the latter has 2n = 44 and FN = 88. Hence, these karyotypes consist exclusively of bi-armed chromosomes. Such karyotypes were previously found in D.affinis and in many lineages of Neotropical species of another suborder of Characiformes, Characoidei. In contrast, the karyotypes dominated by uni-armed elements are typical for a number of phylogenetically basal lineages of Afrotropical and Neotropical Characoidei. We discuss the importance of our data on Citharinidae for the understanding of the karyotype evolution within the order Characiformes.

Detecting and removing sample contamination in phylogenomic data: an example and its implications for Cicadidae phylogeny (Insecta: Hemiptera)

Contamination of a genetic sample with DNA from one or more non-target species is a continuing concern of molecular phylogenetic studies, both Sanger sequencing studies and Next-Generation Sequencing (NGS) studies. We developed an automated pipeline for identifying and excluding likely cross-contaminated loci based on detection of bimodal distributions of patristic distances across gene trees. When the contamination occurs between samples within a dataset, comparisons between a contaminated sample and its contaminant taxon will yield bimodal distributions with one peak close to zero patristic distance. This new method does not rely on a priori knowledge of taxon relatedness nor does it determine the causes(s) of the contamination. Exclusion of putatively contaminated loci from a dataset generated for the insect family Cicadidae showed that these sequences were affecting some topological patterns and branch supports, although the effects were sometimes subtle, with some contamination-influenced relationships exhibiting strong bootstrap support. Long tip branches and outlier values for one anchored phylogenomic pipeline statistic (AvgNHomologs) were correlated with the presence of contamination. While the AHE markers used here, which target hemipteroid taxa, proved effective in resolving deep and shallow level Cicadidae relationships in aggregate, individual markers contained inadequate phylogenetic signal, in part probably due to short length. The cleaned dataset, consisting of 429 loci, from 90 genera representing 44 of 56 current Cicadidae tribes, supported three of the four sampled Cicadidae subfamilies in concatenated-matrix maximum likelihood (ML) and multispecies coalescent-based species tree analyses, with the fourth subfamily weakly supported in the ML trees. No well-supported patterns from previous family-level Sanger sequencing studies of Cicadidae phylogeny were contradicted. One taxon (Aragualna plenalinea) did not fall with its current subfamily in the genetic tree, and this genus and its tribe Aragualnini is reclassified to Tibicininae following morphological re-examination. Only subtle differences were observed in trees after removal of loci for which divergent base frequencies were detected. Greater success may be achieved by increased taxon sampling and developing a probe set targeting a more recent common ancestor and longer loci. Searches for contamination are an essential step in phylogenomic analyses of all kinds and our pipeline is an effective solution.

New phylogenetic insights into the African catfish families Mochokidae and Austroglanididae

Several hundred catfish species (order: Siluriformes) belonging to 11 families inhabit Africa, of which at least six families are endemic to the continent. Although four of those families are well-known to belong to the 'Big-Africa clade', no previous study has addressed the phylogenetic placement of the endemic African catfish family Austroglanididae in a comprehensive framework with molecular data. Furthermore, interrelationships within the 'Big-Africa clade', including the most diverse family Mochokidae, remain unclear. This study was therefore designed to help reconstruct inter- and intrarelationships of all currently valid mochokid genera, to infer their position within the 'Big Africa clade' and to establish a first molecular phylogenetic hypothesis of the relationships of the enigmatic Austroglanididae within the Siluriformes. We assembled a comprehensive mitogenomic dataset comprising all protein coding genes and representing almost all recognized catfish families (N = 33 of 39) with carefully selected species (N = 239). We recovered the monophyly of the previously identified multifamily clades 'Big Asia' and 'Big Africa' and determined Austroglanididae to be closely related to Pangasiidae, Ictaluroidea and Ariidae. Mochokidae was recovered as the sister group to a clade encompassing Auchenoglanididae, Claroteidae, Malapteruridae and the African Schilbeidae, albeit with low statistical support. The two mochokid subfamilies Mochokinae and Chiloglanidinae as well as the chiloglanid tribe Atopochilini were recovered as reciprocally monophyletic. The genus Acanthocleithron forms the sister group of all remaining Mochokinae, although with low support. The genus Atopodontus is the sister group of all remaining Atopochilini. In contrast to morphological reconstructions, the monophyly of the genus Chiloglanis was strongly supported in our analysis, with Chiloglanis macropterus nested within a Chiloglanis sublineage encompassing only other taxa from the Congo drainage. This is an important result because the phylogenetic relationships of C. macropterus have been controversial in the past, and because we and other researchers assumed that this species would be resolved as sister to most or all other members of Chiloglanis. The apparent paraphyly of Synodontis with respect to Microsynodontis provided an additional surprise, with Synodontis punu turning out to be the sister group of the latter genus.

Tracking the evolutionary pathways among Brazilian Lebiasina species (Teleostei: Lebiasinidae): a chromosomal and genomic comparative investigation

Abstract Despite several difficulties in chromosomal analyses of small-sized fishes, the cytogenetics of the Lebiasinidae was largely improved in the last years, showing differential patterns in the chromosomal evolution inside the family. In this context, it has been shown that genus Lebiasina preserves its karyotypic macrostructure, composed of 2n = 36 chromosomes, whereas the other genera generally present higher 2n. This study focused on the comparative cytogenetics of three Lebiasina species, one of them analyzed here for the first time, using conventional and molecular procedures. The results reinforced the differentiated evolutionary path of the genus Lebiasina while, at the same time, highlighted the genomic particularities that have accompanied the evolution of each species. In this sense, the repetitive components of the genome played a significant role in the differentiation of each species. It is also notable that L. minuta and L. melanoguttata, the two species that occur exclusively in the Brazilian territory, show greater chromosomal similarities to each other than to the trans-Andean sister species, L. bimaculata.

Phylogenetic relationships and description of two new species of Diapoma (Characidae: Stevardiinae) from the La Plata River basin

Abstract Herein we describe two new species of Diapoma, one from the Negro River, a tributary of the Uruguay River in Brazil and Uruguay, and one from the Iguaçu River, in Brazil and Argentina. The new species from the Negro River basin is distinguished from its congeners by the following combination of characters: a black narrow and conspicuous line restricted to the body horizontal septum, incomplete lateral line, tricuspid teeth in the inner series of the premaxilla, and a lower body depth at vertical through the dorsal-fin origin (29.3–32.8% SL in males and 27.7–33.3% SL in females). The new species from the Iguaçu River basin is distinguished from its congeners by the following combination of characters: a discontinuous lateral line, adipose fin hyaline, longer anal-fin base (26.5–32.4% SL), and a longitudinal black stripe along the median region of caudal-fin rays. Additionally, we updated the molecular phylogeny of the genus, including new sequences from these two new species and Diapoma thauma. An identification key for species of Diapoma is presented, modified from previous study.

Divergence times of the Rhoadsia clade (Characiformes: Characidae)

Abstract The family Characidae is the most diverse group of fishes in the Neotropics with challenging systematics. The three genera Carlana, Parastremma, and Rhoadsia, formerly considered the subfamily Rhoadsiinae, are now included in the subfamily Stethaprioninae. Previous phylogenetic analyses did not include all genera of Rhoadsiinae, specifically Parastremma. Here, we estimated the phylogenetic relationships and divergence times of the genera of Rhoadsiinae (the Rhoadsia clade) relative to the most representative genera of the Characidae. We used six molecular markers from the mitochondrial and nuclear genome to estimate the phylogeny and divergence times. We confirmed the monophyly of the Rhoadsia clade. Furthermore, we estimated that the Central American genus Carlana and the western Colombian genus Parastremma diverged approximately 13 Mya (95% HPD 8.36–18.11), consistent with the early-closure estimates of the Isthmus of Panama (~15 Mya). The genus Rhoadsia, endemic to Western Ecuador and Northern Peru, was estimated to originate at around 20 Mya (95% HPD 14.35–25.43), consistent with the Andean uplift (~20 Mya).

Tracking the Evolutionary Trends Among Small-Size Fishes of the Genus Pyrrhulina (Characiforme, Lebiasinidae): New Insights From a Molecular Cytogenetic Perspective

Miniature fishes have always been a challenge for cytogenetic studies due to the difficulty in obtaining chromosomal preparations, making them virtually unexplored. An example of this scenario relies on members of the family Lebiasinidae which include miniature to medium-sized, poorly known species, until very recently. The present study is part of undergoing major cytogenetic advances seeking to elucidate the evolutionary history of lebiasinids. Aiming to examine the karyotype diversification more deeply in Pyrrhulina, here we combined classical and molecular cytogenetic analyses, including Giemsa staining, C-banding, repetitive DNA mapping, comparative genomic hybridization (CGH), and whole chromosome painting (WCP) to perform the first analyses in five Pyrrhulina species (Pyrrhulina aff. marilynae, Pyrrhulina sp., P. obermulleri, P. marilynae and Pyrrhulina cf. laeta). The diploid number (2n) ranged from 40 to 42 chromosomes among all analyzed species, but P. marilynae is strikingly differentiated by having 2n = 32 chromosomes and a karyotype composed of large meta/submetacentric chromosomes, whose plesiomorphic status is discussed. The distribution of microsatellites does not markedly differ among species, but the number and position of the rDNA sites underwent significant changes among them. Interspecific comparative genome hybridization (CGH) found a moderate divergence in the repetitive DNA content among the species’ genomes. Noteworthy, the WCP reinforced our previous hypothesis on the origin of the X₁X₂Y multiple sex chromosome system in P. semifasciata. In summary, our data suggest that the karyotype differentiation in Pyrrhulina has been driven by major structural rearrangements, accompanied by high dynamics of repetitive DNAs.

Interrogating Phylogenetic Discordance Resolves Deep Splits in the Rapid Radiation of Old World Fruit Bats (Chiroptera: Pteropodidae)

The family Pteropodidae (Old World fruit bats) comprises $>$200 species distributed across the Old World tropics and subtropics. Most pteropodids feed on fruit, suggesting an early origin of frugivory, although several lineages have shifted to nectar-based diets. Pteropodids are of exceptional conservation concern with $>$50% of species considered threatened, yet the systematics of this group has long been debated, with uncertainty surrounding early splits attributed to an ancient rapid diversification. Resolving the relationships among the main pteropodid lineages is essential if we are to fully understand their evolutionary distinctiveness, and the extent to which these bats have transitioned to nectar-feeding. Here we generated orthologous sequences for $>$1400 nuclear protein-coding genes (2.8 million base pairs) across 114 species from 43 genera of Old World fruit bats (57% and 96% of extant species- and genus-level diversity, respectively), and combined phylogenomic inference with filtering by information content to resolve systematic relationships among the major lineages. Concatenation and coalescent-based methods recovered three distinct backbone topologies that were not able to be reconciled by filtering via phylogenetic information content. Concordance analysis and gene genealogy interrogation show that one topology is consistently the best supported, and that observed phylogenetic conflicts arise from both gene tree error and deep incomplete lineage sorting. In addition to resolving long-standing inconsistencies in the reported relationships among major lineages, we show that Old World fruit bats have likely undergone at least seven independent dietary transitions from frugivory to nectarivory. Finally, we use this phylogeny to identify and describe one new genus. [Chiroptera; coalescence; concordance; incomplete lineage sorting; nectar feeder; species tree; target enrichment.].

Comparative cytogenetics of Serrasalmidae (Teleostei: Characiformes): The relationship between chromosomal evolution and molecular phylogenies

Serrasalmidae has high morphological and chromosomal diversity. Based on molecular hypotheses, the family is currently divided into two subfamilies, Colossomatinae and Serrasalminae, with Serrasalminae composed of two tribes: Myleini (comprising most of pacus species) and Serrasalmini (represented by Metynnis, Catoprion, and remaining piranha’s genera). This study aimed to analyze species of the tribes Myleini (Myloplus asterias, M. lobatus, M. rubripinnis, M. schomburgki, and Tometes camunani) and Serrasalmini (Metynnis cuiaba, M. hypsauchen, and M. longipinnis) using classical and molecular cytogenetic techniques in order to understand the chromosomal evolution of the family. The four species of the genus Myloplus and T. camunani presented 2n = 58 chromosomes, while the species of Metynnis presented 2n = 62 chromosomes. The distribution of heterochromatin occurred predominantly in pericentromeric regions in all species. Tometes camunani and Myloplus spp. presented only one site with 5S rDNA. Multiple markers of 18S rDNA were observed in T. camunani, M. asterias, M. lobatus, M. rubripinnis, and M. schomburgkii. For Metynnis, however, synteny of the 18S and 5S rDNA was observed in the three species, in addition to an additional 5S marker in M. longipinnis. These data, when superimposed on the phylogeny of the family, suggest a tendency to increase the diploid chromosome number from 54 to 62 chromosomes, which occurred in a nonlinear manner and is the result of several chromosomal rearrangements. In addition, the different karyotype formulas and locations of ribosomal sequences can be used as cytotaxonomic markers and assist in the identification of species.

Ultra-Conserved Elements and morphology reciprocally illuminate conflicting phylogenetic hypotheses in Chalcididae (Hymenoptera, Chalcidoidea)

Recent technical advances combined with novel computational approaches have promised the acceleration of our understanding of the tree of life. However, when it comes to hyperdiverse and poorly known groups of invertebrates, studies are still scarce. As published phylogenies will be rarely challenged by future taxonomists, careful attention must be paid to potential analytical bias. We present the first molecular phylogenetic hypothesis for the family Chalcididae, a group of parasitoid wasps, with a representative sampling (144 ingroups and seven outgroups) that covers all described subfamilies and tribes, and 82% of the known genera. Analyses of 538 Ultra-Conserved Elements (UCEs) with supermatrix (RAxML and IQTREE) and gene tree reconciliation approaches (ASTRAL, ASTRID) resulted in highly supported topologies in overall agreement with morphology but reveal conflicting topologies for some of the deepest nodes. To resolve these conflicts, we explored the phylogenetic tree space with clustering and gene genealogy interrogation methods, analyzed marker and taxon properties that could bias inferences and performed a thorough morphological analysis (130 characters encoded for 40 taxa representative of the diversity). This joint analysis reveals that UCEs enable attainment of resolution between ancestry and convergent/divergent evolution when morphology is not informative enough, but also shows that a systematic exploration of bias with different analytical methods and a careful analysis of morphological features is required to prevent publication of artifactual results. We highlight a GC content bias for maximum-likelihood approaches, an artifactual mid-point rooting of the ASTRAL tree and a deleterious effect of high percentage of missing data (>85% missing UCEs) on gene tree reconciliation methods. Based on the results we propose a new classification of the family into eight subfamilies and ten tribes that lay the foundation for future studies on the evolutionary history of Chalcididae.

Comparing Ultraconserved Elements and Exons for Phylogenomic Analyses of Middle American Cichlids: When Data Agree to Disagree

Choosing among types of genomic markers to be used in a phylogenomic study can have a major influence on the cost, design, and results of a study. Yet few attempts have been made to compare categories of next-generation sequence markers limiting our ability to compare the suitability of these different genomic fragment types. Here, we explore properties of different genomic markers to find if they vary in the accuracy of component phylogenetic trees and to clarify the causes of conflict obtained from different data sets or inference methods. As a test case, we explore the causes of discordance between phylogenetic hypotheses obtained using a novel data set of ultraconserved elements (UCEs) and a recently published exon data set of the cichlid tribe Heroini. Resolving relationships among heroine cichlids has historically been difficult, and the processes of colonization and diversification in Middle America and the Greater Antilles are not yet well understood. Despite differences in informativeness and levels of gene tree discordance between UCEs and exons, the resulting phylogenomic hypotheses generally agree on most relationships. The independent data sets disagreed in areas with low phylogenetic signal that were overwhelmed by incomplete lineage sorting and nonphylogenetic signals. For UCEs, high levels of incomplete lineage sorting were found to be the major cause of gene tree discordance, whereas, for exons, nonphylogenetic signal is most likely caused by a reduced number of highly informative loci. This paucity of informative loci in exons might be due to heterogeneous substitution rates that are problematic to model (i.e., computationally restrictive) resulting in systematic errors that UCEs (being less informative individually but more uniform) are less prone to. These results generally demonstrate the robustness of phylogenomic methods to accommodate genomic markers with different biological and phylogenetic properties. However, we identify common and unique pitfalls of different categories of genomic fragments when inferring enigmatic phylogenetic relationships.

Chromosomal Analysis of Ctenolucius hujeta Valenciennes, 1850 (Characiformes): A New Piece in the Chromosomal Evolution of the Ctenoluciidae

Ctenoluciidae (Characiformes), a family of freshwater fishes, comprises 2 genera, Ctenolucius and Boulengerella, with 7 recognized species. Up to now, only species of the genus Boulengerella have been subjected to cytogenetic studies. Here, we investigated the karyotype and other cytogenetic features of pike characin, Ctenolucius hujeta, using conventional (Giemsa staining, C-banding, Ag-NOR staining) and molecular (rDNA, telomeric sequences, and fiber-FISH mapping) procedures. This species has a diploid chromosome number of 2n = 36, and a karyotype composed of 12m + 20sm + 4a and FN = 68, similar to that found in Boulengerella species. However, differences regarding the number and distribution of several chromosomal markers support a distinct generic status. Colocalization of the 18S and 5S rDNA genes is an exclusive characteristic of the C. hujeta genome, with an interspersed distribution in the chromosomal fiber, an unusual phenomenon among eukaryotes. Additionally, our results support the view that Ctenoluciidae and Lebiasinidae families are closely related.

Total evidence phylogenetic analysis reveals polyphyly of Anostomoides and uncovers an unexpectedly ancient genus of Anostomidae fishes (Characiformes)

The nearly 150 species of Anostomidae comprise one of the most diverse and taxonomically dynamic families of Neotropical freshwater fishes. A recent revision of the enigmatic and poorly diagnosed genus Anostomoides demonstrated that it contains two valid species, each with complicated taxonomic histories; however, that study did not address their phylogenetic placement. Herein, we integrate molecular and morphological data to demonstrate their distant evolutionary relationship, and thus the polyphyly of Anostomoides. While we reconstruct one of the species in a previously hypothesized placement within a clade also containing Laemolyta, Rhytiodus and Schizodon, the other represents a morphologically and genetically distinctive lineage that diverged early in the history of the family. We describe and illustrate the osteology of this remarkable species, discuss the evolutionary implications of its unique suite of features, and use those characteristics to diagnose a new genus that evolved independently of all other known members of the family for approximately 37 Myr.

Accelerated Diversification Explains the Exceptional Species Richness of Tropical Characoid Fishes

The Neotropics harbor the most species-rich freshwater fish fauna on the planet, but the timing of that exceptional diversification remains unclear. Did the Neotropics accumulate species steadily throughout their long history, or attain their remarkable diversity recently? Biologists have long debated the relative support for these museum and cradle hypotheses, but few phylogenies of megadiverse tropical clades have included sufficient taxa to distinguish between them. We used 1,288 ultraconserved element loci (UCE) spanning 293 species, 211 genera and 21 families of characoid fishes to reconstruct a new, fossil-calibrated phylogeny and infer the most likely diversification scenario for a clade that includes a third of Neotropical fish diversity. This phylogeny implies paraphyly of the traditional delimitation of Characiformes because it resolves the largely Neotropical Characoidei as the sister lineage of Siluriformes (catfishes), rather than the African Citharinodei. Time-calibrated phylogenies indicate an ancient origin of major characoid lineages and reveal a much more recent emergence of most characoid species. Diversification rate analyses infer increased speciation and decreased extinction rates during the Oligocene at around 30 million years ago (Ma) during a period of mega-wetland formation in the proto-Orinoco-Amazonas. Three species-rich and ecomorphologically diverse lineages (Anostomidae, Serrasalmidae, and Characidae) that originated more than 60 Ma in the Paleocene experienced particularly notable bursts of Oligocene diversification and now account collectively for 68% of the approximately 2,150 species of Characoidei. In addition to paleogeographic changes, we discuss potential accelerants of diversification in these three lineages. While the Neotropics accumulated a museum of ecomorphologically diverse characoid lineages long ago, this geologically dynamic region also cradled a much more recent birth of remarkable species-level diversity.

Development and phylogenetic utilities of a new set of single-/low-copy nuclear genes in Senecioneae (Asteraceae), with new insights into the tribal position and the relationships within subtribe Tussilagininae

The tribe Senecioneae is one of the largest tribes in Asteraceae, with a nearly cosmopolitan distribution. Despite great efforts devoted to elucidate the evolution of Senecioneae, many questions still remain concerning the systematics of this group, from the tribal circumscription and position to species relationships in many genera. The hybridization-based target enrichment method of next-generation sequencing has been accepted as a promising approach to resolve phylogenetic problems. We herein develop a set of single-/low-copy genes for Senecioneae, and test their phylogenetic utilities. Our results demonstrate that these genes work highly efficiently for Senecioneae, with a high average gene recovery of 98.8% across the tribe and recovering robust phylogenetic hypotheses at different levels. In particular, the delimitation of the Senecioneae has been confirmed to include Abrotanella and exclude Doronicum, with the former sister to core Senecioneae and the latter shown to be more closely related to Calenduleae. Moreover, Doronicum and Calenduleae are inferred to be the closest relatives of Senecioneae, which is a new hypothesis well supported by statistical topology tests, morphological evidence, and the profile of pyrrolizidine alkaloids, a special kind of chemical characters generally used to define Senecioneae. Furthermore, this study suggests a complex reticulation history in the diversification of Senecioneae, accounting for the prevalence of polyploid groups in the tribe. With subtribe Tussilagininae s.str. as a case study showing a more evident pattern of gene duplication, we further explored reconstructing the phylogeny in the groups with high ploidy levels. Our results also demonstrate that tree topologies based on sorted paralogous copies are stable across different methods of phylogenetic inference, and more congruent with the morphological evidence and the results of previous phylogenetic studies.

New insights into patterns and rates of tooth replacement in serrasalmid and characid fishes, with implications for the subsistence fishery of Peru's remote ribereños villages

We explored patterns, rates and unexpected socio-ecological consequences of tooth replacement in serrasalmids and characids of the Peruvian Amazon using microcomputed tomography. Of 24 specimens collected in February 2019, representing a mix of red-bellied piranha Pygocentrus nattereri, redeye piranha Serrasalmus rhombeus, silver dollar fish Ctenobrycon hauxwellianus and mojara Astyanax abramis, six individuals possessed edentulous jaw quadrants. On average, 22.9% of fish collected per day from these species featured incomplete dentition, a value three to five times higher than anticipated based on replacement rates estimated from captive fish, differences that may be driven by ontogeny, seasonality or environmental quality.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Testing the utility of alternative metrics of branch support to address the ancient evolutionary radiation of tunas, stromateoids, and allies (Teleostei: Pelagiaria)

The use of high-throughput sequencing technologies to produce genome-scale datasets was expected to settle some long-standing controversies across the Tree of Life, particularly in areas where short branches occur at deep timescales. Instead, these datasets have often yielded many well-supported but conflicting topologies, and highly variable gene-tree distributions. A variety of branch-support metrics beyond the nonparametric bootstrap are now available to assess how robust a phylogenetic hypothesis may be, as well as new methods to quantify gene-tree discordance. We applied multiple branch support metrics to an ancient group of marine fishes (Teleostei: Pelagiaria) whose interfamilial relationships have proven difficult to resolve due to a rapid accumulation of lineages very early in its history. We analyzed hundreds of loci including published UCE data and newly generated exonic data along with their flanking regions to represent all 16 extant families for more than 150 out of 284 valid species in the group. Branch support was lower for interfamilial relationships (except the SH-like aLRT and aBayes methods) regardless of the type of marker used. Several nodes that were highly supported with bootstrap had very low site and gene-tree concordance, revealing underlying conflict. Despite this conflict, we were able to identify four consistent interfamilial clades, each comprised of two or three families. Combining exons with their flanking regions also produced increased branch lengths in the deep branches of the pelagiarian tree. Our results demonstrate the limitations of employing current metrics of branch support and species-tree estimation when assessing the confidence of ancient evolutionary radiations and emphasize the necessity to embrace alternative measurements to explore phylogenetic uncertainty and discordance in phylogenomic datasets.

Partitioned Gene-Tree Analyses and Gene-Based Topology Testing Help Resolve Incongruence in a Phylogenomic Study of Host-Specialist Bees (Apidae: Eucerinae)

Incongruence among phylogenetic results has become a common occurrence in analyses of genome-scale data sets. Incongruence originates from uncertainty in underlying evolutionary processes (e.g., incomplete lineage sorting) and from difficulties in determining the best analytical approaches for each situation. To overcome these difficulties, more studies are needed that identify incongruences and demonstrate practical ways to confidently resolve them. Here, we present results of a phylogenomic study based on the analysis 197 taxa and 2,526 ultraconserved element (UCE) loci. We investigate evolutionary relationships of Eucerinae, a diverse subfamily of apid bees (relatives of honey bees and bumble bees) with >1,200 species. We sampled representatives of all tribes within the group and >80% of genera, including two mysterious South American genera, Chilimalopsis and Teratognatha. Initial analysis of the UCE data revealed two conflicting hypotheses for relationships among tribes. To resolve the incongruence, we tested concatenation and species tree approaches and used a variety of additional strategies including locus filtering, partitioned gene-trees searches, and gene-based topological tests. We show that within-locus partitioning improves gene tree and subsequent species-tree estimation, and that this approach, confidently resolves the incongruence observed in our data set. After exploring our proposed analytical strategy on eucerine bees, we validated its efficacy to resolve hard phylogenetic problems by implementing it on a published UCE data set of Adephaga (Insecta: Coleoptera). Our results provide a robust phylogenetic hypothesis for Eucerinae and demonstrate a practical strategy for resolving incongruence in other phylogenomic data sets.

The wall lizards of the Balkan peninsula: Tackling questions at the interface of phylogenomics and population genomics

Wall lizards of the genus Podarcis (Sauria, Lacertidae) are the predominant reptile group in southern Europe, including 24 recognized species. Mitochondrial DNA data have shown that, with the exception of P. muralis, the Podarcis species distributed in the Balkan peninsula form a species group that is further sub-divided into two subgroups: the one of "P. tauricus" consisting of P. tauricus, P. milensis, P. gaigeae, and P. melisellensis, and the other of "P. erhardii" comprising P. erhardii, P. levendis, P. cretensis, and P. peloponnesiacus. In an attempt to explore the Balkan Podarcis phylogenomic relationships, assess the levels of genetic structure and to re-evaluate the number of extant species, we employed phylogenomic and admixture approaches on ddRADseq (double digested Restriction site Associated DNA sequencing) genomic data. With this efficient Next Generation Sequencing approach, we were able to obtain a large number of genomic loci randomly distributed throughout the genome and use them to resolve the previously obscure phylogenetic relationships among the different Podarcis species distributed in the Balkans. The obtained phylogenomic relationships support the monophyly of both aforementioned subgroups and revealed several divergent lineages within each subgroup, stressing the need for taxonomic re-evaluation of Podarcis' species in Balkans. The phylogenomic trees and the species delimitation analyses confirmed all recently recognized species (P. levendis, P. cretensis, and P. ionicus) and showed the presence of at least two more species, one in P. erhardii and the other in P. peloponnesiacus.

Differences in developmental potential predict the contrasting patterns of dental diversification in characiform and cypriniform fishes

Morphological diversification during adaptive radiation may depend on factors external or internal to the lineage. We provide evidence for the latter in characiform fishes (tetras and piranhas), which exhibit extensive dental diversity. Phylogenetic character mapping supported regain of lost teeth as contributing to this diversity. To test for latent potential for dentition that would facilitate its evolutionary expansion, we overexpressed a tooth initiation signal, the tumour necrosis factor pathway ligand ectodysplasin, in a model characiform, the Mexican tetra (Astyanax mexicanus). This manipulation resulted in extensive ectopic dentition, in contrast with its previously reported limited effect in the zebrafish (Danio rerio). Tooth location in the order Cypriniformes, to which the zebrafish belongs, is much more restricted than in characiforms, a pattern that may be explained by differences in the retention of ancestral developmental potential. Our results suggest that differences in evolvability between lineages may lead to contrasting patterns of diversification.

Adaptation to herbivory and detritivory drives the convergent evolution of large abdominal cavities in a diverse freshwater fish radiation (Otophysi: Characiformes)

Convergent evolution is often interpreted as evidence of natural selection favoring an optimal phenotype during adaptation. Morphological convergence is frequently found among lineages that converge on diet, but most studies have focused on morphological traits that relate exclusively to food handling and processing. In vertebrates, there is a strong inverse relationship between intestine length and trophic level. However, little is known about whether adaptation to a low trophic level influences the evolution of abdominal cavities that can accommodate larger intestines. Here, I reconstruct the evolutionary history of trophic ecology and examine abdominal cavity shape across 157 species of the fish order Characiformes to determine whether adaptation to an herbivorous-detritivorous diet drives convergent evolution of large abdominal cavities. Herbivorous-detritivorous species evolved significantly larger abdominal cavities than other trophic groups and repeatedly converged on a similar abdominal cavity morphology. Other trophic groups evolved abdominal cavity morphologies either stochastically or by selective pressures from an untested ecological character. These findings demonstrate that the selective demands of a larger intestinal tract promote the repeated convergence of a large abdominal cavity within herbivorous-detritivorous characiform fishes, while allowing other lineages to evolve randomly or adapt in response to other selection pressures, contributing to the overall body shape diversity of the order.

A Long-Term Conserved Satellite DNA That Remains Unexpanded in Several Genomes of Characiformes Fish Is Actively Transcribed

Eukaryotic genomes contain large amounts of repetitive DNA sequences, such as tandemly repeated satellite DNAs (satDNAs). These sequences are highly dynamic and tend to be genus- or species-specific due to their particular evolutionary pathways, although there are few unusual cases of conserved satDNAs over long periods of time. Here, we used multiple approaches to reveal that an satDNA named CharSat01-52 originated in the last common ancestor of Characoidei fish, a superfamily within the Characiformes order, ∼140-78 Ma, whereas its nucleotide composition has remained considerably conserved in several taxa. We show that 14 distantly related species within Characoidei share the presence of this satDNA, which is highly amplified and clustered in subtelomeric regions in a single species (Characidium gomesi), while remained organized as small clusters in all the other species. Defying predictions of the molecular drive of satellite evolution, CharSat01-52 shows similar values of intra- and interspecific divergence. Although we did not provide evidence for a specific functional role of CharSat01-52, its transcriptional activity was demonstrated in different species. In addition, we identified short tandem arrays of CharSat01-52 embedded within single-molecule real-time long reads of Astyanax paranae (536 bp-3.1 kb) and A. mexicanus (501 bp-3.9 kb). Such arrays consisted of head-to-tail repeats and could be found interspersed with other sequences, inverted sequences, or neighbored by other satellites. Our results provide a detailed characterization of an old and conserved satDNA, challenging general predictions of satDNA evolution.

Mitochondrial Genome Structures and Phylogenetic Analyses of Two Tropical Characidae Fishes

The Characidae family contains the largest number of tropical fish species. Morphological similarities make species identification difficult within this family. Here, the complete mitogenomes of two Characidae fish were determined and comparatively analyzed with those of nine other Characidae fish species. The two newly sequenced complete mitogenomes are circular DNA molecules with sizes of 16,701 bp (Hyphessobrycon amandae; MT484069) and 16,710 bp (Hemigrammus erythrozonus; MT484070); both have a highly conserved structure typical of Characidae, with the start codon ATN (ATG/ATT) and stop codon TAR (TAA/TAG) or an incomplete T--/TA-. Most protein-coding genes of the 11 Characidae mitogenomes showed significant codon usage bias, and the protein-coding gene cox1 was found to be a comparatively slow-evolving gene. Phylogenetic analyses via the maximum likelihood and Bayesian inference methods confirmed that H. amandae and H. erythrozonus belong to the family Characidae. In all Characidae species studied, one genus was well supported; whereas other two genera showed marked differentiation. These findings provide a phylogenetic basis for improved classification of the family Characidae. Determining the mitogenomes of H. erythrozonus and H. amandae improves our understanding of the phylogeny and evolution of fish species.

Cryptic speciation in populations of the genus Aphyocharax (Characiformes: Characidae) from eastern Amazon coastal river drainages and surroundings revealed by single locus species delimitation methods

Abstract Recent studies in eastern Amazon coastal drainages and their surroundings have revealed new fish species that sometimes exhibit little morphological differentiation (cryptic species). Thus, we used a DNA-based species delimitation approach to test if populations showing the morphotype and typical character states of the Aphyocharax avary holotype correspond either to A. avary or A. brevicaudatus, two known species from the region, or if they form independent lineages, indicating cryptic speciation. WP and GMYC analyses recovered five lineages (species) in the ingroup, while a bPTP analysis delimited three lineages. ABGD analyses produced two possible results: one corroborating the WP and GMYC methods and another corroborating the bPTP method. All methods indicate undescribed cryptic species in the region and show variation from at least 1 to 4 species in the ingroup, depending on the approach, corroborating previous studies, and revealing this region as a possible hotspot for discovering undescribed fish species.

Chromosome analysis in Saccodon wagneri (Characiformes) and insights into the karyotype evolution of Parodontidae

ABSTRACT Parodontidae is a relatively small group of Neotropical characiform fishes consisting of three genera (Apareiodon, Parodon, and Saccodon) with 32 valid species. A vast cytogenetic literature is available on Apareiodon and Parodon, but to date, there is no cytogenetic data about Saccodon, a genus that contains only three species with a trans-Andean distribution. In the present study the karyotype of S. wagneri was described, based on both conventional (Giemsa staining, Ag-NOR, C-bands) and molecular (repetitive DNA mapping by fluorescent in situ hybridization) methods. A diploid chromosome number of 2n = 54 was observed in both sexes, and the presence of heteromorphic sex chromosomes of the ZZ/ZW type was detected. The W chromosome has a terminal heterochromatin band that occupies approximately half of the long arm, being this band approximately half the size of the Z chromosome. The FISH assay showed a synteny of the 18S-rDNA and 5S-rDNA genes in the chromosome pair 14, and the absence of interstitial telomeric sites. Our data reinforce the hypothesis of a conservative karyotype structure in Parodontidae and suggest an ancient origin of the sex chromosomes in the fishes of this family.

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.

Revisiting the phylogeny of the genus Lolliguncula Steenstrup 1881 improves understanding of their biogeography and proves the validity of Lolliguncula argus Brakoniecki & Roper, 1985

The biogeography of American loliginid squids has been improved in recent years, but certain key taxa have been missing. Given that the most accurate phylogenies and estimates of divergence times of common ancestors depend heavily on good taxonomic coverage we have reanalyzed the genus Lolliguncula in light of new samples that increase the geographic and taxonomic coverage. New sequences were produced using standard methods to update an existing dataset for COI, 16S and Rhodopsin markers. Data was analyzed using various species delimitation methods, rigorous phylogenetic analyses and estimates of divergence times between clades. Within Lolliguncula we recover five monophyletic lineages that relate to the known species L. argus, L. diomedeae, L. panamensis, L. brevis North Atlantic and L. brevis South Atlantic. Except when using low divergence thresholds in ABGD, species delimitation methods only identify four of these lineages as distinct species, grouping L. argus and L. diomedeae as a single species. However, considering the reciprocal monophyly, recent divergence time estimate and morphological diagnoses we refrain from synonymizing L. argus within L. diomedeae, considering them very recently diverged species. The biogeography of the American loliginids is discussed, wherein basal cladogenesis in both Lolliguncula and Doryteuthis occur between the Atlantic and Pacific about 45 mya, with subsequent speciation around 20 mya associated with seafloor changes during the formation of the Caribbean. The recent speciation between L. argus and L. diomedeae is associated to oceanic environmental changes associated with glaciation, deep sea cooling and tropical upwelling.

Phylogenomics of Piranhas and Pacus (Serrasalmidae) Uncovers How Dietary Convergence and Parallelism Obfuscate Traditional Morphological Taxonomy

The Amazon and neighboring South American river basins harbor the world's most diverse assemblages of freshwater fishes. One of the most prominent South American fish families is the Serrasalmidae (pacus and piranhas), found in nearly every continental basin. Serrasalmids are keystone ecological taxa, being some of the top riverine predators as well as the primary seed dispersers in the flooded forest. Despite their widespread occurrence and notable ecologies, serrasalmid evolutionary history and systematics are controversial. For example, the sister taxon to serrasalmids is contentious, the relationships of major clades within the family are inconsistent across different methodologies, and half of the extant serrasalmid genera are suggested to be non-monophyletic. We analyzed exon capture to reexamine the evolutionary relationships among 63 (of 99) species across all 16 serrasalmid genera and their nearest outgroups, including multiple individuals per species to account for cryptic lineages. To reconstruct the timeline of serrasalmid diversification, we time-calibrated this phylogeny using two different fossil-calibration schemes to account for uncertainty in taxonomy with respect to fossil teeth. Finally, we analyzed diet evolution across the family and comment on associated changes in dentition, highlighting the ecomorphological diversity within serrasalmids. We document widespread non-monophyly of genera within Myleinae, as well as between Serrasalmus and Pristobrycon, and propose that reliance on traits like teeth to distinguish among genera is confounded by ecological homoplasy, especially among herbivorous and omnivorous taxa. We clarify the relationships among all serrasalmid genera, propose new subfamily affiliations, and support hemiodontids as the sister taxon to Serrasalmidae. [Characiformes; exon capture; ichthyochory; molecular time-calibration; piscivory.].

An Insight into the Chromosomal Evolution of Lebiasinidae (Teleostei, Characiformes)

Lebiasinidae fishes have been historically neglected by cytogenetical studies. Here we present a genomic comparison in eleven Lebiasinidae species, in addition to a review of the ribosomal DNA sequences distribution in this family. With that, we develop ten sets of experiments in order to hybridize the genomic DNA of representative species from the genus Copeina, Copella, Nannostomus, and Pyrrhulina in metaphase plates of Lebiasina melanoguttata. Two major pathways on the chromosomal evolution of these species can be recognized: (i) conservation of 2n = 36 bi-armed chromosomes in Lebiasininae, as a basal condition, and (ii) high numeric and structural chromosomal rearrangements in Pyrrhulininae, with a notable tendency towards acrocentrization. The ribosomal DNA (rDNA) distribution also revealed a marked differentiation during the chromosomal evolution of Lebiasinidae, since both single and multiple sites, in addition to a wide range of chromosomal locations can be found. With some few exceptions, the terminal position of 18S rDNA appears as a common feature in Lebiasinidae-analyzed species. Altogether with Ctenoluciidae, this pattern can be considered a symplesiomorphism for both families. In addition to the specific repetitive DNA content that characterizes the genome of each particular species, Lebiasina also keeps inter-specific repetitive sequences, thus reinforcing its proposed basal condition in Lebiasinidae.

First Record on Sex Chromosomes in a Species of the Family Cynodontidae: Cynodon gibbus (Agassiz, 1829)

The fish family Cynodontidae belongs to the superfamily Curimatoidea, together with the Hemiodontidae, Serrasalmidae, Parodontidae, Prochilodontidae, Chilodontidae, Curimatidae, and Anostomidae. The majority of the species of this superfamily that have been analyzed to date have a diploid chromosome number of 2n = 54. Differentiated sex chromosomes (with female heterogamety) have been observed only in the Prochilodontidae, Parodontidae, and Anostomidae. The present study provides the first description of differentiated sex chromosomes in the cynodontid species Cynodon gibbus, which has a ZZ/ZW system, and shows that repetitive DNA has played a fundamental role in the differentiation of these sex chromosomes.

Centric Fusions behind the Karyotype Evolution of Neotropical Nannostomus Pencilfishes (Characiforme, Lebiasinidae): First Insights from a Molecular Cytogenetic Perspective

Lebiasinidae is a Neotropical freshwater family widely distributed throughout South and Central America. Due to their often very small body size, Lebiasinidae species are cytogenetically challenging and hence largely underexplored. However, the available but limited karyotype data already suggested a high interspecific variability in the diploid chromosome number (2n), which is pronounced in the speciose genus Nannostomus, a popular taxon in ornamental fish trade due to its remarkable body coloration. Aiming to more deeply examine the karyotype diversification in Nannostomus, we combined conventional cytogenetics (Giemsa-staining and C-banding) with the chromosomal mapping of tandemly repeated 5S and 18S rDNA clusters and with interspecific comparative genomic hybridization (CGH) to investigate genomes of four representative Nannostomus species: N. beckfordi, N. eques, N. marginatus, and N. unifasciatus. Our data showed a remarkable variability in 2n, ranging from 2n = 22 in N. unifasciatus (karyotype composed exclusively of metacentrics/submetacentrics) to 2n = 44 in N. beckfordi (karyotype composed entirely of acrocentrics). On the other hand, patterns of 18S and 5S rDNA distribution in the analyzed karyotypes remained rather conservative, with only two 18S and two to four 5S rDNA sites. In view of the mostly unchanged number of chromosome arms (FN = 44) in all but one species (N. eques; FN = 36), and with respect to the current phylogenetic hypothesis, we propose Robertsonian translocations to be a significant contributor to the karyotype differentiation in (at least herein studied) Nannostomus species. Interspecific comparative genome hybridization (CGH) using whole genomic DNAs mapped against the chromosome background of N. beckfordi found a moderate divergence in the repetitive DNA content among the species’ genomes. Collectively, our data suggest that the karyotype differentiation in Nannostomus has been largely driven by major structural rearrangements, accompanied by only low to moderate dynamics of repetitive DNA at the sub-chromosomal level. Possible mechanisms and factors behind the elevated tolerance to such a rate of karyotype change in Nannostomus are discussed.

Checklist of the fishes from Jamari River basin, in areas under influence of dams, Rondônia, Brazil

Abstract: The Madeira River is the most extensive tributary of the Amazon River and has the largest diversity of fishes in the world. On its right bank, the Madeira River receives the Jamari River, in which the first hydroelectric power plant (HPP) in State of Rondônia, Samuel HPP, was built. Besides this, other dams were built in the Jamari River and its tributaries, however, the available information in the scientific literature about the ichthyofaunistic diversity of this basin is rare. This work aims to provide an ichthyofaunistic inventory in a region of the Jamari river basin, in the State of Rondônia, where three small hydropower plants (SHPs) were implemented. The ichthyofauna was sampled in 16 expeditions between August 2015 and December 2018. Gill nets and seine nets were used with different meshes, as well as longlines and cast nets at different times of the day. Additionally, 81 INPA lots of species from the Samuel HPP area of influence were reanalyzed. Fish were identified according to the specialized literature, as well as in consultations with experts of various taxonomic groups. Voucher specimens of the species were cataloged and deposited in the ichthyological collection of the Núcleo de Pesquisas em Limnologia, Ictiologia e Aquicultura (Nupélia) of the Universidade Estadual de Maringá. A total of 230 species were recorded, of which 22 were putative new species, 117 were added to the Jamari River basin and 28 to the Madeira River basin. The continuation of the studies in this section of the Jamari river basin is fundamental for analysis of local impact due to the presence of dams. Moreover, the addition of putative new species to the Madeira River basin indicates gaps in the knowledge of Neotropical ichthyofauna.