Comparative analysis of a sex chromosome from the blackchin tilapia, Sarotherodon melanotheron

Background Inversions and other structural polymorphisms often reduce the rate of recombination between sex chromosomes, making it impossible to fine map sex-determination loci using traditional genetic mapping techniques. Here we compare distantly related species of tilapia that each segregate an XY system of sex-determination on linkage group 1. We use whole genome sequencing to identify shared sex-patterned polymorphisms, which are candidates for the ancestral sex-determination mutation. Results We found that Sarotherodon melanotheron segregates an XY system on LG1 in the same region identified in Oreochromis niloticus. Both species have higher densities of sex-patterned SNPs, as well as elevated number of ancestral copy number variants in this region when compared to the rest of the genome, but the pattern of differentiation along LG1 differs between species. The number of sex-patterned SNPs shared by the two species is small, but larger than expected by chance, suggesting that a novel Y-chromosome arose just before the divergence of the two species. We identified a shared sex-patterned SNP that alters a Gata4 binding site near Wilms tumor protein that might be responsible for sex-determination. Conclusions Shared sex-patterned SNPs, insertions and deletions suggest an ancestral sex-determination system that is common to both S. melanotheron and O. niloticus. Functional analyses are needed to evaluate shared SNPs near candidate genes that might play a role in sex-determination of these species. Interspecific variation in the sex chromosomes of tilapia species provides an excellent model system for understanding the evolution of vertebrate sex chromosomes. (Résumé d'auteur

Species-Specific Marker Discovery in Tilapia

Tilapias (family Cichlidae) are of importance in aquaculture and fisheries. Hybridisation and introgression are common within tilapia genera but are difficult to analyse due to limited numbers of species-specific genetic markers. We tested the potential of double digested restriction-site associated DNA (ddRAD) sequencing for discovering single nucleotide polymorphism (SNP) markers to distinguish between 10 tilapia species. Analysis of ddRAD data revealed 1,371 shared SNPs in the de novo-based analysis and 1,204 SNPs in the reference-based analysis. Phylogenetic trees based on these two analyses were very similar. A total of 57 species-specific SNP markers were found among the samples analysed of the 10 tilapia species. Another set of 62 species-specific SNP markers was identified from a subset of four species which have often been involved in hybridisation in aquaculture: 13 for Oreochromis niloticus, 23 for O. aureus, 12 for O. mossambicus and 14 for O. u. hornorum. A panel of 24 SNPs was selected to distinguish among these four species and validated using 91 individuals. Larger numbers of SNP markers were found that could distinguish between the pairs of species within this subset. This technique offers potential for the investigation of hybridisation and introgression among tilapia species in aquaculture and in wild populations

Structure and Sequence of the Sex Determining Locus in Two Wild Populations of Nile Tilapia

In domesticated strains of the Nile tilapia, phenotypic sex has been linked to genetic variants on linkage groups 1, 20 and 23. This diversity of sex-loci might reflect a naturally polymorphic sex determination system in Nile tilapia, or it might be an artefact arising from the process of domestication. Here, we searched for sex-determiners in wild populations from Kpandu, Lake Volta (Ghana-West Africa), and from Lake Koka (Ethiopia-East Africa) that have not been subjected to any genetic manipulation. We analysed lab-reared families using double-digest Restriction Associated DNA sequencing (ddRAD) and analysed wild-caught males and females with pooled whole-genome sequencing (WGS). Strong sex-linked signals were found on LG23 in both populations, and sex-linked signals with LG3 were observed in Kpandu samples. WGS uncovered blocks of high sequence coverage, suggesting the presence of B chromosomes. We confirmed the existence of a tandem amh duplication in LG23 in both populations and determined its breakpoints between the oaz1 and dot1l genes. We found two common deletions of ~5 kb in males and confirmed the presence of both amhY and amh∆Y genes. Males from Lake Koka lack both the previously reported 234 bp deletion and the 5 bp frameshift-insertion that creates a premature stop codon in amh∆Y

A high-resolution map of the Nile tilapia genome: a resource for studying cichlids and other percomorphs

Background: The Nile tilapia (Oreochromis niloticus) is the second most farmed fish species worldwide. It is also an important model for studies of fish physiology, particularly because of its broad tolerance to an array of environments. It is a good model to study evolutionary mechanisms in vertebrates, because of its close relationship to haplochromine cichlids, which have undergone rapid speciation in East Africa. The existing genomic resources for Nile tilapia include a genetic map, BAC end sequences and ESTs, but comparative genome analysis and maps of quantitative trait loci (QTL) are still limited. Results: We have constructed a high-resolution radiation hybrid (RH) panel for the Nile tilapia and genotyped 1358 markers consisting of 850 genes, 82 markers corresponding to BAC end sequences, 154 microsatellites and 272 single nucleotide polymorphisms (SNPs). From these, 1296 markers could be associated in 81 RH groups, while 62 were not linked. The total size of the RH map is 34,084 cR3500 and 937,310 kb. It covers 88% of the entire genome with an estimated inter-marker distance of 742 Kb. Mapping of microsatellites enabled integration to the genetic map. We have merged LG8 and LG24 into a single linkage group, and confirmed that LG16-LG21 are also merged. The orientation and association of RH groups to each chromosome and LG was confirmed by chromosomal in situ hybridizations (FISH) of 55 BACs. Fifty RH groups were localized on the 22 chromosomes while 31 remained small orphan groups. Synteny relationships were determined between Nile tilapia, stickleback, medaka and pufferfish. Conclusion:The RH map and associated FISH map provide a valuable gene-ordered resource for gene mapping and QTL studies. All genetic linkage groups with their corresponding RH groups now have a corresponding chromosome which can be identified in the karyotype. Placement of conserved segments indicated that multiple inter-chromosomal rearrangements have occurred between Nile tilapia and the other model fishes. These maps represent a valuable resource for organizing the forthcoming genome sequence of Nile tilapia, and provide a foundation for evolutionary studies of East African cichlid fishes.Additional co-authors: Thomas D Kocher, Catherine Ozouf-Costaz, Jean Francois Baroiller and Francis Galiber

Comparative physical maps derived from BAC end sequences of tilapia (Oreochromis niloticus)

Background: The Nile tilapia is the second most important fish in aquaculture. It is an excellent laboratory model, and is closely related to the African lake cichlids famous for their rapid rates of speciation. A suite of genomic resources has been developed for this species, including genetic maps and ESTs. Here we analyze BAC endsequences to develop comparative physical maps, and estimate the number of genome rearrangements, between tilapia and other model fish species. Results: We obtained sequence from one or both ends of 106,259 tilapia BACs. BLAST analysis against the genome assemblies of stickleback, medaka and pufferfish allowed identification of homologies for approximately 25,000 BACs for each species. We calculate that rearrangement breakpoints between tilapia and these species occur about every 3 Mb across the genome. Analysis of 35,000 clones previously assembled into contigs by restriction fingerprints allowed identification of longer-range syntenies. Conclusions: Our data suggest that chromosomal evolution in recent teleosts is dominated by alternate loss of gene duplicates, and by intra-chromosomal rearrangements (~one per million years). These physical maps are a useful resource for comparative positional cloning of traits in cichlid fishes. The paired BAC end sequences from these clones will be an important resource for scaffolding forthcoming shotgun sequence assemblies of the tilapia genome. (Résumé d'auteur

Genetic, abiotic and social influences on sex differentiation in cichlid fishes and the evolution of sequential hermaphroditism

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/73799/1/j.1467-2979.2005.00184.x.pd

Impact of Daily Thermocycles on Hatching Rhythms, Larval Performance and Sex Differentiation of Zebrafish

In the wild, water temperature cycles daily: it warms up after sunrise, and cools rapidly after sunset. Surprisingly, the impact of such daily thermocycles during the early development of fish remains neglected. We investigated the influence of constant vs daily thermocycles in zebrafish, from embryo development to sexual differentiation, by applying four temperature regimens: two constant (24°C and 28°C) and two daily thermocycles: 28:24°C, TC (thermophase coinciding with daytime, and cryophase coinciding with night-time) and 24:28°C, CT (opposite to TC) in a 12:12 h light:dark cycle (LD). Embryo development was temperature-dependent but enhanced at 28°C and TC. Hatching rhythms were diurnal (around 4 h after lights on), but temperature- and cycle-sensitive, since hatching occurred sooner at 28°C (48 hours post fertilization; hpf) while it was delayed at 24°C (96 hpf). Under TC, hatching occurred at 72 hpf, while under CT hatching displayed two peaks (at 70 hpf and 94 hpf). In constant light (LL) or darkness (DD), hatching rhythms persisted with tau close to 24 h, suggesting a clock-controlled "gating" mechanism. Under 28°C or TC, larvae showed the best performance (high growth and survival, and low malformations). The sex ratio was strongly influenced by temperature, as the proportion of females was higher in CT and TC (79 and 83% respectively), contrasting with 28°C and 24°C, which led to more males (83 and 76%). Ovarian aromatase (cyp19a) expression in females was highest in TC and CT (6.5 and 4.6 fold higher than at 28°C, respectively); while anti-müllerian hormone (amh) expression in males increased in testis at 24°C (3.6 fold higher compared to TC) and particularly at 28°C (14.3 fold increase). Taken together, these findings highlight the key role of environmental cycles during early development, which shaped the daily rhythms in fish embryo and larvae, and ultimately influenced sex differentiation

A high quality assembly of the Nile Tilapia (Oreochromis niloticus) genome reveals the structure of two sex determination regions

Background  Tilapias are the second most farmed fishes in the world and a sustainable source of food. Like many other fish, tilapias are sexually dimorphic and sex is a commercially important trait in these fish. In this study, we developed a significantly improved assembly of the tilapia genome using the latest genome sequencing methods and show how it improves the characterization of two sex determination regions in two tilapia species.  Results  A homozygous clonal XX female Nile tilapia (Oreochromis niloticus) was sequenced to 44X coverage using Pacific Biosciences (PacBio) SMRT sequencing. Dozens of candidate de novo assemblies were generated and an optimal assembly (contig NG50 of 3.3Mbp) was selected using principal component analysis of likelihood scores calculated from several paired-end sequencing libraries. Comparison of the new assembly to the previous O. niloticus genome assembly reveals that recently duplicated portions of the genome are now well represented. The overall number of genes in the new assembly increased by 27.3%, including a 67% increase in pseudogenes. The new tilapia genome assembly correctly represents two recentvasagene duplication events that have been verified with BAC sequencing. At total of 146Mbp of additional transposable element sequence are now assembled, a large proportion of which are recent insertions. Large centromeric satellite repeats are assembled and annotated in cichlid fish for the first time. Finally, the new assembly identifies the long-range structure of both a ~9Mbp XY sex determination region on LG1 in O. niloticus, and a ~50Mbp WZ sex determination region on LG3 in the related species O. aureus.  Conclusions  This study highlights the use of long read sequencing to correctly assemble recent duplications and to characterize repeat-filled regions of the genome. The study serves as an example of the need for high quality genome assemblies and provides a framework for identifying sex determining genes in tilapia and related fish species

Sex Reversal in Zebrafish fancl Mutants Is Caused by Tp53-Mediated Germ Cell Apoptosis

The molecular genetic mechanisms of sex determination are not known for most vertebrates, including zebrafish. We identified a mutation in the zebrafish fancl gene that causes homozygous mutants to develop as fertile males due to female-to-male sex reversal. Fancl is a member of the Fanconi Anemia/BRCA DNA repair pathway. Experiments showed that zebrafish fancl was expressed in developing germ cells in bipotential gonads at the critical time of sexual fate determination. Caspase-3 immunoassays revealed increased germ cell apoptosis in fancl mutants that compromised oocyte survival. In the absence of oocytes surviving through meiosis, somatic cells of mutant gonads did not maintain expression of the ovary gene cyp19a1a and did not down-regulate expression of the early testis gene amh; consequently, gonads masculinized and became testes. Remarkably, results showed that the introduction of a tp53 (p53) mutation into fancl mutants rescued the sex-reversal phenotype by reducing germ cell apoptosis and, thus, allowed fancl mutants to become fertile females. Our results show that Fancl function is not essential for spermatogonia and oogonia to become sperm or mature oocytes, but instead suggest that Fancl function is involved in the survival of developing oocytes through meiosis. This work reveals that Tp53-mediated germ cell apoptosis induces sex reversal after the mutation of a DNA–repair pathway gene by compromising the survival of oocytes and suggests the existence of an oocyte-derived signal that biases gonad fate towards the female developmental pathway and thereby controls zebrafish sex determination