Combined ATAC-seq, RNA-seq, and GWAS analysis reveals glycogen metabolism regulatory network in Jinjiang oyster ( Crassostrea ariakensis).

Glycogen serves as the principal energy reserve for metabolic processes in aquatic shellfish and substantially contributes to the flavor and quality of oysters. The Jinjiang oyster ( Crassostrea ariakensis) is an economically and ecologically important species in China. In the present study, RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) were performed to investigate gene expression and chromatin accessibility variations in oysters with different glycogen contents. Analysis identified 9 483 differentially expressed genes (DEGs) and 7 215 genes with significantly differential chromatin accessibility (DCAGs) were obtained, with an overlap of 2 600 genes between them. Notably, a significant proportion of these genes were enriched in pathways related to glycogen metabolism, including "Glycogen metabolic process" and "Starch and sucrose metabolism". In addition, genome-wide association study (GWAS) identified 526 single nucleotide polymorphism (SNP) loci associated with glycogen content. These loci corresponded to 241 genes, 63 of which were categorized as both DEGs and DCAGs. This study enriches basic research data and provides insights into the molecular mechanisms underlying the regulation of glycogen metabolism in C. ariakensis.

Molecular phylogenetic relationships based on mitochondrial genomes of novel deep-sea corals (Octocorallia: Alcyonacea): Insights into slow evolution and adaptation to extreme deep-sea environments.

A total of 10 specimens of Alcyonacea corals were collected at depths ranging from 905 m to 1 633 m by the manned submersible Shenhai Yongshi during two cruises in the South China Sea (SCS). Based on mitochondrial genomic characteristics, morphological examination, and sclerite scanning electron microscopy, the samples were categorized into four suborders (Calcaxonia, Holaxonia, Scleraxonia, and Stolonifera), and identified as 9 possible new cold-water coral species. Assessments of GC-skew dissimilarity, phylogenetic distance, and average nucleotide identity (ANI) revealed a slow evolutionary rate for the octocoral mitochondrial sequences. The nonsynonymous ( Ka) to synonymous ( Ks) substitution ratio ( Ka/ Ks) suggested that the 14 protein-coding genes (PCGs) were under purifying selection, likely due to specific deep-sea environmental pressures. Correlation analysis of the median Ka/ Ks values of five gene families and environmental factors indicated that the genes encoding cytochrome b (cyt b) and DNA mismatch repair protein ( mutS) may be influenced by environmental factors in the context of deep-sea species formation. This study highlights the slow evolutionary pace and adaptive mechanisms of deep-sea corals.

Genome-wide DNA methylation and gene expression patterns of androgenetic haploid tiger pufferfish (Takifugu rubripes) provide insights into haploid syndrome.

Androgenesis is an important chromosome set manipulation technique used in sex control in aquaculture. Haploid embryos exhibit haploid syndrome with body abnormalities and even die during early embryonic development. In this study, we used whole genome bisulfite sequencing (WGBS) to investigate the genome-wide DNA methylation profiles in haploid females (1n-X) and males (1n-Y), and diploid females (2n-XX) and males (2n-XY) of tiger pufferfish (Takifugu rubripes), an economically important fish in China. A total of 96.32 Gb clean data was produced. Differentially methylated regions (DMRs) were found between haploids and diploids, which may be related to abnormal development and early embryonic death in haploids. There were 3,641 hyper-methylated differentially methylated genes (DMGs) and 2,179 hypo-methylated DMGs in haploid vs. diploid comparisons in both females and males. These DMGs were mainly related to genomic stability maintenance and cell cycle regulation. slf1, actr8, gas2, and pbrm1 genes were selected to validate the methylation sequencing. After combining the methylation data with the corresponding transcriptome data, we identified several genes, including guca2a, myoc, fezf2, rprml, telo2, s100a1, and marveld1, which exhibited differential expression levels modulated by DNA methylation. In conclusion, our study revealed different methylation and expression profiles between haploid and diploid T. rubripes for the first time. Several DMGs were identified between different ploidy levels, which may be related to haploid syndrome formation. The results expand the understanding of the effects of ploidy on the early development of teleosts and provide knowledge about target genes and networks to improve the survival rate of haploids.

Identification of global alternative splicing and sex-specific splicing via comparative transcriptome analysis of gonads of Chinese tongue sole ( Cynoglossus semilaevis).

Chinese tongue sole ( Cynoglossus semilaevis) is an economically important marine fish species with a ZZ/ZW sex determination mechanism, which can be influenced by temperature. Alternative splicing (AS) is an important mechanism regulating the expression of genes related to sex determination and gonadal differentiation, but has rarely been reported in fish. In this study, to explore the molecular regulatory mechanisms of sex determination and gonadal differentiation, we combined isoform and RNA sequencing (Iso-Seq and RNA-Seq) to perform transcriptome profiling of male and female gonads in C. semilaevis. In total, 81 883 and 32 341 full-length transcripts were obtained in males and females, respectively. A total of 8 279 AS genes were identified, including 2 639 genes showing differential AS (DAS) between males and females. Many intersecting DAS genes and differentially expressed genes (DEGs) were enriched in the meiotic cell cycle pathway, and genes related to gonadal differentiation, such as esrrb and wt1a, were found to have sex-specific isoforms. Thus, this study revealed AS events in the gonadal transcriptomes of male and female C. semilaevis, described the characteristics of active transcription in the testes, and identified candidate genes for studying the regulatory mechanisms of AS during gonadal differentiation.

Draft genomes of female and male turbot Scophthalmus maximus.

Turbot (Scophthalmus maximus) is a commercially important flatfish species in aquaculture. It has a drastic sexual dimorphism, with females growing faster than males. In the present study, we sequenced and de novo assembled female and male turbot genomes. The assembled female genome was 568 Mb (scaffold N50, 6.2 Mb, BUSCO 97.4%), and the male genome was 584 Mb (scaffold N50, 5.9 Mb, BUSCO 96.6%). Using two genetic maps, we anchored female scaffolds representing 535 Mb onto 22 chromosomes. Annotation of the female anchored genome identified 87.8 Mb transposon elements and 20,134 genes. We identified 17,936 gene families, of which 369 gene families were flatfish specific. Phylogenetic analysis showed that the turbot, Japanese flounder and Chinese tongue sole form a clade that diverged from other teleosts approximately 78 Mya. This report of female and male turbot draft genomes and annotated genes provides a new resource for identifying sex determination genes, elucidating the evolution of adaptive traits in flatfish and developing genetic techniques to increase the sustainability of turbot aquaculture.

Changes in DNA methylation during epigenetic-associated sex reversal under low temperature in Takifugu rubripes.

DNA methylation has frequently been implicated in sex determination and differentiation in teleost species. In order to detect the DNA methylation patterns established during sexual differentiation in tiger pufferfish T. rubripes, we performed comprehensive whole genome methylation sequencing and analyses of the gonads of male, female, and pseudo male. We obtained a total of 33.12, 32.44, and 31.60 Gb clean data for male, female, and pseudo male, with a sequencing depth of 66.44×, 60.47× and 54.86×, respectively. The methylation level of cytosine (C) residues in the genomic DNA from gonads was 11.016%, 10.428%, and 11.083% in male, female, and pseudo male, respectively. More than 65% of C methylation was at CpG sites, and less than 1% was at CHG and CHH sites. In each regulatory element, there were low methylation levels on both sides of the transcription start site, and higher methylation levels in exons, introns, and downstream of genes. The highest mCpG was on chromosome 8 and the lowest mCpG was on chromosome 5. Comparisons of whole-genome DNA methylation between pairs of samples revealed that there were 3,173 differentially methylated regions (DMRs) between female and male, and 3,037 DMRs between male and pseudo male, corresponding to 0.232% and 0.223% of the length of the genome, respectively. There were only 1,635 DMRs between female and pseudo male, representing 0.127% of the length of the genome. A number of differentially methylated genes (DMGs) related to sex determination and differentiation were selected, such as amhr2 and pfcyp19a. After Bisulfite Sequencing PCR (BSP) verification, amhr2 was exhibited low methylation level in normal males and pseudo male, and high methylation level in normal females but pfcyp19a showed low methylation level in normal females and high methylation level in normal males and pseudo males. These results provide information about the molecular epigenetic mechanisms of DNA methylation during low-temperature induced masculinization of tiger pufferfish, and increase our understanding of the mechanisms of sex determination and differentiation in this important aquaculture fish species.

Molecular characterization and functional divergence of two Gadd45g homologs in sex determination in half-smooth tongue sole (Cynoglossus semilaevis).

The growth arrest and DNA-damage-inducible protein 45 gamma (Gadd45g) is known to play a major role in embryonic development and sex determination. In this study, two Gadd45g genes were isolated from half-smooth tongue sole (Cynoglossus semilaevis). Using chromosomal fluorescence in situ hybridization (FISH), Gadd45g1 and Gadd45g2 were located on the W and Z chromosomes, respectively. The full-length cDNA sequences of Gadd45g1 (1270bp) and Gadd45g2 (1181bp) were predicted to contain a 480-bp coding sequence that could encode a protein of 159 amino acids residues. A phylogenetic tree showed that the predicted Gadd45g1 and Gadd45g2 amino acid sequences clustered closely in one branch. It is proposed that Gadd45g1 and Gadd45g2 are paralogous genes derived from the divergence of the sex chromosome. Ka/Ks ratios indicated that Gadd45g1 and Gadd45g2 may have undergone a high number of mutations and have a divergence time of only about 68,000years, although Gadd45g homologs are highly conserved. The qRT-PCR demonstrated that Gadd45g1 and Gadd45g2 were highly expressed in ovary, and negligibly expressed in testis of male and neo-male. During development of the ovary (from 80 to 150days), the expression levels of both genes reached high levels. Gadd45g1 was also highly expressed at 50days, the stage just before gonad differentiation in C. semilaevis. All these findings imply functional divergence of the two Gadd45g homologs; Gadd45g1 may be necessary for sex differentiation in the early stage of gonad development, and then Gadd45g1 and Gadd45g2 maintain ovary development and the female character of half-smooth tongue sole.

[Spatial distribution of methane in surface water and sediment of Jiulongjiang estuary and the effect environment factors of it].

Distribution of methane in surface water and sediment of Jiulongjiang Estuary was investigated during July, 2009 through head-space method. The concentration of methane varies from 10.7 to 456.7 nmol x L(-1) in the surface water at 56 sampled stations, and supersaturates relative to equilibrium with atmospheric methane. The concentration of methane decreases rapidly from estuarine upside margin to the open coastal ocean, resulting from mixing between high CH4-containing fresh water and low CH4-containing seawater. The sediment cores are situated in the upper estuarine coast and seaward boundary along the estuarine salinity gradient, representing the freshwater, half-brackish and marine water environment. Distribution of methane in porewater is consistent with that of surface water, which decreases rapidly from B1, B2, B3 to B4 stations, from 2 212 micromol x L(-1) to 5 micromol x L(-1). The concentration of sulfate in porewater increases gradually from B1, B2, B3 to B4 stations, with average value of 0.13, 0.64, 5.3 and 16.3 mmol x L(-1) respectively. The trends of methane in surface water and porewater have illustrated a large amount of methane is generated via the process of organic matter degradation mediated by methanogens, moved across sediment-water interface, and entered to overlying water. In seaward boundary sediment with an abundance of sulphate in sediment, and sulphate in porewater inhibits the methanogenesis, the methane input from the sediment rapidly decreases. Depth profiles of methane in porewater B2 and B3 stations show an increase in concentration from 43 and 10 micromol x L(-1) near the sediment-water interface to about 1 051 and 57 micromol L(-1) at core end. According to the vertical profile of methane, total organic carbon (TOC) and sulfate trend, a large amount of methane is depleted via anoxic oxidation in methane-sulfate transition. The methane released from the low concentration of sulfate sediment intertidal wetland situated in upper estuarine could be the most important source in Jiulongjiang estuary.

Induction of mitogynogenetic diploids and identification of WW super-female using sex-specific SSR markers in half-smooth tongue sole (Cynoglossus semilaevis).

The half-smooth tongue sole (Cynoglossus semilaevis) is an important cultured marine fish as well as a promising model fish for the study of sex determination mechanisms. In the present study, a protocol for artificial induction of mitogynogenesis by hydrostatic pressure using heterologous sperm was developed in half-smooth tongue sole in order to assess homozygosity of gynogens and to identify WW super-female. The optimal initiation time for pressure shock of mitogynogenetic embryos was determined to be 21.5 min after insemination when water temperature is at 22-23°C, while the optimal pressure and treatment duration were determined to be 70 MPa for 4 min. About 1,500 mitogynogenetic diploid larvae were obtained. Ten tongue sole microsatellite markers were used for homozygosity analysis of 24 mitogynogenetic larvae. Among the 24 larvae, the percentage of homozygosity ranged from 73.91% to 87.50% with an average homozygosity of 80.54%. Sex-specific simple sequence repeat (SSR) markers, CseF-SSR1, were isolated and used for identifying WW super-female mitogynogens in the tongue sole. The amplification of genomic DNA using the sex-specific SSR marker produced one DNA band of 206 bp in ZZ males, two DNA bands of 206 and 218 bp in ZW females, and one DNA band of 218 bp in WW super-females. Four WW "super-female" gynogens were observed in 39 mitogynogenetic diploids, indicating a ZW sex determination mechanism in the tongue sole. Thus, a protocol for the induction of artificial mitogynogenesis has been developed for the first time in half-smooth tongue sole, and the WW super-female diploids were identified in the mitogynogens by sex-specific SSR markers. These findings lay the foundation and provide important tool for the elaboration of sex determination mechanism, generation of WW super-females, and development of clone line and breeding of all-female stock in the half-smooth tongue sole.

Comparison of chromosome preparation methods for the different developmental stages of the half-smooth tongue sole, Cynoglossus semilaevis.

Here we present a comparison of the different chromosome preparation methods, including the method of the single embryo or larva, the method of juvenile-swimming and the method of phytohemagglutinin injection, used for the different developmental stages of the half-smooth tongue sole, Cynoglossus semilaevis. The mean index of mitosis for the three methods was 0.79%, 0.09% and 0.15%, respectively. From the well-spread metaphase chromosomes obtained, it was determined that C. semilaevis has a diploid chromosome number of 42 and heterotypic sex chromosomes. Also, the metaphase chromosomes produced could be used for karyotype preparation and banding studies. The three methods are also successfully used for chromosome preparation in Japanese flounder (Paralichthys olivaceus) and summer flounder (Paralichthys dentatus).

Construction of two BAC libraries from half-smooth tongue sole Cynoglossus semilaevis and identification of clones containing candidate sex-determination genes.

Half-smooth tongue sole (Cynoglossus semilaevis) is an increasingly important aquaculture species in China. It is also a tractable model to study sex chromosome evolution and to further elucidate the mechanism of sex determination in teleosts. Two bacterial artificial chromosome (BAC) libraries for C. semilaevis, with large, high-quality inserts and deep coverage, were constructed in the BamHI and HindIII sites of the vector pECBAC1. The two libraries contain a total of 55,296 BAC clones arrayed in 144 384-well microtiter plates and correspond to 13.36 haploid genome equivalents. The combined libraries have a greater than 99% probability of containing any single-copy sequence. Screening high-density arrays of the libraries with probes for female-specific markers and sex-related genes generated between 4-46 primary positive clones per probe. Thus, the two BAC libraries of C. semilaevis provided a readily useable platform for genomics research, illustrated by the isolation of sex determination gene(s).

Construction of a genetic linkage map and mapping of a female-specific DNA marker in half-smooth tongue sole (Cynoglossus semilaevis).

The half-smooth tongue sole (Cynoglossus semilaevis, hereafter, "tongue sole") is a marine flatfish with great commercial importance for fisheries and aquaculture in China. It has also been a promising model for the study of sex determination mechanisms in fish. Here, we report the construction of a genetic linkage map for the tongue sole, based on 137 markers including 103 AFLP markers, 33 microsatellite markers, and one female-specific DNA marker. Twenty-six linkage groups (LGs) were found. The total map length was 934.6 cM (Kosambi), with an average spacing of 8.4 cM, covering 64.4% of the estimated genome size. Furthermore, a female-specific SCAR marker, CseF-382, was mapped on LG5. This study represents the first genetic linkage map in the tongue sole. This map has great potential in the identification of quantitative traits loci and sex-related genes and marker-assisted selection in the tongue sole. Meanwhile, the new set of polymorphic microsatellite markers developed in this study is not only useful for genetic mapping but also of critical importance for studies on genetic diversity and broodstock management in tongue sole.

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

Half-smooth tongue sole (Cynoglossus semilaevis) is an important cultured marine fish as well as a promising model fish for the study of sex determination mechanisms. In the present study, a protocol for artificial gynogenesis of half-smooth tongue sole was developed in order to identify the sex determination mechanism and to generate all-female stock. The optimal UV-irradiation dose for genetically inactivating sea perch spermatozoa was determined to be > or =30 mJ/cm(2). The optimal initiation time for cold shock of gynogenetic embryos was determined to be 5 min after fertilization, while the optimal temperature and treatment duration were determined to be 20-25 min at 5 degrees C. Chromosomes from common diploids, gynogenetic haploids, and diploids were analyzed. WW chromosomes were discovered in some of the gynogenetic diploids. The microsatellite marker was applied to analyze gynogenetic diploid fry. Among the 30 gynogenetic diploid fry, 11 fry contained only one allele, while 19 contained two alleles, which had the same genotype as their mother. The female-specific DNA marker was observed in four individuals out of ten gynogenetic diploid fry. Ploidy analysis of 20 putative gynogenetic fry showed them all to be diploid. Thus, a protocol for the induction of artificial gynogenesis has been developed for the first time in half smooth tongue sole, and the sex determination mechanism in the tongue sole was determined to be female heterogametic with the ZW chromosome.