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.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data.

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.

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.

Independent Origin of XY and ZW Sex Determination Mechanisms in Mosquitofish Sister Species.

Fish are known for the outstanding variety of their sex determination mechanisms and sex chromosome systems. The western (Gambusia affinis) and eastern mosquitofish (G. holbrooki) are sister species for which different sex determination mechanisms have been described: ZZ/ZW for G. affinis and XX/XY for G. holbrooki Here, we carried out restriction-site associated DNA (RAD-) and pool sequencing (Pool-seq) to characterize the sex chromosomes of both species. We found that the ZW chromosomes of G. affinis females and the XY chromosomes of G. holbrooki males correspond to different linkage groups, and thus evolved independently from separate autosomes. In interspecific hybrids, the Y chromosome is dominant over the W chromosome, and X is dominant over Z. In G. holbrooki, we identified a candidate region for the Y-linked melanic pigmentation locus, a rare male phenotype that constitutes a potentially sexually antagonistic trait and is associated with other such characteristics, e.g., large body size and aggressive behavior. We developed a SNP-based marker in the Y-linked allele of GIPC PDZ domain containing family member 1 (gipc1), which was linked to melanism in all tested G. holbrooki populations. This locus represents an example for a color locus that is located in close proximity to a putative sex determiner, and most likely substantially contributed to the evolution of the Y.

Analysis of the putative tumor suppressor gene cdkn2ab in pigment cells and melanoma of Xiphophorus and medaka.

In humans, the CDKN2A locus encodes two transcripts, INK4A and ARF. Inactivation of either one by mutations or epigenetic changes is a frequent signature of malignant melanoma and one of the most relevant entry points for melanomagenesis. To analyze whether cdkn2ab, the fish ortholog of CDKN2A, has a similar function as its human counterpart, we studied its action in fish models for human melanoma. Overexpression of cdkn2ab in a Xiphophorus melanoma cell line led to decreased proliferation and induction of a senescence-like phenotype, indicating a melanoma-suppressive function analogous to mammals. Coexpression of Xiphophorus cdkn2ab in medaka transgenic for the mitfa:xmrk melanoma-inducing gene resulted in full suppression of melanoma development, whereas CRISPR/Cas9 knockout of cdkn2ab resulted in strongly enhanced tumor growth. In summary, this provides the first functional evidence that cdkn2ab acts as a potent tumor suppressor gene in fish melanoma models.

The Colorful Sex Chromosomes of Teleost Fish.

Teleost fish provide some of the most intriguing examples of sexually dimorphic coloration, which is often advantageous for only one of the sexes. Mapping studies demonstrated that the genetic loci underlying such color patterns are frequently in tight linkage to the sex-determining locus of a species, ensuring sex-specific expression of the corresponding trait. Several genes affecting color synthesis and pigment cell development have been previously described, but the color loci on the sex chromosomes have mostly remained elusive as yet. Here, we summarize the current knowledge about the genetics of such color loci in teleosts, mainly from studies on poeciliids and cichlids. Further studies on these color loci will certainly provide important insights into the evolution of sex chromosomes.

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.

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.