Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape.

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.

Dynamic regulation of gonadal transposon control across the lifespan of the naturally short-lived African turquoise killifish.

Although germline cells are considered to be functionally "immortal," both the germline and supporting somatic cells in the gonad within an organism experience aging. With increased age at parenthood, the age-related decline in reproductive success has become an important biological issue for an aging population. However, molecular mechanisms underlying reproductive aging across sexes in vertebrates remain poorly understood. To decipher molecular drivers of vertebrate gonadal aging across sexes, we perform longitudinal characterization of the gonadal transcriptome throughout the lifespan in the naturally short-lived African turquoise killifish (Nothobranchius furzeri). By combining mRNA-seq and small RNA-seq from 26 individuals, we characterize the aging gonads of young-adult, middle-aged, and old female and male fish. We analyze changes in transcriptional patterns of genes, transposable elements (TEs), and piRNAs. We find that testes seem to undergo only marginal changes during aging. In contrast, in middle-aged ovaries, the time point associated with peak female fertility in this strain, PIWI pathway components are transiently down-regulated, TE transcription is elevated, and piRNA levels generally decrease, suggesting that egg quality may already be declining at middle-age. Furthermore, we show that piRNA ping-pong biogenesis declines steadily with age in ovaries, whereas it is maintained in aging testes. To our knowledge, this data set represents the most comprehensive transcriptomic data set for vertebrate gonadal aging. This resource also highlights important pathways that are regulated during reproductive aging in either ovaries or testes, which could ultimately be leveraged to help restore aspects of youthful reproductive function.

Molecular mechanisms of embryonic tail development in the self-fertilizing mangrove killifish Kryptolebias marmoratus.

Using the self-fertilizing mangrove killifish, we characterized two mutants, shorttail (stl) and balltail (btl). These mutants showed abnormalities in the posterior notochord and muscle development. Taking advantage of a highly inbred isogenic strain of the species, we rapidly identified the mutated genes, noto and msgn1 in the stl and btl mutants, respectively, using a single lane of RNA sequencing without the need of a reference genome or genetic mapping techniques. Next, we confirmed a conserved morphant phenotype in medaka and demonstrate a crucial role of noto and msgn1 in cell sorting between the axial and paraxial part of the tail mesoderm. This novel system could substantially accelerate future small-scale forward-genetic screening and identification of mutations. Therefore, the mangrove killifish could be used as a complementary system alongside existing models for future molecular genetic studies.

Niche conservatism and spread explain introgression between native and invasive fish.

Hybridisation can be an important driver of evolutionary change, but hybridisation with invasive species can have adverse effects on native biodiversity. While hybridisation has been documented across taxa, there is limited understanding of ecological factors promoting patterns of hybridisation and the spatial distribution of hybrid individuals. We combined the results of ecological niche modelling (ENM) and restriction site-associated DNA sequencing to test theories of niche conservatism and biotic resistance on the success of invasion, admixture, and extent of introgression between native and non-native fishes. We related Maxent predictions of habitat suitability based on the native ranges of invasive Eastern Banded Killifish (Fundulus diaphanus diaphanus Lesueur 1817) and native Western Banded Killifish (Fundulus diaphanus menona Jordan and Copeland 1877) to admixture indices of individual Banded Killifish. We found that Eastern Banded Killifish predominated at sites predicted as suitable from their ENM, consistent with niche conservatism. Admixed individuals were more common as Eastern Banded Killifish habitat suitability declined. We also found that Eastern Banded Killifish were most common at sites closest to the presumed source of this invasion, whereas the proportion of admixed individuals increased with distance from the source of invasion. Lastly, we found little evidence that habitat suitability for Western Banded Killifish provides biotic resistance from either displacement by, or admixture with, invasive Eastern Banded Killifish. Our study demonstrates that ENMs can inform conservation-relevant outcomes between native and invasive taxa while emphasising the importance of protecting isolated Western Banded Killifish populations from invasive conspecifics.

Centrality to the metapopulation is more important for population genetic diversity than habitat area or fragmentation.

Drift and gene flow affect genetic diversity. Given that the strength of genetic drift increases as population size decreases, management activities have focused on increasing population size through preserving habitats to preserve genetic diversity. Few studies have empirically evaluated the impacts of drift and gene flow on genetic diversity. Kryptolebias marmoratus, henceforth 'rivulus', is a small killifish restricted to fragmented New World mangrove forests with gene flow primarily associated with ocean currents. Rivulus form distinct populations across patches, making them a well-suited system to test the extent to which habitat area, fragmentation and connectivity are associated with genetic diversity. Using over 1000 individuals genotyped at 32 microsatellite loci, high-resolution landcover data and oceanographic simulations with graph theory, we demonstrate that centrality (connectivity) to the metapopulation is more strongly associated with genetic diversity than habitat area or fragmentation. By comparing models with and without centrality standardized by the source population's genetic diversity, our results suggest that metapopulation centrality is critical to genetic diversity regardless of the diversity of adjacent populations. While we find evidence that habitat area and fragmentation are related to genetic diversity, centrality is always a significant predictor with a larger effect than any measure of habitat configuration.

Impacts on Atlantic Killifish from Neurotoxicants: Genes, Behavior, and Population-Relevant Outcomes.

Molecular, cellular, and organismal alterations are important descriptors of toxic effects, but our ability to extrapolate and predict ecological risks is limited by the availability of studies that link measurable end points to adverse population relevant outcomes such as cohort survival and growth. In this study, we used laboratory gene expression and behavior data from two populations of Atlantic killifish Fundulus heteroclitus [one reference site (SCOKF) and one PCB-contaminated site (NBHKF)] to inform individual-based models simulating cohort growth and survival from embryonic exposures to environmentally relevant concentrations of neurotoxicants. Methylmercury exposed SCOKF exhibited brain gene expression changes in the si:ch211-186j3.6, si:dkey-21c1.4, scamp1, and klhl6 genes, which coincided with changes in feeding and swimming behaviors, but our models simulated no growth or survival effects of exposures. PCB126-exposed SCOKF had lower physical activity levels coinciding with a general upregulation in nucleic and cellular brain gene sets (BGS) and downregulation in signaling, nucleic, and cellular BGS. The NBHKF, known to be tolerant to PCBs, had altered swimming behaviors that coincided with 98% fewer altered BGS. Our models simulated PCB126 decreased growth in SCOKF and survival in SCOKF and NBHKF. Overall, our study provides a unique demonstration linking molecular and behavioral data to develop quantitative, testable predictions of ecological risk.

Fast satellite DNA evolution in Nothobranchius annual 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.

Persistent interactive effects of developmental salinity and temperature in Atlantic killifish (Fundulus heteroclitus).

Climate change alters multiple abiotic environmental factors in aquatic environments but relatively little is known about their interacting impacts, particularly in developing organisms where these exposures have the potential to cause long-lasting effects. To explore these issues, we exposed developing killifish embryos (Fundulus heteroclitus) to 26 °C or 20 °C and 20 ppt or 3 ppt salinity in a fully-factorial design. After hatching, fish were transferred to common conditions of 20 °C and 20 ppt to assess the potential for persistent developmental plasticity. Warm temperature increased hatching success and decreased hatch time, whereas low salinity negatively affected hatching success, but this was only significant in fish developed at 20 °C. Temperature, salinity, or their interaction affected mRNA levels of genes typically associated with thermal and hypoxia tolerance (hif1a, hsp90b, hsp90a, hsc70, and hsp70.2) across multiple developmental timepoints. These differences were persistent into the juvenile stage, where the fish that developed at 26 °C had higher expression of hif1a, hsp90b, hsp90a, and hsp70.2 than fish developed at 20 °C, and this was particularly evident for the group developed at both high temperature and salinity. There were also long-lasting effects of developmental treatments on body size after four months of rearing under common conditions. Fish developed at low salinity or temperature were larger than fish developed at high temperature or salinity, but there was no interaction between the two factors. These data highlight the complex nature of the developmental effects of interacting stressors which has important implications for predicting the resilience of fishes in the context of climate change.

Validation of qPCR reference genes in the endangered annual killifish Austrolebias charrua considering different tissues, gender and environmental conditions.

The annual killifish Austrolebias charrua is an endangered species, endemic to the southern region of South America, which inhabits temporary ponds that emerges in the rainy season. The main anthropogenic threat driving the extinction of A. charrua stems from extensive agriculture, primarily due to the widrespread use of glyphosate-based herbicides near their habitats. Annual killifishes have been used as models for ecotoxicological studies but, up to now, there are no studies about reference genes in any Austrolebias species. This represents an obstacle to the use of qPCR-based technologies, the standard method for gene expression quantification. The present study aimed to select and validate potential reference genes for qPCR normalization in the annual killifish Austrolebias charrua considering different tissues, gender and environmental conditions. The candidate reference genes 18 s, actb, gapdh, ef1a, shox, eif3g, and the control gene atp1a1 were evaluated in male and female individuals in three different tissues (brain, liver, and gills) under two experimental conditions (control and acute exposition to Roundup Transorb®). The collected tissues were submitted to RNA extraction, followed by cDNA synthesis, cloning, sequencing, and qPCR. Overall, 18 s was the most stable reference gene, and 18 s and ef1a were the most stable combination. Otherwise, considering all variables, gapdh and shox were the least stable candidate genes. Foremost, suitable reference genes were validated in A. charrua, facilitating accurate mRNA quantification in this species, which might be useful for developing molecular tools of ecotoxicological assessment based on gene expression analysis for environmental monitoring of annual killifish.

Differential DNA methylation and metabolite profiling of Atlantic killifish (Fundulus heteroclitus) from the New Bedford Harbor Superfund site.

Atlantic killifish (Fundulus heteroclitus) is a valuable model in evolutionary toxicology to study how the interactions between genetic and environmental factors serve the adaptive ability of organisms to resist chemical pollution. Killifish populations inhabiting environmental toxicant-contaminated New Bedford Harbor (NBH) show phenotypes tolerant to polychlorinated biphenyls (PCBs) and differences at the transcriptional and genomic levels. However, limited research has explored epigenetic alterations and metabolic effects in NBH killifish. To identify the involvement of epigenetic and metabolic regulation in the adaptive response of killifish, we investigated tissue- and sex-specific differences in global DNA methylation and metabolomic profiles of NBH killifish populations, compared to sensitive populations from a non-polluted site, Scorton Creek (SC). The results revealed that liver-specific global DNA hypomethylation and differential metabolites were evident in fish from NBH compared with those from SC. The sex-specific differences were not greater than the tissue-specific differences. We demonstrated liver-specific enriched metabolic pathways (e.g., amino acid metabolic pathways converged into the urea cycle and glutathione metabolism), suggesting possible crosstalk between differential metabolites and DNA hypomethylation in the livers of NBH killifish. Additional investigation of methylated gene regions is necessary to understand the functional role of DNA hypomethylation in the regulation of enzyme-encoding genes associated with metabolic processes and physiological changes in NBH populations.

A Pilot Study of Gene Expression Modulation from Antioxidant System of Killifish Austrolebias charrua After Exposure to Roundup Transorb®.

Roundup Transorb® (RDT) is the most popular glyphosate-based herbicide (GHB) used in agriculture, and its impact extends to non-target organisms. The annual killifish Austrolebias charrua is an endangered species endemic to southern South America and inhabits temporary ponds. This study evaluates the effects of RDT concentrations (0.065 and 5 mg/L GAE) on A. charrua exposed for 96 h. Gene expression of cat, sod2, gstα, gclc, and ucp1 was evaluated on the liver and gills. Highlighting that even at low concentrations permitted by Brazilian legislation, the RDT can have adverse effects on A. charrua.

RNAseq Analysis of Brain Aging in Wild Specimens of Short-Lived Turquoise Killifish: Commonalities and Differences With Aging Under Laboratory Conditions.

A vast body of studies is available that describe age-dependent gene expression in relation to aging in a number of different model species. These data were obtained from animals kept in conditions with reduced environmental challenges, abundant food, and deprivation of natural sensory stimulation. Here, we compared wild- and captive aging in the short-lived turquoise killifish (Nothobranchius furzeri). These fish inhabit temporary ponds in the African savannah. When the ponds are flooded, eggs hatch synchronously, enabling a precise timing of their individual and population age. We collected the brains of wild fish of different ages and quantified the global age-dependent regulation of transcripts using RNAseq. A major difference between captive and wild populations is that wild populations had unlimited access to food and hence grew to larger sizes and reached asymptotic size more rapidly, enabling the analysis of age-dependent gene expression without the confounding effect of adult brain growth. We found that the majority of differentially expressed genes show the same direction of regulation in wild and captive populations. However, a number of genes were regulated in opposite direction. Genes downregulated in the wild and upregulated in captivity were enriched for terms related to neuronal communication. Genes upregulated in the wild and downregulated in captive conditions were enriched in terms related to DNA replication. Finally, the rate of age-dependent gene regulation was higher in wild animals, suggesting a phenomenon of accelerated aging.

Evidence of prezygotic isolation, but not assortative mating, between locally adapted populations of Fundulus heteroclitus across a salinity gradient.

Selection along environmental gradients can drive reproductive isolation and speciation. Among fishes, salinity is a major factor limiting species distributions, and despite its importance in generating species diversity, speciation events between marine and freshwater are rare. Here, we tested for mechanisms of reproductive isolation between locally adapted freshwater and brackish water-native populations of killifish, Fundulus heteroclitus, from either side of a hybrid zone along a salinity gradient. There was evidence for pre-zygotic endogenous reproductive isolation with reduced fertilization success between crosses of freshwater-native males and brackish water-native females. Exogenous pre-zygotic isolation was also present where females had highest fertilization in their native salinity. We used a replicated mass spawning design to test for mate choice in both brackish and fresh water. After genotyping 187 parents and 2523 offspring at 2347 SNPs across the genome, 85% of offspring were successfully assign to their parents. However, no reinforcing mate choice was observed. These results therefore demonstrate emerging, yet limited, reproductive isolation and incipient speciation across a marine to freshwater salinity gradient and suggest that both endogenous and exogenous mechanisms, but not assortative mating, contribute to divergence.

Reproductive isolating mechanisms contributing to asymmetric hybridization in Killifishes (Fundulus spp.).

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.

Sex chromosome differentiation via changes in the Y chromosome repeat landscape in African annual killifishes Nothobranchius furzeri and N. kadleci.

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.

Genome sequencing and transcriptomic analysis of the Andean killifish Orestias ascotanensis reveals adaptation to high-altitude aquatic life.

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.

Fish models for investigating nutritional regulation of embryonic development.

In recent decades, biologist have focused on the spatiotemporal regulation and function of genes to understand embryogenesis. It is clear that maternal diet impacts fetal development but how nutrients, like lipids and vitamins, modify developmental programs is not completely understood. Fish are useful research organisms for such investigations. Most species of fish produce eggs that develop outside the mother, dependent on a finite amount of yolk to form and grow. The developing embryo is a closed system that can be readily biochemically analyzed, easily visualized, and manipulated to understand the role of nutrients in tissue specification, organogenesis, and growth. Natural variation in yolk composition observed across fish species may be related to unique developmental strategies. In this review, we discuss the reasons that teleost fishes are powerful models to understand nutritional control of development and highlight three species that are particularly valuable for future investigations: the zebrafish, Danio rerio, the African Killifish, Nothobranchius furzeri, and the Mexican tetra, Astyanax mexicanus. This review is a part of a special issue on nutritional, hormonal, and metabolic drivers of development.

Lifespan and telomere length variation across populations of wild-derived African killifish.

Telomeres and telomerase prevent the continuous erosion of chromosome-ends caused by lifelong cell division. Shortened telomeres are associated with age-related pathologies. While short telomere length is positively correlated with increased lethality at the individual level, in comparisons across species short telomeres are associated with long (and not short) lifespans. Here, we tested this contradiction between individual and evolutionary patterns in telomere length using African annual killifish. We analysed lifespan and telomere length in a set of captive strains derived from well-defined wild populations of Nothobranchius furzeri and its sister species, N. kadleci, from sites along a strong gradient of aridity which ultimately determines maximum natural lifespan. Overall, males were shorter-lived than females, and also had shorter telomeres. Male lifespan (measured in controlled laboratory conditions) was positively associated with the amount of annual rainfall in the site of strain origin. However, fish from wetter climates had shorter telomeres. In addition, individual fish which grew largest over the juvenile period possessed shorter telomeres at the onset of adulthood. This demonstrates that individual condition and environmentally-driven selection indeed modulate the relationship between telomere length and lifespan in opposite directions, validating the existence of inverse trends within a single taxon. Intraindividual heterogeneity of telomere length (capable to detect very short telomeres) was not associated with mean telomere length, suggesting that the shortest telomeres are controlled by regulatory pathways other than those that determine mean telomere length. The substantial variation in telomere length between strains from different environments identifies killifish as a powerful system in understanding the adaptive value of telomere length.

Phylogenomics reveals extensive introgression and a case of mito-nuclear discordance in the killifish genus Kryptolebias.

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.

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

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