Quantitative trait loci on LGs 9 and 14 affect the reproductive interaction between two Oreochromis species, O. niloticus and O. aureus

Absence of Figla-like Gene Is Concordant with Femaleness in Cichlids Harboring the LG1 Sex-Determination System

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.

Gene Variant of Barrier to Autointegration Factor 2 (Banf2w) Is Concordant with Female Determination in Cichlids

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.

All-male production by marker-assisted selection for sex determining loci of admixed Oreochromis niloticus and Oreochromis aureus stocks

Crossing Oreochromis niloticus (On) females with Oreochromis aureus (Oa) males results in all-male 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 sex-determination (SD) of the respective purebred species. Here, we test this approach in a data set of 160 On × 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 × 10^-32 ) and for LG3-UNH168 × LG23-UNH898 interaction (P < 1 × 10^-17 ; R²  = 0.98). Furthermore, a maternal effect for LG3-UNH168 (P < 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.

High-quality chromosome-level genomes of two tilapia species reveal their evolution of repeat sequences and sex chromosomes

Tilapias are one of the most farmed fishes that are coined as "aquatic chicken" by the food industry. Nile tilapia and blue tilapia exhibit very recent transition of sex chromosome systems since their divergence approximately five million years ago, making them a great model for elucidating the molecular and evolutionary mechanisms of sex chromosome turnovers. Studies of their sex-determining pathways are also critical for developing genetic sex control in aquaculture. We report here the newly produced genomes of Nile tilapia and blue tilapia that integrate long-read sequencing and chromatin conformation data. The two nearly complete genomes have anchored over 97% of the sequences into linkage groups (LGs), and assembled majorities of complex repetitive regions including telomeres, centromeres and rDNA clusters. In particular, we inferred two episodes of repeat expansion at LG3 respectively in the ancestor of cichlids and that of tilapias. The consequential large heterochromatic region concentrated at one end of LG3 comprises tandem arrays of mRNA and small RNA genes, among which we have identified a candidate female determining gene Paics in blue tilapia. Paics shows female-specific patterns of single-nucleotide variants, copy numbers and expression patterns in gonads during early gonadogenesis. Our work provides a very important genomic resource for functional studies of cichlids, and suggested that unequal distribution of repeat content that impacts the local recombination rate might make some chromosomes more likely to become sex chromosomes.

A duplication of the Anti-Müllerian hormone gene is associated with genetic sex determination of different Oreochromis niloticus strains

Sex determination (SD) mechanisms are ancient and conserved, yet much diversity is exhibited in primary sex-determining signals that trigger male or female development. In O. niloticus, SD is associated with a male-specific locus on linkage group (LG) 23 which harbors the Y-linked Anti-Müllerian hormone (amh) gene, and a truncated duplication, denoted amhΔy. We have evaluated the possible role of identified indels and SNPs in the amh gene on SD, based on conservation in different O. niloticus strains. A fluorescent assay for the detection of a 5 bp insertion in amhΔy exon VI, efficiently discriminated between XX, XY, and YY genotypes. Concordance rate between amhΔy and sex varied in six Oreochromis strains, from 100% (Ghana) through 90% (Swansea) to 85% (Thai-Chitralada). The association of amhΔy with sex was found to be conserved in all tested O. niloticus strains, and thus supports its key role in SD. However, the previously identified missense SNP (C/T) in amh exon II was found only in the Swansea strain, thus excluding its candidacy for the causal variation of SD across all strains. Effects of markers on LGs 1, 3, and 23 (amhΔy) fully explained sex distribution in one Thai-Chitralada family (R2 = 1.0), whereas in another family only the major effect of LG23 (amhΔy) was significant (R2 = 0.37). Thus, amhΔy on LG23 is associated with genetic SD, either as a single causal gene in different O. niloticus strains, or in combination with segregating genes on LGs 1 and 3 in the Thai-Chitralada hybrid strain.