Mapping of DNA sex-specific markers and genes related to sex differentiation in turbot (Scophthalmus maximus)

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.

Whole-Transcriptome Analysis Identifies Gender Dimorphic Expressions of Mrnas and Non-Coding Rnas in Chinese Soft-Shell Turtle (Pelodiscus sinensis)

In aquaculture, the Chinese soft-shelled turtle (Pelodiscus sinensis) is an economically important species with remarkable gender dimorphism in its growth patterns. However, the underlying molecular mechanisms of this phenomenon have not been elucidated well. Here, we conducted a whole-transcriptome analysis of the female and male gonads of P. sinensis. Overall, 7833 DE mRNAs, 619 DE lncRNAs, 231 DE circRNAs, and 520 DE miRNAs were identified. Some "star genes" associated with sex differentiation containing dmrt1, sox9, and foxl2 were identified. Additionally, some potential genes linked to sex differentiation, such as bmp2, ran, and sox3, were also isolated in P. sinensis. Functional analysis showed that the DE miRNAs and DE ncRNAs were enriched in the pathways related to sex differentiation, including ovarian steroidogenesis, the hippo signaling pathway, and the calcium signaling pathway. Remarkably, a lncRNA/circRNA-miRNA-mRNA interaction network was constructed, containing the key genes associated with sex differentiation, including fgf9, foxl3, and dmrta2. Collectively, we constructed a gender dimorphism profile of the female and male gonads of P. sinensis, profoundly contributing to the exploration of the major genes and potential ncRNAs involved in the sex differentiation of P. sinensis. More importantly, we highlighted the potential functions of ncRNAs for gene regulation during sex differentiation in P. sinensis as well as in other turtles.

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.

QTL Fine Mapping for Sex Determination Region in Bighead Carp (Hypophthalmichthys nobilis) and Comparison with Silver Carp (Hypophthalmichthys molitrix)

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix) are genetically close aquaculture fish in the Cyprinidae, which have been confirmed to hold XX/XY sex determination. However, genomic locations of potential sex-related loci in these two fishes are still unknown. In this study, a high-resolution genetic linkage map was constructed by using 2976 SNP and 924 microsatellite markers in a F1 full-sib family of bighead carp, the length of which spanned 2022.34 cM with an average inter-marker distance of 0.52 cM. Comparative genomics revealed a high level of genomic synteny between bighead carp and zebrafish as well as grass carp. QTL fine mapping for sex trait was performed based on this linkage map of bighead carp and an unpublished linkage map of silver carp. A map distance of 3.863 cM (69.787-73.650 cM) on LG19 of bighead carp and 4.705 cM (79.096-83.801 cM) on LG21 of silver carp was significantly associated with sex phenotypes, and these two LGs are homologous between two fish species. Fourteen markers harboring in these regions were in strong linkage disequilibrium with the sex phenotype variance explained (PVE) varying from 89 to 100%. Two common markers were mapped on the QTL regions of bighead carp and silver carp, suggesting that these two carp species may have similar genetic bases for sex determination. Eleven potentially sex-related genes were identified within or near the sex QTL markers in two species. This study provided insights into elucidating mechanisms and evolution of sex determination in cyprinid fishes.

Characterization of gonadal transcriptomes from the turbot (Scophthalmus maximus)

The mechanisms underlying sexual reproduction and sex ratio determination remains unclear in turbot, a flatfish of great commercial value. And there is limited information in the turbot database regarding genes related to the reproductive system. Here, we conducted high-throughput transcriptome profiling of turbot gonad tissues to better understand their reproductive functions and to supply essential gene sequence information for marker-assisted selection programs in the turbot industry. In this study, two gonad libraries representing sex differences in Scophthalmus maximus yielded 453 818 high-quality reads that were assembled into 24 611 contigs and 33 713 singletons by using 454 pyrosequencing, 13 936 contigs and singletons (CS) of which were annotated using BLASTx. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses revealed that various biological functions and processes were associated with many of the annotated CS. Expression analyses showed that 510 genes were differentially expressed in males versus females; 80% of these genes were annotated. In addition, 6484 and 6036 single nucleotide polymorphisms (SNPs) were identified in male and female libraries, respectively. This transcriptome resource will serve as the foundation for cDNA or SNP microarray construction, gene expression characterization, and sex-specific linkage mapping in turbot.

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.

Comprehensive transcriptomic analysis of the process of gonadal sex differentiation in the turbot (Scophthalmus maximus)

The turbot is a flatfish with a ZW/ZZ sex determination system but with a still unknown sex determining gene(s), and with a marked sexual growth dimorphism in favor of females. To better understand sexual development in turbot we sampled young turbot encompassing the whole process of gonadal differentiation and conducted a comprehensive transcriptomic study on its sex differentiation using a validated custom oligomicroarray. Also, the expression profiles of 18 canonical reproduction-related genes were studied along gonad development. The expression levels of gonadal aromatase cyp19a1a alone at three months of age allowed the accurate and early identification of sex before the first signs of histological differentiation. A total of 56 differentially expressed genes (DEG) that had not previously been related to sex differentiation in fish were identified within the first three months of age, of which 44 were associated with ovarian differentiation (e.g., cd98, gpd1 and cry2), and 12 with testicular differentiation (e.g., ace, capn8 and nxph1). To identify putative sex determining genes, ∼4.000 DEG in juvenile gonads were mapped and their positions compared with that of previously identified sex- and growth-related quantitative trait loci (QTL). Although no genes mapped to the previously identified sex-related QTLs, two genes (foxl2 and 17βhsd) of the canonical reproduction-related genes mapped to growth-QTLs in linkage group (LG) 15 and LG6, respectively, suggesting that these genes are related to the growth dimorphism in this species.

Mapping QTL for Resistance Against Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), leads to mass mortality in mariculture. However, phenotypic selection for resistance against VNN is very difficult. To facilitate marker-assisted selection (MAS) for resistance against VNN and understanding of the genetic architecture underlying the resistance against this disease, we mapped quantitative trait loci (QTL) for resistance against VNN in Asian seabass. We challenged fingerlings at 37 days post-hatching (dph), from a single back-cross family, with NNV at a concentration of 9 × 10(6) TCID50/ml for 2 h. Daily mortalities were recorded and collected. A panel of 330 mortalities and 190 surviving fingerlings was genotyped using 149 microsatellites with 145 successfully mapped markers covering 24 linkage groups (LGs). Analysis of QTL for both resistance against VNN and survival time was conducted using interval mapping. Five significant QTL located in four LGs and eight suggestive QTL in seven LGs were identified for resistance. Another five significant QTL in three LGs and five suggestive QTL in three LGs were detected for survival time. One significant QTL, spanning 3 cM in LG20, was identified for both resistance and survival time. These QTL explained 2.2-4.1% of the phenotypic variance for resistance and 2.2-3.3% of the phenotypic variance for survival time, respectively. Our results suggest that VNN resistance in Asian seabass is controlled by many loci with small effects. Our data provide information for fine mapping of QTL and identification of candidate genes for a better understanding of the mechanism of disease resistance.

Gene expression analysis at the onset of sex differentiation in turbot (Scophthalmus maximus)

BACKGROUND

Controlling sex ratios is essential for the aquaculture industry, especially in those species with sex dimorphism for relevant productive traits, hence the importance of knowing how the sexual phenotype is established in fish. Turbot, a very important fish for the aquaculture industry in Europe, shows one of the largest sexual growth dimorphisms amongst marine cultured species, being all-female stocks a desirable goal for the industry. Although important knowledge has been achieved on the genetic basis of sex determination (SD) in this species, the master SD gene remains unknown and precise information on gene expression at the critical stage of sex differentiation is lacking. In the present work, we examined the expression profiles of 29 relevant genes related to sex differentiation, from 60 up to 135 days post fertilization (dpf), when gonads are differentiating. We also considered the influence of three temperature regimes on sex differentiation.

RESULTS

The first sex-related differences in molecular markers could be observed at 90 days post fertilization (dpf) and so we have called that time the onset of sex differentiation. Three genes were the first to show differential expression between males and females and also allowed us to sex turbot accurately at the onset of sex differentiation (90 dpf): cyp19a1a, amh and vasa. The expression of genes related to primordial germ cells (vasa, gsdf, tdrd1) started to increase between 75-90 dpf and vasa and tdrd1 later presented higher expression in females (90-105 dpf). Two genes placed on the SD region of turbot (sox2, fxr1) did not show any expression pattern suggestive of a sex determining function. We also detected changes in the expression levels of several genes (ctnnb1, cyp11a, dmrt2 or sox6) depending on culture temperature.

CONCLUSION

Our results enabled us to identify the first sex-associated genetic cues (cyp19a1a, vasa and amh) at the initial stages of gonad development in turbot (90 dpf) and to accurately sex turbot at this age, establishing the correspondence between gene expression profiles and histological sex. Furthermore, we profiled several genes involved in sex differentiation and found specific temperature effects on their expression.

Co-localization of growth QTL with differentially expressed candidate genes in rainbow trout

We tested whether genes differentially expressed between large and small rainbow trout co-localized with familial QTL regions for body size. Eleven chromosomes, known from previous work to house QTL for weight and length in rainbow trout, were examined for QTL in half-sibling families produced in September (1 XY male and 1 XX neomale) and December (1 XY male). In previous studies, we identified 108 candidate genes for growth expressed in the liver and white muscle in a subset of the fish used in this study. These gene sequences were BLASTN aligned against the rainbow trout and stickleback genomes to determine their location (rainbow trout) and inferred location based on synteny with the stickleback genome. Across the progeny of all three males used in the study, 63.9% of the genes with differential expression appear to co-localize with the QTL regions on 6 of the 11 chromosomes tested in these males. Genes that co-localized with QTL in the mixed-sex offspring of the two XY males primarily showed up-regulation in the muscle of large fish and were related to muscle growth, metabolism, and the stress response.

Identification of Sex-Linked SNPs and Sex-Determining Regions in the Yellowtail Genome

Unlike the conservation of sex-determining (SD) modes seen in most mammals and birds, teleost fishes exhibit a wide variety of SD systems and genes. Hence, the study of SD genes and sex chromosome turnover in fish is one of the most interesting topics in evolutionary biology. To increase resolution of the SD gene evolutionary trajectory in fish, identification of the SD gene in more fish species is necessary. In this study, we focused on the yellowtail, a species widely cultivated in Japan. It is a member of family Carangidae in which no heteromorphic sex chromosome has been observed, and no SD gene has been identified to date. By performing linkage analysis and BAC walking, we identified a genomic region and SNPs with complete linkage to yellowtail sex. Comparative genome analysis revealed the yellowtail SD region ancestral chromosome structure as medaka-fugu. Two inversions occurred in the yellowtail linage after it diverged from the yellowtail-medaka ancestor. An association study using wild yellowtails and the SNPs developed from BAC ends identified two SNPs that can reasonably distinguish the sexes. Therefore, these will be useful genetic markers for yellowtail breeding. Based on a comparative study, it was suggested that a PDZ domain containing the GIPC protein might be involved in yellowtail sex determination. The homomorphic sex chromosomes widely observed in the Carangidae suggest that this family could be a suitable marine fish model to investigate the early stages of sex chromosome evolution, for which our results provide a good starting point.

High-density genetic linkage mapping in turbot (Scophthalmus maximus L.) based on SNP markers and major sex- and growth-related regions detection

This paper describes the development of a high density consensus genetic linkage map of a turbot (Scophthalmus maximus L.) family composed of 149 mapping individuals using Single Nucleotide Polymorphisms (SNP) developed using the restriction-site associated DNA (RAD) sequencing technique with the restriction enzyme, PstI. A total of 6,647 SNPs were assigned to 22 linkage groups, which is equal to the number of chromosome pairs in turbot. For the first time, the average marker interval reached 0.3958 cM, which is equal to approximately 0.1203 Mb of the turbot genome. The observed 99.34% genome coverage indicates that the linkage map was genome-wide. A total of 220 Quantitative Traits Locus (QTLs) associated with two body length traits, two body weight traits in different growth periods and sex determination were detected with an LOD > 5.0 in 12 linkage groups (LGs), which explained the corresponding phenotypic variance (R²), ranging from 14.4–100%. Among them, 175 overlapped with linked SNPs, and the remaining 45 were located in regions between contiguous SNPs. According to the QTLs related to growth trait distribution and the changing of LGs during different growth periods, the growth traits are likely controlled by multi-SNPs distributed on several LGs; the effect of these SNPs changed during different growth periods. Most sex-related QTLs were detected at LG 21 with a linkage span of 70.882 cM. Additionally, a small number of QTLs with high feasibility and a narrow R² distribution were also observed on LG7 and LG14, suggesting that multi LGs or chromosomes might be involved in sex determination. High homology was recorded between LG21 in Cynoglossus semilaevis and turbot. This high-saturated turbot RAD-Seq linkage map is undoubtedly a promising platform for marker assisted selection (MAS) and flatfish genomics research.

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

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

Genetic basis and biotechnological manipulation of sexual dimorphism and sex determination in fish

Aquaculture has made an enormous contribution to the world food production, especially to the sustainable supply of animal proteins. The utility of diverse reproduction strategies in fish, such as the exploiting use of unisexual gynogenesis, has created a typical case of fish genetic breeding. A number of fish species show substantial sexual dimorphism that is closely linked to multiple economic traits including growth rate and body size, and the efficient development of sex-linked genetic markers and sex control biotechnologies has provided significant approaches to increase the production and value for commercial purposes. Along with the rapid development of genomics and molecular genetic techniques, the genetic basis of sexual dimorphism has been gradually deciphered, and great progress has been made in the mechanisms of fish sex determination and identification of sex-determining genes. This review summarizes the progress to provide some directive and objective thinking for further research in this field.

Combinatorial immunotherapy of sorafenib and blockade of programmed death-ligand 1 induces effective natural killer cell responses against hepatocellular carcinoma

Sorafenib, a multi-tyrosine kinase inhibitor, is a standard treatment for advanced hepatocellular carcinoma (HCC). Herein, we report that the combinatorial therapy of sorafenib and anti-programmed death-ligand 1 (PD-L1) monoclonal antibody (mAb) can be implemented with good results for HCC. Cancer mouse models were used to evaluate therapeutic efficacy and examine the immunologic mechanisms of the sorafenib/anti-PD-L1 mAb therapy. The combined administration of sorafenib and anti-PD-L1 mAb into tumor-bearing mice generated potent immune responses resulting in the complete eradication or remarkable reduction of tumor growth. In some instances, the sorafenib/anti-PD-L1 mAb therapy induced long-lasting protection against tumor rechallenges. The results indicate that NK cells but not CD4T cells or CD8 cells mediated the therapeutic efficacy of this combinatorial therapy. The overall results suggest that immunotherapy consisting of the combination of sorafenib/anti-PD-L1 mAb could be a promising new approach for treating patients with HCC.

Fine mapping and evolution of the major sex determining region in turbot (Scophthalmus maximus)

Fish sex determination (SD) systems are varied, suggesting evolutionary changes including either multiple evolution origins of genetic SD from nongenetic systems (such as environmental SD) and/or turnover events replacing one genetic system by another. When genetic SD is found, cytological differentiation between the two members of the sex chromosome pair is often minor or undetectable. The turbot (Scophthalmus maximus), a valuable commercial flatfish, has a ZZ/ZW system and a major SD region on linkage group 5 (LG5), but there are also other minor genetic and environmental influences. We here report refined mapping of the turbot SD region, supported by comparative mapping with model fish species, to identify the turbot master SD gene. Six genes were located to the SD region, two of them associated with gonad development (sox2 and dnajc19). All showed a high association with sex within families (P = 0), but not at the population level, so they are probably partially sex-linked genes, but not SD gene itself. Analysis of crossovers in LG5 using two families confirmed a ZZ/ZW system in turbot and suggested a revised map position for the master gene. Genetic diversity and differentiation for 25 LG5 genetic markers showed no differences between males and females sampled from a wild population, suggesting a recent origin of the SD region in turbot. We also analyzed associations with markers of the most relevant sex-related linkage groups in brill (S. rhombus), a closely related species to turbot; the data suggest that an ancient XX/XY system in brill changed to a ZZ/ZW mechanism in turbot.

SOX10 is a novel oncogene in hepatocellular carcinoma through Wnt/β-catenin/TCF4 cascade

SOX (high mobility group) genes play an important role in a number of developmental processes. Potential roles of SOXs have been demonstrated in various neoplastic tissues as tumor suppressors or promoters depending on tumor status and types. The aim of this study was to investigate the function role of SOXs in the human hepatocellular carcinoma (HCC). The gene expression changes of SOXs in HCC tissues compared with those in noncancerous hepatic tissues were detected using real-time quantitative reverse transcriptase polymerase chain reaction (QRT-PCR) analysis and immunohistochemistry. In addition, we identified the gene SOX10 that was significantly upregulated in HCC by QRT-PCR analysis and immunohistochemistry. Furthermore, we discovered that SOX10 promoted cancer cell proliferation in vitro, and SOX10 expression correlated with elevated β-catenin levels in HCC, and β-catenin function was required for SOX10's oncogenic effects. Mechanistically, SOX10 facilitates TCF4 to bind to β-catenin and form a stable SOX10/TCF4/β-catenin complex and trans-activate its downstream target gene. SOX10 mutations that disrupt the SOX10-β-catenin interaction partially prevent its function in tumor cells. All in all, SOX10 is a commonly activated tumor promoter that activates Wnt/β-catenin signaling in cancer cells of HCC.

Consolidation of the genetic and cytogenetic maps of turbot (Scophthalmus maximus) using FISH with BAC clones

Bacterial artificial chromosomes (BAC) have been widely used for fluorescence in situ hybridization (FISH) mapping of chromosome landmarks in different organisms, including a few in teleosts. In this study, we used BAC-FISH to consolidate the previous genetic and cytogenetic maps of the turbot (Scophthalmus maximus), a commercially important pleuronectiform. The maps consisted of 24 linkage groups (LGs) but only 22 chromosomes. All turbot LGs were assigned to specific chromosomes using BAC probes obtained from a turbot 5× genomic BAC library. It consisted of 46,080 clones with inserts of at least 100 kb and <5 % empty vectors. These BAC probes contained gene-derived or anonymous markers, most of them linked to quantitative trait loci (QTL) related to productive traits. BAC clones were mapped by FISH to unique marker-specific chromosomal positions, which showed a notable concordance with previous genetic mapping data. The two metacentric pairs were cytogenetically assigned to LG2 and LG16, and the nucleolar organizer region (NOR)-bearing pair was assigned to LG15. Double-color FISH assays enabled the consolidation of the turbot genetic map into 22 linkage groups by merging LG8 with LG18 and LG21 with LG24. In this work, a first-generation probe panel of BAC clones anchored to the turbot linkage and cytogenetical map was developed. It is a useful tool for chromosome traceability in turbot, but also relevant in the context of pleuronectiform karyotypes, which often show small hardly identifiable chromosomes. This panel will also be valuable for further integrative genomics of turbot within Pleuronectiformes and teleosts, especially for fine QTL mapping for aquaculture traits, comparative genomics, and whole-genome assembly.

A sex-associated sequence identified by RAPD screening in gynogenetic individuals of turbot (Scophthalmus maximus)

Understanding the genetic basis of sex determination mechanisms is essential for improving the productivity of farmed aquaculture fish species like turbot (Scophthalmus maximus). In culture conditions turbot males grow slower than females starting from eight months post-hatch, and this differential growth rate is maintained until sexual maturation is reached, being mature females almost twice as big as males of the same age. The goal of this study was to identify sex-specific DNA markers in turbot using comparative random amplified polymorphism DNA (RAPD) profiles in males and females to get new insights of the genetic architecture related to sex determination. In order to do this, we analyzed 540 commercial 10-mer RAPD primers in male and female pools of a gynogenetic family because of its higher inbreeding, which facilitates the detection of associations across the genome. Two sex-linked RAPD markers were identified in the female pool and one in the male pool. After the analysis of the three markers on individual samples of each pool and also in unrelated individuals, only one RAPD showed significant association with females. This marker was isolated, cloned and sequenced, containing two sequences, a microsatellite (SEX01) and a minisatellite (SEX02), which were mapped in the turbot reference map. From this map position, through a comparative mapping approach, we identified Foxl2, a relevant gene related to initial steps of sex differentiation, and Wnt4, a gene related with ovarian development, close to the microsatellite and minisatellite markers, respectively. The position of Foxl2 and Wnt4 was confirmed by linkage mapping in the reference turbot map.

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

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

Male-biased genes in catfish as revealed by RNA-Seq analysis of the testis transcriptome

Background

Catfish has a male-heterogametic (XY) sex determination system, but genes involved in gonadogenesis, spermatogenesis, testicular determination, and sex determination are poorly understood. As a first step of understanding the transcriptome of the testis, here, we conducted RNA-Seq analysis using high throughput Illumina sequencing.

Methodology/Principal Findings

A total of 269.6 million high quality reads were assembled into 193,462 contigs with a N50 length of 806 bp. Of these contigs, 67,923 contigs had hits to a set of 25,307 unigenes, including 167 unique genes that had not been previously identified in catfish. A meta-analysis of expressed genes in the testis and in the gynogen (double haploid female) allowed the identification of 5,450 genes that are preferentially expressed in the testis, providing a pool of putative male-biased genes. Gene ontology and annotation analysis suggested that many of these male-biased genes were involved in gonadogenesis, spermatogenesis, testicular determination, gametogenesis, gonad differentiation, and possibly sex determination.

Conclusion/Significance

We provide the first transcriptome-level analysis of the catfish testis. Our analysis would lay the basis for sequential follow-up studies of genes involved in sex determination and differentiation in catfish.

A combined strategy involving Sanger and 454 pyrosequencing increases genomic resources to aid in the management of reproduction, disease control and genetic selection in the turbot (Scophthalmus maximus)

Background

Genomic resources for plant and animal species that are under exploitation primarily for human consumption are increasingly important, among other things, for understanding physiological processes and for establishing adequate genetic selection programs. Current available techniques for high-throughput sequencing have been implemented in a number of species, including fish, to obtain a proper description of the transcriptome. The objective of this study was to generate a comprehensive transcriptomic database in turbot, a highly priced farmed fish species in Europe, with potential expansion to other areas of the world, for which there are unsolved production bottlenecks, to understand better reproductive- and immune-related functions. This information is essential to implement marker assisted selection programs useful for the turbot industry.

Results

Expressed sequence tags were generated by Sanger sequencing of cDNA libraries from different immune-related tissues after several parasitic challenges. The resulting database (“Turbot 2 database”) was enlarged with sequences generated from a 454 sequencing run of brain-hypophysis-gonadal axis-derived RNA obtained from turbot at different development stages. The assembly of Sanger and 454 sequences generated 52,427 consensus sequences (“Turbot 3 database”), of which 23,661 were successfully annotated. A total of 1,410 sequences were confirmed to be related to reproduction and key genes involved in sex differentiation and maturation were identified for the first time in turbot (AR, AMH, SRY-related genes, CYP19A, ZPGs, STAR FSHR, etc.). Similarly, 2,241 sequences were related to the immune system and several novel key immune genes were identified (BCL, TRAF, NCK, CD28 and TOLLIP, among others). The number of genes of many relevant reproduction- and immune-related pathways present in the database was 50–90% of the total gene count of each pathway. In addition, 1,237 microsatellites and 7,362 single nucleotide polymorphisms (SNPs) were also compiled. Further, 2,976 putative natural antisense transcripts (NATs) including microRNAs were also identified.

Conclusions

The combined sequencing strategies employed here significantly increased the turbot genomic resources available, including 34,400 novel sequences. The generated database contains a larger number of genes relevant for reproduction- and immune-associated studies, with an excellent coverage of most genes present in many relevant physiological pathways. This database also allowed the identification of many microsatellites and SNP markers that will be very useful for population and genome screening and a valuable aid in marker assisted selection programs.