Draft genomes of female and male turbot Scophthalmus maximus

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

Analysis of the implementation effects of ecological restoration projects based on carbon storage and eco-environmental quality: A case study of the Yellow River Delta, China

As an important means to address global climate change and land-use/land-cover (LULC) change, ecological restoration projects (ERPs) have a large effect on carbon storage functions and eco-environmental quality. However, the various ERPs carried out in the Yellow River Delta region have important implications for ecological security strategies in China. Therefore, based on land-use data and remote sensing image data, with the help of ArcGIS and Google Earth Engine (GEE) platforms, this study uses the Integrated Valuation of Ecosystem Services and Trade-offs (InVEST) model, an improved remote sensing ecological index (RSEI) model and other methods to deeply examine the evolutionary trends of eco-environmental quality and carbon storage during the implementation of ERPs in the Yellow River Delta and selects key implementation areas for in-depth analysis to determine the implementation effects of ERPs. Our findings suggested that the RSEI and carbon storage levels in the study area had opposite evolutionary trends from 2001 to 2020. Among them, the RSEI showed a fluctuating upwards trend (0.4461 (2001) and 0.5185 (2020)), while the total carbon stock showed a fluctuating downwards trend (30.67 Tg (2001) and 26.40 Tg (2020)). However, from 2015 to 2020, the RSEI and carbon storage were at a relatively stable level, which indirectly indicated that the ERPs carried out during the period from 2015 to 2020 had achieved a good comprehensive implementation effect. In addition, the areas with better improvement effects from 2015 to 2020 were primarily located in the mouth of the Yellow River Delta (Areas C and D), and their RSEI and the total carbon stock showed a certain upwards trend. This research can promote the formulation of the management strategy of ERPs in the Yellow River Delta, which is of tremendous importance to the ecological environmental preservation and high-quality development of the Yellow River Basin.

Genome-wide identification and functional characterization of inhibitor of nuclear factor-κB (IκB) kinase (IKK) in turbot (Scophthalmus maximus)

The inhibitor of nuclear factor-κB (IκB) kinase (IKK) is involved in a variety of intracellular cell signaling pathways and is an important component of the NF-κB signaling pathway. IKK genes have been suggested to play important roles in the innate immune response to pathogen infection in both vertebrates and invertebrates. However, little information is available about IKK genes in turbot (Scophthalmus maximus). In this study, six IKK genes were identified including SmIKKα, SmIKKα2, SmIKKβ, SmIKKε, SmIKKγ, and SmTBK1. The IKK genes of turbot showed the highest identity and similarity with Cynoglossus semilaevis. Then, phylogenetic analysis showed that the IKK genes of turbot were most closely related to C. semilaevis. In addition, IKK genes were widely expressed in all the examined tissues. Meanwhile, the expression patterns of IKK genes were investigated by QRT-PCR after Vibrio anguillarum and Aeromonas salmonicida infection. The results showed that IKK genes had varying expression patterns in mucosal tissues after bacteria infection, indicating that they may play key roles in maintaining the integrity of the mucosal barrier. Subsequently, protein and protein interaction (PPI) network analysis showed that most proteins interacting with IKK genes were located in the NF-κB signaling pathway. Finally, the double luciferase report and overexpression experiments showed that SmIKKα/SmIKKα2/SmIKKβ involved in the activation of NF-κB in turbot. In summary, our results suggested that IKK genes of turbot played important roles in the innate immune response of teleost, and provide valuable information for further study of the function of IKK genes.

Identification of stress-related genes by co-expression network analysis based on the improved turbot genome

Turbot (Scophthalmus maximus), commercially important flatfish species, is widely cultivated in Europe and China. With the continuous expansion of the intensive breeding scale, turbot is exposed to various stresses, which greatly impedes the healthy development of turbot industry. Here, we present an improved high-quality chromosome-scale genome assembly of turbot using a combination of PacBio long-read and Illumina short-read sequencing technologies. The genome assembly spans 538.22 Mb comprising 27 contigs with a contig N50 size of 25.76 Mb. Annotation of the genome assembly identified 104.45 Mb repetitive sequences, 22,442 protein-coding genes and 3,345 ncRNAs. Moreover, a total of 345 stress responsive candidate genes were identified by gene co-expression network analysis based on 14 published stress-related RNA-seq datasets consisting of 165 samples. Significantly improved genome assembly and stress-related candidate gene pool will provide valuable resources for further research on turbot functional genome and stress response mechanism, as well as theoretical support for the development of molecular breeding technology for resistant turbot varieties.

Measurement(s)	whole genome sequencing
Technology Type(s)	PacBio long-read and Illumina short-read sequencing technologies

Molecular cloning, expression analysis of the IgT gene and detection of IgT+ B cells in the half-smooth tongue sole (Cynoglossus semilaevis)

IgT is a specific Ig isotype in teleosts, which plays extremely important roles in the mucosal immunity of fish. In the present study, the membrane-bound and secretory IgT of the half-smooth tongue sole (Cynoglossus semilaevis) were identified for the first time. The V-D-J-C structure of two forms of csIgT are translated by the same Cτ gene, and the secretory tail and transmembrane domain were encoded through alternative splicing at the 3' end of the Cτ4. The CH regions of csIgT had high similarity with that of other flatfish (P. olivaceus and S. maximus). In healthy C. semilaevis, sIgT and mIgT were mainly expressed in mucus related tissues such as skin, intestine and gill. The transcript levels of sIgT and mIgT mRNA showed a significant induction in the immune-related tissues upon Vibrio Harveyi infection. A polyclonal rabbit anti-csIgT was successfully prepared using the csIgT heavy chain recombinant protein. Using this antibody, we detected the native IgT with the molecular mass at 220 kDa in skin total protein under non-reducing SDS-PAGE condition. Immunofluorescence analysis indicated that IgT⁺ B lymphocytes were intensively located in the skin, gill, intestine, and head kidney of C. semilaevis. These results suggest that IgT may participate in the immune response of C. semilaevis, which will facilitate the investigations of the immunoglobulins of marine fish.

Chromosome anchoring in Senegalese sole (Solea senegalensis) reveals sex-associated markers and genome rearrangements in flatfish

The integration of physical and high-density genetic maps is a very useful approach to achieve chromosome-level genome assemblies. Here, the genome of a male Senegalese sole (Solea senegalensis) was de novo assembled and the contigs were anchored to a high-quality genetic map for chromosome-level scaffolding. Hybrid assembled genome was 609.3 Mb long and contained 3403 contigs with a N50 of 513 kb. The linkage map was constructed using 16,287 informative SNPs derived from ddRAD sequencing in 327 sole individuals from five families. Markers were assigned to 21 linkage groups with an average number of 21.9 markers per megabase. The anchoring of the physical to the genetic map positioned 1563 contigs into 21 pseudo-chromosomes covering 548.6 Mb. Comparison of genetic and physical distances indicated that the average genome-wide recombination rate was 0.23 cM/Mb and the female-to-male ratio 1.49 (female map length: 2,698.4 cM, male: 2,036.6 cM). Genomic recombination landscapes were different between sexes with crossovers mainly concentrated toward the telomeres in males while they were more uniformly distributed in females. A GWAS analysis using seven families identified 30 significant sex-associated SNP markers located in linkage group 18. The follicle-stimulating hormone receptor appeared as the most promising locus associated with sex within a region with very low recombination rates. An incomplete penetrance of sex markers with males as the heterogametic sex was determined. An interspecific comparison with other Pleuronectiformes genomes identified a high sequence similarity between homologous chromosomes, and several chromosomal rearrangements including a lineage-specific Robertsonian fusion in S. senegalensis.

A genome-wide association study, supported by a new chromosome-level genome assembly, suggests sox2 as a main driver of the undifferentiatiated ZZ/ZW sex determination of turbot (Scophthalmus maximus)

BACKGROUND

Understanding sex determination (SD) across taxa is a major challenge for evolutionary biology. The new genomic tools are paving the way to identify genomic features underlying SD in fish, a group frequently showing limited sex chromosome differentiation and high SD evolutionary turnover. Turbot (Scophthalmus maximus) is a commercially important flatfish with an undifferentiated ZW/ZZ SD system and remarkable sexual dimorphism. Here we describe a new long-read turbot genome assembly used to disentangle the genetic architecture of turbot SD by combining genomics and classical genetics approaches.

RESULTS

The new turbot genome assembly consists of 145 contigs (N50 = 22.9 Mb), 27 of them representing >95% of its estimated genome size. A genome wide association study (GWAS) identified a ~ 6.8 Mb region on chromosome 12 associated with sex in 69.4% of the 36 families analyzed. The highest associated markers flanked sox2, the only gene in the region showing differential expression between sexes before gonad differentiation. A single SNP showed consistent differences between Z and W chromosomes. The analysis of a broad sample of families suggested the presence of additional genetic and/or environmental factors on turbot SD.

CONCLUSIONS

The new chromosome-level turbot genome assembly, one of the most contiguous fish assemblies to date, facilitated the identification of sox2 as a consistent candidate gene putatively driving SD in this species. This chromosome SD system barely showed any signs of differentiation, and other factors beyond the main QTL seem to control SD in a certain proportion of families.

Development of whole-genome multiplex assays and construction of an integrated genetic map using SSR markers in Senegalese sole

The Senegalese sole (Solea senegalensis) is an economically important flatfish species. In this study, a genome draft was analyzed to identify microsatellite (SSR) markers for whole-genome genotyping. A subset of 224 contigs containing SSRs were preselected and validated by using a de novo female hybrid assembly. Overall, the SSR density in the genome was 886.7 markers per megabase of genomic sequences and the dinucleotide motif was the most abundant (52.4%). In silico comparison identified a set of 108 SSRs (with di-, tetra- or pentanucleotide motifs) widely distributed in the genome and suitable for primer design. A total of 106 markers were structured in thirteen multiplex PCR assays (with up to 10-plex) and the amplification conditions were optimized with a high-quality score. Main genetic diversity statistics and genotyping reliability were assessed. A subset of 40 high polymorphic markers were selected to optimize four supermultiplex PCRs (with up to 11-plex) for pedigree analysis. Theoretical exclusion probabilities and real parentage allocation tests using parent–offspring information confirmed their robustness and effectiveness for parental assignment. These new SSR markers were combined with previously published SSRs (in total 229 makers) to construct a new and improved integrated genetic map containing 21 linkage groups that matched with the expected number of chromosomes. Synteny analysis with respect to C. semilaevis provided new clues on chromosome evolution in flatfish and the formation of metacentric and submetacentric chromosomes in Senegalese sole.