Swordtail fish hybrids reveal that genome evolution is surprisingly predictable after initial hybridization.

Over the past 2 decades, biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common-not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past approximately 40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni × X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis.

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Genome biology of the darkedged splitfin, Girardinichthys multiradiatus, and the evolution of sex chromosomes and placentation.

Viviparity evolved independently about 150 times in vertebrates and more than 20 times in fish. Several lineages added to the protection of the embryo inside the body of the mother, the provisioning of nutrients, and physiological exchange. This often led to the evolution of a placenta. Among fish, one of the most complex systems serving the function of the placenta is the embryonal trophotaenia/ovarian luminal epithelium of the goodeid fishes. For a better understanding of this feature and others of this group of fishes, high-quality genomic resources are essential. We have sequenced the genome of the darkedged splitfin, Girardinichthys multiradiatus The assembly is chromosome level and includes the X and Y Chromosomes. A large male-specific region on the Y was identified covering 80% of Chromosome 20, allowing some first inferences on the recent origin and a candidate male sex determining gene. Genome-wide transcriptomics uncovered sex-specific differences in brain gene expression with an enrichment for neurosteroidogenesis and testis genes in males. The expression signatures of the splitfin embryonal and maternal placenta showed overlap with homologous tissues including human placenta, the ovarian follicle epithelium of matrotrophic poeciliid fish species and the brood pouch epithelium of the seahorse. Our comparative analyses on the evolution of embryonal and maternal placenta indicate that the evolutionary novelty of maternal provisioning development repeatedly made use of genes that already had the same function in other tissues. In this way, preexisting modules are assembled and repurposed to provide the molecular changes for this novel trait.

Predictability and parallelism in the contemporary evolution of hybrid genomes.

Hybridization between species is widespread across the tree of life. As a result, many species, including our own, harbor regions of their genome derived from hybridization. Despite the recognition that this process is widespread, we understand little about how the genome stabilizes following hybridization, and whether the mechanisms driving this stabilization tend to be shared across species. Here, we dissect the drivers of variation in local ancestry across the genome in replicated hybridization events between two species pairs of swordtail fish: Xiphophorus birchmanni × X. cortezi and X. birchmanni × X. malinche. We find unexpectedly high levels of repeatability in local ancestry across the two types of hybrid populations. This repeatability is attributable in part to the fact that the recombination landscape and locations of functionally important elements play a major role in driving variation in local ancestry in both types of hybrid populations. Beyond these broad scale patterns, we identify dozens of regions of the genome where minor parent ancestry is unusually low or high across species pairs. Analysis of these regions points to shared sites under selection across species pairs, and in some cases, shared mechanisms of selection. We show that one such region is a previously unknown hybrid incompatibility that is shared across X. birchmanni × X. cortezi and X. birchmanni × X. malinche hybrid populations.

Paternal source of germ plasm determinants in the viviparous teleost, Gambusia holbrooki; dads do matter.

The identity of germ cells, the progenitors of life, is thought to be acquired by two modes; either by maternal signals (preformed) or induced de novo from pluripotent cells (epigenesis) in the developing embryos. However, paternal roles seem enshrouded or completely overlooked in this fundamental biological process. Hence, we investigated the presence of germplasm transcripts in the sperm of Gambusia holbrooki, a live-bearing fish, demonstrating their presence and suggesting paternal contributions. Interestingly, not all germplasm markers were present (nanos1 and tdrd6) in the sperm, but some were conspicuous (dazl, dnd-α, piwi II, and vasa), indicating that the latter is required for establishing germ cell identity in the progeny, with a possible parent-specific role. Furthermore, there were also spatial differences in the distribution of these determinants, suggesting additional roles in sperm physiology and/or fertility. Our results support the hypothesis that dads also play a vital role in establishing the germ cell identity, especially in G. holbrooki, which shares elements of both preformation and induction modes of germline determination. This, coupled with its life history traits, makes G. holbrooki an excellent system for dissecting evolutionary relationships between the two germline determination modes, their underpinning mechanisms and ultimately the perpetuity of life.

Genomes of the Orestias pupfish from the Andean Altiplano shed light on their evolutionary history and phylogenetic relationships within Cyprinodontiformes.

BACKGROUND: To unravel the evolutionary history of a complex group, a comprehensive reconstruction of its phylogenetic relationships is crucial. This requires meticulous taxon sampling and careful consideration of multiple characters to ensure a complete and accurate reconstruction. The phylogenetic position of the Orestias genus has been estimated partly on unavailable or incomplete information. As a consequence, it was assigned to the family Cyprindontidae, relating this Andean fish to other geographically distant genera distributed in the Mediterranean, Middle East and North and Central America. In this study, using complete genome sequencing, we aim to clarify the phylogenetic position of Orestias within the Cyprinodontiformes order. RESULTS: We sequenced the genome of three Orestias species from the Andean Altiplano. Our analysis revealed that the small genome size in this genus (~ 0.7 Gb) was caused by a contraction in transposable element (TE) content, particularly in DNA elements and short interspersed nuclear elements (SINEs). Using predicted gene sequences, we generated a phylogenetic tree of Cyprinodontiformes using 902 orthologs extracted from all 32 available genomes as well as three outgroup species. We complemented this analysis with a phylogenetic reconstruction and time calibration considering 12 molecular markers (eight nuclear and four mitochondrial genes) and a stratified taxon sampling to consider 198 species of nearly all families and genera of this order. Overall, our results show that phylogenetic closeness is directly related to geographical distance. Importantly, we found that Orestias is not part of the Cyprinodontidae family, and that it is more closely related to the South American fish fauna, being the Fluviphylacidae the closest sister group. CONCLUSIONS: The evolutionary history of the Orestias genus is linked to the South American ichthyofauna and it should no longer be considered a member of the Cyprinodontidae family. Instead, we submit that Orestias belongs to the Orestiidae family, as suggested by Freyhof et al. (2017), and that it is the sister group of the Fluviphylacidae family, distributed in the Amazonian and Orinoco basins. These two groups likely diverged during the Late Eocene concomitant with hydrogeological changes in the South American landscape.

Genomic Signatures Associated with Transitions to Viviparity in Cyprinodontiformes.

The transition from oviparity to viviparity has occurred independently over 150 times across vertebrates, presenting one of the most compelling cases of phenotypic convergence. However, whether the repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common signature in genomic divergence remains largely unknown. Although recent investigations into the evolution of viviparity have demonstrated striking similarity among the genes and molecular pathways involved across disparate vertebrate groups, quantitative tests for genome-wide convergent have provided ambivalent answers. Here, we investigate the potential role of molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo genomes and utilized publicly available genomes of viviparous and oviparous species to test for molecular convergence across both coding and noncoding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and in rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, both statistical power and biological confounds could constrain our ability to detect significant correlated evolution. We therefore identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. Motif enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development, and immunity occurred alongside the evolution of viviparity. Overall, however, our findings indicate that independent transitions to viviparity in these fish are not strongly associated with an excess of molecular convergence, but a few genes show convincing evidence of convergent evolution.

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

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

Evidence for gene flow and trait reversal during radiation of Mexican Goodeid fish.

Understanding the phylogeographic history of a group and identifying the factors contributing to speciation is an important challenge in evolutionary biology. The Goodeinae are a group of live-bearing fishes endemic to Mexico. Here, we develop genomic resources for species within the Goodeinae and use phylogenomic approaches to characterise their evolutionary history. We sequenced, assembled and annotated the genomes of four Goodeinae species, including Ataeniobius toweri, the only matrotrophic live-bearing fish without a trophotaenia in the group. We estimated timings of species divergence and examined the extent and timing of introgression between the species to assess if this may have occurred during an early radiation, or in more recent episodes of secondary contact. We used branch-site models to detect genome-wide positive selection across Goodeinae, and we specifically asked whether this differs in A. toweri, where loss of placental viviparity has recently occurred. We found evidence of gene flow between geographically isolated species, suggesting vicariant speciation was supplemented by limited post-speciation gene flow, and gene flow may explain previous uncertainties about Goodeid phylogeny. Genes under positive selection in the group are likely to be associated with the switch to live-bearing. Overall, our studies suggest that both volcanism-driven vicariance and changes in reproductive mode influenced radiation in the Goodeinae.

Warming, stochastic diel thermal fluctuations affect physiological performance and gill plasticity in an amphibious mangrove fish.

Natural temperature variation in many marine ecosystems is stochastic and unpredictable, and climate change models indicate that this thermal irregularity is likely to increase. Temperature acclimation may be more challenging when conditions are highly variable and stochastic, and there is a need for empirical physiological data in these thermal environments. Using the hermaphroditic, amphibious mangrove rivulus (Kryptolebias marmoratus), we hypothesized that compared with regular, warming diel thermal fluctuations, stochastic warm fluctuations would negatively affect physiological performance. To test this, we acclimated fish to: (1) non-stochastic and (2) stochastic thermal fluctuations with a similar thermal load (27-35°C), and (3) a stable/consistent control temperature at the low end of the cycle (27°C). We determined that fecundity was reduced in both cycles, with reproduction ceasing in stochastic thermal environments. Fish acclimated to non-stochastic thermal cycles had growth rates lower than those of control fish. Exposure to warm, fluctuating cycles did not affect emersion temperature, and only regular diel cycles modestly increased critical thermal tolerance. We predicted that warm diel cycling temperatures would increase gill surface area. Notably, fish acclimated to either thermal cycle had a reduced gill surface area and increased intralamellar cell mass when compared with control fish. This decreased gill surface area with warming contrasts with what is observed for exclusively aquatic fish and suggests a preparatory gill response for emersion in these amphibious fish. Collectively, our data reveal the importance of considering stochastic thermal variability when studying the effects of temperature on fishes.

Fish microbiota repel ovipositing mosquitoes.

The mere presence of predators causes prey organisms to display predation-avoidance strategies. Predator presence is often communicated through predator-released chemical signals. Ovipositing female mosquitoes of several species are repelled by unknown signals released from larvivorous fish. It was previously suggested that in many cases, a predator's microbiota plays an important role in the release of these signals; however, this mechanism is still poorly understood. In this study, we looked into the effects of the microbiota originating from the larvivorous Gambusia affinis (Baird and Girard) on the oviposition behaviour of gravid female mosquitoes. We used fish with altered microbiota and bacterial isolates in a set of outdoor mesocosm experiments to address this aim. We show that interference with the fish microbiota significantly reduces fish's repellent effect. We further show that the bacterium Pantoea pleuroti, isolated from the skin of the fish, repels oviposition of Culex laticinctus Edwards and Culiseta longiareolata Macquart mosquitoes similarly to the way in which live fish repel them. Our results highlight the importance of bacteria in the interspecies interactions of their hosts. Furthermore, this finding may lead to the development of an ecologically friendly mosquito repellent, that may reduce the use of larvivorous fish for mosquito control.

Habitat structural complexity predicts cognitive performance and behaviour in western mosquitofish.

Urban stream syndrome alters stream habitat complexity. We define habitat complexity as the degree of variation in physical habitat structure, with increasing variation equating to higher complexity. Habitat complexity affects species composition and shapes animal ecology, physiology, behaviour and cognition. We used a delayed detour test to measure whether cognitive processes (motor self-regulation) and behaviour (risk-taking) of female Western mosquitofish, Gambusia affinis, varied with habitat structural complexity (low, moderate and high) that was quantified visually for nine populations. We predicted that motor self-regulation and risk-taking behaviour would increase with increasing habitat complexity, yet we found support for the opposite. Lower complexity habitats offer less refuge potentially leading to higher predation pressure and selecting for greater risk-taking by fish with higher motor self-regulation. Our findings provide insight into how habitat complexity can shape cognitive processes and behaviour and offers a broader understanding of why some species may tolerate conditions of urbanized environments.

Pharmaceutical Pollution Alters the Structure of Freshwater Communities and Hinders Their Recovery from a Fish Predator.

Freshwater ecosystems are under threat from rising pharmaceutical pollution. While such pollutants are known to elicit biological effects on organisms, we have limited knowledge on how these effects might cascade through food-webs, disrupt ecological processes, and shape freshwater communities. In this study, we used a mesocosm experiment to explore how the community impacts of a top-order predator, the eastern mosquitofish (Gambusia holbrooki), are mediated by exposure to environmentally relevant low (measured concentration: ∼10 ng/L) and high concentrations (∼110 ng/L) of the pervasive pharmaceutical pollutant fluoxetine. We found no evidence that exposure to fluoxetine altered the consumptive effects of mosquitofish on zooplankton. However, once mosquitofish were removed from the mesocosms, zooplankton abundance recovered to a greater extent in control mesocosms compared to both low and high fluoxetine-exposed mesocosms. By the end of the experiment, this resulted in fundamental differences in community structure between the control and fluoxetine-treated mesocosms. Specifically, the control mesocosms were characterized by higher zooplankton abundances and lower algal biomass, whereas mesocosms exposed to either low or high concentrations of fluoxetine had lower zooplankton abundances and higher algal biomass. Our results suggest that fluoxetine, even at very low concentrations, can alter aquatic communities and hinder their recovery from disturbances.

Different patterns of coupled predator-prey dynamics when the same species interact in different locations.

Individual predator and prey species exhibit coupled population dynamics in simple laboratory systems and simple natural communities. It is unclear how often such pairwise coupling occurs in more complex communities, in which an individual predator species might feed on several prey species and an individual prey species might be attacked by several predators. To examine this problem, we applied multivariate autoregressive state-space (MARSS) models to 5-year time-series of monthly surveys of a predatory fish, the eastern mosquitofish (Gambusia holbrooki), and its littoral zone prey species, the least killifish (Heterandria formosa), in three locations in north Florida. The MARSS models were consistent with coupled predator-prey dynamics at two of the three locations. In one of these two locations, the estimated densities of the two species displayed classic predator-prey oscillations. In the third location, there was a positive effect of killifish density on mosquitofish density but no detectable effect of mosquitofish density on killifish density. In all three locations, increased submergent vegetation cover was associated with increased prey density but not increased predator density. Eigenvalues analyses for the joint predator-prey dynamics indicated that one of the cyclic locations had more stable dynamics than the other locations. The three different patterns demonstrate that the dynamics of a pairwise predator-prey interaction emerge not only from the characteristics of the prey and the predator, but also those of the habitat and trophic web in which the predator and prey are embedded.

Cortisol enhances aerobic metabolism and locomotor performance during the transition to land in an amphibious fish.

Amphibious fishes on land encounter higher oxygen (O2) availability and novel energetic demands, which impacts metabolism. Previous work on the amphibious mangrove killifish (Kryptolebias marmoratus) has shown that cortisol becomes elevated in response to air exposure, suggesting a possible role in regulating metabolism as fish move into terrestrial environments. We tested the hypothesis that cortisol is the mechanism by which oxidative processes are upregulated during the transition to land in amphibious fishes. We used two groups of fish, treated fish (+metyrapone, a cortisol synthesis inhibitor) and control (-metyrapone), to determine the impact of cortisol during air exposure (0 and 1 h, 7 days) on O2 consumption, terrestrial locomotion, the phenotype of red skeletal muscle, and muscle lipid concentration. Metyrapone-treated fish had an attenuated elevation in O2 consumption rate during the water to air transition and an immediate reduction in terrestrial exercise performance relative to control fish. In contrast, we found no short- (0 h) or long-term (7 days) differences between treatments in the oxidative phenotype of red muscles, nor in muscle lipid concentrations. Our results suggest that cortisol stimulates the necessary increase in aerobic metabolism needed to fuel the physiological changes that amphibious fishes undergo during the acclimation to air, although further studies are required to determine specific mechanisms of cortisol regulation.

Hormonal dynamics of matrotrophy vs. lecithotrophy in live-bearing fish reproduction.

We explored the relationship between gestational states, fecundity, and steroid hormone levels in three species of live-bearing fish with different maternal provisioning strategies. We studied two lecithotrophic species, Gambusia affinis and Xiphophorus couchianus, where embryos feed exclusively on yolk stored in the eggs, and one matrotrophic species, Heterandria formosa, which actively transfers nutrients to embryos through a follicular placenta. We measured water-borne cortisol, estradiol, and progesterone along with brood size (fecundity) and gestational stage(s). We examined the physiological costs of both maternal provisioning modes. Matrotrophy likely imposes energetic demands due to active nutrient transfer, while lecithotrophy may incur costs from carrying many large embryos. We hypothesized that fecundity, gestational stage, and hormones would covary differently in lecithotrophic vs. matrotrophic species. We found no relationships between hormones and fecundity or gestational stage in any species. However, in H. formosa, we found a positive relationship between estradiol levels and female mass, and a negative relationship between progesterone levels and female mass indicating a change in the circulating levels of both hormones as females grow. We observed differences in average hormone levels among species: the matrotrophic species had higher progesterone and lower estradiol compared to lecithotrophic species. Higher estradiol in lecithotrophic species may relate to egg yolk formation, while placental structures could play a role in progesterone production in matrotrophic species. Elevated cortisol in H. formosa suggests either higher energetic costs or a preparative role for reproduction. Our findings highlight progesterone's importance in maintaining gestation in matrotrophic species, like other placental species.

Season, size, and sex: factors influencing monogenean prevalence and intensity on Gambusia affinis in New Zealand.

A number of studies have been conducted on monogenean seasonality, though primarily in continental regions with wide annual temperatures ranges. We investigated seasonal changes in the prevalence and intensity of Salsuginus seculus infesting sexually dimorphic western mosquitofish (Gambusia affinis) in New Zealand. This represents the first examination of seasonality for this species globally, and the first seasonal assessment of any monogenean population in New Zealand, a temperate country with a mild oceanic climate. Prevalence and intensity of S. seculus with respect to fish size and sex was also examined. Prevalence of S. seculus changed temporally, peaking in summer, and was strongly positively correlated with algal concentrations. This relationship may be associated with increasing food levels, leading to an increase in fish courting and mating, resulting in high numbers and close physical associations of G. affinis individuals, facilitating transmission of the monogeneans. Thus, biotic factors may be important in determining temporal changes in S. seculus prevalence in New Zealand. Female G. affinis had a significantly higher prevalence and mean intensity of S. seculus than males. Longer fish had a higher mean intensity and prevalence of S. seculus. Female G. affinis likely host disproportionately more monogeneans as they are larger than males. Alternatively, females may have a compromised immune response during reproductive periods. Overall, seasonal change was observed in S. seculus prevalence and intensity under New Zealand's mild climatic conditions, and the larger female G. affinis in this dimorphic species supported a greater prevalence and intensity of infestation than males.

Massive expansion of sex-specific SNPs, transposon-related elements, and neocentromere formation shape the young W-chromosome from the mosquitofish Gambusia affinis.

BACKGROUND: The Western mosquitofish, Gambusia affinis, is a model for sex chromosome organization and evolution of female heterogamety. We previously identified a G. affinis female-specific marker, orthologous to the aminomethyl transferase (amt) gene of the related platyfish (Xiphophorus maculatus). Here, we have analyzed the structure and differentiation of the G. affinis W-chromosome, using a cytogenomics and bioinformatics approach. RESULTS: The long arm of the G. affinis W-chromosome (Wq) is highly enriched in dispersed repetitive sequences, but neither heterochromatic nor epigenetically silenced by hypermethylation. In line with this, Wq sequences are highly transcribed, including an active nucleolus organizing region (NOR). Female-specific SNPs and evolutionary young transposable elements were highly enriched and dispersed along the W-chromosome long arm, suggesting constrained recombination. Wq copy number expanded elements also include female-specific transcribed sequences from the amt locus with homology to TE. Collectively, the G. affinis W-chromosome is actively differentiating by sex-specific copy number expansion of transcribed TE-related elements, but not (yet) by extensive sequence divergence or gene decay. CONCLUSIONS: The G. affinis W-chromosome exhibits characteristic genomic properties of an evolutionary young sex chromosome. Strikingly, the observed sex-specific changes in the genomic landscape are confined to the W long arm, which is separated from the rest of the W-chromosome by a neocentromere acquired during sex chromosome evolution and may thus have become functionally insulated. In contrast, W short arm sequences were apparently shielded from repeat-driven differentiation, retained Z-chromosome like genomic features, and may have preserved pseudo-autosomal properties.

Predictability and conceptual repeatability of the predator-associated burst speed ecophenotype in western mosquitofish (Gambusia affinis).

Predation exerts strong selection pressure on morphological traits and is often studied in freshwater fishes. A common morphological ecophenotype driven by predation from pursuit predators is the predator-associated burst speed (PABS) ecophenotype. This ecophenotype, characterized by a smaller head, smaller body, and larger caudal region, is commonly found in western mosquitofish (Gambusia affinis) residing in environments with sunfish (family Centrarchidae) predators. However, the repeatability and transferability of the PABS ecophenotype across populations have not been tested. The four objectives of this study were to (1) assess the repeatability of methods to confirm the presence of the PABS ecophenotype in independent populations, (2) test whether morphology is correlated with predator abundance, (3) assess the repeatability of utilizing a discriminant function analysis (DFA) to accurately classify fish to predator status, and (4) use a DFA to conduct a cross-validation test of the PABS model using previously studied populations of G. affinis to predict predator status of newly studied populations. There was consistency in the occurrence of the PABS ecophenotype among populations with predators. Permuted MANOVA (PMANOVA) models fit to each dataset revealed a significant effect of collection site nested within predator status (i.e., presence or absence of predator), suggesting location-specific variations of the PABS ecophenotype. In addition, a PMANOVA model fit to body shape as a function of predator abundance (i.e., 0, 0.1-1.0, or >1.0 predators per minute of electrofishing) revealed a significant effect of predation level in the newly studied populations, suggesting morphology differs among populations with varying predator abundances. Cross-study validation of the DFA revealed low between-study accuracy compared to within-study accuracy, but elongation of the caudal region in the presence of predators was consistent across studies. Our findings ultimately suggest that although the PABS ecophenotype at studied locations provides partial predictive capacity at unstudied locations, the nature and magnitude of the PABS ecophenotype depend on environmental settings, encounter histories with predators, level of abundance of predators, or other spatially structured mechanisms.

Shoals in troubled waters? The impact of rising temperatures on metabolism, foraging, and shoaling behavior in mixed-species shoals.

Rising water temperatures across aquatic habitats, in the current global climate change scenario, can directly affect metabolism and food intake in fish species. This can potentially alter their physiological, behavioral, and shoaling properties. In the current study, we examined the effects of high temperatures on metabolism, foraging, and shoaling in tropical fish. Mixed-species (comprising flying barbs, zebrafish, and gambusia) and single-species (flying barbs and zebrafish) shoals were conditioned for 45 days to three kinds of temperature regimes: the current temperature regime (CTR), in which shoals were maintained at water temperature of 24°C (i.e., the current mean temperature of their habitat), the predicted temperature regime (PTR) at 31°C (i.e., simulating conditions projected for their habitat in 2100), and the dynamic temperature regime (DTR), which experienced daily temperature fluctuations between 24 and 31°C (i.e., resembling rapid temperature changes expected in their natural environments). We found species-specific responses to these temperature regimes. Flying barbs exhibited significantly lower body weight at PTR but maintained consistent muscle glycogen content across all temperature regimes. In contrast, zebrafish and gambusia displayed significantly elevated muscle glycogen content at PTR, with similar body weights across all three temperature regimes. Cohesion within flying barb shoals and cohesion/polarization in mixed-species shoals decreased significantly at PTR. Shoals exposed to DTR exhibited intermediate characteristics between those conditioned to CTR and PTR, suggesting that shoals may be less impacted by dynamic temperatures compared to prolonged high temperatures. This study highlights species-specific metabolic responses to temperature changes and their potential implications for larger-scale shoal properties.