Disentangling historical relationships within Poeciliidae (Teleostei: Cyprinodontiformes) using ultraconserved elements

Poeciliids (Cyprinodontiformes: Poeciliidae), commonly known as livebearers, are popular fishes in the aquarium trade (e.g., guppies, mollies, swordtails) that are widely distributed in the Americas, with 274 valid species in 27 genera. This group has undergone various taxonomic changes recently, spurred by investigations using traditional genetic markers. Here we used over 1,000 ultraconserved loci to infer the relationships within Poeciliidae in the first attempt at understanding their diversification based on genome-scale data. We explore gene tree discordance and investigate potential incongruence between concatenation and coalescent inference methods. Our aim is to examine the influence of incomplete lineage sorting and reticulate evolution on the poeciliids' evolutionary history and how these factors contribute to the observed gene tree discordace. Our concatenated and coalescent phylogenomic inferences recovered four major clades within Poeciliidae. Most supra-generic level relationships we inferred were congruent with previous molecular studies, but we found some disagreements; the Middle American taxa Phallichthys and Poecilia (Mollienesia) were recovered as non-monophyletic, and unlike other recent molecular studies, we recovered Brachyrhaphis as monophyletic. Our study is the first to provide signatures of reticulate evolution in Poeciliidae at the family level; however, continued finer-scale investigations are needed to understand the complex evolutionary history of the family along with a much-needed taxonomic re-evaluation.

Investigating the utility of Anchored Hybrid Enrichment data to investigate the relationships among the Killifishes (Actinopterygii: Cyprinodontiformes), a globally distributed group of fishes

The Killifishes (Cyprinodontiformes) are a diverse and well-known group of fishes that contains sixteen families inclusive of Anablepidae, Aphaniidae Aplocheilidae, Cubanichthyidae, Cyprinodontidae, Fluviphylacidae, Fundulidae, Goodeidae, Nothobranchiidae, Orestiidae, Pantanodontidae, Poeciliidae, Procatopodidae, Profundulidae, Rivulidae, and Valenciidae and more than 1,200 species that are globally distributed in tropical and temperate, freshwater and estuarine habitats. The evolutionary relationships among the families within the group, based on different molecular and morphological data sets, have remained uncertain. Therefore, the objective of this study was to use a targeted approach, anchored hybrid enrichment, to investigate the phylogenetic relationships among the families within the Cyprindontiformes. This study included more than 100 individuals, representing all sixteen families within the Cyprinodontiformes, including many recently diagnosed families. We recovered an average of 244 loci per individual. These data were submitted to phylogenetic analyses (RaxML and ASTRAL) and although we recovered many of the same relationships as in previous studies of the group, several novel sets of relationships for other families also were recovered. In addition, two well-established clades (Suborders Cyprinodontoidei and Aplocheilodei) were recovered as monophyletic and are in agreement with most previous studies. We also assessed the degree of gene tree discordance in our dataset to evaluate support for alternative topological hypotheses for interfamilial relationships within the Cyprinodontiformes using a variety of different analyses. The results from this study will provide a robust, historical framework needed to investigate a plethora of biogeographic, taxonomic, ecological, and physiological questions for this group of fishes.

Genomic data support the taxonomic validity of Middle American livebearers Poeciliopsis gracilis and Poeciliopsis pleurospilus (Cyprinodontiformes: Poeciliidae)

Poeciliopsis (Cyprinodontiformes: Poeciliidae) is a genus comprised of 25 species of freshwater fishes. Several well-known taxonomic uncertainties exist within the genus, especially in relation to the taxonomic status of Poeciliopsis pleurospilus and P. gracilis. However, to date, no studies have been conducted to specifically address the taxonomic status of these two species. The goal of this study was to examine the taxonomic validity of P. pleurospilus and P. gracilis using genomic data (ddRADseq) in phylogenetic, population genetic, and species delimitation frameworks. Multiple analyses support the recognition of both taxa as distinct species and also permits us to revise their respective distributions. A species delimitation analysis indicates that P. pleurospilus and P. gracilis are distinct species, each of which consists of two distinct lineages that are geographically structured. Phylogenetic and population genetic analyses provide clear evidence that individuals of P. gracilis are distributed north and west of the Isthmus of Tehuantepec in both Pacific and Atlantic river systems in Mexico, whereas individuals of P. pleurospilus are distributed in both Atlantic and Pacific river systems south and east of the Isthmus of Tehuantepec, from southern Mexico to Honduras.

Molecular systematics of the Awaous banana complex (River gobies; Teleostei: Oxudercidae)

Diadromous fishes can exhibit interesting evolutionary and population-level patterns given their use of freshwater and marine environments as part of their life histories. The River goby genus Awaous are prominent members of riverine ichthyofaunas and occur throughout Atlantic and Pacific slopes of the Americas from the southern United States to Ecuador and Brazil. Here we study the widespread and polymorphic Awaous banana complex to assess phylogeographic patterns and test previous hypotheses that all populations of this species in the Americas belong to the same species. Analysis of sequence data based on the mitochondrial cytochrome oxidase I gene shows multiple clades within the Atlantic and Pacific basins, which correspond to previously described species. Additionally, haplotype analysis demonstrates unique and unconnected networks between these species. Within these clades we document biogeographic patterns that are congruent with results of other co-occurring diadromous species, as well as a novel biogeographic pattern for the region. Our results support the recognition of distinct species of Awaous in the Atlantic (A. banana and A. tajasica) and Pacific (A. transandeanus) basins. These results are concordant with previously established morphological characters permitting the separation of these species.

Genetic differentiation among populations of the blackfin goodea Goodea atripinnis (Cyprinodontiformes: Goodeidae): implications for its evolutionary history

Central Mexico is characterized by a complex topography that is the result of historic and contemporary tectonic and climatic factors. These events have influenced the evolutionary history of numerous freshwater fishes in the region. Nonetheless, recent studies have shown that life-history traits and ecological characteristics of species may influence dispersal capabilities and the degree of genetic connectivity. Goodea (Cyprinodontiformes: Goodeidae) is one of the most widely distributed and environmentally tolerant genera of goodeids. In this study, the authors analysed variation in the mitochondrial cytochrome b gene to evaluate the phylogeographic relationships, genetic structure, genetic diversity and demographic history of Goodea from across its distribution range. They found low genetic differentiation and identified shared haplotypes among several regions. Geographic segregation was found in samples southwest and northeast of the Lower Lerma region, with some internal isolated groups showing phylogeographic differentiation and unique haplotypes. The AMOVA best explained genetic structure when grouped by haplogroups rather than when grouped by recognized biogeographic regions. Several regions showed null genetic diversity, raising the possibility of dispersal mediated by humans. Finally, Bayesian Skyline Plot analysis showed a population expansion for the Southwest haplogroup, except for the Armería population and sub-group II of the Northeast haplogroup. All this suggests a recent colonization of Goodea atripinnis throughout some of the biogeographic regions currently inhabited by this species.

Distribution and current conservation status of the Mexican Goodeidae (Actinopterygii, Cyprinodontiformes)

The current distribution and abundance of the 40 species of Goodeidae fishes known from Mexico are described, and a total of 84 Evolutionarily Significant Units (ESUs) is designated within these species. Two species and four ESUs are likely extinct with no captive populations, and three species and eight ESUs are probably extinct in the wild but have at least one captive population in Mexico, the United States, or Europe. Of the 35 extant species, the analyses indicate that nine should be considered as critically endangered, 14 as endangered, nine as vulnerable, and only three as least concern. Twenty-seven of these species have experienced substantial declines in distribution or abundance or both since 2000, and only eight appear to have remained relatively stable. Of the 72 extant ESUs, our analyses indicate that 29 should be considered as critically endangered, 21 as endangered, 18 as vulnerable, and only four as least concern. Brief summaries of the historic and current distributions and abundance of each species are provided, as well as ESU. Three strategies are recommended to conserve Mexican goodeids: protect the best-quality remaining habitats where goodeids still persist, restore degraded habitat and re-introduce species or ESUs where practical, and establish captive populations to ensure continued survival of the many species and ESUs that will almost inevitably go extinct in the coming years. Limited resources require cooperation and collaboration between scientists, conservationists, and aquarium hobbyists for successful captive maintenance.

Species or population? Systematic status of Vieja coatlicue (Teleostei: Cichlidae)

ABSTRACT Pacific and Atlantic-slope rivers of the Isthmus of Tehuantepec in Mexico contain two cichlid species of the genus Vieja relevant to the current study, V. zonata and V. coatlicue, respectively. The Atlantic-slope species was initially considered a population of uncertain taxonomic status and recently described as a distinct species based on three putatively diagnostic morphological characters. The objective of this study was to combine new and existing genetic data, along with reanalysis of morphological characters, to test the distinctiveness of V. coatlicue. Little genetic variability and no geographic structuring were recovered. Additionally, putatively diagnostic morphological character states were present across both species, failing to separate the forms. The synthesis of molecular and morphological data supports the recognition of V. coatlicue as a junior synonym of V. zonata.

Disentangling the drivers of diversification in an imperiled group of freshwater fishes (Cyprinodontiformes: Goodeidae)

BACKGROUND

One of the most perplexing questions in evolutionary biology is why some lineages diversify into many species, and others do not. In many cases, ecological opportunity has played an important role, leading to diversification along trophic or habitat-based axes. The Goodeidae (Teleostomi: Cyprinodontiformes) are a family of freshwater fishes with two subfamilies: Goodeinae (42 species, viviparous, heterogeneous habitats, Mesa Central of Mexico) and Empetrichthyinae (4 species, oviparous, homogeneous habitats, Great Basin of the United States). These discrepant sets of characteristics and their sister-group relationship make the goodeids amenable to a comparative study of diversification. We gathered lateral body images from more than 1600 specimens of all extant species in the family. Geometric morphometric, and phylogenetic comparative analyses were used to address whether higher species diversity correlates with higher rates of morphological shape evolution and whether there are differences in functional/habitat modules between the two subfamilies.

RESULTS

This study recovered a higher rate of overall body shape evolution in the Goodeinae that is nearly double in magnitude compared to the Empetrichthyinae. A modularity test indicated that the Goodeinae displayed elevated rates of morphological evolution in comparison to the Empetrichthyinae when only trunk (locomotor) regions were compared between subfamilies. No significant differences in evolutionary shape rates were recovered when the trophic (head) regions were compared between subfamilies.

DISCUSSION

These results support the hypothesis that Mexican goodeids radiated via an ecological opportunity scenario into a wide-array of novel habitats in the island-like Mesa Central as evidenced by their high rate of shape evolution, relative to the Empetrichthyinae. This study quantitatively unraveled the drivers of evolution and eliminated trophic specialization as a driving force within the Goodeidae.

CONCLUSIONS

A combination of phylogenetic and morphometric data, and phylogenetic comparative analyses were used to examine body shape rate evolution within the Goodeidae. Results support the hypothesis that species in the subfamily Goodeinae on the central Mexican plateau had a higher rate of body shape evolution relative to its sister subfamily Empetrichthyinae in the Great Basin suggesting that the Goodeinae diversified via an ecological opportunity scenario along habitat, rather than trophic axes.

Assessing population-level morphometric variation of the Mountain Mullet Agonostomus monticola (Teleostei: Mugilidae) across its Middle American distribution

ABSTRACT Population-level morphometric variation of the Mountain Mullet (Agonostomus monticola) was assessed in 419 adult specimens from 25 sample sites (river basins) across its Middle American distribution (Pacific and Atlantic-Caribbean drainages). This analysis was based on 36 standardized linear measurements and 19 landmarks on geometric morphometrics approach. Discriminant function analysis (DFA) revealed 19 linear morphological characters with significant variation among groups. Geometrically, the most notable changes were associated to the curvature of the frontal region of the head, the anterior and posterior insertion of the first dorsal and anal fins. The resulting grouping based on the DFA and geometric morphometrics techniques (Pacific-A, Pacific-B and NE México-Caribbean) were similar to those previously recovered by genetic techniques, where the Pacific-B (Ayuquila river basin) was the most different group. Our results provide morphological evidence for considering Agonostomus monticola as a complex of evolutionary entities, represented by two forms in the Pacific Ocean and another in the Atlantic Ocean.

Let's jump in: A phylogenetic study of the great basin springfishes and poolfishes, Crenichthys and Empetrichthys (Cyprinodontiformes: Goodeidae)

North America’s Great Basin has long been of interest to biologists due to its high level of organismal endemicity throughout its endorheic watersheds. One example of such a group is the subfamily Empetricthyinae. In this paper, we analyzed the relationships of the Empetrichtyinae and assessed the validity of the subspecies designations given by Williams and Wilde within the group using concatenated phylogenetic tree estimation and species tree estimation. Samples from 19 populations were included covering the entire distribution of the three extant species of Empetricthyinae–Crenichthys nevadae, Crenichthys baileyi and Empetricthys latos. Three nuclear introns (S8 intron 4, S7 intron 1, and P0 intron 1) and one mitochondrial gene (Cytb) were sequenced for phylogenetic analysis. Using these sequences, we generated two separate hypotheses of the evolutionary relationships of Empetrichtyinae- one based on the mitochondrial data and one based on the nuclear data using Bayesian phylogenetics. Haplotype networks were also generated to look at the relationships of the populations within Empetrichthyinae. After comparing the two phylogenetic hypotheses, species trees were generated using *BEAST with the nuclear data to further test the validity of the subspecies within Empetrichthyinae. The mitochondrial analyses supported four lineages within C. baileyi and 2 within C. nevadae. The concatenated nuclear tree was more conserved, supporting one clade and an unresolved polytomy in both species. The species tree analysis supported the presence of two species within both C. baileyi and C. nevadae. Based on the results of these analyses, the subspecies designations of Williams and Wilde are not valid, rather a conservative approach suggests there are two species within C. nevadae and two species within C. baileyi. No structure was found for E. latos or the populations of Empetricthyinae. This study represents one of many demonstrating the invalidity of subspecies and their detriment to species identification, conservation, and understanding.

Shaping up: a geometric morphometric approach to assemblage ecomorphology

This study adopts an ecomorphological approach to test the utility of body shape as a predictor of niche relationships among a stream fish assemblage of the Tickfaw River (Lake Pontchartrain Basin) in southeastern Louisiana, U.S.A. To examine the potential influence of evolutionary constraints, analyses were performed with and without the influence of phylogeny. Fish assemblages were sampled throughout the year, and ecological data (habitat and tropic guild) and body shape (geometric morphometric) data were collected for each fish specimen. Multivariate analyses were performed to examine relationships and differences between body shape and ecological data. Results indicate that a relationship exists between body shape and trophic guild as well as flow regime, but no significant correlation between body shape and substratum was found. Body shape was a reliable indicator of position within assemblage niche space.

Stable isotope evidence for trophic overlap of sympatric Mexican Lake Chapala silversides (Teleostei: Atherinopsidae: Chirostoma spp.)

We explore the trophic role that a diverse sympatric group of fishes in the genus Chirostoma play in a large, shallow lake in central Mexico, Lake Chapala. We use δ13C and δ15N stable isotope - based food web analyses to explore how they relate to other components of the Lake Chapala ecosystem. We find five Chirostoma species in top trophic levels of the Chapala food web compared to other fishes, relying on a combination of zooplankton, fish and benthic resources as energy sources. Food web metric analyses showed generally overlapping trophic niches for members of Chirostoma, especially in terms of δ13C. However, C. jordani had lower mean δ15N isotopic values than C. promelas. As a group, "pescados blancos" (C. sphyraena and C. promelas) also had higher δ15N signatures than "charales" (C. consocium, C. jordani and C. labarcae) reflecting greater piscivory, but these differences were not strong for all food web metrics used. Trophic overlap among species of Chirostoma in Lake Chapala raises questions about the forces that might have led to a morphologically diverse but functionally similar and monophyletic group of species.

Black liquor and the hangover effect: fish assemblage recovery dynamics following a pulse disturbance

Anthropogenic perturbations impact aquatic systems causing wide-ranging responses, from assemblage restructuring to assemblage recovery. Previous studies indicate the duration and intensity of disturbances play a role in the dynamics of assemblage recovery. In August 2011, the Pearl River, United States, was subjected to a weak black liquor spill from a paper mill which resulted in substantial loss of fish in a large stretch of the main channel. We quantified resilience and recovery of fish assemblage structure in the impacted area following the event. We compared downstream (impacted) assemblages to upstream (unimpacted) assemblages to determine initial impacts on structure. Additionally, we incorporated historic fish collections (1988–2011) to examine impacts on assemblage structure across broad temporal scales. Based on NMDS, upstream and downstream sites generally showed similar assemblage structure across sample periods with the exception of the 2 months postdischarge, where upstream and downstream sites visually differed. Multivariate analysis of variance (PERMANOVA) indicated significant seasonal variation among samples, but found no significant interaction between impacted and unimpacted assemblages following the discharge event. However, multivariate dispersion (MVDISP) showed greater variance among assemblage structure following the discharge event. These results suggest that 2 months following the disturbance represent a time period of stochasticity in regard to assemblage structure dynamics, and this was followed by rapid recovery. We term this dynamic the “hangover effect” as it represents the time frame from the cessation of the perturbation to the assemblage's return to predisturbance conditions. The availability and proximity of tributaries and upstream refugia, which were not affected by the disturbance, as well as the rapid recovery of abiotic parameters likely played a substantial role in assemblage recovery. This study not only demonstrates rapid recovery in an aquatic system, but further demonstrates the value of continuous, long-term, data collections which enhance our understanding of assemblage dynamics.

Scale-stacking in the live-bearer Belonesox belizanus (Poeciliidae)

An incidence of scale-stacking was found in a museum collection of Belonesox belizanus (Poeciliidae) from the Río Tesechocán, Veracruz, México. This phenomenon is thought to be unique to lepidophagous East African cichlids, although no studies have addressed the behaviour in any detail. This finding in B. belizanus, the first in Cyprinodontiformes, indicates Belonesox is able to manipulate scales into stacks in the buccal cavity.

It's a family matter: molecular phylogenetics of Atheriniformes and the polyphyly of the surf silversides (family: Notocheiridae)

Phylogenetic relationships among families of Atheriniformes have long been problematic. The affinities of one of the most enigmatic lineages, surf silversides (Notocheiridae), have proven particularly elusive due to this taxon's unique morphology and rarity in museum collections. In this study, we use mitochondrial and nuclear sequence data to generate a phylogeny for seven of the eight families of Atheriniformes. Our results reveal that four families within Atheriniformes (Atherinopsidae, Notocheiridae, Atherinidae, Melanotaeniidae) are not monophyletic. Most notably, Notocheiridae is polyphyletic, with Notocheirus hubbsi nested within New World silversides (Atherinopsidae), while members of Iso are sister to all other Old World Atheriniforms. These data suggest that the unique morphology of Notocheirus and Iso is a result of adaptive convergent evolution to the high-energy surf habitat where these species live.

Permanent genetic resources added to Molecular Ecology Resources Database 1 October 2010-30 November 2010

This article documents the addition of 277 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Ascochyta rabiei, Cambarellus chapalanus, Chionodraco hamatus, Coptis omeiensis, Cynoscion nebulosus, Daphnia magna, Gerbillus nigeriae, Isurus oxyrinchus, Lates calcarifer, Metacarcinus magister, Oplegnathus fasciatus, Pachycondyla verenae, Phaethon lepturus, Pimelodus grosskopfii, Rotylenchulus reniformis, Scomberomorus niphonius, Sepia esculenta, Terapon jarbua, Teratosphaeria cryptica and Thunnus obesus. These loci were cross-tested on the following species: Austropotamobius italicus, Cambarellus montezumae, Cambarellus puer, Cambarellus shufeldtii, Cambarellus texanus, Chionodraco myersi, Chionodraco rastrospinosus, Coptis chinensis, Coptis chinensis var. brevisepala, Coptis deltoidea, Coptis teeta, Orconectes virilis, Pacifastacus leniusculus, Pimelodus bochii, Procambarus clarkii, Pseudopimelodus bufonius, Rhamdia quelen, Sepia andreana, Sepiella maindroni, Thunnus alalunga, Thunnus albacares, Thunnus maccoyii, Thunnus orientalis, Thunnus thynnus and Thunnus tonggol.