Artificial gynogenesis and sex determination in half-smooth tongue sole (Cynoglossus semilaevis)

The Chinese medaka (Oryzias sinensis) dmrt1 gene converts females to males in medaka (Oryzias latipes)

BACKGROUND

Chinese medaka (Oryzias sinensis) is widely distributed in freshwater rivers in China. Similar to the medaka (Oryzias latipes), Chinese medaka has the characteristics of small size, rapid reproductive cycle, and strong adaptability, which makes it suitable as a model organism for studies in basic biology and environmental toxicology. Chinese medaka exhibits distinct sexual dimorphism. However, due to the lack of complete genomic information, the regulation of sex determination and differentiation-related genes in Chinese medaka remains unclear.

METHODS

Chinese medaka dmrt1 (Osdmrt1) was cloned by PCR, and transgenic individuals of medaka [Tg(CMV:Osdmrt1)] overexpressing Osdmrt1 were generated to investigate the role of Osdmrt1 in sex determination. Western blot was used to validate the integration of the Osdmrt1 into the medaka genome. Tissue sectioning and HE staining were used to identify Tg(CMV:Osdmrt1) physiological gender and phenotype. qRT-PCR was used to analyze the expression of gonad-specific genes.

RESULTS

Osdmrt1 was cloned and identified, and it shared similar evolutionary relationships with medaka dmrt1. Tg(CMV:Osdmrt1) exhibited partial sex reversal from female to male in the F2 generation, with genetically female individuals developing testes and producing functional sperm. Additionally, the secondary sexual characteristics of the transgenic females also changed to males.

CONCLUSION

The Chinese medaka dmrt1 gene could convert females to males in medaka.

GENERAL SIGNIFICANCE

These results not only elucidate the function of Chinese medaka dmrt1, but also accumulate knowledge for studying the function of economically important fish genes in model fish by transgenic technology.

Genome-wide analysis of BMP/GDF family and DAP-seq of YY1 suggest their roles in Cynoglossus semilaevis sexual size dimorphism

Sexual size dimorphism (SSD) characterized by different body size between females and males have been reported in various animals. Gonadectomy experiments have implied important regulatory roles of the gonad in SSD. Among multiple factors from the gonad, TGF-β superfamily (especially BMP/GDF family) attracted our interest due to its pleiotropy in growth and reproduction regulations. Thus, whether BMP/GDF family members serve as crucial regulators for SSD was studied in a typically female-biased SSD flatfish named Chinese tongue sole (Cynoglossus semilaevis). Firstly, a total of 26 BMP/GDF family members were identified. The PPI network analysis showed that they may interact with ACVR2a, ACVR2b, ACVR1, BMPR2, SMAD3, BMPR1a, and other proteins. Subsequently, DAP-seq was employed to reveal the binding sites for yin yang 1 (yy1), a transcription factor involved in gonad function and cell growth partly by regulating TGF-β superfamily. The results revealed that two yy1 homologues yy1a and yy1b in C. semilaevis could regulate Hippo signaling pathway, mTOR signaling pathway, and AMPK signaling pathway. Moreover, BMP/GDF family genes including bmp2, bmp4, bmp5, gdf6a, and gdf6b were important components of Hippo pathway. In future, the crosstalk among yy1a, yy1b, and TGF-β family would provide more insight into sexual size dimorphism in C. semilaevis.

Quality of fish eggs and production of androgenetic and gynogenetic doubled haploids (DHs)

Induced development of haploid embryos (H) with only paternal (androgenesis) or maternal (gynogenesis) chromosomes requires irradiation of eggs before fertilization or activation of eggs with irradiated spermatozoa, respectively. To provide doubled haploids (DHs), androgenetic and gynogenetic haploid zygotes need to be subjected to the thermal or high hydrostatic pressure (HHP) shock to suppress the first mitotic cleavage and to double paternal or maternal haploid set of chromosomes. Androgenesis and mitotic gynogenesis (mito-gynogenesis) result in the generation of fully homozygous individuals in a single generation. DHs have been utilized in selective breeding programs, in studies concerning the phenotypic consequences of recessive alleles and to evaluate the impact of sex chromosomes on the early ontogeny. Moreover, the use of DHs for the NGS approach radically improves de novo the assembly of the genomes. However, reduced survival of the doubled haploids limits the wide application of androgenotes and gynogenotes. The high mortality of DHs may be only partly explained by the expression of recessive traits. Observed inter-clutch variation in the survival of DHs developing in eggs originating from different females make it necessary to take a closer look at the quality of the eggs used during induced androgenesis and gynogenesis. Moreover, the developmental competence of eggs that are subjected to irradiation before fertilization in order to deactivate maternal chromosomes when undergoing induced androgenesis and exposed to the physical shock after fertilization that leads to the duplication of the zygotes in both mito-gynogenesis and androgenesis may be also altered as irradiation and sublethal values of temperatures and hydrostatic pressure are considered as harmful for the cell organelles and biomolecules. Here, recently provided results concerning the morphological, biochemical, genomic, and transcriptomic characteristics of fish eggs showing high and low competence for androgenesis and mito-gynogenesis are reviewed.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

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.

The Role of the Heat Shock Cognate Protein 70 Genes in Sex Determination and Differentiation of Chinese Tongue Sole (Cynoglossus semilaevis)

Fish sex determination can be affected by environmental temperature. This process relies on temperature-sensitive proteins such as heat shock proteins (HSPs). Our previous work found that heat shock cognate proteins (HSCs) may participate in high-temperature associated sex reversal of Chinese tongue sole (Cynoglossus semilaevis). However, the role of hsc genes in responding to high temperature and affecting sex determination/differentiation remains unclear. Here, by using C. semilaevis as model, we identified hsc70 and hsc70-like. hsc70 was abundant in the gonads with a testicular-higher expression at all gonadal development stages except for 6 months post fertilization (mpf). Intriguingly, hsc70-like showed higher expression in testes from 6 mpf on. Both long-term heat treatment during the temperature-sensitive sex-determining period and short-term heat stress at the end of this period caused different expression of hsc70/hsc70-like between sexes. The dual-luciferase assay results also suggested that these genes can respond to high temperature rapidly in vitro. Heat treatment of C. semilaevis testis cells overexpressed with hsc70/hsc70-like could affect the expression of sex-related genes sox9a and cyp19a1a. Our results indicated that hsc70 and hsc70-like were key regulators linking external high-temperature signals with sex differentiation in vivo and provide a new idea for understanding the mechanism by which high temperature affects sex determination/differentiation in teleosts.

Transcriptomic Analysis Revealed Candidate Genes Involved in Pseudomale Sperm Abnormalities in Chinese Tongue Sole (Cynoglossus semilaevis)

Chinese tongue sole (Cynoglossus semilaevis) has a ZZ/ZW sex determination system, but the genotypic female (ZW) can be sex-reversed into phenotypic males, namely, pseudomales. Pseudomale fish can produce only Z-type sperm but not W sperm. However, the molecular mechanism is unclear. To screen the key genes involved in pseudomale sperm abnormalities, we analysed the transcriptomic profiles of pseudomale and male sperm. In comparison to male sperm, 592 differentially expressed genes (DEGs) were identified in pseudomale sperm, including 499 upregulated and 93 downregulated genes. KEGG analysis indicated that the FoxO signalling pathway, especially the foxo3a and foxo6-like genes, may play an important role in spermatogenesis. The DEGs were mainly distributed on sex chromosomes, with 158 downregulated genes on the Z chromosome and 41 upregulated genes on the W chromosome. A specific area (14-15 M) on the Z chromosome was identified, which enriched eight DEGs inside the ~1 M region. In addition, there were five gene alleles on the sex chromosomes, which showed the opposite transcription pattern (upregulated for the W allele, downregulated for the Z allele). This study has provided valuable data for screening candidate genes involved in the pseudomale sperm abnormality.

Identification and Expression Pattern of cyp26b1 Gene in Gonad of the Chinese Tongue Sole (Cynoglossus semilaevis)

Simple Summary

In fish, it is obvious that the asynchronous development of the gonads and sexual dimorphism limit the development of aquaculture, so the research into sex-differentiation and gonadal growth is very important. Due to the sexual reversal phenomenon (genetic females becoming phenotypic males), the Chinese tongue sole (Cynoglossus semilaevis) is a great model for investigating sex-differentiation. Herein, we report one gene involved in sex-differentiation and gonadal growth of the Chinese Tongue Sole. The gene cyp26b1 (cytochrome P450 family 26 subfamily b member 1) is a metabolizing Retinoic Acid (RA) enzyme. Since it regulates RA to control sex determination and differentiation, cyp26b1 is considered a critical part of mammals’ ovary-antagonizing and testis-determining downstream passageway of Sry (sex-determining region Y) and Sox9 (sry-box transcription factor 9). In fish, the related research is reported only on the Japanese flounder (Paralichthys olivaceus) and zebrafish (Danio rerio). In the current investigation, the identification and expression pattern of the cyp26b1 gene in the Chinese tongue sole suggested that cyp26b1 might impact sex-differentiation and gonadal development.

Abstract

As an RA-metabolizing enzyme, cyp26b1 has a substantial impact on RA-signaling pathways. The cyp26b1 gene from the Chinese tongue sole was cloned and identified in this investigation. The cyp26b1 ORF was 1536 bp in length and encoded a 512 amino acid protein. A quantitative real-time PCR (qPCR) indicated that the cyp26b1 expression is no significant sexual dimorphism in the gonads at the 80 days post-hatching (dph) stages. After 4 months post-hatching (mph), the expression of cyp26b1 showed sexual dimorphism and lower level of expression in the ovaries than in the testes. An in situ hybridization demonstrated that cyp26b1 mRNA was primarily located in the testis. Interestingly, the cyp26b1 mRNA probe was also detected in the ovaries. These results suggested that cyp26b1 participates in the sex-differentiation and gonadal development of the Chinese tongue sole.

An updated review of cold shock and cold stress in fish

Temperature is critical in regulating virtually all biological functions in fish. Low temperature stress (cold shock/stress) is an often-overlooked challenge that many fish face as a result of both natural events and anthropogenic activities. In this study, we present an updated review of the cold shock literature based on a comprehensive literature search, following an initial review on the subject by M.R. Donaldson and colleagues, published in a 2008 volume of this journal. We focus on how knowledge on cold shock and fish has evolved over the past decade, describing advances in the understanding of the generalized stress response in fish under cold stress, what metrics may be used to quantify cold stress and what knowledge gaps remain to be addressed in future research. We also describe the relevance of cold shock as it pertains to environmental managers, policymakers and industry professionals, including practical applications of cold shock. Although substantial progress has been made in addressing some of the knowledge gaps identified a decade ago, other topics (e.g., population-level effects and interactions between primary, secondary and tertiary stress responses) have received little or no attention despite their significance to fish biology and thermal stress. Approaches using combinations of primary, secondary and tertiary stress responses are crucial as a research priority to better understand the mechanisms underlying cold shock responses, from short-term physiological changes to individual- and population-level effects, thereby providing researchers with better means of quantifying cold shock in laboratory and field settings.

Identification of crucial factors involved in Cynoglossus semilaevis sexual size dimorphism by GWAS and demonstration of zbed1 regulatory network by DAP-seq

Sexual size dimorphism (SSD), whereby females and males exhibit different body sizes, are widely documented in animals. To explore crucial regulators implicated in female-biased SSD of Chinese tongue sole (Cynoglossus semilaevis), GWAS was conducted on 350 females and 59 males. Twenty SNPs and 25 genes including zbed1, nsd3, cdc45, klhl29, and smad4 with -log(p) > 7 were screened, mainly mapping to sex chromosome. The chromosome W-linked gene zbed1 attracted particular attention because it is a master key for cell proliferation. Thus, the regulatory network of zbed1 in C. semilaevis was explored by DAP-seq and 1352 peaks were discovered in the female brain. Moreover, zbed1 potentially regulated hippo signaling pathway, cell cycle, translation, and PI3k-Akt signaling pathway in C. semilaevis. These findings identify crucial SNPs and genes associated with female-biased SSD in C. semilaevis, also provide the first genome-wide survey for the zbed1 regulatory network in fish species.

Advantages, Factors, Obstacles, Potential Solutions, and Recent Advances of Fish Germ Cell Transplantation for Aquaculture-A Practical Review

Simple Summary

This review aims to provide practical information and viewpoints regarding fish germ cell transplantation for enhancing its commercial applications. We reviewed and summarized the data from more than 70 important studies and described the advantages, obstacles, recent advances, and future perspectives of fish germ cell transplantation. We concluded and proposed the critical factors for achieving better success and various options for germ cell transplantation with their pros and cons. Additionally, we discussed why this technology has not actively been utilized for commercial purposes, what barriers need to be overcome, and what potential solutions can advance its applications in aquaculture.

Abstract

Germ cell transplantation technology enables surrogate offspring production in fish. This technology has been expected to mitigate reproductive barriers, such as long generation time, limited fecundity, and complex broodstock management, enhancing seed production and productivity in aquaculture. Many studies of germ cell transplantation in various fish species have been reported over a few decades. So far, surrogate offspring production has been achieved in many commercial species. In addition, the knowledge of fish germ cell biology and the related technologies that can enhance transplantation efficiency and productivity has been developed. Nevertheless, the commercial application of this technology still seems to lag behind, indicating that the established models are neither beneficial nor cost-effective enough to attract potential commercial users of this technology. Furthermore, there are existing bottlenecks in practical aspects such as impractical shortening of generation time, shortage of donor cells with limited resources, low efficiency, and unsuccessful surrogate offspring production in some fish species. These obstacles need to be overcome through further technology developments. Thus, we thoroughly reviewed the studies on fish germ cell transplantation reported to date, focusing on the practicality, and proposed potential solutions and future perspectives.

Integration of Transcriptome and Methylome Highlights the Roles of Cell Cycle and Hippo Signaling Pathway in Flatfish Sexual Size Dimorphism

Sexual size dimorphism (SSD) is the difference in segments or body size between sexes prevalent in various species. Understanding the genetic architecture of SSD has remained a significant challenge owing to the complexity of growth mechanisms and the sexual influences among species. The Chinese tongue sole (Cynoglossus semilaevis), which exhibits a female-biased SSD and sex reversal from female to pseudomale, is an ideal model for exploring SSD mechanism at the molecular level. The present study aimed to integrate transcriptome and methylome analysis to unravel the genetic and epigenetic changes in female, male, and pseudomale C. semilaevis. The somatotropic and reproductive tissues (brain, liver, gonad, and muscle) transcriptomes were characterized by RNA-seq technology. Transcriptomic analysis unravelled numerous differentially expressed genes (DEGs) involved in cell growth and death-related pathways. The gonad and muscle methylomes were further employed for screening differentially methylated genes (DMGs). Relatively higher DNA methylation levels were observed in the male and pseudomale individuals. In detail, hypermethylation of the chromosome W was pronounced in the pseudomale group than in the female group. Furthermore, weighted gene co-expression network analysis showed that turquoise and brown modules positively and negatively correlated with the female-biased SSD, respectively. A combined analysis of the module genes and DMGs revealed the female-biased mRNA transcripts and hypomethylated levels in the upstream and downstream regions across the cell cycle-related genes. Moreover, the male and pseudomale-biased gene expression in the hippo signaling pathway were positively correlated with their hypermethylation levels in the gene body. These findings implied that the activation of the cell cycle and the inhibition of the hippo signaling pathway were implicated in C. semilaevis female-biased SSD. In addition, the dynamic expression pattern of the epigenetic regulatory factors, including dnmt1, dnmt3a, dnmt3b, and uhrf1, among the different sexes correspond with their distinct DNA methylation levels. Herein, we provide valuable clues for understanding female-biased SSD in C. semilaevis.

A Z-Linked E3 Ubiquitin Ligase Cs-rchy1 Is Involved in Gametogenesis in Chinese Tongue Sole, Cynoglossus semilaevis

Simple Summary

The sexual growth dimorphism prevails in animals and this phenomenon is even more obvious in marine fish, so understanding the mechanism of gonadal development and gametogenesis is of great importance for sex control, thus increased productivity in aquaculture. In mammal, ubiquitin ligase plays a versatile role in gonadal development and spermatogenesis, whereas its function in fish is little reported. Using Cynoglossus semilaevis (one-year-old female individual usually grows 2–4 times bigger than male) as the fish model, a Z-chromosome linked ubiquitin ligase neurl3 was previously identified and characterized, which suggested its involvement in spermatogenesis. However, in this study, characterization of another Z-chromosome linked ubiquitin ligase Cs-rchy1 suggested it might function both in spermatogenesis and oogenesis, as well as the potential role in growth. These data may provide the genetic resource for gene editing or marker exploration in future.

Abstract

Ubiquitin ligase (E3) plays a versatile role in gonadal development and spermatogenesis in mammals, while its function in fish is little reported. In this study, a Z-chromosome linked ubiquitin ligase rchy1 in C. semilaevis (Cs-rchy1) was cloned and characterized. The full-length cDNA was composed of 1962 bp, including 551 bp 5′UTR, 736 bp 3′UTR, and 675 bp ORF encoding a 224-amino-acid (aa) protein. Cs-rchy1 was examined among seven different tissues and found to be predominantly expressed in gonads. In testis, Cs-rchy1 could be detected from 40 days post hatching (dph) until 3 years post hatching (yph), but there was a significant increase at 6 months post hatching (mph). In comparison, the expression levels in ovary were rather stable among different developmental stages. In situ hybridization showed that Cs-rchy1 was mainly localized in germ cells, that is, spermatid and spermatozoa in testis and stage I, II and III oocytes in ovary. In vitro RNA interference found that Cs-rchy1 knockdown resulted in the decline of sox9 and igf1 in ovarian cell line and down-regulation of cyp19a in the testicular cell line. These data suggested that Cs-rchy1 might participate in gonadal differentiation and gametogenesis, via regulating steroid hormone synthesis.

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.

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.

Cloning and functional analysis of c/ebpα as negative regulator of dmrt1 in Chinese tongue sole (Cynoglossus semilaevis)

c/ebpα is a member of the C/EBP family of transcription factors, which are involved in cell growth and differentiation and have a conserved basic leucine zipper (bZIP) domain. However, little is known about its function in sex determination and differentiation. In the present study, c/ebpα was cloned from the gonads of Chinese tongue sole (Cynoglossus semilaevis). The full-length cDNA of c/ebpα was 1583 bp, with a 198-bp 5' UTR, a 446-bp 3' UTR, and a 939-bp open reading frame encoding a 312-amino acid peptide. qRT-PCR revealed that c/ebpα was predominantly expressed in undifferentiated gonads of male C. semilaevis at 30 dpf and 60 dpf and peaked at 60 dpf. Expression levels of c/ebpα in the testis were constantly higher than those in ovaries at all developmental stages. Moreover, a dual-luciferase assay revealed that c/ebpα could negatively regulate the male-determining gene dmrt1 in vitro. These results provide fundamental information indicating that C. semilaevis c/ebpa might be involved in early gonadal differentiation and functions as a negative regulator of dmrt1 by repressing its transcription.

Comparative transcriptome analyses and identification of candidate genes involved in vertebral abnormality of bighead carp Hypophthalmichthys nobilis

Body deformity occurs both in wild and farmed fishes, which is one of the most challenging problems for aquaculture industry. In most cases, such body deformities are linked to skeletal deformities. Currently, very limited information is available on skeletal deformities of farmed fish species which may be caused by genetic factor. In this study, we performed muscle and vertebra transcriptome analyses in body deformity and normality of bighead carp Hypophthalmichthys nobilis (from one meiotic gynogenesis family) using RNA-Seq. A total of 43,923 and 44,416 unigenes were predicted in muscles and vertebrae, respectively. Based on these data, we further explored the gene expression profiles in gynogenetic normal and abnormal bighead carp. No differentially expressed gene (DEG) was found in transcriptome data of muscles. Totally, 20 key DEGs were identified in transcriptome data of vertebrae, such as low density lipoprotein-related protein 2 (lrp2), bone morphogenetic protein 2B (bmp2b) and collagen alpha-1(IV) (col4a1). 12 potential pathways were also identified in vertebra transcriptome data, which were mainly involved in development, growth, cytoskeleton and energy metabolism, such as MAPK signaling pathway, regulation of actin cytoskeleton and TGF-beta signaling pathway. Results of this study will be informative for the understanding of genetic mechanisms for body shape formation and also provide potential candidate genes for selection program involved in body shape and skeletal development in H. nobilis.

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.

Further evidence for paternal DNA transmission in gynogenetic grass carp

Gynogenesis is an important breeding method in aquaculture and has been widely applied to many fish species. If gynogenetic progenies are to inherit paternal partial genomic DNA, this will increase genetic variation and will provide a useful outcome for breeding. In this study, we investigated the genetic variation in homeobox (Hox) gene clusters (HoxA4a, HoxA9a, HoxA11b, HoxB1b, HoxC4a, HoxC6b, and HoxD10a) among koi carp (Cyprinus carpio haematopterus, KOC; the stimulation sperm source), grass carp (Ctenopharyngodon idellus), and gynogenetic grass carp (GGC). We found paternal DNA (a special DNA fragment and HoxC6b) derived from KOC and a recombinant gene belonging to HoxC6b in GGC. We are the first to report the recombinant HoxC6b in GGC. Our study provides further evidence for paternal DNA transmission to gynogenetic progenies, which is a finding with great significance for the genetic breeding of fish.

Large scale SNP unearthing and genetic architecture analysis in sea-captured and cultured populations of Cynoglossus semilaevis

Knowledge on population structure and genetic diversity is a focal point for association mapping studies and genomic selection. Genotyping by sequencing (GBS) represents an innovative method for large scale SNP detection and genotyping of genetic resources. Here we used the GBS approach for the genome-wide identification of SNPs in a collection of Cynoglossus semilaevis and for the assessment of the level of genetic diversity in C. semilaevis genotypes. GBS analysis generated a total of 55.12 Gb high-quality sequence data, with an average of 0.63 Gb per sample. The total number of SNP markers was 563, 109. In order to explore the genetic diversity of C. semilaevis and to select a minimal core set representing most of the total genetic variation with minimum redundancy, C. semilaevis sequences were analyzed using high quality SNPs. Based on hierarchical clustering, it was possible to divide the collection into 2 clusters. The marine fishing populations were clustered and clearly separated from the cultured populations, and the cultured populations from Hebei was also distinct from the other two local populations. These analyses showed that genotypes were clustered based on species-related features. Differential significant SNPs were also captured and validated by GBS and SNaPshot, with linkage disequilibrium and haplotype analysis, seven SNPs have been confirmed to have obvious differentiation in two populations, which may be used as the characteristic evaluation sites of sea-captured and cultured Cynoglossus semilaevis populations. And SNP markers and information on population structure developed in this study will undoubtedly support genome-wide association mapping studies and marker-assisted selection programs. These differential SNPs could be also employed as the characteristic evaluation sites of sea-captured and cultured Cynoglossus semilaevis populations in future.

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.

Integrative analysis reveals pathways associated with sex reversal in Cynoglossus semilaevis

Sex reversal is a complex biological phenomenon exhibited by Cynoglossus semilaevis. Some genetic females may irreversibly convert to pseudomales, thus increasing aquaculture costs because males grow much more slowly than females. In this study, an integrative analysis of transcriptome and proteome was performed to compare differences in gene and protein expression in females and pseudomales after gonad differentiation in C. semilaevis. Based on RNA-Seq results, 1893 genes showed differences in expression at the transcript level between females and pseudomales. Of these differentially expressed genes (DEGs), zona pellucida sperm-binding protein 4-like (LOC103393374 , ZP4), zona pellucida sperm-binding protein 4-like (LOC103396071, ZP4) and forkhead box L2 (foxl2) were highly expressed in females and doublesex and mab-3 related transcription factor 1(dmrt1) and doublesex and mab-3 related transcription factor 3 (dmrt3) were highly expressed in pseudomales. GO enrichment analysis results indicate that wnt signaling pathways and oocyte maturation are two terms enriched in female. At the protein level, Tandem Mass Tags analysis revealed that 324 proteins differed in their relative abundance between pseudomales and females. KEGG analysis found that pseudo-highly expressed proteins were enriched in the ubiquitin mediated proteolysis pathway. For integrative analysis, the Spearman correlation coefficient between the transcriptome and proteome was 0.59. Among 52 related genes, 46 DEGs (88%) were well matched in their levels of change in protein abundance. These findings reveal major active pathways in female and pseudomale gonads after sex reversal and provide new insights into molecular mechanisms associated with sex reversal regulatory network.

Genome-wide identification and characterization of toll-like receptor genes in spotted sea bass (Lateolabrax maculatus) and their involvement in the host immune response to Vibrio harveyi infection

Toll-like receptor (TLR) genes are the earliest reported pathogen recognition receptors (PRRs) and have been extensively studied. These genes play pivotal roles in the innate immune defense against pathogen invasion. In this study, a total of 16 tlr genes were identified and characterized in spotted sea bass (Lateolabrax maculatus). The tlr genes of spotted sea bass were classified into five subfamilies (tlr1-subfamily, tlr3-subfamily, tlr5-subfamily, tlr7-subfamily, and tlr11-subfamily) according to the phylogenetic analysis, and their annotations were confirmed by a syntenic analysis. The protein domain analysis indicated that most tlr genes had the following three major TLR protein domains: a leucine-rich repeat (LRR) domain, a transmembrane region (TM) and a Toll/interleukin-1 receptor (TIR) domain. The tlr genes in spotted sea bass were distributed in 11 of 24 chromosomes. The mRNA expression levels of 16 tlr genes in response to Vibrio harveyi infection were quantified in the head kidney. Most genes were downregulated following V. harveyi infection, while only 5 tlr genes, including tlr1-1, tlr1-2, tlr2-2, tlr5, and tlr7, were significantly upregulated. Collectively, these results help elucidate the crucial roles of tlr genes in the immune response of spotted sea bass and may supply valuable genomic resources for future studies investigating fish disease management.

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

Dmrt2 is a member of the dmrt gene family with a conserved zinc finger-like DNA-binding motif (DM domain). In the present study, CS-dmrt2 was cloned from the gonads of Chinese tongue sole (Cynoglossus semilaevis). The full-length cDNA of CS-dmrt2 is 2834 bp in length, with a 251 bp 5'-untranslated region (UTR), a 1086 bp 3'-UTR and a 1503 bp open reading frame (ORF) that encodes a 501-amino-acid peptide. qPCR revealed that CS-dmrt2 was mainly expressed in C. semilaevis testes. In situ hybridization (ISH) showed CS-dmrt2 expression throughout early gonadal development (36 days after hatching (dah) and 86 dah), but the expression was higher in male gonads than in female gonads. CS-dmrt2 mRNA was highly expressed in male germ cells. Comparison of methylation levels between females and males demonstrated hypo-methylated levels of the CS-dmrt2 promoter in the male gonads, which is consistent with the high mRNA expression. These results suggest that the CS-dmrt2 gene may play a functional role in gonadal differentiation/development and germ cell maturation in the testis.

The autosomal Gsdf gene plays a role in male gonad development in Chinese tongue sole (Cynoglossus semilaevis)

Gsdf is a key gene for testicular differentiation in teleost. However, little is known about the function of Gsdf in Chinese tongue sole (Cynoglossus semilaevis). In this study, we obtained the full-length Gsdf gene (CS-Gsdf), and functional characterization revealed its potential participation during germ cell differentiation in testes. CS-Gsdf transcription was predominantly detected in gonads, while the levels in testes were significantly higher than those in ovaries. During the different developmental stages in male gonads, the mRNA level was significantly upregulated at 86 dph, and a peak appeared at 120 dph; then, the level decreased at 1 and 2 yph. In situ hybridization revealed that CS-Gsdf mRNA was mainly localized in the Sertoli cells, spermatogonia, and spermatids in mature testes. After CS-Gsdf knockdown in the male testes cell line by RNA interference, a series of sex-related genes was influenced, including several sex differentiation genes, CS-Wnt4a, CS-Cyp19a1a and CS-Star. Based on these data, we speculated that CS-Gsdf may play a positive role in germ differentiation and proliferation via influencing genes related to sex differentiation.

Characterization and expression pattern of r-spondin1 in Cynoglossus semilaevis

r-spondin1 (rspo1) encodes a secreted protein that is involved in the determination and differentiation of the mammalian ovary. However, little information is yet available for teleosts. Here, we identified a homologue of rspo1 in Cynoglossus semilaevis. The full-length cDNA of rspo1 had a length of 2,703 bp with an open reading frame of 834 bp, encoding a protein with a length of 277 amino-acids. rspo1 expression was detected via qRT-PCR in various tissues, and significant sexually dimorphic expression was observed in the gonads. Furthermore, ISH located rspo1 in germ cells such as spermatogonia, spermatocytes, spermatids, spermatozoa, and oocytes, as well as in somatic cells of the gonads. Following knockdown of rspo1 in an ovarian cell line, the expressions of wnt4a, β-catenin, foxl2, and StAR were highly affected; wnt4a and β-catenin were significantly downregulated, whereas foxl2 and StAR were significantly upregulated. In summary, these data suggest that rspo1 may be involved in the regulation of ovarian development and differentiation through a conserved pathway, while the function of the gene in the testis remains elusive.

Environmental epigenetics in zebrafish

It is widely accepted that the epigenome can act as the link between environmental cues, both external and internal, to the organism and phenotype by converting the environmental stimuli to phenotypic responses through changes in gene transcription outcomes. Environmental stress endured by individual organisms can also enforce epigenetic variations in offspring that had never experienced it directly, which is termed transgenerational inheritance. To date, research in the environmental epigenetics discipline has used a wide range of both model and non-model organisms to elucidate the various epigenetic mechanisms underlying the adaptive response to environmental stimuli. In this review, we discuss the advantages of the zebrafish model for studying how environmental toxicant exposures affect the regulation of epigenetic processes, especially DNA methylation, which is the best-studied epigenetic mechanism. We include several very recent studies describing the state-of-the-art knowledge on this topic in zebrafish, together with key concepts in the function of DNA methylation during vertebrate embryogenesis.

Genome editing reveals dmrt1 as an essential male sex-determining gene in Chinese tongue sole (Cynoglossus semilaevis)

Chinese tongue sole is a marine fish with ZW sex determination. Genome sequencing suggested that the Z-linked dmrt1 is a putative male determination gene, but direct genetic evidence is still lacking. Here we show that TALEN of dmrt1 efficiently induced mutations of this gene. The ZZ dmrt1 mutant fish developed ovary-like testis, and the spermatogenesis was disrupted. The female-related genes foxl2 and cyp19a1a were significantly increased in the gonad of the ZZ dmrt1 mutant. Conversely, the male-related genes Sox9a and Amh were significantly decreased. The dmrt1 deficient ZZ fish grew much faster than ZZ male control. Notably, we obtained an intersex ZW fish with a testis on one side and an ovary on the other side. This fish was chimeric for a dmrt1 mutation in the ovary, and wild-type dmrt1 in the testis. Our data provide the first functional evidence that dmrt1 is a male determining gene in tongue sole.

Single Locus Maintains Large Variation of Sex Reversal in Half-Smooth Tongue Sole (Cynoglossus semilaevis)

Sex determination is a fundamental biological process for individual sex development and population sex ratios. However, for some species, the primary sex might be altered during development, and individuals can develop into the opposite sex. Sex reversal may happen in insects, reptiles, amphibians, and fishes. In half-smooth tongue sole (Cynoglossus semilaevis), some genetically female fish irreversibly reverse to pseudomales, resulting in higher costs in aquaculture owing to a lower growth rate of male fish during a 2-yr growth period. Here, we identified a locus with large controlling effect on sex reversal in the half-smooth tongue sole through genome-wide association study with high-density single nucleotide polymorphisms (SNPs). This SNP is located at the third intron of the F-box and leucine rich repeat protein 17 (FBXL17) gene on the Z chromosome, and it has two alleles, A and T. Genetic females with Z^AW genotypes will never reverse into phenotypic males, but those with Z^TW genotypes can sometimes undergo sex reversal. This SNP explains 82.7% of the genetic variation, or 58.4% of the phenotypic variation. Based on our results, a reproductive management program could be developed to improve the phenotypic female ratio in aquaculture, and elucidate the mechanism of sex reversal in half-smooth tongue sole. We expect that these findings will have a substantial impact on the population management in many harvested species where sex reversal occurs.

Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production

Flatfish have a high market acceptance thus representing a profitable aquaculture production. The main farmed species is the turbot (Scophthalmus maximus) followed by Japanese flounder (Paralichthys olivaceous) and tongue sole (Cynoglossus semilaevis), but other species like Atlantic halibut (Hippoglossus hippoglossus), Senegalese sole (Solea senegalensis) and common sole (Solea solea) also register an important production and are very promising for farming. Important genomic resources are available for most of these species including whole genome sequencing projects, genetic maps and transcriptomes. In this work, we integrate all available genomic information of these species within a common framework, taking as reference the whole assembled genomes of turbot and tongue sole (>210× coverage). New insights related to the genetic basis of productive traits and new data useful to understand the evolutionary origin and diversification of this group were obtained. Despite a general 1:1 chromosome syntenic relationship between species, the comparison of turbot and tongue sole genomes showed huge intrachromosomic reorganizations. The integration of available mapping information supported specific chromosome fusions along flatfish evolution and facilitated the comparison between species of previously reported genetic associations for productive traits. When comparing transcriptomic resources of the six species, a common set of ~2500 othologues and ~150 common miRNAs were identified, and specific sets of putative missing genes were detected in flatfish transcriptomes, likely reflecting their evolutionary diversification.

Two Figla homologues have disparate functions during sex differentiation in half-smooth tongue sole (Cynoglossus semilaevis)

Figla is a germ-cell-specific transcription factor associated with ovary development and differentiation. In vertebrates, one transcriptional form of Figla is commonly found. However, besides the common form of this gene (named Figla_tv1), a new transcriptional form (named Figla_tv2) was identified in half-smooth tongue sole (Cynoglossus semilaevis). The full-length cDNA of Figla_tv1 was 1057 bp long with a 591-bp open reading frame encoding a predicted 196 amino acid protein, while Figla_tv2 encoded a 125 amino acid protein. Figla_tv1 and Figla_tv2 expression in various tissues was detected by qRT-PCR. Figla_tv1 was expressed mainly in ovary, skin and liver, while Figla_tv2 was expressed in all examined tissues. In the gonads, Figla_tv1 was expressed in ovary, while Figla_tv2 was predominately expressed in testis of pseudomales. Further, in situ hybridization located Figla_tv1 only in oocytes and Figla_tv2 mainly in germ cells of pseudomale testis. After knocking down Figla_tv2 in a pseudomale testis cell line, the expression of two steroid hormone-encoding genes, StAR and P450scc, was significantly up-regulated (P < 0.05). Our findings suggest that Figla_tv1 has a conserved function in folliculogenesis, as in other vertebrates, and that Figla_tv2 may have a role in the spermatogenesis of pseudomales by regulating the synthesis and metabolism of steroid hormones.

Integrated gene mapping and synteny studies give insights into the evolution of a sex proto-chromosome in Solea senegalensis

The evolution of genes related to sex and reproduction in fish shows high plasticity and, to date, the sex determination system has only been identified in a few species. Solea senegalensis has 42 chromosomes and an XX/XY chromosome system for sex determination, while related species show the ZZ/ZW system. Next-generation sequencing (NGS), multi-color fluorescence in situ hybridization (mFISH) techniques, and bioinformatics analysis have been carried out, with the objective of revealing new information about sex determination and reproduction in S. senegalensis. To that end, several bacterial artificial chromosome (BAC) clones that contain candidate genes involved in such processes (dmrt1, dmrt2, dmrt3, dmrt4, sox3, sox6, sox8, sox9, lh, cyp19a1a, amh, vasa, aqp3, and nanos3) were analyzed and compared with the same region in other related species. Synteny studies showed that the co-localization of dmrt1-dmrt2-drmt3 in the largest metacentric chromosome of S. senegalensis is coincident with that found in the Z chromosome of Cynoglossus semilaevis, which would potentially make this a sex proto-chromosome. Phylogenetic studies show the close proximity of S. senegalensis to Oryzias latipes, a species with an XX/XY system and a sex master gene. Comparative mapping provides evidence of the preferential association of these candidate genes in particular chromosome pairs. By using the NGS and mFISH techniques, it has been possible to obtain an integrated genetic map, which shows that 15 out of 21 chromosome pairs of S. senegalensis have at least one BAC clone. This result is important for distinguishing those chromosome pairs of S. senegalensis that are similar in shape and size. The mFISH analysis shows the following co-localizations in the same chromosomes: dmrt1-dmrt2-dmrt3, dmrt4-sox9-thrb, aqp3-sox8, cyp19a1a-fshb, igsf9b-sox3, and lysg-sox6.

Establishment and characterization of a gonad cell line from half-smooth tongue sole Cynoglossus semilaevis pseudomale

A new cell line was established from half-smooth tongue sole Cynoglossus semilaevis pseudomale gonad (CSPMG). Primary culture was initiated from gonad tissues pieces, and the CSPMG cells were cultured at 24 °C in Dulbecco's modified Eagle medium/F12 medium (1:1) (pH7.0), supplemented with 20 % fetal bovine serum, basic fibroblast growth factor, epidermal growth factor, insulin-like growth factor-I, 2-mercaptoethanol, penicillin and streptomycin. The cultured CSPMG cells, in fibroblast shape, proliferated to 100 % confluency 10 days later and had been subcultured to passage 109. Chromosome analyses indicated that the CSPMG cells exhibited chromosomal aneuploidy with a modal chromosome number of 42, which displayed the normal diploid karyotype of half-smooth tongue sole (2n = 42t, NF = 42). Reverse transcription polymerase chain reaction revealed CSPMG cells could express gonad somatic cell functional genes Sox9a, Wt1a and weakly germ cell marker gene Vasa, but not male specific gene Dmrt1. Transfection experiment demonstrated that CSPMG cells transfected with pEGFP-N3 plasmid and small RNA could express green and red fluorescence signals with high transfection efficiency. In conclusion, a continuous CSPMG cell line has been established successfully. The cell line might serve as a valuable tool for studies on the mechanism of sex determination, sex reversal and gonad development in flatfish.

Sperm of the giant grouper: cryopreservation, physiological and morphological analysis and application in hybridizations with red-spotted grouper

In order to develop excellent germplasm resources for giant grouper (Epinephelus lanceolatus), cryopreservation of giant grouper sperm was examined in the present study. Firstly, 13 kinds of sperm dilution (ELS1-3, EM1-2, TS-2, MPRS, ELRS0-6) were prepared with physiological salt, sucrose, glucose and fetal bovine serum. The physiological parameters of ELRS3 (ratio of fast motion, ratio of slow motion, time of fast motion, time of slow motion, lifespan and motility) and ELS3 (sperm ratio of slow motion, time of slow motion and motility) were significantly higher than those of the other dilutions (P < 0.05). Secondly, after adding 15% DMSO and 10% FBS to ELRS3 and ELS3, most physiological parameters of frozen sperm were also significantly higher than the other gradients (P < 0.05), and sperm motility was as high as 63.68 ± 4.16% to74.75 ± 12.71% (fresh sperm motility, 80.70 ± 1.37% to 80.71 ± 1.49%). Mixed with the above dilutions, a final volume of 105 ml semen was cryopreserved. Finally, the sperm of giant grouper cryopreserved with cryoprotectants (ELRS3 + 15% DMSO + 10% FBS) was used for electron-microscopic observation and crossbreeding with red-spotted groupers (Epinephelus akaara). The electron-microscopic observation revealed that part of the frozen-thawed sperm was cryodamaged, e.g., flagellum fracturing and mitochondria falling out, while the ultrastructure of sperm membrane, mitochondria and flagellum remained intact. Also, the fertilization and hatchability rates of giant grouper frozen sperm and red-spotted grouper eggs were as high as 94.56% and 75.56%, respectively. Thus, a technique for cryopreservation of giant grouper sperm was successfully developed and applied to crossbreeding with red-spotted grouper eggs.

Cloning and characterization of tesk1, a novel spermatogenesis-related gene, in the tongue sole (Cynoglossus semilaevis)

Testis-specific protein kinase 1 (Tesk1) is a serine/threonine kinase with unique structural features. In the present study, we cloned and characterized the tesk1 gene of tongue sole, Cynoglossus semilaevis. The full-length tesk1 cDNA consists of 1,672 nucleotides, encoding a 331 amino acid polypeptide with a characteristic structure composed of an N-terminal kinase domain and a C-terminal proline-rich domain. The tesk1 genomic sequence contains eight exons and seven introns. Real-time quantitative PCR revealed that tesk1 mRNA is expressed predominantly in the testis, though the level of expression varied throughout development. We used in situ hybridization to show that tesk1 mRNA is expressed in the spermatids of males and pseudo-males, but not in triploid males. Our results suggest that tongue sole Tesk1 may play a role in spermatogenesis.

Chromosomal manipulation in Senegalese sole (Solea senegalensis Kaup, 1858): induction of triploidy and gynogenesis

In this study we have developed protocols for induced triploidy and gynogenesis of Senegalese sole (Solea senegalensis), a promising flatfish species for marine aquaculture, in order to: 1) identify the sex-determination mechanism; and 2) to improve its production by generating a) sterile fish, avoiding problems related with sexual maturation, and b) all-female stocks, of higher growth rate. Triploidy was induced by means of a cold shock. Gynogenesis was induced by activating eggs with UV-irradiated sperm, and to prompt diploid gynogenesis, a cold-shock step was also used. Ploidy of putative triploid larvae and gynogenetic embryos were determined by means of karyotyping and microsatellite analysis. Haploid gynogenetic embryos showed the typical "haploid syndrome". As expected, triploid and gynogenetic groups showed lower fertilization, hatching, and survival rates than in the diploid control group. Survival rate, calculated 49 days after hatching, for haploid and diploid gynogenetic groups was similar to those observed in other fish species (0% and 62.5%, respectively), whereas triploids showed worse values (45%). Sex was determined macroscopically and by histological procedures, revealing that all the diploid gynogenetic individuals were females. In conclusion, we have successfully applied chromosomal-manipulation techniques in the flatfish species Senegalese sole in order to produce triploid, haploid, and diploid gynogenetic progenies.

Molecular method for sex identification of half-smooth tongue sole (Cynoglossus semilaevis) using a novel sex-linked microsatellite marker

Half-smooth tongue sole (Cynoglossus semilaevis) is one of the most important flatfish species for aquaculture in China. To produce a monosex population, we attempted to develop a marker-assisted sex control technique in this sexually size dimorphic fish. In this study, we identified a co-dominant sex-linked marker (i.e., CyseSLM) by screening genomic microsatellites and further developed a novel molecular method for sex identification in the tongue sole. CyseSLM has a sequence similarity of 73%–75% with stickleback, medaka, Fugu and Tetraodon. At this locus, two alleles (i.e., A244 and A234) were amplified from 119 tongue sole individuals with primer pairs CyseSLM-F1 and CyseSLM-R. Allele A244 was present in all individuals, while allele A234 (female-associated allele, FAA) was mostly present in females with exceptions in four male individuals. Compared with the sequence of A244, A234 has a 10-bp deletion and 28 SNPs. A specific primer (CyseSLM-F2) was then designed based on the A234 sequence, which amplified a 204 bp fragment in all females and four males with primer CyseSLM-R. A time-efficient multiplex PCR program was developed using primers CyseSLM-F2, CyseSLM-R and the newly designed primer CyseSLM-F3. The multiplex PCR products with co-dominant pattern could be detected by agarose gel electrophoresis, which accurately identified the genetic sex of the tongue sole. Therefore, we have developed a rapid and reliable method for sex identification in tongue sole with a newly identified sex-linked microsatellite marker.

A first generation BAC-based physical map of the half-smooth tongue sole (Cynoglossus semilaevis) genome

Background

Half-smooth tongue sole (Cynoglossus semilaevis Günther) has been exploited as a commercially important cultured marine flatfish, and female grows 2–3 times faster than male. Genetic studies, especially on the chromosomal sex-determining system of this species, have been carried out in the last decade. Although the genome of half-smooth tongue sole was relatively small (626.9 Mb), there are still some difficulties in the high-quality assembly of the next generation genome sequencing reads without the assistance of a physical map, especially for the W chromosome of this fish due to abundance of repetitive sequences. The objective of this study is to construct a bacterial artificial chromosome (BAC)-based physical map for half-smooth tongue sole with the method of high information content fingerprinting (HICF).

Results

A physical map of half-smooth tongue sole was constructed with 30, 294 valid fingerprints (7.5 × genome coverage) with a tolerance of 4 and an initial cutoff of 1e-60. A total of 29,709 clones were assembled into 1,485 contigs with an average length of 539 kb and a N50 length of 664 kb. There were 394 contigs longer than the N50 length, and these contigs will be a useful resource for future integration with linkage map and whole genome sequence assembly. The estimated physical length of the assembled contigs was 797 Mb, representing approximately 1.27 coverage of the half-smooth tongue sole genome. The largest contig contained 410 BAC clones with a physical length of 3.48 Mb. Almost all of the 676 BAC clones (99.9%) in the 21 randomly selected contigs were positively validated by PCR assays, thereby confirming the reliability of the assembly.

Conclusions

A first generation BAC-based physical map of half-smooth tongue sole was constructed with high reliability. The map will promote genetic improvement programs of this fish, especially integration of physical and genetic maps, fine-mappings of important gene and/or QTL, comparative and evolutionary genomics studies, as well as whole genome sequence assembly.

Epigenetic modification and inheritance in sexual reversal of fish

Environmental sex determination (ESD) occurs in divergent, phylogenetically unrelated taxa, and in some species, co-occurs with genetic sex determination (GSD) mechanisms. Although epigenetic regulation in response to environmental effects has long been proposed to be associated with ESD, a systemic analysis on epigenetic regulation of ESD is still lacking. Using half-smooth tongue sole (Cynoglossus semilaevis) as a model—a marine fish that has both ZW chromosomal GSD and temperature-dependent ESD—we investigated the role of DNA methylation in transition from GSD to ESD. Comparative analysis of the gonadal DNA methylomes of pseudomale, female, and normal male fish revealed that genes in the sex determination pathways are the major targets of substantial methylation modification during sexual reversal. Methylation modification in pseudomales is globally inherited in their ZW offspring, which can naturally develop into pseudomales without temperature incubation. Transcriptome analysis revealed that dosage compensation occurs in a restricted, methylated cytosine enriched Z chromosomal region in pseudomale testes, achieving equal expression level in normal male testes. In contrast, female-specific W chromosomal genes are suppressed in pseudomales by methylation regulation. We conclude that epigenetic regulation plays multiple crucial roles in sexual reversal of tongue sole fish. We also offer the first clues on the mechanisms behind gene dosage balancing in an organism that undergoes sexual reversal. Finally, we suggest a causal link between the bias sex chromosome assortment in the offspring of a pseudomale family and the transgenerational epigenetic inheritance of sexual reversal in tongue sole fish.

Cloning, genomic structure and expression analysis of ubc9 in the course of development in the half-smooth tongue sole (Cynoglossus semilaevis)

The small ubiquitin-like modifier (SUMO) pathway is an essential biological process in eukaryote, and Ubc9 is an important E2 conjugating enzyme (UBE2) for SUMO pathway and plays a critical role in cellular differentiation, development and sex modification in various species. However, the relationship between Ubc9 and sex modification and development in fish remains elusive. To elucidate the impact of Ubc9 on sex modification and development, the full length of the cDNA and genomic sequence was cloned from half-smooth tongue sole, Cynoglossus semilaevis. Real-time quantitative RT-PCR demonstrated that ubc9 was ubiquitously expressed in different tissues, and the expression levels varied in the different stages of embryonic and gonadal development. In addition, the expression level was significantly higher in the temperature-treated females than the normal females and males. Moreover, the PET-32-Ubc9 plasmid was constructed and the recombinant protein was expressed in Escherichia coli. Follistatin gene expression was initially up-regulated and FSE genes (cyp19a1a, ctnnb1, foxl2) were initially down-regulated after the injection of Ubc9 protein, prior to 96 h eventually recovered to normal levels. Taken together, the results show that Ubc9 is involved in embryogenesis, gametogenesis and sex modification, and exerts an effect on gene expression.