Effective farming of tilapia requires all-male culture, characterized by uniformity and high growth rate. Males of O. aureus (Oa) and females of O. niloticus (On) produce all-male offspring, but there is a behavioral reproductive barrier between the two species that prevents mass production. In crosses between Oa and On broodstocks, few hybrid females are attracted to the Oa male nests (denoted responders), and if they harbor the On alleles for the sex determination (SD) sites on linkage groups (LGs) 1, 3, and 23, all-male progeny are produced. Yet, without controlling for the alleles underlying SD, the parental stocks gradually lose their capability for all-male production. Hypothesizing that marker-assisted selection for female responders would allow production of sustainable broodstocks, we applied genotyping-by-sequencing to generate 4983 informative SNPs from 13 responding and 28 non-responding females from two full-sib families. Accounting for multiple comparisons in a genome-wide association study, seven SNPs met a false discovery rate of 0.061. Lowest nominal probabilities were on LGs 9 and 14, for which microsatellite DNA markers were designed within the candidate genes PTGDSL and CASRL, respectively. By increasing the sample size to 22 responders and 47 non-responders and by genotyping additional established microsatellites, we confirmed the association of these LGs with female responsiveness. The combined effects of microsatellites GM171 and CARSL-LOC100690618 on LGs 9 and 14 explained 37% of the phenotypic variance of reproductive interaction (p < 0.0001). Based on these findings, we propose a strategy for mass production of all-male tilapia hybrids through selection for genomic loci affecting SD and female responsiveness.

Genetic parameters and selection responses were obtained for harvest body weight of blue tilapia (Oreochromis aureus) from data collected over three generations in a selected population. A total of 18 194 records representing 186 sires and 201 dams were used in the analysis. Within generation heritability estimates for harvest body weight ranged from 0.18 to 0.58. When data from more than one generation were included in the analysis, heritability estimates became more stable (0.33-0.40) and it was 0.33 when all data were included in the analysis. The common full-sib effect accounted for 10% of the phenotypic variance in the full data set. Heritability for survival from stocking to harvest was estimated at 0.01 and 0.09 in actual units (fitting an animal model) and in the logit (sire model) scale respectively. The genetic correlation between harvest body weight and survival was 0.22 and not significantly different from zero. The total selection response for harvest body weight over the three generations of selection measured as the difference between least-squares means of selected and control lines was 17.7%. The corresponding figure when response was measured as the difference between mean breeding values of selected and control lines was 19.6%. The average inbreeding coefficient was 0.003 after three generations of selection. These results indicate that there are good prospects for the genetic improvement of harvest body weight in blue tilapia.

Crossing Oreochromis niloticus (On) females with Oreochromis aureus (Oa) males results in allmale progeny that are essential for effective tilapia farming. However, a reproductive barrier between these species limits mating and mass-fry production. One approach to overcoming this barrier is to select parental stocks of mixed genetic backgrounds, which allow interspecific reproductive recognition, while closely maintaining the genetic profiles for sexdetermination (SD) of the respective purebred species. Here, we test this approach in a data set of 160 On 9 Oa spawns of 109 male and 100 female parents randomly collected from admixed stocks, and genotyped for microsatellite markers representing the known SD loci on linkage groups (LGs) 1, 3, and 23. Following crossbreeding, the most significant paternal effects on male proportions in progeny were found for LG1-BYL018 (P < 2 9 10 À32 ) and for LG3-UNH168 9 LG23-UNH898 interaction (P < 1 9 10 À17 ; R 2 = 0.98). Furthermore, a maternal effect for LG3-UNH168 (P < 9 9 10 À7 ) was associated with low female proportions in progeny (<7%), indicating a non-Mendelian effect on SD. Eighty-four males (77%) and 30 females (30%) were selected as parents, based on their genetic profiles for the SD loci that were associated with male production. Of these, 51 of 53 crosses produced all-male progeny, while two crosses had low female proportions in their progeny (<4%). This suggests that selection could be improved using the causative sequence variation underlying SD on LG3, since the large non-recombining block of the SD region in purebred Oa readily breaks down in hybrids. Nevertheless, marker-assisted selection for sex determining loci of admixed parental stocks may be used for all-male production.

Oreochromis fishes exhibit variability of sex-determination (SD) genes whose characterization contributes to understanding of the sex differentiation network, and to effective tilapia farming, which requires all-male culture. However, O. niloticus (On) amh is the only master-key regulator (MKR) of SD that has been mapped (XY/XX SD-system on LG23). In O. aureus (Oa), LG3 controls a WZ/ZZ SD-system that has recently been delimited to 9.2 Mbp, with an embedded interval rich with female-specific variation, harboring two paics genes and banf2. Developing genetic markers within this interval and using a hybrid Oa stock that demonstrates no recombination repression in LG3, we mapped the critical SD region to 235 Kbp on the orthologous On physical map (p < 1.5 × 10−26). DNA-seq assembly and peak-proportion analysis of variation based on Sanger chromatograms allowed the characterization of copy-number variation (CNV) of banf2. Oa males had three exons capable of encoding 90-amino-acid polypeptides, yet in Oa females, we found an extra copy with an 89-amino-acid polypeptide and three non-conservative amino acid substitutions, designated as banf2w. CNV analysis suggested the existence of two to five copies of banf2 in diploidic Cichlidae. Disrupting the Hardy–Weinberg equilibrium (p < 4.2 × 10−3), banf2w was concordant with female determination in Oa and in three cichlids with LG3 WZ/ZZ SD-systems (O. tanganicae, O. hornorum and Pelmatolapia mariae). Furthermore, exclusive RNA-seq expression in Oa females strengthened the candidacy of banf2w as the long-sought LG3 SD MKR. As banf genes mediate nuclear assembly, chromatin organization, gene expression and gonad development, banf2w may play a fundamental role inducing female nucleus formation that is essential for WZ/ZZ SD.

Oreochromis niloticus has been used as a reference genome for studies of tilapia sex determination (SD) revealing segregating genetic loci on linkage groups (LGs) 1, 3, and 23. The master key regulator genes (MKR) underlying the SD regions on LGs 3 and 23 have been already found. To identify the MKR in fish that segregate for the LG1 XX/XY SD-system, we applied short variant discovery within the sequence reads of the genomic libraries of the Amherst hybrid stock, Coptodon zillii and Sarotherodon galilaeus, which were aligned to a 3-Mbp-region of the O. aureus genome. We obtained 66,372 variants of which six were concordant with the XX/XY model of SD and were conserved across these species, disclosing the male specific figla-like gene. We further validated this observation in O. mossambicus and in the Chitralada hybrid stock. Genome alignment of the 1252-bp transcript showed that the figla-like gene’s size was 2664 bp, and that its three exons were capable of encoding 99 amino acids including a 45-amino-acid basic helix–loop–helix domain that is typical of the ovary development regulator—factor-in-the-germline-alpha (FIGLA). In Amherst gonads, the figla-like gene was exclusively expressed in testes. Thus, the figla-like genomic presence determines male fate by interrupting the female developmental program. This indicates that the figla-like gene is the long-sought SD MKR on LG1.