Search

The Turquoise Killifish is an important vertebrate for the study of aging and age-related diseases due to its short lifespan. Within Nothobranchiidae, species possess annual, semi-annual, or non-annual life-histories. We took a comparative approach and examined gene expression profiles (QuantSeq) from 62 individuals from eleven nothobranchid species that span three life-histories. Our results show significant differences in differentially expressed genes (DEGs) across life-histories with non-annuals and semi-annuals being most similar, and annuals being the most distinct. At finer scales, we recovered significant differences in DEGs for DNA repair genes and show that non-annual and semi-annuals share similar gene expression profiles, while annuals are distinct. Most of the GO terms enriched in annuals are related to metabolic processes. However, GO terms, including translation, protein transport, and DNA replication initiation also are enriched in annuals. Non-annuals are enriched in Notch signaling pathway genes and downregulated in the canonical Wnt signaling pathway compared to annual species, which suggests that non-annuals have stronger regulation in cellular processes. This study provides support for congruency in DEGs involved in these life-histories and provides strong evidence that a particular set of candidate genes may be worthy of study to investigate their role in the aging process., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Leow, Piller. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Repeated evolution of novel life histories that are correlated with ecological variables offers opportunities to study convergence in genetic, developmental, and metabolic features. Nearly half of the 800 species of Aplocheiloid killifishes, a clade of teleost fishes with a circumtropical distribution, are "annual" or seasonal species that survive in ephemeral bodies of water that desiccate and are unfeasible for growth, reproduction, or survival for weeks to months every year. But the repeated evolution of adaptations that are key features of the annual life history among these fishes remains poorly known without a robust phylogenetic framework. We present a large-scale phylogenomic reconstruction of aplocheiloid killifishes evolution using newly sequenced transcriptomes obtained from a diversity of killifish lineages representing putative independent origins of annualism. Ancestral state estimation shows that developmental dormancy (diapause), a key trait of the killifish annual life cycle, may have originated up to seven times independently among African and South American lineages. To further explore the genetic basis of this unique trait, we measure changes in evolutionary rates among orthologous genes across the killifish tree of life by quantifying codon evolution using dN/dS ratios. We show that some genes have higher dN/dS ratios in lineages leading to species with annual life history. Many of them constitute key developmental genes or nuclear-encoded metabolic genes that control oxidative phosphorylation. Lastly, we compare these genes with higher ω to genes previously associated to developmental dormancy and metabolic shifts in killifishes and other vertebrates, and thereby identify molecular evolutionary signatures of repeated transitions to extreme environments., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Three nominal species of the killifish genus Aplocheilus are reported from the lowlands of Sri Lanka. Two of these, Aplocheilus dayi and Aplocheilus werneri, are considered endemic to the island, whereas Aplocheilus parvus is reported from both Sri Lanka and Peninsular India. Here, based on a collection from 28 locations in Sri Lanka, also including a dataset of Asian Aplocheilus downloaded from GenBank, we present a phylogeny constructed from the mitochondrial cytochrome b (cytb), mitochondrial cytochrome c oxidase subunit 1 (cox1), and nuclear recombination activating protein 1 (rag1), and investigate the interrelationships of the species of Aplocheilus in Sri Lanka. The endemic Sri Lankan aplocheilid clade comprising A. dayi and A. werneri is recovered as the sister group to the clade comprising A. parvus from Sri Lanka and Aplocheilus blockii from Peninsular India. The reciprocal monophyly of A. dayi and A. werneri is not supported in our molecular phylogeny. A. dayi and A. werneri display strong sexual dimorphism, but species-level differences are subtle, explained mostly by pigmentation patterns. Their phenotypes exhibit a parapatric distribution and may represent locally adapted forms of a single species. Alternatively, the present study does not rule out the possibility that A. dayi and A. werneri may represent an incipient species pair or that they have undergone introgression or hybridization in their contact zones. We provide evidence that the Nilwala-Gin region of southwestern Sri Lanka may have acted as a drought refugium for these fishes., (© 2024 Fisheries Society of the British Isles.)

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species., Competing Interests: Declaration of interests A.B. is a scientific advisory board member of Calico. M.P.S. is a co-founder and the scientific advisory board member of Personalis, Qbio, January AI, SensOmics, Filtricine, Protos, Mirvie, Onza, Marble Therapeutics, Iollo, and NextThought AI. He is also on the scientific advisory board of Jupiter, Applied Cognition, Neuvivo, Mitrix, and Enovone. K.C. is currently an AstraZeneca employee., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish ( Cyprinodon tularosa ) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global F ST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose., Competing Interests: Competing interests statement:The authors declare no competing interest.

Kryptolebias marmoratus (Kmar), a teleost fish of the order Cyprinodontiformes, has a suite of unique phenotypes and behaviors not observed in other fishes. Many of these phenotypes are discrete and highly plastic-varying over time within an individual, and in some cases reversible. Kmar and its interfertile sister species, K. hermaphroditus, are the only known self-fertile vertebrates. This unusual sexual mode has the potential to provide unique insights into the regulation of vertebrate sexual development, and also lends itself to genetics. Kmar is easily adapted to the lab and requires little maintenance. However, its internal fertilization and small clutch size limits its experimental use. To support Kmar as a genetic model, we compared alternative husbandry techniques to maximize recovery of early cleavage-stage embryos. We find that frequent egg collection enhances yield, and that protease treatment promotes the greatest hatching success. We completed a forward mutagenesis screen and recovered several mutant lines that serve as important tools for genetics in this model. Several will serve as useful viable recessive markers for marking crosses. Importantly, the mutant kissylips lays embryos at twice the rate of wild-type. Combining frequent egg collection with the kissylips mutant background allows for a substantial enhancement of early embryo yield. These improvements were sufficient to allow experimental analysis of early development and the successful mono- and bi-allelic targeted knockout of an endogenous tyrosinase gene with CRISPR/Cas9 nucleases. Collectively, these tools will facilitate modern developmental genetics in this fascinating fish, leading to future insights into the regulation of plasticity., (© 2023 The Authors. Journal of Experimental Zoology Part B: Molecular and Developmental Evolution published by Wiley Periodicals LLC.)

The Deepwater Horizon disaster of April 2010 was the largest oil spill in U.S. history and exerted catastrophic effects on several ecologically important fish species in the Gulf of Mexico (GoM). Within fish, the microbiome plays a key symbiotic role in maintaining host health and aids in acquiring nutrients, supporting immune function, and modulating behavior. The aim of this study was to examine if exposure to weathered oil might produce significant shifts in fish gut-associated microbial communities as determined from taxa and genes known for hydrocarbon degradation, and whether foraging behavior was affected. The gut microbiome (16S rRNA and shotgun metagenomics) of sheepshead minnow ( Cyprinodon variegatus ) was characterized after fish were exposed to oil in High Energy Water Accommodated Fractions (HEWAF; tPAH = 81.1 ± 12.4 µg/L) for 7 days. A foraging behavioral assay was used to determine feeding efficiency before and after oil exposure. The fish gut microbiome was not significantly altered in alpha or beta diversity. None of the most abundant taxa produced any significant shifts as a result of oil exposure, with only rare taxa showing significant shifts in abundance between treatments. However, several bioindicator taxa known for hydrocarbon degradation were detected in the oil treatment, primarily Sphingomonas and Acinetobacter . Notably, the genus Stenotrophomonas was detected in high abundance in 16S data, which previously was not described as a core member of fish gut microbiomes. Data also demonstrated that behavior was not significantly affected by oil exposure. Potential low bioavailability of the oil may have been a factor in our observation of minor shifts in taxa and no behavioral effects. This study lays a foundation for understanding the microbiome of captive sheepshead minnows and indicates the need for further research to elucidate the responses of the fish gut-microbiome under oil spill conditions.

Hypsolebia antenori is a seasonal killifish considered to be broadly distributed along the Caatingas basins of northeastern Brazil, with records in the coastal drainages of the rio Pacoti, rio Jaguaribe, crrego Virglio and rio Apodi-Mossor basins. Based on morphological characters and molecular phylogenetic reconstructions, it was possible to diagnose two new species and restrict H. antenori to the rio Jaguaribe basin. Hypsolebias gongobira new species is described from a temporary pool in the rio Pacoti basin in Cear State, in syntopy with the species threatened H. longignatus. Hypsolebias bonita new species occurs in the floodplains of the rio Apodi-Mossor basin and in the crrego Virglio microbasin in the Furna Feia National Park, in Rio Grande do Norte State. We also discuss the conservation status of these new species and how the So Francisco Interbasin Water transfer, and the agricultural and urban expansion in the semiarid may be adversely affecting these seasonal killifishes.

Satellite DNA (satDNA) is a rapidly evolving class of tandem repeats, with some monomers being involved in centromere organization and function. To identify repeats associated with (peri)centromeric regions, we investigated satDNA across Southern and Coastal clades of African annual killifishes of the genus Nothobranchius. Molecular cytogenetic and bioinformatic analyses revealed that two previously identified satellites, designated here as NkadSat01-77 and NfurSat01-348, are associated with (peri)centromeres only in one lineage of the Southern clade. NfurSat01-348 was, however, additionally detected outside centromeres in three members of the Coastal clade. We also identified a novel satDNA, NrubSat01-48, associated with (peri)centromeres in N. foerschi, N. guentheri, and N. rubripinnis. Our findings revealed fast turnover of satDNA associated with (peri)centromeres and different trends in their evolution in two clades of the genus Nothobranchius., (© 2023. The Author(s).)

Understanding the genetic basis of novel adaptations in new species is a fundamental question in biology. Here we demonstrate a new role for galr2 in vertebrate craniofacial development using an adaptive radiation of trophic specialist pupfishes endemic to San Salvador Island, Bahamas. We confirmed the loss of a putative Sry transcription factor binding site upstream of galr2 in scale-eating pupfish and found significant spatial differences in galr2 expression among pupfish species in Meckel's cartilage using in situ hybridization chain reaction (HCR). We then experimentally demonstrated a novel role for Galr2 in craniofacial development by exposing embryos to Garl2-inhibiting drugs. Galr2-inhibition reduced Meckel's cartilage length and increased chondrocyte density in both trophic specialists but not in the generalist genetic background. We propose a mechanism for jaw elongation in scale-eaters based on the reduced expression of galr2 due to the loss of a putative Sry binding site. Fewer Galr2 receptors in the scale-eater Meckel's cartilage may result in their enlarged jaw lengths as adults by limiting opportunities for a circulating Galr2 agonist to bind to these receptors during development. Our findings illustrate the growing utility of linking candidate adaptive SNPs in non-model systems with highly divergent phenotypes to novel vertebrate gene functions.

Methylmercury (MeHg) is a ubiquitous bioaccumulative neurotoxicant present in aquatic ecosystems. It is known to alter behaviors, sensory functions and learning abilities in fish and other vertebrates. Developmental and early-life stages exposure to MeHg can lead to brain damage with immediate consequences on larvae behavior, but may also induce long term effects in adults after a detoxification period. However, very little is known about developmental origin of behavioral impairment in adults due to early exposure to MeHg. The aim of this study is to assess whether early-life MeHg exposure induces immediate and/or delayed effects on behaviors, related genes expression and DNA methylation (one of epigenetic mechanisms). To reach this goal, newly hatched larvae of mangrove rivulus fish, Kryptolebias marmoratus, were exposed to two sub-lethal concentrations of MeHg (90 μg/L and 135 µg/L) for 7 days, and immediate and delayed effects were assessed respectively in 7 dph (days post-hatching) and 90 dph fish. This species naturally produces isogenic lineages due to its self-fertilizing reproduction system, which is unique among vertebrates. It allows to study how environment stressors can influence organism's phenotype while minimizing genetic variability. As results, both MeHg exposures are associated with a decreased foraging efficiency and thigmotaxis, and a dose-dependent reduction in larvae locomotor activity. Regarding molecular analysis in larvae whole bodies, both MeHg exposures induced significant decreased expression of DNMT3a, MAOA, MeCP2 and NIPBL, and significant increase of GSS, but none of those genes underwent methylation changes in targeted CpGs. None of significant behavioral and molecular impairments observed in 7-dph larvae were found in 90-dph adults, which highlight a distinction between immediate and delayed effects of developmental MeHg exposure. Our results suggest implications of aminergic system and its neurotransmitters, redox/methylation trade-off and possibly other epigenetic mechanisms in MeHg neurotoxicity underlying behavioral alterations in rivulus., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper, (Copyright © 2023 Elsevier B.V. All rights reserved.)

Background: The Spanish toothcarp (Aphanius iberus Valenciennes, 1846) is a small fish endemic to the eastern coastline of the Iberian Peninsula and is currently listed as "Endangered" (category IUCN: EN). It mainly inhabits brackish waters which can exhibit large fluctuations in temperature and salinity throughout the year. The genetics of A. iberus are not well-known since most studies have only evaluated the genetic structure of the species under a conservation framework in order to identify its potential conservation units. Different phylogenetic relationships of Aphanius have been published based on some particular genes. In the present study, the entire mitochondrial genome of A. iberus was obtained for the first time in the context of an A. iberus reference genome and a hypothesis regarding its phylogenetic position was considered., Methods and Results: The mitogenome (a circular doble-stranded DNA sequence of 16,708 bp) was reconstructed and aligned against 83 Cyprinodontiformes and two outgroup taxa to identify the phylogenetic position of A. iberus. PartitionFinder was first used to test for the best evolutionary model and the phylogenetic analyses were performed using two methods: Maximun-Likelihood Approximation (IQ-Tree) and Bayesian inference (MrBayes). Our results show that A. iberus forms a sister group with Orestias ascotanensis, a cyprinodontiform species native to South America., Conclusions: The results were congruent with the traditional morphometric reconstructed trees and with a geological vicariant hypothesis involving Cyprinodontiformes where Aphaniidae is shown as a monophyletic family separated from the family Cyprinodontidae. The information gathered from this study is not only valuable for improving our understanding of the evolutionary history of A. iberus, but for future genomic studies involving the species., (© 2023. The Author(s).)

When species hybridize, one F1 hybrid cross type often predominates. Such asymmetry can arise from differences in a variety of reproductive barriers, but the relative roles and concordance of pre-mating, post-mating prezygotic, and post-zygotic barriers in producing these biases in natural animal populations have not been widely investigated. Here, we study a population of predominantly F1 hybrids between two killifish species (Fundulus heteroclitus and F. diaphanus) in which >95% of F1 hybrids have F. diaphanus mothers and F. heteroclitus fathers (D♀ × H♂). To determine why F. heteroclitus × F. diaphanus F1 hybrids (H♀ × D♂) are so rare, we tested for asymmetry in pre-mating reproductive barriers (female preference and male aggression) at a common salinity (10 ppt) and post-mating, pre-zygotic (fertilization success) and post-zygotic (embryonic development time and hatching success) reproductive barriers at a range of ecologically relevant salinities (0, 5, 10, and 15 ppt). We found that F. heteroclitus females preferred conspecific males, whereas F. diaphanus females did not, matching the observed cross bias in the wild. Naturally rare H♀ × D♂ crosses also had lower fertilization success than all other cross types, and a lower hatching success than the prevalent D♀ × H♂ crosses at the salinity found in the hybrid zone centre (10 ppt). Furthermore, the naturally predominant D♀ × H♂ crosses had a higher hatching success than F. diaphanus crosses at 10 ppt, which may further increase their relative abundance. The present study suggests that a combination of incomplete mating, post-mating pre-zygotic and post-zygotic reproductive isolating mechanisms act in concert to produce hybrid asymmetry in this system., (© 2023 European Society for Evolutionary Biology.)

Homomorphic sex chromosomes and their turnover are common in teleosts. We investigated the evolution of nascent sex chromosomes in several populations of two sister species of African annual killifishes, Nothobranchius furzeri and N. kadleci, focusing on their under-studied repetitive landscape. We combined bioinformatic analyses of the repeatome with molecular cytogenetic techniques, including comparative genomic hybridization, fluorescence in situ hybridization with satellite sequences, ribosomal RNA genes (rDNA) and bacterial artificial chromosomes (BACs), and immunostaining of SYCP3 and MLH1 proteins to mark lateral elements of synaptonemal complexes and recombination sites, respectively. Both species share the same heteromorphic XY sex chromosome system, which thus evolved prior to their divergence. This was corroborated by sequence analysis of a putative master sex determining (MSD) gene gdf6Y in both species. Based on their divergence, differentiation of the XY sex chromosome pair started approximately 2 million years ago. In all populations, the gdf6Y gene mapped within a region rich in satellite DNA on the Y chromosome long arms. Despite their heteromorphism, X and Y chromosomes mostly pair regularly in meiosis, implying synaptic adjustment. In N. kadleci, Y-linked paracentric inversions like those previously reported in N. furzeri were detected. An inversion involving the MSD gene may suppress occasional recombination in the region, which we otherwise evidenced in the N. furzeri population MZCS-121 of the Limpopo clade lacking this inversion. Y chromosome centromeric repeats were reduced compared with the X chromosome and autosomes, which points to a role of relaxed meiotic drive in shaping the Y chromosome repeat landscape. We speculate that the recombination rate between sex chromosomes was reduced due to heterochiasmy. The observed differences between the repeat accumulations on the X and Y chromosomes probably result from high repeat turnover and may not relate closely to the divergence inferred from earlier SNP analyses., (© 2022. The Author(s), under exclusive licence to Springer Nature B.V.)

Introgression is a widespread evolutionary process leading to phylogenetic inconsistencies among distinct parts of the genomes, particularly between mitochondrial and nuclear-based phylogenetic reconstructions (e.g., mito-nuclear discordances). Here, we used mtDNA and genome-wide nuclear sites to provide the first phylogenomic-based hypothesis on the evolutionary relationships within the killifish genus Kryptolebias. In addition, we tested for evidence of past introgression in the genus given the multiple reports of undergoing hybridization between its members. Our mtDNA phylogeny generally agreed with the relationships previously proposed for the genus. However, our reconstruction based on nuclear DNA revealed an unknown lineage - Kryptolebias sp. 'ESP' - as the sister group of the self-fertilizing mangrove killifishes, K. marmoratus and K. hermaphroditus. All individuals sequenced of Kryptolebias sp. 'ESP' had the same mtDNA haplotype commonly observed in K. hermaphroditus, demonstrating a clear case of mito-nuclear discordance. Our analysis further confirmed extensive history of introgression between Kryptolebias sp. 'ESP' and K. hermaphroditus. Population genomics analyses indicate no current gene flow between the two lineages, despite their current sympatry and history of introgression. We also confirmed introgression between other species pairs in the genus that have been recently reported to form hybrid zones. Overall, our study provides a phylogenomic reconstruction covering most of the Kryptolebias species, reveals a new lineage hidden in a case of mito-nuclear discordance, and provides evidence of multiple events of ancestral introgression in the genus. These findings underscore the importance of investigating different genomic information in a phylogenetic framework, particularly in taxa where introgression is common as in the sexually diverse mangrove killifishes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Adaptive radiations offer an excellent opportunity to understand the eco-evolutionary dynamics of gut microbiota and host niche specialization. In a laboratory common garden, we compared the gut microbiota of two novel derived trophic specialist pupfishes, a scale-eater and a molluscivore, to closely related and distant outgroup generalist populations, spanning both rapid trophic evolution within 10 kya and stable generalist diets persisting over 11 Mya. We predicted an adaptive and highly divergent microbiome composition in the trophic specialists reflecting their rapid rates of craniofacial and behavioral diversification. We sequenced 16S rRNA amplicons of gut microbiomes from lab-reared adult pupfishes raised under identical conditions and fed the same high protein diet. In contrast to our predictions, gut microbiota largely reflected phylogenetic distance among species, rather than generalist or specialist life history, in support of phylosymbiosis. However, we did find significant enrichment of Burkholderiaceae bacteria in replicated lab-reared scale-eater populations. These bacteria sometimes digest collagen, the major component of fish scales, supporting an adaptive shift. We also found some enrichment of Rhodobacteraceae and Planctomycetia in lab-reared molluscivore populations, but these bacteria target cellulose. Overall phylogenetic conservation of microbiome composition contrasts with predictions of adaptive radiation theory and observations of rapid diversification in all other trophic traits in these hosts, including craniofacial morphology, foraging behavior, aggression, and gene expression, suggesting that the functional role of these minor shifts in microbiota will be important for understanding the role of the microbiome in trophic diversification., Competing Interests: The authors have declared that no competing interests exist.

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., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scale-eater) endemic to San Salvador Island, Bahamas. We leveraged whole-genome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 single nucleotide polymorphisms (SNPs) that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene ( mettl21e ) was also misexpressed in lab-reared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scale-eating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation., Competing Interests: AP, ER, KG, LB, CM No competing interests declared, (© 2022, Patton et al.)

Adaptive radiations involve astounding bursts of phenotypic, ecological and species diversity. However, the microevolutionary processes that underlie the origins of these bursts are still poorly understood. We report the discovery of an intermediate C. sp. 'wide-mouth' scale-eating ecomorph in a sympatric radiation of Cyprinodon pupfishes, illuminating the transition from a widespread algae-eating generalist to a novel microendemic scale-eating specialist. We first show that this ecomorph occurs in sympatry with generalist C. variegatus and scale-eating specialist C. desquamator on San Salvador Island, Bahamas, but is genetically differentiated, morphologically distinct and often consumes scales. We then compared the timing of selective sweeps on shared and unique adaptive variants in trophic specialists to characterize their adaptive walk. Shared adaptive regions swept first in both the specialist desquamator and the intermediate 'wide-mouth' ecomorph, followed by unique sweeps of introgressed variation in 'wide-mouth' and de novo variation in desquamator . The two scale-eating populations additionally shared 9% of their hard selective sweeps with the molluscivore C. brontotheroides , despite no single common ancestor among specialists. Our work provides a new microevolutionary framework for investigating how major ecological transitions occur and illustrates how both shared and unique genetic variation can provide a bridge for multiple species to access novel ecological niches.

Orestias ascotanensis (Cyprinodontidae) is a teleost pupfish endemic to springs feeding into the Ascotan saltpan in the Chilean Altiplano (3,700 m.a.s.l.) and represents an opportunity to study adaptations to high-altitude aquatic environments. We have de novo assembled the genome of O. ascotanensis at high coverage. Comparative analysis of the O. ascotanensis genome showed an overall process of contraction, including loss of genes related to G-protein signaling, chemotaxis and signal transduction, while there was expansion of gene families associated with microtubule-based movement and protein ubiquitination. We identified 818 genes under positive selection, many of which are involved in DNA repair. Additionally, we identified novel and conserved microRNAs expressed in O. ascotanensis and its closely-related species, Orestias gloriae. Our analysis suggests that positive selection and expansion of genes that preserve genome stability are a potential adaptive mechanism to cope with the increased solar UV radiation to which high-altitude animals are exposed to., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

The White Sands pupfish (Cyprinodon tularosa), endemic to New Mexico in Southwestern North America, is of conservation concern due in part to invasive species, chemical pollution, and groundwater withdrawal. Herein, we developed a draft reference genome and use it to provide biological insights into the evolution and conservation of C. tularosa. We used our assembly to localize microsatellite markers previously used to demarcate evolutionary significant units (ESU), quantified genomic divergence and transposable element profiles between species, and compared C. tularosa genomic diversity related species. Our de novo assembly of PacBio Sequel II error-corrected reads resulted in a 1.08 Gb draft genome with a contig N50 of 1.4 Mb and 25,260 annotated protein coding genes, including 95% of the expected Actinopterygii conserved complete single-copy orthologues. Many of the C. tularosa microsatellite markers used for conservation assessments fell within, or near, genes and exhibited a pattern of increased heterozygosity near genic areas compared to those in intergenic regions. Nuclear alignments between these two species revealed 193 genes contained in rapidly diverging tracts; transposable element profiles were largely concordant and suggest a shared, rapid expansion of LINE and Gypsy elements. Genome-wide heterozygosity was markedly lower in C. tularosa compared to estimates from other related species, probably because of smaller long-term effective population sizes constrained by their isolated and limited habitat. Overall, these inferences provide new insights into C. tularosa that should help inform future management efforts., (© 2021 John Wiley & Sons Ltd.)

Increased nutrient loading has led to eutrophication of coastal shelf waters which has resulted in increased prevalence of persistent hypoxic zones - areas in which the dissolved oxygen content of the water drops below 2 mg/L. The northern Gulf of Mexico, fed primarily by the Mississippi River watershed, undergoes annual establishment of one of the largest hypoxic zones in the world. Exposure to hypoxia can induce physiological impacts in fish cardiac systems that include bradycardia, changes in stroke volume, and altered cardiovascular vessel development. While these impacts have been addressed at the functional level, there is little information regarding the molecular basis for these changes. This study used transcriptomic analysis techniques to interrogate the effects of hypoxia exposure on the developing cardiovascular system in newly hatched larvae of two estuarine species that occupy the same ecological niche - the sheepshead minnow (Cyprinodon variegatus) and the Gulf killifish (Fundulus grandis). Results suggest that while differential gene expression is largely distinct between the two species, downstream impacts on pathways and functional responses such as reduced cardiac hypertrophy, modulation of blood pressure, and increased incidence of apoptosis appear to be conserved. Further, differences in the magnitude of these conserved responses may suggest that the length of embryonic development could impart a level of resiliency to hypoxic perturbation in early life stage fish., (Copyright © 2021 Elsevier Inc. All rights reserved.)

To investigate the origins and stages of vertebrate adaptive radiation, we reconstructed the spatial and temporal histories of adaptive alleles underlying major phenotypic axes of diversification from the genomes of 202 Caribbean pupfishes. On a single Bahamian island, ancient standing variation from disjunct geographic sources was reassembled into new combinations under strong directional selection for adaptation to the novel trophic niches of scale-eating and molluscivory. We found evidence for two longstanding hypotheses of adaptive radiation: hybrid swarm origins and temporal stages of adaptation. Using a combination of population genomics, transcriptomics, and genome-wide association mapping, we demonstrate that this microendemic adaptive radiation of novel trophic specialists on San Salvador Island, Bahamas experienced twice as much adaptive introgression as generalist populations on neighboring islands and that adaptive divergence occurred in stages. First, standing regulatory variation in genes associated with feeding behavior ( prlh , cfap20 , and rmi1 ) were swept to fixation by selection, then standing regulatory variation in genes associated with craniofacial and muscular development ( itga5 , ext1 , cyp26b1 , and galr2 ) and finally the only de novo nonsynonymous substitution in an osteogenic transcription factor and oncogene ( twist1 ) swept to fixation most recently. Our results demonstrate how ancient alleles maintained in distinct environmental refugia can be assembled into new adaptive combinations and provide a framework for reconstructing the spatiotemporal landscape of adaptation and speciation., Competing Interests: The authors declare no competing interest.

This paper reports a phylogeny of the African killifishes (Genus Nothobranchius, Order Cyprinodontiformes) informed by five genetic markers (three nuclear, two mitochondrial) of 80 taxa (seven undescribed and 73 of the 92 recognized species). These short-lived annual fishes occupy seasonally wet habitats in central and eastern Africa, and their distribution coincides largely with the East African Rift System (EARS). The fossil dates of sister clades used to constrain a chronometric tree of all sampled Nothobranchius recovered the origin of the genus at ~13.27 Mya. It was followed by the radiations of six principal clades through the Neogene. An ancestral area estimation tested competing biogeographical hypotheses to constrain the ancestral origin of the genus to the Nilo-Sudan Ecoregion, which seeded a mid-Miocene dispersal event into the Coastal ecoregion, followed closely (~10 Mya) by dispersals southward across the Mozambique coastal plain into the Limpopo Ecoregion. Extending westwards across the Tanzanian plateau, a pulse of radiations through the Pliocene were associated with dispersals and fragmentation of wetlands across the Kalahari and Uganda Ecoregions. We interpret this congruence of drainage rearrangements with dispersals and cladogenic events of Nothobranchius to reflect congruent responses to recurrent uplift and rifting. The coevolution of these freshwater fishes and wetlands is attributed to ultimate control by tectonics, as the EARS extended southwards during the Neogene. Geobiological consilience of the combined evidence supports a tectonic hypothesis for the evolution of Nothobranchius., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

The Austrolebias adloffi species group encompasses a diverse lineage of annual killifishes that occurs along the Laguna dos Patos/Lagoa Mirim system, in both Brazilian and Uruguayan territories. We herein employ an integrative taxonomy approach to describe two new species of the group, inferring their phylogenetic relationships and evaluating their conservation status. Austrolebias cheffei sp. nov. and Austrolebias lourenciano sp. nov. are herein described from the western portion of the Laguna dos Patos system. Austrolebias cheffei is distinguished from the remaining species of the A. adloffi species group by presenting a yellowish green or yellowish blue dorsal fin, with wide black to dark brown bars extending from the base to the middle portion of the dorsal and anal fins in the males. Austrolebias lourenciano is distinguished from the remaining species of the A. adloffi species group by presenting a yellowish green dorsal fin, with light yellow or light bluish bars forming small triangles, interspersed with small dark brown rows of blotches in the dorsal fin base, and greenish blue anal fin, sometimes with lighter elongated yellowish iridescent blotches, limited to the basal region. According to mitochondrial cytb sequences, both species are reciprocally monophyletic relative to other species of the A. adloffi species group, and present positive barcoding gap values. Interestingly, both new species form a grade that is closely related to Austrolebias aff. minuano 1, an undescribed species that occurs at the opposite margin of the Laguna dos Patos. Among the other evaluated species, A. bagual, A. aff. minuano 1, A. nigrofasciatus, A. pelotapes, A. pongondo, A. arachan, and A. viarius also revealed to be reciprocally monophyletic, whereas A. minuano and A. adloffi revealed to be paraphyletic in regard to A. charrua and A. aff. minuano 2, respectively, and A. nachtigalli is subdivided in two clades, one of which including A. reicherti, which points to the need of a taxonomic review of the group. In addition, we discussed the conservation status of the new species, corrected the type locality of A. pongondo, and provided a dichotomous identification key of the A. adloffi species group.

Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists-a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species-only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The Deepwater Horizon oil spill released 3.19 million barrels of crude oil into the Gulf of Mexico, making it the largest oil spill in U.S. history. Weathering and the application of dispersants can alter the solubility of compounds within crude oil, thus modifying the acute toxicity of the crude oil to aquatic life. The primary aim of our study was to determine the lasting impact of early-life stage sheepshead minnow (Cyprinodon variegatus variegatus) exposure to weathered, unweathered and dispersed crude oil on prey capture, male aggression, novel object interaction and global DNA methylation. Embryos were exposed from 1 to 10 dpf to water accommodations of crude oil and were raised to adulthood in artificial seawater. Our results suggest exposure to crude oil did not result in lasting impairment of complex behavioral responses of male sheepshead minnow. Exposure to dispersed weathered oil, however, decreased border dwelling in response to a novel object (i.e. decreased anxiety). Principal component analysis revealed that exposure to weathered oil had no overarching effect, but that unweathered crude oil increased variability in exploratory behaviors but decreased variability in anxiety-associated behaviors. Further work is needed to understand the effects of oil exposure on fish behavior and the potential ecological impact of subtle behavioral changes in fishes., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Vertebrates vary in their ability to regenerate, and the genetic mechanisms underlying such disparity remain elusive. Comparative epigenomic profiling and single-cell sequencing of two related teleost fish uncovered species-specific and evolutionarily conserved genomic responses to regeneration. The conserved response revealed several regeneration-responsive enhancers (RREs), including an element upstream to inhibin beta A ( inhba ), a known effector of vertebrate regeneration. This element activated expression in regenerating transgenic fish, and its genomic deletion perturbed caudal fin regeneration and abrogated cardiac regeneration altogether. The enhancer is present in mammals, shares functionally essential activator protein 1 (AP-1)-binding motifs, and responds to injury, but it cannot rescue regeneration in fish. This work suggests that changes in AP-1-enriched RREs are likely a crucial source of loss of regenerative capacities in vertebrates., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species., (© 2020 John Wiley & Sons Ltd.)

Understanding the genetic basis for phenotypic differences is fundamental to the study of macroevolutionary patterns of biological diversity. While technological advances in DNA sequencing have made researching genetic variation in wild taxa routine, fully understanding how these variants affect phenotype requires taking the next step to investigate how genetic changes alter cell and tissue interactions that ultimately produce phenotypes. In this article, we investigate a role for cell proliferation as a developmental source of craniofacial diversity in a radiation of 3 species of Cyprinodon from San Salvador Island, Bahamas. Patterns of cell proliferation in the heads of hatching-age fish differ among species of Cyprinodon, and correlate with differences in allometric growth rate among the jaws of 3 distinct species. Regional patterns of cell proliferation in the head are complex, resulting in an unintuitive result in which lower levels of cell proliferation in the posterior head region are associated with the development of relatively larger jaws in one species. We combine these data with previously published morphological and genomic data to show how studying the mechanisms generating phenotype at the cellular and tissue levels of biological organization can help mechanistically link genomic studies with classic morphological studies., (© The American Genetic Association 2019. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Explaining how populations adapt to environments is among the foremost objectives of evolutionary theory. Over generations, natural selection impels the phenotypic distribution of a population based on individual variation in phenotype and fitness. However, environmental conditions can also shape how individuals develop within their lifetime to influence which phenotypes are expressed in a population. It has been proposed that such environmentally-initiated phenotypic variation - also termed developmental plasticity - may enable adaptive evolution under some scenarios. As dynamic regulators of development and phenotypic expression, hormones are important physiological mediators of developmental plasticity. Patterns of hormone secretion, hormone transport, and the sensitivity of tissues to hormones can each be altered by environmental conditions, and understanding how endocrine regulation shapes phenotypic development in an ecologically-relevant context has much to contribute toward clarifying the role of plasticity in evolutionary adaptation. This article explores how the environmental sensitivity of endocrine regulation may facilitate 'plasticity-first' evolution by generating phenotypic variants that precede adaptation to altered or novel environments. Predictions arising from 'plasticity-first' evolution are examined in the context of thyroid hormone mediation of morphological plasticity in Cyprinodon pupfishes from the Death Valley region of California and Nevada, USA. This clade of extremophile fishes diversified morphologically over the last ~20,000 years, and observations that some populations experienced contemporary phenotypic differentiation under recent habitat change provide evidence that hormone-mediate plasticity preceded genetic assimilation of morphology in one of the region's species: the Devils Hole pupfish, Cyprinodon diabolis. This example illustrates how conceptualizing hormones not only as regulators of homeostasis, but also as developmental intermediaries between environment conditions and phenotypic variation at the individual-, population-, and species-levels can enrich our understanding of endocrine regulation both as a facilitator of phenotypic change under shifting environments, and as important proximate mechanisms that may initiate 'plasticity-first' evolutionary adaptation., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Background: African annual killifishes (Nothobranchius spp.) are adapted to seasonally desiccating habitats (ephemeral pools), surviving dry periods as dormant eggs. Given their peculiar life history, geographic aspects of their diversity uniquely combine patterns typical for freshwater taxa (river basin structure and elevation gradient) and terrestrial animals (rivers acting as major dispersal barriers). However, our current knowledge on fine-scale inter-specific and intra-specific genetic diversity of African annual fish is limited to a single, particularly dry region of their distribution (subtropical Mozambique). Using a widespread annual killifish from coastal Tanzania and Kenya, we tested whether the same pattern of genetic divergence pertains to a wet equatorial region in the centre of Nothobranchius distribution., Results: In populations of Nothobranchius melanospilus species group across its range, we genotyped a part of mitochondrial cytochrome oxidase subunit 1 (COI) gene (83 individuals from 22 populations) and 10 nuclear microsatellite markers (251 individuals from 16 populations). We found five lineages with a clear phylogeographic structure but frequent secondary contact. Mitochondrial lineages were largely congruent with main population genetic clusters identified on microsatellite markers. In the upper Wami basin, populations are isolated as a putative Nothobranchius prognathus, but include also a population from a periphery of the middle Ruvu basin. Other four lineages (including putative Nothobranchius kwalensis) coexisted in secondary contact zones, but possessed clear spatial pattern. Main river channels did not form apparent barriers to dispersal. The most widespread lineage had strong signal of recent population expansion., Conclusions: We conclude that dispersal of a Nothobranchius species from a wet part of the genus distribution (tropical lowland) is not constrained by main river channels and closely related lineages frequently coexist in secondary contact zones. We also demonstrate contemporary connection between the Ruvu and Rufiji river basins. Our data do not provide genetic support for existence of recently described cryptic species from N. melanospilus complex, but cannot resolve this issue.

Samples of Austrolebias nigrofasciatus (n = 103), an endangered species of annual fish endemic to a small area of the Patos-Mirim lagoon system encompassing the São Gonçalo Channel lowlands, were collected from eight isolated temporary ponds, four located at the known distribution range of the species and four located along the Piratini River lowlands, where morphologically different individuals were found. In the laboratory, fragments of the mitochondrial cytochrome c oxidase I (coI), cytochrome b (cytb) and nuclear rhodopsin (rho) genes were amplified, purified and sequenced for 100, 99 and 58 of these individuals, respectively. Samples were further analysed using phylogenetic and phylogeographic methods to evaluate the patterns of genetic diversity and differentiation presented within and between populations, while assessing their evolutionary history, in order to guide the application of further conservation strategies. We found that the four new populations from the Piratini River lowlands encompass a different lineage of A. nigrofasciatus that diverged from that encountered in the São Gonçalo Channel at approximately 0.165 M years before present, during a population expansion and did not yet attain reciprocal monophyly. This divergence was associated with a glacial event that was preceded by an interglacial period putatively associated with the dispersal. Moreover, significant levels of genetic differentiation and a high number of exclusive haplotypes could be encountered even in micro-geographical scales, as in the comparisons between populations located within the same major lineage, indicating each of them may encompass independent management units. Conservation actions are certainly urgent, especially in the face of signs of a recent bottleneck., (© 2019 The Fisheries Society of the British Isles.)

The 2010 Deepwater Horizon oil spill released approximately 780 million liters of crude oil contaminating coastal habitats from Texas to Florida which are important habitats for many fish species during early life stages. These diverse habitats are also prone to rapid fluctuations in water quality, such as dissolved oxygen concentration and salinity. The consequence of combined exposure to crude oil and suboptimal environmental conditions during early life stage development of fish is still largely unknown. The objective of this project was to investigate the impacts of exposure to crude oil in combination with varying environmental stressors on developing Cyprinodon variegatus survival, growth, and gene expression. Three life stages (embryonic, post-hatch, and post-larval) were exposed to four nominal concentrations (6.25%, 12.5%, 50% and 100% with actual polycyclic aromatic hydrocarbon (PAH) concentrations ranging from 0 to 512 μg/L) of high energy water accommodated fractions (HEWAF) under different oxic (2.0 or >5.0 mg/L) and salinity (10 or 30 ppt) regimes at 30 °C for 48 h. We found that the post-larval developmental stage was the most sensitive to oil toxicity. Median lethal concentrations during the post-larval exposures followed a treatment-dependent pattern with the highest mortality observed under hypoxic-high salinity conditions (64.55 μg/L). Real-time PCR analysis identified down regulation of target genes, encoding cytochrome P450-1α (cyp1a1), erythropoietin (epo), and the aryl hydrocarbon receptor nuclear translocator (arnt1) only when oil exposure occurred under hypoxic-high salinity conditions in treatments with PAH concentrations greater than 226 μg/L. The target genes measured in this experiment are involved in the aryl hydrocarbon receptor signaling pathway which modulates metabolism of PAHs (a major component of crude oil), and the hypoxia inducible 1-α signaling pathway which is responsible for resilience to hypoxic stress, and it is known that disruption of these pathways can lead to an array of acute and chronic effects. Our results indicated that sheepshead minnow are most sensitive to oil exposure during the post-larval developmental stage. Survival data from this age-stage also indicate that oil toxicity response is exacerbated in hypoxic and high salinity environments. The increased mortality observed during the post-larval developmental stage might be attributed to the suppression of the aryl hydrocarbon receptor signaling and the hypoxia inducible 1-α signaling pathways which is evident in by the down-regulated expression of cyp1a1, epo, and arnt1. These findings provide more information about interactions between oil and abiotic factors which enable us to make better assumptions of the ecological impacts of DWH on coastal estuaries., (Copyright © 2019 Elsevier B.V. All rights reserved.)

To uncover the selective forces shaping life-history trait evolution across species, we investigate the genomic basis underlying adaptations to seasonal habitat desiccation in African killifishes, identifying the genetic variants associated with positive and relaxed purifying selection in 45 killifish species and 231 wild individuals distributed throughout sub-Saharan Africa. In annual species, genetic drift led to the expansion of nuclear and mitochondrial genomes and caused the accumulation of deleterious genetic variants in key life-history modulating genes such as mtor, insr, ampk, foxo3, and polg. Relaxation of purifying selection is also significantly associated with mitochondrial function and aging in human populations. We find that relaxation of purifying selection prominently shapes genomes and is a prime candidate force molding the evolution of lifespan and the distribution of genetic variants associated with late-onset diseases in different species. VIDEO ABSTRACT., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Genetic incompatibilities constitute the final stages of reproductive isolation and speciation, but little is known about incompatibilities that occur within recent adaptive radiations among closely related diverging populations. Crossing divergent species to form hybrids can break up coadapted variation, resulting in genetic incompatibilities within developmental networks shaping divergent adaptive traits. We crossed two closely related sympatric Cyprinodon pupfish species-a dietary generalist and a specialized molluscivore-and measured expression levels in their F1 hybrids to identify regulatory variation underlying the novel craniofacial morphology found in this recent microendemic adaptive radiation. We extracted mRNA from eight day old whole-larvae tissue and from craniofacial tissues dissected from 17-20 day old larvae to compare gene expression between a total of seven F1 hybrids and 24 individuals from parental species populations. We found 3.9% of genes differentially expressed between generalists and molluscivores in whole-larvae tissues and 0.6% of genes differentially expressed in craniofacial tissue. We found that 2.1% of genes were misregulated in whole-larvae hybrids whereas 19.1% of genes were misregulated in hybrid craniofacial tissues, after correcting for sequencing biases. We also measured allele specific expression across 15,429 heterozygous sites to identify putative compensatory regulatory mechanisms underlying differential expression between generalists and molluscivores. Together, our results highlight the importance of considering misregulation as an early indicator of genetic incompatibilities in the context of rapidly diverging adaptive radiations and suggests that compensatory regulatory divergence drives hybrid gene misregulation in developing tissues that give rise to novel craniofacial traits., Competing Interests: The authors have declared that no competing interests exist.

The nonapeptide arginine vasotocin (AVT) regulates osmotic balance in teleost fishes, but its mechanisms of action are not fully understood. Recently, it was discovered that nonapeptide receptors in teleost fishes are differentiated into two V1a-type, several V2-type, and two isotocin (IT) receptors, but it remains unclear which receptors mediate AVT's effects on gill osmoregulation. Here, we examined the role of nonapeptide receptors in the gill of the euryhaline Amargosa pupfish ( Cyprinodon nevadensis amargosae ) during osmotic acclimation. Transcripts for the teleost V1a-type receptor v1a2 were upregulated over fourfold in gill 24 h after transferring pupfish from 7.5 ppt to seawater (35 ppt) or hypersaline (55 ppt) conditions and downregulated after transfer to freshwater (0.3 ppt). Gill transcripts for the nonapeptide degradation enzyme leucyl-cystinyl aminopeptidase (LNPEP) also increased in fish acclimating to 35 ppt. To test whether the effects of AVT on the gill might be mediated by a V1a-type receptor, we administered AVT or a V1-type receptor antagonist (Manning compound) intraperitoneally to pupfish before transfer to 0.4 ppt or 35 ppt. Pupfish transferred to 35 ppt exhibited elevated gill mRNA abundance for cystic fibrosis transmembrane conductance regulator ( cftr ), but that upregulation diminished under V1-receptor inhibition. AVT inhibited the increase in gill Na + /Cl - cotransporter 2 ( ncc2 ) transcript abundance that occurs following transfer to hypoosmotic environments, whereas V1-type receptor antagonism increased ncc2 mRNAs even without a change in salinity. These findings indicate that AVT acts via a V1-type receptor to regulate gill Cl - transport by inhibiting Cl - uptake and facilitating Cl - secretion during seawater acclimation.

Pupfishes (genus Cyprinodon) evolved some of the broadest salinity tolerances of teleost fishes, with some taxa surviving in conditions from freshwater to nearly 160 ppt. In this study, we examined transcriptional dynamics of ion transporters and aquaporins in the gill of the desert Amargosa pupfish (Cyprinodon nevadensis amargosae) during rapid salinity change. Pupfish acclimated to 7.5 ppt were exposed to freshwater (0.3 ppt), seawater (35 ppt), or hypersaline (55 ppt) conditions over 4 h and sampled at these salinities over 14 d. Plasma osmolality and Cl - concentration became elevated 8 h after the start of exposure to 35 or 55 ppt but returned to baseline levels after 14 d. Osmolality recovery was paralleled by increased gill Na + /K + -ATPase activity and higher relative levels of messenger RNAs (mRNAs) encoding cystic fibrosis transmembrane conductance regulator (cftr) and Na + /K + /2Cl - cotransporter-1 (nkcc1). Transcripts encoding one Na + -HCO 3 - cotransporter-1 isoform (nbce1.1) also increased in the gills at higher salinities, while a second isoform (nbce1.2) increased expression in freshwater. Pupfish in freshwater also had lower osmolality and elevated gill mRNAs for Na + /H + exchanger isoform-2a (nhe2a) and V-type H + -ATPase within 8 h, followed by increases in Na + /H + exchanger-3 (nhe3), carbonic anhydrase 2 (ca2), and aquaporin-3 (aqp3) within 1 d. Gill mRNAs for Na + /Cl - cotransporter-2 (ncc2) also were elevated 14 d after exposure to 0.3 ppt. These results offer insights into how coordinated transcriptional responses for ion transporters in the gill facilitate reestablishment of osmotic homeostasis after changes in environmental salinity and provide evidence that the teleost gill expresses two Na + -HCO 3 - cotransporter-1 isoforms with different roles in freshwater and seawater acclimation.

The Danakil Depression in northeastern Africa represents one of the harshest arid environments on Earth, yet two genera of fishes, Danakilia (Cichlidae) and Aphanius (Cyprinodontidae), share its sparse aquatic habitats. The evolutionary history of these fishes is investigated here in the context of genetic, geological and paleoenvironmental information. We collected samples from seven sites and assessed phylogeographic relationships using concatenated COI and cytb mtDNA genes. Danakilia morphospecies show low differentiation at mitochondrial markers, but variation is partitioned between a northern cluster containing D. dinicolai plus three undescribed riverine populations, and a southern cluster including two creek populations of D. franchettii separated by the hypersaline waters of Lake Afrera. Aphanius displayed four genetically distinct clades (A. stiassnyae in Lake Afrera; one distributed across the entire area; one in Lake Abaeded; and one in the Shukoray River), but without clear large-scale geographic structure. However, Danakil Aphanius are clearly differentiated from A. dispar sensu stricto from the Sinai Peninsula. Geological evidence suggests that after the Late Pleistocene closure of the Danakil-Red Sea connection, increased post-glacial groundwater availability caused the formation of a brackish paleo-lake flooding the entire region below the -50 m contour. Fish populations previously isolated in coastal oases during glaciation were able to mix in the paleo-lake. Subsequently, in a more arid phase starting ∼7300 BP, paleo-lake regression isolated fishes in separate drainages, triggering their still ongoing diversification.

Annual killifishes are members of the Aplocheiloidea and live in ephemeral habitats that desiccate regularly during the dry season and refill during the rainy season. Populations of these fishes survive the dry season by producing drought-resistant diapausing eggs that are buried in the substrate. When the pool refills during the rainy season the juveniles hatch, grow rapidly and reproduce until the pool desiccates again during the next dry season. The association with such unpredictable habitats has led to the evolution to a variety of developmental adaptations such as a dispersed/reaggregation phase of the deep blastomeres, three possible diapause stages, extreme tolerance to high salinity and anoxia, an efficient DNA repair system and an extremely short life span. Here, we review the course of the dispersed/reaggregation phase, its evolution and phylogenetic distribution and diversity within the Aplocheiloidea. The phenomenon of blastomere dispersion/reaggregation in these fishes was first described in the 1960s and 70s. Blastomeres of most teleost fishes segregate into three groups that give rise to the enveloping cell layer, the yolk syncytial layer and the deep blastomeres that will form the embryo itself. When epiboly commences, the deep blastomeres form a more or less coherent cell sheet with a so called embryonic shield at it marginal zone marking the area where gastrulation takes place. In annual killifishes, the deep blastomeres segregate when epiboly starts and disperse when epiboly commences. After epiboly has been completed, the deep blastomeres are randomly distributed and migrate all over the enveloping cell layer. After several days they start to reaggregate and form the actual embryo that starts gastrulation. The evolutionary origin and mechanism behind this peculiar developmental pathway have puzzled developmental biologists for almost 50 years. However, several of these annual killifishes (Nothobranchius furzeri, Austrofundulus limnaeus, Austrolebias charrua and Austrolebias bellottii) have become model organisms in studies on developmental physiology, aging and stress tolerance. This has led to the establishment of modern genetic techniques such as transgenesis and cell fate mapping that are now used to tackle questions about the origin and mechanisms behind the dispersal/reaggregation phase., (Copyright © 2018. Published by Elsevier Inc.)

Secondary freshwater fish species inhabiting fluctuating and extreme environments are susceptible to changes in dispersion, effective population size, and genetic structure. The Iberian toothcarp Aphanius iberus is an endemic cyprinodontid of the Iberian Peninsula restricted to brackish water of salt marshes and coastal lagoons on the eastern Spanish Mediterranean coast. In this study, we analyzed mitochondrial cytochrome b (cyt b) DNA and microsatellite variation to evaluate ways in which the processes of extinction, dispersal, and colonization of A. iberus across its geographic distribution have affected its population genetic structure over time and space. The A. iberus network reconstruction indicated subtle levels of phylogeographic structuring. This, combined with substantial mitochondrial DNA (mtDNA) genetic diversity, suggests that Pleistocene glaciations had a lesser effect on the demographic structure of its populations than was the case for Iberian freshwater species with a similar distribution. Haplotype network, hierarchical analysis of molecular variance, and pairwise ΦST comparisons involving some Levantine samples showed a relatively high degree of mtDNA differentiation, which could be explained by historical isolation of the Villena Lagoon population. Conversely, significant genetic differentiation that follows an isolation-by-distance pattern, and a reduction in Ne though time was detected with microsatellites, suggesting extensive habitat fragmentation on the Mediterranean coast of the Iberian Peninsula over the past hundreds of years. At a smaller geographical scale (Mar Menor Lagoon), habitat fragmentation, probably due to human activity, appears to have resulted in substantially reduced migration and increased genetic drift, as shown by expanded genetic differentiation of populations.

Background: The annual killifish Austrofundulus limnaeus inhabits ephemeral ponds in northern Venezuela, South America, and is an emerging extremophile model for vertebrate diapause, stress tolerance, and evolution. Embryos of A. limnaeus regularly experience extended periods of desiccation and anoxia as a part of their natural history and have unique metabolic and developmental adaptations. Currently, there are limited genomic resources available for gene expression and evolutionary studies that can take advantage of A. limnaeus as a unique model system., Results: We describe the first draft genome sequence of A. limnaeus. The genome was assembled de novo using a merged assembly strategy and was annotated using the NCBI Eukaryotic Annotation Pipeline. We show that the assembled genome has a high degree of completeness in genic regions that is on par with several other teleost genomes. Using RNA-seq and phylogenetic-based approaches, we identify several candidate genes that may be important for embryonic stress tolerance and post-diapause development in A. limnaeus. Several of these genes include heat shock proteins that have unique expression patterns in A. limnaeus embryos and at least one of these may be under positive selection., Conclusion: The A. limnaeus genome is the first South American annual killifish genome made publicly available. This genome will be a valuable resource for comparative genomics to determine the genetic and evolutionary mechanisms that support the unique biology of annual killifishes. In a broader context, this genome will be a valuable tool for exploring genome-environment interactions and their impacts on vertebrate physiology and evolution.

The Caatinga is the largest nucleus of seasonally dry tropical forests in South America, but little is known about the evolutionary history and biogeography of endemic organisms. Evolutionary diversification and distribution of terrestrial vertebrates endemic to the Caatinga have been explained by palaeogeographical Neogene episodes, mostly related to changes in the course of the São Francisco River, the largest river in the region. Our objective is to estimate the timing of divergence of two endemic groups of short-lived seasonal killifishes inhabiting all ecoregions of the Caatinga, testing the occurrence of synchronic events of spatial diversification in light of available data on regional palaeogeography. We performed independent time-calibrated phylogenetic molecular analyses for two clades of sympatric and geographically widespread seasonal killifishes endemic to the Caatinga, the Hypsolebias antenori group and the Cynolebias alpha-clade. Our results consistently indicate that species diversification took place synchronically in both groups, as well as it is contemporary to diversification of other organisms adapted to life in the semi-arid Caatinga, including lizards and small mammals. Both groups originated during the Miocene, but species diversification started between the Late Miocene and Early Pliocene, when global cooling probably favoured the expansion of semi-arid areas. Synchronic diversification patterns found are chronologically related to Tertiary palaeogeographical reorganizations associated to continental drift and to Quaternary climatic changes, corroborating the recent proposal that South American biodiversity has been continuously shaped between the Late Paleogene and Pleistocene.

Background: The mangrove killifish Kryptolebias marmoratus is the only vertebrate that reproduces by self-fertilizing and is an important model species in genetics and marine ecotoxicology. Using whole-genome and transcriptome sequences, we identified all members of the cytochrome P450 (CYP) family in this model teleost and compared them with those of other teleosts., Results: A total of 74 cytochrome P450 genes and one pseudogene were identified in K. marmoratus. Phylogenetic analysis indicated that the CYP genes in clan 2 were most expanded, while synteny analysis with other species showed orthologous relationships of CYP subfamilies among teleosts. In addition to the CYP2K expansions, five tandem duplicated gene copies of CYP5A were observed. These features were unique to K. marmoratus., Conclusions: These results shed a light on CYP gene evolution, particularly the co-localized CYP2K, CYP5A, and CYP46A subfamilies in fish. Future studies of CYP expression could identify specific endogenous and exogenous environmental factors that triggered the evolution of tandem CYP duplication in K. marmoratus.

The rich biological diversity of South America has motivated a series of studies associating evolution of endemic taxa with the dramatic geologic and climatic changes that occurred during the Cainozoic. The organism here studied is the killifish tribe Cynolebiini, a group of seasonal fishes uniquely inhabiting temporary pools formed during the rainy seasons. The Cynolebiini are found in open vegetation areas inserted in the main tropical and subtropical South American phytogeographical regions east of the Andes. Here, we present the first molecular phylogeny sampling all the eight genera of the Cynolebiini, using fragments of two mitochondrial and four nuclear genes for 35 species of Cynolebiini plus 19 species as outgroups. The dataset, 4448bp, was analysed under Bayesian and maximum likelihood approaches, providing a relatively well solved tree, which retrieves high support values for the Cynolebiini and most included clades. The resulting tree was used to estimate the time of divergence in included lineages using two cyprinodontiform fossils to calibrate the tree. We further investigated historical biogeography through the likelihood-based DEC model. Our estimates indicate that divergence between the clades comprising New World and Old World aplocheiloids occurred during the Eocene, about 50Mya, much more recent than the Gondwanan fragmentation scenario assumed in previous studies. This estimation is nearly synchronous to estimated splits involving other South American and African vertebrate clades, which have been explained by transoceanic dispersal through an ancient Atlantic island chain during the Palaeogene. We estimate that Cynolebiini split from its sister group Cynopoecilini in the Oligocene, about 25Mya and that Cynolebiini started to diversify giving origin to the present genera during the Miocene, about 20-14Mya. The Cynolebiini had an ancestral origin in the Atlantic Forest and probably were not present in the open vegetation formations of central and northeastern South America until the Middle Miocene, when expansion of dry open vegetation was favoured by cool temperatures and strike seasonality. Initial splitting between the genera Cynolebias and Simpsonichthys during the Miocene (about 14Mya) is attributed to the uplift of the Central Brazilian Plateau., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Much of the current research on longevity focuses on the aging process within a single species. Several molecular players (e.g. IGF1 and MTOR), pharmacological compounds (e.g. rapamycin and metformin), and dietary approaches (e.g. calorie restriction and methionine restriction) have been shown to be important in regulating and modestly extending lifespan in model organisms. On the other hand, natural lifespan varies much more significantly across species. Within mammals alone, maximum lifespan differs more than 100 fold, but the underlying regulatory mechanisms remain poorly understood. Recent comparative studies are beginning to shed light on the molecular signatures associated with exceptional longevity. These include genome sequencing of microbats, naked mole rat, blind mole rat, bowhead whale and African turquoise killifish, and comparative analyses of gene expression, metabolites, lipids and ions across multiple mammalian species. Together, they point towards several putative strategies for lifespan regulation and cancer resistance, as well as the pathways and metabolites associated with longevity variation. In particular, longevity may be achieved by both lineage-specific adaptations and common mechanisms that apply across the species. Comparing the resulting cross-species molecular signatures with the within-species lifespan extension strategies will improve our understanding of mechanisms of longevity control and provide a starting point for novel and effective interventions., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Rapid diversification often involves complex histories of gene flow that leave variable and conflicting signatures of evolutionary relatedness across the genome. Identifying the extent and source of variation in these evolutionary relationships can provide insight into the evolutionary mechanisms involved in rapid radiations. Here we compare the discordant evolutionary relationships associated with species phenotypes across 42 whole genomes from a sympatric adaptive radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas and several outgroup pupfish species in order to understand the rarity of these trophic specialists within the larger radiation of Cyprinodon. 82% of the genome depicts close evolutionary relationships among the San Salvador Island species reflecting their geographic proximity, but the vast majority of variants fixed between specialist species lie in regions with discordant topologies. Top candidate adaptive introgression regions include signatures of selective sweeps and adaptive introgression of genetic variation from a single population in the northwestern Bahamas into each of the specialist species. Hard selective sweeps of genetic variation on San Salvador Island contributed 5 times more to speciation of trophic specialists than adaptive introgression of Caribbean genetic variation; however, four of the 11 introgressed regions came from a single distant island and were associated with the primary axis of oral jaw divergence within the radiation. For example, standing variation in a proto-oncogene (ski) known to have effects on jaw size introgressed into one San Salvador Island specialist from an island 300 km away approximately 10 kya. The complex emerging picture of the origins of adaptive radiation on San Salvador Island indicates that multiple sources of genetic variation contributed to the adaptive phenotypes of novel trophic specialists on the island. Our findings suggest that a suite of factors, including rare adaptive introgression, may be necessary for adaptive radiation in addition to ecological opportunity.

A new genus and species, Pseudorestias lirimensis, is described from the southern part of the Chilean Altiplano. While sharing several characters that clearly align the new species with Orestias, this new fish is characterized by numerous autapomorphies: the Meckel cartilage is a continuous cartilage that broadly expands posteriorly (in large specimens, it keeps its anterior part and is resorbed posteriorly), the basibranchials are fused into one long element, the second pharyngobranchial is not displaced dorsally over pharyngobranchial tooth plate 3+4, but they are aligned, the anterior and posterior ceratohyals are closely articulated keeping a scarce amount of cartilage between both bones and ventral to them, ossified middle and distal dorsal radials are present in females as well as ossified middle and distal anal radials. Pseudorestias lirimensis presents strong sexual dimorphism associated to size. Females are almost twice as large and long than males, neuromast lines are absent in males, a mesethmoid is present in males, squamation on head is reduced in males, and ossified middle and distal radial of dorsal fin are cartilaginous in males. Pseudorestias and Orestias are suggested as the sole members of the tribe Orestiini. A list of characters diagnosing the tribe is provided. The presence of the new genus is interpreted as a possible result of the ecosystem isolation where the fish is living from surrounding basins-as early as possibly from the Miocene-Pliocene times-and its physical and chemical characteristics. Small populations, living conditions, small habitat, and reduced distribution make this species a strong candidate to be considered critically endangered, a situation already established for all other Chilean species living in the Altiplano. There is high probability it will become extinct due to water demands and climate change in the region.

Vertebrates exhibit substantial variation in eye size. Eye size correlates positively with visual capacity and behaviors that enhance fitness, such as predator avoidance. This foreshadows a connection between predation and eye size evolution. Yet, the conditions that favor evolutionary shifts in eye size, besides the well-known role for light availability, are unclear. We tested the influence of predation on the evolution of eye size in Trinidadian killifish, Rivulus hartii. Rivulus are located across a series of communities where they coexist with visually oriented piscivores ("high predation" sites), and no predators ("Rivulus-only" sites). Wild-caught Rivulus from high predation sites generally exhibited a smaller relative eye size than communities that lack predators. Yet, such differences were inconsistent across rivers. Second-generation common garden reared fish revealed repeatable decreases in eye size in Rivulus from high predation sites. We performed additional experiments that tested the importance of light and resources on eye size evolution. Sites that differ in light or resource availability did not differ in eye size. Our results argue that differences in predator-induced mortality underlie genetically-based shifts in vertebrate eye size. We discuss the drivers of eye size evolution in light of the nonparallel trends between the phenotypic and common garden results., (© 2017 The Author(s). Evolution © 2017 The Society for the Study of Evolution.)