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

Extreme Y chromosome polymorphism corresponds to five male reproductive morphs of a freshwater fish

Loss of recombination between sex chromosomes often depletes Y chromosomes of functional content and genetic variation, which might limit their potential to generate adaptive diversity. Males of the freshwater fish Poecilia parae occur as one of five discrete morphs, all of which shoal together in natural populations where morph frequency has been stable for over 50 years. Each morph uses a different complex reproductive strategy and morphs differ dramatically in colour, body size and mating behaviour. Morph phenotype is passed perfectly from father to son, indicating there are five Y haplotypes segregating in the species, which encode the complex male morph characteristics. Here, we examine Y diversity in natural populations of P. parae. Using linked-read sequencing on multiple P. parae females and males of all five morphs, we find that the genetic architecture of the male morphs evolved on the Y chromosome after recombination suppression had occurred with the X. Comparing Y chromosomes between each of the morphs, we show that, although the Ys of the three minor morphs that differ in colour are highly similar, there are substantial amounts of unique genetic material and divergence between the Ys of the three major morphs that differ in reproductive strategy, body size and mating behaviour. Altogether, our results suggest that the Y chromosome is able to overcome the constraints of recombination loss to generate extreme diversity, resulting in five discrete Y chromosomes that control complex reproductive strategies.

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.

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.

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.

Comprehensive analysis of Sichuan white geese (Anser cygnoides) transcriptome

High-throughput RNA sequencing was performed for comprehensively analyzing the transcriptome of geese. A total of 28,803,759 bp of raw sequence data was generated by 454 GS Flx+. After removal of adaptor sequences, 28,730,361 bp remained and 117,279 reads were obtained, with an average length of 244 bases. Simultaneously, complementary DNA samples from two different reproductive stages of goose ovarian, hypothalamus and pituitary tissue were sequenced separately using Illumina MiSeq platform. A total of 12 688 673 148 bp of raw sequence data were generated by Illumina MiSeq. After removal of adaptor sequences, 8 198 126 562 bp remained and 60 382 786 clean reads were obtained, with an average length of 135 bases. Assembly of all the reads from both 454 Flx+ and Illumina platforms formed 56,839 contigs. The sequence size ranges from 38 to 28,206 bp in size, with an average size of 2584 bp and an N50 of 4624. The assembly produced a substantial number of large contigs: 35,545 (62.5%) were longer than 1 kb, of which 8850 (15.6%) were longer than 5 kb. The sequencing depth was 85 X on average. We performed comprehensive function annotations on unigenes including protein sequence similarity, gene ontology (GO) term classification, and Kyoto Encylcopedia of Genes and Genomes (KEGG) pathway enrichment. GO analysis showed that approximately 63% of the contigs had annotation information, among the 35,953 annotated isotigs in Nr database, 24,783 (68.9%) sequences were assigned with one or more GO terms. There were 14,634 (40.7%) isotigs for biological processes, 10,557(29.3%) isotigs for cellular component, 22,607 (62.9%) isotigs for molecular function. The result of KEGG pathway mapping 8926 sequences had the pathway annotation, and took part in 477 pathways. Additionally, 10,685 simple sequence repeat (SSR) markers were identified from the assembled sequences. The most frequent repeat motifs were trinucleotides, which accounted for 53.03% of all SSRs, followed by dinucleotides (39.9%), tetranucleotides (5.08%), pentanucleotides (1.68%) and hexanucleotides (0.32%). Transcriptome sequencing on mixture issue of the geese yielded substantial transcriptional sequences and potentially useful SSR markers which provide an important data source for geese research.

Mining antimicrobial peptides from small open reading frames in Ciona intestinalis

Though being able to encode various kinds of bioactive peptides, small open reading frames (sORFs) are poorly annotated in many genomic data. The present study was conducted to evaluate the potential of sORFs in encoding antimicrobial peptides (AMPs) in the basal chordate model Ciona intestinalis. About 4.8 m genomic sequence was first retrieved for sORFs mining by the program sORF finder, then the sORFs were translated into amino acid sequences for AMP prediction via CAMP server, and thereafter, ten putative AMPs were selected for expression and antimicrobial activity validation. In total, over 180 peptides deduced from the sORFs were predicted to be AMPs. Among the ten tested peptides, six were found to have significant expressed sequence tag matches, providing strong evidence for gene expression; five were proved to be active against the bacterial strains. These results indicate that many sORFs in C. intestinalis genome contain AMP information. This work can serve as an important initial step to investigate the role of sORFs in the innate defense of C. intestinalis.

Application of single nucleotide polymorphisms to non-model species: a technical review

Single nucleotide polymorphisms (SNPs) have gained wide use in humans and model species and are becoming the marker of choice for applications in other species. Technology that was developed for work in model species may provide useful tools for SNP discovery and genotyping in non-model organisms. However, SNP discovery can be expensive, labour intensive, and introduce ascertainment bias. In addition, the most efficient approaches to SNP discovery will depend on the research questions that the markers are to resolve as well as the focal species. We discuss advantages and disadvantages of several past and recent technologies for SNP discovery and genotyping and summarize a variety of SNP discovery and genotyping studies in ecology and evolution.

Genomic resources for a model in adaptation and speciation research: characterization of the Poecilia mexicana transcriptome

Background

Elucidating the genomic basis of adaptation and speciation is a major challenge in natural systems with large quantities of environmental and phenotypic data, mostly because of the scarcity of genomic resources for non-model organisms. The Atlantic molly (Poecilia mexicana, Poeciliidae) is a small livebearing fish that has been extensively studied for evolutionary ecology research, particularly because this species has repeatedly colonized extreme environments in the form of caves and toxic hydrogen sulfide containing springs. In such extreme environments, populations show strong patterns of adaptive trait divergence and the emergence of reproductive isolation. Here, we used RNA-sequencing to assemble and annotate the first transcriptome of P. mexicana to facilitate ecological genomics studies in the future and aid the identification of genes underlying adaptation and speciation in the system.

Description

We provide the first annotated reference transcriptome of P. mexicana. Our transcriptome shows high congruence with other published fish transcriptomes, including that of the guppy, medaka, zebrafish, and stickleback. Transcriptome annotation uncovered the presence of candidate genes relevant in the study of adaptation to extreme environments. We describe general and oxidative stress response genes as well as genes involved in pathways induced by hypoxia or involved in sulfide metabolism. To facilitate future comparative analyses, we also conducted quantitative comparisons between P. mexicana from different river drainages. 106,524 single nucleotide polymorphisms were detected in our dataset, including potential markers that are putatively fixed across drainages. Furthermore, specimens from different drainages exhibited some consistent differences in gene regulation.

Conclusions

Our study provides a valuable genomic resource to study the molecular underpinnings of adaptation to extreme environments in replicated sulfide spring and cave environments. In addition, this study adds to the increasing number of genomic resources in the family Poeciliidae, which are widely used in comparative analyses of behavior, ecology, evolution, and medical genetics.

Molecular analysis of the sex chromosomes of the platyfish Xiphophorus maculatus: Towards the identification of a new type of master sexual regulator in vertebrates

In contrast to mammals and birds, fish display an amazing diversity of genetic sex determination systems, with frequent changes during evolution possibly associated with the emergence of new sex chromosomes and sex-determining genes. To better understand the molecular and evolutionary mechanisms driving this diversity, several fish models are studied in parallel. Besides the medaka (Oryzias latipes Temminck and Schlegel, 1846) for which the master sex-determination gene has been identified, one of the most advanced models for studying sex determination is the Southern platyfish (Xiphophorus maculatus, Günther 1966). Xiphophorus maculatus belongs to the Poeciliids, a family of live-bearing freshwater fish, including platyfish, swordtails and guppies that perfectly illustrates the diversity of genetic sex-determination mechanisms observed in teleosts. For X. maculatus, bacterial artificial chromosome contigs covering the sex-determination region of the X and Y sex chromosomes have been constructed. Initial molecular analysis demonstrated that the sex-determination region is very unstable and frequently undergoes duplications, deletions, inversions and other rearrangements. Eleven gene candidates linked to the master sex-determining gene have been identified, some of them corresponding to pseudogenes. All putative genes are present on both the X and the Y chromosomes, suggesting a poor degree of differentiation and a young evolutionary age for platyfish sex chromosomes. When compared with other fish and tetrapod genomes, syntenies were detected only with autosomes. This observation supports an independent origin of sex chromosomes, not only in different vertebrate lineages but also between different fish species.

Sequencing and characterization of the guppy (Poecilia reticulata) transcriptome

Background

Next-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.

Results

We 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.

Conclusions

We 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.

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.

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.