Identification of a sex-determining region in Nile tilapia (Oreochromis niloticus) using bulked segregant analysis

Diversity of Sex Chromosomes in Vertebrates: Six Novel Sex Chromosomes in Basal Haplochromines (Teleostei: Cichlidae)

African cichlid fishes are known for their high rates of phenotypic evolution. A rapid rate of diversification is apparent also in the diversity of their sex chromosomes. To date, sex determiners have been identified on 18 of 22 chromosomes in the standard karyotype. Here, we use whole-genome sequencing to characterize the sex chromosomes of seven populations of basal haplochromines, focusing on the genus Pseudocrenilabrus. We identify six new sex chromosome systems, including the first report of a cichlid sex-determining system on linkage group 12. We then quantify the rates and patterns of sex chromosome turnover in this clade. Finally, we test whether some autosomes become sex chromosomes in East African cichlids more often than expected by chance.

Effect of high temperatures on sex ratio and differential expression analysis (RNA-seq) of sex-determining genes in Oreochromis niloticus from different river basins in Benin

The high temperature sex reversal process leading to functional phenotypic masculinization during development has been widely described in Nile tilapia (Oreochromis n iloticus) under laboratory or aquaculture conditions and in the wild. In this study, we selected five wild populations of O. niloticus from different river basins in Benin and produced twenty full-sib families of mixed-sex (XY and XX) by natural reproduction. Progenies were exposed to room temperature or high (36.5°C) temperatures between 10 and 30 days post-fertilization (dpf). In control groups, we observed sex ratios from 40% to 60% males as expected, except for 3 families from the Gobé region which showed a bias towards males. High temperature treatment significantly increased male rates in each family up to 88%. Transcriptome analysis was performed by RNA-sequencing (RNA-seq) on brains and gonads from control and treated batches of six families at 15 dpf and 40 dpf. Analysis of differentially expressed genes, differentially spliced genes, and correlations with sex reversal was performed. In 40 dpf gonads, genes involved in sex determination such as dmrt1, cyp11c1, amh, cyp19a1b, ara, and dax1 were upregulated. In 15 dpf brains, a negative correlation was found between the expression of cyp19a1b and the reversal rate, while at 40 dpf a negative correlation was found between the expression of foxl2, cyp11c1, and sf1 and positive correlation was found between dmrt1 expression and reversal rate. Ontology analysis of the genes affected by high temperatures revealed that male sex differentiation processes, primary male sexual characteristics, autophagy, and cilium organization were affected. Based on these results, we conclude that sex reversal by high temperature treatment leads to similar modifications of the transcriptomes in the gonads and brains in offspring of different natural populations of Nile tilapia, which thus may activate a common cascade of reactions inducing sex reversal in progenies.

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.

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.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Candidate powdery mildew resistance gene in wheat landrace cultivar Hongyoumai discovered using SLAF and BSR-seq

BACKGROUND

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is an important disease affecting wheat production. Planting resistant cultivars is an effective, safe, and economical method to control the disease. Map construction using next-generation sequencing facilitates gene cloning based on genetic maps and high-throughput gene expression studies. In this study, specific-locus amplified fragment sequencing (SLAF) was used to analyze Huixianhong (female parent), Hongyoumai (male parent) and two bulks (50 homozygous resistant and 50 susceptible F_2:3 segregating population derived from Huixianhong × Hongyoumai to determine a candidate gene region for resistance to powdery mildew on the long arm of chromosome 7B in wheat landrace Hongyoumai. Gene expressions of candidate regions were obtained using bulked segregant RNA-seq in 10 homozygous resistant and 10 susceptible progeny inoculated by Bgt.. Candidate genes were obtained using homology-based cloning in two parents.

RESULTS

A 12.95 Mb long candidate region in chromosome 7BL was identified, and five blocks in SLAF matched the scaffold of the existing co-segregation marker Xmp1207. In the candidate region, 39 differentially expressed genes were identified using RNA-seq, including RGA4 (Wheat_Chr_Trans_newGene_16173)-a disease resistance protein whose expression was upregulated in the resistant pool at 16 h post inoculation with Bgt. Quantitative reverse transcription (qRT)-PCR was used to further verify the expression patterns in Wheat_Chr_Trans_newGene_16173 that were significantly different in the two parents Hongyoumai and Huixianhong. Two RGA4 genes were cloned based on the sequence of Wheat_Chr_Trans_newGene_16173, respectively from two parent and there was one amino acid mutation: S to G in Huixianhong on 510 loci.

CONCLUSION

The combination of SLAF and BSR-seq methods identified a candidate region of pmHYM in the chromosome 7BL of wheat landrace cultivar Hongyoumai. Comparative analysis between the scaffold of co-segregating marker Xmp1207 and SLAF-seq showed five matching blocks. qRT-PCR showed that only the resistant gene Wheat_Chr_Trans_newGene_16173 was significantly upregulated in the resistant parent Hongyoumai after inoculation with Bgt, and gene cloning revealed a difference in one amino acid between the two parent genes, indicating it was involved in the resistance response and may be the candidate resistance gene pmHYM.

A Chromosome-Level Genome Assembly of Mozambique Tilapia (Oreochromis mossambicus) Reveals the Structure of Sex Determining Regions

The Mozambique tilapia (Oreochromis mossambicus) is a fascinating taxon for evolutionary and ecological research. It is an important food fish and one of the most widely distributed tilapias. Because males grow faster than females, genetically male tilapia are preferred in aquaculture. However, studies of sex determination and sex control in O. mossambicus have been hindered by the limited characterization of the genome. To address this gap, we assembled a high-quality genome of O. mossambicus, using a combination of high coverage of Illumina and Nanopore reads, coupled with Hi-C and RNA-Seq data. Our genome assembly spans 1,007 Mb with a scaffold N50 of 11.38 Mb. We successfully anchored and oriented 98.6% of the genome on 22 linkage groups (LGs). Based on re-sequencing data for male and female fishes from three families, O. mossambicus segregates both an XY system on LG14 and a ZW system on LG3. The sex-patterned SNPs shared by two XY families narrowed the sex determining regions to ∼3 Mb on LG14. The shared sex-patterned SNPs included two deleterious missense mutations in ahnak and rhbdd1, indicating the possible roles of these two genes in sex determination. This annotated chromosome-level genome assembly and identification of sex determining regions represents a valuable resource to help understand the evolution of genetic sex determination in tilapias.

Origin of a Giant Sex Chromosome

Chromosome size and morphology vary within and among species, but little is known about the proximate or ultimate causes of these differences. Cichlid fish species in the tribe Oreochromini share an unusual giant chromosome that is ∼3 times longer than the other chromosomes. This giant chromosome functions as a sex chromosome in some of these species. We test two hypotheses of how this giant sex chromosome may have evolved. The first hypothesis proposes that it evolved by accumulating repetitive elements as recombination was reduced around a dominant sex determination locus, as suggested by canonical models of sex chromosome evolution. An alternative hypothesis is that the giant sex chromosome originated via the fusion of an autosome with a highly repetitive B chromosome, one of which carried a sex determination locus. We test these hypotheses using comparative analysis of chromosome-scale cichlid and teleost genomes. We find that the giant sex chromosome consists of three distinct regions based on patterns of recombination, gene and transposable element content, and synteny to the ancestral autosome. The WZ sex determination locus encompasses the last ∼105 Mb of the 134-Mb giant chromosome. The last 47 Mb of the giant chromosome shares no obvious homology to any ancestral chromosome. Comparisons across 69 teleost genomes reveal that the giant sex chromosome contains unparalleled amounts of endogenous retroviral elements, immunoglobulin genes, and long noncoding RNAs. The results favor the B chromosome fusion hypothesis for the origin of the giant chromosome.

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.

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.

Sex determination in the GIFT strain of tilapia is controlled by a locus in linkage group 23

Background

Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The ‘genetically improved farmed tilapia’ (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins.

Results

A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284).

Conclusions

Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

AMH and AMHR2 mutations: A spectrum of reproductive phenotypes across vertebrate species

Anti-Müllerian hormone (AMH) is a member of the Transforming Growth Factor-β family of secreted signaling proteins. AMH is expressed in Sertoli cells of the fetal and adult testes and granulosa cells of the postnatal ovary. AMH is required for the regression of the Müllerian ducts in mammalian fetuses during male differentiation. AMH signals through its Type II receptor, AMHR2. AMHR2 is expressed in mesenchyme adjacent to the Müllerian ducts, and in Sertoli, Leydig, and granulosa cells. Although AMH and AMHR2 genes have been identified in numerous vertebrate species, spontaneous or engineered mutations or variants have been found or created in only a few mammals and teleost fishes. AMH or AMHR2 mutations in mammals lead to the development of Persistent Müllerian Duct Syndrome (PMDS), a recessive condition in which affected males are fully virilized but retain Müllerian duct-derived tissues, including a uterus and oviducts, and in human and dog, undescended testes. Amh mutant female mice had accelerated ovarian primordial follicle recruitment, suggesting a role for AMH in regulating germ cells. amh and amhr2 mutations have also been experimentally generated in various teleost fishes. Depending on the fish species, loss of AMH signaling results in infertility, germ cell tumors, or male-to-female sex reversal. Here we compare the spectrum of phenotypes caused by AMH and AMHR2 mutations in a variety of vertebrate species. There are both common and unique phenotypes between species, highlighting the range of biological processes regulated by AMH signaling.

Novel Sex Chromosomes in 3 Cichlid Fishes from Lake Tanganyika

African cichlids are well known for their adaptive radiations, but it is now apparent that they also harbor an extraordinary diversity of sex chromosome systems. In this study, we sequenced pools of males and females from species in 3 different genera of cichlids from Lake Tanganyika. We then searched for regions that were differentiated following the patterns expected for sex chromosomes. We report 2 novel sex chromosomes systems, an XY system on LG19 in Tropheus sp. "black" and a ZW system on LG7 in Hemibates stenosoma. We also identify a ZW system on LG5 in Cyprichromis leptosoma that may be convergent with a system previously described in Lake Malawi cichlids. Our data also identify candidate single nucleotide polymorphisms for the blue/yellow tail color polymorphism observed among male C. leptosoma.

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.

QTL Mapping for Red Blotches in Malaysia Red Tilapia (Oreochromis spp.)

Body color is an interesting economic trait in fish. Red tilapia with red blotches may decrease its commercial values. Conventional selection of pure red color lines is a time-consuming and labor-intensive process. To accelerate selection of pure lines through marker-assisted selection, in this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) technology was applied to genotype a full-sib mapping family of Malaysia red tilapia (Oreochromis spp.) (N = 192). Genome-wide significant quantitative trait locus (QTL)-controlling red blotches were mapped onto two chromosomes (chrLG5 and chrLG15) explaining 9.7% and 8.2% of phenotypic variances by a genome-wide association study (GWAS) and linkage-based QTL mapping. Six SNPs from the chromosome chrLG5 (four), chrLG15 (one), and unplaced supercontig GL831288-1 (one) were significantly associated to the red blotch trait in GWAS analysis. We developed nine microsatellite markers and validated significant correlations between genotypes and blotch data (p < 0.05). Our study laid a foundation for exploring a genetic mechanism of body colors and carrying out genetic improvement for color quality in tilapia.

Constructing High-Density Genetic Maps and Developing Sexing Markers in Northern Snakehead (Channa argus)

High-density genetic maps are essential for mapping QTL, improving genome assembly, comparative genomics, and studying sex chromosome evolution. The northern snakehead (Channa argus) is an economically important foodfish species with significant sexual dimorphism, where the males grow much faster and bigger than the females. However, to date, the sex determination pattern is still not clear, limiting identification of sex chromosomes, even sex determination genes and development of monosex populations that are valuable for both sex evolution of vertebrates and aquaculture practices. Here, a sex-averaged map and two sex-specific genetic maps were constructed with 2974, 2323, and 2338 SNPs, respectively. Little difference was observed in the pattern of sex-specific recombination between female- and male-specific genetic maps. Genome scan identified a major locus for sex determination at LG16. Females and males are, respectively, homogametic and heterogametic, suggesting an XY sex determination system for this species. By resequencing genomes, InDels in the sex-associated QTL region were discovered and used for developing sex-specific PCR assays for fast sexing of snakehead. These high-density genetic maps provide useful resources for future genomic studies in snakehead and its related species. The PCR assays for sexing are of importance in developing all male populations for aquaculture.

Construction of a High-Density Linkage Map and QTL Fine Mapping for Growth- and Sex-Related Traits in Channel Catfish (Ictalurus punctatus)

A high-density genetic linkage map is of particular importance in the fine mapping for important economic traits and whole genome assembly in aquaculture species. The channel catfish (Ictalurus punctatus), a species native to North America, is one of the most important commercial freshwater fish in the world. Outside of the United States, China has become the major producer and consumer of channel catfish after experiencing rapid development in the past three decades. In this study, based on restriction site associated DNA sequencing (RAD-seq), a high-density genetic linkage map of channel catfish was constructed by using single nucleotide polymorphisms (SNPs) in a F₁ family composed of 156 offspring and their two parental individuals. A total of 4,768 SNPs were assigned to 29 linkage groups (LGs), and the length of the linkage map reached 2,480.25 centiMorgans (cM) with an average distance of 0.55 cM between loci. Based on this genetic linkage map, 223 genomic scaffolds were anchored to the 29 LGs of channel catfish, and a total length of 704.66 Mb was assembled. Quantitative trait locus (QTL) mapping and genome-wide association analysis identified 10 QTLs of sex-related and six QTLs of growth-related traits at LG17 and LG28, respectively. Candidate genes associated with sex dimorphism, including spata2, spata5, sf3, zbtb38, and fox, were identified within QTL intervals on the LG17. A sex-linked marker with simple sequence repeats (SSR) in zbtb38 gene of the LG17 was validated for practical verification of sex in the channel catfish. Thus, the LG17 was considered as a sex-related LG. Potential growth-related genes were also identified, including important regulators such as megf9, npffr1, and gas1. In a word, we constructed the high-density genetic linkage map and developed the sex-linked marker in channel catfish, which are important genetic resources for future marker-assisted selection (MAS) of this economically important teleost.

Marker-assisted selection of YY supermales from a genetically improved farmed tilapia-derived strain

Genetically improved farmed tilapia (GIFT) and GIFT-derived strains account for the majority of farmed tilapia worldwide. As male tilapias grow much faster than females, they are often considered more desirable in the aquacultural industry. Sex reversal of females to males using the male sex hormone 17-α-methyltestosterone (MT) is generally used to induce phenotypic males during large-scale production of all male fingerlings. However, the widespread use of large quantities of sex reversal hormone in hatcheries may pose a health risk to workers and ecological threats to surrounding environments. Breeding procedures to produce genetically all-male tilapia with limited or no use of sex hormones are therefore urgently needed. In this study, by applying marker-assisted selection (MAS) for the selection of YY supermales from a GIFT-derived strain, we identified 24 XY pseudofemale and 431 YY supermale tilapias. Further performance evaluation on the progenies of the YY supermales resulted in male rates of 94.1%, 99.5% and 99.6%, respectively, in three populations, and a daily increase in body weight of 1.4 g at 3 months (n=997). Our study established a highly effective MAS procedure in the selection of YY supermales from a GIFT-derived strain. Furthermore, the development of MAS-selected YY supermales will help reduce the utilization of hormones for controlling sex in the tilapia aquaculture.

A novel PDZ domain-containing gene is essential for male sex differentiation and maintenance in yellow catfish (Pelteobagrus fulvidraco)

The sex-determining genes are found to be variable among different fish species. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish species in China with XX/XY sex-determining type. Recently, YY super-male yellow catfish has been successfully produced by combining hormonal-induced sex reversal method with sex chromosome-linked markers. Here, we identified a novel PDZ domain-containing gene in yellow catfish designated as pfpdz1, in whose intron the sex-linked marker was located. The coding sequence of pfpdz1 in Y chromosome was identical to that in X chromosome except a missense SNP (A/T) that changes an amino acid (E8V) in the N-terminal region. Pfpdz1 displayed male-specific expression during sex differentiation. Overexpression of pfpdz1 using additive transgenesis induces XX ovary to differentiate into testis-like tissue, while the targeted inactivation of pfpdz1 in Y chromosome using CRISPR/Cas9-mediated mutagenesis triggers ovarian differentiation. Furthermore, we demonstrated that pfpdz1 initiates testicular differentiation through upregulating expression of amh, dmrt1 and sox9a1, as well as downregulating expression of cyp19a1, foxl2 and wnt4. Our data provide functional evidence that pfpdz1 is significant for male differentiation and maintenance in yellow catfish.

Unusual Diversity of Sex Chromosomes in African Cichlid Fishes

African cichlids display a remarkable assortment of jaw morphologies, pigmentation patterns, and mating behaviors. In addition to this previously documented diversity, recent studies have documented a rich diversity of sex chromosomes within these fishes. Here we review the known sex-determination network within vertebrates, and the extraordinary number of sex chromosomes systems segregating in African cichlids. We also propose a model for understanding the unusual number of sex chromosome systems within this clade.

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

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.

An NGS-based approach for the identification of sex-specific markers in snakehead (Channa argus)

We described a next generation sequencing (NGS)-based approach to identify sex-specific markers and subsequently determine whether a species has male or female heterogamety. To test the accuracy of this technique, we examined the snakehead (Channa argus), which is economically important freshwater fish in China. Males grow faster than females, and there is significant interest in developing methods to skew breeding towards all-males to increase biomass yields. NGS was conducted on DNAs of individual female and male, the male reads were spitted into 60 bp K-mers and aligned to the female reference genome assembled by female reads, unaligned male K-mers-60 were kept in next filter process. Meanwhile, DNA sample of 48 females was pooled and sequenced, this data was further used to filter out the previous unaligned male K-mers-60. Hence, numbers of candidate Y chromosome-specific sequences were screened out, their sex-specificity were validated in wild snakeheads through PCR amplification. Finally, three Y chromosome-specific fragments (Contig-275834, Contig-359642, and Contig-418354) were identified, and specific primers were obtained to distinguish the sex of snakehead. Additionally, a pair of primers of Contig-275834 (275834X/Y-F and 275834X/Y-R) was exploited to distinguish XX females, XY males, and YY super-males, whose amplification products of different lengths were produced for different sexes. Therefore, our work demonstrated the ability of NGS data in identification of sex-specific markers, and the pipeline adopted in our study could be applied in any species of sex differentiation. Furthermore, the sex-specific markers have tremendous potential for improving the efficiency of all-male breeding practices in snakehead.

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

BACKGROUND

In Nile tilapia sex determination is governed by a male heterogametic system XX/XY either on LG1 or LG23. The latter carries a Y-specific duplicate of the amh gene, which is a testis-determining factor. Allelic variants in the amh gene demonstrated to be major triggers for autosomal and temperature-dependent sex reversal. Further, QTL on LG23 and LG20 show a temperature-responsiveness with influence on the phenotypic sex relative to the sex chromosomes. Here we present a ddRADseq based approach to identify genomic regions that show unusual large differentiation in terms of fixation index (F_ST) between temperature-treated pseudomales and non-masculinized females using a comparative genome-scan. Genome-wide associations were identified for the temperature-dependent sex using a genetically all-female population devoid of amh-ΔY.

RESULTS

Twenty-two thousand three hundred ninety-two SNPs were interrogated for the comparison of temperature-treated pseudomales and females, which revealed the largest differentiation on LG23. Outlier F_ST-values (0.35-0.44) were determined for six SNPs in the genomic interval (9,190,077-11,065,693) harbouring the amh gene (9,602,693-9,605,808), exceeding the genome-wide low F_ST of 0.013. Association analysis with a set of 9104 selected SNPs confirmed that the same genomic region on LG23 exerts a significant effect on the temperature-dependent sex.

CONCLUSIONS

This study highlights the role of LG23 in sex determination, harbouring major determinants for temperature-dependent sex reversal in Nile tilapia. Furthermore F_ST outlier detection proves a powerful tool for detection of sex-determining regions in fish genomes.

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 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.

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.

Electronic supplementary material

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

Screening and identification of a microsatellite marker associated with sex in Wami tilapia, Oreochromis urolepis hornorum

In this study, primer pairs of 15 microsatellite markers associated with sex determination of tilapia were selected and amplified in Wami tilapia, Oreochromis urolepis hornorum. While one marker, UNH168, on linkage group 3 (LG3) was associated (P <0.001) with the phenotypic sex in the experimental population, nine genotypes were detected in both sexes. Only 99-bp allele was detected in the female samples, while 141, 149 and 157-bp alleles were present in both male and female samples. UNH168 was localized by fluorescence in situ hybridization (FISH) on the long arm of the largest tilapia chromosome pair (chromosome 1, equivalent to LG3). This sex-linked microsatellite marker could potentially be used for marker-assisted selection in tilapia breeding programmes to produce monosex male tilapia.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

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.

A Tandem Duplicate of Anti-Müllerian Hormone with a Missense SNP on the Y Chromosome Is Essential for Male Sex Determination in Nile Tilapia, Oreochromis niloticus

Variation in the TGF-β signaling pathway is emerging as an important mechanism by which gonadal sex determination is controlled in teleosts. Here we show that amhy, a Y-specific duplicate of the anti-Müllerian hormone (amh) gene, induces male sex determination in Nile tilapia. amhy is a tandem duplicate located immediately downstream of amhΔ-y on the Y chromosome. The coding sequence of amhy was identical to the X-linked amh (amh) except a missense SNP (C/T) which changes an amino acid (Ser/Leu92) in the N-terminal region. amhy lacks 5608 bp of promoter sequence that is found in the X-linked amh homolog. The amhΔ-y contains several insertions and deletions in the promoter region, and even a 5 bp insertion in exonVI that results in a premature stop codon and thus a truncated protein product lacking the TGF-β binding domain. Both amhy and amhΔ-y expression is restricted to XY gonads from 5 days after hatching (dah) onwards. CRISPR/Cas9 knockout of amhy in XY fish resulted in male to female sex reversal, while mutation of amhΔ-y alone could not. In contrast, overexpression of Amhy in XX fish, using a fosmid transgene that carries the amhy/amhΔ-y haplotype or a vector containing amhy ORF under the control of CMV promoter, resulted in female to male sex reversal, while overexpression of AmhΔ-y alone in XX fish could not. Knockout of the anti-Müllerian hormone receptor type II (amhrII) in XY fish also resulted in 100% complete male to female sex reversal. Taken together, these results strongly suggest that the duplicated amhy with a missense SNP is the candidate sex determining gene and amhy/amhrII signal is essential for male sex determination in Nile tilapia. These findings highlight the conserved roles of TGF-β signaling pathway in fish sex determination.

Unlike mammals, the identity of the master sex-determining gene varies among fish species, and it is not yet clear if there is a common molecular pathway regulating gonadal sex determination across teleosts. Here we show that a Y-linked duplicate of the anti-Mullerian hormone (amhy) is essential for male sex determination in tilapia. Mutation of amhy resulted in male to female sex reversal, while overexpression of it resulted in female to male sex reversal. A missense single nucleotide polymorphisms (SNP) (C/T) in the open reading frame (ORF) of amhy might contribute to male sex determination in tilapia. Knockout of the anti-Müllerian hormone receptor type II (amhrII) also resulted in male to female sex reversal. Taken the amhy in Patagonian pejerrey, amhrII in Takifugu rubripes, gsdf^Y in Oryzias luzonensis into consideration, these data highlight an important role for TGF-β signaling in teleost sex determination.

The role of Amh signaling in teleost fish--Multiple functions not restricted to the gonads

This review summarizes the important role of Anti-Müllerian hormone (Amh) during gonad development in fishes. This Tgfβ-domain bearing hormone was named after one of its known functions, the induction of the regression of Müllerian ducts in male mammalian embryos. Later in development it is involved in male and female gonad differentiation and extragonadal expression has been reported in mammals as well. Teleosts lack Müllerian ducts, but they have amh orthologous genes. amh expression is reported from 21 fish species and possible regulatory interactions with further factors like sex steroids and gonadotropic hormones are discussed. The gonadotropin Fsh inhibits amh expression in all fish species studied. Sex steroids show no consistent influence on amh expression. Amh is produced in male Sertoli cells and female granulosa cells and inhibits germ cell proliferation and differentiation as well as steroidogenesis in both sexes. Therefore, Amh might be a central player in gonad development and a target of gonadotropic Fsh. Furthermore, there is evidence that an Amh-type II receptor is involved in germ cell regulation. Amh and its corresponding type II receptor are also present in brain and pituitary, at least in some teleosts, indicating additional roles of Amh effects in the brain-pituitary-gonadal axis. Unraveling Amh signaling is important in stem cell research and for reproduction as well as for aquaculture and in environmental science.

A novel sex-determining QTL in Nile tilapia (Oreochromis niloticus)

BACKGROUND

Fish species often exhibit significant sexual dimorphism for commercially important traits. Accordingly, the control of phenotypic sex, and in particular the production of monosex cultures, is of particular interest to the aquaculture industry. Sex determination in the widely farmed Nile tilapia (Oreochromis niloticus) is complex, involving genomic regions on at least three chromosomes (chromosomes 1, 3 and 23) and interacting in certain cases with elevated early rearing temperature as well. Thus, sex ratios may vary substantially from 50%.

RESULTS

This study focused on mapping sex-determining quantitative trait loci (QTL) in families with skewed sex ratios. These included four families that showed an excess of males (male ratio varied between 64% and 93%) when reared at standard temperature (28°C) and a fifth family in which an excess of males (96%) was observed when fry were reared at 36°C for ten days from first feeding. All the samples used in the current study were genotyped for two single-nucleotide polymorphisms (rs397507167 and rs397507165) located in the expected major sex-determining region in linkage group 1 (LG 1). The only misassigned individuals were phenotypic males with the expected female genotype, suggesting that those offspring had undergone sex-reversal with respect to the major sex-determining locus. We mapped SNPs identified from double digest Restriction-site Associated DNA (ddRAD) sequencing in these five families. Three genetic maps were constructed consisting of 641, 175 and 1,155 SNPs from the three largest families. QTL analyses provided evidence for a novel genome-wide significant QTL in LG 20. Evidence was also found for another sex-determining QTL in the fifth family, in the proximal region of LG 1.

CONCLUSIONS

Overall, the results from this study suggest that these previously undetected QTLs are involved in sex determination in the Nile tilapia, causing sex reversal (masculinisation) with respect to the XX genotype at the major sex-determining locus in LG 1.

A microsatellite-based genetic linkage map and putative sex-determining genomic regions in Lake Victoria cichlids

Cichlid fishes in East Africa have undergone extensive adaptive radiation, which has led to spectacular diversity in their morphology and ecology. To date, genetic linkage maps have been constructed for several tilapias (riverine), Astatotilapia burtoni (Lake Tanganyika), and hybrid lines of Lake Malawi cichlids to facilitate genome-wide comparative analyses. In the present study, we constructed a genetic linkage map of the hybrid line of Lake Victoria cichlids, so that maps of cichlids from all the major areas of East Africa will be available. The genetic linkage map shown here is derived from the F2 progeny of an interspecific cross between Haplochromis chilotes and Haplochromis sauvagei and is based on 184 microsatellite and two single-nucleotide polymorphism (SNP) markers. Most of the microsatellite markers used in the present study were originally designed for other genetic linkage maps, allowing us to directly compare each linkage group (LG) among different cichlid groups. We found 25 LGs, the total length of which was 1133.2cM with an average marker spacing of about 6.09cM. Our subsequent linkage mapping analysis identified two putative sex-determining loci in cichlids. Interestingly, one of these two loci is located on cichlid LG5, on which the female heterogametic ZW locus and several quantitative trait loci (QTLs) related to adaptive evolution have been reported in Lake Malawi cichlids. We also found that V1R1 and V1R2, candidate genes for the fish pheromone receptor, are located very close to the recently detected sex-determining locus on cichlid LG5. The genetic linkage map study presented here may provide a valuable foundation for studying the chromosomal evolution of East African cichlids and the possible role of sex chromosomes in generating their genomic diversity.

Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish

The traits of cultured fish must continually be genetically improved to supply high-quality animal protein for human consumption. Economically important fish traits are controlled by multiple gene quantitative trait loci (QTL), most of which have minor effects, but a few genes may have major effects useful for molecular breeding. In this review, we chose relevant studies on some of the most intensively cultured fish and concisely summarize progress on identifying and verifying QTLs for such traits as growth, disease and stress resistance and sex in recent decades. The potential applications of these major-effect genes and their associated markers in marker-assisted selection and molecular breeding, as well as future research directions are also discussed. These genetic and genomic analyses will be valuable for elucidating the mechanisms modulating economically important traits and to establish more effective molecular breeding techniques in fish.

A syntenic region conserved from fish to Mammalian x chromosome

Sex chromosomes bearing the sex-determining gene initiate development along the male or female pathway, no matter which sex is determined by XY male or ZW female heterogamety. Sex chromosomes originate from ancient autosomes but evolved rapidly after the acquisition of sex-determining factors which are highly divergent between species. In the heterogametic male system (XY system), the X chromosome is relatively evolutionary silent and maintains most of its ancestral genes, in contrast to its Y counterpart that has evolved rapidly and degenerated. Sex in a teleost fish, the Nile tilapia (Oreochromis niloticus), is determined genetically via an XY system, in which an unpaired region is present in the largest chromosome pair. We defined the differences in DNA contents present in this chromosome with a two-color comparative genomic hybridization (CGH) and the random amplified polymorphic DNA (RAPD) approach in XY males. We further identified a syntenic segment within this region that is well conserved in several teleosts. Through comparative genome analysis, this syntenic segment was also shown to be present in mammalian X chromosomes, suggesting a common ancestral origin of vertebrate sex chromosomes.

Structure and decay of a proto-Y region in Tilapia, Oreochromis niloticus

Background

Sex-determination genes drive the evolution of adjacent chromosomal regions. Sexually antagonistic selection favors the accumulation of inversions that reduce recombination in regions adjacent to the sex-determination gene. Once established, the clonal inheritance of sex-linked inversions leads to the accumulation of deleterious alleles, repetitive elements and a gradual decay of sex-linked genes. This in turn creates selective pressures for the evolution of mechanisms that compensate for the unequal dosage of gene expression. Here we use whole genome sequencing to characterize the structure of a young sex chromosome and quantify sex-specific gene expression in the developing gonad.

Results

We found an 8.8 Mb block of strong differentiation between males and females that corresponds to the location of a previously mapped sex-determiner on linkage group 1 of Oreochromis niloticus. Putatively disruptive mutations are found in many of the genes within this region. We also found a significant female-bias in the expression of genes within the block of differentiation compared to those outside the block of differentiation. Eight candidate sex-determination genes were identified within this region.

Conclusions

This study demonstrates a block of differentiation on linkage group 1, suggestive of an 8.8 Mb inversion encompassing the sex-determining locus. The enrichment of female-biased gene expression inside the proposed inversion suggests incomplete dosage compensation. This study helps establish a model for studying the early-to-intermediate stages of sex chromosome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-975) contains supplementary material, which is available to authorized users.

Identification of male-specific amh duplication, sexually differentially expressed genes and microRNAs at early embryonic development of Nile tilapia (Oreochromis niloticus)

Background

The probable influence of genes and the environment on sex determination in Nile tilapia suggests that it should be regarded as a complex trait. Detection of sex determination genes in tilapia has both scientific and commercial importance. The main objective was to detect genes and microRNAs that were differentially expressed by gender in early embryonic development.

Results

Artificial fertilization of Oreochromis niloticus XX females with either sex-reversed ΔXX males or genetically-modified YY ‘supermales’ resulted in all-female and all-male embryos, respectively. RNA of pools of all-female and all-male embryos at 2, 5 and 9 dpf were used as template for a custom Agilent eArray hybridization and next generation sequencing. Fifty-nine genes differentially expressed between genders were identified by a false discovery rate of p < 0.05. The most overexpressed genes were amh and tspan8 in males, and cr/20β-hsd, gpa33, rtn4ipl and zp3 in females (p < 1 × 10⁻⁹). Validation of gene expression using qPCR in embryos and gonads indicated copy number variation in tspan8, gpa33, cr/20β-hsd and amh. Sequencing of amh identified a male-specific duplication of this gene, denoted amhy, differing from the sequence of amh by a 233 bp deletion on exonVII, hence lacking the capability to encode the protein motif that binds to the transforming growth factor beta receptor (TGF-β domain). amh and amhy segregated in the mapping family in full concordance with SD-linked marker on LG23 signifying the QTL for SD. We discovered 831 microRNAs in tilapia embryos of which nine had sexually dimorphic expression patterns by a false discovery rate of p < 0.05. An up-regulated microRNA in males, pma-mir-4585, was characterized with all six predicted target genes including cr/20β-hsd, down-regulated in males.

Conclusions

This study reports the first discovery of sexually differentially expressed genes and microRNAs at a very early stage of tilapia embryonic development, i.e. from 2 dpf. Genes with sexually differential expression patterns are enriched for copy number variation. A novel male-specific duplication of amh, denoted amhy, lacking the TGF-β domain was identified and mapped to the QTL region on LG23 for SD, thus indicating its potential role in SD.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-774) contains supplementary material, which is available to authorized users.

Allelic variant in the anti-Müllerian hormone gene leads to autosomal and temperature-dependent sex reversal in a selected Nile tilapia line

Owing to the demand for sustainable sex-control protocols in aquaculture, research in tilapia sex determination is gaining momentum. The mutual influence of environmental and genetic factors hampers disentangling the complex sex determination mechanism in Nile tilapia (Oreochromis niloticus). Previous linkage analyses have demonstrated quantitative trait loci for the phenotypic sex on linkage groups 1, 3, and 23. Quantitative trait loci for temperature-dependent sex reversal similarly reside on linkage group 23. The anti-Müllerian hormone gene (amh), located in this genomic region, is important for sexual fate in higher vertebrates, and shows sexually dimorphic expression in Nile tilapia. Therefore this study aimed at detecting allelic variants and marker-sex associations in the amh gene. Sequencing identified six allelic variants. A significant effect on the phenotypic sex for SNP ss831884014 (p<0.0017) was found by stepwise logistic regression. The remaining variants were not significantly associated. Functional annotation of SNP ss831884014 revealed a non-synonymous amino acid substitution in the amh protein. Consequently, a fluorescence resonance energy transfer (FRET) based genotyping assay was developed and validated with a representative sample of fish. A logistic linear model confirmed a highly significant effect of the treatment and genotype on the phenotypic sex, but not for the interaction term (treatment: p<0.0001; genotype: p<0.0025). An additive genetic model proved a linear allele substitution effect of 12% in individuals from controls and groups treated at high temperature, respectively. Moreover, the effect of the genotype on the male proportion was significantly higher in groups treated at high temperature, giving 31% more males on average of the three genotypes. In addition, the groups treated at high temperature showed a positive dominance deviation (+11.4% males). In summary, marker-assisted selection for amh variant ss831884014 seems to be highly beneficial to increase the male proportion in Nile tilapia, especially when applying temperature-induced sex reversal.

Construction of a microsatellites-based linkage map for the white grouper (Epinephelus aeneus)

The white grouper (Epinephelus aeneus) is a promising candidate for domestication and aquaculture due to its fast growth, excellent taste, and high market price. A linkage map is an essential framework for mapping quantitative trait loci for economic traits and the study of genome evolution. DNA of a single individual was deep-sequenced, and microsatellite markers were identified in 177 of the largest scaffolds of the sequence assembly. The success rate of developing polymorphic homologous markers was 94.9% compared with 63.1% of heterologous markers from other grouper species. Of the 12 adult mature fish present in the broodstock tank, two males and two females were identified as parents of the assigned offspring by parenthood analysis using 34 heterologous markers. A single full-sib family of 48 individuals was established for the construction of first-generation linkage maps based on genotyping data of 222 microsatellites. The markers were assigned to 24 linkage groups in accordance to the 24 chromosomal pairs. The female and male maps consisting of 203 and 202 markers spanned 1053 and 886 cM, with an average intermarker distance of 5.8 and 5.0 cM, respectively. Mapping of markers to linkage groups ends was enriched by using markers originating from scaffolds harboring telomeric repeat-containing RNA. Comparative mapping showed high synteny relationships among the white grouper, kelp grouper (E. bruneus), orange-spotted grouper (E. coioides), and Nile tilapia (Oreochromis niloticus). Thus, it would be useful to integrate the markers that were developed for different groupers, depending on sharing of sequence data, into a comprehensive consensus map.

Bulk segregant RNA-seq reveals expression and positional candidate genes and allele-specific expression for disease resistance against enteric septicemia of catfish

Background

The application of RNA-seq has accelerated gene expression profiling and identification of gene-associated SNPs in many species. However, the integrated studies of gene expression along with SNP mapping have been lacking. Coupling of RNA-seq with bulked segregant analysis (BSA) should allow correlation of expression patterns and associated SNPs with the phenotypes.

Results

In this study, we demonstrated the use of bulked segregant RNA-seq (BSR-Seq) for the analysis of differentially expressed genes and associated SNPs with disease resistance against enteric septicemia of catfish (ESC). A total of 1,255 differentially expressed genes were found between resistant and susceptible fish. In addition, 56,419 SNPs residing on 4,304 unique genes were identified as significant SNPs between susceptible and resistant fish. Detailed analysis of these significant SNPs allowed differentiation of significant SNPs caused by genetic segregation and those caused by allele-specific expression. Mapping of the significant SNPs, along with analysis of differentially expressed genes, allowed identification of candidate genes underlining disease resistance against ESC disease.

Conclusions

This study demonstrated the use of BSR-Seq for the identification of genes involved in disease resistance against ESC through expression profiling and mapping of significantly associated SNPs. BSR-Seq is applicable to analysis of genes underlining various performance and production traits without significant investment in the development of large genotyping platforms such as SNP arrays.

Genetic differentiation and efficient sex-specific marker development of a pair of Y- and X-linked markers in yellow catfish

Pf62-Y and Pf62-X is a pair of allelic Y chromosome-linked and X chromosome-linked markers, and have been used to identify YY super-males, XY males and XX females for commercial production of all-male populations in yellow catfish (Pelteobagrus fulvidraco). However, the SCAR primers used previously have only two nucleotide difference, which restricts the wide utility because of nucleotide polymorphism. In this study, a continuous 8102 bp Pf62-Y sequence and a 5362 bp Pf62-X sequence have been cloned by genome walking, and significant genetic differentiation has been revealed between the corresponding X and Y chromosome allele sequences. Moreover, three pairs of primers were designed to efficiently identify YY super-males, XY males and XX females in an artificial breeding population, and to distinguish XY males and XX females in various wild populations. Together, the three new sex-specific genetic markers develop a highly stable and efficient method for genetic sex identification and sex control application in sustainable aquaculture of all-male yellow catfish.

Mapping and validation of the major sex-determining region in Nile tilapia (Oreochromis niloticus L.) Using RAD sequencing

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.

Characterization of gonadal transcriptomes from Nile tilapia (Oreochromis niloticus) reveals differentially expressed genes

Four pairs of XX and XY gonads from Nile tilapia were sequenced at four developmental stages, 5, 30, 90, and 180 days after hatching (dah) using Illumina Hiseq(TM) technology. This produced 28 Gb sequences, which were mapped to 21,334 genes. Of these, 259 genes were found to be specifically expressed in XY gonads, and 69 were found to be specific to XX gonads. Totally, 187 XX- and 1,358 XY-enhanced genes were identified, and 2,978 genes were found to be co-expressed in XX and XY gonads. Almost all steroidogenic enzymes, including cyp19a1a, were up-regulated in XX gonads at 5 dah; but in XY gonads these enzymes, including cyp11b2, were significantly up-regulated at 90 dah, indicating that, at a time critical to sex determination, the XX fish produced estrogen and the XY fish did not produce androgens. The most pronounced expression of steroidogenic enzyme genes was observed at 30 and 90 dah for XX and XY gonads, corresponding to the initiation of germ cell meiosis in the female and male gonads, respectively. Both estrogen and androgen receptors were found to be expressed in XX gonads, but only estrogen receptors were expressed in XY gonads at 5 dah. This could explain why exogenous steroid treatment induced XX and XY sex reversal. The XX-enhanced expression of cyp19a1a and cyp19a1b at all stages suggests an important role for estrogen in female sex determination and maintenance of phenotypic sex. This work is the largest collection of gonadal transcriptome data in tilapia and lays the foundation for future studies into the molecular mechanisms of sex determination and maintenance of phenotypic sex in non-model teleosts.

Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish

BACKGROUND

The availability of a large number of recently sequenced vertebrate genomes opens new avenues to integrate cytogenetics and genomics in comparative and evolutionary studies. Cytogenetic mapping can offer alternative means to identify conserved synteny shared by distinct genomes and also to define genome regions that are still not fine characterized even after wide-ranging nucleotide sequence efforts. An efficient way to perform comparative cytogenetic mapping is based on BAC clones mapping by fluorescence in situ hybridization. In this report, to address the knowledge gap on the genome evolution in cichlid fishes, BAC clones of an Oreochromis niloticus library covering the linkage groups (LG) 1, 3, 5, and 7 were mapped onto the chromosomes of 9 African cichlid species. The cytogenetic mapping data were also integrated with BAC-end sequences information of O. niloticus and comparatively analyzed against the genome of other fish species and vertebrates.

RESULTS

The location of BACs from LG1, 3, 5, and 7 revealed a strong chromosomal conservation among the analyzed cichlid species genomes, which evidenced a synteny of the markers of each LG. Comparative in silico analysis also identified large genomic blocks that were conserved in distantly related fish groups and also in other vertebrates.

CONCLUSIONS

Although it has been suggested that fishes contain plastic genomes with high rates of chromosomal rearrangements and probably low rates of synteny conservation, our results evidence that large syntenic chromosome segments have been maintained conserved during evolution, at least for the considered markers. Additionally, our current cytogenetic mapping efforts integrated with genomic approaches conduct to a new perspective to address important questions involving chromosome evolution in fishes.

First evidence for family-specific QTL for temperature-dependent sex reversal in Nile tilapia (Oreochromis niloticus)

This study for the first time screens microsatellite markers for associations with the temperature-dependent sex of Oreochromis niloticus. Previous studies revealed markers on linkage groups (LG) 1, 3, and 23 to be linked to the phenotypic sex of Oreochromis spp. at normal rearing temperatures. Moreover, candidate genes for sex determination and differentiation have been mapped to these linkage groups. Here, 6 families of a temperature-treated genetically all-female (XX) F(1)-population were genotyped for 21 microsatellites on the 3 LGs. No population-wide QTL (quantitative trait loci) or marker trait associations could be detected. However, family-specific QTL were found on LG 1 flanked by UNH995 and UNH104, on LG 3 at the position of GM213, and on LG 23 next to GM283. Moreover, family-specific single marker associations for UNH995 and UNH104 on LG 1, GM213 on LG 3, as well as for UNH898 and GM283 on LG 23 were detected. Yet, marker trait associations could not explain the temperature-dependent sex of all fish in the respective families. The molecular cue for the temperature-dependent sex in Nile tilapia might partially coincide with allelic variants at major and minor genetic sex determining factors. Moreover, additional QTL contributing to variable liabilities towards temperature might persist on other LGs.

Microsatellites Cross-Species Amplification across Some African Cichlids

The transfer of the genomic resources developed in the Nile tilapia, Oreochromis niloticus, to other Tilapiines sensu lato and African cichlid would provide new possibilities to study this amazing group from genetics, ecology, evolution, aquaculture, and conservation point of view. We tested the cross-species amplification of 32 O. niloticus microsatellite markers in a panel of 15 species from 5 different African cichlid tribes: Oreochromines (Oreochromis, Sarotherodon), Boreotilapiines (Tilapia), Chromidotilapines, Hemichromines, and Haplochromines. Amplification was successfully observed for 29 markers (91%), with a frequency of polymorphic (P₉₅) loci per species around 70%. The mean number of alleles per locus and species was 3.2 but varied from 3.7 within Oreochromis species to 1.6 within the nontilapia species. The high level of cross-species amplification and polymorphism of the microsatellite markers tested in this study provides powerful tools for a wide range of molecular genetic studies within tilapia species as well as for other African cichlids.

Linkage and Physical Mapping of Sex Region on LG23 of Nile Tilapia (Oreochromis niloticus)

Evidence supports that sex determination (SD) in tilapia is controlled by major genetic factors that may interact with minor genetic as well as environmental factors, thus implying that SD should be analyzed as a quantitative trait. Quantitative trait loci (QTL) for SD in Oreochromis niloticus were previously detected on linkage groups (LG) 1 and 23. Twenty-one short single repeats (SSR) of >12 TGs and one single nucleotide polymorphism were identified using the unpublished tilapia genome sequence on LG23. All markers showed two segregating alleles in a mapping family that was obtained by a cross between O. niloticus male (XY) and sex-reversed female (ΔXY) yielding 29 females (XX) and 61 males (XY and YY). Interval mapping analysis mapped the QTL peak between SSR markers ARO172 and ARO177 with a maximum F value of 78.7 (P < 7.6 × 10(-14)). Twelve adjacent markers found in this region were homozygous in females and either homozygous for the alternative allele or heterozygous in males. This segment was defined as the sex region (SR). The SR encompasses 1.5 Mbp on a single tilapia scaffold (no. 101) harboring 51 annotated genes. Among 10 candidate genes for SD that were tested for gene expression, anti-Müllerian hormone (Amh), which is located in the center of the SR, showed the highest overexpression in male vs. female embryos at 3 to 7 days postfertilization.