- Research article
- Open Access
Structure and decay of a proto-Y region in Tilapia, Oreochromis niloticus
© Gammerdinger et al.; licensee BioMed Central Ltd. 2014
- Received: 16 July 2014
- Accepted: 29 October 2014
- Published: 17 November 2014
Sex-determination genes drive the evolution of adjacent chromosomal regions. Sexually antagonistic selection favors the accumulation of inversions that reduce recombination in regions adjacent to the sex-determination gene. Once established, the clonal inheritance of sex-linked inversions leads to the accumulation of deleterious alleles, repetitive elements and a gradual decay of sex-linked genes. This in turn creates selective pressures for the evolution of mechanisms that compensate for the unequal dosage of gene expression. Here we use whole genome sequencing to characterize the structure of a young sex chromosome and quantify sex-specific gene expression in the developing gonad.
We found an 8.8 Mb block of strong differentiation between males and females that corresponds to the location of a previously mapped sex-determiner on linkage group 1 of Oreochromis niloticus. Putatively disruptive mutations are found in many of the genes within this region. We also found a significant female-bias in the expression of genes within the block of differentiation compared to those outside the block of differentiation. Eight candidate sex-determination genes were identified within this region.
This study demonstrates a block of differentiation on linkage group 1, suggestive of an 8.8 Mb inversion encompassing the sex-determining locus. The enrichment of female-biased gene expression inside the proposed inversion suggests incomplete dosage compensation. This study helps establish a model for studying the early-to-intermediate stages of sex chromosome evolution.
- Linkage Group
- Intermediate Frequency
- Dosage Compensation
- Therian Mammal
- Antagonistic Selection
The classic model of sex chromosome evolution begins with the emergence of a new sex- determining gene on an autosome . The new sex-determiner may be linked with genes experiencing sexually antagonistic selection. Selection favors mechanisms, such as chromosomal inversions, that reduce recombination between the sex-determination locus and sexually antagonistic genes [2, 3]. The human sex chromosomes have undergone at least four such inversions, which may have limited recombination between the sex-determination locus and nearby sexually antagonistic genes .
Inversions create a clonally inherited chromosomal segment with a relatively small effective population size when compared to the rest of the genome . Sex chromosomes therefore become a haven for deleterious mutations and repetitive elements that are difficult to purge. These deleterious mutations accumulate via Muller’s Ratchet, as well as by hitchhiking with advantageous mutations [6, 7]. Degradation of functional genes on the Y- or W-chromosome leaves the homogametic sex carrying two functional copies of a particular gene, while the heterogametic sex carries only one functional copy. Therefore, mechanisms are needed to maintain appropriate expression of dosage-sensitive genes on emerging sex chromosomes [8–10]. In mammals, global dosage compensation is accomplished through X-inactivation [10, 11]. However, in many species, dosage compensation is partial and the expression of many genes is not compensated .
Some sex-determining genes are conserved for long periods of time. An example is Sry, a gene that has controlled sex-determination in therian mammals for approximately 180 million years [7, 12–14]. Other sex-determination genes hold sway for much shorter periods of time. There have been at least five transitions in the mechanism controlling sex-determination in rice fish (genus Oryzias) during the last 20 million years . Similarly rapid rates of sex chromosome evolution have been identified among sticklebacks (Family Gasterosteidae) .
The evolution of new sex-determination genes may have contributed to the rapid radiation of African cichlid fishes . Among the closely related haplochromine cichlids of Lake Malawi, sex-determination regions have been localized to linkage groups 3 (ZW), 5 (ZW), 7 (XY) and 20 (XY) [18–20]. Among tilapia cichlids, sex-determination regions have been localized to linkage groups 1 (XY), 3 (ZW) and 23 (XY) [21, 22]. Multiple sex-determination genes often segregate within a single species [18, 23]. The blue tilapia, Oreochromis aureus, segregates both an XY system on linkage group 1 and a ZW system on linkage group 3 . Some strains of the Nile tilapia, O. niloticus, have an XY system on linkage group 1, while others segregate an XY system on linkage group 23 [22, 24, 25].
The goal of this study was to characterize the sex-determination locus on linkage group 1 in O. niloticus. We took a family-based strategy, separately pooling males and females from two crosses, and performing whole genome sequencing on the pooled DNAs. We cataloged the density and frequency of single nucleotide polymorphisms (SNPs) and assessed their functional impact. We identified an 8.8 Mb block of differentiation suggestive of a Y-linked inversion on linkage group 1. We found high densities of functionally significant SNPs within this differentiated block. Analysis of gonadal transcriptomes demonstrated an enrichment of female-biased gene models within the inversion, which suggests that dosage compensation is incomplete in this strain of O. niloticus.
Sequencing of male and female DNA pools
We obtained ~202 million reads from the pool of male DNA and ~219 million reads from the pool of female DNA. 90.12% of the male and 90.67% of the female reads were aligned to the O. niloticus reference genome. Genome-wide coverage was slightly lower in males (32.97, standard deviation = 24.41 alignments per site), compared to females (36.68, standard deviation = 31.39 alignments per site).
Large block of divergence on linkage group 1
We also counted the number of positions per 10 kb window that were fixed in female pools and had a SNP in intermediate frequency in male pools, as would be consistent with females having two X chromosomes and males having an X and a Y chromosome, using Sex_SNP_finder_now.pl. There were 40,514 of these SNPs found across the genome. 18,277 (2,076.932/Mb) lay inside the differentiated block and 22,237 (24.197/Mb) lay outside. Among the 300 non-overlapping 10 kb windows with the highest frequency of these SNPs, 290 were found within the differentiated block on linkage group 1. The mean number of such SNPs per window was 21.81 (standard deviation = 13.84) within the differentiated block and only 0.33 (standard deviation = 1.29) outside of this region (Figures 1c, and 2c). The elevated FST, along with the abundance of intermediate frequency SNPs in males that are fixed in females, suggests that this region has limited, if any, recombination between the X and Y alleles.
We considered the possibility that this block of differentiation is an artifact of the process by which we selected individuals for sequencing. We initially screened individuals by genotyping two sex-linked microsatellites in order to confirm family identity and sex. We required males to demonstrate a heterogametic pattern and females to demonstrate a homogametic pattern for both markers. Five male and five female individuals were excluded by these criteria and may represent naturally sex-reversed individuals. The sharply defined edges of the block lie 4.22 Mb upstream and 3.37 Mb downstream of the microsatellites we genotyped (Figures 1 and 2), which would normally represent approximately 5 cM of genetic distance in this species . However, there is no evidence of an exponential decay of FST in the flanking regions as would be expected if there was recombination between the markers and the edges of the block. We also considered the possibility that the high level of differentiation might be due to an 8.8 Mb duplication on the Y. However, the depth of read coverage is relatively consistent across this entire linkage group. Additionally, cytogenetic studies have not revealed any evidence of heteromorphy in this chromosome pair as would arise from a translocation . The sum of the evidence suggests that this block of differentiation most likely reflects an 8.8 Mb inversion on the Y-chromosome.
The relatively small size of the putative inversion, and its location in the middle of the chromosome, make it challenging to characterize using standard cytogenetic techniques. Ideally, we would characterize the breakpoints, but we were unable to identify anomalous Illumina mate pairs near the ends of the inversion in our short insert libraries. Longer reads or more widely spaced mate pairs will be needed to characterize the breakpoints of the proposed inversion.
Functionally significant SNPs
Putative functional mutations in the proposed inversion
SNP location on LG1
Effect on the Y
Effect on the X
Frequency of SNP in males
Frequency of SNP in females
Ras-related protein R-Ras2 (LOC100693950)
AMP deaminase 3 (LOC100694225)
Zinc finger protein 821 (LOC100712266)
Zinc finger protein 821 (LOC100712266)
SAFB-like transcription modulator (LOC100711186)
Hepatic lipase (Lipc)
Ammonium transporter Rh type C 2 (LOC100706367)
AFG3-like protein 1 (LOC100702885)
CUB and sushi domain-containing protein 1 (LOC100698036)
CUB and sushi domain-containing protein 1 (LOC100698036)
Neuromedin-K receptor (LOC100693904)
Splice Site Acceptor Lost
Protein FAM176A (LOC100700039)
GC-rich sequence DNA-binding factor (LOC100700589)
Splice Site Acceptor Gain
BTB/POZ domain-containing protein KCTD3 (LOC100703295)
Hypothetical protein (LOC100705710)
Hypothetical protein (LOC100705710)
Nuclear factor of activated T-cells
The elevated density of high impact SNPs within the proposed inversion leads us to believe that deleterious alleles have begun to accumulate on this proto-Y. This is in accordance with the canonical model of heterogametic sex-chromosome evolution [2, 32] and empirical observations of the therian mammal Y-chromosome, Silene, Drosophila and tongue sole [4, 33–35].
Localization of the sex-determining gene
The multifactorial nature of sex-determination in this species causes difficulties for genetic mapping studies. An XX individual may develop as a male due to other genetic factors, or environmental effects on differentiation. These individuals would appear to be recombinant in the sex interval. We previously claimed to exclude Wilm’s tumor protein homolog (Wt1b) as the sex-determining gene on the basis of two recombinant individuals , but this conclusion is now in doubt. Conversely, the absence of recombination within the proposed inversion may preclude any further fine-mapping of the gene responsible for sex-determination.
Differences in gene expression
The block of differentiation on linkage group 1 comprises just more than 1% of the assembled genome and contains 257 RefSeq annotated genes. Cufflinks predicted 234 gene models within the block of differentiation and predicted 22,411 gene models across the entire transcriptome. Of the gene models that showed an FPKM of >0.05 in at least one sex, 7,977 gene models (37.4%) showed higher expression in males, while 13,375 (59.7%) gene models showed higher expression in females. Furthermore, within the inverted region, only 68 of these gene models (29.6%) showed a male bias (Additional file 3), while 162 of these gene models (69.2%) showed a female bias (Additional file 4). The enrichment of female biased gene models within the proposed inversion is significant (χ2 = 5.58, p <0.05). These data suggest that this sex chromosome is at an early-to-intermediate evolutionary phase where the degradation of a proto-Y has begun and expression of Y-linked genes in males is reduced. However, mechanisms for complete dosage compensation have yet to take hold.
Candidate sex determiners
Candidate genes in the proposed inversion
Gene location on LG1 (Mb)
SNP location on LG1
Pool frequency in males
Pool frequency in females
Transcription factor SOX-6 (LOC100694759)
Ras-related protein R-Ras2 (LOC100693950)
Suppression of tumorigenicity 5 protein (LOC100693420)
Ras association domain-containing protein 10 (LOC100693148)
AFG3-like protein 1 (LOC100702885)
Wilms tumor protein homolog (LOC100701078)
Estrogen-related receptor gamma (LOC100704106)
Growth regulation by estrogen in breast cancer 1 (GREB1)
First, we analyzed all SNPs that SnpEff classified as high impact mutations (Table 1). One prominent candidate within the proposed inversion is Ras-related protein R-Ras2 (10.49 Mb-10.51 Mb), which is part of the Ras-MEK-ERK pathway within the TGF-β signaling network . Alterations to the TGF-β network have been suggested as the mechanism for sex-determination in several fish species . Ras2 has been implicated as particularly important in the proliferation of cells  and is expressed during early development in a hermaphroditic fish, Kryptolebias marmoratus. Ras-related protein R-ras2 has a stop codon gain in intermediate frequency in males that is absent in females. Disruption of R-ras2 could lead to decreased cell proliferation of primordial germ cells, resulting in increased likelihood of maleness [15, 39].
Next, we evaluated SNPs that SnpEff categorized as missense mutations (Additional file 2). The first of these candidate genes is Wilms tumor protein homolog, Wt1b (14.86 Mb-14.88 Mb), which has been implicated in gonadal development and acts directly upstream of AMH, the sex-determination gene in Odontesthes hatcheri. Wt1b has also been demonstrated to bind to DNA and upregulate the sex-determination gene Sry in mammals. There is an A237V missense mutation in Wt1b that is absent in females and in intermediate frequency in males. Although our previous paper rejected Wt1b on the basis of two recombinant individuals , in light of the proposed inversion, we now believe that these individuals represented instances of natural sex reversal, not recombination.
A third candidate gene is estrogen-related receptor gamma, ERRγ (17.05-17.11 Mb). It has a R172H missense mutation within a predicted DNA-binding domain . ERRγ has been shown to be a transcriptional activator of DAX-1, and DAX-1 has been implicated as having an antagonistic effect to Sry in mammals . Therefore, a mutation in the DNA-binding domain of ERRγ could reduce DAX-1 transcription and thus have a masculinizing effect.
Growth regulation by estrogen in breast cancer 1 (GREB1) is another candidate gene (17.41-17.42 Mb) with a missense mutation. The R1775C mutation alters the side chain from a basic side chain to a polar side chain. GREB1 has been shown to be predominantly expressed within ovaries of young mice . Additionally, GREB1 has been demonstrated to be a coactivator of estrogen receptor-α . Therefore, the missense mutation in GREB1 could downregulate the expression of estrogen receptor-α, resulting in a masculinizing effect on the developing embryo.
Another potential sex-determination gene is transcription factor SOX-6 (10.22 Mb-10.30 Mb). There is a T789K missense mutation in intermediate frequency in males that is fixed in females and changes a polar side chain into a basic one. SOX-6 protein is localized to the same nuclear speckles as Sry and it has been suggested that it might play a role in sex-specific splicing in mammals .
We also evaluated gene models showing differential expression between males and females (Additional files 3 and 4). AFG3(ATPase Family Gene 3)-like protein 1 (13.72 Mb-13.73 Mb) has over a nine-fold male-biased expression. It is also on the list of SNPs with high impact coding alterations with a stop codon gain. However, a clear tie to sex-determination has yet to be elucidated.
Suppression of tumorigenicity 5 protein (10.80 Mb-10.85 Mb) and Ras association domain-containing protein 10 (11.40-11.41 Mb) were also identified for having over a three-fold male-biased expression pattern. Ras association domain family proteins have been implicated as tumor suppressors [47–49]. Therefore, upregulation of these genes could suppress primordial germ cell proliferation leading to maleness.
Lastly, it is possible that there could be Y-specific genes that were not captured in our study, because the reference genome that the reads were aligned to is a homozygous clonal XX individual.
Inversions have been well-documented in sex-chromosome evolution and are one possible mechanism for resolving sexually antagonistic selection near the novel sex-determiner through a reduction in recombination . This study revealed an 8.8 Mb block of differentiation between males and females. The variety of evidence presented here is most consistent with the presence of an inversion. The decay of genes and overall level of differentiation indicate that this region has substantially reduced recombination. We have also documented an accumulation of SNPs causing functional alterations within this region, as would be expected for a genomic region suffering both the deleterious effects of Muller’s Ratchet and accumulation of deleterious alleles hitchhiking to fixation with advantageous alleles. The transcriptome data indicates that genes inside the proposed inversion show significant enrichment for female-biased expression. These data suggest that O. niloticus has not yet evolved complete dosage compensation. Future functional studies are needed to identify the master sex-determination gene(s) within this region. Further research on cichlid sex determination will help unravel the underlying sex-determination network that underlies the rapid turnover of sex-determination mechanisms within teleosts.
Genomic DNA pools
All animal procedures were conducted in accordance with University of Maryland IACUC Protocol #R-10-73. The fish sequenced are 3rd generation descendants of fish collected from a commercial tilapia farm in Amherst, Massachusetts USA. Individuals from two related lab-raised families were sacrificed and visually inspected for the presence testes or ovaries to determine the sex of each fish. Fish with ambiguous or immature gonads were excluded from the study. DNA was extracted from fin clips using a standard phenol/chloroform protocol. To confirm the family identity we genotyped each individual for two sex-linked microsatellite markers selected from the Broad anchored tilapia assembly on linkage group 1 (MS1045 at 14.32 Mb and MS1141 at 15.53 Mb). 33 males and 20 females from family BYL078 and 25 males and 13 females from family BYL084 were selected for pooling. DNA from each individual was then quantified by Picogreen fluorescence on a BioTek FLx800 spectrophotometer and appropriate dilutions were made to ensure equal representation of each individual in the pooled samples. The pooled male (or female) DNA from each family was sheared to 500 bp using a Covaris shearer and indexed separately during library construction. Paired-end libraries for each family/sex were constructed for Illumina sequencing using the Illumina TruSeq DNA Sample Prep Kit (Illumina, San Diego, CA). The male (or female) libraries from each family were combined and each sex was sequenced in a separate lane on an Illumina HiSeq 2000. The male and female reads were deposited to NCBI with the accession numbers SRR1606298 and SRR1606304, respectively. Only reads passing the Illumina CASAVA filtering were retained.
Read qualities was checked with FASTQC . Alignments to the O. niloticus anchored reference assembly  were performed with Bowtie 2  using the --very-sensitive setting (Additional file 5). The mean alignment rate was 90.12% in males and 90.67% in females (Separate values for each family are given in Additional file 5). Read alignments were filtered for a minimum mapping quality (MAPQ) of 20 before further analysis. Insert sizes were analyzed using Picard tools CollectInsertSizeMetrics package . The aligned mean insert size was 188.76 bp (standard deviation = 44.81 bp) for males and 167.62 bp (standard deviation = 37.29 bp) for females. Variants were called using GATK .
Popoolation2  was used to calculate FST and Fisher’s exact test on allele frequency differences between the male and female pools. Initial FST results from the individually adapter-indexed families were very similar, so all subsequent analyses were performed on the combined male or female pool, including unassigned reads from the male and female lane which could not be assigned to a particular family.
A custom Perl script, Sex_SNP_finder_now.pl (available at https://github.com/Gammerdinger/sex-SNP-finder), was used to identify SNPs at intermediate frequencies in the male pools and were fixed or nearly fixed in female pools at the same position. Intermediate SNPs were defined as SNPs with a frequency between 0.3 and 0.7 within the male pool. Fixed or nearly fixed sites required a frequency less than or equal to 0.1 or greater than or equal to 0.9 within the female pool. We used a non-overlapping window of 10 kb to determine the density of these SNPs. The non-overlapping window did not include positions with coverage less than 10 reads in both sexes. The Sex_SNP_finder_now.pl script outputs a tab-delimited file with the number of SNPs per window along with an Integrative Genomics Viewer file  that lists all SNPs that were fixed or nearly fixed in one designated pool and in intermediate frequency in the other.
We used SnpEff  to identify variants predicted to alter gene function. The SnpEff output was filtered to consider only the SNPs found using Sex_SNP_finder_now.pl. SnpSift  was used to extract out SNPs with similar effects and impacts. A complete list of genes within the proposed inversion can be found in Additional file 6.
Gonads were dissected from individual larvae 28 days post-fertilization. The sex of each larvae was determined by genotyping microsatellite markers highly associated with sex. RNA from approximately 20 male or 20 female larvae was pooled and cDNA libraries were constructed using the Illumina TruSeq DNA Sample Prep Kit. Sequencing of these libraries yielded ~392 million reads for each male and female pool. Reads were aligned to the O. niloticus reference sequence with TopHat2 . NCBI RefSeq annotations were used to guide the Cufflinks  assembly (-g) and Cuffdiff was used to was used to determine FPKM values for those gene models. The results were subsequently filtered to exclude gene models whose FPKM value was less than 0.05 in both males and females. Additionally, when comparisons between FPKM of the two sexes was carried out, if the FPKM value exceeded 0.05 in one sex and was zero in the other sex, it was considered an undefined bias favoring the sex with expression. Female-biased and male-biased gene models from inside and outside the proposed inversion were counted and statistical significance was looked for using χ2 with Yates’ correction on a 2x2 contingency table. These male and female reads from the RNA-Seq experiment were deposited to NCBI with the accession numbers SRR1606274 and SRR1606273, respectively.
We wish to thank the National Science Foundation for providing computational space and support through Indiana University Pervasive Technology Institute on the National Center for Genome Analysis Support (NCGAS) Mason Computing Cluster. We also wish to thank Karen Carleton for her insightful comments and suggestions on the manuscript. This work was supported in part by the National Science Foundation under Grant Number DEB-1143920. Funding for Open Access provided by the UMD Libraries Open Access Publishing Fund.
- Charlesworth B: The evolution of sex chromosomes. Science. 1991, 251: 1030-1033. 10.1126/science.1998119.PubMedView ArticleGoogle Scholar
- Rice WR: The accumulation of sexually antagonistic genes as a selective agent promoting the evolution of reduced recombination between primitive sex chromosomes. Evolution. 1987, 41: 911-914. 10.2307/2408899.View ArticleGoogle Scholar
- Rice WR: Evolution of the Y sex in animals: Y chromosomes evolve through the degeneration of autosomes. Bioscience. 1996, 46: 331-343. 10.2307/1312947.View ArticleGoogle Scholar
- Lahn BT, Page DC: Four evolutionary strata on the human X chromosome. Science. 1999, 286: 964-967. 10.1126/science.286.5441.964.PubMedView ArticleGoogle Scholar
- Hedrick PW: Genetic drift and effective population size. Genetics of Populations. 2011, Sudbury, MA: Jones and Bartlett Publishers, 187-245. 4Google Scholar
- Green M: Muller’s Ratchet and the evolution of supernumerary chromosomes. Genome. 1990, 33: 818-824. 10.1139/g90-123.View ArticleGoogle Scholar
- Bachtrog D: Y-chromosome evolution: emerging insights into processes of Y-chromosome degeneration. Nat Rev Genetics. 2013, 14: 113-124.View ArticleGoogle Scholar
- Mank JE: Sex chromosome dosage compensation: definitely not for everyone. Trends Genet. 2013, 29: 677-683. 10.1016/j.tig.2013.07.005.PubMedView ArticleGoogle Scholar
- Vicoso B, Bachtrog D: Progress and prospects toward our understanding of the evolution of dosage compensation. Chromosome Res. 2009, 17: 585-602. 10.1007/s10577-009-9053-y.PubMed CentralPubMedView ArticleGoogle Scholar
- Heard E, Disteche CM: Dosage compensation in mammals: fine-tuning the expression of the X chromosome. Genes Dev. 2006, 20: 1848-1867. 10.1101/gad.1422906.PubMedView ArticleGoogle Scholar
- Lyon MF: Sex chromatin and gene action in the mammalian X-chromosome. Am J Hum Genet. 1962, 14: 135-148.PubMed CentralPubMedGoogle Scholar
- Bellott DW, Hughes JF, Skaletsky H, Brown LG, Pyntikova T, Cho T, Koutseva N, Zaghlul S, Graves T, Rock S, Kremitzki C, Fulton RS, Dugan S, Ding Y, Morton D, Khan Z, Lewis L, Buhay C, Wang Q, Watt J, Holder M, Lee S, Nazareth L, Rozen S, Muzny DM, Warren WC, Gibbs RA, Wilson RK, Page DC: Mammalian Y chromosomes retain widely expressed dosage-sensitive regulators. Nature. 2014, 508: 494-499. 10.1038/nature13206.PubMed CentralPubMedView ArticleGoogle Scholar
- Cortez D, Marin R, Toledo-Flores D, Froidevaux L, Liechti A, Waters P, Grützner F, Kaessman H: Origins and functional evolution of Y chromosomes across mammals. Nature. 2014, 508: 488-493. 10.1038/nature13151.PubMedView ArticleGoogle Scholar
- Koopman P, Gubbay J, Vivian N, Goodfellow P, Lovell-Badge R: Male development of chromosomally female mice transgenic for Sry. Nature. 1991, 351: 117-121. 10.1038/351117a0.PubMedView ArticleGoogle Scholar
- Kikuchi K, Hamaguchi S: Novel sex-determining genes in fish and sex chromosome evolution. Dev Dyn. 2013, 242: 339-353. 10.1002/dvdy.23927.PubMedView ArticleGoogle Scholar
- Ross JA, Urton JR, Boland J, Shapiro MD, Peichel CL: Turnover of sex chromosomes in the stickleback fishes (Gasterosteidae). PLoS Genet. 2009, 5: e1000391-10.1371/journal.pgen.1000391.PubMed CentralPubMedView ArticleGoogle Scholar
- Kocher TD: Adaptive evolution and explosive speciation: the cichlid fish model. Nat Rev Genet. 2004, 5: 288-298. 10.1038/nrg1316.PubMedView ArticleGoogle Scholar
- Ser JR, Roberts RB, Kocher TD: Multiple interacting loci control sex determination in lake Malawi cichlid fish. Evolution. 2010, 64: 486-501. 10.1111/j.1558-5646.2009.00871.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Roberts RB, Ser JR, Kocher TD: Sexual conflict resolved by invasion of a novel sex determiner in Lake Malawi cichlid fishes. Science. 2009, 326: 998-1001. 10.1126/science.1174705.PubMed CentralPubMedView ArticleGoogle Scholar
- Parnell NF, Streelman JT: Genetic interactions controlling sex and color establish the potential for sexual conflict in Lake Malawi cichlid fishes. Heredity. 2013, 110: 239-246. 10.1038/hdy.2012.73.PubMed CentralPubMedView ArticleGoogle Scholar
- Cnaani A, Lee BY, Zilberman N, Ozouf-Costaz C, Hulata G, Ron M, D’Hout A, Baroiller JF, D’Cotta H, Penman DJ, Tomasino E, Coutanceau JP, Pepey E, Shirak A, Kocher TD: Genetics of sex determination in tilapiine species. Sexual Dev. 2008, 2: 43-54. 10.1159/000117718.View ArticleGoogle Scholar
- Eshel O, Shirak A, Weller JI, Hulata G, Ron M: Linkage and physical mapping of sex region on LG23 of Nile tilapia (Oreochromis niloticus). G3. 2012, 2: 35-42. 2012.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee BY, Hulata G, Kocher TD: Two unlinked loci controlling the sex of blue tilapia (Oreochromis aureus). Heredity. 2004, 92: 543-549. 10.1038/sj.hdy.6800453.PubMedView ArticleGoogle Scholar
- Lee BY, Penman DJ, Kocher TD: Identification of a sex-determining region in Nile tilapia (Oreochromis niloticus) using bulked segregant analysis. Anim Genet. 2003, 34: 1-5. 10.1046/j.1365-2052.2003.00920.x.View ArticleGoogle Scholar
- Lee BY, Countanceau JP, Ozouf-Costaz C, D’Cotta H, Baroiller JF, Kocher TD: Genetic and physical mapping of sex-linked AFLP markers in Nile tilapia (Oreochromis niloticus). Marine Biotechnol. 2011, 13: 557-562. 10.1007/s10126-010-9326-7.View ArticleGoogle Scholar
- Ezaz MT, Harvey SC, Boonphakdee C, Teale AJ, McAndrew BJ, Penman DJ: Isolation and physical mapping of sex-linked AFLP markers in Nile tilapia (Oreochromis niloticus L.). Marine Biotechnol. 2004, 6: 435-445. 10.1007/s10126-004-3004-6.View ArticleGoogle Scholar
- Palaiokostas C, Bekaert M, Kahn MGQ, Taggart JB, Gharbi K, McAndrew BJ, Penman DJ: Mapping and validation of the major sex-determining region in Nile tilapia (Oreochromis niloticus L.) using RAD sequencing. PLoS One. 2013, 8: e68389-10.1371/journal.pone.0068389.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee BY, Lee WJ, Streelman T, Carleton KL, Howe AE, Hulata G, Slettan A, Stern JE, Terai Y, Kocher TD: A second-generation genetic linkage map of tilapia (Oreochromis spp.). Genetics. 2005, 170: 237-244. 10.1534/genetics.104.035022.PubMed CentralPubMedView ArticleGoogle Scholar
- Oliveira C, Wright JM: Molecular cytogenetic analysis of heterochromatin in the chromosomes of tilapia, Oreochromis niloticus (Teleostei: Cichlidae). Chromosome Res. 1998, 6: 205-211. 10.1023/A:1009211701829.PubMedView ArticleGoogle Scholar
- Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Ruden DM, Lu X: A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012, 6: 80-92. 10.4161/fly.19695.PubMed CentralPubMedView ArticleGoogle Scholar
- Cingolani P, Patel VM, Coon M, Nguyen T, Land SJ, Ruden DM, Lu X: Using Drosophila melanogaster as a model for genotoxic chemical mutational studies with a new program, SnpSift. Front Genet. 2012, 3: 35-PubMed CentralPubMedView ArticleGoogle Scholar
- Fisher RA: The evolution of dominance. Biol Rev. 1931, 6: 345-368. 10.1111/j.1469-185X.1931.tb01030.x.View ArticleGoogle Scholar
- Chibalina MV, Filatov DA: Plant Y chromosome degeneration is retarded by haploid purifying selection. Curr Biol. 2011, 21: 1475-1479. 10.1016/j.cub.2011.07.045.PubMedView ArticleGoogle Scholar
- Zhou Q, Bachtrog D: Chromosome-wide gene silencing initiates Y degeneration in Drosophila. Curr Biol. 2012, 22: 522-525. 10.1016/j.cub.2012.01.057.PubMedView ArticleGoogle Scholar
- Chen S, Guojie Z, Shao C, Huang Q, Liu G, Zhang P, Song W, An N, Chalopin D, Volff JN, Hong Y, Li Q, Sha Z, Zhou H, Xie M, Yu Q, Liu Y, Xiang H, Wang N, Wu K, Yang C, Zhou Q, Liao X, Yang L, Hu Q, Zhang J, Meng L, Jin L, Tian Y, Lian J, et al: Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle. Nat Genet. 2014, 46: 253-260. 10.1038/ng.2890.PubMedView ArticleGoogle Scholar
- Mair GC, Abucay JS, Skibinski DOF, Abella TA, Beardmore JA: Genetic manipulation of sex ratio for the large-scale production of all-male tilapia, Oreochromis niloticus. Can J Fish Aquat Sci. 1997, 54: 396-404. 10.1139/f96-282.View ArticleGoogle Scholar
- Lee BY, Kocher TD: Exclusion of Wilms tumour (WT1b) and ovarian cytochrome P450 aromatase (CYP19A1) as candidates for sex determination genes in Nile tilapia (Oreochromis niloticus). Anim Genet. 2007, 38: 85-86. 10.1111/j.1365-2052.2007.01563.x.PubMedView ArticleGoogle Scholar
- Massagué J: How cells read TGF-β signals. Nat Rev Mol Cell Biol. 2000, 1: 169-178. 10.1038/35043051.PubMedView ArticleGoogle Scholar
- Patmore DM, Welch S, Fulkerson PC, Wu J, Choi K, Eaves D, Kordich JJ, Collins MH, Cripe TP, Ratner N: In vivo regulation of TGF-β by R-Ras2 revealed through loss of the RasGAP protein NF1. Cancer Res. 2012, 72: 5317-5327. 10.1158/0008-5472.CAN-12-1972.PubMed CentralPubMedView ArticleGoogle Scholar
- Rhee JS, Lee YM, Raisuddin S, Lee JS: Expression of R-ras oncogenes in the hermaphroditic fish Kryptolebias marmoratus, exposed to endocrine disrupting chemicals. Comp Biochem Physiol Part C Toxicol Pharmacol. 2009, 149: 433-439. 10.1016/j.cbpc.2008.10.102.View ArticleGoogle Scholar
- Trukhina AV, Lukina NA, Wackerow-Kouzova ND, Smirnov AF: The variety of vertebrate mechanisms of sex determination. Biomed Res Int. 2013, 2013: 587460-PubMed CentralPubMedView ArticleGoogle Scholar
- MOTIF. http://www.genome.jp/tools/motif/,
- Swain A, Narvaez V, Burgoyne P, Camerino G, Lovell-Badge R: Dax1 antagonizes Sry action in mammalian sex determination. Nature. 1998, 391: 761-767. 10.1038/35799.PubMedView ArticleGoogle Scholar
- Gallardo TD, John GB, Shirley L, Contreras CM, Akbay EA, Haynie JM, Ward SE, Shidler MJ, Castrillon DH: Genomewide discovery and classification of candidate ovarian fertility genes in the mouse. Genetics. 2007, 177: 179-194. 10.1534/genetics.107.074823.PubMed CentralPubMedView ArticleGoogle Scholar
- Mohammed H, D’Santos C, Serandour AA, Ali HR, Brown GD, Atkins A, Rueda OM, Holmes KA, Theodorou V, Robinson JLL, Zwart W, Saadi A, Ross-Innes CS, Chin SF, Menon S, Stingl J, Palmieri C, Caldas C, Carroll JS: Endogenous purification reveals GREB1 as a key estrogen receptor regulatory factor. Cell Reports. 2013, 3: 342-349. 10.1016/j.celrep.2013.01.010.PubMedView ArticleGoogle Scholar
- Ohe K, Lalli E, Sassone Corsi P: A direct role of SRY and SOX proteins in pre-mRNA splicing. Proc Natl Acad Sci U S A. 2002, 99: 1146-1151. 10.1073/pnas.022645899.PubMed CentralPubMedView ArticleGoogle Scholar
- Dammann R, Li C, Yoon JH, Chin PL, Bates S, Pfeifer GP: Epigenetic inactivation of a RAS association domain family protein from the lung tumour suppressor locus 3p21.3. Nat Genet. 2000, 25: 315-319. 10.1038/77083.PubMedView ArticleGoogle Scholar
- Burbee DG, Forgacs E, Zöchbauer-Müller S, Shivakumar L, Fong K, Gao B, Randle D, Kondo M, Virmani A, Bader S, Sekido Y, Latif F, Milchgrub S, Toyooka S, Gazdar AF, Lerman MI, Zabarovsky E, White M, Minna JD: Epigenetic inactivation of RASSF1A in lung and breast cancers and malignant phenotype suppression. J Natl Cancer Inst. 2001, 93: 691-699. 10.1093/jnci/93.9.691.PubMed CentralPubMedView ArticleGoogle Scholar
- Shivakumar L, Minna J, Sakamaki T, Pestell R, White MA: The RASSF1A tumor suppressor blocks cell cycle progression and inhibits cyclin D1 accumulation. Mol Cellular Biol. 2002, 22: 4309-4318. 10.1128/MCB.22.12.4309-4318.2002.View ArticleGoogle Scholar
- FASTQC. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/,
- Brawand D, Wagner CE, Li YI, Malinsky M, Keller I, Fan S, Simakov O, Ng AY, Lim ZW, Bezault E, Turner-Maier J, Johnson J, Alcazar R, Noh HJ, Russell P, Aken B, Alföldi J, Amemiya C, Azzouzi N, Baroiller JF, Barloy-Hubler F, Berlin A, Bloomquist R, Carleton KL, Conte MA, D’Cotta HD, Eshel O, Gaffney L, Galibert F, Gante HF, et al: The genomic substrate for adaptive radiation in African cichlid fish. Nature. 2014, 513: 375-381. 10.1038/nature13726.PubMed CentralPubMedView ArticleGoogle Scholar
- Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2013, 9: 357-359.View ArticleGoogle Scholar
- Picard. http://broadinstitute.github.io/picard/,
- DePristo MA, Banks E, Poplin RE, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, McKenna A, Fennell TJ, Kernytsky AM, Sivachenko AY, Cibulskis K, Gabriel SB, Altshuler D, Daly MJ: A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011, 43: 491-498. 10.1038/ng.806.PubMed CentralPubMedView ArticleGoogle Scholar
- Kofler R, Pandey RV, Schlötterer C: PoPoolation2: identifying differentiation between populations using sequencing of pooled DNA samples (Pool-Seq). Bioinformatics. 2011, 27: 3435-3436. 10.1093/bioinformatics/btr589.PubMed CentralPubMedView ArticleGoogle Scholar
- Thorvaldsdóttir H, Robinson JT, Mesirov JP: Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2012, 14: 178-192.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14: R36-10.1186/gb-2013-14-4-r36.PubMed CentralPubMedView ArticleGoogle Scholar
- Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L: Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 31: 46-53.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.