Improved Reference Genome Uncovers Novel Sex-Linked Regions in the Guppy (Poecilia reticulata).

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 favor 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 color 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 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.

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.

Adenylate cyclase 5 is required for melanophore and male pattern development in the guppy (Poecilia reticulata).

Guppies (Poecilia reticulata) are colorful fish that have attracted the attention of pigmentation researchers for almost a century. Here, we report that the blond phenotype of the guppy is caused by a spontaneous mutation in the guppy ortholog of adenylate cyclase 5 (adcy5). Using double digest restriction site-associated DNA sequencing (ddRADseq) and quantitative trait locus (QTL) mapping, we linked the blond phenotype to a candidate region of 118 kb, in which we subsequently identified a 2-bp deletion in adcy5 that alters splicing and leads to a premature stop codon. We show that adcy5, which affects life span and melanoma growth in mouse, is required for melanophore development and formation of male orange pigmentation traits in the guppy. We find that some components of the male orange pattern are particularly sensitive to loss of Adcy5 function. Our work thus reveals a function for Adcy5 in patterning of fish color ornaments.

Population genomics of natural and experimental populations of guppies (Poecilia reticulata).

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.

Transcriptome assemblies for studying sex-biased gene expression in the guppy, Poecilia reticulata.

BACKGROUND: Sexually 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. RESULTS: To 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. CONCLUSIONS: Our 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.

Sex chromosome polymorphism in guppies.

Sex chromosomes differ from autosomes by dissimilar gene content and, at a more advanced stage of their evolution, also in structure and size. This is driven by the divergence of the Y or W from their counterparts, X and Z, due to reduced recombination and the resulting degeneration as well as the accumulation of sex-specific and sexually antagonistic genes. A paradigmatic example for Y-chromosome evolution is found in guppies. In these fishes, conflicting data exist for a morphological and molecular differentiation of sex chromosomes. Using molecular probes and the previously established linkage map, we performed a cytogenetic analysis of sex chromosomes. We show that the Y chromosome has a very large pseudoautosomal region, which is followed by a heterochromatin block (HCY) separating the subtelomeric male-specific region from the rest of the chromosome. Interestingly, the size of the HCY is highly variable between individuals from different population. The largest HCY was found in one population of Poecilia wingei, making the Y almost double the size of the X and the largest chromosome of the complement. Comparative analysis revealed that the Y chromosomes of different guppy species are homologous and share the same structure and organization. The observed size differences are explained by an expansion of the HCY, which is due to increased amounts of repetitive DNA. In one population, we observed also a polymorphism of the X chromosome. We suggest that sex chromosome-linked color patterns and other sexually selected genes are important for maintaining the observed structural polymorphism of sex chromosomes.

Multiple pigment cell types contribute to the black, blue, and orange ornaments of male guppies (Poecilia reticulata).

The fitness of male guppies (Poecilia reticulata) highly depends on the size and number of their black, blue, and orange ornaments. Recently, progress has been made regarding the genetic mechanisms underlying male guppy pigment pattern formation, but we still know little about the pigment cell organization within these ornaments. Here, we investigate the pigment cell distribution within the black, blue, and orange trunk spots and selected fin color patterns of guppy males from three genetically divergent strains using transmission electron microscopy. We identified three types of pigment cells and found that at least two of these contribute to each color trait. Further, two pigment cell layers, one in the dermis and the other in the hypodermis, contribute to each trunk spot. The pigment cell organization within the black and orange trunk spots was similar between strains. The presence of iridophores in each of the investigated color traits is consistent with a key role for this pigment cell type in guppy color pattern formation.

Pigment pattern formation in the guppy, Poecilia reticulata, involves the Kita and Csf1ra receptor tyrosine kinases.

Males of the guppy (Poecilia reticulata) vary tremendously in their ornamental patterns, which are thought to have evolved in response to a complex interplay between natural and sexual selection. Although the selection pressures acting on the color patterns of the guppy have been extensively studied, little is known about the genes that control their ontogeny. Over 50 years ago, two autosomal color loci, blue and golden, were described, both of which play a decisive role in the formation of the guppy color pattern. Orange pigmentation is absent in the skin of guppies with a lesion in blue, suggesting a defect in xanthophore development. In golden mutants, the development of the melanophore pattern during embryogenesis and after birth is affected. Here, we show that blue and golden correspond to guppy orthologs of colony-stimulating factor 1 receptor a (csf1ra; previously called fms) and kita. Most excitingly, we found that both genes are required for the development of the black ornaments of guppy males, which in the case of csf1ra might be mediated by xanthophore-melanophore interactions. Furthermore, we provide evidence that two temporally and genetically distinct melanophore populations contribute to the adult camouflage pattern expressed in both sexes: one early appearing and kita-dependent and the other late-developing and kita-independent. The identification of csf1ra and kita mutants provides the first molecular insights into pigment pattern formation in this important model species for ecological and evolutionary genetics.

Estimates of genetic differentiation measured by F(ST) do not necessarily require large sample sizes when using many SNP markers.

Population genetic studies provide insights into the evolutionary processes that influence the distribution of sequence variants within and among wild populations. F(ST) is among the most widely used measures for genetic differentiation and plays a central role in ecological and evolutionary genetic studies. It is commonly thought that large sample sizes are required in order to precisely infer F(ST) and that small sample sizes lead to overestimation of genetic differentiation. Until recently, studies in ecological model organisms incorporated a limited number of genetic markers, but since the emergence of next generation sequencing, the panel size of genetic markers available even in non-reference organisms has rapidly increased. In this study we examine whether a large number of genetic markers can substitute for small sample sizes when estimating F(ST). We tested the behavior of three different estimators that infer F(ST) and that are commonly used in population genetic studies. By simulating populations, we assessed the effects of sample size and the number of markers on the various estimates of genetic differentiation. Furthermore, we tested the effect of ascertainment bias on these estimates. We show that the population sample size can be significantly reduced (as small as n = 4-6) when using an appropriate estimator and a large number of bi-allelic genetic markers (k>1,000). Therefore, conservation genetic studies can now obtain almost the same statistical power as studies performed on model organisms using markers developed with next-generation sequencing.

Reasons for the invasive success of a guppy (Poecilia reticulata) population in Trinidad.

The introduction of non-native species into new habitats poses a major threat to native populations. Of particular interest, though often overlooked, are introductions of populations that are not fully reproductively isolated from native individuals and can hybridize with them. To address this important topic we used different approaches in a multi-pronged study, combining the effects of mate choice, shoaling behaviour and genetics. Here we present evidence that behavioural traits such as shoaling and mate choice can promote population mixing if individuals do not distinguish between native and foreign conspecifics. We examined this in the context of two guppy (Poecilia reticulata) populations that have been subject to an introduction and subsequent population mixing event in Trinidad. The introduction of Guanapo River guppies into the Turure River more than 50 years ago led to a marked reduction of the original genotype. In our experiments, female guppies did not distinguish between shoaling partners when given the choice between native and foreign individuals. Introduced fish are therefore likely to benefit from the protection of a shoal and will improve their survival chances as a result. The additional finding that male guppies do not discriminate between females on the basis of origin will further increase the process of population mixing, especially if males encounter mixed shoals. In a mesocosm experiment, in which the native and foreign populations were allowed to mate freely, we found, as expected on the basis of these behavioural interactions, that the distribution of offspring genotypes could be predicted from the proportions of the two types of founding fish. This result suggests that stochastic and environmental processes have reinforced the biological ones to bring about the genetic dominance of the invading population in the Turure River. Re-sampling the Turure for genetic analysis using SNP markers confirmed the population mixing process and showed that it is an on-going process in this river and has led to the nearly complete disappearance of the original genotype.

ShapePheno: unsupervised extraction of shape phenotypes from biological image collections.

MOTIVATION: Accurate large-scale phenotyping has recently gained considerable importance in biology. For example, in genome-wide association studies technological advances have rendered genotyping cheap, leaving phenotype acquisition as the major bottleneck. Automatic image analysis is one major strategy to phenotype individuals in large numbers. Current approaches for visual phenotyping focus predominantly on summarizing statistics and geometric measures, such as height and width of an individual, or color histograms and patterns. However, more subtle, but biologically informative phenotypes, such as the local deformation of the shape of an individual with respect to the population mean cannot be automatically extracted and quantified by current techniques. RESULTS: We propose a probabilistic machine learning model that allows for the extraction of deformation phenotypes from biological images, making them available as quantitative traits for downstream analysis. Our approach jointly models a collection of images using a learned common template that is mapped onto each image through a deformable smooth transformation. In a case study, we analyze the shape deformations of 388 guppy fish (Poecilia reticulata). We find that the flexible shape phenotypes our model extracts are complementary to basic geometric measures. Moreover, these quantitative traits assort the observations into distinct groups and can be mapped to polymorphic genetic loci of the sample set. AVAILABILITY: Code is available under: http://bioweb.me/GEBI CONTACT: theofanis.karaletsos@tuebingen.mpg.de; oliver.stegle@tuebingen.mpg.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Paired-end RAD-seq for de novo assembly and marker design without available reference.

MOTIVATION: Next-generation sequencing technologies have facilitated the study of organisms on a genome-wide scale. A recent method called restriction site associated DNA sequencing (RAD-seq) allows to sample sequence information at reduced complexity across a target genome using the Illumina platform. Single-end RAD-seq has proven to provide a large number of informative genetic markers in reference as well as non-reference organisms. RESULTS: Here, we present a method for de novo assembly of paired-end RAD-seq data in order to produce extended contigs flanking a restriction site. We were able to reconstruct one-tenth of the guppy genome represented by 200-500 bp contigs associated to EcoRI recognition sites. In addition, these contigs were used as reference allowing the detection of thousands of new polymorphic markers that are informative for mapping and population genetic studies in the guppy. AVAILABILITY: A perl and C++ implementation of the method demonstrated in this article is available under http://guppy.weigelworld.org/weigeldatabases/radMarkers/ as package RApiD. CONTACT: christine.dreyer@tuebingen.mpg.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A major QTL controls susceptibility to spinal curvature in the curveback guppy.

BACKGROUND: Understanding the genetic basis of heritable spinal curvature would benefit medicine and aquaculture. Heritable spinal curvature among otherwise healthy children (i.e. Idiopathic Scoliosis and Scheuermann kyphosis) accounts for more than 80% of all spinal curvatures and imposes a substantial healthcare cost through bracing, hospitalizations, surgery, and chronic back pain. In aquaculture, the prevalence of heritable spinal curvature can reach as high as 80% of a stock, and thus imposes a substantial cost through production losses. The genetic basis of heritable spinal curvature is unknown and so the objective of this work is to identify quantitative trait loci (QTL) affecting heritable spinal curvature in the curveback guppy. Prior work with curveback has demonstrated phenotypic parallels to human idiopathic-type scoliosis, suggesting shared biological pathways for the deformity. RESULTS: A major effect QTL that acts in a recessive manner and accounts for curve susceptibility was detected in an initial mapping cross on LG 14. In a second cross, we confirmed this susceptibility locus and fine mapped it to a 5 cM region that explains 82.6% of the total phenotypic variance. CONCLUSIONS: We identify a major QTL that controls susceptibility to curvature. This locus contains over 100 genes, including MTNR1B, a candidate gene for human idiopathic scoliosis. The identification of genes associated with heritable spinal curvature in the curveback guppy has the potential to elucidate the biological basis of spinal curvature among humans and economically important teleosts.

Gene duplication and divergence of long wavelength-sensitive opsin genes in the guppy, Poecilia reticulata.

Female preference for male orange coloration in the genus Poecilia suggests a role for duplicated long wavelength-sensitive (LWS) opsin genes in facilitating behaviors related to mate choice in these species. Previous work has shown that LWS gene duplication in this genus has resulted in expansion of long wavelength visual capacity as determined by microspectrophotometry (MSP). However, the relationship between LWS genomic repertoires and expression of LWS retinal cone classes within a given species is unclear. Our previous study in the related species, Xiphophorus helleri, was the first characterization of the complete LWS opsin genomic repertoire in conjunction with MSP expression data in the family Poeciliidae, and revealed the presence of four LWS loci and two distinct LWS cone classes. In this study we characterized the genomic organization of LWS opsin genes by BAC clone sequencing, and described the full range of cone cell types in the retina of the colorful Cumaná guppy, Poecilia reticulata. In contrast to X. helleri, MSP data from the Cumaná guppy revealed three LWS cone classes. Comparisons of LWS genomic organization described here for Cumaná to that of X. helleri indicate that gene divergence and not duplication was responsible for the evolution of a novel LWS haplotype in the Cumaná guppy. This lineage-specific divergence is likely responsible for a third additional retinal cone class not present in X. helleri, and may have facilitated the strong sexual selection driven by female preference for orange color patterns associated with the genus Poecilia.

Genome-wide single nucleotide polymorphisms reveal population history and adaptive divergence in wild guppies.

Adaptation of guppies (Poecilia reticulata) to contrasting upland and lowland habitats has been extensively studied with respect to behaviour, morphology and life history traits. Yet population history has not been studied at the whole-genome level. Although single nucleotide polymorphisms (SNPs) are the most abundant form of variation in many genomes and consequently very informative for a genome-wide picture of standing natural variation in populations, genome-wide SNP data are rarely available for wild vertebrates. Here we use genetically mapped SNP markers to comprehensively survey genetic variation within and among naturally occurring guppy populations from a wide geographic range in Trinidad and Venezuela. Results from three different clustering methods, Neighbor-net, principal component analysis (PCA) and Bayesian analysis show that the population substructure agrees with geographic separation and largely with previously hypothesized patterns of historical colonization. Within major drainages (Caroni, Oropouche and Northern), populations are genetically similar, but those in different geographic regions are highly divergent from one another, with some indications of ancient shared polymorphisms. Clear genomic signatures of a previous introduction experiment were seen, and we detected additional potential admixture events. Headwater populations were significantly less heterozygous than downstream populations. Pairwise F(ST) values revealed marked differences in allele frequencies among populations from different regions, and also among populations within the same region. F(ST) outlier methods indicated some regions of the genome as being under directional selection. Overall, this study demonstrates the power of a genome-wide SNP data set to inform for studies on natural variation, adaptation and evolution of wild populations.

Genetic linkage map of the guppy, Poecilia reticulata, and quantitative trait loci analysis of male size and colour variation.

We report construction of a genetic linkage map of the guppy genome using 790 single nucleotide polymorphism markers, integrated from six mapping crosses. The markers define 23 linkage groups (LGs), corresponding to the known haploid number of guppy chromosomes. The map, which spans a genetic length of 899 cM, includes 276 markers linked to expressed genes (expressed sequence tag), which have been used to derive broad syntenic relationships of guppy LGs with medaka chromosomes. This combined linkage map should facilitate the advancement of genetic studies for a wide variety of complex adaptive phenotypes relevant to natural and sexual selection in this species. We have used the linkage data to predict quantitative trait loci for a set of variable male traits including size and colour pattern. Contributing loci map to the sex LG for many of these traits.

Linkage analysis reveals the independent origin of Poeciliid sex chromosomes and a case of atypical sex inheritance in the guppy (Poecilia reticulata).

Among different teleost fish species, diverse sex-determining mechanisms exist, including environmental and genetic sex determination, yet chromosomal sex determination with male heterogamety (XY) prevails. Different pairs of autosomes have evolved as sex chromosomes among species in the same genus without evidence for a master sex-determining locus being identical. Models for evolution of Y chromosomes predict that male-advantageous genes become linked to a sex-determining locus and suppressed recombination ensures their co-inheritance. In the guppy, Poecilia reticulata, a set of genes responsible for adult male ornaments are linked to the sex-determining locus on the incipient Y chromosome. We have identified >60 sex-linked molecular markers to generate a detailed map for the sex linkage group of the guppy and compared it with the syntenic autosome 12 of medaka. We mapped the sex-determining locus to the distal end of the sex chromosome. We report a sex-biased distribution of recombination events in female and male meiosis on sex chromosomes. In one mapping cross, we observed sex ratio and male phenotype deviations and propose an atypical mode of genetic sex inheritance as its basis.

Natural variation of male ornamental traits of the guppy, Poecilia reticulata.

Male ornamental traits of the guppy, Poecilia reticulata, provide an outstanding example of natural variation in sex-linked male-advantageous traits that are shaped by both sexual and environmental selection. A substantial fraction of the underlying genes is known to be genetically linked to the sex-determining region on the differentiating Y-chromosome. Intercrosses between parental populations originating from geographically distant locations in East Trinidad and Cumaná (Venezuela) were used to study segregation of ornamental traits in male progeny. In addition, we performed backcrosses to compare segregation of ornaments in presence or absence of prominent traits linked to the Y-chromosome. Another backcross strategy involving XY females from the laboratory strain zebrinus maculatus allowed studying additive and dominant effects of alleles on two different Y-chromosomes on pattern formation. For genetic mapping, we have previously developed nuclear SNP markers linked to expressed genes, including several genes known to be important for pattern formation in other species. Of these candidate genes 15 were placed on 11 different linkage groups. Our phenotypic and genotypic analysis of progeny from mapping crosses and backcrosses suggests several genetic mechanisms that enhance natural variation, namely, additive effects of codominant alleles, suppressive actions of dominant alleles, and a complex interplay between sex-linked and autosomal cofactors.

ESTs and EST-linked polymorphisms for genetic mapping and phylogenetic reconstruction in the guppy, Poecilia reticulata.

BACKGROUND: The guppy, Poecilia reticulata, is a well-known model organism for studying inheritance and variation of male ornamental traits as well as adaptation to different river habitats. However, genomic resources for studying this important model were not previously widely available. RESULTS: With the aim of generating molecular markers for genetic mapping of the guppy, cDNA libraries were constructed from embryos and different adult organs to generate expressed sequence tags (ESTs). About 18,000 ESTs were annotated according to BLASTN and BLASTX results and the sequence information from the 3' UTRs was exploited to generate PCR primers for re-sequencing of genomic DNA from different wild type strains. By comparison of EST-linked genomic sequences from at least four different ecotypes, about 1,700 polymorphisms were identified, representing about 400 distinct genes. Two interconnected MySQL databases were built to organize the ESTs and markers, respectively. A robust phylogeny of the guppy was reconstructed, based on 10 different nuclear genes. CONCLUSION: Our EST and marker databases provide useful tools for genetic mapping and phylogenetic studies of the guppy.

Opsin gene duplication and diversification in the guppy, a model for sexual selection.

Identification of genes that control variation in adaptive characters is a prerequisite for understanding the processes that drive sexual and natural selection. Male coloration and female colour perception play important roles in mate choice in the guppy (Poecilia reticulata), a model organism for studies of natural and sexual selection. We examined a potential source for the known variation in colour perception, by analysing genomic and complementary DNA sequences of genes that code for visual pigment proteins. We find high sequence variability, both within and between populations, and expanded copy number for long-wave sensitive (LWS) opsin genes. Alleles with non-synonymous changes that suggest dissimilar spectral tuning properties occur in the same population and even in the same individual, and the high frequency of non-synonymous substitutions argues for diversifying selection acting on these proteins. Therefore, variability in tuning amino acids is partitioned within individuals and populations of the guppy, in contrast to variability for LWS at higher taxonomic levels in cichlids, a second model system for differentiation owing to sexual selection. Since opsin variability parallels the extreme male colour polymorphism within guppy populations, we suggest that mate choice has been a major factor driving the coevolution of opsins and male ornaments in this species.