Multi-locus fossil-calibrated phylogeny of Atheriniformes (Teleostei, Ovalentaria)

Differential expression analyses and detection of SNP loci associated with environmental variables: Are salinity and temperature factors involved in population differentiation and speciation in Odontesthes?

Environmental factors play a key role in individual adaptation to different local conditions. Because of this, studies about the physiological and genetic responses of individuals exposed to different natural environments offer clues about mechanisms involved in population differentiation, and as a subsequent result, speciation. Marine environments are especially suited to survey this kind of phenomena because they commonly harbor species adapted to different local conditions along a geographic continuum. Silversides belonging to Odontesthes are commonly distributed in tropical and temperate regions of South America and exhibit noticeable phenotypic plasticity, which allows them to adapt to contrasting environments. In this study, the genetic expression of O. argentinensis sampled along the Uruguayan Atlantic coast and estuarine adjacent areas was investigated. In addition, the correlation between individual genotypes and environmental variables was also analysed in O. argentinensis and O. bonariensis. Results obtained suggest a differential expression pattern of low magnitude among individuals from the different areas sampled and a correlation between several SNP loci and environmental variables. The analyses carried out did not show a clear differentiation among individuals sampled along different salinity regimens, but enriched GOTerms seem to be driven by water oxygen content. On the other hand, a total of 46 SNPs analysed in O. argentinensis and O. bonariensis showed a correlation with salinity and temperature. Although none of the correlated SNPs and corresponding genes from our both analyses were directly associated with hypoxia, genes related to the cardiovascular system and muscle cell differentiation were found. All these genes are interesting candidates for future studies since they are closely related to the differentially expressed genes. Although salinity was also mentioned as an important parameter limiting introgression between O. argentinensis and O. bonariensis, it was found that salinity does not drive differential expression in O. argentinensis, but rather oxygen levels. Moreover, salinity does not directly affect the structure and genetic divergence of the populations, they appear to be structured based on their degree of isolation and geographical distance between them. Further studies, like genome-wide analyses, could help to elucidate additional genes adapted to the different environments in these silverside species.

Cenozoic colonisation of the Indian Ocean region by the Australian freshwater originating glassperch family Ambassidae (Teleostei)

We examined the phylogeny and biogeography of the glassperch family Ambassidae (Teleostei), which is widely distributed in the freshwater, brackish and marine coastal habitats across the Indo-West Pacific region. We first built a comprehensive time-calibrated phylogeny of Ambassidae using five genes. We then used this tree to reconstruct the evolution of the salinity preference and ancestral areas. Our results indicate that the two largest genera of Ambassidae, Ambassis and Parambassis, are each not monophyletic. The most recent common ancestor of Ambassidae was freshwater adapted and lived in Australia about 56 million years ago. Three independent freshwater-to-marine transitions are inferred, but no marine-to-freshwater ones. To explain the distribution of ambassids, we hypothesise two long-distance marine dispersal events from Australia. A first event was towards Southeast Asia during the early Cenozoic, followed by a second one towards Africa during mid-Cenozoic. The phylogenetic signal associated with the salinity adaptation of these events was not detected, possibly because of the selective extinction of intermediate marine lineages. The Ambassidae shares two characteristics with other freshwater fish groups distributed in continental regions surrounding the Indian Ocean: They are too young to support the hypothesis that their distribution is the result of the fragmentation of Gondwana, but they did not retain the phylogenetic signal of their marine dispersal.

Modulation of miR-429 during osmotic stress in the silverside Odontesthes humensis

Silverside fish inhabit marine coastal waters, coastal lagoons, and estuarine regions in southern South America. Although silversides are not fully adapted to freshwater, they can tolerate a wide range of salinity variations. MicroRNAs (miRNAs) are a class of ∼22 nucleotide noncoding RNAs, which are crucial regulators of gene expression at post-transcriptional level. Current data indicate that miRNAs biogenesis is altered by situations of environmental stress, thereby altering the expression of target mRNAs. Foremost, the silversides were acutely exposed to 30 g.L−1 of salt to reveal in which tissue miR-429 could be differentially expressed. Thus, fish were acclimated to freshwater (0 g.L−1) and to brackish water (10 g.L−1), and then exposed to opposite salinity treatment. Here, we reveal that miR-429, a gill-enriched miRNA, emerges as a prime osmoregulator in silversides. Taken together, our findings suggest that miR-429 is an endogenous regulator of osmotic stress, which may be developed as a biomarker to assist silverside aquaculture.

Patterns of Phenotypic Evolution Associated with Marine/Freshwater Transitions in Fishes

Evolutionary transitions between marine and freshwater ecosystems have occurred repeatedly throughout the phylogenetic history of fishes. The theory of ecological opportunity predicts that lineages that colonize species-poor regions will have greater potential for phenotypic diversification than lineages invading species-rich regions. Thus, transitions between marine and freshwaters may promote phenotypic diversification in trans-marine/freshwater fish clades. We used phylogenetic comparative methods to analyze body size data in nine major fish clades that have crossed the marine/freshwater boundary. We explored how habitat transitions, ecological opportunity, and community interactions influenced patterns of phenotypic diversity. Our analyses indicated that transitions between marine and freshwater habitats did not drive body size evolution, and there are few differences in body size between marine and freshwater lineages. We found that body size disparity in freshwater lineages is not correlated with the number of independent transitions to freshwaters. We found a positive correlation between body size disparity and overall species richness of a given area, and a negative correlation between body size disparity and diversity of closely related species. Our results indicate that the diversity of incumbent freshwater species does not restrict phenotypic diversification, but the diversity of closely related taxa can limit body size diversification. Ecological opportunity arising from colonization of novel habitats does not seem to have a major effect in the trajectory of body size evolution in trans-marine/freshwater clades. Moreover, competition with closely related taxa in freshwaters has a greater effect than competition with distantly related incumbent species.

Morphological characters in light of new molecular phylogenies: the caudal-fin skeleton of Ovalentaria

The Ovalentaria is a taxon of teleosts that has been proposed based on molecular analyses only. Previously widely separated families are assembled in this taxon. For the first time, the Ovalentaria are analysed using a comparative morphological approach. The caudal-fin skeleton of 355 species covering all 48 ovalentarian families are examined in cleared and stained specimens, µCT datasets and X-ray images as well as from the literature. A total of 38 morphological characters are evaluated and used for ancestral character state reconstructions and phylogenetic analyses. Results provide hypotheses for a scenario of the evolution of the caudal-fin skeleton and its ground plan in Ovalentaria. An evolutionary trend towards the reduction of skeletal elements in the caudal fin is observed. Connections between the evolution of the caudal-fin skeleton and modes of locomotion found in ovalentarian taxa are discussed. Phylogenetic analyses based on the caudal-fin morphology provide topologies for intra-ovalentarian relationships that largely agree with molecular hypotheses.

Where the Ends Meet: An Overview of Sex Determination in Atheriniform Fishes

Atheriniform fishes have recently emerged as attractive models for evolutionary, ecological, and molecular/physiological studies on sex determination. Many species in this group have marked temperature-dependent sex determination (TSD) and yet many species also have a sex determinant gene that provides a strong drive for male differentiation. Thus, in these species the 2 forms of sex determination that were once considered to be mutually exclusive, environmental (ESD) and genotypic (GSD) sex determination, can coexist at environmentally relevant conditions. Here, we review the current knowledge on sex determination in atheriniform fishes with emphasis on the molecular and physiological mechanisms of ESD and GSD, the coexistence and cross-talk between these 2 mechanisms, the possibility of extragonadal transduction of environmental information and/or extragonadal onset of sex determination, and the results of field studies applying novel tools such as otolith increment analysis and molecular markers of genetic sex developed for selected New World and Old World atheriniform species. We also discuss the existence of molecular and histological mechanisms to prevent the discrepant differentiation in parts of the gonads because of ambiguous or conflicting environmental and genetic signals and particularly the possibility that the female is the default state in these species.

A fossil-calibrated time-tree of all Australian freshwater fishes

Australian freshwater fishes are a relatively species-poor assemblage, mostly comprising groups derived from older repeated freshwater invasions by marine ancestors, plus a small number of Gondwanan lineages. These taxa are both highly endemic and highly threatened, but a comprehensive phylogeny for Australian freshwater fishes is lacking. This has hampered efforts to study their phylogenetic diversity, distribution of extinction risk, speciation rates, and rates of trait evolution. Here, we present a comprehensive dated phylogeny of 412 Australian fishes. We include all formally recognized freshwater species plus a number of genetically distinct subpopulations, species awaiting formal description, and predominantly brackish-water species that sometimes enter fresh water. The phylogeny was inferred using maximum-likelihood analysis of a multilocus data set comprising six mitochondrial and three nuclear genes from 326 taxa. We inferred the evolutionary timescale using penalized likelihood, then used a statistical approach to add 86 taxa for which no molecular data were available. The time-tree inferred in our study will provide a useful resource for macroecological studies of Australian freshwater fishes by enabling corrections for phylogenetic non-independence in evolutionary and ecological comparative analyses.

Mitochondrial genome architecture and phylogenetic relationships of Odontesthes argentinensis within Atherinomorpha

Silversides are a widely distributed group across South America, with several species occupying marine, freshwater and estuarine environments. Several authors suggest main transitions among these environments took place during Pleistocene, and were accompanied with rapid speciation events. This scenario produced very limited genetic and morphological differentiation among the species. However, most of these surveys have an incomplete coverage of the intraspecific genetic diversity of the taxa studied. In this work, we reconstructed six mitochondrial genomes of O. argentinensis using transcriptomic data, and used them-in combination with several nuclear markers retrieved from the same transcriptomes-to explore the effect of additional coverage of intraspecific diversity of this species in phylogenetic reconstructions. Unlike previous works, phylogenetic analyses failed to identify O. argentinensis as a monophyletic group in relation with closely related taxa. Our results suggest that several species of the genus, especially those related to O. argentinensis, need further taxonomic revision.

"Colonization routes and demographic history of Chirostoma humboldtianum in the central Mexican plateau"

In the present study we evaluate the population structure and potential colonization routes of the silverside Chirostoma humboldtianum through approximate Bayesian computations. Six microsatellite loci were amplified in a total of 288 individuals from six different locations covering the complete geographic distribution of the species. Additionally, two mitochondrial DNA markers, a D loop control region and cytochrome b were amplified in a subset of 107 individuals. The results found with microsatellites allow recovering well-structured populations that have experienced a drastic reduction in the effective population size. On the other hand, mtDNA sequences showed a moderate phylogeographic structure with shared haplotypes between geographic localities and signalsof a slight increase in the effective population size. Finally, the approximate Bayesian computation analysis performed with both datasets suggested a west-to-east colonization route for the species in Central Mexico.

Evidence of Interspecific Chromosomal Diversification in Rainbowfishes (Melanotaeniidae, Teleostei)

Rainbowfishes (Melanotaeniidae) are the largest monophyletic group of freshwater fishes occurring in Australia and New Guinea, with 112 species currently recognised. Despite their high taxonomic diversity, rainbowfishes remain poorly studied from a cytogenetic perspective. Using conventional (Giemsa staining, C banding, chromomycin A3 staining) and molecular (fluorescence in situ hybridisation with ribosomal DNA (rDNA) and telomeric probes) cytogenetic protocols, karyotypes and associated chromosomal characteristics of five species were examined. We covered all major lineages of this group, namely, Running River rainbowfish Melanotaenia sp., red rainbowfish Glossolepisincisus, threadfin rainbowfish Iriatherina werneri, ornate rainbowfish Rhadinocentrus ornatus, and Cairns rainbowfish Cairnsichthys rhombosomoides. All species had conserved diploid chromosome numbers 2n = 48, but karyotypes differed among species; while Melanotaenia sp., G. incisus, and I. werneri possessed karyotypes composed of exclusively subtelo/acrocentric chromosomes, the karyotype of R. ornatus displayed six pairs of submetacentric and 18 pairs of subtelo/acrocentric chromosomes, while C. rhombosomoides possessed a karyotype composed of four pairs of submetacentric and 20 pairs of subtelo/acrocentric chromosomes. No heteromorphic sex chromosomes were detected using conventional cytogenetic techniques. Our data indicate a conserved 2n in Melanotaeniidae, but morphologically variable karyotypes, rDNA sites, and heterochromatin distributions. Differences were observed especially in taxonomically divergent species, suggesting interspecies chromosome rearrangements.

Calcified gill filaments increase respiratory function in fishes

The morphology of fish gills is closely linked to aerobic capacity and tolerance of environmental stressors such as hypoxia. The importance of gill surface area is well studied, but little is known about how the mechanical properties of gill tissues determine function. In some fishes, the bases of the gill filaments are surrounded by a calcified 'sheath' of unknown function. We tested two non-exclusive hypotheses: (i) calcified gill filaments enhance water flow through the gill basket, improving aquatic respiratory function, and (ii) in amphibious fishes, calcification provides support for gills out of water. In a survey of more than 100 species of killifishes and related orders, we found filament calcification was widespread and thus probably arose before the evolution of amphibious lifestyles in killifishes. Calcification also did not differ between amphibious and fully aquatic species, but terrestrial acclimation caused calcium deposition on the filaments of the killifish Kryptolebias marmoratus, suggesting a possible structural role when out of water. We found strong evidence supporting a role for filament calcification in enhancing aquatic respiratory function. First, acclimation to increased respiratory demands (hypoxia, elevated temperatures) induced calcium deposition on the filaments of K. marmoratus. Next, gentle removal of filament calcification decreased branchial resistance to water flow, indicating disruption of gill basket positioning. Thus, the mechanical properties of the gill filaments appear to play an important and previously unappreciated role in determining fish respiratory function.

A comprehensive non-redundant reference transcriptome for the Atlantic silverside Menidia menidia

The Atlantic silverside (Menidia menidia) has been the focus of extensive research efforts in ecology, evolutionary biology, and physiology over the past three decades, but lack of genomic resources has so far hindered examination of the molecular basis underlying the remarkable patterns of phenotypic variation described in this species. We here present the first reference transcriptome for M. menidia. We sought to capture a single representative sequence from as many genes as possible by first using a combination of Trinity and the CLC Genomics Workbench to de novo assemble contigs based on RNA-seq data from multiple individuals, tissue types, and life stages. To reduce redundancy, we passed the combined raw assemblies through a stringent filtering pipeline based both on sequence similarity to related species and computational predictions of transcript quality, condensing an initial set of >480,000 contigs to a final set of 20,998 representative contigs, amounting to a total length of 53.3 Mb. In this final assembly, 91% of the contigs were functionally annotated with putative gene function and gene ontology (GO) terms and/or InterProScan identifiers. The assembly contains complete or nearly complete copies of >95% of 248 highly conserved core genes present in low copy number across higher eukaryotes, and partial copies of another 3.8%, suggesting that our assembly provides relatively comprehensive coverage of the M. menidia transcriptome. The assembly provided here will be an important resource for future research.

Evidence of a landlocked reproducing population of the marine pejerrey Odontesthes argentinensis (Actinopterygii; Atherinopsidae)

In South America, the order Atheriniformes includes the monophyletic genus Odontesthes with 20 species that inhabit freshwater, estuarine and coastal environments. Pejerrey Odontesthes argentinensis is widely distributed in coastal and estuarine areas of the Atlantic Ocean and is known to foray into estuaries of river systems, particularly in conditions of elevated salinity. However, to our knowledge, a landlocked self-sustaining population has never been recorded. In this study, we examined the pejerrey population of Salada de Pedro Luro Lake (south-east of Buenos Aires Province, Argentina) to clarify its taxonomic identity. An integrative taxonomic analysis based on traditional meristic, landmark-based morphometrics and genetic techniques suggests that the Salada de Pedro Luro pejerrey population represents a novel case of physiological and morphological adaptation of a marine pejerrey species to a landlocked environment and emphasises the environmental plasticity of this group of fishes.

Origins of Afrotropical freshwater fishes

The Afrotropics house a diverse freshwater ichthyofauna with > 3000 species, almost all of which are endemic. Recent progress in dated phylogenetics and palaeontology of several groups of Afrotropical freshwater fishes (AFFs) has allowed the testing of palaeoecology- and palaeogeography-based hypotheses explaining their early presence in Africa. Seven hypotheses were tested for 37 most-inclusive monophyletic groups of AFFs. Results indicated that ten lineages originated from direct, but asynchronous, marine-to-freshwater shifts. These lineages contribute < 2% to the current AFF species richness. Eleven lineages colonized the Afrotropics from the Orient after the Afro-Arabian plate collided with Eurasia in the early Oligocene. These lineages contribute ~20% to the total diversity. There are seven sister relationships between Afrotropical and Neotropical taxa. For only three of them (4% of the species diversity), the continental drift vicariance hypothesis was not rejected. Distributions of the other four younger trans-Atlantic lineages are better explained by post-drifting long-distance dispersal. In those cases, I discuss the possibility of dispersal through the Northern Hemisphere as an alternative to direct trans-Atlantic dispersal. The origins of ten AFF lineages, including the most species-rich Pseudocrenilabrinae (> 1100 species), are not yet established with confidence.

The Duplicated Y-specific amhy Gene Is Conserved and Linked to Maleness in Silversides of the Genus Odontesthes

Sex-determining genes have been successively isolated in several teleosts. In Odontesthes hatcheri and O. bonariensis, the amhy gene has been identified as a master sex-determining gene. However, whether this gene is conserved along related species is still unknown. In this study, the presence of amhy and its association with phenotypic sex was analyzed in 10 species of Odontesthes genus. The primer sets from O. hatcheri that amplify both amhs successfully generated fragments that correspond to amha and amhy in all species. The full sequences of amhy and amha isolated for four key species revealed higher identity values among presumptive amhy, including the 0.5 Kbp insertion in the third intron and amhy-specific insertions/deletions. Amha was present in all specimens, regardless of species and sex, whereas amhy was amplified in most but not all phenotypic males. Complete association between amhy-homologue with maleness was found in O. argentinensis, O. incisa, O. mauleanum, O. perugiae, O. piquava, O. regia, and O. smitti, whereas O. humensis, O. mirinensis, and O. nigricans showed varied degrees of phenotypic/genotypic sex mismatch. The conservation of amhy gene in Odontesthes provide an interesting framework to study the evolution and the ecological interactions of genotypic and environmental sex determination in this group.

Contrasting genomic shifts underlie parallel phenotypic evolution in response to fishing

Humans cause widespread evolutionary change in nature, but we still know little about the genomic basis of rapid adaptation in the Anthropocene. We tracked genomic changes across all protein-coding genes in experimental fish populations that evolved pronounced shifts in growth rates due to size-selective harvest over only four generations. Comparisons of replicate lines show parallel allele frequency shifts that recapitulate responses to size-selection gradients in the wild across hundreds of unlinked variants concentrated in growth-related genes. However, a supercluster of genes also rose rapidly in frequency and dominated the evolutionary dynamic in one replicate line but not in others. Parallel phenotypic changes thus masked highly divergent genomic responses to selection, illustrating how contingent rapid adaptation can be in the face of strong human-induced selection.

Characterization of the karyotype and accumulation of repetitive sequences in Australian Darling hardyhead Craterocephalus amniculus (Atheriniformes, Teleostei)

Belonging to the order Atheriniformes, Craterocephalus is one of the most widespread genera of freshwater fishes in Australia, spanning along the northern coast from central Western Australia to central New South Wales and across the Murray-Darling and Lake Eyre basins. In this study, both conventional cytogenetic techniques (Giemsa, C-banding, CMA₃/DAPI staining), and fluorescence in situ hybridization (FISH) with telomeric DNA and rDNA probes were used to examine the karyotypes and other chromosomal characteristics of Darling hardyhead (Craterocephalus amniculus) from New South Wales, Australia. We identified a diploid chromosome number 2n = 48 (NF = 58) in all studied individuals. FISH with rDNA probes showed a nonsyntenic pattern, with signals on one pair of subtelocentric chromosomes for 5S rDNA and one pair of submetacentric chromosomes for 28S rDNA. C-banding displayed the accumulation of constitutive heterochromatin in the centromeric regions of approximately 40 chromosomes. CMA₃/DAPI fluorescence staining revealed extremely GC-rich signals in the pericentromeric region of one submetacentric chromosomal pair with size polymorphism. We detected telomeric signals at the end of all chromosomes and no interstitial signals.

Inferring boundaries among fish species of the new world silversides (Atherinopsidae; genus Odontesthes): new evidences of incipient speciation between marine and brackish populations of Odontesthes argentinensis

Species of new world silversides (Actinopterygii; Atherinopsidae; genus Odontesthes) possess economic relevance, biological interest and ecological importance. In the present paper we: (A) investigate the molecular diversity in marine species of Odontesthes from the South West Atlantic Ocean (SWAO), and analyse their interspecific relationships and divergence by means of DNA Barcoding, including its freshwater congeners, as well. (B) Explore the suitability of DNA Barcoding to analyse the diversity and distribution of haplotypes in Odontesthes argentinensis, the only well documented marine species from the SWAO that exhibit putative estuarine and marine populations. Molecular analysis revealed 100% of agreement between morphological identification and molecular identity. Odontesthes argentinensis, Odontesthes platensis, Odontesthes smitti, Odontesthes nigricans and Odontesthes incisa were assigned to five different barcode index numbers (BINs). Maximum-likelihood analysis showed that all marine species of Odontesthes clustered separately in a unique monophyletic phylogroup, comprising five well defined haplogroups, with genetic divergence between groups ranging from 2.75 to 7.11%. The genetic analysis including freshwater congeners showed that O. incisa clustered alone occupying a basal position. The F_st pairwise comparisons within O. argentinensis support the existence of three population groups: one conformed by Mar Chiquita Lagoon (MCh) specimens, and the others by Mar del Plata/Mar Chiquita coast and San Blas Bay coastal specimens, respectively. The AMOVA showed significant overall differentiation (F_st = 0.238; p = 0.00001) for the entire data set. The previous/present evidence is discussed, and strongly suggests that incipient speciation is occurring in O. argentinensis argentinean populations, and specimens from MCh would be considered at present as the leading candidate of a marine to freshwater incipient speciation event.

Roundup® Herbicide Decreases Quality Parameters of Spermatozoa of Silversides Odontesthes Humensis

The silverside (Odontesthes humensis) is a very interesting model for toxicological studies due its high sensitivity and need for good water quality. The aim of this study was to evaluate the effects of Roundup on spermatozoa of O. humensis, after acute exposure. The fish were exposed to 0 and 7.8 mg L-1 (a.e.) of glyphosate, respectively. Through computer-assisted sperm analysis, a significant decrease in concentration, total and progressive motility, average path distance, straight line distance, path average velocity, curved line velocity, straight line velocity linearity, wobble, amplitude of lateral head displacement, cross beat frequency, and motility period of silverside spermatozoa exposed to Roundup was observed. Also, increase in membrane fluidity, ROS production and lipid peroxidation and a decrease in the mitochondrial functionality was observed in spermatozoa of Roundup exposed silversides. It was demonstrated that Roundup exposure in a concentration that can be achieve in natural water bodies soon after its application in fields is able to cause losses in several sperm quality parameters, consequently decreasing the fertilization potential of O. humensis spermatozoa.

Gene and Blood Analysis Reveal That Transfer from Brackish Water to Freshwater Is More Stressful to the Silverside Odontesthes humensis

Silversides are fish that inhabit marine coastal waters, coastal lagoons, and estuarine regions in southern South America. The freshwater (FW) silversides have the ability to tolerate salinity variations. Odontesthes humensis have similar habitats and biological characteristics of congeneric O. bonariensis, the most studied silverside species and with great economic importance. Studies revealed that O. bonariensis is not fully adapted to FW, despite inhabiting hyposmotic environments in nature. However, there is little information about stressful environments for cultivation of silverside O. humensis. Thus, the aim of this study was to evaluate the stress and osmoregulation responses triggered by the osmotic transfers on silverside O. humensis. Silversides were acclimated to FW (0 ppt) and to brackish water (BW, 10 ppt) and then they were exposed to opposite salinity treatment. Silverside gills and blood were sampled on pre-transfer (D0) and 1, 7, and 15 days (D1, D7, and D15) after changes in environmental salinity, the expression levels of genes atp1a3a, slc12a2b, kcnh1, and hspa1a were determined by quantitative reverse transcription-PCR for evaluation of osmoregulatory and stress responses. Furthermore, glycemia, hematocrit, and osmolality were also evaluated. The expression of atp1a3a was up- and down-regulated at D1 after the FW-BW and BW-FW transfers, respectively. Slc12a2b was up-regulated after FW-BW transfer. Similarly, kcnh1 and hspa1a were up-regulated at D1 after the BW-FW transfer. O. humensis blood osmolality decreased after the exposure to FW. It remained stable after exposure to BW, indicating an efficient hyposmoregulation. The glycemia had a peak at D1 after BW-FW transfer. No changes were observed in hematocrit. The return to the pre-transfer levels at D7 after the significant increases in responses of almost all evaluated molecular and blood parameters indicated that this period is enough for acclimation to the experimental conditions. In conclusion, our results suggest that BW-FW transfer is more stressful to O. humensis than FW-BW transfer and the physiology of O. humensis is only partially adapted to FW.

Insight into the population structure of hardhead silverside, Atherinomorus stipes (Teleostei: Atherinidae), in Belize and the Florida Keys using nd2

Little is known about the natural history, biology, and population genetic structure of the Hardhead Silverside, Atherinomorus stipes, a small schooling fish found around islands throughout the Caribbean. Our field observations of A. stipes in the cays of Belize and the Florida Keys found that populations tend to be in close association with the shoreline in mangrove habitats. Due to this potential island‐based population structuring, A. stipes represents an ideal system to examine questions about gene flow and isolation by distance at different geographic scales. For this study, the mitochondrial gene nd2 was amplified from 394 individuals collected from seven different Belizean Cays (N = 175) and eight different Floridian Keys (N = 219). Results show surprisingly high haplotype diversity both within and between island‐groups, as well as a high prevalence of unique haplotypes within each island population. The results are consistent with models that require gene flow among populations as well as in situ evolution of rare haplotypes. There was no evidence for an isolation by distance model. The nd2 gene tree consists of two well‐supported monophyletic groups: a Belizean‐type clade and a Floridian‐type clade, indicating potential species‐level differentiation.

Pulsed evolution shaped modern vertebrate body sizes

The diversity of forms found among animals on Earth is striking. Despite decades of study, it has been difficult to reconcile the patterns of diversity seen between closely related species with those observed when studying single species on ecological timescales. We propose a set of models, called Lévy processes, to attempt to reconcile rapid evolution between species with the relatively stable distributions of phenotypes seen within species. These models, which have been successfully used to model stock market data, allow for long periods of stasis followed by bursts of rapid change. We find that many vertebrate groups are well fitted by Lévy models compared with models for which traits evolve toward a stationary optimum or evolve in an incremental and wandering manner.

The relative importance of different modes of evolution in shaping phenotypic diversity remains a hotly debated question. Fossil data suggest that stasis may be a common mode of evolution, while modern data suggest some lineages experience very fast rates of evolution. One way to reconcile these observations is to imagine that evolution proceeds in pulses, rather than in increments, on geological timescales. To test this hypothesis, we developed a maximum-likelihood framework for fitting Lévy processes to comparative morphological data. This class of stochastic processes includes both an incremental and a pulsed component. We found that a plurality of modern vertebrate clades examined are best fitted by pulsed processes over models of incremental change, stationarity, and adaptive radiation. When we compare our results to theoretical expectations of the rate and speed of regime shifts for models that detail fitness landscape dynamics, we find that our quantitative results are broadly compatible with both microevolutionary models and observations from the fossil record.

Phylogenetic classification of bony fishes

BACKGROUND

Fish classifications, as those of most other taxonomic groups, are being transformed drastically as new molecular phylogenies provide support for natural groups that were unanticipated by previous studies. A brief review of the main criteria used by ichthyologists to define their classifications during the last 50 years, however, reveals slow progress towards using an explicit phylogenetic framework. Instead, the trend has been to rely, in varying degrees, on deep-rooted anatomical concepts and authority, often mixing taxa with explicit phylogenetic support with arbitrary groupings. Two leading sources in ichthyology frequently used for fish classifications (JS Nelson's volumes of Fishes of the World and W. Eschmeyer's Catalog of Fishes) fail to adopt a global phylogenetic framework despite much recent progress made towards the resolution of the fish Tree of Life. The first explicit phylogenetic classification of bony fishes was published in 2013, based on a comprehensive molecular phylogeny ( www.deepfin.org ). We here update the first version of that classification by incorporating the most recent phylogenetic results.

RESULTS

The updated classification presented here is based on phylogenies inferred using molecular and genomic data for nearly 2000 fishes. A total of 72 orders (and 79 suborders) are recognized in this version, compared with 66 orders in version 1. The phylogeny resolves placement of 410 families, or ~80% of the total of 514 families of bony fishes currently recognized. The ordinal status of 30 percomorph families included in this study, however, remains uncertain (incertae sedis in the series Carangaria, Ovalentaria, or Eupercaria). Comments to support taxonomic decisions and comparisons with conflicting taxonomic groups proposed by others are presented. We also highlight cases were morphological support exist for the groups being classified.

CONCLUSIONS

This version of the phylogenetic classification of bony fishes is substantially improved, providing resolution for more taxa than previous versions, based on more densely sampled phylogenetic trees. The classification presented in this study represents, unlike any other, the most up-to-date hypothesis of the Tree of Life of fishes.

Transcriptomic differentiation underlying marine-to-freshwater transitions in the South American silversides Odontesthes argentinensis and O. bonariensis (Atheriniformes)

Salinity gradients are critical habitat determinants for freshwater organisms. Silverside fishes in the genus Odontesthes have recently and repeatedly transitioned from marine to freshwater habitats, overcoming a strong ecological barrier. Genomic and transcriptomic changes involved in this kind of transition are only known for a few model species. We present new data and analyses of gene expression and microbiome composition in the gills of two closely related silverside species, marine O. argentinensis and freshwater O. bonariensis and find more than three thousand transcripts differentially expressed, with osmoregulatory/ion transport genes and immune genes showing very different expression patterns across species. Interspecific differences also involve more than one thousand transcripts with nonsynonymous SNPs in the coding sequences, most of which were not differentially expressed. In addition to characterizing gill transcriptomes from wild‐caught marine and freshwater fishes, we test experimentally the response to salinity increases by O. bonariensis collected from freshwater habitats. Patterns of expression in gill transcriptomes of O. bonariensis exposed to high salinity do not resemble O. argentinensis mRNA expression, suggesting lack of plasticity for adaptation to marine conditions in this species. The diversity of functions associated with both the differentially expressed set of transcripts and those with sequence divergence plus marked microbiome differences suggest that multiple abiotic and biotic factors in marine and freshwater habitats are driving transcriptomic differences between these species.

Pinpointing cryptic borders: Fine-scale phylogeography and genetic landscape analysis of the Hormogaster elisae complex (Oligochaeta, Hormogastridae)

Spatial and temporal aspects of the evolution of cryptic species complexes have received less attention than species delimitation within them. The phylogeography of the cryptic complex Hormogaster elisae (Oligochaeta, Hormogastridae) lacks knowledge on several aspects, including the small-scale distribution of its lineages or the palaeogeographic context of their diversification. To shed light on these topics, a dense specimen collection was performed in the center of the Iberian Peninsula - resulting in 28 new H. elisae collecting points, some of them as close as 760m from each other- for a higher resolution of the distribution of the cryptic lineages and the relationships between the populations. Seven molecular regions were amplified: mitochondrial subunit 1 of cytochrome c oxidase (COI), 16S rRNA and tRNA Leu, Ala, and Ser (16S t-RNAs), one nuclear ribosomal gene (a fragment of 28S rRNA) and one nuclear protein-encoding gene (histone H3) in order to infer their phylogenetic relationships. Different representation methods of the pairwise divergence in the cytochrome oxidase I sequence (heatmap and genetic landscape graphs) were used to visualize the genetic structure of H. elisae. A nested approach sensu Mairal et al. (2015) (connecting the evolutionary rates of two datasets of different taxonomic coverage) was used to obtain one approximation to a time-calibrated phylogenetic tree based on external Clitellata fossils and a wide molecular dataset. Our results indicate that limited active dispersal ability and ecological or biotic barriers could explain the isolation of the different cryptic lineages, which never co-occur. Rare events of long distance dispersal through hydrochory appear as one of the possible causes of range expansion.

Lessons from the first dorsal fin in atheriniforms-A new mode of dorsal fin development and its phylogenetic implications

The median fins in extant actinopterygians are the product of millions of years of evolution. During this time, different developmental patterns for the dorsal and anal fins emerged leading to a high variation in median fin morphology and ontogeny. In this study, the development of anal and dorsal fins in atheriniforms is described and its consequences for the current phylogenetic hypothesis are discussed. Developmental series of five atheriniform species were investigated using clearing and staining as well as antibody staining. The skeletal elements of the second dorsal fin and the anal fin emerge in a bidirectional pattern. The first dorsal fin, however, arises separately in front of the second dorsal fin after this one is almost completely formed. The pterygiophores of the first dorsal fin, including the interdorsal pterygiophores, develop from caudal to rostral, but the fin‐spines of the first dorsal fin form in the opposite direction. This new mode of fin development has been found in all examined atheriniform species with two dorsal fins. Several morphological characters of atheriniforms, including interdorsal pterygiophores, are also found in one other taxon: the Mugiliformes. Thus, several dorsal fin characteristics may provide evidence for a closer relationship of these two taxa.

First DNA Barcode Reference Library for the Identification of South American Freshwater Fish from the Lower Paraná River

Valid fish species identification is essential for biodiversity conservation and fisheries management. Here, we provide a sequence reference library based on mitochondrial cytochrome c oxidase subunit I for a valid identification of 79 freshwater fish species from the Lower Paraná River. Neighbour-joining analysis based on K2P genetic distances formed non-overlapping clusters for almost all species with a ≥99% bootstrap support each. Identification was successful for 97.8% of species as the minimum genetic distance to the nearest neighbour exceeded the maximum intraspecific distance in all these cases. A barcoding gap of 2.5% was apparent for the whole data set with the exception of four cases. Within-species distances ranged from 0.00% to 7.59%, while interspecific distances varied between 4.06% and 19.98%, without considering Odontesthes species with a minimum genetic distance of 0%. Sequence library validation was performed by applying BOLDs BIN analysis tool, Poisson Tree Processes model and Automatic Barcode Gap Discovery, along with a reliable taxonomic assignment by experts. Exhaustive revision of vouchers was performed when a conflicting assignment was detected after sequence analysis and BIN discordance evaluation. Thus, the sequence library presented here can be confidently used as a benchmark for identification of half of the fish species recorded for the Lower Paraná River.