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.     

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.     

Demonstration of viability and fertility and development of a molecular tool to identify YY supermales in a fish with both genotypic and environmental sex determination

The pejerrey possesses a genotypic sex determination system driven by the amhy gene and yet shows marked temperature‐dependent sex determination. Sex‐reversed XY females have been found in a naturally breeding population established in Lake Kasumigaura, Japan. These females could mate with normal XY males and generate YY “supermale” individuals that, if viable and fertile, would sire only genotypic male offspring. This study was conducted to verify the viability, gender, and fertility of YY pejerrey and to develop a molecular method for their identification. Production of YY fish was attempted by crossing a thermally sex‐reversed XY female and an XY male, and rearing the progeny until sexual maturation. To identify the presumable YY individuals, we first conducted a PCR analysis using amhy‐specific primers to screen only amhy‐positive (XY and YY) fish. This screening showed that 60.6% of the progeny was amhy‐positive, which suggested the presence of YY fish. We then conducted a second screening by qPCR in order to identify the individuals with two amhy copies in their genome. This screening revealed 13 individuals, all males, with values twice higher than the other 30 amhy‐positive fishes, suggesting they have a YY complement. This assumption as well as the viability, fertility, and “supermale” nature of these individuals was confirmed in progeny tests with XX females that yielded 100% amhy‐positive offspring. These results demonstrate that qPCR can obviate progeny test as a means to identify the genotypic sex and therefore may be useful for the survey of all three possible genotypes in wild populations.     

A General Model to Explain Repeated Turnovers of Sex Determination in the Salicaceae

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex-determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.     

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.     

Identification of scuticociliates (Pseudocohnilembus persalinus,P. longisetus,Uronema marinum and Miamiensis avidus) based on the cox1 sequence

Scuticociliatosis is characterized as highly histophagous, causing systemic tissue destruction and high mortality in cultured marine fish. Some of the scuticociliates have been implicated as the causative agents of scuticociliatosis. Here, we describe our study to differentially identify various species in complex animal-sourced samples, namely olive flounder Paralichthys olivaceus and black rockfish Sebastes schlegelii suffering from scuticociliatosis. The mitochondrial cytochrome c oxidase 1 (cox1) gene from the scuticociliates was amplified and sequenced. The divergence percentage of small subunit ribosomal DNA sequence between average scuticociliate species was found to be low (8.3%) but the genetic divergence of cox1 sequence reached 23.5%, suggesting that a hyper-variable region of the cox1 gene could be used as a diagnostic DNA barcoding region. Thus, we developed species-specific primers for use in multiplex PCR of complex (pooled) samples. The primers yielded species-specific fragments (of distinct size) that allowed for simple, rapid, and effective identification and differentiation of multiple species present in a single sample.     

Inferring the Genetic Basis of Sex Determination from the Genome of a Dioecious Nightshade

Dissecting the genetic mechanisms underlying dioecy (i.e., separate female and male individuals) is critical for understanding the evolution of this pervasive reproductive strategy. Nonetheless, the genetic basis of sex determination remains unclear in many cases, especially in systems where dioecy has arisen recently. Within the economically important plant genus Solanum (∼2,000 species), dioecy is thought to have evolved independently at least 4 times across roughly 20 species. Here, we generate the first genome sequence of a dioecious Solanum and use it to ascertain the genetic basis of sex determination in this species. We de novo assembled and annotated the genome of Solanum appendiculatum (assembly size: ∼750 Mb scaffold N50: 0.92 Mb; ∼35,000 genes), identified sex-specific sequences and their locations in the genome, and inferred that males in this species are the heterogametic sex. We also analyzed gene expression patterns in floral tissues of males and females, finding approximately 100 genes that are differentially expressed between the sexes. These analyses, together with observed patterns of gene-family evolution specific to S. appendiculatum, consistently implicate a suite of genes from the regulatory network controlling pectin degradation and modification in the expression of sex. Furthermore, the genome of a species with a relatively young sex-determination system provides the foundational resources for future studies on the independent evolution of dioecy in this clade.     

Short-Term Fidelity,Habitat Use and Vertical Movement Behavior of the Black Rockfish Sebastes schlegelii as Determined by Acoustic Telemetry

The recent miniaturization of acoustic tracking devices has allowed fishery managers and scientists to collect spatial and temporal data for sustainable fishery management. The spatial and temporal dimensions of fish behavior (movement and/or vertical migrations) are particularly relevant for rockfishes (Sebastes spp.) because most rockfish species are long-lived and have high site fidelity, increasing their vulnerability to overexploitation. In this study, we describe the short-term (with a tracking period of up to 46 d) spatial behavior, as determined by acoustic tracking, of the black rockfish Sebastes schlegelii, a species subject to overexploitation in the Yellow Sea of China. The average residence index (the ratio of detected days to the total period from release to the last detection) in the study area was 0.92 ± 0.13, and most of the tagged fish were detected by only one region of the acoustic receiver array, suggesting relatively high site fidelity to the study area. Acoustic tracking also suggested that this species is more frequently detected during the day than at night in our study area. However, the diel detection periodicity (24 h) was only evident for certain periods of the tracking time, as revealed by a continuous wavelet transform. The habitat selection index of tagged S. schlegelii suggested that S. schlegelii preferred natural reefs, mixed sand/artificial reef bottoms and mixed bottoms of boulder, cobble, gravel and artificial reefs. The preference of this species for the artificial reefs that were recently deployed in the study area suggests that artificial seascapes may be effective management tools to attract individuals. The vertical movement of tagged S. schlegelii was mostly characterized by bottom dwelling behavior, and there was high individual variability in the vertical migration pattern. Our results have important implications for S. schlegelii catchability, the implementation of marine protected areas, and the identification of key species habitats, and our study provides novel information for future studies on the sustainability of this important marine resource in eastern China.     

Gene duplication in the evolution of sexual dimorphism

Males and females share most of the same genes, so selection in one sex will typically produce a correlated response in the other sex. Yet, the sexes have evolved to differ in a multitude of behavioral, morphological, and physiological traits. How did this sexual dimorphism evolve despite the presence of a common underlying genome? We investigated the potential role of gene duplication in the evolution of sexual dimorphism. Because duplication events provide extra genetic material, the sexes each might use this redundancy to facilitate sex‐specific gene expression, permitting the evolution of dimorphism. We investigated this hypothesis at the genome‐wide level in Drosophila melanogaster, using the presence of sex‐biased expression as a proxy for the sex‐specific specialization of gene function. We expected that if sexually antagonistic selection is a potent force acting upon individual genes, duplication will result in paralog families whose members differ in sex‐biased expression. Gene members of the same duplicate family can have different expression patterns in males versus females. In particular, duplicate pairs containing a male‐biased gene are found more frequently than expected, in agreement with previous studies. Furthermore, when the singleton ortholog is unbiased, duplication appears to allow one of the paralog copies to acquire male‐biased expression. Conversely, female‐biased expression is not common among duplicates; fewer duplicate genes are expressed in the female‐soma and ovaries than in the male‐soma and testes. Expression divergence exists more in older than in younger duplicates pairs, but expression divergence does not correlate with protein sequence divergence. Finally, genomic proximity may have an effect on whether paralogs differ in sex‐biased expression. We conclude that the data are consistent with a role of gene duplication in fostering male‐biased, but not female‐biased, gene expression, thereby aiding the evolution of sexual dimorphism.     

Phenotypic/genotypic sex mismatches and temperature‐dependent sex determination in a wild population of an Old World atherinid,the cobaltcap silverside Hypoatherina tsurugae

Several New World atheriniforms have been recognized as temperature‐dependent sex determined (TSD) and yet possess a genotypic sex determinant (amhy) which is primarily functional at mid‐range temperatures. In contrast, little is known about the sex determination in Old World atheriniforms, even though such knowledge is crucial to understand the evolution of sex determination mechanisms in fishes and to model the effects of global warming and climate change on their populations. This study examined the effects of water temperature on sex determination of an Old World atheriniform, the cobaltcap silverside Hypoatherina tsurugae, in which we recently described an amhy homologue. We first assessed the occurrence of phenotypic/genotypic sex mismatches in wild specimens from Tokyo Bay for three years (2014–2016) and used otolith analysis to estimate their birth dates and approximate thermal history during the presumptive period of sex determination. Phenotypic sex ratios became progressively biased towards males (47.3%–78.2%) during the period and were associated with year‐to‐year increases in the frequency of XX‐males (7.3%–52.0%) and decreases in XY/YY‐females (14.5%–0%). The breeding season had similar length but was delayed by about 1 month per year between 2014 and 2016, causing larvae to experience higher temperatures during the period of sex determination from year to year. Larval rearing experiments confirmed increased likelihood of feminization and masculinization at low and high temperatures, respectively. The results suggest that cobaltcap silverside has TSD, or more specifically the coexistence of genotypic and environmental sex determinants, and that it affects sex ratios in wild populations.     

Occurrence and Molecular Identification of Microcotyle sebastis Isolated from Fish Farms of the Korean Rockfish,Sebastes schlegelii

Microcotyle sebastis is a gill monogenean ectoparasite that causes serious problems in the mariculture of the Korean rockfish, Sebastes schlegelii. In this study, we isolated the parasite from fish farms along the coasts of Tongyeong, South Korea in 2016, and characterized its infection, morphology and molecular phylogeny. The prevalence of M. sebastis infection during the study period ranged from 46.7% to 96.7%, and the mean intensity was 2.3 to 31.4 ind./fish, indicating that the fish was constantly exposed to parasitic infections throughout the year. Morphological observations under light and scanning electron microscopes of the M. sebastis isolates in this study showed the typical characteristics of the anterior prohaptor and posterior opisthaptor of monogenean parasites. In phylogenetic trees reconstructed using the nuclear 28S ribosomal RNA gene and the mitochondrial cytochrome c oxidase I gene (cox1), they consistently clustered together with their congeneric species, and showed the closest phylogenetic relationships to M. caudata and M. kasago in the cox1 tree.     

Germ line control of female sex determination in zebrafish

A major transition during development of the gonad is commitment from an undifferentiated “bi-potential” state to ovary or testis fate. In mammals, the oogonia of the developing ovary are known to be important for folliculogenesis. An additional role in promoting ovary fate or female sex determination has been suggested, however it remains unclear how the germ line might regulate this process. Here we show that the germ line is required for the ovary versus testis fate choice in zebrafish. When the germ line is absent, the gonad adopts testis fate. These germ line deficient testes have normal somatic structures indicating that the germ line influences fate determination of surrounding somatic tissues. In germ line deficient animals the expression of the ovary specific gene cyp19a1a fails to be maintained whereas the testis genes sox9a and amh remain expressed. Furthermore, we observed decreased levels of the ovary specific genes cyp19a1a and foxL2 in germ line deficient animals prior to morphological sex differentiation of the gonad. We propose that the germ line has a common role in female sex determination in fish and mammals. Additionally, we show that testis specification is sufficient for masculinization of the fish pointing to a direct role of hormone signaling from the gonad in directing sex differentiation of non-gonadal tissues.     

Genomic characterization of sex‐identification markers in Sebastes carnatus and Sebastes chrysomelas rockfishes

Fish have evolved a variety of sex‐determining (SD) systems including male heterogamy (XY), female heterogamy (ZW) and environmental SD. Little is known about SD mechanisms of Sebastes rockfishes, a highly speciose genus of importance to evolutionary and conservation biology. Here, we characterize the sex determination system in the sympatrically distributed sister species Sebastes chrysomelas and Sebastes carnatus. To identify sex‐specific genotypic markers, double digest restriction site – associated DNA sequencing (ddRAD‐seq) of genomic DNA from 40 sexed individuals of both species was performed. Loci were filtered for presence in all of the individuals of one sex, absence in the other sex and no heterozygosity. Of the 74 965 loci present in all males, 33 male‐specific loci met the criteria in at least one species and 17 in both. Conversely, no female‐specific loci were detected, together providing evidence of an XY sex determination system in both species. When aligned to a draft reference genome from Sebastes aleutianus, 26 sex‐specific loci were interspersed among 1168 loci that were identical between sexes. The nascent Y chromosome averaged 5% divergence from the X chromosome and mapped to reference Sebastes genome scaffolds totalling 6.9Mbp in length. These scaffolds aligned to a single chromosome in three model fish genomes. Read coverage differences were also detected between sex‐specific and autosomal loci. A PCR‐RFLP assay validated the bioinformatic results and correctly identified sex of five additional individuals of known sex. A sex‐determining gene in other teleosts gonadal soma‐derived factor (gsdf) was present in the model fish chromosomes that spanned our sex‐specific markers.