Characterization of sex chromosomes in rainbow trout and coho salmon using fluorescence in situ hybridization (FISH)

The Gene Toolkit Implicated in Functional Sex in Sparidae Hermaphrodites: Inferences From Comparative Transcriptomics

Sex-biased gene expression is the mode through which sex dimorphism arises from a nearly identical genome, especially in organisms without genetic sex determination. Teleost fishes show great variations in the way the sex phenotype forms. Among them, Sparidae, that might be considered as a model family displays a remarkable diversity of reproductive modes. In this study, we sequenced and analyzed the sex-biased transcriptome in gonads and brain (the tissues with the most profound role in sexual development and reproduction) of two sparids with different reproductive modes: the gonochoristic common dentex, Dentex dentex, and the protandrous hermaphrodite gilthead seabream, Sparus aurata. Through comparative analysis with other protogynous and rudimentary protandrous sparid transcriptomes already available, we put forward common male and female-specific genes and pathways that are probably implicated in sex-maintenance in this fish family. Our results contribute to the understanding of the complex processes behind the establishment of the functional sex, especially in hermaphrodite species and set the groundwork for future experiments by providing a gene toolkit that can improve efforts to control phenotypic sex in finfish in the ever-increasingly important field of aquaculture.

Evolutionary and developmental dynamics of sex-biased gene expression in common frogs with proto-Y chromosomes

BACKGROUND

The patterns of gene expression on highly differentiated sex chromosomes differ drastically from those on autosomes, due to sex-specific patterns of selection and inheritance. As a result, X chromosomes are often enriched in female-biased genes (feminization) and Z chromosomes in male-biased genes (masculinization). However, it is not known how quickly sexualization of gene expression and transcriptional degeneration evolve after sex-chromosome formation. Furthermore, little is known about how sex-biased gene expression varies throughout development.

RESULTS

We sample a population of common frogs (Rana temporaria) with limited sex-chromosome differentiation (proto-sex chromosome), leaky genetic sex determination evidenced by the occurrence of XX males, and delayed gonadal development, meaning that XY individuals may first develop ovaries before switching to testes. Using high-throughput RNA sequencing, we investigate the dynamics of gene expression throughout development, spanning from early embryo to froglet stages. Our results show that sex-biased expression affects different genes at different developmental stages and increases during development, reaching highest levels in XX female froglets. Additionally, sex-biased gene expression depends on phenotypic, rather than genotypic sex, with similar expression in XX and XY males; correlates with gene evolutionary rates; and is not localized to the proto-sex chromosome nor near the candidate sex-determining gene Dmrt1.

CONCLUSIONS

The proto-sex chromosome of common frogs does not show evidence of sexualization of gene expression, nor evidence for a faster rate of evolution. This challenges the notion that sexually antagonistic genes play a central role in the initial stages of sex-chromosome evolution.

Tissue Specificity and Dynamics of Sex-Biased Gene Expression in a Common Frog Population with Differentiated, Yet Homomorphic, Sex Chromosomes

Sex-biased genes are central to the study of sexual selection, sexual antagonism, and sex chromosome evolution. We describe a comprehensive de novo assembled transcriptome in the common frog Rana temporaria based on five developmental stages and three adult tissues from both sexes, obtained from a population with karyotypically homomorphic but genetically differentiated sex chromosomes. This allows the study of sex-biased gene expression throughout development, and its effect on the rate of gene evolution while accounting for pleiotropic expression, which is known to negatively correlate with the evolutionary rate. Overall, sex-biased genes had little overlap among developmental stages and adult tissues. Late developmental stages and gonad tissues had the highest numbers of stage- or tissue-specific genes. We find that pleiotropic gene expression is a better predictor than sex bias for the evolutionary rate of genes, though it often interacts with sex bias. Although genetically differentiated, the sex chromosomes were not enriched in sex-biased genes, possibly due to a very recent arrest of XY recombination. These results extend our understanding of the developmental dynamics, tissue specificity, and genomic localization of sex-biased genes.

Utility of pooled sequencing for association mapping in nonmodel organisms

High-density genome-wide sequencing increases the likelihood of discovering genes of major effect and genomic structural variation in organisms. While there is an increasing availability of reference genomes across broad taxa, the greatest limitation to whole-genome sequencing of multiple individuals continues to be the costs associated with sequencing. To alleviate excessive costs, pooling multiple individuals with similar phenotypes and sequencing the homogenized DNA (Pool-Seq) can achieve high genome coverage, but at the loss of individual genotypes. Although Pool-Seq has been an effective method for association mapping in model organisms, it has not been frequently utilized in natural populations. To extend bioinformatic tools for rapid implementation of Pool-Seq data in nonmodel organisms, we developed a pipeline called PoolParty and illustrate its effectiveness in genetic association mapping. Alignment expectations based on five pooled Chinook salmon (Oncorhynchus tshawytscha) libraries showed that approximately 48% genome coverage per library could be achieved with reasonable sequencing effort. We additionally examined male and female O. tshawytscha libraries to illustrate how Pool-Seq techniques can successfully map known genes associated with functional differences among sexes such as growth hormone 2. Finally, we compared pools of individuals of different spawning ages for each sex to discover novel genes involved with age at maturity in O. tshawytscha such as opsin4 and transmembrane protein19. While not appropriate for every system, Pool-Seq data processed by the PoolParty pipeline is a practical method for identifying genes of major effect in nonmodel organisms when high genome coverage is necessary and cost is a limiting factor.

The transcriptomic signature of different sexes in two protogynous hermaphrodites: Insights into the molecular network underlying sex phenotype in fish

Sex differentiation is a puzzling problem in fish due to the variety of reproductive systems and the flexibility of their sex determination mechanisms. The Sparidae, a teleost family, reflects this remarkable diversity of sexual mechanisms found in fish. Our aim was to capture the transcriptomic signature of different sexes in two protogynous hermaphrodite sparids, the common pandora Pagellus erythrinus and the red porgy Pagrus pagrus in order to shed light on the molecular network contributing to either the female or the male phenotype in these organisms. Through RNA sequencing, we investigated sex-specific differences in gene expression in both species' brains and gonads. The analysis revealed common male and female specific genes/pathways between these protogynous fish. Whereas limited sex differences found in the brain indicate a sexually plastic tissue, in contrast, the great amount of sex-biased genes observed in gonads reflects the functional divergence of the transformed tissue to either its male or female character. Α common "crew" of well-known molecular players is acting to preserve either sex identity of the gonad in these fish. Lastly, this study lays the ground for a deeper understanding of the complex process of sex differentiation in two species with an evolutionary significant reproductive system.

Large-scale transcriptome sequencing reveals novel expression patterns for key sex-related genes in a sex-changing fish

Background

Teleost fishes exhibit remarkably diverse and plastic sexual developmental patterns. One of the most astonishing is the rapid socially controlled female-to-male (protogynous) sex change observed in bluehead wrasses (Thalassoma bifasciatum). Such functional sex change is widespread in marine fishes, including species of commercial importance, yet its underlying molecular basis remains poorly explored.

Methods

RNA sequencing was performed to characterize the transcriptomic profiles and identify genes exhibiting sex-biased expression in the brain (forebrain and midbrain) and gonads of bluehead wrasses. Functional annotation and enrichment analysis were carried out for the sex-biased genes in the gonad to detect global differences in gene products and genetic pathways between males and females.

Results

Here we report the first transcriptomic analysis for a protogynous fish. Expression comparison between males and females reveals a large set of genes with sex-biased expression in the gonad, but relatively few such sex-biased genes in the brain. Functional annotation and enrichment analysis suggested that ovaries are mainly enriched for metabolic processes and testes for signal transduction, particularly receptors of neurotransmitters and steroid hormones. When compared to other species, many genes previously implicated in male sex determination and differentiation pathways showed conservation in their gonadal expression patterns in bluehead wrasses. However, some critical female-pathway genes (e.g., rspo1 and wnt4b) exhibited unanticipated expression patterns. In the brain, gene expression patterns suggest that local neurosteroid production and signaling likely contribute to the sex differences observed.

Conclusions

Expression patterns of key sex-related genes suggest that sex-changing fish predominantly use an evolutionarily conserved genetic toolkit, but that subtle variability in the standard sex-determination regulatory network likely contributes to sexual plasticity in these fish. This study not only provides the first molecular data on a system ideally suited to explore the molecular basis of sexual plasticity and tissue re-engineering, but also sheds some light on the evolution of diverse sex determination and differentiation systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-015-0044-8) contains supplementary material, which is available to authorized users.

Chromosome landmarks and autosome-sex chromosome translocations in Rumex hastatulus, a plant with XX/XY1Y2 sex chromosome system

Rumex hastatulus is the North American endemic dioecious plant with heteromorphic sex chromosomes. It is differentiated into two chromosomal races: Texas (T) race characterised by a simple XX/XY sex chromosome system and North Carolina (NC) race with a polymorphic XX/XY1Y2 sex chromosome system. The gross karyotype morphology in NC race resembles the derived type, but chromosomal changes that occurred during its evolution are poorly understood. Our C-banding/DAPI and fluorescence in situ hybridization (FISH) experiments demonstrated that Y chromosomes of both races are enriched in DAPI-positive sequences and that the emergence of polymorphic sex chromosome system was accompanied by the break of ancestral Y chromosome and switch in the localization of 5S rDNA, from autosomes to sex chromosomes (X and Y2). Two contrasting domains were detected within North Carolina Y chromosomes: the older, highly heterochromatinised, inherited from the original Y chromosome and the younger, euchromatic, representing translocated autosomal material. The flow-cytometric DNA estimation showed ∼3.5 % genome downsizing in the North Carolina race. Our results are in contradiction to earlier reports on the lack of heterochromatin within Y chromosomes of this species and enable unambiguous identification of autosomes involved in the autosome-heterosome translocation, providing useful chromosome landmarks for further studies on the karyotype and sex chromosome differentiation in this species.

Y-specific sequences and polymorphisms in rainbow trout and Chinook salmon

Improved methods for genetically sexing salmonids and for characterization of Y-chromosome homologies between species can contribute to understanding the evolution of sex chromosomes and sex-determining mechanisms. In this study we have explored 12.5 kb of Y-chromosome-specific sequence flanking the previously described OtY2 locus in Chinook salmon (Oncorhynchus tshawytscha) and 21 kb of homologous rainbow trout (Oncorhynchus mykiss) Y-chromosome-specific sequence. This is the first confirmed Y-specific sequence for rainbow trout. New Y-specific markers are described for Chinook salmon (OtY3) and rainbow trout (OmyY1), which are readily detected by PCR assays and are advantageous because they also produce autosomal control amplification products. Additionally, AFLP analysis of Chinook salmon yielded another potential Y-chromosome marker. These descriptions will facilitate genotypic sexing and should be useful for population studies of Y-chromosome polymorphisms and for future studies to characterize what appears to be a common sex-determining mechanism between these species.

Cytogenetic analysis of three Italian populations of Coregonus lavaretus (Pisces, Salmoniformes) with chromosomal localization of major and minor ribosomal genes, and telomeric repeats

The European whitefish, Coregonus lavaretus, widely distributed in freshwater of northern Europe and introduced into the major lakes of northern Italy, has been restocked in central Italian lakes. In accordance with current managing practices, a reduced number of spawners contribute to reproduction within each lake and a certain degree of isolation is to be expected between populations from different lakes, resulting in the rapid fixing of chromosomal changes. A detailed survey of three populations from different lakes was carried out using classical and molecular cytogenetic techniques, to verify if specific chromosomal markers are present in the distinct populations. The comparative analysis revealed intraspecific variability of NORs and fixed differences in their number in the three populations. A co-localization of major and minor rRNA genes on one chromosome site was also observed. The original data regarding the chromosome mapping of the (TTAGGG)(n) telomeric repeat obtained in this study, demonstrated their exclusively terminal distribution, and a conspicuous inter-chromosomal variation in the number of repeats. The results are compared with data available for populations from native geographic ranges.

Identification of the sex-determining locus of Atlantic salmon (Salmo salar) on chromosome 2

We have integrated data from linkage mapping, physical mapping and karyotyping to gain a better understanding of the sex-determining locus, SEX, in Atlantic salmon (Salmo salar). SEX has been mapped to Atlantic salmon linkage group 1 (ASL1) and is associated with several microsatellite markers. We have used probes designed from the flanking regions of these sex-linked microsatellite markers to screen a bacterial artificial chromosome (BAC) library, representing an 11.7x coverage of the Atlantic salmon genome, which has been HindIII fingerprinted and assembled into contigs. BACs containing sex-linked microsatellites and their related contigs have been identified and representative BACs have been placed on the Atlantic salmon chromosomes by fluorescent in situ hybridization (FISH). This identified chromosome 2, a large metacentric, as the sex chromosome. By positioning several BACs on this chromosome by FISH, it was possible to orient ASL1 with respect to chromosome 2. The region containing SEX appears to lie on the long arm between marker Ssa202DU and a region of heterochromatin identified by DAPI staining. BAC end-sequencing of clones within sex-linked contigs revealed five hitherto unmapped genes along the sex chromosome. We are using an in silico approach coupled with physical probing of the BAC library to extend the BAC contigs to provide a physical map of ASL1, with a view to sequencing chromosome 2 and, in the process, identifying the sex-determining gene.

Polymorphism and differentiation of rainbow trout Y chromosomes

Most fish species show little morphological differentiation in the sex chromosomes. We have coupled molecular and cytogenetic analyses to characterize the male-determining region of the rainbow trout (Oncorhynchus mykiss) Y chromosome. Four genetically diverse male clonal lines of this species were used for genetic and physical mapping of regions in the vicinity of the sex locus. Five markers were genetically mapped to the Y chromosome in these male lines, indicating that the sex locus was located on the same linkage group in each of the lines. We also confirmed the presence of a Y chromosome morphological polymorphism among these lines, with the Y chromosomes from two of the lines having the more common heteromorphic Y chromosome and two of the lines having Y chromosomes morphologically similar to the X chromosome. The fluorescence in situ hybridization (FISH) pattern of two probes linked to sex suggested that the sex locus is physically located on the long arm of the Y chromosome. Fishes appear to be an excellent group of organisms for studying sex chromosome evolution and differentiation in vertebrates because they show considerable variability in the mechanisms and (or) patterns involved in sex determination.Key words: sex chromosomes, sex markers, cytogenetics, rainbow trout, fish.

Does differential selection on the 5S rDNA explain why the rainbow trout sex chromosome heteromorphism is not linked to the SEX locus?

Many but not all rainbow trout strains have morphologically distinguishable sex chromosomes. In these strains, the short arm of the X has multiple copies of 5S rDNA and a bright DAPI band near the centromere, both of which are missing from the Y chromosome, which has a very small short arm. We examined the presence of these markers using fluorescence in situ hybridization (FISH) in four different YY clonal lines derived from different strains and compared the results with sexed fish of the Donaldson strain with the normal X/Y heteromorphism. The Y chromosome in two of the YY clonal lines (Arlee and Swanson) is indistinguishable from the X chromosome and it is positive for 5S rDNA and the DAPI bright band. On the other hand, both 5S rDNA sequences and the DAPI band were not found on the Y chromosome in Hot Creek and Clearwater which have the normal Y. Thus the presence of these two cytogenetic markers may account for the size difference between the short arm of the X and Y chromosome found in most rainbow trout strains. In fishes the expression of one type of 5S rRNA is restricted to oocytes and previous work suggests that although XX males are fairly common, XY females are rare, implying a selective disadvantage for XY females. A hypothesis is presented to explain why this sex chromosome heteromorphism is not closely linked to the SEX locus, which is found on the long arm of the Y chromosome in rainbow trout.

A consolidated linkage map for rainbow trout (Oncorhynchus mykiss)

Androgenetic doubled haploid progeny produced from a cross between the Oregon State University and Arlee clonal rainbow trout (Oncorhynchus mykiss) lines, used for a previous published rainbow trout map, were used to update the map with the addition of more amplified fragment length polymorphic (AFLP) markers, microsatellites, type I and allozyme markers. We have added more than 900 markers, bringing the total number to 1359 genetic markers and the sex phenotype including 799 EcoRI AFLPs, 174 PstI AFLPs, 226 microsatellites, 72 VNTR, 38 SINE markers, 29 known genes, 12 minisatellites, five RAPDs, and four allozymes. Thirty major linkage groups were identified. Synteny of linkage groups in our map with the outcrossed microsatellite map has been established for all except one linkage group in this doubled haploid cross. Putative homeologous relationships among linkage groups, resulting from the autotetraploid nature of the salmonid genome, have been revealed based on the placement of duplicated microsatellites and type I loci.

Comparative genome analysis of the primary sex-determining locus in salmonid fishes

We compared the Y-chromosome linkage maps for four salmonid species (Arctic charr, Salvelinus alpinus; Atlantic salmon, Salmo salar; brown trout, Salmo trutta; and rainbow trout, Oncorhynchus mykiss) and a putative Y-linked marker from lake trout (Salvelinus namaycush). These species represent the three major genera within the subfamily Salmoninae of the Salmonidae. The data clearly demonstrate that different Y-chromosomes have evolved in each of the species. Arrangements of markers proximal to the sex-determining locus are preserved on homologous, but different, autosomal linkage groups across the four species studied in detail. This indicates that a small region of DNA has been involved in the rearrangement of the sex-determining region. Placement of the sex-determining region appears telomeric in brown trout, Atlantic salmon, and Arctic charr, whereas an intercalary location for SEX may exist in rainbow trout. Three hypotheses are proposed to account for the relocation: translocation of a small chromosome arm; transposition of the sex-determining gene; or differential activation of a primary sex-determining gene region among the species.