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.

Theory predicts that the sexes can achieve greater fitness if loci with sexually antagonistic polymorphisms become linked to the sex determining loci, and this can favour the spread of reduced recombination around sex determining regions. Given that sex-linked regions are frequently repetitive and highly heterozygous, few complete Y chromosome assemblies are available to test these ideas. The guppy system (Poecilia reticulata) has long been invoked as an example of sex chromosome formation resulting from sexual conflict. Early genetics studies revealed that male colour patterning genes are mostly but not entirely Y-linked, and that X-linkage may be most common in low predation populations. More recent population genomic studies of guppies have reached varying conclusions about the size and placement of the Y-linked region. However, this previous work used a reference genome assembled from short-read sequences from a female guppy. Here, we present a new guppy reference genome assembly from a male, using long-read PacBio single-molecule real-time sequencing (SMRT) and chromosome contact information. Our new assembly sequences across repeat- and GC-rich regions and thus closes gaps and corrects mis-assemblies found in the short-read female-derived guppy genome. Using this improved reference genome, we then employed broad population sampling to detect sex differences across the genome. We identified two small regions that showed consistent male-specific signals. Moreover, our results help reconcile the contradictory conclusions put forth by past population genomic studies of the guppy sex chromosome. Our results are consistent with a small Y-specific region and rare recombination in male guppies.

MEGASAT is software that enables genotyping of microsatellite loci using next-generation sequencing data. Microsatellites are amplified in large multiplexes, and then sequenced in pooled amplicons. MEGASAT reads sequence files and automatically scores microsatellite genotypes. It uses fuzzy matches to allow for sequencing errors and applies decision rules to account for amplification artefacts, including nontarget amplification products, replication slippage during PCR (amplification stutter) and differential amplification of alleles. An important feature of MEGASAT is the generation of histograms of the length-frequency distributions of amplification products for each locus and each individual. These histograms, analogous to electropherograms traditionally used to score microsatellite genotypes, enable rapid evaluation and editing of automatically scored genotypes. MEGASAT is written in Perl, runs on Windows, Mac OS X and Linux systems, and includes a simple graphical user interface. We demonstrate MEGASAT using data from guppy, Poecilia reticulata. We genotype 1024 guppies at 43 microsatellites per run on an Illumina MiSeq sequencer. We evaluated the accuracy of automatically called genotypes using two methods, based on pedigree and repeat genotyping data, and obtained estimates of mean genotyping error rates of 0.021 and 0.012. In both estimates, three loci accounted for a disproportionate fraction of genotyping errors; conversely, 26 loci were scored with 0-1 detected error (error rate ≤0.007). Our results show that with appropriate selection of loci, automated genotyping of microsatellite loci can be achieved with very high throughput, low genotyping error and very low genotyping costs.

BackgroundSexually dimorphic phenotypes are generally associated with differential gene expression between the sexes. The study of molecular evolution and genomic location of these differentially expressed, or sex-biased, genes is important for understanding inter-sexual divergence under sex-specific selection pressures. Teleost fish provide a unique opportunity to examine this divergence in the presence of variable sex-determination mechanisms of recent origin. The guppy, Poecilia reticulata, displays sexual dimorphism in size, ornaments, and behavior, traits shaped by natural and sexual selection in the wild.ResultsTo gain insight into molecular mechanisms underlying the guppy’s sexual dimorphism, we assembled a reference transcriptome combining genome-independent as well as genome-guided assemblies and analyzed sex-biased gene expression between different tissues of adult male and female guppies. We found tissue-associated sex-biased expression of genes related to pigmentation, signal transduction, and spermatogenesis in males; and growth, cell-division, extra-cellular matrix organization, nutrient transport, and folliculogenesis in females. While most sex-biased genes were randomly distributed across linkage groups, we observed accumulation of ovary-biased genes on the sex linkage group, LG12. Both testis-biased and ovary-biased genes showed a significantly higher rate of non-synonymous to synonymous substitutions (dN/dS) compared to unbiased genes. However, in somatic tissues only female-biased genes, including those co-expressed in multiple tissues, showed elevated ratios of non-synonymous substitutions.ConclusionsOur work identifies a set of annotated gene products that are candidate factors affecting sexual dimorphism in guppies. The differential genomic distribution of gonad-biased genes provides evidence for sex-specific selection pressures acting on the nascent sex chromosomes of the guppy. The elevated rates of evolution of testis-biased and female-biased genes indicate differing evolution under distinct selection pressures on the reproductive versus non-reproductive tissues.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-400) contains supplementary material, which is available to authorized users.

Convergent evolution represents one of the best lines of evidence for adaptation, but few cases of phenotypic convergence are understood at the genetic level. Guppies inhabiting the Northern Mountain Range of Trinidad provide a classic example of phenotypic convergent evolution, where adaptation to low or high predation environments has been found for a variety of traits. A major advantage of this system is the possibility of long-term experimental studies in nature, including transplantation from high to low predation sites. We used genome scans of guppies from three natural high and low predation populations and from two experimentally established populations and their sources to examine whether phenotypic convergent evolution leaves footprints at the genome level. We used population-genetic modelling approaches to reconstruct the demographic history and migration among sampled populations. Naturally colonized low predation populations had signatures of increased effective population size since colonization, while introduction populations had signatures of decreased effective population size. Only a small number of regions across the genome had signatures of selection in all natural populations. However, the two experimental populations shared many genomic regions under apparent selection, more than expected by chance. This overlap coupled with a population decrease since introduction provides evidence for convergent selection occurring in the two introduced populations. The lack of genetic convergence in the natural populations suggests that convergent evolution is lacking in these populations or that the effects of selection become difficult to detect after a long-time period.

BackgroundNext-generation sequencing is providing researchers with a relatively fast and affordable option for developing genomic resources for organisms that are not among the traditional genetic models. Here we present a de novo assembly of the guppy (Poecilia reticulata) transcriptome using 454 sequence reads, and we evaluate potential uses of this transcriptome, including detection of sex-specific transcripts and deployment as a reference for gene expression analysis in guppies and a related species. Guppies have been model organisms in ecology, evolutionary biology, and animal behaviour for over 100 years. An annotated transcriptome and other genomic tools will facilitate understanding the genetic and molecular bases of adaptation and variation in a vertebrate species with a uniquely well known natural history.ResultsWe generated approximately 336 Mbp of mRNA sequence data from male brain, male body, female brain, and female body. The resulting 1,162,670 reads assembled into 54,921 contigs, creating a reference transcriptome for the guppy with an average read depth of 28×. We annotated nearly 40% of this reference transcriptome by searching protein and gene ontology databases. Using this annotated transcriptome database, we identified candidate genes of interest to the guppy research community, putative single nucleotide polymorphisms (SNPs), and male-specific expressed genes. We also showed that our reference transcriptome can be used for RNA-sequencing-based analysis of differential gene expression. We identified transcripts that, in juveniles, are regulated differently in the presence and absence of an important predator, Rivulus hartii, including two genes implicated in stress response. For each sample in the RNA-seq study, >50% of high-quality reads mapped to unique sequences in the reference database with high confidence. In addition, we evaluated the use of the guppy reference transcriptome for gene expression analyses in a congeneric species, the sailfin molly (Poecilia latipinna). Over 40% of reads from the sailfin molly sample aligned to the guppy transcriptome.ConclusionsWe show that next-generation sequencing provided a reliable and broad reference transcriptome. This resource allowed us to identify candidate gene variants, SNPs in coding regions, and sex-specific gene expression, and permitted quantitative analysis of differential gene expression.

Understanding the selective forces influencing genetic diversity is a fundamental goal of evolutionary ecology. The genes of the major histocompatibility complex (MHC) play a key role in the adaptive immune response of vertebrates and thus provide an excellent opportunity to examine the agents of selection on a functionally important gene. Here we examine the genetic architecture of the MHC class IIB genes in 10 wild populations of guppies (Poecilia reticulata) in Northern Trinidad. We have previously shown that these populations are significantly less diverged at the class IIB locus than expected based on neutral (microsatellite) loci. We now survey infection by Gyrodactylus turnbulli and G. bullatarudis, common parasitic worms that infect guppies, as a potential agent of homogenizing selection. We used a genetic algorithm to partition both additive and non-additive genetic effects of the five most common MHC allele types as well as a rare allele category. Although we found no evidence for non-additive effects, across the populations we found that one allele type (the a-type) had a significant negative additive effect on parasite load. Thus, individuals who had more copies of the a-type allele were infected with fewer gyrodactylus than individuals with fewer copies of the allele. These results not only link parasite infection with MHC genotype, they provide a mechanism of homogenizing selection across these otherwise disparate populations.

A major goal of modern evolutionary biology is to understand the causes and consequences of phenotypic plasticity, the ability of a single genotype to produce multiple phenotypes in response to variable environments. While ecological and quantitative genetic studies have evaluated models of the evolution of adaptive plasticity, some long-standing questions about plasticity require more mechanistic approaches. Here, we address two of those questions: does plasticity facilitate adaptive evolution? And do physiological costs place limits on plasticity? We examine these questions by comparing genetically and plastically regulated behavioural variation in sailfin mollies (Poecilia latipinna), which exhibit striking variation in plasticity for male mating behaviour. In this species, some genotypes respond plastically to a change in the social environment by switching between primarily courting and primarily sneaking behaviour. In contrast, other genotypes have fixed mating strategies (either courting or sneaking) and do not display plasticity. We found that genetic and plastic variation in behaviour were accompanied by partially, but not completely overlapping changes in brain gene expression, in partial support of models that predict that plasticity can facilitate adaptive evolution. We also found that behavioural plasticity was accompanied by broader and more robust changes in brain gene expression, suggesting a substantial physiological cost to plasticity. We also observed that sneaking behaviour, but not courting, was associated with upregulation of genes involved in learning and memory, suggesting that sneaking is more cognitively demanding than courtship.