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

Gender Impacted Gut Microbiota and Growth Performance in the Blotched Snakehead (Channa maculata)

The blotched snakehead Channa maculata is an important economical freshwater species in East Asia. However, there has been relatively little research conducted on the correlation between gender and gut microbes. In this study, 36 of 1000 blotched snakeheads were randomly selected for growth performance measurement and gut microbiota high-throughput sequencing. Results showed that microbial diversity, composition, and metabolic functions were altered by gender and growth performance except the microbial network. In our study, Proteobacteria were the most abundant phylum, with Fusobacteria showing enrichment in males and Bacteroidetes in females. Notably, phylum Deinococcus-Thermus was identified as a significant biomarker. The Cetobacterium was the most abundant genus-level taxon. Furthermore, gut microbes specializing in the production of gut-healthy substances, such as coenzymes and vitamins, were identified as biomarkers in the fast-growing group. Our investigation highlighted the impact of gender on the composition and abundance of gut microbial biomarkers in both males and females, thereby influencing differential growth performance through the modulation of specific metabolic functions.

Identification of differentially expressed genes and SNPs linked to harvest body weight of genetically improved rohu carp, Labeo rohita

In most of the aquaculture selection programs, harvest body weight has been a preferred performance trait for improvement. Molecular interplay of genes linked to higher body weight is not elucidated in major carp species. The genetically improved rohu carp with 18% average genetic gain per generation with respect to harvest body weight is a promising candidate for studying genes' underlying performance traits. In the present study, muscle transcriptome sequencing of two groups of individuals, with significant difference in breeding value, belonging to the tenth generation of rohu carp was performed using the Illumina HiSeq 2000 platform. A total of 178 million paired-end raw reads were generated to give rise to 173 million reads after quality control and trimming. The genome-guided transcriptome assembly and differential gene expression produced 11,86,119 transcripts and 451 upregulated and 181 downregulated differentially expressed genes (DEGs) between high-breeding value and low-breeding value (HB & LB) groups, respectively. Similarly, 39,158 high-quality coding SNPs were identified with the Ts/Tv ratio of 1.23. Out of a total of 17 qPCR-validated transcripts, eight were associated with cellular growth and proliferation and harbored 13 SNPs. The gene expression pattern was observed to be positively correlated with RNA-seq data for genes such as myogenic factor 6, titin isoform X11, IGF-1 like, acetyl-CoA, and thyroid receptor hormone beta. A total of 26 miRNA target interactions were also identified to be associated with significant DETs (p-value < 0.05). Genes such as Myo6, IGF-1-like, and acetyl-CoA linked to higher harvest body weight may serve as candidate genes in marker-assisted breeding and SNP array construction for genome-wide association studies and genomic selection.

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.

Validation of a male-specific DNA marker confirms XX/XY-type sex determination in several Hungarian strains of African catfish (Clarias gariepinus)

African catfish (Clarias gariepinus) is a promising food fish species with significant potential and growing mass of production in freshwater aquaculture. Male African catfish possess improved production characteristics over females, therefore the use of monosex populations could be advantageous for aquaculture production. However, our knowledge about the sex determination mechanism of this species is still limited and controversial. A previously isolated male-specific DNA marker (CgaY1) was validated using offspring groups from targeted crosses (n = 630) and it was found to predict the sex of 608 individuals correctly (96.43% accuracy). Using the proportion of recombinants, we estimated the average genetic distance between the potential sex determination locus and the sex-specific marker to be 3.57 cM. As an earlier study suggested that both XX/XY and ZZ/ZW systems coexist in this species, we tested the applicability of their putative 'moderately sex-linked loci' and found that no sex-specific amplification could be detected for any of them. In addition, temperature-induced masculinization suggested by others was also tested, but no such effect was detected in our stocks when the published parameters were used for heat treatment. Altogether, our results support an exclusive XX/XY sex determination system in our African catfish stock and indicate a good potential for the future use of this male-specific DNA marker in research and commercial production.

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.

Sexually Dimorphic Gene Expression in X and Y Sperms Instructs Sexual Dimorphism of Embryonic Genome Activation in Yellow Catfish (Pelteobagrus fulvidraco)

Paternal factors play an important role in embryonic morphogenesis and contribute to sexual dimorphism in development. To assess the effect of paternal DNA on sexual dimorphism of embryonic genome activation, we compared X and Y sperm and different sexes of embryos before sex determination. Through transcriptome sequencing (RNA-seq) and whole-genome bisulfite sequencing (WGBS) of X and Y sperm, we found a big proportion of upregulated genes in Y sperm, supported by the observation that genome-wide DNA methylation level is slightly lower than in X sperm. Cytokine-cytokine receptor interaction, TGF-beta, and toll-like receptor pathways play important roles in spermatogenesis. Through whole-genome re-sequencing (WGRS) of parental fish and RNA-seq of five early embryonic stages, we found the low-blastocyst time point is a key to maternal transcriptome degradation and zygotic genome activation. Generally, sexual differences emerged from the bud stage. Moreover, through integrated analysis of paternal SNPs and gene expression, we evaluated the influence of paternal inheritance on sexual dimorphism of genome activation. Besides, we screened out gata6 and ddx5 as potential instructors for early sex determination and gonad development in yellow catfish. This work is meaningful for revealing the molecular mechanisms of sex determination and sexual dimorphism of fish species.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

High-temperature stress will put the thermo-sensitive teleost yellow catfish (Tachysurus fulvidraco) in danger through reducing reproductivity

Recently, concerns for species that sex differentiation is influenced by temperature in the context of global warming have increased because disrupted operational sex ratios could threaten population maintenance. In contrast, little attention has been given to the reproductive ability of populations that experienced elevated temperatures. In this study, we demonstrated that high temperature (HT) would decrease population size via three different aspects of reproductive ability for the first time. We show that, in a thermo-sensitive teleost yellow catfish, a short period of HT (+3 °C) exposure during the critical period of sex differentiation leads to a different percentage of masculinization of XX genotypic females (1-23%) in wet-lab and natural water bodies. Combining the results of gonadal appearance, histology, sperm parameters, and fertilization rate, we found that XX pseudo-males induced by HT display significantly discounted fertility and reproductive performance compared to XY normal males. We demonstrate that the survival of the XY genotype is lower than XX genotype under environmental stress, including HT, hypoxia, and parasite infection, and the differential survival seems unrelated to male-biased sexual size dimorphism. The mathematical model predicts that the phenotypic female percent will be stabilized at 50% and the population will be sustainably maintained when masculinizing force is less than 0.5, while HT will put the population in danger when the masculinizing force exceeds 0.5. However, when we combine the real-world data of reproductive ability and mathematic model, our results suggest the population size decreases and the long-term survival of the studied species are threatened under the projected pace of increasing temperature. These findings will be useful for understanding the long-term effects of increasing temperature on sex ratio, reproduction and population maintenance in teleost.

Insights into chromosomal evolution and sex determination of Pseudobagrus ussuriensis (Bagridae, Siluriformes) based on a chromosome-level genome

Pseudobagrus ussuriensis is an aquaculture catfish with significant sexual dimorphism. In this study, a chromosome-level genome with a size of 741.97 Mb was assembled for female P. ussuriensis. A total of 26 chromosome-level contigs covering 97.34% of the whole-genome assembly were obtained with an N50 of 28.53 Mb and an L50 of 11. A total of 24,075 protein-coding genes were identified, with 91.54% (22,039) genes being functionally annotated. Based on the genome assembly, four chromosome evolution clusters of catfishes were identified and the formation process of P. ussuriensis chromosomes was predicted. A total of 55 sex-related quantitative trait loci (QTLs) with a phenotypic variance explained value of 100% were located on chromosome 8 (chr08). The QTLs and other previously identified sex-specific markers were located in a sex-determining region of 16.83 Mb (from 6.90 to 23.73 Mb) on chr08, which was predicted as the X chromosome. The sex-determining region comprised 554 genes, with 135 of which being differently expressed between males and females/pseudofemales, and 16 candidate sex-determining genes were screened out. The results of this study provided a useful chromosome-level genome for genetic, genomic and evolutionary studies of P. ussuriensis, and also be useful for further studies on sex-determination mechanism analysis and sex-control breeding of this fish.

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.

A time-course transcriptome analysis of gonads from yellow catfish (Pelteobagrus fulvidraco) reveals genes associated with gonad development

Background

The yellow catfish, Pelteobagrus fulvidraco, is a commercially important fish species. It is widely distributed in the fresh water areas of China, including rivers, lakes, and reservoirs. Like many other aquaculture fish species, people have observed significant size dimorphism between male and female yellow catfish and it shows a growth advantage in males.

Results

Here, at the first time, the time-course transcriptome was used to explore the various expression profiles of genes in different gonad developmental stages and genders. A total of 2696 different expression genes (DEGs) were identified from different stages. Based on these DEGs, 13 gonad development related genes were identified which showed time-specific or sex biased expression patterns.

Conclusion

This study will provide the crucial information on the molecular mechanism of gonad development of female and male yellow catfish. Especially, during the different gonad development stages, these 13 gonad development related genes exhibit various expression patterns in female and male individual respectively. These results could inspire and facilitate us to understanding the various roles of these genes play in different gonad development stages and genders.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08651-0.

Chromosome-level genome assembly of Asian yellow pond turtle (Mauremys mutica) with temperature-dependent sex determination system

Knowledge of sex determination has important implications in physiology, ecology and genetics, but the evolutionary mechanisms of sex determination systems in turtles have not been fully elucidated, due to a lack of reference genomes. Here, we generate a high-quality genome assembly of Asian yellow pond turtle (Mauremys mutica) using continuous long-read (PacBio platform), Illumina, and high-throughput chromatin conformation capture (Hi-C) technologies. The M. mutica haplotype has a genome size of 2.23 Gb with a contig N50 of 8.53 Mb and scaffold N50 of 141.98 Mb. 99.98% sequences of the total assembly are anchored to 26 pseudochromosomes. Comparative genomics analysis indicated that the lizard-snake-tuatara clade diverged from the bird-crocodilian-turtle clade at approximately 267.0-312.3 Mya. Intriguingly, positive selected genes are mostly enriched in the calcium signaling pathway and neuroactive ligand-receptor interaction, which are involved in the process of temperature-dependent sex determination. These findings provide important evolutionary insights into temperature-dependent sex determination system.

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.

An ancient truncated duplication of the anti‐Mullerian hormone receptor type 2 gene is a potential conserved master sex determinant in the Pangasiidae catfish family

The evolution of sex determination (SD) in teleosts is amazingly dynamic, as reflected by the variety of different master sex-determining genes identified. Pangasiids are economically important catfishes in South Asian countries, but little is known about their SD system. Here, we generated novel genomic resources for 12 Pangasiids and characterized their SD system. Based on a Pangasianodon hypophthalmus chromosome-scale genome assembly, we identified an anti-Müllerian hormone receptor type Ⅱ gene (amhr2) duplication, which was further characterized as being sex-linked in males and expressed only in testes. These results point to a Y chromosome male-specific duplication (amhr2by) of the autosomal amhr2a. Sequence annotation revealed that the P. hypophthalmus Amhr2by is truncated in its N-terminal domain, lacking the cysteine-rich extracellular part of the receptor that is crucial for ligand binding, suggesting a potential route for its neofunctionalization. Reference-guided assembly of 11 additional Pangasiids, along with sex-linkage studies, revealed that this truncated amhr2by duplication is a male-specific conserved gene in Pangasiids. Reconstructions of the amhr2 phylogeny suggested that amhr2by arose from an ancient duplication/insertion event at the root of the Siluroidei radiation that is dated to ~100 million years ago. Together these results bring multiple lines of evidence supporting that amhr2by is an ancient and conserved master sex-determining gene in Pangasiids, a finding that highlights the recurrent use of the transforming growth factor β pathway, which is often used for the recruitment of teleost master SD genes, and provides another empirical case towards firther understanding of dynamics of SD systems.

Screening and identification of female-specific DNA sequences in octaploid sturgeon using comparative genomics with high-throughput sequencing

In this study, six candidate female-specific DNA sequences of octaploid Amur sturgeon (Acipenser schrenckii) were identified using comparative genomic approaches with high-throughput sequencing data. Their specificity was confirmed by traditional PCR. Two of these sex-specific sequences were also validated as female-specific in other eight sturgeon species and two hybrid sturgeons. The identified female-specific DNA fragments suggest that the family Acipenseridae has a ZZ/ZW sex-determining system. However, one of the two DNA sequences has been deleted in some sturgeons such as Sterlet sturgeon (Acipenser ruthenus), Beluga (Huso huso) and Kaluga (H. dauricus). The difference of sex-specific sequences among sturgeons indicates that there are different sex-specific regions among species of sturgeon. This study not only provided the sex-specific DNA sequences for management, conservation and studies of sex-determination mechanisms in sturgeons, but also confirmed the capability of the workflow to identify sex-specific DNA sequences in the polyploid species with complex genomes.

Arginase plays an important role in ammonia detoxification of yellow catfish Pelteobagrus fulvidraco

A two-stage study was carried out to test the mechanism of arginase in ammonia detoxification of yellow catfish. At stage 1, fish was injected lethal half concentration ammonium acetate and 0.9% sodium chloride respectively every 12 h in six replicates for 72 h. The result found that no significant different in serum ammonia contents of fish in ammonium acetate group at hours 12, 24, 36, 48, 60 and 72. At stage 2, ammonium acetate group was split in two, one continued to injected with ammonium acetate (NH₃ group) and the other with ammonium acetate and valine (an inhibitor of arginase; Val group); Sodium chloride group also was split in two, one continued to injected with sodium chloride (NaCl group) and the other with sodium chloride and valine (NaCl + Val group). The experiment continued for 12 h. Serum ammonia and liver arginine contents of fish in Val group were higher than those of fish in NH₃ group; Compared with NaCl group, arginase activity and ARG 1 expression in liver of fish in Val group were lower; Fish in NaCl and NaCl + Val groups had the lowest serum superoxide dismutase activities, malondialdehyde, tumor necrosis factor-α, interleukin 1 and 8 contents, TNF-α, IL-1 and IL-8 expressions than fish in NH₃ and Val groups, and had the higher lysozyme activities, complement 3 and 4 contents. This study indicates that ammonia poisoning would lead to oxidative damage, immunosuppression and inflammation in yellow catfish; Arginase may be an important target of ammonia toxicity in yellow catfish; Exogenous arginine supplementation might alleviate the symptoms of ammonia poisoning in yellow catfish.

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

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

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.

A high-density SNP-based genetic map and several economic traits-related loci in Pelteobagrus vachelli

BACKGROUND

A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection in aquaculture species. Pelteobagrus vachelli is a very popular commercial species in Asia. However, some specific characters hindered achievement of the traditional selective breeding based on phenotypes, such as lack of large-scale genomic resource and short of markers tightly associated with growth, sex determination and hypoxia tolerance related traits.

RESULTS

By making use of 5059 ddRAD markers in P. vachelli, a high-resolution genetic linkage map was successfully constructed. The map' length was 4047.01 cM by using an interval of 0.11 cm, which is an average marker standard. Comparative genome mapping revealed that a high proportion (83.2%) of markers with a one-to-one correspondence were observed between P. vachelli and P. fulvidraco. Based on the genetic map, 8 significant genome-wide QTLs for 4 weight, 1 body proportion, 2 sex determination, and 1 hypoxia tolerance related traits were detected on 4 LGs. Some SNPs from these significant genome-wide QTLs were observably associated with these phenotypic traits in other individuals by Kompetitive Allele Specific PCR. In addition, two candidate genes for weight, Sipa1 and HSD11B2, were differentially expressed between fast-, medium- and slow-growing P. vachelli. Sema7a, associated with hypoxia tolerance, was induced after hypoxia exposure and reoxygenation.

CONCLUSIONS

We mapped a set of suggestive and significant QTLs as well as candidate genes for 12 growth, 1 sex determination and 1 hypoxia tolerance related traits based on a high-density genetic linkage map by making use of SNP markers for P. fulvidraco. Our results have offered a valuable method about the much more efficient production of all-male, fast growth and hypoxia tolerance P. vachelli for the aquaculture industry.

Synergistic Combination of Exogenous Hormones to Improve the Spawning and Post-spawning Survival of Female Yellow Catfish

Multiple repeat spawners make large contributions to long-term population stability and aquaculture breeding programs. A high percentage of female yellow catfish (Pelteobagrus fulvidraco) died for spawning failure or incomplete spawning after artificial spawning by traditional synthetic hormones including human chorionic gonadotropin (hCG), luteinizing hormone releasing hormone (LHRH), and domperidone (DOM). The present study was designed to compare the efficacy of different combinations of exogenous hormones for inducing ovulation in yellow catfish using hCG, LHRH, DOM, and carp pituitary extraction (CPE). We found a optimal strategy for exogenous hormones administration, the mixture of LHRH/CPE for the first injection and LHRH/CPE/DOM for the second injection, could greatly improve the rates of spawning success, weight of ovulated eggs and survival rate after spawning. Interestingly, a population of female yellow catfish with defective reproductive duct could not spawn and showed high mortality after induced by a combination of hCG/LHRH/DOM, whereas a synergistic combination of hCG, LHRH, DOM, and CPE could efficiently induce spawning and reduce mortality in the defective yellow catfish, in which a significant decrease of Vitellin and E2 levels. Altogether, our findings provide an effective combination of exogenous hormones to improve spawning and post-spawning survival of female yellow catfish.

Identifying difference in primordial germ cells between XX female and XY male yellow catfish embryos

Primordial germ cells (PGCs) are singled out from somatic cells very early during embryogenesis, then they migrate towards the genital ridge and differentiate into gametes through oogenesis or spermatogenesis. Labeling PGCs with Localized RNAexpression (LRE) technique by fluorescent proteins has been widely applied among teleost species to study the germ cell development and gonad differentiation. In this study, we first cloned and characterized the 3' untranslated regions (3'UTRs) of nanos homolog 1-like (nos1l), dead end (dnd), and vasa in yellow catfish (Pelteobagrus fulvidraco), and then synthesized the GFP-nos1l/dnd/vasa 3'UTR mRNAs. Each of these three 3'UTRs could label PGCs in yellow catfish embryos, of which, vasa 3'UTR exhibited the highest labeling efficiency. To identify the differences in PGCs at embryonic stage, XX all-female and XY all-male yellow catfish embryos were produced and injected with GFP-vasa 3'UTR mRNA. We observed the PGC migration route in these two monosex embryos from 24 hpf to 7 dpf, and found there was no difference between them. Besides, the PGC number was counted at 48 hpf, and the result showed that the average PGC number in XX females (11.3) was significantly larger than that in XY males (8.1).These findings provide an insight into the development of PGCs in yellow catfish embryos and the relationship between embryonicPGCnumberand thelatergonaddifferentiation.

Cloning, expression pattern, and potential role of MAPKp38 and NF-κBp65 in response to lipopolysaccharide in yellow catfish (Pelteobagrus fulvidraco)

Mitogen-activated protein kinases (MAPK) and nuclear factor kappa B (NF-κB) pathways are considered to be two crucial intracellular signaling cascades in pro-inflammatory responses. In this study, we reported the coding sequences (CDS) of MAPKp38 and NF-κBp65 from yellow catfish. We also investigated the gene structure, expression patterns and functional role in yellow catfish. The CDS of MAPKp38 is 1,086 bp encoding 361 amino acids (AA). The MAPKp38 protein has a long highly conserved serine/threonine protein kinases catalytic domain. The NF-κBp65 CDS is 1,794 bp, and the gene encodes 597 AA, with a Rel homology domain (RHD) which consists of a RHD-DNA-binding domain and an Ig-like, plexins, transcription factors (IPT) domain. Moreover, MAPKp38 and NF-κBp65 protein of bony fish and other vertebrates have a single clade. Quantitative real-time PCR analysis revealed the presence of the MAPKp38 and NF-κBp65 transcript in 12 tissues of healthy yellow catfish. The highest expression levels of MAPKp38 and NF-κBp65 were detected in the heart and liver, respectively. Upon stimulation with an intraperitoneal injection of lipopolysaccharide (LPS), the expression levels of MAPKp38 and NF-κBp65 were up-regulated in the intestine. These results indicated that MAPKp38 and NF-κBp65 play important roles in mediating the response protection against LPS in yellow catfish.

Comparative Transcriptome Analysis Reveals Differentially Expressed Genes and Signaling Pathways Between Male and Female Red-Tail Catfish (Mystus wyckioides)

Sexual dimorphism is widespread in fish species. The red-tail catfish (Mystus wyckioides) is a commercially important catfish in the lower reaches of the Lancang River and the Mekong basin, and it shows a growth advantage in males. Here, RNA-seq was for the first time used to explore the gene expression difference between the sexes in the hypothalamus and pituitary of red-tail catfish, respectively. In the hypothalamus, 5732 and 271 unigenes have significantly higher and lower expressions, respectively, in males compared with females. KEGG analysis showed that 212 DEGs were annotated to 216 signaling pathways, and enrichment analysis suggested different levels of cAMP and glutamatergic synapse signaling between male and female hypothalami and some of the DEGs appear involved in gonad development and growth. In the pituitary, we found only 19 differentially expressed unigenes, which were annotated to 32 signaling pathways, most of which play important roles in gonad development.

Development and Validation of Sex-Specific Markers in Pelodiscus Sinensis Using Restriction Site-Associated DNA Sequencing

The sex of an animal influences its economic traits, especially in species displaying sexual dimorphism. The Chinese soft-shelled turtle, Pelodiscus sinensis, is an economically important aquatic species that shows significant male sexual dimorphism, with a large body size, faster growth, a thick and wide calipash, and lower body fat. In this study, ten male and ten female turtles were subjected to restriction site-associated DNA sequencing (RAD-seq) using the Hi-Seq 4000 sequencing platform to isolate female-specific DNA fragments. We identified 5967 bp and 6532 bp fragments using genome walking. Three female-specific markers designed from these two fragments were confirmed to separate the sexes of Pelodiscus sinensis perfectly. One of the female-specific markers showed dosage association in female and male individuals. Individuals from different populations (n = 296) were used to validate that the female-specific markers could identify the genetic sex of Pelodiscus sinensis with 100% accuracy. The results of the present study demonstrated that RAD-seq was useful to develop sex-related markers in animals, and verified that the sex determination system of Pelodiscus sinensis belonged to the ZZ/ZW heterogametic system. Importantly, the developed markers could lead to a method for sex-controlled breeding in the Chinese soft-shelled turtle.

Construction of a Genomic Bacterial Artificial Chromosome (BAC) Library for the Prawn Macrobrachium rosenbergii and Initial Analysis of ZW Chromosome-Derived BAC Inserts

Knowledge on sex determination has proven valuable for commercial production of the prawn Macrobrachium rosenbergii due to sex dimorphism of the male and female individuals. Previous studies indicated that prawn sex is determined by a ZW-ZZ chromosomal system, but no genomic information is available for the sex chromosome. Herein, we constructed a genomic bacterial artificial chromosome (BAC) library and identified the ZW-derived BAC clones for initial analysis of the sex chromosomal DNA sequence. The arrayed BAC library contains 200,448 clones with average insert size of 115.4 kb, corresponding to ∼ 4× coverage of the estimated 5.38 Gb genome. Based on a short female-specific marker, a Z- and a W-fragment were retrieved with the genomic walking method. Screening the BAC library using a ZW-specific marker as probe resulted in 12 positive clones. From these, a Z-derived (P331M17) and a W-derived (P122G2) BAC clones were randomly selected and sequenced by PacBio method. We report the construction of a large insert, deep-coverage, and high-quality BAC library for M. rosenbergii that provides a useful resource for positional cloning of target genes, genomic organization, and comparative genomics analysis. Our study not only confirmed the ZW/ZZ system but also discovered sex-linked genes on ZW chromosomes for the first time, contributing to a comprehensive understanding of the genomic structure of sex chromosomes in M. rosenbergii.

Identification of male-specific SNP markers and development of PCR-based genetic sex identification technique in crucifix crab (Charybdis feriatus) with implication of an XX/XY sex determination system

In this study, we first identified male-specific SNP markers using restriction site-associated DNA sequencing, and further developed a PCR-based sex identification technique for Charybdis feriatus. A total of 296.96 million clean reads were obtained, with 114.95 and 182.01 million from females and males. After assembly and alignment, 10 SNP markers were identified being heterozygous in males but homozygous in females. Five markers were further confirmed to be male-specific in a large number of individuals. Moreover, two male-specific sense primers and a common antisense primer were designed, using which, a PCR-based genetic sex identification method was successfully developed and used to identify the sex of 103 individuals, with a result of 49 females and 54 males. The presence of male-specific SNP markers suggests an XX/XY sex determination system for C. feriatus. These findings should be helpful for better understanding sex determination mechanism, and drafting artificial breeding program in crustaceans.

Whole-Genome Sequencing of Chinese Yellow Catfish Provides a Valuable Genetic Resource for High-Throughput Identification of Toxin Genes

Naturally derived toxins from animals are good raw materials for drug development. As a representative venomous teleost, Chinese yellow catfish (Pelteobagrus fulvidraco) can provide valuable resources for studies on toxin genes. Its venom glands are located in the pectoral and dorsal fins. Although with such interesting biologic traits and great value in economy, Chinese yellow catfish is still lacking a sequenced genome. Here, we report a high-quality genome assembly of Chinese yellow catfish using a combination of next-generation Illumina and third-generation PacBio sequencing platforms. The final assembly reached 714 Mb, with a contig N50 of 970 kb and a scaffold N50 of 3.65 Mb, respectively. We also annotated 21,562 protein-coding genes, in which 97.59% were assigned at least one functional annotation. Based on the genome sequence, we analyzed toxin genes in Chinese yellow catfish. Finally, we identified 207 toxin genes and classified them into three major groups. Interestingly, we also expanded a previously reported sex-related region (to ≈6 Mb) in the achieved genome assembly, and localized two important toxin genes within this region. In summary, we assembled a high-quality genome of Chinese yellow catfish and performed high-throughput identification of toxin genes from a genomic view. Therefore, the limited number of toxin sequences in public databases will be remarkably improved once we integrate multi-omics data from more and more sequenced species.

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

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

Diverse and variable sex determination mechanisms in vertebrates

Sex is prevalent in nature and sex determination is one of the most fundamental biological processes, while the way of initiating female and male development exhibits remarkable diversity and variability across vertebrates. The knowledge on why and how sex determination mechanisms evolve unusual plasticity remains limited. Here, we summarize sex determination systems, master sex-determining genes and gene-regulatory networks among vertebrates. Recent research advancements on sex determination system transition are also introduced and discussed in some non-model animals with multiple sex determination mechanisms. This review will provide insights into the origin, transition and evolutionary adaption of different sex determination strategies in vertebrates, as well as clues for future perspectives in this field.

Chromosomal-level assembly of yellow catfish genome using third-generation DNA sequencing and Hi-C analysis

Background

The yellow catfish, Pelteobagrus fulvidraco, belonging to the Siluriformes order, is an economically important freshwater aquaculture fish species in Asia, especially in Southern China. The aquaculture industry has recently been facing tremendous challenges in germplasm degeneration and poor disease resistance. As the yellow catfish exhibits notable sex dimorphism in growth, with adult males about two- to three-fold bigger than females, the way in which the aquaculture industry takes advantage of such sex dimorphism is another challenge. To address these issues, a high-quality reference genome of the yellow catfish would be a very useful resource.

Findings

To construct a high-quality reference genome for the yellow catfish, we generated 51.2 Gb short reads and 38.9 Gb long reads using Illumina and Pacific Biosciences (PacBio) sequencing platforms, respectively. The sequencing data were assembled into a 732.8 Mb genome assembly with a contig N50 length of 1.1 Mb. Additionally, we applied Hi-C technology to identify contacts among contigs, which were then used to assemble contigs into scaffolds, resulting in a genome assembly with 26 chromosomes and a scaffold N50 length of 25.8 Mb. Using 24,552 protein-coding genes annotated in the yellow catfish genome, the phylogenetic relationships of the yellow catfish with other teleosts showed that yellow catfish separated from the common ancestor of channel catfish ∼81.9 million years ago. We identified 1,717 gene families to be expanded in the yellow catfish, and those gene families are mainly enriched in the immune system, signal transduction, glycosphingolipid biosynthesis, and fatty acid biosynthesis.

Conclusions

Taking advantage of Illumina, PacBio, and Hi-C technologies, we constructed the first high-quality chromosome-level genome assembly for the yellow catfish P. fulvidraco. The genomic resources generated in this work not only offer a valuable reference genome for functional genomics studies of yellow catfish to decipher the economic traits and sex determination but also provide important chromosome information for genome comparisons in the wider evolutionary research community.

Origin and transition of sex determination mechanisms in a gynogenetic hexaploid fish

Most vertebrates reproduce sexually, and plastic sex determination mechanisms including genotypic sex determination (GSD) and environmental sex determination (ESD) have been extensively revealed. However, why sex determination mechanisms evolve diversely and how they correlate with diverse reproduction strategies remain largely unclear. Here, we utilize the superiority of a hexaploid gibel carp (Carassius gibelio) that is able to reproduce by unisexual gynogenesis and contains a rare but diverse proportion of males to investigate these puzzles. A total of 2248 hexaploid specimens were collected from 34 geographic wild populations throughout mainland China, in which 24 populations were revealed to contain 186 males with various incidences ranging from 1.2 to 26.5%. Subsequently, the proportion of temperature-dependent sex determination (TSD) was revealed to be positively correlated to average annual temperature in wild populations, and male incidence in lab gynogenetic progenies was demonstrated to increase with the increasing of larval rearing temperature. Meanwhile, extra microchromosomes were confirmed to play genotypic male determination role as previously reported. Thereby, GSD and TSD were found to coexist in gibel carp, and the proportions of GSD were observed to be much higher than that of TSD in sympatric wild populations. Our findings uncover a potential new mechanism in the evolution of sex determination system in polyploid vertebrates with unisexual gynogenesis ability, and also reveal a possible association of sex determination mechanism transition between TSD and GSD and reproduction mode transition between unisexual gynogenesis and bisexual reproduction.

Stat5b Regulates Sexually Dimorphic Gene Expression in Zebrafish Liver

Sexual size dimorphism is an interesting phenomenon occurred in many fish species. Wildtype zebrafish exhibits a significant sexual dimorphism in body size at the adult stage. Previous studies indicated that sexual size dimorphism was eliminated in stat5b-mutated zebrafish. Herein, the comparative transcriptome analysis was conducted to observe the genes and pathways involved in sexual size dimorphism. The number of male-biased and female-biased genes was much less in the liver of stat5b mutant zebrafish than in wildtype. Gene ontology (GO) enrichment and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that multiple pathways related to metabolism were affected upon loss of stat5b function. qRT-PCR results also validated that sexually dimorphic expression of a set of genes was lost when stat5b was mutated. Furthermore, the weighted correlation network analysis (WGCNA) detected many candidate genes related to the growth traits and stat5b function, such as greb1, lepr, and igf2b. Our data suggest that stat5b should regulate the sexually dimorphic gene expression in zebrafish liver and add in understanding of the molecular mechanisms underlying sexual size dimorphism in fish species.

Sex-specific markers developed by next-generation sequencing confirmed an XX/XY sex determination system in bighead carp (Hypophthalmichehys nobilis) and silver carp (Hypophthalmichthys molitrix)

Sex-specific markers are powerful tools for identifying sex-determination system in various animals. Bighead carp (Hypophthalmichehys nobilis) and silver carp (Hypophthalmichthys molitrix) are two of the most important edible fish in Asia, which have a long juvenility period that can lasts for 4-5 years. In this study, we found one sex-specific marker by next-generation sequencing together with bioinformatics analysis in bighead carp. The male-specific markers were used to perform molecular sexing in the progenies of artificial gynogenetic diploids and found all progenies (n = 160) were females. Meanwhile, around 1 : 1 sex ratio was observed in a total of 579 juvenile offspring from three other families. To further extend the male-specific region, we performed genome walking and got a male-specific sequence of 8,661 bp. Five pairs of primers were designed and could be used to efficiently distinguish males from females in bighead carp and silver carp. The development of these male-specific markers and results of their molecular sexing in different populations provide strong evidence for a sex determination system of female homogametry or male heterogametry (XX/XY) in bighead carp and silver carp. To the best of our knowledge, this is the first report of effective sex-specific markers in these two large carp species.

Potential Contributions of miR-200a/-200b and Their Target Gene-Leptin to the Sexual Size Dimorphism in Yellow Catfish

Sexual size dimorphism is the consequence of differential expression of sex-biased genes related to feeding and growth. Leptin is known to regulate energy balance by regulating food intake. In order to investigate the molecular mechanism of sexual size dimorphism in yellow catfish (Pelteobagrus fulvidraco), the expression of leptin (lep) and its functional receptor (lepr) were detected during larval development. Both lep and lepr have lower expression in males than in females during 1–4 weeks post hatching. 17a-Methyltestosterone (MT) treatment resulted in decreased expression of lep and lepr in both male and female larval fish. Interestingly, the mRNA levels of lep and lepr in juvenile male were significantly decreased compared with juvenile female during short-term fasting periods. Lep was predicted to be a potential target of miR-200a and miR-200b that had an opposite expression pattern to lep in male and female larvas. The results of luciferase reporter assay suggested that lep is a target of miR-200a/-200b. Subsequently, male hormone and fasting treatment have opposite effects on the expression of miR-200a/-200b and lep between males and females. In summary, our results suggest that sexual size dimorphism in fish species is probably caused by the sexually dimorphic expression of leptin, which could be negatively regulated by miR-200a/-200b.

Essential roles of stat5.1/stat5b in controlling fish somatic growth

Signal transducer and activator of transcription 5b (STAT5b) has been identified as a key downstream mediator of growth hormone (GH) signaling in somatic growth of mammalian. However, the corresponding homologue gene of Stat5b is unknown in fish species. In this study, we generated loss-of-function mutants in stat5.1 and stat5.2, two stat5 homologues existing in zebrafish. In stat5.1-deficient zebrafish, a significant reduction of body length and body weight was detected in the embryos/larvae and adults compared with the wild-type control fish, and sexual size dimorphism in adult zebrafish was also eliminated. However, the stat5.2-deficient zebrafish displayed a normal developmental phenotype during all lifespan. Chromatin immunoprecipitation combined with deep sequencing (ChIP-seq) method was adopted to further investigate the potential transcriptional targets of Stat5 protein and cast much light upon the biological function of Stat5. We identified more than 800 genes as transcriptional targets of Stat5 during zebrafish embryogenesis. KEGG analysis indicated that the Stat5 target gene network is predominantly linked to the metabolic pathways, neuroactive ligand-receptor interaction and JAK-STAT signaling pathways. Further validation studies suggested that Stat5.1 protein could directly regulate the expression of gh1, and stat5.1-mutated zebrafish showed a reduction of gh1 mRNA level. In the present study, stat5.1 was revealed as the corresponding homologue gene of Stat5b in fish species. Additionally, we found a novel molecular interaction between Stat5.1/Stat5b and GH, and unraveled a positive feedback loop Stat5.1-GH-Stat5.1 which is necessary for somatic growth and body development in zebrafish.

Deciphering transcriptome profile of the yellow catfish (Pelteobagrus fulvidraco) in response to Edwardsiella ictaluri

Edwardsiella ictaluri is one of the most important pathogens posing a serious threat for yellow catfish (Pelteobagrus fulvidraco), a highly valuable fish species of increasing commercial interest in China. Here, a transcriptomic strategy was undertaken to investigate the yellow catfish gene expression profile against infection by the bacterial pathogen E. ictaluri. Comparison of the transcriptome profiles between the infected and uninfected samples showed that a massive gene expression change occurred in yellow catfish following bacterial exposure. A total of 5527 differentially expressed genes (DEGs) were detected, of which 2265 showed up-regulation and 3262 down-regulation. Gene set enrichment analysis revealed the presence of canonical pathways directly linked to innate and adaptive immune response, such as pattern recognition receptor (PRR) signaling pathways, complement and coagulation cascades, as well as T-cell receptor (TCR) and B-cell receptor (BCR) signaling pathways. Additionally, 47,526 putative EST-liked simple sequence repeats (SSRs) markers were retrieved for use in genetic studies. This study establishes the first molecular clues to understand the potential mechanisms of yellow catfish resistance to E. ictaluri, thus enabling future efforts on disease control programs in this valuable aquaculture species.

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

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

Gonadal morphogenesis and sex differentiation in cultured Ussuri catfish Tachysurus ussuriensis

The objective of this study was to investigate the optimal developmental time to perform sex reversal in Ussuri catfish Tachysurus ussuriensis, to develop monosex breeding in aquaculture. Systematic observations of gonadal sex differentiation of P. ussiriensis were conducted. The genital ridge formed at 9 days post fertilization (dpf) and germ cells begin to proliferate at 17 dpf. The ovarian cavity began forming on 21 dpf and completed by 25 dpf while presumptive testis remained quiescent. The primary oocytes were at the chromatin nucleolus stage by 30 dpf, the peri-nucleolus stage by 44 dpf and the cortical alveoli stage by 64 dpf. The germinal vesicle migrated towards the animal pole (polarization) at 120 dpf. In presumptive testis, germ cells entered into mitosis and blood vessels appeared in the proximal gonad on 30 dpf. The efferent duct anlage appeared on 36 dpf and formation of seminal lobules with spermatogonia and lobules interstitium occurred at 120 dpf. Therefore, gonadal sex differentiation occurred earlier in females than in males, with the histological differentiation preceding cytologic differentiation in T. ussuriensis. This indicates that undifferentiated gonads directly differentiate into ovary or testis between 17 and 21 dpf and artificial induction of sexual reversal by oral steroid administration must be conducted before 17 dpf.

Potential contributions of heat shock proteins and related genes in sexual differentiation in yellow catfish (Pelteobagrus fulvidraco)

Sex determination and differentiation in ectotherms are very complicated affairs and usually affected by both genetic and environmental factors. Because of their temperature-sensitive expression, heat shock proteins (HSPs) are good candidates for temperature-dependent sex determination (TSD). Similar to most thermosensitive fish species, the male to female ratio increases with temperature in yellow catfish (Pelteobagrus fulvidraco). Yellow catfish is also a type of sexual size dimorphic fish, and the male individuals grow much faster than females of the same age. Therefore, research of sex differentiation in yellow catfish is important in aquiculture. In this attempt, a total of seven HSPs and related genes were identified from transcriptomes of yellow catfish by 454 pyrosequencing and Solexa sequencing that we did previously, including five genes with complete open reading frame (ORF). Phylogenetically, all these genes were compared with their counterparts from other vertebrates. All these genes were sex-biased expressed in gonads. Hspa5, Hip, and Cdc37 were expressed more highly in ovary than in testis, whereas Hsp90α, Hspb2, Hspb8, and Hspbp1 were expressed more highly in testis than in ovary. Additionally, the expression of these genes was assessed after 17α-methyltestosterone (MT) and 17α-ethinylestradiol (EE2) treatment, respectively. Our result showed that working as co-chaperones, these HSPs and related genes may regulate sex steroid receptor activities to influence gonad development in yellow catfish. Our work would help in the understanding of the mechanism of sexual differentiation in teleosts.

Molecular cloning and characterization of a C-type lectin in yellow catfish Tachysurus fulvidraco

This study represents the first report of a C-type lectin (ctl) in yellow catfish Tachysurus fulvidraco. The complete sequence of ctl complementary (c)DNA consisted of 685 nucleotides. The open reading frame potentially encoded a protein of 177 amino acids with a calculated molecular mass of c.y 20.204 kDa. The deduced amino-acid sequence contained a signal peptide and a single carbohydrate recognition domain with four cysteine residues and GlnProAsp (QPD) and TrpAsnAsp (WND) motifs. Ctl showed the highest identity (56.0%) to the predicted lactose binding lectin from channel catfish Ictalurus punctatus. Quantitative real-time (qrt)-PCR analysis showed that ctl messenger (m)RNA was constitutively expressed in all examined tissues in normal fish, with high expression in trunk kidney and head kidney, which was increased following Aeromonas hydrophila challenge in a duration-dependent manner. Purified recombinant Ctl (rCtl) from Escherichia coli BL21 was able to bind and agglutinate Gram-positive and Gram-negative bacteria in a calcium-dependent manner. These results suggested that Ctl might be a C-type lectin of T. fulvidraco involved in innate immune responses as receptors (PRR).

Effects of astaxanthin and emodin on the growth, stress resistance and disease resistance of yellow catfish (Pelteobagrus fulvidraco)

Yellow catfish (Pelteobagrus fulvidraco) has become a commercially important fish species in China and eastern Asia. High-density aquaculture has led to congestion and excessive stress and contributed to bacterial infection outbreaks that have caused high mortality. We investigated the effects of dietary supplementation with astaxanthin and emodin alone and in combination on the growth and stress resistance of yellow catfish. After 60 days of feeding, each group of fish (control, astaxanthin, emodin, and astaxanthin plus emodin (combination) groups) was exposed to acute crowding stress for 24 h, and a subsample of fish from the four groups was challenged with the bacterial septicemia pathogen Proteus mirabilis after the end of the crowding stress experiment. Compared with the control, the astaxanthin and emodin groups showed increases in serum total protein (TP), hepatic superoxide dismutase (SOD) activity and hepatic heat shock proteins 70 (HSP70) mRNA levels at 12 and 24 h after the initiation of crowding stress. The combination group exhibited increases in alanine aminotransferase (ALT) activity, aspartate aminotransferase (AST) activity, serum TP, hepatic SOD activity and hepatic HSP70 mRNA levels within 24 h after the initiation of crowding stress. However, decreases relative to the control were observed in the serum cortisol and glucose contents in the three treatment groups at 12 and 24 h after the initiation of crowding stress, in ALT and AST activity in the astaxanthin and emodin group at 24 h after the initiation of crowding stress, and in the serum lysozyme activity, serum alkaline phosphatase (ALP) activity, and hepatic catalase (CAT) and malondialdehyde (MDA) activity in the combination group at 24 h after the initiation of crowding stress. Additionally, the cumulative mortality after P. mirabilis infection was lower in all three treatment groups (57.00%-70.33%) than in the control (77.67%). Dietary supplementation with astaxanthin and emodin decreased the specific growth rate (SGR) and weight gain (WG) of healthy yellow catfish, although significant differences in mortality were not observed. These results indicate that dietary supplementation with 80 mg/kg astaxanthin and 150 mg/kg emodin can improve the anti-oxidative capabilities, hepatic HSP70 levels, and resistance to acute crowding stress of yellow catfish. Finally, an appropriate strategy for enhance yellow catfish stress resistance and disease resistance is proposed.

Extra Microchromosomes Play Male Determination Role in Polyploid Gibel Carp

Sex is generally determined by sex chromosomes in vertebrates, and sex chromosomes exhibit the most rapidly-evolving traits. Sex chromosome evolution has been revealed previously in numerous cases, but the association between sex chromosome origin and the reproduction mode transition from unisexual to sexual reproduction remains unclear. Here, we have isolated a male-specific sequence via analysis of amplified fragment length polymorphism from polyploid gibel carp (Carassius gibelio), a species that not only has the ability to reproduce unisexually but also contains males in wild populations. Subsequently, we have found through FISH analysis that males have several extra microchromosomes with repetitive sequences and transposable elements when compared to females. Moreover, we produced sex-reversed physiological females with a male-specific marker by using estradiol hormone treatment, and two gynogenetic families were established from them. In addition, the male incidence rates of two gynogenetic families were revealed to be closely associated with the extra microchromosome number of the sex-reversed physiological females. These results suggest that the extra microchromosomes in males might resemble a common feature of sex chromosomes and might play a significant role in male determination during the evolutionary trajectory of the reproduction mode transition from unisexual to sexual reproduction in the polyploid fish.

Sex biased expression of ghrelin and GHSR associated with sexual size dimorphism in yellow catfish

Sexual size dimorphism has been observed in many cultivable fish species including yellow catfish, in which male fish grow much faster than female fish. Ghrelin is a potent stimulator of pituitary growth hormone (GH) release and known to potentially promote food intake and body weight gain. In order to investigate the molecular mechanism of sexual size dimorphism in yellow catfish (Pelteobagrus fulvidraco), ghrelin and its functional receptor, growth hormone secretagogue receptor (GHSR) cDNAs were cloned. Real-time PCR indicated that both ghrelin and GHSR were more highly expressed in hypothalamus and gut of male fish than female. During normal larval development, expression of ghrelin and GHSR genes was significantly higher in males than in females. 17a-Methyltestosterone (MT) treatment enhanced the expression of ghrelin in female larval fish and GHSR in both sexes, whereas the expression of ghrelin in male larval fish increased in the beginning, then decreased as the treatment time prolonged. Furthermore, the expression of ghrelin and GHSR in male juvenile was significantly increased compared with female juvenile, in short and long term fasting periods, suggesting that male fish may have a better appetite than female during fasting. Our results demonstrate that sex difference in the expression of ghrelin and GHSR may be involved in sexual size dimorphism by regulating feeding and GH/IGF signaling in yellow catfish.

Development of PCR-Based Markers to Determine the Sex of Kelps

Sex discriminating genetic markers are commonly used to facilitate breeding programs in economically and ecologically important animal and plant species. However, despite their considerable economic and ecological value, the development of sex markers for kelp species has been very limited. In this study, we used the recently described sequence of the sex determining region (SDR) of the brown algal model Ectocarpus to develop novel DNA-based sex-markers for three commercially relevant kelps: Laminaria digitata, Undaria pinnatifida and Macrocystis pyrifera. Markers were designed within nine protein coding genes of Ectocarpus male and female (U/V) sex chromosomes and tested on gametophytes of the three kelp species. Seven primer pairs corresponding to three loci in the Ectocarpus SDR amplified sex-specific bands in the three kelp species, yielding at least one male and one female marker for each species. Our work has generated the first male sex-specific markers for L. digitata and U. pinnatifida, as well as the first sex markers developed for the genus Macrocystis. The markers and methodology presented here will not only facilitate seaweed breeding programs but also represent useful tools for population and demography studies and provide a means to investigate the evolution of sex determination across this largely understudied eukaryotic group.