The tree of life of fishes is in a state of flux because we still lack a comprehensive phylogeny that includes all major groups. The situation is most critical for a large clade of spiny-finned fishes, traditionally referred to as percomorphs, whose uncertain relationships have plagued ichthyologists for over a century. Most of what we know about the higher-level relationships among fish lineages has been based on morphology, but rapid influx of molecular studies is changing many established systematic concepts. We report a comprehensive molecular phylogeny for bony fishes that includes representatives of all major lineages. DNA sequence data for 21 molecular markers (one mitochondrial and 20 nuclear genes) were collected for 1410 bony fish taxa, plus four tetrapod species and two chondrichthyan outgroups (total 1416 terminals). Bony fish diversity is represented by 1093 genera, 369 families, and all traditionally recognized orders. The maximum likelihood tree provides unprecedented resolution and high bootstrap support for most backbone nodes, defining for the first time a global phylogeny of fishes. The general structure of the tree is in agreement with expectations from previous morphological and molecular studies, but significant new clades arise. Most interestingly, the high degree of uncertainty among percomorphs is now resolved into nine well-supported supraordinal groups. The order Perciformes, considered by many a polyphyletic taxonomic waste basket, is defined for the first time as a monophyletic group in the global phylogeny. A new classification that reflects our phylogenetic hypothesis is proposed to facilitate communication about the newly found structure of the tree of life of fishes. Finally, the molecular phylogeny is calibrated using 60 fossil constraints to produce a comprehensive time tree. The new time-calibrated phylogeny will provide the basis for and stimulate new comparative studies to better understand the evolution of the amazing diversity of fishes.

Sequencing reduced-representation libraries of restriction site-associated DNA (RADseq) to identify single nucleotide polymorphisms (SNPs) is quickly becoming a standard methodology for molecular ecologists. Because of the scale of RADseq data sets, putative loci cannot be assessed individually, making the process of filtering noise and correctly identifying biologically meaningful signal more difficult. Artefacts introduced during library preparation and/or bioinformatic processing of SNP data can create patterns that are incorrectly interpreted as indicative of population structure or natural selection. Therefore, it is crucial to carefully consider types of errors that may be introduced during laboratory work and data processing, and how to minimize, detect and remove these errors. Here, we discuss issues inherent to RADseq methodologies that can result in artefacts during library preparation and locus reconstruction resulting in erroneous SNP calls and, ultimately, genotyping error. Further, we describe steps that can be implemented to create a rigorously filtered data set consisting of markers accurately representing independent loci and compare the effect of different combinations of filters on four RAD data sets. At last, we stress the importance of publishing raw sequence data along with final filtered data sets in addition to detailed documentation of filtering steps and quality control measures.

High species richness and evenness in structurally complex habitats has been hypothesized to be associated with niche partitioning. To test this idea, relationships between habitat structural complexity in river littoral-zone habitats and morphological diversity of tropical fishes were examined in the Cinaruco River, Venezuela. Six habitat attributes were quantified in 45 sites spanning a range of structural complexity. Fishes were collected during day and night to estimate species density and relative abundances at each site. Twenty-two morphological variables were measured for each species. Principal components analysis (PCA) of physical habitat data yielded two axes that modeled >80% of variation across sites. The first two axes from PCA of fish morphological variables modeled >70% of variation. Species density during both day and night was negatively associated with flow velocity and positively associated with habitat complexity. Similarity of day and night samples from the same site was significantly greater for sites with high habitat complexity and low flow. In general, mean local assemblage morphological PC scores were not significantly associated with habitat PC scores. Average, maximum, and standard deviation of morphological Euclidean distances of local assemblages revealed positive associations with structural complexity and negative associations with flow. These relationships held even when the positive relationship of species density was statistically removed from assemblage morphological patterns. Findings suggest that both species niche compression and assemblage niche space increase when habitat complexity is greater and flow velocity is lower in this tropical lowland river.