Phylogenetic analyses of the subgenus Mollienesia (Poecilia, Poeciliidae, Teleostei) reveal taxonomic inconsistencies, cryptic biodiversity, and spatio-temporal aspects of diversification in Middle America

Biogeography of Greater Antillean freshwater fishes, with a review of competing hypotheses

In biogeography, vicariance and long-distance dispersal are often characterised as competing scenarios. However, they are related concepts, both relying on collective geological, ecological, and phylogenetic evidence. This is illustrated by freshwater fishes, which may immigrate to islands either when freshwater connections are temporarily present and later severed (vicariance), or by unusual means when ocean gaps are crossed (long-distance dispersal). Marine barriers have a strong filtering effect on freshwater fishes, limiting immigrants to those most capable of oceanic dispersal. The roles of vicariance and dispersal are debated for freshwater fishes of the Greater Antilles. We review three active hypotheses [Cretaceous vicariance, Greater Antilles-Aves Ridge (GAARlandia), long-distance dispersal] and propose long-distance dispersal to be an appropriate model due to limited support for freshwater fish use of landspans. Greater Antillean freshwater fishes have six potential source bioregions (defined from faunal similarity): Northern Gulf of México, Western Gulf of México, Maya Terrane, Chortís Block, Eastern Panamá, and Northern South America. Faunas of the Greater Antilles are composed of taxa immigrating from many of these bioregions, but there is strong compositional disharmony between island and mainland fish faunas (>90% of Antillean species are cyprinodontiforms, compared to <10% in Northern Gulf of México and Northern South America, and ≤50% elsewhere), consistent with a hypothesis of long-distance dispersal. Ancestral-area reconstruction analysis indicates there were 16 or 17 immigration events over the last 51 million years, 14 or 15 of these by cyprinodontiforms. Published divergence estimates and evidence available for each immigration event suggests they occurred at different times and by different pathways, possibly with rafts of vegetation discharged from rivers or washed to sea during storms. If so, ocean currents likely provide critical pathways for immigration when flowing from one landmass to another. On the other hand, currents create dispersal barriers when flowing perpendicularly between landmasses. In addition to high salinity tolerance, cyprinodontiforms collectively display a variety of adaptations that could enhance their ability to live with rafts (small body size, viviparity, low metabolism, amphibiousness, diapause, self-fertilisation). These adaptations likely also helped immigrants establish island populations after arrival and to persist long term thereafter. Cichlids may have used a pseudo bridge (Nicaragua Rise) to reach the Greater Antilles. Gars (Lepisosteidae) may have crossed the Straits of Florida to Cuba, a relatively short crossing that is not a barrier to gene flow for several cyprinodontiform immigrants. Indeed, widespread distributions of Quaternary migrants (Cyprinodon, Gambusia, Kryptolebias), within the Greater Antilles and among neighbouring bioregions, imply that long-distance dispersal is not necessarily inhibitory for well-adapted species, even though it appears to be virtually impossible for all other freshwater fishes.

Nascent transcription reveals regulatory changes in extremophile fishes inhabiting hydrogen sulfide-rich environments

Regulating transcription allows organisms to respond to their environment, both within a single generation (plasticity) and across generations (adaptation). We examined transcriptional differences in gill tissues of fishes in the Poecilia mexicana species complex (family Poeciliidae), which have colonized toxic springs rich in hydrogen sulfide (H2S) in southern Mexico. There are gene expression differences between sulfidic and non-sulfidic populations, yet regulatory mechanisms mediating this gene expression variation remain poorly studied. We combined capped-small RNA sequencing (csRNA-seq), which captures actively transcribed (i.e. nascent) transcripts, and messenger RNA sequencing (mRNA-seq) to examine how variation in transcription, enhancer activity, and associated transcription factor binding sites may facilitate adaptation to extreme environments. csRNA-seq revealed thousands of differentially initiated transcripts between sulfidic and non-sulfidic populations, many of which are involved in H2S detoxification and response. Analyses of transcription factor binding sites in promoter and putative enhancer csRNA-seq peaks identified a suite of transcription factors likely involved in regulating H2S-specific shifts in gene expression, including several key transcription factors known to respond to hypoxia. Our findings uncover a complex interplay of regulatory processes that reflect the divergence of extremophile populations of P. mexicana from their non-sulfidic ancestors and suggest shared responses among evolutionarily independent lineages.

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.

Continental fishes of Nicaragua: diversity, distribution and conservation status; with an annotated and illustrated checklist of species and an identification guide to families

Based on a combination of intensive literature review, electronic database searches, re-identification of museum specimens, and fieldwork, we provide an updated checklist of the continental fishes of Nicaragua. This checklist, systematically arranged at the ordinal and familial level, includes nomenclatural revisions, distributional information, conservation status, and when appropriate, cross-references to previous lists and revisionary works on the fish fauna of the country and region (Central America); illustrations and descriptions to the family level and color photographs of some representative species are also included. According to our results, the native Nicaraguan continental fish fauna is composed of 244 species, divided into 111 genera, 45 families, 24 orders, and two classes; increasing by 60 (i.e., 32.6%) the number of species originally reported by Villa (1982; i.e., the most comprehensive work published to date). The majority of the Nicaraguan native continental fish species, according to their supposed tolerance to salinity, are peripheral (68.9%), followed by secondary freshwater fishes (23.8%), and primary freshwater fishes (only 7.4%). Fish diversity in the 19 major river basins ranged between 21 and 90 species (mean = 58 species, median = 60); the Escondido (Es; 87 species), in the Atlantic, and the Pacfico de Nicaragua (PN; 90), in the Pacific, presented the highest species richness. Most species (188) are restricted to elevations between 0 and 100 meters above sea level (masl); only a few species (18), in contrast, were found at elevations greater than 500 masl (up to 780 masl). Fourteen species in this revision are listed as endemic to Nicaragua. Regarding their conservation status, 22 species (9.0%) are categorized, according to the International Union for Conservation of Nature (IUCN), as Data Deficient, 186 (76.2%) as Least Concern, and seven (2.9%) as Near Threatened; while 24 species (9.8%) are listed as Threatened, thirteen (5.3%) as Vulnerable, two (0.8%) as Endangered, and nine (3.7%) as Critically Endangered. Five species (2.0%) have not yet been evaluated. Six exotic species are also reported. Our results provide a framework for future taxonomic and biogeographic works on fishes from this country and region (Central America), particularly by providing up-to-date knowledge on nomenclature and distributions. In addition, the data provided in this revision will help to monitor fish distributional changes in the future due to human introductions and global change and will aid in conservation decisions concerning Nicaraguan continental fishes.

Sperm specificity and potential paternal effects in gynogenesis in the Amazon Molly (Poecilia formosa)

The Amazon Molly (Poecilia formosa) reproduces by gynogenesis, a relatively rare form of asexual reproduction where sperm is required to trigger embryogenesis, but male genes are not incorporated into the genome of the embryo. Studying gynogenesis could isolate paternal non-genetic effects on reproduction. This study explored which of eleven related species can produce sperm to trigger gynogenesis through natural mating in P. formosa, and whether sympatry affects reproductive success in P. formosa. Reproductive outcomes measured were relative reproductive output (number of offspring in the first brood divided by female standard length), relative embryo output (number of embryos in the first brood divided by female standard length) and combined relative reproductive output (sum of relative reproductive output and relative embryo output). For large (>4 cm) P. formosa, combined relative reproductive output was higher with sympatric Atlantic Molly (Poecilia mexicana) males than with allopatric P. mexicana males. P. formosa produced live offspring or late-stage embryos with all species tested in the genera Poecilia and Limia but did not produce offspring or embryos with males from the genera Gambusia, Girardinus, Heterandria, Poeciliopsis, or Xiphophorus. This information, as well as the limitations characterized in this study, will set a foundation for use of P. formosa as a model for paternal effects and the species specificity of sperm on fertilization, embryogenesis, and reproductive success.

Biogeography of the West Indies: A complex scenario for species radiations in terrestrial and aquatic habitats

Studies of the biogeography of the West Indies are numerous but not all taxonomic groups have received the same attention. Many of the contributions to this field have historically focused on terrestrial vertebrates from a perspective closely linked to the classical theory of island biogeography. However, some recent works have questioned whether some of the assumptions of this theory are too simplistic. In this review, we compiled information about the West Indies biogeography based on an extensive and rigorous literature search. While we offer some background of the main hypotheses that explain the origin of the Caribbean biota, our main purpose here is to highlight divergent diversification patterns observed in terrestrial versus aquatic groups of the West Indian biota and also to shed light on the unbalanced number of studies covering the biogeography of these groups of organisms. We use an objective method to compile existing information in the field and produce a rigorous literature review. Our results show that most of the relevant literature in the field is related to the study of terrestrial organisms (mainly vertebrates) and only a small portion covers aquatic groups. Specifically, livebearing fishes show interesting deviations from the species-area relationship predicted by classical island biogeography theory. We found that species richness on the Greater Antilles is positively correlated with island size but also with the presence of elevations showing that not only island area but also mountainous relief may be an important factor determining the number of freshwater species in the Greater Antilles. Our findings shed light on mechanisms that may differently drive speciation in aquatic versus terrestrial environments suggesting that ecological opportunity could outweigh the importance of island size in speciation. Investigations into freshwater fishes of the West Indies offer a promising avenue for understanding origins and subsequent diversification of the Caribbean biota.

microRNA expression variation as a potential molecular mechanism contributing to adaptation to hydrogen sulphide

microRNAs (miRNAs) are post-transcriptional regulators of gene expression and can play an important role in modulating organismal development and physiology in response to environmental stress. However, the role of miRNAs in mediating adaptation to diverse environments in natural study systems remains largely unexplored. Here, we characterized miRNAs and their expression in Poecilia mexicana, a species of small fish that inhabits both normal streams and extreme environments in the form of springs rich in toxic hydrogen sulphide (H₂ S). We found that P. mexicana has a similar number of miRNA genes as other teleosts. In addition, we identified a large population of mature miRNAs that were differentially expressed between locally adapted populations in contrasting habitats, indicating that miRNAs may contribute to P. mexicana adaptation to sulphidic environments. In silico identification of differentially expressed miRNA-mRNA pairs revealed, in the sulphidic environment, the downregulation of miRNAs predicted to target mRNAs involved in sulphide detoxification and cellular homeostasis, which are pathways essential for life in H₂ S-rich springs. In addition, we found that predicted targets of upregulated miRNAs act in the mitochondria (16.6% of predicted annotated targets), which is the main site of H₂ S toxicity and detoxification, possibly modulating mitochondrial function. Together, the differential regulation of miRNAs between these natural populations suggests that miRNAs may be involved in H₂ S adaptation by promoting functions needed for survival and reducing functions affected by H₂ S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H₂ S. Overall, this study provides a critical stepping-stone towards a comprehensive understanding of the regulatory mechanisms underlying the adaptive variation in gene expression in a natural system.

A new livebearing fish of the genus Limia (Cyprinodontiformes: Poeciliidae) from Lake Miragoane, Haiti

Limia islai, a new species of livebearing fish, is described from Lake Miragoane in south-western Haiti on Hispaniola. The new species has a conspicuous barred pattern consisting of several (4-12) black bars along the body, ray 4p serrae of the gonopodium in males with 10 segments and origin of dorsal fin in females slightly behind the origin of the anal fin. Although the new species colour pattern is similar to that of the humpbacked limia Limia nigrofasciata Regan 1913, L. islai sp. nov. has exclusive morphological features, such as slender body, lack of hump anterior to dorsal fin in males and presence of specific features in the gonopodial suspensory, which allow an unambiguous diagnosis from L. nigrofasciata. L. islai further differs from L. nigrofasciata in reproductive behaviour since L. islai males rely on sneak copulations and gonopodial thrusting, whereas L. nigrofasciata display an elaborate courtship behaviour. The new species is also genetically distinct in both nuclear (Rh, Myh6) and mitochondrial (12S, ND2, D-loop, Cytb) genes from other species in the genus showing reciprocal monophyly. The description of this new Limia species from Lake Miragoane confirms this lake as an important centre of endemism for the genus, with a total of eight endemic species described so far.

Ichthyological Differentiation and Homogenization in the Pánuco Basin, Mexico

Species introductions and extirpations are key aspects of aquatic ecosystem change that need to be examined at large geographic and temporal scales. The Pánuco Basin (Eastern Mexico) has high ichthyological diversity and ecological heterogeneity. However, freshwater fish (FWF) introductions and extirpations since the mid-1900s have modified species range and distribution. We examine changes in FWF species composition in and among four sub-basins of the Pánuco by comparing fish collection records pre-1980 to 2018. Currently, the FWF of the Pánuco includes 95 species. Fishes in the Poeciliidae, Cyprinidae, and Cichlidae, respectively, comprised most records over time. Significant differences in species composition were found between the first (pre-1980) and last (2011–2018) study periods, but not for periods in-between. Eight independent species groups were key for explaining changes in Pánuco river ichthyofauna; one group was dominated by invasive species, and saw increases in the number of records across study periods (faunal homogenization). Another group was formed by species with conservation concern with a declining number of records over time. Thirteen (2 native and 11 non-native) species were responsible for temporal turnover. These results strongly suggest high rates of differentiation over time (via native species loss) following widespread non-native species introductions.

Molecular clocks, biogeography and species diversity in Herichthys with evaluation of the role of Punta del Morro as a vicariant brake along the Mexican Transition Zone in the context of local and global time frame of cichlid diversification

Using molecular dated phylogenies and biogeographic reconstructions, the species diversity, biogeography and time frame of evolution of the genus Herichthys were evaluated. In particular, we test the role of Punta del Morro (PdM) as a vicariant brake along the Mexican Transition Zone in the context of local and global time frame of cichlid diversification using several sets of calibrations. Species diversity in Herichthys is complex and the here employed dating methods suggest young age and rapid divergence for many species while species delimitation methods did not resolve these young species including both sympatric species pairs. Based on our molecular clock dating analyses, Herichthys has colonized its present distribution area significantly prior to the suggested vicariance by PdM (10-17.1 Ma vs. 5 to 7.5 Ma). The PdM constraint is in conflict with all other paleogeographic and fossil constraints including novel ones introduced in this study that are, however, congruent among each other. Our study demonstrates that any cichlid datings significantly older or younger than the bounds presented by our analyses and discussion have to be taken as highly questionable from the point of view of Middle American paleogeography and cichlid biogeography unless we allow the option that cichlid biogeography is completely independent from ecological and geological constraints.

Placing the hybrid origin of the asexual Amazon molly (Poecilia formosa) based on historical climate data

Asexual hybrids are important model organisms for addressing questions in evolution and ecology, especially for understanding the role of hybridization in speciation. They are rare in nature and several hypotheses have been suggested to explain this. We use an asexual fish, the Amazon molly (Poecilia formosa), to establish the area in which it was formed via hybridization 125 000 years ago. Using species distribution models and climate models for the Last Interglacial (LIG) we found that model projections to the LIG show a similar map to the present climate model and parental species potentially overlapped in a relatively small area near Tampico, Mexico. This makes P. formosa one of a few hybrid species for which we know the parental species, the time of hybridization, and likely the place of hybridization. Based on the small area of overlap, our data is in agreement with the idea that asexual hybrids may be rare not because they are evolutionary dead ends but are formed rarely.

Freshwater-to-marine transitions may explain the evolution of herbivory in the subgenus Mollienesia (genus Poecilia, mollies and guppies)

The ability of organisms to cross ecosystem boundaries is an important catalyst of evolutionary diversification. The genus Poecilia (mollies and guppies) is an excellent system for studying ecosystem transitions because species display a range of salinity and dietary preferences, with herbivory concentrated in the subgenus Mollienesia. We reconstructed ancestral habitats and diets across a phylogeny of the genus Poecilia, evaluated diversification rates and used phylogenetically independent contrasts to determine whether diet evolved in response to habitat transition in this group. The results suggest that ancestors of subgenus Mollienesia were exclusively herbivorous, whereas ancestral diets of other Poecilia included animals. We found that transitions across euryhaline boundaries occurred at least once in this group, probably after the divergence of the subgenus Mollienesia. Furthermore, increased salinity affiliation explained 24% of the decrease in animals in the gut, and jaw morphology was associated with the percentage of animals in the gut, but not with the percentage of species occupying saline habitats. These findings suggest that in the genus Poecilia, herbivory evolved in association with transitions from fresh to euryhaline habitats, and jaw morphology evolved in response to the appearance of herbivory. These results provide a rare example of increased diet diversification associated with the transition from freshwater to euryhaline habitats.

Local ancestry analysis reveals genomic convergence in extremophile fishes

The molecular basis of convergent phenotypes is often unknown. However, convergence at a genomic level is predicted when there are large population sizes, gene flow among diverging lineages or strong genetic constraints. We used whole-genome resequencing to investigate genomic convergence in fishes ( Poecilia spp.) that have repeatedly colonized hydrogen sulfide (H₂S)-rich environments in Mexico. We identified genomic similarities in both single nucleotide polymorphisms (SNPs) and structural variants (SVs) among independently derived sulfide spring populations, with approximately 1.2% of the genome being shared among sulfidic ecotypes. We compared these convergent genomic regions to candidate genes for H₂S adaptation identified from transcriptomic analyses and found that a significant proportion of these candidate genes (8%) were also in regions where sulfidic individuals had similar SNPs, while only 1.7% were in regions where sulfidic individuals had similar SVs. Those candidate genes included genes involved in sulfide detoxification, the electron transport chain (the main toxicity target of H₂S) and other processes putatively important for adaptation to sulfidic environments. Regional genomic similarity across independent populations exposed to the same source of selection is consistent with selection on standing variation or introgression of adaptive alleles across divergent lineages. However, combined with previous analyses, our data also support that adaptive changes in mitochondrially encoded subunits arose independently via selection on de novo mutations. Pressing questions remain on what conditions ultimately facilitate the independent rise of adaptive alleles at the same loci in separate populations, and thus, the degree to which evolution is repeatable or predictable. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

The role of Central American barriers in shaping the evolutionary history of the northernmost glassfrog, Hyalinobatrachium fleischmanni (Anura: Centrolenidae)

The complex geological history of Central America has been useful for understanding the processes influencing the distribution and diversity of multiple groups of organisms. Anurans are an excellent choice for such studies because they typically exhibit site fidelity and reduced movement. The objective of this work was to identify the impact of recognized geographic barriers on the genetic structure, phylogeographic patterns and divergence times of a wide-ranging amphibian species, Hyalinobatrachium fleischmanni. We amplified three mitochondrial regions, two coding (COI and ND1) and one ribosomal (16S), in samples collected from the coasts of Veracruz and Guerrero in Mexico to the humid forests of Chocó in Ecuador. We examined the biogeographic history of the species through spatial clustering analyses (Geneland and sPCA), Bayesian and maximum likelihood reconstructions, and spatiotemporal diffusion analysis. Our data suggest a Central American origin of H. fleischmanni and two posterior independent dispersals towards North and South American regions. The first clade comprises individuals from Colombia, Ecuador, Panama and the sister species Hyalinobatrachium tatayoi; this clade shows little structure, despite the presence of the Andes mountain range and the long distances between sampling sites. The second clade consists of individuals from Costa Rica, Nicaragua, and eastern Honduras with no apparent structure. The third clade includes individuals from western Honduras, Guatemala, and Mexico and displays deep population structure. Herein, we synthesize the impact of known geographic areas that act as barriers to glassfrog dispersal and demonstrated their effect of differentiating H. fleischmanni into three markedly isolated clades. The observed genetic structure is associated with an initial dispersal event from Central America followed by vicariance that likely occurred during the Pliocene. The southern samples are characterized by a very recent population expansion, likely related to sea-level and climatic oscillations during the Pleistocene, whereas the structure of the northern clade has probably been driven by dispersal through the Isthmus of Tehuantepec and isolation by the Motagua–Polochic–Jocotán fault system and the Mexican highlands.

Female Choice Undermines the Emergence of Strong Sexual Isolation between Locally Adapted Populations of Atlantic Mollies (Poecilia mexicana)

Divergent selection between ecologically dissimilar habitats promotes local adaptation, which can lead to reproductive isolation (RI). Populations in the Poecilia mexicana species complex have independently adapted to toxic hydrogen sulfide and show varying degrees of RI. Here, we examined the variation in the mate choice component of prezygotic RI. Mate choice tests across drainages (with stimulus males from another drainage) suggest that specific features of the males coupled with a general female preference for yellow color patterns explain the observed variation. Analyses of male body coloration identified the intensity of yellow fin coloration as a strong candidate to explain this pattern, and common-garden rearing suggested heritable population differences. Male sexual ornamentation apparently evolved differently across sulfide-adapted populations, for example because of differences in natural counterselection via predation. The ubiquitous preference for yellow color ornaments in poeciliid females likely undermines the emergence of strong RI, as female discrimination in favor of own males becomes weaker when yellow fin coloration in the respective sulfide ecotype increases. Our study illustrates the complexity of the (partly non-parallel) pathways to divergence among replicated ecological gradients. We suggest that future work should identify the genomic loci involved in the pattern reported here, making use of the increasing genomic and transcriptomic datasets available for our study system.

Molecular evolution and expression of oxygen transport genes in livebearing fishes (Poeciliidae) from hydrogen sulfide rich springs

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S—along with related lineages from non-sulfidic environments—to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H₂ S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H₂ S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H₂ S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H₂ S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H₂ S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

Genome-scale data reveal that endemic Poecilia populations from small sulphidic springs display no evidence of inbreeding

Populations with limited ranges can be highly vulnerable to changes in their environment and are, thus, of high conservation concern. Populations that experience human-induced range reductions are often highly inbred and lack genetic diversity, but it is unknown whether this is also the case for populations with naturally small ranges. The fishes Poecilia sulphuraria (listed as critically endangered) and Poecilia thermalis, which are endemic to small hydrogen sulphide-rich springs in southern Mexico, are examples of such populations with inherently small habitats. We used geometric morphometrics and population genetics to quantify phenotypic and genetic variation within and among two populations of P. sulphuraria and one population of P. thermalis. Principal component analyses revealed phenotypic and genetic differences among the populations. Evidence for inbreeding was low compared to populations that have undergone habitat reduction. The genetic data were also used to infer the demographic history of these populations to obtain estimates for effective population sizes and migration rates. Effective population sizes were large given the small habitats of these populations. Our results imply that these three endemic extremophile populations should each be considered separately for conservation purposes. Additionally, this study suggests that populations in naturally small habitats may have lower rates of inbreeding and higher genetic diversity than expected, and therefore may be better equipped to handle environmental perturbations than anticipated. We caution, however, that the inferred lack of inbreeding and the large effective population sizes could potentially be a result of colonization by genetically diverse ancestors.

The origin and biogeographic diversification of fishes in the family Poeciliidae

The fish subfamily Poeciliinae (sensu Parenti, 1981) is widely distributed across the Western Hemisphere and a dominant component of the fish communities of Central America. Poeciliids have figured prominently in previous studies on the roles of dispersal and vicariance in shaping current geographic distributions. Most recently, Hrbek et al. combined a DNA-based phylogeny of the family with geological models to provide a biogeographic perspective that emphasized the role of both vicariance and dispersal. Here we expand on that effort with a database enlarged in the quantity of sequence represented per species, in the number of species included, and in an enlarged and more balanced representation of the order Cyprinodontiformes. We combine a robust timetree based upon multiple fossil calibrations with enhanced biogeographic analyses that include ancestral area reconstructions to provide a detailed biogeographic history of this clade. Key features of our results are that the family originated in South America, but its major diversification dates to a later colonization of Central America. We also resolve additional colonizations among South, Central and North America and the Caribbean and consider how this reconstruction contributes to our understanding of the mechanisms of dispersal.

Checklist of the continental fishes of the state of Chiapas, Mexico, and their distribution

An updated checklist of the distribution of fishes that inhabit the continental waters of the Mexican state of Chiapas is presented. The state was compartmentalized into 12 hydrological regions for the purpose of understanding the distribution of fish fauna across a state with large physiographic variance. The ichthyofauna of Chiapas is represented by 311 species distributed in two classes, 26 orders, 73 families, and 182 genera, including 12 exotic species. The families with the highest number of species were Cichlidae, Poeciliidae, Sciaenidae, Carangidae, Ariidae, Gobiidae, and Haemulidae. This study attempts to close gaps in knowledge of the distribution of ichthyofauna in the diverse hydrological regions of Chiapas, Mexico.