Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

We identified DNA markers linked to sex determining genes in six closely related species of tilapiine fishes. The mode of sex determination differed among species. In Oreochromis karongae and Tilapia mariae the sex-determining locus is on linkage group (LG) 3 and the female is heterogametic (WZ-ZZ system). In O. niloticus and T. zillii the sex-determining locus is on LG1 and the male is heterogametic (XX-XY system). A more complex pattern was observed in O. aureus and O. mossambicus, in which markers on both LG1 and LG3 were associated with sex. We found evidence for sex-linked lethal effects on LG1, as well as interactions between loci in the two linkage groups. Comparison of genetic and physical maps demonstrated a broad region of recombination suppression harboring the sex-determining locus on LG3. Sex-specific recombination suppression was found in the female heterogametic sex. Sequence analysis showed the accumulation of repetitive elements in this region. Phylogenetic analysis suggests that at least two transitions in the mode of sex determination have occurred in this clade. This variation in sex determination mechanisms among closely related species makes tilapias an excellent model system for studying the evolution of sex chromosomes in vertebrates.

BackgroundDense single nucleotide polymorphism (SNP) genotyping arrays provide extensive information on polymorphic variation across the genome of species of interest. Such information can be used in studies of the genetic architecture of quantitative traits and to improve the accuracy of selection in breeding programs. In Atlantic salmon (Salmo salar), these goals are currently hampered by the lack of a high-density SNP genotyping platform. Therefore, the aim of the study was to develop and test a dense Atlantic salmon SNP array.ResultsSNP discovery was performed using extensive deep sequencing of Reduced Representation (RR-Seq), Restriction site-Associated DNA (RAD-Seq) and mRNA (RNA-Seq) libraries derived from farmed and wild Atlantic salmon samples (n = 283) resulting in the discovery of > 400 K putative SNPs. An Affymetrix Axiom® myDesign Custom Array was created and tested on samples of animals of wild and farmed origin (n = 96) revealing a total of 132,033 polymorphic SNPs with high call rate, good cluster separation on the array and stable Mendelian inheritance in our sample. At least 38% of these SNPs are from transcribed genomic regions and therefore more likely to include functional variants. Linkage analysis utilising the lack of male recombination in salmonids allowed the mapping of 40,214 SNPs distributed across all 29 pairs of chromosomes, highlighting the extensive genome-wide coverage of the SNPs. An identity-by-state clustering analysis revealed that the array can clearly distinguish between fish of different origins, within and between farmed and wild populations. Finally, Y-chromosome-specific probes included on the array provide an accurate molecular genetic test for sex.ConclusionsThis manuscript describes the first high-density SNP genotyping array for Atlantic salmon. This array will be publicly available and is likely to be used as a platform for high-resolution genetics research into traits of evolutionary and economic importance in salmonids and in aquaculture breeding programs via genomic selection.

Sex in Oreochromis niloticus (Nile tilapia) is principally determined by an XX/XY locus but other genetic and environmental factors also influence sex ratio. Restriction Associated DNA (RAD) sequencing was used in two families derived from crossing XY males with females from an isogenic clonal line, in order to identify Single Nucleotide Polymorphisms (SNPs) and map the sex-determining region(s). We constructed a linkage map with 3,802 SNPs, which corresponded to 3,280 informative markers, and identified a major sex-determining region on linkage group 1, explaining nearly 96% of the phenotypic variance. This sex-determining region was mapped in a 2 cM interval, corresponding to approximately 1.2 Mb in the O. niloticus draft genome. In order to validate this, a diverse family (4 families; 96 individuals in total) and population (40 broodstock individuals) test panel were genotyped for five of the SNPs showing the highest association with phenotypic sex. From the expanded data set, SNPs Oni23063 and Oni28137 showed the highest association, which persisted both in the case of family and population data. Across the entire dataset all females were found to be homozygous for these two SNPs. Males were heterozygous, with the exception of five individuals in the population and two in the family dataset. These fish possessed the homozygous genotype expected of females. Progeny sex ratios (over 95% females) from two of the males with the “female” genotype indicated that they were neomales (XX males). Sex reversal induced by elevated temperature during sexual differentiation also resulted in phenotypic males with the “female” genotype. This study narrows down the region containing the main sex-determining locus, and provides genetic markers tightly linked to this locus, with an association that persisted across the population. These markers will be of use in refining the production of genetically male O. niloticus for aquaculture.

BackgroundTilapias 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.ResultsA 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 recent vasa gene 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. ConclusionsThis 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-017-3723-5) contains supplementary material, which is available to authorized users.

No abstract

A brain aromatase gene was identified from the Nile tilapia Oreochromis niloticus. The cDNA sequence of this gene differed from that of the ovarian aromatase gene previously reported from this species. Tissue specific expression for both brain and ovarian aromatase genes was examined in the tissues of adult tilapia. Brain aromatase mRNA was expressed in the brain, kidney, eye, ovary, and testis, but not in the liver and spleen. Ovarian aromatase mRNA was expressed in the brain, spleen, ovary, and testis but not in the eye, kidney, and liver. Differential aromatase gene expression between the sexes was investigated in all-male (XY) and all-female (XX) groups of tilapia fry from fertilisation throughout the sexual differentiation period. Semi-quantitative RT-PCR analysis revealed that the initiation of expression of both aromatase genes lay between 3 and 4 dpf (days post fertilisation) in both sexes. The level of brain aromatase mRNA gradually increased throughout the period studied with little difference between the sexes. This contrasted with marked sexual dimorphism of ovarian aromatase mRNA expression. In females, the expression level was maintained or increased gradually throughout ontogeny, while the level in males was dramatically down-regulated between 15 and 27 dpf. Subsequently, the level of ovarian aromatase mRNA expression fluctuated slightly in both sexes, with the expression in females always being higher than in males. These findings clearly suggest that ovarian aromatase plays a decisive role in sexual differentiation in this species and that this is achieved by down-regulation of the expression of this gene in males. Mol. Reprod. Dev. 59: 359-370, 2001.

Sex determination in the Nile tilapia (Oreochromis niloticus) is thought to be an XX-XY (male heterogametic) system controlled by a major gene. We searched for DNA markers linked to this major locus using bulked segregant analysis. Ten microsatellite markers belonging to linkage group 8 were found to be linked to phenotypic sex. The putative Y-chromosome alleles correctly predict the sex of 95% of male and female individuals in two families. Our results suggest a major sex-determining locus within a few centimorgans of markers UNH995 and UNH104. A third family from the same population showed no evidence for linkage of this region with phenotypic sex, indicating that additional genetic and/or environmental factors regulate sex determination in some families. These markers have immediate utility for studying the strength of different Y chromosome alleles, and for identifying broodstock carrying one or more copies of the Y haplotype.