Isolation of female-specific AFLP markers and molecular identification of genetic sex in half-smooth tongue sole (Cynoglossus semilaevis)

Development of DNA Insertion-specific Markers Based on the Intergenic Region of Oplegnathus punctatus Cdkn1/srsf3 for Sex Identification

Spotted knifejaw (Oplegnathus punctatus) is a marine economic fish with high food and ecological value, and its growth process has obvious male and female sexual dimorphism, with males growing significantly faster than females. However, the current sex identification technology is not yet mature, which will limit the growth rate of O. punctatus aquaculture and the efficiency of separate sex breeding, so the development of efficient sex molecular markers is imperative. This study identified a 926 bp DNA insertion fragment in the cdkn1/srsf3 intergenic region of O. punctatus males through whole-genome scanning, comparative genomics, and structural variant analysis. A pair of primers was designed based on the insertion information of the Y chromosome intergenic region in male individuals. Agarose gel electrophoresis revealed the amplification of two DNA fragments, 1118 bp and 192 bp, in male O. punctatus individuals. The 926 bp fragment was identified as the insertion in the intergenic region of cdkn1/srsf3 in males, while only a single 192 bp DNA fragment was amplified in females. The biological sex of the individuals identified in this manner was consistent with their known phenotypic sex. In this study, we developed a method to detect DNA insertion variants in the intergenic region of O. punctatus. Additionally, we introduced a new DNA marker for the rapid identification of the sex of O. punctatus, which enhances detection efficiency. The text has important reference significance and application value in sex identification, all-male breeding, and lineage selection. It provides new insights into the regulation of variation in the intergenic region of cdkn1/srsf3 genes and the study of RNA shearing.

Multi-organ transcriptomic profiles and gene-regulation network underlying vibriosis resistance in tongue sole

Vibrio spp. are major pathogens responsible for mortality and disease in various marine aquaculture organisms. Effective disease control and genetic breeding strategies rely heavily on understanding host vibriosis resistance mechanisms. The Chinese tongue sole (Cynoglossus semilaevis) is economically vital but suffers from substantial mortalities due to vibriosis. Through continuous selective breeding, we have successfully obtained vibriosis-resistant families of this species. In this study, we conducted RNA-seq analysis on three organs, including liver, spleen and intestine from selected resistant and susceptible tongue soles. Additionally, we integrated these data with our previously published RNA-seq datasets of skin and gill, enabling the construction of organ-specific transcriptional profiles and a comprehensive gene co-expression network elucidating the differences in vibriosis resistance. Furthermore, we identified 12 modules with organ-specific functional implications. Overall, our findings provide a valuable resource for investigating the molecular basis of vibriosis resistance in fish, offering insights into target genes and pathways essential for molecular selection and genetic manipulation to enhance vibriosis resistance in fish breeding programs.

Unveiling Gene Expression Dynamics during Early Embryogenesis in Cynoglossus semilaevis: A Transcriptomic Perspective

Commencing with sperm-egg fusion, the early stages of metazoan development include the cleavage and formation of blastula and gastrula. These early embryonic events play a crucial role in ontogeny and are accompanied by a dramatic remodeling of the gene network, particularly encompassing the maternal-to-zygotic transition. Nonetheless, the gene expression dynamics governing early embryogenesis remain unclear in most metazoan lineages. We conducted transcriptomic profiling on two types of gametes (oocytes and sperms) and early embryos (ranging from the four-cell to the gastrula stage) of an economically valuable flatfish-the Chinese tongue sole Cynoglossus semilaevis (Pleuronectiformes: Cynoglossidae). Comparative transcriptome analysis revealed that large-scale zygotic genome activation (ZGA) occurs in the blastula stage, aligning with previous findings in zebrafish. Through the comparison of the most abundant transcripts identified in each sample and the functional analysis of co-expression modules, we unveiled distinct functional enrichments across different gametes/developmental stages: actin- and immune-related functions in sperms; mitosis, transcription inhibition, and mitochondrial function in oocytes and in pre-ZGA embryos (four- to 1000-cell stage); and organ development in post-ZGA embryos (blastula and gastrula). These results provide insights into the intricate transcriptional regulation of early embryonic development in Cynoglossidae fish and expand our knowledge of developmental constraints in vertebrates.

Production of sterile mono-sex triploid yellow drum (Nibea albiflora): genotypic females and sex-reversed phenotypic males with emphasis on utilization as surrogate broodstock

Production of sterile mono-sex fish is of great significance for sustainable aquaculture as well as germ cell transplantation. In this study, we aimed to produce mono-sex triploid yellow drum, including genotypic females (XXX female) and sex-reversed phenotypic males (XXX male). Firstly, the mono-female triploids were produced through cold-shock treatment on eggs fertilized with sperm from neo-males. Then, the mono-male triploids were produced by the sex reversal of mono-female triploids with oral administration of letrozole (LZ). We comparatively investigated the growth and gonadal development in the mono-sex triploids. The results showed that the triploids displayed similar growth performance to their diploids throughout their first year, but had impaired gonadosomatic index and gametogenesis. No mature gametes were produced in the triploids during their first spawning season. Meanwhile, we analyzed the process of gametogenesis in the both sex of triploids. Ultrastructure of gametogenesis showed that the germ cells arrested at abnormal metaphase 1 in females, while males had irregular meiotic divisions, variable-sized spermatid and degenerated cells. The expression levels of meiosis-related genes (i.e., sycp3 and rec8) confirmed the abnormal meiosis in the triploids. Furthermore, the gonadal development was also determined by the expression patterns of vasa, dmrt1 and cyp19a1a. Abnormal expression of vasa mRNA and protein were detected in triploids. High cyp19a1a expression levels suggested the sex steroid hormones production might be at least partially functional in triploid females. In addition, high dmrt1 expression levels confirmed the masculinization and testicular development of sex-reversed triploid males by LZ. Our findings provide an efficient protocol to produce sterile mono-sex triploid yellow drum and provide new insights into the mechanism of gonadal sterility of triploid fish.

Phosphoproteomics Reveal New Candidates in Abnormal Spermatogenesis of Pseudomales in Cynoglossus semilaevis

Phosphorylation is a post-translational modification that contributes to versatile protein functions in spermatogenesis, and the variations they generate usually results in abnormal spermatogenesis or sperm dysfunction. The sex-reversal phenomenon exists in Chinese tongue sole under certain conditions such that individuals with a ZW genotype can acquire a male phenotype and are thus called pseudomales. Pseudomale tongue sole can reach sexual maturity but produce only Z-type sperm, and the Z sperm carries paternal epigenetic information. Whether phosphorylation plays a role in the sperm abnormality of pseudomales is unknown. In this study, a phosphoproteomic analysis was performed to compare protein phosphorylation profiles between pseudomale and male testes. Altogether, we identified 14,253 phosphopeptides matching with 4843 proteins, with 1329 differentially phosphorylated peptides corresponding to 1045 differentially phosphorylated proteins (DPPs). Phosphorylation at 781 sites was upregulated and at 548 sites was downregulated. Four motifs were identified among differentially phosphorylated peptides, which were "SP", "SD", "RxxS", and "TP". Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that the cell cycle and DNA/RNA processing were significantly enriched with the genes encoding DPPs. To analyze DPP function in depth, a protein-protein interaction network was constructed, and Ran-binding protein 2 was found to play a central role in spermatogenesis by regulating several processes such as the cell cycle, eukaryotic translation, ubiquitination, and minichromosome maintenance. In kinase-associated network analyses, two "mitogen-activated protein kinase (Mapk)-centered" clusters were identified that may account for abnormal spermatogenesis in pseudomales. One cluster was centered on Mapk6, which predominantly regulated the cell cycle by interacting with several cyclin-dependent kinases, and the other was centered on the "testis-expressed kinase 1-like (Tesk1l)/Pim1l-Mapk4l- testis-expressed 14 (Tex14)" kinase cascade, which might contribute to spermatogenesis by regulating β-catenin. Taken together, these data suggested the new candidates involved in pseudomale sperm abnormalities and provided clues to discover the phosphorylated regulatory mechanism underlying tongue sole spermatogenesis.

A chromosome-level genome assembly of Hong Kong catfish (Clarias fuscus) uncovers a sex-determining region

BACKGROUND

Hong Kong catfish (Clarias fuscus) is an ecologically and economically important species that is widely distributed in freshwater regions of southern China. Hong Kong catfish has significant sexual growth dimorphism. The genome assembly of the Hong Kong catfish would facilitate study of the sex determination and evolution mechanism of the species.

RESULTS

The first high-quality chromosome-level genome of the Hong Kong catfish was constructed. The total genome was 933.4 Mb, with 416 contigs and a contig N50 length of 8.52 Mb. Using high-throughput chromosome conformation capture (Hi-C) data, the genome assembly was divided into 28 chromosomes with a scaffold N50 length of 36.68 Mb. A total of 23,345 protein-coding genes were predicted in the genome, and 94.28% of the genes were functionally annotated in public databases. Phylogenetic analysis indicated that C. fuscus and Clarias magur diverged approximately 63.7 million years ago. The comparative genome results showed that a total of 60 unique, 353 expanded and 851 contracted gene families were identified in Hong Kong catfish. A sex-linked quantitative trait locus identified in a previous study was located in a sex-determining region of 30.26 Mb (0.02 to 30.28 Mb) on chromosome 13 (Chr13), the predicted Y chromosome. This QTL region contained 785 genes, of which 18 were identified as sex-related genes.

CONCLUSIONS

This study is the first to report the chromosome-level genome assembly of Hong Kong catfish. The study provides an excellent genetic resource that will facilitate future studies of sex determination mechanisms and evolution in fish.

Gonadal transcriptomes reveal sex-biased expression genes associated with sex determination and differentiation in red-tail catfish (Hemibagrus wyckioides)

BACKGROUND

Red-tail catfish (Hemibagrus wyckioides) is an important commercially farmed catfish in southern China. Males of red-tail catfish grow faster than females, suggesting that all-male catfish will produce more significant economic benefits in aquaculture practice. However, little research has been reported on sex determination and gonadal development in red-tail catfish.

RESULTS

In this study, we performed the first transcriptomic analysis of male and female gonads at four developmental stages at 10, 18, 30, and 48 days post hatching (dph) using RNA-seq technology. A total of 23,588 genes were screened in 24 sequenced samples, of which 28, 213, 636, and 1381 differentially expressed genes (DEGs) were detected at four developmental stages, respectively. Seven candidate genes of sex determination and differentiation were further identified. Real-time quantitative PCR (RT-qPCR) further confirmed that anti-Mullerian hormone (amh), growth differentiation factor 6a (gdf6a), testis-specific gene antigen 10 (tsga10), and cytochrome P450 family 17 subfamily A (cyp17a) were highly expressed mainly in the male, while cytochrome P450 family 19 subfamily A polypeptide 1b (cyp19a1b), forkhead box L2 (foxl2), and hydroxysteroid 17-beta dehydrogenase 1 (hsd17b1) were highly expressed in the female. The KEGG pathway enrichment data showed that these identified DEGs were mainly involved in steroid hormone biosynthesis and TGF-β signaling pathways.

CONCLUSIONS

Based on RNA-seq data of gonads at the early developmental stages, seven DEGs shared by the four developmental stages were identified, among which amh and gdf6a may be the male-biased expression genes, while foxl2, cyp19a1b and hsd17b1 may be the female-biased expression genes in red-tail catfish. Our study will provide crucial genetic information for the research on sex control in red-tail catfish, as well as for exploring the evolutionary processes of sex determination mechanisms in fish.

Screening and characterization of sex-specific sequences through 2b-RAD sequencing in American shad (Alosa sapidissima)

American shad (Alosa sapidissima), introduced from the United States, has become one of the most expensive farmed fish in the aquatic product market of China. The shad reveals significant sexual dimorphism in growth and behaviors. For the study, five male-specific tags were identified in two-generation breeding populations of Alosa sapidissima and were verified by PCR amplification. Averages of 10,245,091 and 8,685,704 raw and enzyme reads were obtained by high-throughput sequencing of the 2b-RAD library, respectively. 301,022 unique tags were obtained from the sequences of twenty samples with sequencing depths of 0 to 500. Finally, 274,324 special tags and 29,327 SNPs were selected with a sequencing depth of 3 to 500. Eleven preliminary screening male-specific tags and three male heterogametic SNP loci were isolated. After verification by PCR amplification, five male-specific sequences of 27 bp located on chromosome 3 were screened out. Chromosome 3 could be assumed to be the sex chromosome of Alosa sapidissima. Sex-specific markers will provide invaluable and systematic animal germplasm resources to allow for the precise identification of neo-males for the all-female breeding of Alosa sapidissima in commercial aquaculture.

Characterization and potential function of 7-dehydrocholesterol reductase (dhcr7) and lathosterol 5-desaturase (sc5d) in Cynoglossus semilaevis sexual size dimorphism

The typical sexual size dimorphism (SSD) phenomenon of Chinese tongue sole (Cynoglossus semilaevis) seriously restricts the sustainable development of the fishing industry. Previous transcriptome analysis has found a close relationship between the steroid biosynthesis and C. semilaevis SSD. The 7-dehydrocholesterol reductase (dhcr7) and lathosterol 5-desaturase (sc5d) are two genes in the steroid biosynthesis pathway, playing important roles in lipid synthesis, cellular metabolism, and growth. The present study assessed their roles in the mechanism of C. semilaevis SSD. The quantitative polymerase chain reaction (qPCR) results showed that C. semilaevis dhcr7 was mainly expressed in female livers, and C. semilaevis sc5d was highly expressed in female livers and gonads. Dual-luciferase experiment showed that dhcr7 and sc5d promoters had strong transcriptional activity. The transcription factors E2F transcription factor 1 (E2F1), and CCAAT enhancer binding protein alpha (C/EBPα) significantly regulated the transcriptional activity of dhcr7 and sc5d promoters, respectively. Furthermore, small interfering RNA (siRNA) knockdown results showed that expression levels of several genes [SREBF chaperone (scap), membrane-bound transcription factor peptidase, site 1 (mbtps1), fatty acid synthase (fasn), sonic hedgehog (shh), bone morphogenetic protein 2b (bmp2b) and AKT serine/threonine kinase 1 (akt1)] were suppressed. Protein subcellular localization results indicated that Dhcr7 and Sc5d were both specifically distributed in the cytoplasm, with co-localization been observed. The present study provides evidence that dhcr7 and sc5d might regulate C. semilaevis sexual size dimorphism by involving in energy homeostasis and cell cycle, or by affecting PI3K-Akt and Shh signaling pathways. The detailed roles of these steroid biosynthesis genes regulating C. semilaevis SSD needed more information.

Comparative transcriptome profiles of four sexually size dimorphic fish

Sexual size dimorphism is widespread in fish species. Although sex growth differences in multiple species have been studied successively, the commonalities of regulatory mechanisms across sexually dimorphic species are unknown. In this study, we performed RNA-seq analysis of four representative fish (loach, half-smooth tongue sole, yellow catfish, and Nile tilapia) with significant growth differences between females and males. Clean reads were identified from four fish species, ranging from 45,718,052 to 57,733,120. Following comparison transcriptome analysis, there were 1,132 and 1,108, 1,290 and 1,102, 4,732 and 4,266, 748 and 192 differentially expressed genes (DEGs) in the brain and muscle of loach, half-smooth tongue sole, yellow catfish, and Nile tilapia, respectively. Furthermore, the expression levels were validated by quantitative real-time PCR (qRT-PCR). Comparative transcriptome profiles of four fish described here will provide fundamental information for further studies on the commonalities of sexually size dimorphic fish in regulating growth differences between females and males.

Genomics and transcriptomics reveal new molecular mechanism of vibriosis resistance in fish

Infectious diseases have caused dramatic production decline and economic loss for fish aquaculture. However, the poor understanding of fish disease resistance severely hampered disease prevention. Chinese tongue sole (Cynoglossus semilaevis) is an important economic flatfish suffering from vibriosis. Here we used genomic, transcriptomic and experimental approaches to investigate the molecular genetic mechanisms underlying fish vibriosis resistance. A genome-wide comparison revealed that the genes under selective sweeps were enriched for glycosaminoglycan (GAG) chondroitin sulfate (CS)/dermatan sulfate (DS) metabolism. Transcriptomic analyses prioritized synergic gene expression patterns in this pathway, which may lead to an increased CS/DS content in the resistant family. Further experimental evidence showed that carbohydrate sulfotransferases 12 (Chst12), a key enzyme for CS/DS biosynthesis, has a direct antibacterial activity. To the best of our knowledge, this is the first report that the chst12 gene has a bactericidal effect. In addition, CS/DS is a major component of the extracellular matrix (ECM) and the selection signatures and fine-tuned gene expressions of ECM-receptor interaction genes indicated a modification in the ECM structure with an enhancement of the barrier function. Furthermore, functional studies conducted on Col6a2, encoding a collagen gene which constitutes the ECM, pointed to that it may act as a cellular receptor for Vibrio pathogens, thus plays an important role for the Vibrio invasion. Taken together, these findings provide new insights into the molecular protective mechanism underlying vibriosis resistance in fish, which offers crucial genomic resources for the resistant germplasm breeding and infectious disease control in fish culturing.

Intestinal microbiome-mediated resistance against vibriosis for Cynoglossus semilaevis

BACKGROUND

Infectious diseases have caused huge economic loss and food security issues in fish aquaculture. Current management and breeding strategies heavily rely on the knowledge of regulative mechanisms underlying disease resistance. Though the intestinal microbial community was linked with disease infection, there is little knowledge about the roles of intestinal microbes in fish disease resistance. Cynoglossus semilaevis is an economically important and widely cultivated flatfish species in China. However, it suffers from outbreaks of vibriosis, which results in huge mortalities and economic loss.

RESULTS

Here, we used C. semilaevis as a research model to investigate the host-microbiome interactions in regulating vibriosis resistance. The resistance to vibriosis was reflected in intestinal microbiome on both taxonomic and functional levels. Such differences also influenced the host gene expressions in the resistant family. Moreover, the intestinal microbiome might control the host immunological homeostasis and inflammation to enhance vibriosis resistance through the microbe-intestine-immunity axis. For example, Phaeobacter regulated its hdhA gene and host cyp27a1 gene up-expressed in bile acid biosynthesis pathways, but regulated its trxA gene and host akt gene down-expressed in proinflammatory cytokines biosynthesis pathways, to reduce inflammation and resist disease infection in the resistant family. Furthermore, the combination of intestinal microbes and host genes as biomarkers could accurately differentiate resistant family from susceptible family.

CONCLUSION

Our study uncovered the regulatory patterns of the microbe-intestine-immunity axis that may contribute to vibriosis resistance in C. semilaevis. These findings could facilitate the disease control and selective breeding of superior germplasm with high disease resistance in fish aquaculture. Video Abstract.

Screening of Genes Related to Sex Determination and Differentiation in Mandarin Fish (Siniperca chuatsi)

Mandarin fish has an XX/XY sex-determination system. The female mandarin fish is typically larger than the male. Sex identification and the discovery of genes related to sex determination in mandarin fish have important theoretical significance in the elucidation of the regulation and evolutionary mechanism of animal reproductive development. In this study, the chromosome-level genome of a female mandarin fish was assembled, and we found that LG24 of the genome was an X chromosome. A total of 61 genes on the X chromosome showed sex-biased expression. Only six gonadal genes (LG24G00426, LG24G003280, LG24G003300, LG24G003730, LG24G004200, and LG24G004770) were expressed in the testes, and the expression of the other gene LG24G003870 isoform 1 in the ovaries was significantly higher than that in the testes (p < 0.01). Five (except LG24G003280 and LG24G003300) of the seven aforementioned genes were expressed at the embryonic development stage, suggesting their involvement in early sex determination. The expression of LG24G004770 (encoding HS6ST 3-B-like) was also significantly higher in female muscles than in male muscles (p < 0.01), indicating other functions related to female growth. ZP3 encoded by LG24G003870 isoform 1 increased the C-terminal transmembrane domain, compared with that encoded by other fish zp3 isoforms, indicating their different functions in sex determination or differentiation. This study provides a foundation for the identification of sex-determining genes in mandarin fish.

De novo gonad transcriptome analysis of elongate loach (Leptobotia elongata) provides novel insights into sex-related genes

Elongate loach (Leptobotia elongata) is endemic in middle and upper reaches of the Yangtze River in China. Because of many anthropogenic factors such as overfishing and dam construction, the loach has become a highly endangered species. So far, the genomic resources which benefit for species conservation and utilization are still lacking in elongate loach. Therefore, the first gonad transcriptome analysis of the loach was conducted in this study, providing novel insights into sex-related genes. A total of 286,800,660 clean reads with a total length of 42.02 Gb were obtained. 18,975 differentially expressed genes (DEGs) were identified, where 12,976 DEGs, especially Sox9a, Sox9b, Igf2 and Fgfr2, were upregulated in the testis, and 5999 DEGs, especially Zp3, Eg2, Plk1, Ccnb1, Cdc20 and Mos, were upregulated in the ovary. Meanwhile, some testis-specific genes (i.e. Cald1, Atp1a, Muc2 and Ca2) and ovary-specific genes (i.e. Ca4, Tuba, Acp5, Ccna, Larp6 and Nop4) were identified and verified. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs, we found a series of enrichment pathways related to reproduction in elongate loach, such as the MAPK signaling pathway, oxytocin signaling pathway and oocyte meiosis pathway. Twelve DEGs were randomly selected to verify RNA-seq results by qPCR. In conclusion, this study provides a data source to study the molecular characteristics and regulatory mechanisms of sex-related genes in elongate loach, which has a potential to improve the resource protection and aquaculture production of the loach.

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.

The chromosome-level genome assembly of the Japanese yellowtail jack Seriola aureovittata provides insights into genome evolution and efficient oxygen transport

Fishes of the genus Seriola are widely farmed and highly valued in global aquaculture production. To further understand their economically important traits and help improve aquaculture product quality and sustainability, we performed a chromosome-level genome construction for Seriola aureovittata. Combining two technologies, PacBio and BGISEQ-500, we assembled 649.86 Mb S. aureovittata genome sequences with a contig N50 of 22.21 Mb, and 98% of BUSCO genes were detected in total. The initial assembly was then further scaffolded into 24 pseudochromosomes using Hi-C data, indicating the high quality of the genome. Genome evolution analysis showed that many genes related to fatty acid metabolism and oxygen binding, or transport were expanded, which provided insights into the metabolic characteristics of fatty acids and efficient oxygen transport. Based on the genome data, we confirmed the evolutionary relationship of S. aureovittata, S. dorsalis and S. lalandi and identified chr12 as the putative sex chromosome of S. aureovittata. Our chromosome-level genome assembly provides a genetic foundation for the phylogenetic and taxonomic investigation of different Seriola species. Moreover, the genome will provide an important genomic resource for further biological and aquaculture studies of S. aureovittata.

Sexual dimorphism in Odontobutis sinensis brain-pituitary-gonad axis and liver highlighted by histological and transcriptomic approach

In this study, sexual dimorphism in Chinese dark sleeper (Odontobutis sinensis) brain-pituitary-gonad axis and liver was highlighted by histological and transcriptomic approach. The results showed that there were two significant differences between males and females. Firstly, males grew larger and faster than females. Transcriptomic analysis and qPCR validation indicated that two key growth genes, insulin-like growth factor (igf) and 25-hydroxyvitamin D-1 alpha hydroxylase (cyp27b), were more highly detected in male liver than that in female liver. Secondly, histological analysis displayed that the liver in males showed an obvious ivory fatty phenomenon with more fat vacuoles and lipid droplet aggregation compared to that in females. Transcriptomic analysis indicated that the transcript level of vitellogenin (vtg) in male liver were significantly lower than that in female liver. After 17β-estradiol (E2) treatment of primary cultured hepatocytes, the vtg mRNA expression was induced significantly, while dihydrotestosterone (DHT) treatment had little effect on it. Generally, this study will provide some ideas for further exploring the mechanism of sexual dimorphism in Odontobutis sinensis.

Comparative piRNAs Profiles Give a Clue to Transgenerational Inheritance of Sex-Biased piRNAs in Cynoglossus semilaevis

Piwi interacting RNAs (piRNAs) are involved in the epigenetic and post-transcriptional gene silencing of retrotransposons in germ line cells, especially in spermatogenesis. There are many related reports on model organisms, such as flies and mice. In fish, however, there are few studies on piRNAs. Cynoglossus semilaevis, a benthic warm water flatfish, with remarkable sexual dimorphism, especially the "pseudo males" with sex reversal, mating with normal females to produce viable offspring, is an ideal material for the study of sex development. Here, sperm piwi-interacting RNAs profiles of Cynoglossus semilaevis were characterized, comparing between male and pseudomale groups. Differential piRNAs were identified with their predicted and annotated targets. Attention was then focused on candidate piRNAs associated with sex development and methylation. We continued to compare the expression levels of 10 candidates differentially expressed piRNAs in F1 spermatozoa. Quantitative RT-PCR demonstrated that five of the ten piRNAs showed sex bias consistent with parental sequencing results, with four significantly higher expression level in sperm of five males offspring than that of pseudomales, while one piRNAs showed the opposite expression profile. The five signature piRNAs (piR-mmu-49600337, piR-mmu-95849, piR-xtr-7474223, piR-xtr-1790334, and piR-mmu-4491546) could be employed as male-specific molecular biomarkers for C. semilaevis. Besides, this study also implied the possibility of transgenerational inheritance of sex-biased piRNAs exiting in sperm of Cynoglossus semilaevis.

Transcriptomic analysis revealed gene expression profiles during the sex differentiation of Chinese tongue sole (Cynoglossus semilaevis)

Sex differentiation in aquatic fish is important both for theoretical study and practical production, as growth dimorphism frequently appears in different sexes, especially in marine fish. The deciphered genome, identification of the male-determining gene dmrt1 and established genotypic sex screening method make Chinese tongue sole (Cynoglossus semilaevis) an ideal model to study sex differentiation in fish. In this study, comparative gonadal transcriptomic analyses were conducted for genetic females and males at 48, 68, and 108 days post hatching (dph), representing pre-, during- and post-gonadal differentiation stages, although the gonad is not completely differentiated and isolable in 48 and 68 dph individuals, while it is in 108 dph individuals. Altogether, 28 libraries were constructed, and a mean of 46.64 M clean reads was obtained. Differentially expressed gene (DEG) analysis revealed that 179 genes had similar expression patterns in males and females in all three stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that the enriched pathways included ubiquitin-mediated proteolysis, lysosomes, and RNA degradation. Moreover, weighted gene coexpression network analyses (WGCNA) identified 14 modules, one of which was closely correlated with female differentiation, exhibiting female-biased expression in all three stages (48, 68, 108 dph). An illustrated core gene interaction network of this module identified 50 genes, most of which are on W chromosomes. Six genes, including two ubiquitin conjugating enzymes, were selected for further investigation, and their female-biased expression was confirmed in even earlier stages, at 10 and 30 dph. These data facilitate our understanding of sex differentiation in fish and provide a genomic rationale for screening candidate genes (preferentially W-linked genes) that could be involved in the female differentiation process.

Characterization of Chinese tongue sole (Cynoglossus semilaevis) 24-dehydrocholesterol reductase: Expression profile, epigenetic modification, and its knock-down effect

The sexual size dimorphism of the Chinese tongue sole (Cynoglossus semilaevis) has greatly obstructed its sustainable development; however, the underlying mechanism remains unclear. Based on C. semilaevis transcriptomic information, 24-dehydrocholesterol reductase (dhcr24) was identified in steroid biosynthesis, showing female-liver-biased expression. Dhcr24 has been reported to participate in various processes, such as cholesterol synthesis, oxidative stress response, neuroprotection, and cell survival. The present study assessed its role in the sexual size dimorphism in fish. First, detailed expression pattern analysis showed that dhcr24 mRNAs were extensively expressed in tissues and the highest levels were found in the liver and gonads of females. Analysis of the dhcr24 promoter region demonstrated different DNA methylation statuses in female, male, and pseudomale gonads with higher epigenetic modification in males. The confirmation of transcription activity of the dhcr24 promoter and putative transcription factors (e.g., ER, AR, SREBP, and POU1F1a) provides the foundation for studying its regulatory mechanism. Finally, dhcr24-siRNA mediated knock-down assay using C. semilaevis liver cells showed that steroid biosynthesis related genes (e.g., ebp, dhcr7, and sc5d), core component of PI3K/Akt pathway (e.g., pi3k), and igf1r exhibited different expression patterns. Further investigation on the interplay between steroid hormones, dhcr24, PI3K/Akt, and IGF-1 systems will be valuable to better understand the mechanism underlying the sexual size dimorphism in C. semilaevis.

Female-specific genomic regions and molecular sex identification of the clearhead icefish (Protosalanx hyalocranius)

Background

The clearhead icefish, Protosalanx hyalocranius, is an economically important fishery species in China. Since 1980s, P. hyalocranius was widely introduced into lakes and reservoirs of northern China for aquaculture. However, the lack of a rapid and cost-effective sex identification method based on sex specific genetic markers has hindered study on sex determination mechanisms and breeding applications.

Results

Female-specific genomic regions were discovered by comparing whole genome re-sequencing data of both males and females. Two female-specific genomic regions larger than 50 bp were identified, and one (598 bp) contained a putative FOXI gene, which was paralogous to another FOXI gene with sex-associated SNPs. The two FOXI sequences displayed significant length difference with nine deletions of total length of 230 bp. This deletion-type structural variation could be easily and efficiently detected by traditional PCR and agarose gel electrophoresis with one 569 bp band for males and two bands (569 and 339 bp) for females, which were validated in 50 females and 40 males with known phenotypic sexes.

Conclusions

The results provided structural genomic evidence for the ZZ/ZW sex determination system in P. hyalocranius discovered in our previous study with association analysis of SNPs. Moreover, the female-specific markers and rapid and cost-effective PCR-based genetic sex identification method should have applications in further studies of sex determination mechanism for this species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07830-9.

Sex bias miRNAs in Cynoglossus semilaevis could play a role in transgenerational inheritance

Alterations of non-coding RNA profiling in spermatozoa are candidate mechanisms related to changes in paternal environment and progeny. Transgenerational inheritance of sex in pseudomales of Cynoglossus semilaevis, a fish with significant sex dimorphism, is a typical example of non-Mendelian inheritance. In the present study, miRNA profiles of spermatozoa were compared between male and pseudomale of C. semilaevis. Differential miRNAs in sperm from F0 and F1 generation also provides clues for revealing the possible role of non-coding RNA mediated transgenerational inheritance. Four sexual bias miRNAs, dre-miR-26a-5p, dre-miR-27b-3p, dre-miR-125b-5p,pol-199a-5p, were identified and verified in F0 and F1 generation of C. semilaevis. All of them were highly expressed in male sperm compared with pseudomale sperm. Function of target genes indicates that target genes of these differential RNAs are highly correlated with sex differentiation, gametogenesis and maintenance of secondary sexual characteristics. In a word, identification of epigenetic markers in gametes has great prospects in predicting susceptibility and properties in offsprings, and providing an indicator of parentalgenetic property.

A simple and rapid method for fish sex identification based on recombinase-aided amplification and its use in Cynoglossus semilaevis

Fish sex identification is a basic technique of great importance for both fish genetic studies and fisheries. Due to the sexual reversal phenomenon in many fish species, a simple and rapid molecular identification method for fish genetic sex is urgently needed to suit versatile detection scenarios, such as point-of-need applications. In this study, we took Cynoglossus semilaevis as an example, established a recombinase-aided amplification (RAA)-based method for sex identification, and combined the RAA-detection with two result visualization approaches with distinct features, capillary electrophoresis (CE) and lateral flow dipstick (LFD). Specific primers and probe were designed to specifically detect the sex chromosome W of C. semilaevis in order to distinguish the genetic sex between males, pseudo-males and females. To evaluate the performance of our methods, the genetic sex for twenty-eight males, sixty-eight pseudo-males and fifty-four females were examined with the RAA-based method and classical PCR-based genotyping method, demonstrating the consistent results of sex identification between both methods. The RAA-LFD method is operationally simple, rapid (~ 30 min) and holds great potential for point-of-need applications of fish sex identification, including fishery fields. The method presented here could be effective for identifying fish gender with the ZW karyotype.

Involvement of glycolysis activation in flatfish sexual size dimorphism: Insights from transcriptomic analyses of Platichthys stellatus and Cynoglossus semilaevis

The starry flounder (Platichthys stellatus), a flatfish cultured at the margins of the North Pacific, displays an obvious female-biased growth advantage, similar to many other fish species. To reveal the molecular mechanism underlying sexual size dimorphism, a comparative transcriptomic analysis of the somatotropic and reproductive axes was conducted. In total, 156, 67, 3434, and 378 differentially expressed genes (DEGs) between female and male samples were obtained in the brain, liver, gonad, and muscle tissues (q < 0.05). These DEGs were significantly enriched for various GO terms, including ion channel activity, protein binding, lipid transporter activity, and glycolytic process. The significantly enriched KEGG pathways included insulin secretion, axon guidance, and glycolysis/gluconeogenesis. In a detailed analysis of DEGs in these significantly enriched pathways, 35 genes showed higher expression levels in female muscle tissues than in male muscle tissues. A protein-protein interaction network further revealed specific interactions involving the glycolysis related-protein enolase (ENO), triosephosphate isomerase (TPI), Bisphosphoglycerate mutase (BPGM), fructose-bisphosphate aldolase (ALDO), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH). Interestingly, the role of glycolysis/gluconeogenesis was supported by an analysis of common DEGs between P. stellatus and Chinese tongue sole (Cynoglossus semilaevis). These results indicate that the activation of glycolysis in female muscle tissues contributes to flatfish sexual size dimorphism.

A First Insight into the Gonad Transcriptome of Hong Kong Catfish (Clarias fuscus)

Hong Kong catfish (Clarias fuscus) exhibit sexual dimorphism, particularly in body size. Due to the fast growth rate of males, the sexual size dimorphism of Hong Kong catfish has become an economically important trait. However, limited knowledge is known about the molecular mechanisms of sex determination and sex differentiation in this species. In this study, a first de novo transcriptome sequencing analysis of testes and ovaries was performed to identify sex-biased genes in Hong Kong catfish. The results showed that a total of 290,291 circular consensus sequences (CCSs) were obtained, from which 248,408 full-length non-chimeric (FLNC) reads were generated. After non-redundant analysis, a total of 37,305 unigenes were predicted, in which 34,342 unigenes were annotated with multiple public databases. Comparative transcriptomic analysis identified 5750 testis-biased differentially expressed genes (DEGs) and 6991 ovary-biased DEGs. The enrichment analysis showed that DEGs were classified into 783 Gene Ontology (GO) terms and 16 Kyoto Encyclopedia of Gene and Genome (KEGG) pathways. Many DEGs were involved with sex-related GO terms and KEGG pathways, such as oocyte maturation, androgen secretion, gonadal development and steroid biosynthesis pathways. In addition, the expression levels of 23 unigenes were confirmed to validate the transcriptomic data by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first investigation into the transcriptome of Hong Kong catfish testes and ovaries. This study provides an important molecular basis for the sex determination and sex control breeding of Hong Kong catfish.

The linkage of cell cycle and DNA replication with growth difference in female Chinese tongue sole (Cynoglossus semilaevis): Analysis from transcriptomic study and WGCNA

In addition to the typical sexual size dimorphism, considerable size differences within the female population of the Chinese tongue sole (Cynoglossus semilaevis) have become a further bottleneck of the improvement of sole aquaculture. To identify the internal mechanism, transcriptomic analysis and weighted gene co-expression network analysis (WGCNA) were employed simultaneously. Transcriptomic analyses of brain, pituitary gland, liver, gonad, and muscle tissues from two female groups with size differences identified 109, 698, 1325, 2299, and 2141 differentially expressed genes (DEGs), respectively. The results of these enrichment analyses suggest that the up-regulation of neuroactive ligand-receptor interaction, cell cycle, DNA replication, and MAPK signaling pathway in the group with larger females may be involved in the regulation of the observed growth differences. WGCNA of DEGs showed that cell cycle and DNA replication might be crucial pathways for accelerating cell growth in the groups with larger females. Finally, a series of hub genes including 6-phosphofructokinase type C (pfkp), ribosome biogenesis protein (wdr12), bleomycin hydrolase (blmh), and semaphorin-3A (sema3a) were recognized by the illustrated network map of modules. The linkage of cell cycle, DNA replication, and hub genes in the growth regulation of C. semilaevis provides further information for a better understanding of growth differences in fish.

Novel molecular marker for RAA-LFD visual detection of Cynoglossus semilaevis sex

Due to the significant sex dimorphism, Cynoglossus semilaevis has long been a species of research interest in the field of artificial sex manipulation. The existence of pseudo-males both in the natural habitat and aquaculture enterprises also is indicative of the importance for identification of the genetic sex in this species. In the present study, there was elucidation of a novel molecular marker for utilizing the recombinase aided amplification-lateral flow dipstick (RAA-LFD) visual system to identify the genetic sex of C. semilaevis. This 533 bp novel marker is a differential single copy fragment between the Z and W chromosome of C. semilaevis and exists only in the W chromosome. After primer designing and probe labeling, this marker has been utilized in a RAA isothermal amplification system. There were 49 C. semilaevis specimens evaluated for genetic sex identification using both PCR-agarose gel electrophoresis based InDel marker detection and the novel RAA-LFD system. The results from conducting evaluations with the two methods were consistent in all samples. Also, results from sensitivity analysis with use of the RAA-LFD system indicated the detection system was effective and reliable from 10⁸ copy number to 10¹. With use of the RAA reaction, there was only need to utilize a constant temperature of 37 ℃ for specific DNA amplification within 30 min. The combination use of RAA with LFD resulted in more efficient and convenient sex determination with there being a lesser technical threshold.

Application of SNP in Genetic Sex Identification and Effect of Estradiol on Gene Expression of Sex-Related Genes in Strongylocentrotus intermedius

Sea urchin (Strongylocentrotus intermedius) is an economically important mariculture species in Asia, and its gonads are the only edible part. The efficiency of genetic breeding in sea urchins is hampered due to the inability to distinguish gender by appearance. In this study, we first identified a sex-associated single nucleotide polymorphism (SNP) by combining type IIB endonuclease restriction site-associated DNA sequencing (2b-RAD-seq) and genome survey. Importantly, this SNP is located within spata4, a gene specifically expressed in male. Knocking down of spata4 by RNA interference (RNAi) in male individuals led to the downregulation of other conserved testis differentiation-related genes and germ cell marker genes. We also revealed that sex ratio in this validated culture population of S. intermedius is not 1:1. Moreover, after a 58-day feeding experiment with estradiol, the expression levels of several conserved genes that are related to testis differentiation, ovary differentiation, and estrogen metabolism were dynamically changed. Taken together, our results will contribute toward improving breeding efficiency, developing sex-controlled breeding, and providing a solid base for understanding sex determination mechanisms in sea urchins.

Cloning, expression and functional analysis of the desert hedgehog (dhh) gene in Chinese tongue sole (Cynoglossus semilaevis)

Desert hedgehog (dhh) is a gene that is crucial for spermatogenesis and Leydig cell differentiation, but little is known regarding its influence on gonadal differentiation and development in fish. To understand its function, we cloned and characterized the dhh gene from Cynoglossus semilaevis (csdhh). The full length csdhh cDNA was 2473 bp, including a 1386 bp open reading frame (ORF), a 475 bp 5'-UTR, and a 612 bp 3'-UTR, encoding a predicted protein of 461 amino acid residues. Phylogenetic analysis showed that the putative protein belongs to the hedgehog (HH) family, and contains typical HH-N and HH-C domains. Amino acid sequence analysis revealed that CsDhh shares many features with Dhh analogues in other teleost species. Real-time quantitative PCR showed that csdhh was detected in eight different tissues in male and female tongue sole. During early embryonic development, the relative expression of the csdhh was significantly higher in the neural stage than in other embryonic developmental stages (P < 0.05). csdhh was detected at 20 days after hatching (dah) and at the critical period of male gonadal differentiation (80-95 dah), the relative expression of the csdhh was significantly higher in the male gonads than the female gonads. In 5, 8, and 12 month old gonads, the relative expression of the csdhh was significantly higher in male and pseudo-male than in female fish. The in situ hybridization (ISH) results showed that the hybridization signal was strongly expressed in primary and secondary spermatocytes, spermatids, and sertoli cells of the 1-year-old fish testis, with only weak signal expression in the corresponding ovarian tissue. These results suggest that csdhh is highly conserved in evolution and plays an important role in spermatogenesis in males and pseudo-males.

Efficacy of Whole Cell Inactivated Vibrio harveyi Vaccine against Vibriosis in a Marine Red Hybrid Tilapia (Oreochromis niloticus × O. mossambicus) Model

Vibrio harveyi causes vibriosis in various commercial marine fish species. The infection leads to significant economic losses for aquaculture farms, and vaccination is an alternative approach for the prevention and control of fish diseases for aquaculture sustainability. This study describes the use of formalin-killed Vibrio harveyi (FKVh) strain Vh1 as a vaccine candidate to stimulate innate and adaptive immunities against vibriosis in a marine red hybrid tilapia model. Tilapia are fast growing; cheap; resistant to diseases; and tolerant to adverse environmental conditions of fresh water, brackish water, and marine water and because of these advantages, marine red hybrid tilapia is a suitable candidate as a model to study fish diseases and vaccinations against vibriosis. A total of 180 healthy red hybrid tilapias were gradually adapted to the marine environment before being divided into two groups, with 90 fish in each group and were kept in triplicate with 30 fish per tank. Group 1 was vaccinated intraperitoneally with 100 µL of FKVh on week 0, and a booster dose was similarly administered on week 2. Group 2 was similarly injected with PBS. Skin mucus, serum, and gut lavage were collected weekly for enzyme-linked immunosorbent assay (ELISA) and a lysozyme activity assay from a total of 30 fish of each group. On week 4, the remaining 60 fish of Groups 1 and 2 were challenged with 10⁸ cfu/fish of live Vibrio harveyi. The clinical signs were monitored while the survival rate was recorded for 48 h post-challenge. Vaccination with FKVh resulted in a significantly (p < 0.05) higher rate of survival (87%) compared to the control (20%). The IgM antibody titer and lysozyme activities of Group 1 were significantly (p < 0.05) higher than the unvaccinated Groups 2 in most weeks throughout the experiment. Therefore, the intraperitoneal exposure of marine red hybrid tilapia to killed V. harveyi enhanced the resistance and antibody response of the fish against vibriosis.

Isolation of a Male-Specific Molecular Marker and Development of a Genetic Sex Identification Technique in Spotted Knifejaw (Oplegnathus punctatus)

Sex-specific DNA markers are very helpful for identifying genetic sex and studying sex determination mechanisms in fish. To identify the sex-specific markers of spotted knifejaw (Oplegnathus punctatus), we performed a comparative analysis of the female and male genomes. In this study, an 18 bp insertion was identified in the male genome after verification by sequencing depth and PCR. An effective and rapid method based on PCR was then developed to identify the genetic sex. A male-female-shared primer pair and a male-specific primer were designed for PCR amplification to avoid false-negative phenomena. To examine the primers in practice, we utilized hundreds of spotted knifejaw fish from different groups to identify their genetic sex, and the results were consistent with their phenotypic sex. The male-specific DNA marker would be helpful for artificial breeding, Y chromosome assembly and further study of the sex determination mechanism. This study is the first to identify an effective sex-specific marker in spotted knifejaw.

A combination of genome-wide association study screening and SNaPshot for detecting sex-related SNPs and genes in Cynoglossus semilaevis

Chinese tongue sole (Cynoglossus semilaevis) males and females exhibit great differences in growth rate and appearance. The species is heterogametic (ZW/ZZ) and has sex-reversed "pseudomales" that are genetically female and physiologically male. In this study, we identified eight sex-specific single nucleotide polymorphism (SNP) markers for the sex identification of C. semilaevis by using a combination of genome-wide association study (GWAS) screening and SnaPshot validation. Candidate SNPs were screened using genotyping by sequencing to perform GWAS of the differential SNPs between the sexes of C. semilaevis. The SNP loci were amplified using a multiplex PCR system and detected via SNaPshot, which enables multiplexing of up to 30-40 SNPs in a single assay and ensures high accuracy of the results. The molecular markers detected in our study were used to successfully identify normal males and pseudomales from 45 caught and 40 cultured C. semilaevis specimens. Linkage disequilibrium analysis showed that the eight SNP loci were related to each other, with a strong linkage. Moreover, we investigated the expression of prdm6 mRNA containing a missense SNP and confirmed that the gene is differentially expressed in the gonads of the different sexes of C. semilaevis; the expression of prdm6 mRNA was significantly higher in the males than in the females and pseudomales. This means prdm6 may be related to sex differentiation in C. semilaevis.

A newly developed genetic sex marker and its application to understanding chemically induced feminisation in roach (Rutilus rutilus)

Oestrogenic wastewater treatment works (WwTW) effluents discharged into UK rivers have been shown to affect sexual development, including inducing intersex, in wild roach (Rutilus rutilus). This can result in a reduced breeding capability with potential population level impacts. In the absence of a sex probe for roach it has not been possible to confirm whether intersex fish in the wild arise from genetic males or females, or whether sex reversal occurs in the wild, as this condition can be induced experimentally in controlled exposures to WwTW effluents and a steroidal oestrogen. Using restriction site-associated DNA sequencing (RAD-seq), we identified a candidate for a genetic sex marker and validated this marker as a sex probe through PCR analyses of samples from wild roach populations from nonpolluted rivers. We also applied the sex marker to samples from roach exposed experimentally to oestrogen and oestrogenic effluents to confirm suspected phenotypic sex reversal from males to females in some treatments, and also that sex-reversed males are able to breed as females. We then show, unequivocally, that intersex in wild roach populations results from feminisation of males, but find no strong evidence for complete sex reversal in wild roach at river sites contaminated with oestrogens. The discovered marker has utility for studies in roach on chemical effects, wild stock assessments, and reducing the number of fish used where only one sex is required for experimentation. Furthermore, we show that the marker can be applied nondestructively using a fin clip or skin swab, with animal welfare benefits.

Identification and application of piwi-interacting RNAs from seminal plasma exosomes in Cynoglossus semilaevis

BACKGROUND

Piwi-interacting RNAs (piRNAs) have been linked to epigenetic and post-transcriptional gene silencing of retrotransposons in germ line cells, particularly in spermatogenesis. Exosomes are important mediators of vesicle transport, and the piRNAs in exosomes might play an important role in cell communication and signal pathway regulation. Moreover, exosomic piRNAs are promising biomarkers for disease diagnosis and physiological status indication. We used Cynoglossus semilaevis because of its commercial value and its sexual dimorphism, particularly the sex reversed "pseudomales" who have a female karyotype, produce sperm, and copulate with normal females to produce viable offspring.

RESULTS

To determine whether piRNAs from fish germ line cells have similar features, seminal plasma exosomes from half-smooth tongue sole, C. semilaevis, were identified, and their small RNAs were sequenced and analysed. We identified six signature piRNAs as biomarkers in exosomes of seminal plasma from males and pseudomale C. semilaevis. Bioinformatic analysis showed that all six signatures were sex-related, and four were DNA methylation-related and transposition-related piRNAs. Their expression profiles were verified using real-time quantitative PCR. The expression of the signature piRNAs was markedly higher in males than in pseudomales. The signature piRNAs could be exploited as male-specific biomarkers in this fish.

CONCLUSIONS

These signatures provide an effective tool to explore the regulatory mechanism of sex development in C. semilaevis and may provide guidance for future research on the function of piRNAs in the generative mechanism of sex reversed "pseudomales" in C. semilaevis.

Accurate Sex Identification of Ancient Elephant and Other Animal Remains Using Low-Coverage DNA Shotgun Sequencing Data

Sex identification of ancient animal biological remains can benefit our understanding of historical population structure, demography and social behavior. Traditional methods for sex identification (e.g., osteological and morphometric comparisons) may be ineffective when animal remains are not well preserved, when sex distinguishing characteristics have not yet developed, or where organisms do not exhibit sex-associated phenotypic dimorphisms. Here we adapt a method developed for human sex determination so that it can be used to identify the sex of ancient and modern animal taxa. The method identifies sex by calculating the ratio of DNA reads aligning to the X chromosome to DNA reads aligning to autosomes (termed the Rx ratio). We tested the accuracy of this method using low coverage genomes from 15 modern elephants (Loxodonta africana) for which sex was known. We then applied this method to ancient elephant ivory samples for which sex was unknown, and describe how this method can be further adapted to the genomes of other taxa. This method may be especially useful when only low-coverage genomic data are obtainable. Furthermore, because this method relies on only the X and not the Y chromosome, it can be used to determine the sex of organisms for which a reference genome was obtained from a female or for which only the X chromosome is reported. Such taxa include the domestic cat, sheep, goat, and horse; and non-domesticated animals such as the Sumatran orangutan, western lowland gorilla and meerkat.

Genome-Wide Association Mapping and Gene Expression Analyses Reveal Genetic Mechanisms of Disease Resistance Variations in Cynoglossus semilaevis

The sustainable development of aquaculture has been impeded by infectious diseases worldwide. However, the genomic architecture and the genetic basis underlying the disease resistance remain poorly understood, which severely hampers both the understanding of the evolution of fish disease resistance traits and the prevention of these diseases in the aquaculture community. Cynoglossus semilaevis is a representative and commercially-important flatfish species. Here we combined genome-wide association study and Fst and nucleotide diversity filtration to identify loci important for the disease resistance. Based on 1,016,774 single-nucleotide polymorphisms (SNPs) identified from 650 Gb genome resequencing data of 505 individuals, we detected 33 SNPs significantly associated with disease resistance and 79 candidate regions after filtration steps. Both the allele frequencies and genotype frequencies of the associated loci were significantly different between the resistant and susceptible fish, suggesting a role in the genetic basis of disease resistance. The SNP with strongest association with disease resistance was located in Chr 17, at 145 bp upstream of fblx19 gene, and overlapped with the major quantitative trait locus previously identified. Several genes, such as plekha7, nucb2, and fgfr2, were also identified to potentially play roles in the disease resistance. Furthermore, the expression of some associating genes were likely under epigenetic regulations between the bacterial resistant and susceptible families. These results provide insights into the mechanism that enable variation of disease resistance to bacterial pathogen infection. The identified polymorphisms and genes are valuable targets and molecular resources for disease resistance and other traits, and for advanced breeding practice for superior germplasm in fish aquaculture.

Improvement of the Pacific bluefin tuna (Thunnus orientalis) reference genome and development of male-specific DNA markers

The Pacific bluefin tuna, Thunnus orientalis, is a highly migratory species that is widely distributed in the North Pacific Ocean. Like other marine species, T. orientalis has no external sexual dimorphism; thus, identifying sex-specific variants from whole genome sequence data is a useful approach to develop an effective sex identification method. Here, we report an improved draft genome of T. orientalis and male-specific DNA markers. Combining PacBio long reads and Illumina short reads sufficiently improved genome assembly, with a 38-fold increase in scaffold contiguity (to 444 scaffolds) compared to the first published draft genome. Through analysing re-sequence data of 15 males and 16 females, 250 male-specific SNPs were identified from more than 30 million polymorphisms. All male-specific variants were male-heterozygous, suggesting that T. orientalis has a male heterogametic sex-determination system. The largest linkage disequilibrium block (3,174 bp on scaffold_064) contained 51 male-specific variants. PCR primers and a PCR-based sex identification assay were developed using these male-specific variants. The sex of 115 individuals (56 males and 59 females; sex was diagnosed by visual examination of the gonads) was identified with high accuracy using the assay. This easy, accurate, and practical technique facilitates the control of sex ratios in tuna farms. Furthermore, this method could be used to estimate the sex ratio and/or the sex-specific growth rate of natural populations.

Seminal Plasma Exosomes: Promising Biomarkers for Identification of Male and Pseudo-Males in Cynoglossus semilaevis

In mammals, small RNAs enclosed in exosomes have been identified as appropriate signatures for disease diagnosis. However, there is limited information on exosomes derived from seminal plasma, and few studies have reported analyzed the composition of exosomes and enclosed small RNAs in fish. The half-smooth tongue sole (Cynoglossus semilaevis) is an economically important fish for aquaculture, and it exhibits sexual dimorphism: the female gender show higher growth rates and larger body sizes than males. Standard karyotype analysis and artificial gynogenesis tests have revealed that this species uses heterogametic sex determination (ZW/ZZ), and so-called sex-reversed pseudo-males exist. In this study, we successfully identified exosomes in the seminal plasma of C. semilaevis; to the best of our knowledge, this is the first report of exosomes in fish seminal plasma. Analysis of the nucleotide composition showed that miRNAs were dominant in the exosomes, and the miRNAs were sequenced and compared to identify signature miRNAs as sexual biomarkers. Moreover, target genes of the signature miRNAs were predicted by sequence matching and annotation. Finally, four miRNAs (dre-miR-141-3P, dre-miR-10d-5p, ssa-miR-27b-3p, and ssa-miR-23a-3p) with significant differential expression in the males and pseudo-males were selected from the signature candidate miRNAs as markers for sex identification, and their expression profiles were verified using real-time quantitative PCR. Our findings could provide an effective detection method for sex differentiation in fish.

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

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

Satellitome landscape analysis of Megaleporinus macrocephalus (Teleostei, Anostomidae) reveals intense accumulation of satellite sequences on the heteromorphic sex chromosome

The accumulation of repetitive DNA sequences on the sex-limited W or Y chromosomes is a well-known process that is likely triggered by the suppression of recombination between the sex chromosomes, which leads to major differences in their sizes and genetic content. Here, we report an analysis conducted on the satellitome of Megaleporinus macrocephalus that focuses specifically on the satDNAs that have been shown to have higher abundances in females and are putatively located on the W chromosome in this species. We characterized 164 satellite families in M. macrocephalus, which is, by far, the most satellite-rich species discovered to date. Subsequently, we mapped 30 satellites, 22 of which were located on the W chromosome, and 14 were shown to exist only on the W chromosome. Finally, we report two simple, quick and reliable methods that can be used for sex identification in M. macrocephalus individuals using fin clips or scales, which could be applicable to future studies conducted in the field of aquaculture.

Expression analysis and characterization of zglp1 in the Chinese tongue sole (Cynoglossus semilaevis)

Zinc finger GATA like protein-1 (ZGLP1) is a nuclear zinc finger protein that regulates the interaction between somatic cells and germ cells during gonad developmental process in mammals. In this study, the zglp1 of Chinese tongue sole, Cynoglossus semilaevis (cysezglp1), was cloned and characterized for the first time in fish. Cysezglp1 had an open reading frame with five exons and was located to chromosome 9. The open reading frame of cysezglp1 consisted of 1692 nucleotides and encoded a 583 amino acid polypeptide. The predicted protein contained two zinc finger structures (Znf1 and Znf2), one of which was highly homologous to the GATA-type zinc finger domain. Multiple sequence alignment showed that Znf1 was conserved across different species while Znf2 was more divergent. Through quantitative Real-time PCR (qRT-PCR), we found that cysezglp1 was predominantly expressed in gonads, and the expression level of the ovary was significantly higher than that of the testis. We compared expression level in different embryonic stages and found that cysezglp1 mRNAs were mainly expressed in the fertilized egg to the cleavage stage, subsequently declining in the blastula stage. Cysezglp1 expression was not detected from the gastrulation stage onward. In the ovary, cysezglp1 expression was detected at 120 days after hatching and expression gradually increased with the maturation of the ovary. In situ hybridization showed that the cysezglp1 was mainly expressed in oocytes. Taken together, our results suggest that cysezglp1 may play an important role in the process of oogenesis in Chinese tongue sole.

Identification, Characterization, and Mapping of a Novel SNP Associated with Body Color Transparency in Juvenile Red Sea Bream (Pagrus major)

We previously reported a body color deformity in juvenile red sea bream, which shows transparency in the juvenile stage because of delayed chromatophore development compared with normal individuals, and this finding suggested a genetic cause based on parentage assessments. To conduct marker-assisted selection to eliminate broodstock inheriting the causative gene, developing DNA markers associated with the phenotype was needed. We first conducted SNP mining based on AFLP analysis using bulked-DNA from normal and transparent individuals. One SNP was identified from a transparent-specific AFLP fragment, which significantly associated with transparent individuals. Two alleles (A/G) were observed in this locus, and the genotype G/G was dominantly observed in the transparent groups (97.1%) collected from several production lots produced from different broodstock populations. A few normal individuals inherited the G/G genotype (5.0%), but the A/A and A/G genotypes were dominantly observed in the normal groups. The homologs region of the SNP was searched using a medaka genome database, and intron 12 of the Nell2a gene (located on chromosome 6 of the medaka genome) was highly matched. We also mapped the red sea bream Nell2a gene on the previously developed linkage maps, and this gene was mapped on a male linkage group, LG4-M. The newly found SNP was useful in eliminating broodstock possessing the causative gene of the body color transparency observed in juvenile stage of red sea bream.

Liver transcriptome analysis reveals extensive transcriptional plasticity during acclimation to low salinity in Cynoglossus semilaevis

BACKGROUND

Salinity is an important abiotic stress that influences the physiological and metabolic activity, reproduction, growth and development of marine fish. It has been suggested that half-smooth tongue sole (Cynoglossus semilaevis), a euryhaline fish species, uses a large amount of energy to maintain osmotic pressure balance when exposed to fluctuations in salinity. To delineate the molecular response of C. semilaevis to different levels of salinity, we performed RNA-seq analysis of the liver to identify the genes and molecular and biological processes involved in responding to salinity changes.

RESULTS

The present study yielded 330.4 million clean reads, of which 83.9% were successfully mapped to the reference genome of C. semilaevis. One hundred twenty-eight differentially expressed genes (DEGs), including 43 up-regulated genes and 85 down-regulated genes, were identified. These DEGs were highly represented in metabolic pathways, steroid biosynthesis, terpenoid backbone biosynthesis, butanoate metabolism, glycerolipid metabolism and the 2-oxocarboxylic acid metabolism pathway. In addition, genes involved in metabolism, osmoregulation and ion transport, signal transduction, immune response and stress response, and cytoskeleton remodeling were affected during acclimation to low salinity. Genes acat2, fdps, hmgcr, hmgcs1, mvk, pmvk, ebp, lss, dhcr7, and dhcr24 were up-regulated and abat, ddc, acy1 were down-regulated in metabolic pathways. Genes aqp10 and slc6a6 were down-regulated in osmoregulation and ion transport. Genes abat, fdps, hmgcs1, mvk, pmvk and dhcr7 were first reported to be associated with salinity adaptation in teleosts.

CONCLUSIONS

Our results revealed that metabolic pathways, especially lipid metabolism were important for salinity adaptation. The candidate genes identified from this study provide a basis for further studies to investigate the molecular mechanism of salinity adaptation and transcriptional plasticity in marine fish.

Mayfly and fish species identification and sex determination in bleak (Alburnus alburnus) by MALDI-TOF mass spectrometry

Besides food quality control of fish or cephalopods, the novel mass spectrometry (MS) approaches could be effective and beneficial methods for the investigation of biodiversity in ecological research. Our aims were to verify the applicability of MALDI-TOF MS in the rapid identification of closely related species, and to further develop it for sex determination in phenotypically similar fish focusing on the low mass range. For MALDI-TOF MS spectra analysis, ClinProTools software was applied, but our observed classification was also confirmed by Self Organizing Map. For verifying the wide applicability of the method, brains from invertebrate and vertebrate species were used in order to detect the species related markers from two mayflies and eight fish as well as sex-related markers within bleak. Seven Ephemera larvae and sixty-one fish species related markers were observed and nineteen sex-related markers were identified in bleak. Similar patterns were observed between the individuals within one species. In contrast, there were markedly diverse patterns between the different species and sexes visualized by SOMs. Two different Ephemera species and male or female fish were identified with 100% accuracy. The various fish species were classified into 8 species with a high level of accuracy (96.2%). Based on MS data, dendrogram was generated from different fish species by using ClinProTools software. This MS-based dendrogram shows relatively high correspondence with the phylogenetic relationships of both the studied species and orders. In summary, MALDI-TOF MS provides a cheap, reliable, sensitive and fast identification tool for researchers in the case of closely related species using mass spectra acquired in a low mass range to define specific molecular profiles. Moreover, we presented evidence for the first time for determination of sex within one fish species by using this method. We conclude that it is a powerful tool that can revolutionize ecological and environmental research.

Molecular characterization of Pod1 during sex development in Chinese tongue sole (Cynoglossus semilaevis)

Pod1 encodes a Class II bHLH transcription factor involved in the development of a number of tissues such as gonad, spleen, lungs and heart. However, to date, little is known about its function in teleosts. In this study, we cloned and characterized Pod1 gene from Cynoglossus semilaevis. This gene contains three exons and two introns, with the full-length cDNA of 918 nucleotides that encodes a 183 amino acid protein with a conserved bHLH domain. Realtime quantitative PCR revealed that Pod1 was predominantly expressed in the testes of C. semilaevis. In different stages of testes development, Pod1 expression was undetectable up to 120 days after hatching (dah), and then increased at 210 dah and 1 year after hatching (yah). Furthermore, in situ hybridization (ISH) analysis revealed that Pod1 was mainly localized in the germ cells of testes, but was not detected in ovarian cells; which suggested its possible functions in spermatogenesis of C. semilaevis. The methylation profile analysis of Pod1 genomic sequence in the gonads showed that the differences in their putative promoter regions of Pod1 among ovary, male and pseudo-male testes were not obvious. Thus, further research might be needed to evaluate whether Pod1 expression is regulated by epigenetic level.