ddRADseq reveals determinants for temperature-dependent sex reversal in Nile tilapia on LG23

Yes, polygenic sex determination is a thing!

The process of sexual development in animals is modulated by a variety of mechanisms. Some species respond to environmental cues, while, in others, sex determination is thought to be controlled by a single 'master regulator' gene. However, many animals respond to a combination of environmental cues (e.g., temperature) and genetic factors (e.g., sex chromosomes). Even among species in which genetic factors predominate, there is a continuum between monofactorial and polygenic systems. The perception that polygenic systems are rare may result from experiments that lack the statistical power to detect multiple loci. Intellectual biases against the existence of polygenic sex determination (PSD) may further arise from misconceptions about the regulation of developmental processes and a misreading of theoretical results on the stability of polygenic systems of sex determination.

A rapid approach for sex assignment by RAD-seq using a reference genome

Sex identification is a common objective in molecular ecology. While many vertebrates display sexual dimorphism, determining the sex can be challenging in certain situations, such as species lacking clear sex-related phenotypic characteristics or in studies using non-invasive methods. In these cases, DNA analyses serve as valuable tools not only for sex determination but also for validating sex assignment based on phenotypic traits. In this study, we developed a bioinformatic framework for sex assignment using genomic data obtained through GBS, and having an available closely related genome assembled at the chromosome level. Our method consists of two ad hoc indexes that rely on the different properties of the mammalian heteromorphic sex chromosomes. For this purpose, we mapped RAD-seq loci to a reference genome and then obtained missingness and coverage depth values for the autosomes and X and Y chromosomes of each individual. Our methodology successfully determined the sex of 165 fur seals that had been phenotypically sexed in a previous study and 40 sea lions sampled in a non-invasive way. Additionally, we evaluated the accuracy of each index in sequences with varying average coverage depths, with Index Y proving greater reliability and robustness in assigning sex to individuals with low-depth coverage. We believe that the approach presented here can be extended to any animal taxa with known heteromorphic XY/ZW sex chromosome systems and that it can tolerate various qualities of GBS sequencing data.

Tilapia, a good model for studying reproductive endocrinology

The Nile tilapia (Oreochromis niloticus), with a system of XX/XY sex determination, is a worldwide farmed fish with a shorter sexual maturation time than that of most cultured fish. Tilapia show a spawning cycle of approximately 14 days and can be artificially propagated in the laboratory all year round to obtain genetically all female (XX) and all male (XY) fry. Its genome sequence has been opened, and a perfect gene editing platform has been established. With a moderate body size, it is convenient for taking enough blood to measure hormone level. In recent years, using tilapia as animal model, we have confirmed that estrogen is crucial for female development because 1) mutation of star2, cyp17a1 or cyp19a1a (encoding aromatase, the key enzyme for estrogen synthesis) results in sex reversal (SR) due to estrogen deficiency in XX tilapia, while mutation of star1, cyp11a1, cyp17a2, cyp19a1b or cyp11c1 affects fertility due to abnormal androgen, cortisol and DHP levels in XY tilapia; 2) when the estrogen receptors (esr2a/esr2b) are mutated, the sex is reversed from female to male, while when the androgen receptors are mutated, the sex cannot be reversed; 3) the differentiated ovary can be transdifferentiated into functional testis by inhibition of estrogen synthesis, and the differentiated testis can be transdifferentiated into ovary by simultaneous addition of exogenous estrogen and androgen synthase inhibitor; 4) loss of male pathway genes amhy, dmrt1, gsdf causes SR with upregulation of cyp19a1a in XY tilapia. Disruption of estrogen synthesis rescues the male to female SR of amhy and gsdf but not dmrt1 mutants; 5) mutation of female pathway genes foxl2 and sf-1 causes SR with downregulation of cyp19a1a in XX tilapia; 6) the germ cell SR of foxl3 mutants fails to be rescued by estrogen treatment, indicating that estrogen determines female germ cell fate through foxl3. This review also summarized the effects of deficiency of other steroid hormones, such as androgen, DHP and cortisol, on fish reproduction. Overall, these studies demonstrate that tilapia is an excellent animal model for studying reproductive endocrinology of fish.

Polymorphism of Sex Determination Amongst Wild Populations Suggests its Rapid Turnover Within the Nile Tilapia Species

Sex-determining regions have been identified in the Nile tilapia on linkage groups (LG) 1, 20 and 23, depending on the domesticated strains used. Sex determining studies on wild populations of this species are scarce. Previous work on two wild populations, from Lake Volta (Ghana) and from Lake Koka (Ethiopia), found the sex-determining region on LG23. These populations have a Y-specific tandem duplication containing two copies of the Anti-Müllerian Hormone amh gene (named amhY and amhΔY). Here, we performed a whole-genome short-reads analysis using male and female pools on a third wild population from Lake Hora (Ethiopia). We found no association of sex with LG23, and no duplication of the amh gene. Furthermore, we found no evidence of sex linkage on LG1 or on any other LGs. Long read whole genome sequencing of a male from each population confirmed the absence of a duplicated region on LG23 in the Lake Hora male. In contrast, long reads established the structure of the Y haplotype in Koka and Kpandu males and the order of the genes in the duplicated region. Phylogenies constructed on the nuclear and mitochondrial genomes, showed a closer relationship between the two Ethiopian populations compared to the Ghanaian population, implying an absence of the LG23Y sex-determination region in Lake Hora males. Our study supports the hypothesis that the amh region is not the sex-determining region in Hora males. The absence of the Y amh duplication in the Lake Hora population reflects a rapid change in sex determination within Nile tilapia populations. The genetic basis of sex determination in the Lake Hora population remains unknown.

Genetic structure of marine and lake forms of Pacific herring Clupea pallasii

The Pacific herring (Clupea pallasii) is one of the most important species in the commercial fisheries distributed in the North Pacific Ocean and the northeastern European seas. This teleost has marine and lake ecological forms a long its distribution in the Holarctic. However, the level of genetic differentiation between these two forms is not well known. In the present study, we used ddRAD-sequencing to genotype 54 specimens from twelve wild Pacific herring populations from the Kara Sea and the Russian part of the northwestern Pacific Ocean for unveiling the genetic structure of Pacific herring. We found that the Kara Sea population is significantly distinct from Pacific Ocean populations. It was demonstrated that lake populations of Pacific herring differ from one another as well as from marine specimens. Our results show that fresh and brackish water Pacific herring, which inhabit lakes, can be distinguished as a separate lake ecological form. Moreover, we demonstrate that each observed lake Pacific herring population has its own and unique genetic legacy.

Epigenetic Regulation and Environmental Sex Determination in Cichlid Fishes

Studying environmental sex determination (ESD) in cichlids provides a phylogenetic and comparative approach to understand the evolution of the underlying mechanisms, their impact on the evolution of the overlying systems, and the neuroethology of life history strategies. Natural selection normally favors parents who invest equally in the development of male and female offspring, but evolution may favor deviations from this 50:50 ratio when environmental conditions produce an advantage for doing so. Many species of cichlids demonstrate ESD in response to water chemistry (temperature, pH, and oxygen concentration). The relative strengths of and the exact interactions between these factors vary between congeners, demonstrating genetic variation in sensitivity. The presence of sizable proportions of the less common sex towards the environmental extremes in most species strongly suggests the presence of some genetic sex-determining loci acting in parallel with the ESD factors. Sex determination and differentiation in these species does not seem to result in the organization of a final and irreversible sexual fate, so much as a life-long ongoing battle between competing male- and female-determining genetic and hormonal networks governed by epigenetic factors. We discuss what is and is not known about the epigenetic mechanism behind the differentiation of both gonads and sex differences in the brain. Beyond the well-studied tilapia species, the 2 best-studied dwarf cichlid systems showing ESD are the South American genus Apistogramma and the West African genus Pelvicachromis. Both species demonstrate male morphs with alternative reproductive tactics. We discuss the further neuroethology opportunities such systems provide to the study of epigenetics of alternative life history strategies and other behavioral variation.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.

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.

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.

miRTil: An Extensive Repository for Nile Tilapia microRNA Next Generation Sequencing Data

Nile tilapia is the third most cultivated fish worldwide and a novel model species for evolutionary studies. Aiming to improve productivity and contribute to the selection of traits of economic impact, biotechnological approaches have been intensively applied to species enhancement. In this sense, recent studies have focused on the multiple roles played by microRNAs (miRNAs) in the post-transcriptional regulation of protein-coding genes involved in the emergence of phenotypes with relevance for aquaculture. However, there is still a growing demand for a reference resource dedicated to integrating Nile Tilapia miRNA information, obtained from both experimental and in silico approaches, and facilitating the analysis and interpretation of RNA sequencing data. Here, we present an open repository dedicated to Nile Tilapia miRNAs: the “miRTil database”. The database stores data on 734 mature miRNAs identified in 11 distinct tissues and five key developmental stages. The database provides detailed information about miRNA structure, genomic context, predicted targets, expression profiles, and relative 5p/3p arm usage. Additionally, miRTil also includes a comprehensive pre-computed miRNA-target interaction network containing 4936 targets and 19,580 interactions.

Population Structure and Genetic Diversity of Nile Tilapia (Oreochromis niloticus) Strains Cultured in Tanzania

Understanding population structure and genetic diversity within and between local Nile tilapia lines cultured in Tanzania is important for sustainable aquaculture production. This study investigated the genetic structure and diversity among seven Nile tilapia populations in Tanzania (Karanga, Igunga, Ruhila, Fisheries Education and Training Agency, Tanzania Fisheries Research Institute, Kunduchi, and Lake Victoria). Double-digest restriction site-associated DNA (ddRAD) libraries were prepared from 140 individual fish (20 per population) and sequenced using an Illumina HiSeq 4000 resulting in the identification of 2,180 informative single nucleotide polymorphisms (SNPs). Pairwise F_st values revealed strong genetic differentiation between the closely related populations; FETA, Lake Victoria, and Igunga and those from TAFIRI and Karanga with values ranging between 0.45 and 0.55. Population structure was further evaluated using Bayesian model-based clustering (STRUCTURE) and discriminant analysis of principal components (DAPC). Admixture was detected among Karanga, Kunduchi, and Ruhila populations. A cross-validation approach (25% of individual fish from each population was considered of unknown origin) was conducted in order to test the efficiency of the SNP markers to correctly assign individual fish to the population of origin. The cross-validation procedure was repeated 10 times resulting in 77% of the tested individual fish being allocated to the correct population. Overall our results provide a new database of informative SNP markers for both conservation management and aquaculture activities of Nile tilapia strains in Tanzania.

Population Structure and Genetic Diversity of Nile Tilapia (Oreochromis niloticus) Strains Cultured in Tanzania

Understanding population structure and genetic diversity within and between local Nile tilapia lines cultured in Tanzania is important for sustainable aquaculture production. This study investigated the genetic structure and diversity among seven Nile tilapia populations in Tanzania (Karanga, Igunga, Ruhila, Fisheries Education and Training Agency, Tanzania Fisheries Research Institute, Kunduchi, and Lake Victoria). Double-digest restriction site-associated DNA (ddRAD) libraries were prepared from 140 individual fish (20 per population) and sequenced using an Illumina HiSeq 4000 resulting in the identification of 2,180 informative single nucleotide polymorphisms (SNPs). Pairwise Fst values revealed strong genetic differentiation between the closely related populations; FETA, Lake Victoria, and Igunga and those from TAFIRI and Karanga with values ranging between 0.45 and 0.55. Population structure was further evaluated using Bayesian model-based clustering (STRUCTURE) and discriminant analysis of principal components (DAPC). Admixture was detected among Karanga, Kunduchi, and Ruhila populations. A cross-validation approach (25% of individual fish from each population was considered of unknown origin) was conducted in order to test the efficiency of the SNP markers to correctly assign individual fish to the population of origin. The cross-validation procedure was repeated 10 times resulting in 77% of the tested individual fish being allocated to the correct population. Overall our results provide a new database of informative SNP markers for both conservation management and aquaculture activities of Nile tilapia strains in Tanzania.

Fine Mapping Using Whole-Genome Sequencing Confirms Anti-Müllerian Hormone as a Major Gene for Sex Determination in Farmed Nile Tilapia (Oreochromis niloticus L.)

Nile tilapia (Oreochromis niloticus) is one of the most cultivated and economically important species in world aquaculture. Intensive production promotes the use of monosex animals, due to an important dimorphism that favors male growth. Currently, the main mechanism to obtain all-male populations is the use of hormones in feeding during larval and fry phases. Identifying genomic regions associated with sex determination in Nile tilapia is a research topic of great interest. The objective of this study was to identify genomic variants associated with sex determination in three commercial populations of Nile tilapia. Whole-genome sequencing of 326 individuals was performed, and a total of 2.4 million high-quality bi-allelic single nucleotide polymorphisms (SNPs) were identified after quality control. A genome-wide association study (GWAS) was conducted to identify markers associated with the binary sex trait (males = 1; females = 0). A mixed logistic regression GWAS model was fitted and a genome-wide significant signal comprising 36 SNPs, spanning a genomic region of 536 kb in chromosome 23 was identified. Ten out of these 36 genetic variants intercept the anti-Müllerian (Amh) hormone gene. Other significant SNPs were located in the neighboring Amh gene region. This gene has been strongly associated with sex determination in several vertebrate species, playing an essential role in the differentiation of male and female reproductive tissue in early stages of development. This finding provides useful information to better understand the genetic mechanisms underlying sex determination in Nile tilapia.

Mismatches between the genetic and phenotypic sex in the wild Kou population of Nile tilapia Oreochromis niloticus

Sex determination and sex chromosomes can be very diverse between teleost species. The group of tilapias shows a polymorphism in sex determination not only between closely related species but also between domestic strains within a species. In the Nile tilapia, the major effect genes and therefore the Y chromosome have been located on either linkage group 1 (LG1) or LG23 depending on the strains. In a Japanese strain, the sex determinant of LG23 (the amhY gene) has been identified as a duplicated amh (anti-Müllerian hormone) gene, with its gametolog found on the X chromosome (amhX). AmhY is located in tandem with the amhΔY gene (a truncated form) on the Y chromosome. X and Y chromosome markers based on the amh genes have been validated only on a few domestic strains but not in wild populations. Here, we used four of these markers in order to examine (1) the possible variation in sex determination of a wild population of Nile tilapia living in Lake Kou (Burkina Faso), (2) putative polymorphisms for these amh copies and (3) the existence of sex reversed individuals in the wild. Our genotyping of 91 wild Kou individuals with the amh sex-diagnostic markers of LG23 showed that while phenotypic females were all XX, phenotypic males were either XY or XX. Progeny testing of eight of these XX males revealed that one of these males consistently sired all-female progenies, suggesting that it is a wild sex reversed male (which could result from high temperature effects). The other XX males gave balanced sex ratios, suggesting that sex is controlled by another locus (possibly on another LG) which may be epistatically dominant over the LG23 locus. Finally, identification of unexpected amh genotypes was found for two individuals. They produced either balanced or female-biased sex ratios, depending on the breeder with whom they were crossed, suggesting possible recombination between the X and the Y chromosomes.

Genome-wide RAD sequencing to identify a sex-specific marker in Chinese giant salamander Andrias davidianus

Background

Chinese giant salamander Andrias davidianus is an endangered species. The success of artificial breeding provides a useful way to protect this species. However, the method to identify the sex and mechanism of sex determination were unclear which hinder the improvement of the artificial breeding. Detection of a sex specific marker provides an effective approach to identify genetic sex and investigate the sex determination mechanism.

Results

We used restriction-site-associated DNA (RAD) sequencing to isolate a sex-specific genetic marker in A. davidianus to expand knowledge of the sex determination mechanism. Four male and four female specimens were subjected to RAD sequencing, which generated 934,072,989 reads containing approximately 134.4 Gb of sequences. The first round of comparison of the assembled sequence against the opposite sex raw reads revealed 19,097 female and 17,994 male unmatched sequences. Subsequently, 19,097 female sequences were subjected to a BLAST search against male genomic data, which revealed 308 sequences unmapped to the male genome. One hundred of these were randomly selected and validated by PCR in five male and five female specimens, and four putative sex-specific sequences were produced. Further validation was performed by PCR in another 24 females and 24 males, and all female individuals exhibited the expected specific bands, while the males did not. To apply the sex-specific marker, three specimens reversed from genetic female to physiological male were found in a group exposed to elevated temperature, and 13 individuals reversed from genetic male to physiological female were obtained in a 17β-estradiol exposed group.

Conclusion

This is the first report of a sex-specific marker in A. davidianus and may have potential for elucidation of its sex determination mechanism and, hence, its conservation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5771-5) contains supplementary material, which is available to authorized users.

Temperature preference of Nile tilapia (Oreochromis niloticus) juveniles induces spontaneous sex reversal

Nile tilapia (Oreochromis niloticus) is an African freshwater fish that displays a genetic sex determination system (XX|XY) where high temperatures (above 32°C to 36.5°C) induce masculinization. In Nile tilapia, the thermosensitive period was reported from 10 to 30 days post fertilization. In their natural environment, juveniles may encounter high temperatures that are above the optimal temperature for growth (27-30°C). The relevance of the thermal sex reversal mechanism in a natural context remains unclear. The main objective of our study is to determine whether sexually undifferentiated juveniles spontaneously prefer higher, unfavorable temperatures and whether this choice skews the sex ratio toward males. Five full-sib progenies (from 100% XX crosses) were subjected to (1) a horizontal three-compartment thermal step gradient (thermal continuum 28°C- 32°C- 36.5°C) during the thermosensitive period, (2) a control continuum (28°C- 28°C- 28°C) and (3) a thermal control tank (36.5°C). During the first days of the treatment, up to an average of 20% of the population preferred the masculinizing compartment of the thermal continuum (36.5°C) compared to the control continuum. During the second part of the treatment, juveniles preferred the lower, nonmasculinizing 32°C temperature. This short exposure to higher temperatures was sufficient to significantly skew the sex ratio toward males, compared to congeners raised at 28°C (from 5.0 ± 6.7% to 15.6 ± 16.5% of males). The proportion of males was significantly different in the thermal continuum, thermal control tank and control continuum, and it was positively correlated among populations. Our study shows for the first time that Nile tilapia juveniles can choose a masculinizing temperature during a short period of time. This preference is sufficient to induce sex reversal to males within a population. For the first time, behavior is reported as a potential player in the sex determination mechanism of this species.