Trophic niche drives the evolution of craniofacial shape in Trinidadian guppies

Diverse clades of fishes adapted to feeding on the benthos repeatedly converge on steep craniofacial profiles and shorter, wider heads. But in an incipient radiation, to what extent is this morphological evolution measurable and can we distinguish the relative genetic vs. plastic effects? We use the Trinidadian guppy (Poecilia reticulata) to test the repeatability of adaptation and the alignment of genetic and environmental effects shaping poecilid craniofacial morphology. We compare wild-caught and common garden lab-reared fish to quantify the genetic and plastic components of craniofacial morphology across 4 populations from 2 river drainage systems (n = 56 total). We first use micro-computed tomography to capture 3D morphology, then place both landmarks and semilandmarks to perform size-corrected 3D morphometrics and quantify shape space. We find a measurable, significant, and repeatable divergence in craniofacial shape between high-predation invertivore and low-predation detritivore populations. As predicted from previous examples of piscine adaptive trophic divergence, we find increases in head slope and craniofacial compression among the benthic detritivore foragers. Furthermore, the effects of environmental plasticity among benthic detritivores produce exaggerated craniofacial morphological change along a parallel axis to genetic morphological adaptation from invertivore ancestors. Overall, many of the major patterns of benthic-limnetic craniofacial evolution appear convergent among disparate groups of teleost fishes.

Identification of Philaster apodigitiformis in aquaculture and functional characterization of its β-PKA gene and expression analysis of infected Poecilia reticulata

Cyclic adenosine monophosphate (cAMP)-dependent protein kinase A (PKA) is a distinctive member of the serine–threonine protein AGC kinase family and an effective kinase for cAMP signal transduction. In recent years, scuticociliate has caused a lot of losses in domestic fishery farming, therefore, we have carried out morphological and molecular biological studies. In this study, diseased guppies (Poecilia reticulata) were collected from an ornamental fish market, and scuticociliate Philaster apodigitiformis Miao et al., 2009 was isolated. In our prior transcriptome sequencing research, we discovered significant expression of the β-PKA gene in P. apodigitiformis during its infection process, leading us to speculate its involvement in pathogenesis. A complete sequence of the β-PKA gene was cloned, and quantified by quantitative reverse transcription-polymerase chain reaction to analyse or to evaluate the functional characteristics of the β-PKA gene. Morphological identification and phylogenetic analysis based on small subunit rRNA sequence, infection experiments and haematoxylin–eosin staining method were also carried out, in order to study the pathological characteristics and infection mechanism of scuticociliate. The present results showed that: (1) our results revealed that β-PKA is a crucial gene involved in P. apodigitiformis infection in guppies, and the findings provide valuable insights for future studies on scuticociliatosis; (2) we characterized a complete gene, β-PKA, that is generally expressed in parasitic organisms during infection stage and (3) the present study indicates that PKA plays a critical role in scuticociliate when infection occurs by controlling essential steps such as cell growth, development and regulating the activity of the sensory body structures and the irritability system.

Aesthetic preference for artificially selected color variant affects mate choice copying behavior in female Poecilia latipinna

Three experiments were conducted examining whether an artificially selected "gold" color variant in female "models" affects mate choice copying behavior in sailfin mollies (Poecilia latipinna). Experiment I consisted of a pair of female preference assays, first assessing preference for male body size, followed by a mate choice copying assay that paired a model female with the smaller, non-preferred male from the initial preference test. Female subjects were divided into three groups that used either a wildtype female model, an artificially selected "gold" variant (cultivated within the aquarium fish trade) model, or control wherein no model was presented. Results showed females consistently copied the model's choice, switching preferences from the larger to smaller male when paired with a model regardless of color. In the second experiment wildtype females were presented with a pair of size-matched dummy males both of which paired with model females (one gold and the other wild type). Subjects consistently preferred the male previously paired with the gold- over the male with the wildtype-model, suggesting pre-existing sensory/perceptual biases may have affected their mate choice copying behavior. Previous studies have offered evidence for the spread of novel traits in males via sensory exploitation. However, these results indicate such biases may influence courtship behavior in circumstances where the novel trait is expressed in females as well. For the third experiment, wildtype females were presented with a choice between gold vs wildtype dummy males, the results of which revealing significant preferences for gold. In a follow-up assay pairing a wild type model with the non-preferred wildtype male, females maintained their preference for gold males despite the conflicting social driver of mate choice copying. These data offer evidence for the existence of a perceptual/cognitive bias in the context of mate choice copying, favoring the gold phenotype and/or novelty in general.

Identification of Candidate Genes for Red-Eyed (Albinism) Domestic Guppies Using Genomic and Transcriptomic Analyses

Guppies are small tropical fish with brightly colored bodies and variable tail shapes. There are two phenotypes of domestic guppy eye color: red and black. The wild type is black-eyed. The main object of this study was to identify candidate genes for the red-eyed phenotype in domestic guppies. We hope to provide molecular genetic information for the development of new domestic guppy strains. Additionally, the results also contribute to basic research concerning guppies. In this study, 121 domestic guppies were used for genomic analysis (GWAS), and 44 genes were identified. Furthermore, 21 domestic guppies were used for transcriptomic analysis, and 874 differentially expressed genes (DEGs) were identified, including 357 upregulated and 517 downregulated genes. Through GO and KEGG enrichment, we identified some important terms or pathways mainly related to melanin biosynthesis and ion transport. qRT-PCR was also performed to verify the differential expression levels of four important candidate genes (TYR, OCA2, SLC45A2, and SLC24A5) between red-eyed and black-eyed guppies. Based on the results of genomic and transcriptomic analyses, we propose that OCA2 is the most important candidate gene for the red-eyed phenotype in guppies.

Identification of a Tetrahymena species infecting guppies, pathology, and expression of beta-tubulin during infection

Tetrahymenosis is caused by the ciliated protozoan Tetrahymena and is responsible for serious economic losses to the aquaculture industry worldwide. However, information regarding the molecular mechanism leading to tetrahymenosis is limited. In previous transcriptome sequencing work, it was found that one of the two β-tubulin genes in T. pyriformis was significantly expressed in infected fish, we speculated that β-tubulin is involved in T. pyriformis infecting fish. Herein, the potential biological function of the β-tubulin gene in Tetrahymena species when establishing infection in guppies was investigated by cloning the full-length cDNA of this T. pyriformis β-tubulin (BTU1) gene. The full-length cDNA of T. pyriformis BTU1 gene was 1873 bp, and the ORF occupied 1134 bp, whereas 5' UTR 434 bp, and 3' UTR 305 bp whose poly (A) tail contained 12 bases. The predicted protein encoded by T. pyriformis BTU1 gene had a calculated molecular weight of 42.26 kDa and pI of 4.48. Moreover, secondary structure analysis and tertiary structure prediction of BTU1 protein were also conducted. In addition, morphology, infraciliature, phylogeny, and histopathology of T. pyriformis isolated from guppies from a fish market in Harbin were also investigated. Furthermore, qRT-PCR analysis and experimental infection assays indicated that the expression of BTU1 gene resulted in efficient cell proliferation during infection. Collectively, our data revealed that BTU1 is a key gene involved in T. pyriformis infection in guppies, and the findings discussed herein provide valuable insights for future studies on tetrahymenosis.

Has recombination changed during the recent evolution of the guppy Y chromosome?

Genome sequencing and genetic mapping of molecular markers have demonstrated nearly complete Y-linkage across much of the guppy (Poecilia reticulata) XY chromosome pair. Predominant Y-linkage of factors controlling visible male-specific coloration traits also suggested that these polymorphisms are sexually antagonistic (SA). However, occasional exchanges with the X are detected, and recombination patterns also appear to differ between natural guppy populations, suggesting ongoing evolution of recombination suppression under selection created by partially sex-linked SA polymorphisms. We used molecular markers to directly estimate genetic maps in sires from 4 guppy populations. The maps are very similar, suggesting that their crossover patterns have not recently changed. Our maps are consistent with population genomic results showing that variants within the terminal 5 Mb of the 26.5 Mb sex chromosome, chromosome 12, are most clearly associated with the maleness factor, albeit incompletely. We also confirmed occasional crossovers proximal to the male-determining region, defining a second, rarely recombining, pseudo-autosomal region, PAR2. This fish species may therefore have no completely male-specific region (MSY) more extensive than the male-determining factor. The positions of the few crossover events suggest a location for the male-determining factor within a physically small repetitive region. A sex-reversed XX male had few crossovers in PAR2, suggesting that this region's low crossover rate depends on the phenotypic, not the genetic, sex. Thus, rare individuals whose phenotypic and genetic sexes differ, and/or occasional PAR2 crossovers in males can explain the failure to detect fully Y-linked variants.

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.

Estrogen regulates the transcription of guppy isotocin receptors

Estrogen can regulate oxytocin receptor expression, which is mediated through estrogen receptors (ESRs) in mammals, initiating parturition. To further study the reproductive physiological process of ovoviviparous teleosts, guppies (Poecilia reticulata) were employed as the research model in the present study to identify the transcriptional regulation of ESRs on isotocin receptors (itrs). Since guppy embryos develop inside the ovary, in the present study, the levels of itrs in the ovarian stroma of pregnant female guppies treated with estradiol (E2) in vitro were tested. E2 increased only itr2 mRNA levels 3 h post-treatment, with no variation in itr1 mRNA expression levels. In vivo, pregnant guppies were immersed in different concentrations of E2, significantly increasing the relative expression levels of itr1 and itr2 in the ovary. Moreover, based on dual-fluorescence in situ hybridization (ISH), both esrs and itrs mRNAs were localized in the same cells around the embryos in the ovary. To further investigate the regulation of itr transcription by estrogen, a luciferase reporter assay was performed, and the results demonstrated that E2 treatment could induce E2-dependent repression of luciferase activity in cells transfected with ESR1. However, overexpression of ESR2a or ESR2b caused a robust ligand-independent increase in itr2 promoter activity. Deletion analysis of the itr2 promoter indicated that there were novel potential ESR transcription factor-binding sites at -360 bp upstream of the 5' end of the itr2 promoter. Overall, our study provided novel results regarding the ESRs mediating the onset of parturition in ovoviviparous teleosts.

Genetic and context-specific effects on individual inhibitory control performance in the guppy (Poecilia reticulata)

Among-individual variation in cognitive traits, widely assumed to have evolved under adaptive processes, is increasingly being demonstrated across animal taxa. As variation among individuals is required for natural selection, characterizing individual differences and their heritability is important to understand how cognitive traits evolve. Here, we use a quantitative genetic study of wild-type guppies repeatedly exposed to a 'detour task' to test for genetic variance in the cognitive trait of inhibitory control. We also test for genotype-by-environment interactions (GxE) by testing related fish under alternative experimental treatments (transparent vs. semi-transparent barrier in the detour-task). We find among-individual variation in detour task performance, consistent with differences in inhibitory control. However, analysis of GxE reveals that heritable factors only contribute to performance variation in one treatment. This suggests that the adaptive evolutionary potential of inhibitory control (and/or other latent variables contributing to task performance) may be highly sensitive to environmental conditions. The presence of GxE also implies that the plastic response of detour task performance to treatment environment is genetically variable. Our results are consistent with a scenario where variation in individual inhibitory control stems from complex interactions between heritable and plastic components.

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

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

Sperm production is negatively associated with muscle and sperm telomere length in a species subjected to strong sperm competition

Life-history theory suggests that ageing is one of the costs of reproduction. Accordingly, a higher reproductive allocation is expected to increase the deterioration of both the somatic and the germinal lines through enhanced telomere attrition. In most species, males' reproductive allocation mainly regards traits that increase mating and fertilization success, that is sexually selected traits. In this study, we tested the hypothesis that a higher investment in sexually selected traits is associated with a reduced relative telomere length (RTL) in the guppy (Poecilia reticulata), an ectotherm species characterized by strong pre- and postcopulatory sexual selection. We first measured telomere length in both the soma and the sperm over guppies' lifespan to see whether there was any variation in telomere length associated with age. Second, we investigated whether a greater investment in pre- and postcopulatory sexually selected traits is linked to shorter telomere length in both the somatic and the sperm germinal lines, and in young and old males. We found that telomeres lengthened with age in the somatic tissue, but there was no age-dependent variation in telomere length in the sperm cells. Telomere length in guppies was significantly and negatively correlated with sperm production in both tissues and life stages considered in this study. Our findings indicate that telomere length in male guppies is strongly associated with their reproductive investment (sperm production), suggesting that a trade-off between reproduction and maintenance is occurring at each stage of males' life in this species.

Whole Genome Assembly and Annotation of Blackstripe Livebearer Poeciliopsis prolifica

The blackstripe livebearer Poeciliopsis prolifica is a live-bearing fish belonging to the family Poeciliidae with high level of postfertilization maternal investment (matrotrophy). This viviparous matrotrophic species has evolved a structure similarly to the mammalian placenta. Placentas have independently evolved multiple times in Poeciliidae from nonplacental ancestors, which provide an opportunity to study the placental evolution. However, there is a lack of high-quality reference genomes for the placental species in Poeciliidae. In this study, we present a 674 Mb assembly of P. prolifica in 504 contigs with excellent continuity (contig N50 7.7 Mb) and completeness (97.2% Benchmarking Universal Single-Copy Orthologs [BUSCO] completeness score, including 92.6% single-copy and 4.6% duplicated BUSCO score). A total of 27,227 protein-coding genes were annotated from the merged datasets based on bioinformatic prediction, RNA sequencing and homology evidence. Phylogenomic analyses revealed that P. prolifica diverged from the guppy (Poecilia reticulata) ∼19 Ma. Our research provides the necessary resources and the genomic toolkit for investigating the genetic underpinning of placentation.

Selection on standing genetic variation mediates convergent evolution in extremophile fish

Hydrogen sulfide is a toxic gas that disrupts numerous biological processes, including energy production in the mitochondria, yet fish in the Poecilia mexicana species complex have independently evolved sulfide tolerance several times. Despite clear evidence for convergence at the phenotypic level in these fishes, it is unclear if the repeated evolution of hydrogen sulfide tolerance is the result of similar genomic changes. To address this gap, we used a targeted capture approach to sequence genes associated with sulfide processes and toxicity from five sulfidic and five nonsulfidic populations in the species complex. By comparing sequence variation in candidate genes to a reference set, we identified similar population structure and differentiation, suggesting that patterns of variation in most genes associated with sulfide processes and toxicity are due to demographic history and not selection. But the presence of tree discordance for a subset of genes suggests that several loci are evolving divergently between ecotypes. We identified two differentiation outlier genes that are associated with sulfide detoxification in the mitochondria that have signatures of selection in all five sulfidic populations. Further investigation into these regions identified long, shared haplotypes among sulfidic populations. Together, these results reveal that selection on standing genetic variation in putatively adaptive genes may be driving phenotypic convergence in this species complex.

The quest for good genes: Epigamic traits, fitness, MHC and multilocus heterozygosity in the guppy

The 'good genes' hypothesis for the evolution of male secondary sexual traits poses that female preferences for such traits are driven by indirect genetic benefits. However, support for the hypothesis remains ambiguous, and, in particular, the genetic basis for the benefits has rarely been investigated. Here, we use seminatural populations of Trinidadian guppies to investigate whether sexually selected traits (orange, black and iridescent colouration, gonopodium length and body size) predict fitness measured as the number of grandoffspring, a metric that integrates across fitness components and sexes. Furthermore, we tested whether two potential sources of genetic benefits-major histocompatibility complex (MHC) genotypes and multilocus heterozygosity (MLH)-are significant predictors of fitness and of the size of sexually selected traits. We found a significant, nonlinear effect of the area of black pigmentation and male body size on the number of grandoffspring, suggesting stabilizing selection on black area, and nonlinear selection favouring small body size. MLH was heritable (h2  = 0.14) and significantly predicted the number of grandoffspring, indicating the potential for genetic benefits based on heterozygosity. We also found support for local heterozygosity effects, which may reflect a noneven distribution of genetic load across the genome. MHC genotype was not significantly associated with any tested fitness component, or with the load of Gyrodactylus parasites. Neither MHC nor MLH was significant predictor of sexually selected traits. Overall, our results highlight the role of heterozygosity in determining fitness, but do not provide support for male sexually selected traits being indicators of genetic quality.

Extensive variation in germline de novo mutations in Poecilia reticulata

The rate of germline mutation is fundamental to evolutionary processes, as it generates the variation upon which selection acts. The guppy, Poecilia reticulata, is a model of rapid adaptation, however the relative contribution of standing genetic variation versus de novo mutation (DNM) to evolution in this species remains unclear. Here, we use pedigree-based approaches to quantify and characterize germline DNMs in three large guppy families. Our results suggest germline mutation rate in the guppy varies substantially across individuals and families. Most DNMs are shared across multiple siblings, suggesting they arose during early embryonic development. DNMs are randomly distributed throughout the genome, and male-biased mutation rate is low, as would be expected from the short guppy generation time. Overall, our study shows remarkable variation in germline mutation rate and provides insights into rapid evolution of guppies.

Sex chromosome heteromorphism and the Fast-X effect in poeciliids

Fast-X evolution has been observed in a range of heteromorphic sex chromosomes. However, it remains unclear how early in the process of sex chromosome differentiation the Fast-X effect becomes detectible. Recently, we uncovered an extreme variation in sex chromosome heteromorphism across poeciliid fish species. The common guppy, Poecilia reticulata, Endler's guppy, P. wingei, swamp guppy, P. picta and para guppy, P. parae, appear to share the same XY system and exhibit a remarkable range of heteromorphism. Species outside this group lack this sex chromosome system. We combined analyses of sequence divergence and polymorphism data across poeciliids to investigate X chromosome evolution as a function of hemizygosity and reveal the causes for Fast-X effects. Consistent with the extent of Y degeneration in each species, we detect higher rates of divergence on the X relative to autosomes, a signal of Fast-X evolution, in P. picta and P. parae, species with high levels of X hemizygosity in males. In P. reticulata, which exhibits largely homomorphic sex chromosomes and little evidence of hemizygosity, we observe no change in the rate of evolution of X-linked relative to autosomal genes. In P. wingei, the species with intermediate sex chromosome differentiation, we see an increase in the rate of nonsynonymous substitutions on the older stratum of divergence only. We also use our comparative approach to test for the time of origin of the sex chromosomes in this clade. Taken together, our study reveals an important role of hemizygosity in Fast-X evolution.

Evolutionary History of the Poecilia picta Sex Chromosomes

The degree of divergence between the sex chromosomes is not always proportional to their age. In poeciliids, four closely related species all exhibit a male heterogametic sex chromosome system on the same linkage group, yet show a remarkable diversity in X and Y divergence. In Poecilia reticulata and P. wingei, the sex chromosomes remain homomorphic, yet P. picta and P. parae have a highly degraded Y chromosome. To test alternative theories about the origin of their sex chromosomes, we used a combination of pedigrees and RNA-seq data from P. picta families in conjunction with DNA-seq data collected from P. reticulata, P. wingei, P. parae, and P. picta. Phylogenetic clustering analysis of X and Y orthologs, identified through segregation patterns, and their orthologous sequences in closely related species demonstrates a similar time of origin for both the P. picta and P. reticulata sex chromosomes. We next used k-mer analysis to identify shared ancestral Y sequence across all four species, suggesting a single origin to the sex chromosome system in this group. Together, our results provide key insights into the origin and evolution of the poeciliid Y chromosome and illustrate that the rate of sex chromosome divergence is often highly heterogenous, even over relatively short evolutionary time frames.

Effects of 4-n-nonylphenol in liver of male and female viviparous fish (Poecilia vivipara)

4-n-Nonylphenol (NP) is one of the most toxic alkylphenols found in the environment. To evaluate the transcriptional effects of NP in the viviparous fish Poecilia vivipara, a hepatic transcriptome and qPCR analysis of genes were carried out. Guppies separated by sex were injected with two doses of NP (15 μg/g and 150 μg/g) or peanut oil (control). After 24 h, analysis of transcriptional level of Aryl Hydrocarbon Receptor (AhR), Estrogen Nuclear Receptor Alpha (ESR1), Pregnane X Receptor (PXR), Cytochromes P450 (CYP1A, CYP2K1 and CYP3A30), Glutathione S-transferase A3 and Mu 3 (GSTa3 and GSTMu3), SRY-Box Transcription Factor 9 (SOX9), Vitellogenin-1 (VIT), ATP Binding Cassette Subfamily C Member 1 (ABCC1), Multidrug Resistance-Associated Protein 2 (MRP2) and UDP Glucuronosyltransferase Family 1 Member A1 (UGT1A1) was evaluated. 205,046 transcripts were assembled and protein prediction resulted in 203,147 predicted peptides. In females, no significant changes were detected in the transcription of some phase I biotransformation and ABC transporter genes. AhR, PXR, GSTa3 and SOX9 genes where higher in the lower dose group (15 μg/g) compared to control. In male fish, no changes were observed in the transcript levels of the nuclear receptors, in endocrine disruption and phase I biotransformation genes. GSTa3 showed lower transcription in fish treated with both doses. ABCC1 was higher in guppies treated with the lower dose while MRP2 showed less transcripts. This short-term and low-dose exposure to NP caused changes that could serve as early indicators of deleterious processes. These results indicate P. vivipara as a good sentinel in biomonitoring programs.

Triphenyltin exposure induced abnormal morphological colouration in adult male guppies (Poecilia reticulata)

Fish morphological colouration is essential for their survival and reproduction success; however, it is vulnerable to environmental factors, such as pollutants. Triphenyltin (TPT) is widespread in aquatic ecosystems, and its impacts on fish have been problematic. Therefore, the purpose of this study was to investigate the effects of TPT at environment-related concentrations (0, 1, 10 and 100 ng Sn/L) on morphological colouration in male guppies (Poecilia reticulata). The results showed that TPT exposure affected both orange/red and dark morphological colouration in guppies. The faded orange/red colouration might be related to the decrease of coloured pteridine and Pts (6-Pyruvoyltetrahydropterin Synthase) expression. In addition, TPT exposure induced melanogenesis, however, much melanin was distributed diffusely in the skin and did not seem to form a spot pattern, giving the fish a dull appearance. According to the skin transcriptional profiles, the changes of dark morphological colouration might be related to the changes in genes related to the functions of melanosome components (Gpnmb, Slc45a2 and Tyr), construction (Ap3d1, Fig4, Hps3, Hps5, Lyst, Rabggta, Txndc5 and Vps33a), and transport (Rab27a). Additionally, genes related to the regulation of melanogenesis (Atrn and Pomc) and system effects (Atox1, Atp6ap2, Atp6v1f, Atp6v1h, Rpl24, Rps19 and Rps20) might also be involved in the molecular mechanisms of abnormal morphological colouration induced by TPT. The present study provides crucial data on the molecular basis of abnormal morphological colouration in fish exposed to TPT and underscores the importance of toxicological studies of the effects of pollutants in aquatic environments on fish morphological colouration.

Mechanistic revealing of reproductive behavior impairment in male guppy (Poecilia reticulata) induced by environmentally realistic 2,2'-dithiobis-pyridine exposure

Although (PS)₂, the primary degradation product of emerging antifouling biocides metal pyrithiones (MePTs), can disrupt the reproductive behavior of fish at an environmentally relevant ng/L level, the underlying mechanism is still largely unknown. This study exposed sexually mature male guppy (Poecilia reticulata) to 20, 200, and 2000 ng/L (PS)₂ to explore the compromised effect of (PS)₂ on reproductive behavior through a realistic competing scenario. The results showed that (PS)₂ suppressed male guppies' sexual interest to stimulus females, reduced their competitive behavior frequencies toward rival males, and decreased their mating time and frequency. (PS)₂ exposure did not affect male guppies' secondary sexual characteristics or induce estrogenic activity. Whole-brain transcriptome sequencing identified 1070 differentially expressed genes (DEGs) with 872 up-regulated genes, which were functionally enriched into Gene Ontology terms pertaining to extracellular matrix (ECM) and extracellular region. KEGG enrichment for the DEGs uncovered that the activations of ECM-receptor interaction and focal adhesion pathways could be the underlying molecular mechanism implicated in the (PS)₂ induced reproductive behavior impairment. This work would deliver a substantial contribution to the understanding of the ecological safety of MePTs biocides.

How did the guppy Y chromosome evolve?

The sex chromosome pairs of many species do not undergo genetic recombination, unlike the autosomes. It has been proposed that the suppressed recombination results from natural selection favouring close linkage between sex-determining genes and mutations on this chromosome with advantages in one sex, but disadvantages in the other (these are called sexually antagonistic mutations). No example of such selection leading to suppressed recombination has been described, but populations of the guppy display sexually antagonistic mutations (affecting male coloration), and would be expected to evolve suppressed recombination. In extant close relatives of the guppy, the Y chromosomes have suppressed recombination, and have lost all the genes present on the X (this is called genetic degeneration). However, the guppy Y occasionally recombines with its X, despite carrying sexually antagonistic mutations. We describe evidence that a new Y evolved recently in the guppy, from an X chromosome like that in these relatives, replacing the old, degenerated Y, and explaining why the guppy pair still recombine. The male coloration factors probably arose after the new Y evolved, and have already evolved expression that is confined to males, a different way to avoid the conflict between the sexes.

We report new findings concerning the long-studied the guppy XY pair, which has remained somewhat mystifying. We show that it can be understood as a case of a recent sex chromosome turnover event in which an older, highly degenerated Y chromosome was lost, and creation of a new sex chromosome from the ancestral X. This chromosome acquired a male-determining factor, possibly by a mutation in (or a duplication of) a previously X-linked gene, or (less likely) by movement of an ancestral Y-linked maleness factor onto the X. We relate the findings to theoretical models of such events, and argue that the proposed change was free from factors thought to impede such turnovers. The change resulted in the intriguing situation where the X chromosome is old and the Y is much younger, and we discuss some other species where a similar change seems likely to have occurred.

Gill Transcriptome Analysis Revealed the Difference in Gene Expression Between Freshwater and Seawater Acclimated Guppy (Poecilia reticulata)

Guppy (Poecilia reticulata) can adapt to a wide range of salinity changes. To investigate the gene expression changes in the guppy exposed to seawater, we characterized its gill transcriptome using RNA sequencing. Experimental fish were exposed to salinity increase from 0 to 30‰ within 4 days, while control fish were cultured in freshwater (0‰ salinity). Seven days after salinity exposure, the gills were sampled and the mortality within 2 weeks was recorded. No significant difference in the cumulative mortality at the second week was found between the two groups. Transcriptomic analysis identified 3477 differentially expressed genes (DEGs), including 1067 upregulated and 2410 downregulated genes. These DEGs were enriched in several biological processes, including ion transport, ion homeostasis, ATP biosynthetic process, metabolic process, and immune system process. Oxidative phosphorylation was the most activated pathway. DEGs involved in the pathway "endoplasmic reticulum (ER)-mediated phagocytosis," "starch and sucrose metabolism," and "steroid biosynthesis" were mainly downregulated; chemokines and interleukins involved in "cytokine-cytokine receptor interaction" were differentially expressed. The present results suggested that oxidative phosphorylation had essential roles in osmoregulation in the gills of seawater acclimated guppy, during which the decline in the expression of genes encoding V-ATPases and calreticulin had a negative effect on the phagocytosis and immune response. Besides, several metabolic processes including "starch and sucrose metabolism" and "steroid biosynthesis" were affected. This study elucidates transcriptomic changes in osmotic regulation, metabolism, and immunity in seawater acclimated guppy.

Drainage-structuring of ancestral variation and a common functional pathway shape limited genomic convergence in natural high- and low-predation guppies

Studies of convergence in wild populations have been instrumental in understanding adaptation by providing strong evidence for natural selection. At the genetic level, we are beginning to appreciate that the re-use of the same genes in adaptation occurs through different mechanisms and can be constrained by underlying trait architectures and demographic characteristics of natural populations. Here, we explore these processes in naturally adapted high- (HP) and low-predation (LP) populations of the Trinidadian guppy, Poecilia reticulata. As a model for phenotypic change this system provided some of the earliest evidence of rapid and repeatable evolution in vertebrates; the genetic basis of which has yet to be studied at the whole-genome level. We collected whole-genome sequencing data from ten populations (176 individuals) representing five independent HP-LP river pairs across the three main drainages in Northern Trinidad. We evaluate population structure, uncovering several LP bottlenecks and variable between-river introgression that can lead to constraints on the sharing of adaptive variation between populations. Consequently, we found limited selection on common genes or loci across all drainages. Using a pathway type analysis, however, we find evidence of repeated selection on different genes involved in cadherin signaling. Finally, we found a large repeatedly selected haplotype on chromosome 20 in three rivers from the same drainage. Taken together, despite limited sharing of adaptive variation among rivers, we found evidence of convergent evolution associated with HP-LP environments in pathways across divergent drainages and at a previously unreported candidate haplotype within a drainage.

Shifts in bacterial communities and antibiotic resistance genes in surface water and gut microbiota of guppies (Poecilia reticulata) in the upper Rio Uberabinha, Brazil

Anthropogenic activities especially water pollution can affect the diversity and composition of microbial communities and promote the spread of antibiotic resistance genes (ARGs). In this study, water samples and guppies (Poecilia reticulata) were sampled from six sampling sites along the Uberabinha River in southeastern Brazil, both microbial communities and ARGs of surface waters and intestinal microbiota of guppies (Poecilia reticulata) were detected. According to the results of 16S rRNA amplicon sequencing, Proteobacteria, Bacteroidetes, Firmicutes and Actinobacteria were dominant phyla in both water and intestinal microbiota, but the abundance of putative pathogens was higher at heavily polluted sites. Up to 83% of bacteria in intestinal microbiota originated from water microbiota; this proportion was relatively higher in less polluted compared to polluted environments. ARGs providing resistance of tetracyclines and quinolones were dominant in both water and gut microbiota. The relative abundances of class I integrons and ARGs were as high as 1.74 × 10^-1/16S rRNA copies and 3.61 × 10^-1/16S rRNA copies, respectively, at heavily polluted sites. Correlation analysis suggests that integrons and bacteria play key roles in explaining the widespread occurrence of ARGs in the surface, but not in intestinal microbiota. We could rule out the class I integrons a potential intermediary bridge for ARGs between both types of microbiomes. Our results highlight the tight link in microbial communities and ARGs between ambient microbiota of stream ecosystems and intestinal microbiota of fish. Our study could have far-reaching consequences for fisheries and consumer safety and calls for investigations of gut microbiota of target species of both commercial fisheries and recreational (hobby) angling.

Kin Recognition in Guppies Uses Self-Referencing Based on Olfactory Cues

AbstractKin recognition plays an important role in social behavior and evolution, but the proximate mechanisms by which individuals recognize kin remain poorly understood. In many species, individuals form a "kin template" that they compare with conspecifics' phenotypes to assess phenotypic similarity-and, by association, relatedness. Individuals may form a kin template through self-inspection (i.e., self-referencing) and/or by observing their rearing associates (i.e., family referencing). However, despite much interest, few empirical studies have successfully disentangled self-referencing and family referencing. Here, we employ a novel set of breeding crosses using the Trinidadian guppy (Poecilia reticulata) to disentangle referencing systems by manipulating exposure to kin from conception onward. We show that guppies discriminate among their full and maternal half siblings, which can be explained only by self-referencing. Additional behavioral experiments revealed no evidence that guppies incorporate the phenotypes of their broodmates or mother into the kin template. Finally, by manipulating the format of our behavioral tests, we show that olfactory communication is both necessary and sufficient for kin discrimination. These results provide robust evidence that individuals recognize kin by comparing the olfactory phenotypes of conspecifics with their own. This study resolves key questions about the proximate mechanisms underpinning kin recognition, with implications for the ontogeny and evolution of social behavior.

Transcriptome-wide changes associated with the reproductive behaviour of male guppies exposed to 17α-ethinyl estradiol

Although many pharmaceutical compounds (and their metabolites) can induce harmful impacts at the molecular, physiological and behavioural levels, their underlying mechanistic associations have remained largely unexplored. Here, we utilized RNA-Seq to build a whole brain transcriptome profile to examine the impact of a common endocrine disrupting pharmaceutical (17α-ethinyl estradiol, EE2) on reproductive behaviour in wild guppies (Poecilia reticulata). Specifically, we annotated 16,791 coding transcripts in whole brain tissue in relation to the courtship behaviour (i.e. sigmoid display) of EE2 exposed (at environmentally relevant concentration of 8 ng/L for 28-days) and unexposed guppies. Further, we obtained 10,960 assembled transcripts matching in the non-coding orthologous genomes. Behavioural responses were assessed using a standard mate choice experiment, which allowed us to disentangle chemical cues from visual cues. We found that a high proportion of the RNAseq reads aligned back to our de novo assembled transcriptome with 80.59% mapping rate. Behavioural experiments showed that when males were presented only with female visual cues, there was a significant interaction between male treatment and female treatment in the time spent in the preference zone. This is one of the first studies to show that transcriptome-wide changes are associated with the reproductive behaviour of fish: EE2 exposed male guppies that performed high levels of courtship had a gene profile that deviated the most from the other treatment groups, while both non-courting EE2 and control males had similar gene signatures. Using Gene Ontology pathway analysis, our study shows that EE2-exposed males had gene transcripts enriched for pathways associated with altered immunity, starvation, altered metabolism and spermatogenesis. Our study demonstrates that multiple gene networks orchestrate courting behaviour, emphasizing the importance of investigating impacts of pharmaceuticals on gene networks instead of single genes.

Expansion of frozen hybrids in the guppy ectoparasite, Gyrodactylus turnbulli

Hybridization is one of the major factors contributing to the emergence of highly successful parasites. Hybrid vigour can play an important role in this process, but subsequent rounds of recombination in the hybrid population may dilute its effects. Increased fitness of hybrids can, however, be frozen by asexual reproduction. Here, we identify invasion of a 'frozen hybrid' genotype in natural populations of Gyrodactylus turnbulli, a facultatively sexual ectoparasitic flatworm that causes significant damage to its fish host. We resequenced genomes of these parasites infecting guppies from six Trinidad and Tobago populations, and found surprisingly high discrepancy in genome-wide nucleotide diversity between islands. The elevated heterozygosity on Tobago is maintained by predominantly clonal reproduction of hybrids formed from two diverged genomes. Hybridization has been followed by spread of the hybrids across the island, implying a selective advantage compared with native genotypes. Our results thus highlight that a single outcrossing event may be independently sufficient to cause pathogen expansion.

Divergence and Remarkable Diversity of the Y Chromosome in Guppies

The guppy sex chromosomes show an extraordinary diversity in divergence across populations and closely related species. In order to understand the dynamics of the guppy Y chromosome, we used linked-read sequencing to assess Y chromosome evolution and diversity across upstream and downstream population pairs that vary in predator and food abundance in three replicate watersheds. Based on our population-specific genome assemblies, we first confirmed and extended earlier reports of two strata on the guppy sex chromosomes. Stratum I shows significant accumulation of male-specific sequence, consistent with Y divergence, and predates the colonization of Trinidad. In contrast, Stratum II shows divergence from the X, but no Y-specific sequence, and this divergence is greater in three replicate upstream populations compared with their downstream pair. Despite longstanding assumptions that sex chromosome recombination suppression is achieved through inversions, we find no evidence of inversions associated with either Stratum I or Stratum II. Instead, we observe a remarkable diversity in Y chromosome haplotypes within each population, even in the ancestral Stratum I. This diversity is likely due to gradual mechanisms of recombination suppression, which, unlike an inversion, allow for the maintenance of multiple haplotypes. In addition, we show that this Y diversity is dominated by low-frequency haplotypes segregating in the population, suggesting a link between haplotype diversity and female preference for rare Y-linked color variation. Our results reveal the complex interplay between recombination suppression and Y chromosome divergence at the earliest stages of sex chromosome divergence.

Environmental Change, If Unaccounted, Prevents Detection of Cryptic Evolution in a Wild Population

AbstractDetecting contemporary evolution requires demonstrating that genetic change has occurred. Mixed effects models allow estimation of quantitative genetic parameters and are widely used to study evolution in wild populations. However, predictions of evolution based on these parameters frequently fail to match observations. Here, we applied three commonly used quantitative genetic approaches to predict the evolution of size at maturity in a wild population of Trinidadian guppies. Crucially, we tested our predictions against evolutionary change observed in common-garden experiments performed on samples from the same population. We show that standard quantitative genetic models underestimated or failed to detect the cryptic evolution of this trait as demonstrated by the common-garden experiments. The models failed because (1) size at maturity and fitness both decreased with increases in population density, (2) offspring experienced higher population densities than their parents, and (3) selection on size was strongest at high densities. When we accounted for environmental change, predictions better matched observations in the common-garden experiments, although substantial uncertainty remained. Our results demonstrate that predictions of evolution are unreliable if environmental change is not appropriately captured in models.

Co-exposure of iron oxide nanoparticles and glyphosate-based herbicide induces DNA damage and mutagenic effects in the guppy (Poecilia reticulata)

Iron oxide nanoparticles (IONPs) have been tested to remediate aquatic environments polluted by chemicals, such as pesticides. However, their interactive effects on aquatic organisms remain unknown. This study aimed to investigate the genotoxicity and mutagenicity of co-exposure of IONPs (γ-Fe2O3 NPs) and glyphosate-based herbicide (GBH) in the fish Poecilia reticulata. Thus, fish were exposed to citrate-functionalized γ-Fe2O3 NPs (0.3 mg L-1; 5.44 nm) alone or co-exposed to γ-Fe2O3 NPs (0.3 mg L-1) and GBH (65 and 130 μg of glyphosate L-1) during 14 and 21 days. The genotoxicity (DNA damage) was analyzed by comet assay, while the mutagenicity evaluated by micronucleus test (MN test) and erythrocyte nuclear abnormalities (ENA) frequency. The co-exposure induced clastogenic (DNA damage) and aneugenic (nuclear alterations) effects on guppies in a time-dependent pattern. Fish co-exposed to NPs and GBH (130 μg glyphosate L-1) showed high DNA damage when compared to NPs alone and control group, indicating synergic effects after 21 days of exposure. However, mutagenic effects (ENA) were observed in the exposure groups after 14 and 21 days. Results showed the potential genotoxic and mutagenic effects of maghemite NPs and GBH co-exposure to freshwater fish. The transformation and interaction of iron oxide nanoparticles with other pollutants, as herbicides, in the aquatic systems are critical factors in the environmental risk assessment of metal-based NPs.

Locating the Sex Determining Region of Linkage Group 12 of Guppy (Poecilia reticulata)

Despite over 100 years of study, the location of the fully sex-linked region of the guppy (Poecilia reticulata) carrying the male-determining locus, and the regions where the XY pair recombine, remain unclear. Previous population genomics studies to determine these regions used small samples from recently bottlenecked captive populations, which increase the false positive rate of associations between individuals' sexes and SNPs. Using new data from multiple natural populations, we show that a recently proposed candidate for this species' male-determining gene is probably not completely sex-linked, leaving the maleness factor still unidentified. Variants in the chromosome 12 region carrying the candidate gene sometimes show linkage disequilibrium with the sex-determining factor, but no consistently male-specific variant has yet been found. Our genetic mapping with molecular markers spread across chromosome 12 confirms that this is the guppy XY pair. We describe two families with recombinants between the X and Y chromosomes, which confirm that the male-determining locus is in the region identified by all previous studies, near the terminal pseudo-autosomal region (PAR), which crosses over at a very high rate in males. We correct the PAR marker order, and assign two unplaced scaffolds to the PAR. We also detect a duplication, with one copy in the male-determining region, explaining signals of sex linkage in a more proximal region.

Convergent evolution of conserved mitochondrial pathways underlies repeated adaptation to extreme environments

Extreme environments test the limits of life; yet, some organisms thrive in harsh conditions. Extremophile lineages inspire questions about how organisms can tolerate physiochemical stressors and whether the repeated colonization of extreme environments is facilitated by predictable and repeatable evolutionary innovations. We identified the mechanistic basis underlying convergent evolution of tolerance to hydrogen sulfide (H2S)-a toxicant that impairs mitochondrial function-across evolutionarily independent lineages of a fish (Poecilia mexicana, Poeciliidae) from H2S-rich springs. Using comparative biochemical and physiological analyses, we found that mitochondrial function is maintained in the presence of H2S in sulfide spring P. mexicana but not ancestral lineages from nonsulfidic habitats due to convergent adaptations in the primary toxicity target and a major detoxification enzyme. Genome-wide local ancestry analyses indicated that convergent evolution of increased H2S tolerance in different populations is likely caused by a combination of selection on standing genetic variation and de novo mutations. On a macroevolutionary scale, H2S tolerance in 10 independent lineages of sulfide spring fishes across multiple genera of Poeciliidae is correlated with the convergent modification and expression changes in genes associated with H2S toxicity and detoxification. Our results demonstrate that the modification of highly conserved physiological pathways associated with essential mitochondrial processes mediates tolerance to physiochemical stress. In addition, the same pathways, genes, and-in some instances-codons are implicated in H2S adaptation in lineages that span 40 million years of evolution.

Transcriptional effects in the estuarine guppy Poecilia vivipara exposed to sanitary sewage in laboratory and in situ

The urban growth has increased sanitary sewage discharges in coastal ecosystems, negatively affecting the aquatic biota. Mangroves, one of the most human-affected coastal biomes, are areas for reproduction and nursing of several species. In order to evaluate the effects of sanitary sewage effluents in mangrove species, this study assessed the hepatic transcriptional responses of guppy fish Poecilia vivipara exposed to sanitary sewage 33% (v:v), using suppressive subtraction hybridization (SSH), high throughput sequencing of RNA (Ion-proton) and quantification of transcript levels by qPCR of some identified genes in fish kept in a sewage-contaminated environment. Genes identified are related predominantly to xenobiotic biotransformation, immune system and sexual differentiation. The qPCR results confirmed the induction of cytochrome P450 1A (CYP1A), glutathione S transferase A-like (GST A-like) methyltransferase (MET) and UDP glycosyltransferase 1A (UDPGT1A), and repression of complement component C3 (C3), doublesex and mab-3 related transcription factor 1 (DMRT1), and transferrin (TF) in the laboratory experiment. In the field exposure, the transcript levels of CYP1A, DMRT1, MET, GST A-like and UDPGT1A were higher in fishes exposed at the contaminated sites compared to the reference site. Chemical analysis in fish from the laboratory and in situ experiments, and surface sediment from the sewage-contaminated sites revealed relevant levels of polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyl (PCBs) and linear alkylbenzenes (LABs). These data reinforce the use of P. vivipara as a sentinel for monitoring environmental contamination in coastal regions.

Artificial selection on brain size leads to matching changes in overall number of neurons

Neurons are the basic computational units of the brain, but brain size is the predominant surrogate measure of brain functional capacity in comparative and cognitive neuroscience. This approach is based on the assumption that larger brains harbor higher numbers of neurons and their connections, and therefore have a higher information-processing capacity. However, recent studies have shown that brain mass may be less strongly correlated with neuron counts than previously thought. Till now, no experimental test has been conducted to examine the relationship between evolutionary changes in brain size and the number of brain neurons. Here, we provide such a test by comparing neuron number in artificial selection lines of female guppies (Poecilia reticulata) with >15% difference in relative brain mass and numerous previously demonstrated cognitive differences. Using the isotropic fractionator, we demonstrate that large-brained females have a higher overall number of neurons than small-brained females, but similar neuronal densities. Importantly, this difference holds also for the telencephalon, a key region for cognition. Our study provides the first direct experimental evidence that selection for brain mass leads to matching changes in number of neurons and shows that brain size evolution is intimately linked to the evolution of neuron number and cognition.

Red Queen revisited: Immune gene diversity and parasite load in the asexual Poecilia formosa versus its sexual host species P. mexicana

In accordance with the Red Queen hypothesis, the lower genotypic diversity in clonally reproducing species should make them easier targets for pathogen infection, especially when closely related sexually reproducing species occur in close proximity. We analyzed two populations of clonal P. formosa and their sexual parental species P. mexicana by correlating individual parasite infection with overall and immune genotype. Our study revealed lower levels of overall genotypic diversity and marginally fewer MHC class I alleles in P. formosa individuals compared to sexually reproducing P. mexicana. Parasite load, however, differed only between field sites but not between species. We hypothesize that this might be due to slightly higher genotypic diversity in P. formosa at the innate immune system (toll like receptor 8) which is likely due to the species’ hybrid origin. In consequence, it appears that clonal individuals do not necessarily suffer a disadvantage compared to sexual individuals when fighting parasite infection.

Characterization and expression of Xiphophorus maculatus microsatellite Msb069 full sequence in subgenus Poecilia

Msb069 primer pairs encompassed region is believed to be associated with a quantitative trait loci (QTL) of dorsal fin length in subgenus Poecilia. However, detailed investigation on Msb069 which originated from Xiphophorus on subgenus Poecilia remains unexplored. In this study, full sequence of Msb069 was characterized by sequencing bioinformatics analysis and gene expression. The sequence analysis of Msb069 primer pairs encompassed region on three species of Poecilia revealed higher number of microsatellite tandem repeats in Poecilia latipinna (ATG₁₆) compared to P. sphenops (ATG_13-14). There is no notable pattern of ATGtandem repeats discovered in the hybrids. The full sequence of Msb069 is 734 bp in length and showed a 233 bp conserved region between Xiphophorus and Poecilia. BLAST search performed on this sequence revealed no significant similarities. Nonquantitative RT-PCR exhibited the presence of Msb069 transcripts in three different tissues in subgenus Poecilia. Meanwhile, quantitative RTPCR expression on two different tissues showed relatively higher expression of Msb069 transcript in P. latipinna dorsal fin tissues in both male and female fishes, suggesting a repressive function of this transcript with respect to dorsal fin length. However the exact gene expression event of Msb069 is still unknown and requires further investigation.

Variation in expression and activity pattern of mmp2 and mmp9 on different time scales in the regenerating caudal fin of Poecilia latipinna

Alteration in the expression pattern of matrix metalloproteinase (MMP)2 and MMP9 was studied in the regenerating caudal fin of Poecilia latipinna immediately following amputation until the new tissues gained structural integrity. Timed expression pattern of these two MMPs was studied at enzyme, transcript as well as protein levels. Additionally, both the gelatinases were localized in the regenerating caudal fin during three specific stages of regeneration. The results revealed a progressive increase in the expression of MMP2 starting at 1 h post amputation (hpa), indicating its possible role in the remodelling of extracellular matrix early on during caudal-fin regeneration. Nevertheless, a reduction in transcript level expression of MMP2 at 6 hpa and 12 hpa stages, points towards a possible transcriptional regulation, to further moderate its activity. As observed in the case of MMP2, expression of MMP9 too increased from 1 hpa and remained elevated until 5 dpa. However, the active MMP9 revealed its presence only 12 hpa onwards. Moreover, both the gelatinases were localised in the apical epithelial cap and in the progress zone at wound epithelium (1 dpa) and blastema (60 hpa) stages respectively. Further, during early differentiation stage (5 dpa), high intensities of MMP2 and MMP9 were localized in the newly formed actinotrichia as compared with the tissue proximal to it. Based on the results, it could be construed that the controlled up-regulation of MMP2 and MMP9 from 1 hpa until the early differentiation stage ensures a regulated digestion of extracellular matrix, perhaps to facilitate the recruitment, proliferation, morphogenesis and re-patterning of resident stem cells during caudal fin regeneration in P. latipinna.

Clonal polymorphism and high heterozygosity in the celibate genome of the Amazon molly

The extreme rarity of asexual vertebrates in nature is generally explained by genomic decay due to absence of meiotic recombination, thus leading to extinction of such lineages. We explore features of a vertebrate asexual genome, the Amazon molly, Poecilia formosa, and find few signs of genetic degeneration but unique genetic variability and ongoing evolution. We uncovered a substantial clonal polymorphism and, as a conserved feature from its interspecific hybrid origin, a 10-fold higher heterozygosity than in the sexual parental species. These characteristics seem to be a principal reason for the unpredicted fitness of this asexual vertebrate. Our data suggest that asexual vertebrate lineages are scarce not because they are at a disadvantage, but because the genomic combinations required to bypass meiosis and to make up a functioning hybrid genome are rarely met in nature.

Immunogenetic novelty confers a selective advantage in host-pathogen coevolution

The major histocompatibility complex (MHC) is crucial to the adaptive immune response of vertebrates and is among the most polymorphic gene families known. Its high diversity is usually attributed to selection imposed by fast-evolving pathogens. Pathogens are thought to evolve to escape recognition by common immune alleles, and, hence, novel MHC alleles, introduced through mutation, recombination, or gene flow, are predicted to give hosts superior resistance. Although this theoretical prediction underpins host-pathogen "Red Queen" coevolution, it has not been demonstrated in the context of natural MHC diversity. Here, we experimentally tested whether novel MHC variants (both alleles and functional "supertypes") increased resistance of guppies (Poecilia reticulata) to a common ectoparasite (Gyrodactylus turnbulli). We used exposure-controlled infection trials with wild-sourced parasites, and Gyrodactylus-naïve host fish that were F2 descendants of crossed wild populations. Hosts carrying MHC variants (alleles or supertypes) that were new to a given parasite population experienced a 35-37% reduction in infection intensity, but the number of MHC variants carried by an individual, analogous to heterozygosity in single-locus systems, was not a significant predictor. Our results provide direct evidence of novel MHC variant advantage, confirming a fundamental mechanism underpinning the exceptional polymorphism of this gene family and highlighting the role of immunogenetic novelty in host-pathogen coevolution.

Characterisation of the transcriptome of male and female wild-type guppy brains with RNA-Seq and consequences of exposure to the pharmaceutical pollutant, 17α-ethinyl estradiol

Waterways are increasingly being contaminated by chemical compounds that can disrupt the endocrinology of organisms. One such compound is 17α-ethinyl estradiol (EE2), a synthetic estrogen used in the contraceptive pill. Despite considerable research interest in the effects of EE2 on reproduction and gene expression, surprisingly, only a few studies have capitalised on technologies, such as next-generation sequencing (NGS), to uncover the molecular pathways related to EE2 exposure. Accordingly, using high-throughput sequencing technologies, the aim of our study was to explore the effects of EE2 on brain transcriptome in wild-type male and female guppy (Poecilia reticulata). We conducted two sets of experiments, where fish were exposed to EE2 (measured concentrations: 8ng/L and 38ng/L) in a flow-through system for 21days. The effects on the brain transcriptome on both males and females were assessed using Illumina sequencing (MiSeq and HiSeq) platform followed by bioinformatics analysis (edgeR, DESeq2). Here, we report that exposure to EE2 caused both up- and downregulation of specific transcript abundances, and affected transcript abundance in a sex-specific manner. Specifically, we found 773 transcripts, of which 60 were male-specific, 61 female-specific and 285 treatment-specific. EE2 affected expression of 165 transcripts in males, with 88 downregulated and 77 upregulated, while in females, 120 transcripts were affected with 62 downregulated and 58 upregulated. Finally, RT-qPCR validation demonstrated that expression of transcripts related to transposable elements, neuroserpin and heat shock protein were significantly affected by EE2-exposure. Our study is the first to report brain transcriptome libraries for guppies exposed to EE2. Not only does our study provide a valuable resource, it offers insights into the mechanisms underlying the feminizing effects on the brains of organisms exposed to environmentally realistic concentrations of EE2.

Convergent recombination suppression suggests role of sexual selection in guppy sex chromosome formation

Sex chromosomes evolve once recombination is halted between a homologous pair of chromosomes. The dominant model of sex chromosome evolution posits that recombination is suppressed between emerging X and Y chromosomes in order to resolve sexual conflict. Here we test this model using whole genome and transcriptome resequencing data in the guppy, a model for sexual selection with many Y-linked colour traits. We show that although the nascent Y chromosome encompasses nearly half of the linkage group, there has been no perceptible degradation of Y chromosome gene content or activity. Using replicate wild populations with differing levels of sexually antagonistic selection for colour, we also show that sexual selection leads to greater expansion of the non-recombining region and increased Y chromosome divergence. These results provide empirical support for longstanding models of sex chromosome catalysis, and suggest an important role for sexual selection and sexual conflict in genome evolution.

The Genome of the Trinidadian Guppy, Poecilia reticulata, and Variation in the Guanapo Population

For over a century, the live bearing guppy, Poecilia reticulata, has been used to study sexual selection as well as local adaptation. Natural guppy populations differ in many traits that are of intuitively adaptive significance such as ornamentation, age at maturity, brood size and body shape. Water depth, light supply, food resources and predation regime shape these traits, and barrier waterfalls often separate contrasting environments in the same river. We have assembled and annotated the genome of an inbred single female from a high-predation site in the Guanapo drainage. The final assembly comprises 731.6 Mb with a scaffold N50 of 5.3 MB. Scaffolds were mapped to linkage groups, placing 95% of the genome assembly on the 22 autosomes and the X-chromosome. To investigate genetic variation in the population used for the genome assembly, we sequenced 10 wild caught male individuals. The identified 5 million SNPs correspond to an average nucleotide diversity (π) of 0.0025. The genome assembly and SNP map provide a rich resource for investigating adaptation to different predation regimes. In addition, comparisons with the genomes of other Poeciliid species, which differ greatly in mechanisms of sex determination and maternal resource allocation, as well as comparisons to other teleost genera can begin to reveal how live bearing evolved in teleost fish.

Kin Recognition in a Clonal Fish, Poecilia formosa

Relatedness strongly influences social behaviors in a wide variety of species. For most species, the highest typical degree of relatedness is between full siblings with 50% shared genes. However, this is poorly understood in species with unusually high relatedness between individuals: clonal organisms. Although there has been some investigation into clonal invertebrates and yeast, nothing is known about kin selection in clonal vertebrates. We show that a clonal fish, the Amazon molly (Poecilia formosa), can distinguish between different clonal lineages, associating with genetically identical, sister clones, and use multiple sensory modalities. Also, they scale their aggressive behaviors according to the relatedness to other females: they are more aggressive to non-related clones. Our results demonstrate that even in species with very small genetic differences between individuals, kin recognition can be adaptive. Their discriminatory abilities and regulation of costly behaviors provides a powerful example of natural selection in species with limited genetic diversity.

Can multilocus heterozygosity reveal inbreeding depression?

Inbred lines of fish have been widely exploited as model organisms to determine the effect of inbreeding, which is often closely related to fitness such as endurance and productivity compared to morphological traits. Until now, much is unknown about the effects of inbreeding to fish. In the present study, inbred lines of guppies were used to examine the inbreeding effect on morphological traits corresponding to genotype variation. Two strains, called AY and NA1, were selected from the closed culture system. Both strains showed different levels of inbreeding coefficients when compared to microsatellite markers. The AY strain was less inbred as compared to NA1 strain. However, correlation between the standard lengths with multilocus heterozygosity (MLH) at the individual level was observed in the AY, but not in the NA1 strain. This indicated that highly inbred animals have higher similarity in morphological traits as compared to less inbred ones. The inbreeding process showed the importance of heterozygosity, even in laboratory-reared animals. This experiment illustrated the effects of inbreeding towards morphological and genetic changes. ?

Sequence Evolution and Expression of the Androgen Receptor and Other Pathway-Related Genes in a Unisexual Fish, the Amazon Molly, Poecilia formosa, and Its Bisexual Ancestors

The all-female Amazon molly (Poecilia formosa) originated from a single hybridization of two bisexual ancestors, Atlantic molly (Poecilia mexicana) and sailfin molly (Poecilia latipinna). As a gynogenetic species, the Amazon molly needs to copulate with a heterospecific male, but the genetic information of the sperm-donor does not contribute to the next generation, as the sperm only acts as the trigger for the diploid eggs' embryogenesis. Here, we study the sequence evolution and gene expression of the duplicated genes coding for androgen receptors (ars) and other pathway-related genes, i.e., the estrogen receptors (ers) and cytochrome P450, family19, subfamily A, aromatase genes (cyp19as), in the Amazon molly, in comparison to its bisexual ancestors. Mollies possess-as most other teleost fish-two copies of the ar, er, and cyp19a genes, i.e., arα/arβ, erα/erβ1, and cyp19a1 (also referred as cyp19a1a)/cyp19a2 (also referred to as cyp19a1b), respectively. Non-synonymous single nucleotide polymorphisms (SNPs) among the ancestral bisexual species were generally predicted not to alter protein function. Some derived substitutions in the P. mexicana and one in P. formosa are predicted to impact protein function. We also describe the gene expression pattern of the ars and pathway-related genes in various tissues (i.e., brain, gill, and ovary) and provide SNP markers for allele specific expression research. As a general tendency, the levels of gene expression were lowest in gill and highest in ovarian tissues, while expression levels in the brain were intermediate in most cases. Expression levels in P. formosa were conserved where expression did not differ between the two bisexual ancestors. In those cases where gene expression levels significantly differed between the bisexual species, P. formosa expression was always comparable to the higher expression level among the two ancestors. Interestingly, erβ1 was expressed neither in brain nor in gill in the analyzed three molly species, which implies a more important role of erα in the estradiol synthesis pathway in these tissues. Furthermore, our data suggest that interactions of steroid-signaling pathway genes differ across tissues, in particular the interactions of ars and cyp19as.

Selection for brain size impairs innate, but not adaptive immune responses

Both the brain and the immune system are energetically demanding organs, and when natural selection favours increased investment into one, then the size or performance of the other should be reduced. While comparative analyses have attempted to test this potential evolutionary trade-off, the results remain inconclusive. To test this hypothesis, we compared the tissue graft rejection (an assay for measuring innate and acquired immune responses) in guppies (Poecilia reticulata) artificially selected for large and small relative brain size. Individual scales were transplanted between pairs of fish, creating reciprocal allografts, and the rejection reaction was scored over 8 days (before acquired immunity develops). Acquired immune responses were tested two weeks later, when the same pairs of fish received a second set of allografts and were scored again. Compared with large-brained animals, small-brained animals of both sexes mounted a significantly stronger rejection response to the first allograft. The rejection response to the second set of allografts did not differ between large- and small-brained fish. Our results show that selection for large brain size reduced innate immune responses to an allograft, which supports the hypothesis that there is a selective trade-off between investing into brain size and innate immunity.

Poecilia picta, a Close Relative to the Guppy, Exhibits Red Male Coloration Polymorphism: A System for Phylogenetic Comparisons

Studies on the evolution of female preference and male color polymorphism frequently focus on single species since traits and preferences are thought to co-evolve. The guppy, Poecilia reticulata, has long been a premier model for such studies because female preferences and orange coloration are well known to covary, especially in upstream/downstream pairs of populations. However, focused single species studies lack the explanatory power of the comparative method, which requires detailed knowledge of multiple species with known evolutionary relationships. Here we describe a red color polymorphism in Poecilia picta, a close relative to guppies. We show that this polymorphism is restricted to males and is maintained in natural populations of mainland South America. Using tests of female preference we show female P. picta are not more attracted to red males, despite preferences for red/orange in closely related species, such as P. reticulata and P. parae. Male color patterns in these closely related species are different from P. picta in that they occur in discrete patches and are frequently Y chromosome-linked. P. reticulata have an almost infinite number of male patterns, while P. parae males occur in discrete morphs. We show the red male polymorphism in P. picta extends continuously throughout the body and is not a Y-linked trait despite the theoretical prediction that sexually-selected characters should often be linked to the heterogametic sex chromosome. The presence/absence of red male coloration of P. picta described here makes this an ideal system for phylogenetic comparisons that could reveal the evolutionary forces maintaining mate choice and color polymorphisms in this speciose group.

The effect of brain size evolution on feeding propensity, digestive efficiency, and juvenile growth

One key hypothesis in the study of brain size evolution is the expensive tissue hypothesis; the idea that increased investment into the brain should be compensated by decreased investment into other costly organs, for instance the gut. Although the hypothesis is supported by both comparative and experimental evidence, little is known about the potential changes in energetic requirements or digestive traits following such evolutionary shifts in brain and gut size. Organisms may meet the greater metabolic requirements of larger brains despite smaller guts via increased food intake or better digestion. But increased investment in the brain may also hamper somatic growth. To test these hypotheses we here used guppy (Poecilia reticulata) brain size selection lines with a pronounced negative association between brain and gut size and investigated feeding propensity, digestive efficiency (DE), and juvenile growth rate. We did not find any difference in feeding propensity or DE between large- and small-brained individuals. Instead, we found that large-brained females had slower growth during the first 10 weeks after birth. Our study provides experimental support that investment into larger brains at the expense of gut tissue carries costs that are not necessarily compensated by a more efficient digestive system.

Population Genomics of the Euryhaline Teleost Poecilia latipinna

Global climate change and increases in sea levels will affect coastal marine communities. The conservation of these ecologically important areas will be a challenge because of their wide geographic distribution, ecological diversity and species richness. To address this problem, we need to better understand how the genetic variation of the species in these communities is distributed within local populations, among populations and between distant regions. In this study we apply genotyping by sequencing (GBS) and examine 955 SNPs to determine Sailfin molly (Poecilia latipinna) genetic diversity among three geographically close mangrove salt marsh flats in the Florida Keys compared to populations in southern and northern Florida. The questions we are asking are whether there is sufficient genetic variation among isolated estuarine fish within populations and whether there are significant divergences among populations. Additionally, we want to know if GBS approaches agree with previous studies using more traditional molecular approaches. We are able to identify large genetic diversity within each saltmarsh community (π ≈ 36%). Additionally, among the Florida Key populations and the mainland or between southern and northern Florida regions, there are significant differences in allele frequencies seen in population structure and evolutionary relationships among individuals. Surprisingly, even though the cumulative F_ST value using all 955 SNPs within the three Florida Key populations is small, there are 29 loci with significant F_ST values, and 11 of these were outliers suggestive of adaptive divergence. These data suggest that among the salt marsh flats surveyed here, there is significant genetic diversity within each population and small but significant differences among populations. Much of the genetic variation within and among populations found here with GBS is very similar to previous studies using allozymes and microsatellites. However, the meaningful difference between GBS and these previous measures of genetic diversity is the number of loci examined, which allows more precise delineations of population structure as well as facilitates identifying loci with excessive F_ST values that could indicate adaptive divergence.