Whole genome sequencing of turbot (Scophthalmus maximus; Pleuronectiformes): a fish adapted to demersal life

De novo genome hybrid assembly and annotation of the endangered and euryhaline fish Aphanius iberus (Valenciennes, 1846) with identification of genes potentially involved in salinity adaptation

BACKGROUND

The sequencing of non-model species has increased exponentially in recent years, largely due to the advent of novel sequencing technologies. In this study, we construct the Reference Genome of the Spanish toothcarp (Aphanius iberus (Valenciennes, 1846)), a renowned euryhaline fish species. This species is native to the marshes along the Mediterranean coast of Spain and has been threatened with extinction as a result of habitat modification caused by urbanization, agriculture, and its popularity among aquarium hobbyists since the mid-twentieth century. It is also one of the first Reference Genome for Euro-Asian species within the globally distributed order Cyprinodontiformes. Additionally, this effort aims to enhance our comprehension of the species' evolutionary ecology and history, particularly its remarkable adaptations that enable it to thrive in diverse and constantly changing inland aquatic environments.

RESULTS

A hybrid assembly approach was employed, integrating PacBio long-read sequencing with Illumina short-read data. In addition to the assembly, an extensive functional annotation of the genome is provided by using AUGUSTUS, and two different approaches (InterProScan and Sma3s). The genome size (1.15 Gb) is consistent with that of the most closely related species, and its quality and completeness, as assessed with various methods, exceeded the suggested minimum thresholds, thus confirming the robustness of the assembly. When conducting an orthology analysis, it was observed that nearly all genes were grouped in orthogroups that included genes of genetically similar species. GO Term annotation revealed, among others, categories related with salinity regulation processes (ion transport, transmembrane transport, membrane related terms or calcium ion binding).

CONCLUSIONS

The integration of genomic data with predicted genes presents future research opportunities across multiple disciplines, such as physiology, reproduction, disease, and opens up new avenues for future studies in comparative genomic studies. Of particular interest is the investigation of genes potentially associated with salinity adaptation, as identified in this study. Overall, this study contributes to the growing database of Reference Genomes, provides valuable information that enhances the knowledge within the order Cyprinodontiformes, and aids in improving the conservation status of threatened species by facilitating a better understanding of their behavior in nature and optimizing resource allocation towards their preservation.

Unravelling turbot (Scophthalmus maximus) resistance to Aeromonas salmonicida: transcriptomic insights from two full-sibling families with divergent susceptibility

Introduction

Furunculosis, caused by the gram-negative bacterium Aeromonas salmonicida subsp. salmonicida, remains a significant threat to turbot (Scophthalmus maximus) aquaculture. Identifying genetic backgrounds with enhanced disease resistance is critical for improving aquaculture health management, reducing antibiotic dependency, and mitigating economic losses.

Methods

In this study, five full-sibling turbot families were challenged with A. salmonicida, which revealed one family with significantly greater resistance. Transcriptomic analyses (RNA-Seq) were performed on resistant and susceptible families, examining both naïve and 24-h postinfection (hpi) samples from head kidney and liver tissues.

Results

In the absence of infection, differentially expressed genes (DEGs) were identified predominantly in the liver. Following infection, a marked increase in DEGs was observed in the head kidney, with many genes linked to immune functions. Interestingly, the resistant family displayed a more controlled inflammatory response and upregulation of genes related to antigen presentation and T-cell activity in the head kidney at early infection stages, which may have contributed to its increased survival rate. In the liver, transcriptomic differences between the families were associated mainly with cytoskeletal organization, cell cycle regulation, and metabolic processes, including insulin signalling and lipid metabolism, regardless of infection status. Additionally, many DEGs overlapped with previously identified quantitative trait loci (QTLs) associated with resistance to A. salmonicida, providing further insights into the genetic basis of disease resistance.

Discussion

This study represents the first RNA-Seq analysis comparing resistant and susceptible turbot families and contributes valuable knowledge for the development of selective breeding programs targeting disease resistance in turbot and other aquaculture species susceptible to A. salmonicida.

LncRNA (BCO1-AS) regulate inflammatory responses in bacterial infection through caspase-1 in turbot (Scophthalmus maximus)

LncRNA plays key role in several biological processes, including transcriptional regulation, post transcriptional control and epigenetic regulation. However, research on the functional roles of lncRNAs in teleost species remains limited. Here, we discovered a lncRNA (BCO1-AS) with a critical role in antibacterial responses. Briefly, the full length of BCO1-AS was 2005 bp. Subsequently, BCO1-AS was distributed throughout the nucleus, where it may either trans- or cis-regulate the nearby genes. In addition, BCO1-AS was widely expressed in all the examined tissues with the highest expression level in intestine, while the lowest expression level was detected in muscle. Moreover, following Vibrio anguillarum challenge, BCO1-AS was significantly down-regulated in intestine, and up-regulated in gill and skin. In CHIRP experiment, BCO1-AS could effectively enrich RNA and might interact with several immune-related genes. Furthermore, we found that LPS could induce the expression of BCO1-AS. Finally, BCO1-AS could positively regulate caspase-1 at the mRNA and protein level. The BCO1-AS was speculated to inhibit the synthesis of inflammatory components. In summary, these results showed the roles of BCO1-AS in the regulation of inflammatory in turbot, which provided valuable information for further understanding the immune regulation network of lncRNA in teleost fish.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

The Expansion of Sirtuin Gene Family in Gilthead Sea Bream (Sparus aurata)-Phylogenetic, Syntenic, and Functional Insights across the Vertebrate/Fish Lineage

The Sirtuin (SIRT1-7) family comprises seven evolutionary-conserved enzymes that couple cellular NAD availability with health, nutrition and welfare status in vertebrates. This study re-annotated the sirt3/5 branch in the gilthead sea bream, revealing three paralogues of sirt3 (sirt3.1a/sirt3.1b/sirt3.2) and two of sirt5 (sirt5a/sirt5b) in this Perciform fish. The phylogeny and synteny analyses unveiled that the Sirt3.1/Sirt3.2 dichotomy was retained in teleosts and aquatic-living Sarcopterygian after early vertebrate 2R whole genome duplication (WGD). Additionally, only certain percomorphaceae and gilthead sea bream showed a conserved tandem-duplicated synteny block involving the mammalian-clustered sirt3.1 gene (psmd13-sirt3.1a/b-drd4-cdhr5-ctsd). Conversely, the expansion of the Sirt5 branch was shaped by the teleost-specific 3R WGD. As extensively reviewed in the literature, human-orthologues (sirt3.1/sirt5a) showed a high, conserved expression in skeletal muscle that increased as development advanced. However, recent sirt3.2 and sirt5b suffered an overall muscle transcriptional silencing across life, as well as an enhanced expression on immune-relevant tissues and gills. These findings fill gaps in the ontogeny and differentiation of Sirt genes in the environmentally adaptable gilthead sea bream, becoming a good starting point to advance towards a full understanding of its neo-functionalization. The mechanisms originating from these new paralogs also open new perspectives in the study of cellular energy sensing processes in vertebrates.

Transcriptomic analysis of turbot (Scophthalmus maximus) treated with zymosan a reveals that lncRNAs and inflammation-related genes mediate the protection conferred against Aeromonas salmonicida

Aeromonas salmonicida is one of the most harmful pathogens in finfish aquaculture worldwide. Immunostimulants such as β-glucans are used to enhance the immunity of cultured fish. However, their effects on fish physiology are not completely understood. In the present work, we evaluated the effect of a single intraperitoneal (ip) injection of zymosan A on fish survival against A. salmonicida infection. A single administration of this compound protected fish against A. salmonicida challenge and reduce the bacterial load in the head kidney one week after its administration. Transcriptome analyses of head kidney samples revealed several molecular mechanisms involved in the protection conferred by zymosan A and their regulation by long noncoding RNAs. The transcriptome profile of turbot exposed only to zymosan A was practically unaltered one week after ip injection. However, the administration of this immunostimulant induced significant transcriptomic changes once the fish were in contact with the bacteria and increased the survival of the infected turbot. Our results suggest that the restraint of the infection-induced inflammatory response, the management of apoptotic cell death, cell plasticity and cellular processes involving cytoskeleton dynamics support the protective effects of zymosan A. All this information provides insights on the cellular and molecular mechanisms involved in the protective effects of this widely used immunostimulant.

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.

Insights into the Mitochondrial Genetic Makeup and Miocene Colonization of Primitive Flatfishes (Pleuronectiformes: Psettodidae) in the East Atlantic and Indo-West Pacific Ocean

The mitogenomic evolution of the Psettodes flatfishes is still poorly known from their range distribution in eastern Atlantic and Indo-West Pacific Oceans. The study delves into the matrilineal evolutionary pathway of these primitive flatfishes, with a specific focus on the complete mitogenome of the Psettodes belcheri species, as determined through next-generation sequencing. The mitogenome in question spans a length of 16,747 base pairs and comprises a total of 37 genes, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. Notably, the mitogenome of P. belcheri exhibits a bias towards AT base pairs, with a composition of 54.15%, mirroring a similar bias observed in its close relative, Psettodes erumei, which showcases percentages of 53.07% and 53.61%. Most of the protein-coding genes commence with an ATG initiation codon, except for Cytochrome c oxidase I (COI), which initiates with a GTG codon. Additionally, four protein-coding genes commence with a TAA termination codon, while seven others exhibit incomplete termination codons. Furthermore, two protein-coding genes, namely NAD1 and NAD6, terminate with AGG and TAG stop codons, respectively. In the mitogenome of P. belcheri, the majority of transfer RNAs demonstrate the classical cloverleaf secondary structures, except for tRNA-serine, which lacks a DHU stem. Comparative analysis of conserved blocks within the control regions of two Psettodidae species unveiled that the CSB-II block extended to a length of 51 base pairs, surpassing the other blocks and encompassing highly variable sites. A comprehensive phylogenetic analysis using mitochondrial genomes (13 concatenated PCGs) categorized various Pleuronectiformes species, highlighting the basal position of the Psettodidae family and showed monophyletic clustering of Psettodes species. The approximate divergence time (35-10 MYA) between P. belcheri and P. erumei was estimated, providing insights into their separation and colonization during the early Miocene. The TimeTree analysis also estimated the divergence of two suborders, Psettodoidei and Pleuronectoidei, during the late Paleocene to early Eocene (56.87 MYA). The distribution patterns of Psettodes flatfishes were influenced by ocean currents and environmental conditions, contributing to their ecological speciation. In the face of climate change and anthropogenic activities, the conservation implications of Psettodes flatfishes are emphasized, underscoring the need for regulated harvesting and adaptive management strategies to ensure their survival in changing marine ecosystems. Overall, this study contributes to understanding the evolutionary history, genetic diversity, and conservation needs of Psettodes flatfishes globally. However, the multifaceted exploration of mitogenome and larger-scale genomic data of Psettodes flatfish will provide invaluable insights into their genetic characterization, evolutionary history, environmental adaptation, and conservation in the eastern Atlantic and Indo-West Pacific Oceans.

Non-synonymous variation and protein structure of candidate genes associated with selection in farm and wild populations of turbot (Scophthalmus maximus)

Non-synonymous variation (NSV) of protein coding genes represents raw material for selection to improve adaptation to the diverse environmental scenarios in wild and livestock populations. Many aquatic species face variations in temperature, salinity and biological factors throughout their distribution range that is reflected by the presence of allelic clines or local adaptation. The turbot (Scophthalmus maximus) is a flatfish of great commercial value with a flourishing aquaculture which has promoted the development of genomic resources. In this study, we developed the first atlas of NSVs in the turbot genome by resequencing 10 individuals from Northeast Atlantic Ocean. More than 50,000 NSVs where detected in the ~ 21,500 coding genes of the turbot genome, and we selected 18 NSVs to be genotyped using a single Mass ARRAY multiplex on 13 wild populations and three turbot farms. We detected signals of divergent selection on several genes related to growth, circadian rhythms, osmoregulation and oxygen binding in the different scenarios evaluated. Furthermore, we explored the impact of NSVs identified on the 3D structure and functional relationship of the correspondent proteins. In summary, our study provides a strategy to identify NSVs in species with consistently annotated and assembled genomes to ascertain their role in adaptation.

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.

The role of genomic signatures of directional selection and demographic history in the population structure of a marine teleost with high gene flow

Recent studies have uncovered patterns of genomic divergence in marine teleosts where panmixia due to high gene flow has been the general paradigm. These signatures of divergent selection are often impacted by structural variants, acting as "supergenes" facilitating local adaptation. The highly dispersing European plaice (Pleuronectes platessa)-in which putative structural variants (i.e., inversions) have been identified-has successfully colonized the brackish water ecosystem of the Baltic Sea. Thus, the species represents an ideal opportunity to investigate how the interplay of gene flow, structural variants, natural selection, past demographic history, and gene flow impacts on population (sub)structuring in marine systems. Here, we report on the generation of an annotated draft plaice genome assembly in combination with population sequencing data-following the salinity gradient from the Baltic Sea into the North Sea together with samples from Icelandic waters-to illuminate genome-wide patterns of divergence. Neutral markers pointed at large-scale panmixia across the European continental shelf associated with high gene flow and a common postglacial colonization history of shelf populations. However, based on genome-wide outlier loci, we uncovered signatures of population substructuring among the European continental shelf populations, i.e., suggesting signs of ongoing selection. Genome-wide selection analyses (xp-EHH) and the identification of genes within genomic regions of recent selective sweeps-overlapping with the outlier loci-suggest that these represent the signs of divergent selection. Our findings provide support for genomic divergence driven by local adaptation in the face of high gene flow and elucidate the relative importance of demographic history versus adaptive divergence in shaping the contemporary population genetic structure of a marine teleost. The role of the putative inversion(s) in the substructuring-and potentially ongoing adaptation-was seemingly not substantial.

Comprehensive identification and expression profiling of immune-related lncRNAs and their target genes in the intestine of turbot (Scophthalmus maximus L.) in response to Vibrio anguillarum infection

Long non-coding RNA (lncRNA) play vital regulatory roles in various biological processes. Intestine is one of the most sensitive organs to environmental and homeostatic disruptions for fish. However, systematic profiles of lncRNAs in the intestine of teleost in responses to pathogen infections is still limited. Turbot (Scophthalmus maximus L.), an important commercial fish species in China, has been suffering with Vibrio anguillarum infection, resulted in dramatic economic loss. Hereinto, the intestinal tissues of turbot were sampled at 0 h, 2 h, 12 h, and 48 h following V. anguillarum infection. The histopathological analysis revealed that the pathological trauma was mainly present in intestinal tunica mucosal epithelium. After high-throughput sequencing and bioinformatic analysis, a total of 9722 lncRNAs and 21,194 mRNAs were obtained, and the average length and exon number of lncRNAs were both less than those of mRNAs. Among which, a set of 158 lncRNAs and 226 mRNAs were differentially expressed (DE-lncRNAs and DEGs) in turbot intestine at three time points, related to many immune-related genes such as complement, interleukin, chemokine, lysosome, and macrophage, indicating their potential critical roles in immune responses. In addition, 2803 and 1803 GO terms were enriched for DEGs and co-expressed target genes of DE-lncRNAs, respectively. Moreover, 127 and 50 KEGG pathways including cell adhesion molecules (CAMs), phagosome, JAK-STAT signaling pathway, cytokine-cytokine receptor interaction, and intestinal immune network for IgA production, were enriched for DEGs and co-expressed target genes of DE-lncRNAs, respectively. Finally, qRT-PCR was conducted to confirm the reliability of sequencing data. The present study will set the foundation for the future exploration of lncRNA functions in teleost in response to bacterial infection.

A novel full-length transcriptome resource from multiple immune-related tissues in turbot (Scophthalmus maximus) using Pacbio SMART sequencing

Turbot (Scophthalmus maximus) is an important cold-water economic fish. However, the production and development of turbot industry has been constantly hindered by the frequent occurrence of some diseases. Lacking full-length transcriptome for turbot limits immune gene discoveries and gene structures analysis. Therefore, we generated a full-length transcriptome using mixed immune-related tissues of turbot with PacBio Sequel platform. In this study, a total of 31.7 Gb high quality data were generated with the average subreads length of 2618 bp. According to the presence of 5' and 3' primers as well as poly (A) tails, FL (Full-length) and NFL (Non-full-length) isoforms were obtained. Meanwhile, we identified 32,003 non-redundant transcripts, 76.02% of which was novel isoforms of known genes. In addition, 12,176 alternative splicing (AS) events, 6614 polyadenylation (APA) events, 1905 transcription factors, and 2703 lncRNAs were identified. This work is a comprehensive report on the full-length transcriptome of immune-related tissues of turbot, and it also provides valuable molecular resources for future research on the adaptation mechanisms and functional genomics of turbot.

CC and CXC chemokines in turbot (Scophthalmus maximus L.): Identification, evolutionary analyses, and expression profiling after Aeromonas salmonicida infection

Chemokines are a superfamily of structurally related cytokines, which exert essential roles in guiding cell migration in development, homeostasis, and immunity. CC and CXC chemokines are the two major subfamilies in teleost species. In this study, a total of seventeen CC and CXC chemokines, with inclusion of twelve CC and five CXC chemokines, were systematically identified from the turbot genome, making turbot the teleost harboring the least number of CC and CXC chemokines among all teleost species ever reported. Phylogeny, synteny, and genomic organization analyses were performed to annotate these genes, and multiple chemokine genes were identified in the turbot genome, due to the tandem duplications (CCL19 and CCL20), the whole genome duplications (CCL20, CCL25, and CXCL12), and the teleost-specific members (CCL34-36, CCL44, and CXCL18). In addition, chemokines were ubiquitously expressed in nine examined healthy tissues, with high expression levels observed in liver, gill, and spleen. Moreover, most chemokines were significantly differentially expressed in gill and spleen after Aeromonas salmonicida infection, and exhibited tissue-specific and time-dependent manner. Finally, protein-protein interaction network (PPI) analysis indicated that turbot chemokines interacted with a few immune-related genes such as interleukins, cathepsins, stats, and TLRs. These results should be valuable for comparative immunological studies and provide insights for further functional characterization of chemokines in teleost.

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

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Characterization, evolution and expression analysis of Toll-like receptor 7 (TLR7) in turbot (Scophthalmus maximus L.)

The pattern recognition receptors (PRRs) can recognize the conserved molecular structures of pathogens to active the innate immune responses, and subsequently induce the antigen-specific adaptive immune responses for the clearance of infected pathogen. Among the PRRs, Toll-like receptors (TLRs) are the first and best characterized PRRs across all the species. Among the TLR members, TLR7 showed significant conservation across the vertebrates, with the lowest rate of evolution for its LRR domains from primates to fishes. In the current study, one TLR7 (SmTLR7) gene was captured in turbot, with a 3144 bp open reading frame (ORF), that encoding 1047 amino acid residues. Following multiple sequence comparison, SmTLR7 was found to have the highest similarity and identity both to Paralichthys olivaceus with 91.9% and 85.9%, respectively. In phylogenetic analysis, SmTLR7 was firstly clustered with Japanese flounder, and then clustered with fugu, rainbow trout, and zebrafish. In addition, SmTLR7 was widely expressed in all the examined tissues with the highest expression level in spleen, followed by skin, while the lowest expression level was detected in blood. Following both Edwardsiella tarda and Vibrio anguillarum challenge, SmTLR7 was significantly down-regulated in gill and intestine, and up-regulated in skin. Moreover, SmTLR7 was significantly up-regulated in head kidney macrophages following LPS, LTA, PGN and polyI:C stimulation, as well as showed the strongest binding ability to LPS, followed by PGN, LTA, and polyI:C in a dose-dependent manner. Finally, following RNAi of SmTLR7, MyD88 and IL-1β were slightly up-regulated, while TRAF6 and IL-8 were significantly down-regulated. The characterization of TLR7 can expand our understanding of the PRRs in teleost fishes, and eventually aid the exploration of interactions between host and pathogen.

Photoreceptor distributions, visual pigments and the opsin repertoire of Atlantic halibut (Hippoglossus hippoglossus)

Fishes often have cone photoreceptors organized in lattice-like mosaic formations. In flatfishes, these lattices undergo dramatic changes during metamorphosis whereby a honeycomb mosaic of single cones in the larva is replaced by a square mosaic of single and double cones in the adult. The spatio-temporal dynamics of this transition are not well understood. Here, we describe the photoreceptors and mosaic formations that occur during the larva to juvenile transition of Atlantic halibut from the beginning of eye migration to its completion. To gauge the possibility of colour vision, visual pigments in juveniles were measured by microspectrophotometry and the opsin repertoire explored using bioinformatics. At the start of eye migration, the larva had a heterogeneous retina with honeycomb mosaic in the dorsonasal and ventrotemporal quadrants and a square mosaic in the ventronasal and dorsotemporal quadrants. By the end of metamorphosis, the square mosaic was present throughout the retina except in a centrodorsotemporal area where single, double and triple cones occurred randomly. Six cone visual pigments were found with maximum absorbance (λ_max, in nm) in the short [S(431) and S(457)], middle [M(500), M(514) and M(527)], and long [L(550)] wavelengths, and a rod visual pigment with λ_max at 491 nm. These pigments only partially matched the opsin repertoire detected by query of the Atlantic halibut genome. We conclude that the Atlantic halibut undergoes a complex re-organization of photoreceptors at metamorphosis resulting in a multi-mosaic retina adapted for a demersal life style.

Generation of a chromosome‐level genome assembly for Pacific halibut ( Hippoglossus stenolepis ) and characterization of its sex‐determining genomic region

The Pacific halibut (Hippoglossus stenolepis) is a key species in the North Pacific Ocean and Bering Sea ecosystems, where it also supports important fisheries. However, the lack of genomic resources limits our understanding of evolutionary, environmental and anthropogenic forces affecting key life history characteristics of Pacific halibut and prevents the application of genomic tools in fisheries management and conservation efforts. In the present study, we report on the first generation of a high‐quality chromosome‐level assembly of the Pacific halibut genome, with an estimated size of 602 Mb, 24 chromosome‐length scaffolds that contain 99.8% of the assembly and a N₅₀ scaffold length of 27.3 Mb. In the first application of this important resource, we conducted genome‐wide analyses of sex‐specific genetic variation by pool sequencing and characterized a potential sex‐determining region in chromosome 9 with a high density of female‐specific SNPs. Within this region, we identified the bmpr1ba gene as a potential candidate for master sex‐determining (MSD) gene. bmpr1ba is a member of the TGF‐β family that in teleosts has provided the largest number of MSD genes, including a paralogue of this gene in Atlantic herring. The genome assembly constitutes an essential resource for future studies on Pacific halibut population structure and dynamics, evolutionary history and responses to environmental and anthropogenic influences. Furthermore, the genomic location of the sex‐determining region in Pacific halibut has been identified and a putative candidate MSD gene has been proposed, providing further support for the rapid evolution of sex‐determining mechanisms in teleost fish.

Integration of Maps Enables a Cytogenomics Analysis of the Complete Karyotype in Solea senegalensis

The Pleuronectiformes order, which includes several commercially-important species, has undergone extensive chromosome evolution. One of these species is Solea senegalensis, a flatfish with 2n = 42 chromosomes. In this study, a cytogenomics approach and integration with previous maps was applied to characterize the karyotype of the species. Synteny analysis of S. senegalensis was carried out using two flatfish as a reference: Cynoglossus semilaevis and Scophthalmus maximus. Most S. senegalensis chromosomes (or chromosome arms for metacentrics and submetacentrics) showed a one-to-one macrosyntenic pattern with the other two species. In addition, we studied how repetitive sequences could have played a role in the evolution of S. senegalensis bi-armed (3, and 5–9) and acrocentric (11, 12 and 16) chromosomes, which showed the highest rearrangements compared with the reference species. A higher abundance of TEs (Transposable Elements) and other repeated elements was observed adjacent to telomeric regions on chromosomes 3, 7, 9 and 16. However, on chromosome 11, a greater abundance of DNA transposons was detected in interstitial BACs. This chromosome is syntenic with several chromosomes of the other two flatfish species, suggesting rearrangements during its evolution. A similar situation was also found on chromosome 16 (for microsatellites and low complexity sequences), but not for TEs (retroelements and DNA transposons). These differences in the distribution and abundance of repetitive elements in chromosomes that have undergone remodeling processes during the course of evolution also suggest a possible role for simple repeat sequences in rearranged regions.

Characterization and the potential immune role of class A scavenger receptor member 4 (SCARA4) in bacterial infection in turbot (Scophthalmus maximus L.)

The class A scavenger receptors play important roles in innate immunity and are distributed on plasma membrane of macrophages and other cell types. Notably, the class A scavenger receptor 4 (SCARA4) contains a typical C-type (calcium-dependent) lectin domain, which belongs to the collectin family of pattern recognition receptors and is involved in the immune response against infection. Here, one turbot SCARA4 gene was identified with a 2,292 bp open reading frame (ORF) encoding 763 amino acid residues. Multiple sequence analysis and phylogenetic analysis confirmed that SmSCARA4 gene was more close to that of P. olivaceus. Gene structure and syntenic analysis showed conserved exon/intron organization pattern and syntenic pattern across selected vertebrate species. Tissue distribution analysis showed SmSCARA4 was expressed in all the tested healthy tissues with the relative high expression levels in skin, gill and spleen. Following both E. tarda and V. anguillarum challenge in vivo, SmSCARA4 was significantly repressed in gill and intestine. Remarkably, SmSCARA4 showed the strongest binding ability to LPS and strongest upregulation in turbot head kidney macrophages in response to LPS. Knockdown and overexpression of SmSCARA4 revealed its interactions with the two pro-inflammatory cytokines, TNF-α and IL-1β. Finally, repression of SmSCARA4 via combined treatment of LPS and overexpression of SmSCARA4 construct in turbot head kidney macrophages further indicated an inhibitory role of SmSCARA4 in LPS-stimulated inflammation. Taken together, turbot SmSCARA4 plays an important role in turbot immunity, especially in the mucosa-related systems; SmSCARA4 possesses strong binding specificity to LPS, and exerts protective roles in response to LPS infection by reducing the release of pro-inflammatory cytokines. The mechanisms of inhibitory role of SmSCARA4 in LPS-elicited inflammation await further investigation.

Overview and New Insights Into the Diversity, Evolution, Role, and Regulation of Kisspeptins and Their Receptors in Teleost Fish

In the last two decades, kisspeptin (Kiss) has been identified as an important player in the regulation of reproduction and other physiological functions in vertebrates, including several fish species. To date, two ligands (Kiss1, Kiss2) and three kisspeptin receptors (Kissr1, Kissr2, Kissr3) have been identified in teleosts, likely due to whole-genome duplication and loss of genes that occurred early in teleost evolution. Recent results in zebrafish and medaka mutants have challenged the notion that the kisspeptin system is essential for reproduction in fish, in marked contrast to the situation in mammals. In this context, this review focuses on the role of kisspeptins at three levels of the reproductive, brain-pituitary-gonadal (BPG) axis in fish. In addition, this review compiled information on factors controlling the Kiss/Kissr system, such as photoperiod, temperature, nutritional status, sex steroids, neuropeptides, and others. In this article, we summarize the available information on the molecular diversity and evolution, tissue expression and neuroanatomical distribution, functional significance, signaling pathways, and gene regulation of Kiss and Kissr in teleost fishes. Of particular note are recent advances in understanding flatfish kisspeptin systems, which require further study to reveal their structural and functional diversity.

Brains in Metamorphosis: Temporal Transcriptome Dynamics in Hatchery-Reared Flatfishes

Metamorphosis is a captivating process of change during which the morphology of the larva is completely reshaped to face the new challenges of adult life. In the case of fish, this process initiated in the brain has traditionally been considered to be a critical rearing point and despite the pioneering molecular work carried out in other flatfishes, the underlying molecular basis is still relatively poorly characterized. Turbot brain transcriptome of three developmental stages (pre-metamorphic, climax of metamorphosis and post-metamorphic) were analyzed to study the gene expression dynamics throughout the metamorphic process. A total of 1570 genes were differentially expressed in the three developmental stages and we found a specific pattern of gene expression at each stage. Unexpectedly, at the climax stage of metamorphosis, we found highly expressed genes related to the immune response, while the biological pathway enrichment analysis in pre-metamorphic and post-metamorphic were related to cell differentiation and oxygen carrier activity, respectively. In addition, our results confirm the importance of thyroid stimulating hormone, increasing its expression during metamorphosis. Based on our findings, we assume that immune system activation during the climax of metamorphosis stage could be related to processes of larval tissue inflammation, resorption and replacement, as occurs in other vertebrates.

Pattern of New Gene Origination in a Special Fish Lineage, the Flatfishes

Origination of new genes are of inherent interest of evolutionary geneticists for decades, but few studies have addressed the general pattern in a fish lineage. Using our recent released whole genome data of flatfishes, which evolved one of the most specialized body plans in vertebrates, we identified 1541 (6.9% of the starry flounder genes) flatfish-lineage-specific genes. The origination pattern of these flatfish new genes is largely similar to those observed in other vertebrates, as shown by the proportion of DNA-mediated duplication (1317; 85.5%), RNA-mediated duplication (retrogenes; 96; 6.2%), and de novo-origination (128; 8.3%). The emergence rate of species-specific genes is 32.1 per Mya and the whole average level rate for the flatfish-lineage-specific genes is 20.9 per Mya. A large proportion (31.4%) of these new genes have been subjected to selection, in contrast to the 4.0% in primates, while the old genes remain quite similar (66.4% vs. 65.0%). In addition, most of these new genes (70.8%) are found to be expressed, indicating their functionality. This study not only presents one example of systematic new gene identification in a teleost taxon based on comprehensive phylogenomic data, but also shows that new genes may play roles in body planning.

Chromosome-level assembly reveals a putative Y-autosomal fusion in the sex determination system of the Greenland Halibut (Reinhardtius hippoglossoides)

Despite the commercial importance of Greenland Halibut (Reinhardtius hippoglossoides), important gaps still persist in our knowledge of this species, including its reproductive biology and sex determination mechanism. Here, we combined single-molecule sequencing of long reads (Pacific Sciences) with chromatin conformation capture sequencing (Hi-C) data to assemble the first chromosome-level reference genome for this species. The high-quality assembly encompassed more than 598 Megabases (Mb) assigned to 1 594 scaffolds (scaffold N50 = 25 Mb) with 96% of its total length distributed among 24 chromosomes. Investigation of the syntenic relationship with other economically important flatfish species revealed a high conservation of synteny blocks among members of this phylogenetic clade. Sex determination analysis revealed that, similar to other teleost fishes, flatfishes also exhibit a high level of plasticity and turnover in sex-determination mechanisms. A low-coverage whole-genome sequence analysis of 198 individuals revealed that Greenland Halibut possesses a male heterogametic XY system and several putative candidate genes implied in the sex determination of this species. Our study also suggests for the first time in flatfishes that a putative Y-autosomal fusion could be associated with a reduction of recombination typical of the early steps of sex chromosome evolution.

Identification of quantitative trait loci associated with upper temperature tolerance in turbot, Scophthalmus maximus

Temperature tolerance is an important trait from both an economic and evolutionary perspective in fish. Because of difficulties with measurements, genome-wide selection using quantitative trait loci (QTLs) affecting Upper temperature tolerance may be an alternative for genetic improvement. Turbot Scophthalmus maximus (L.) is a cold-water marine fish with high economic value in Europe and Asia. The genetic bases of upper temperature tolerance (UTTs) traits have been rarely studied. In this study, we constructed a genetic linkage map of turbot using simple sequence repeats (SSRs) and single nucleotide polymorphism (SNP) markers. A total of 190 SSR and 8,123 SNP were assigned to 22 linkage groups (LGs) of a consensus map, which spanned 3,648.29 cM of the turbot genome, with an average interval of 0.44 cM. Moreover, we re-anchored genome sequences, allowing 93.8% physical sequences to be clustered into 22 turbot pseudo-chromosomes. A high synteny was observed between two assemblies from the literature. QTL mapping and validation analysis identified thirteen QLTs which are major effect QTLs, of these, 206 linked SNP loci, and two linked SSR loci were considered to have significant QTL effects. Association analysis for UTTs with 129 QTL markers was performed for different families, results showed that eight SNP loci were significantly correlated with UTT, which markers could be helpful in selecting thermal tolerant breeds of turbot. 1,363 gene sequences were genomically annotated, and 26 QTL markers were annotated. We believe these genes could be valuable candidates affecting high temperatures, providing valuable genomic resources for the study of genetic mechanisms regulating thermal stress. Similarly, they may be used in marker-assisted selection (MAS) programs to improve turbot performance.

Comparative analysis of the miRNA-mRNA regulation networks in turbot (Scophthalmus maximus L.) following Vibrio anguillarum infection

MicroRNAs could not only regulate posttranscriptional silencing of target genes in eukaryotic organisms, but also have positive effect on their target genes as well. These microRNAs have been reported to be involved in mucosal immune responses to pathogen infection in teleost. Therefore, we constructed the immune-related miRNA-mRNA networks in turbot intestine following Vibrio anguillarum infection. In our results, 1550 differentially expressed (DE) genes and 167 DE miRNAs were identified. 113 DE miRNAs targeting 89 DE mRNAs related to immune response were used to construct miRNA-mRNA interaction networks. Functional analysis showed that target genes were associated with synthesis and degradation of ketone bodies, mucin type O-Glycan biosynthesis, homologous recombination, biotin metabolism, and intestinal immune network for IgA production that were equivalent to the function of IgT and IgM in fish intestine. Finally, 10 DE miRNAs and 7 DE mRNAs were selected for validating the accuracy of high-throughput sequencing results by qRT-PCR. The results of this study will provide valuable information for the elucidation of the regulation mechanisms of miRNA-mRNA interactions involved in disease resistance in teleost mucosal immune system.

Characterization of the turbot (Scophthalmus maximus) interleukin-18: Identification of splicing variants, phylogeny, synteny and expression analysis

Interleukin-18 (IL-18) is a pro-inflammatory cytokine that belongs to the interleukin-1 (IL-1) family of cytokines. As occurs with IL-1β, it is synthetized as an inactive precursor peptide that is mainly processed by the cysteine protease caspase-1 in the inflammasome complex. In mammals, and in collaboration with IL-12, it has been described as an important cytokine controlling the Th1-mediated immune responses through the induction of IFN-γ. Although its function in mammals is well stablished, the activity of this cytokine in teleost remains to be elucidated. This could be due, among other things, to the absence of this gene in the fish model species zebrafish, but also to its complex regulation. As it was observed for rainbow trout and human, il18 splicing variants were also found in turbot, which could represent a regulatory mechanism of its bioactivity. In the case of turbot, three splicing variants were observed (SV1-3), and one of them showed an insertion of 10 amino acids in the middle of the potential caspase-1 cleavage position, reflecting that this is probably a form resistant to the processing by the inflammasome. Phylogenetic and three-dimensional analyses of turbot Il18 revealed that it is relatively well-conserved in vertebrates, although only a partial conservation of the gene synteny was observed between fish and mammals. As it was expected, turbot il18 splicing variants were mainly expressed in immune tissues under healthy conditions, and their expression was induced by a bacterial challenge, although certain inhibitions were observed after viral and parasitic infections. In the case of the viral challenge, il18 downregulations did not seem to be due to the effect of type I IFNs.

PhylomeDB V5: an expanding repository for genome-wide catalogues of annotated gene phylogenies

PhylomeDB is a unique knowledge base providing public access to minable and browsable catalogues of pre-computed genome-wide collections of annotated sequences, alignments and phylogenies (i.e. phylomes) of homologous genes, as well as to their corresponding phylogeny-based orthology and paralogy relationships. In addition, PhylomeDB trees and alignments can be downloaded for further processing to detect and date gene duplication events, infer past events of inter-species hybridization and horizontal gene transfer, as well as to uncover footprints of selection, introgression, gene conversion, or other relevant evolutionary processes in the genes and organisms of interest. Here, we describe the latest evolution of PhylomeDB (version 5). This new version includes a newly implemented web interface and several new functionalities such as optimized searching procedures, the possibility to create user-defined phylome collections, and a fully redesigned data structure. This release also represents a significant core data expansion, with the database providing access to 534 phylomes, comprising over 8 million trees, and homology relationships for genes in over 6000 species. This makes PhylomeDB the largest and most comprehensive public repository of gene phylogenies. PhylomeDB is available at http://www.phylomedb.org.

The CC and CXC chemokine receptors in turbot (Scophthalmus maximus L.) and their response to Aeromonas salmonicida infection

Chemokines are crucial regulators of cell mobilization for development, homeostasis, and immunity. Chemokines signal through binding to chemokine receptors, a superfamily of seven-transmembrane domain G-coupled receptors. In the present study, eleven CC chemokine receptors (CCRs) and seven CXC chemokine receptors (CXCRs) were identified from turbot genome. Phylogenetic and syntenic analyses were performed to annotate these genes, indicating the closest relationship between the turbot chemokine receptors and their counterparts of Japanese flounders (Paralichthys olivaceus). Evolutionary analyses revealed that the tandem duplications of CCR8 and CXCR3, the whole genome duplications of CCR6, CCR9, CCR12, and CXCR4, and the teleost-specific CCR12 led to the expansion of turbot chemokine receptors. In addition, turbot chemokine receptors were ubiquitously expressed in nine examined healthy tissues, with high expression levels observed in spleen, gill, and head kidney. Moreover, most turbot chemokine receptors were significantly differentially expressed in spleen and gill after Aeromonas salmonicida infection, and exhibited general down-regulations at early time points and then gradually up-regulated. Finally, protein-protein interaction network (PPI) analyses indicated that chemokine receptors interacted with a few immune-related genes such as interleukins, Grk genes, CD genes, etc. These results should be valuable for comparative immunological studies and provide insights for further functional characterization of chemokine receptors in turbots.

Developmental aspects of the hypothalamic-pituitary network related to reproduction in teleost fish

The hypothalamic-pituitary-gonadal axis is the main system that regulates reproduction in vertebrates through a complex network that involves different neuropeptides, neurotransmitters, and pituitary hormones. Considering that this axis is established early on life, the main goal of the present work is to gather information on its development and the actions of its components during early life stages. This review focuses on fish because their neuroanatomical characteristics make them excellent models to study neuroendocrine systems. The following points are discussed: i) developmental functions of the neuroendocrine components of this network, and ii) developmental disruptions that may impact adult reproduction. The importance of the components of this network and their susceptibility to external/internal signals that can alter their specific early functions and/or even the establishment of the reproductive axis, indicate that more studies are necessary to understand this complex and dynamic network.

Characterization of toll-like receptor 1 (TLR1) in turbot (Scophthalmus maximus L.)

TLRs are the first and best-characterized pattern recognition receptors conserved across all the species. Different from mammals, the TLRs in teleost fishes are very diversified due to various evolutionary mechanisms. Here, we characterized one TLR1 gene in turbot, with a 2,415 bp open reading frame (ORF), that encoding 804 amino acid residues, and have the highest similarity and identity both to Paralichthys olivaceus with 88.9% and 79.9%. In phylogenetic analysis, it was firstly clustered with P. olivaceus, and then clustered with Takifugu rubripes. TLR1 was widely expressed in all the examined healthy tissues with the highest expression level in spleen, followed by head-kidney. In addition, it was significantly regulated in gill, skin and intestine following Edwardsiella tarda and Vibrio anguillarum challenge with different expression patterns. In in vitro stimulation with pathogen-associated molecular patterns, TLR1 showed significantly strong and elevated responses to LPS, but only responded to LTA and Poly(I:C) at the highest evaluated concentration, while no response was detected using PGN stimulation. Moreover, in subcellular localization analysis, TLR1 was distributed in the cytoplasm, membrane and nucleus. Taken together, TLR1 played vital roles for host immune response to bacterial infection, only with strong binding ability to LPS and involved in the production of inflammatory cytokines. However, the specific ligand for TLR1 and its functional association with other TLRs should be further characterized in fish species.

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.

Evolutionary ecology of the visual opsin gene sequence and its expression in turbot (Scophthalmus maximus)

Background

As flatfish, turbot undergo metamorphosis as part of their life cycle. In the larval stage, turbot live at the ocean surface, but after metamorphosis they move to deeper water and turn to benthic life. Thus, the light environment differs greatly between life stages. The visual system plays a great role in organic evolution, but reports of the relationship between the visual system and benthic life are rare. In this study, we reported the molecular and evolutionary analysis of opsin genes in turbot, and the heterochronic shifts in opsin expression during development.

Results

Our gene synteny analysis showed that subtype RH2C was not on the same gene cluster as the other four green-sensitive opsin genes (RH2) in turbot. It was translocated to chromosome 8 from chromosome 6. Based on branch-site test and spectral tuning sites analyses, E122Q and M207L substitutions in RH2C, which were found to be under positive selection, are closely related to the blue shift of optimum light sensitivities. And real-time PCR results indicated the dominant opsin gene shifted from red-sensitive (LWS) to RH2B1 during turbot development, which may lead to spectral sensitivity shifts to shorter wavelengths.

Conclusions

This is the first report that RH2C may be an important subtype of green opsin gene that was retained by turbot and possibly other flatfish species during evolution. Moreover, E122Q and M207L substitutions in RH2C may contribute to the survival of turbot in the bluish colored ocean. And heterochronic shifts in opsin expression may be an important strategy for turbot to adapt to benthic life.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01837-2.

Blood Transcriptomics of Turbot Scophthalmus maximus: A Tool for Health Monitoring and Disease Studies

Simple Summary

The analysis of blood gene expression is emerging as a relevant source of information about the health status of an organism. While these investigations are increasingly performed in human and terrestrial animals, their potential is still underexplored in fish pathology. The aim of this work was to analyze the blood transcriptional profile of a commercially important flatfish species, turbot (Scophthalmus maximus), in healthy and diseased specimens. The analysis of the most expressed genes in healthy fish indicated that turbot red blood cells have important immunological functions. In diseased fish, parasitized by a myxozoan, the blood analysis reflected a broad inhibition of the immune response followed by intense inflammatory activation in heavy infections. The results showed that turbot response appears delayed, dysregulated and ineffective in stopping the infection. Particularly, a proper development of the adaptive immune response was lacking. This study points out that blood gene expression profiling is a reliable tool for health monitoring, as well as to advance in the knowledge of fish immunity and diseases.

Abstract

Blood transcriptomics is emerging as a relevant tool to monitor the status of the immune system and assist in diagnosis, prognosis, treatment and pathogenesis studies of diseases. In fish pathology, the potential of transcriptome profiling of blood is still poorly explored. Here, RNA sequencing was applied to analyze the blood transcriptional profile of turbot (Scophthalmus maximus), the most important farmed flatfish. The study was conducted in healthy specimens and specimens parasitized by the myxozoan Enteromyxum scophthalmi, which causes one of the most devastating diseases in turbot aquaculture. The blood of healthy turbot showed a transcriptomic profile mainly related to erythrocyte gas transportation function, but also to antigen processing and presentation. In moderately infected turbot, the blood reflected a broad inhibition of the immune response. Particularly, down-regulation of the B cell receptor signaling pathway was shared with heavily parasitized fish, which showed larger transcriptomic changes, including the activation of the inflammatory response. Turbot response to enteromyxosis proved to be delayed, dysregulated and ineffective in stopping the infection. The study evinces that blood transcriptomics can contribute to a better understanding of the teleost immune system and serve as a reliable tool to investigate the physiopathological status of fish.

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.

Characterization and expression analysis of mitochondrial localization molecule: NOD-like receptor X1 (Nlrx1) in mucosal tissues of turbot (Scophthalmus maximus) following bacterial challenge

The NOD-like receptor X1 (NLRX1) is a member of highly conserved nucleotide-binding domain (NBD)- and leucine-rich-repeat (LRR)-containing family (known as NLR), that localizes to the mitochondrial outer membrane and regulate the innate immunity by interacting with mitochondrial antiviral-signaling protein (MAVS). As one of cytoplasmic PRRs, NLRX1 plays key roles for pathogen recognition, autophagy and regulating of subsequent immune signaling pathways. In this study, we identified the nlrx1 in turbot as well as its expression profiles in mucosal surfaces following bacterial infection. In our results, the full-length nlrx1 transcript consists of an open reading frame (ORF) of 4,886 bp encoding the putative peptide of 966 amino acids. The phylogenetic analysis revealed the SmNlrx1 showed the closest relationship to Cynoglossus semilaevis. In addition, the Nlrx1 mRNA expression could be detected in all the examined tissues, with the most abundant expression level in head kidney, and the lowest expression level in liver. Moreover, Nlrx1 showed similar expression patterns following Vibrio anguillarum and Streptococcus iniae infection, that were both significantly up-regulated following challenge, especially post S. iniae challenge. Finally, fluorescence microscopy unveiled that the SmNlrx1 localized to mitochondria in HEK293T by N-terminal mitochondrial targeting sequence. Characterization of Nlrx1 might have an important implication in bioenergetic adaptation during metabolic stress, oncogenic transformation and innate immunity and will probably contribute to the development of novel intervention strategies for farming turbot.

Biases in Demographic Modeling Affect Our Understanding of Recent Divergence

Testing among competing demographic models of divergence has become an important component of evolutionary research in model and non-model organisms. However, the effect of unaccounted demographic events on model choice and parameter estimation remains largely unexplored. Using extensive simulations, we demonstrate that under realistic divergence scenarios, failure to account for population size (N_e) changes in daughter and ancestral populations leads to strong biases in divergence time estimates as well as model choice. We illustrate these issues reconstructing the recent demographic history of North Sea and Baltic Sea turbots (Scophthalmus maximus) by testing 16 isolation with migration (IM) and 16 secondary contact (SC) scenarios, modeling changes in N_e as well as the effects of linked selection and barrier loci. Failure to account for changes in N_e resulted in selecting SC models with long periods of strict isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for N_e changes suggest recent (<6 kya) divergence with constant gene flow. We further show how interpreting genomic landscapes of differentiation can help discerning among competing models. For example, in the turbot data, islands of differentiation show signatures of recent selective sweeps, rather than old divergence resisting secondary introgression. The results have broad implications for the study of population divergence by highlighting the potential effects of unmodeled changes in N_e on demographic inference. Tested models should aim at representing realistic divergence scenarios for the target taxa, and extreme caution should always be exercised when interpreting results of demographic modeling.

Galectins in turbot (Scophthalmus maximus L.): Characterization and expression profiling in mucosal tissues

Galectins, a family of evolutionary conserved β-galactoside-binding proteins, have been characterized in a wide range of species. Many reports have indicated vital roles of galectins in innate immunity, especially in the mucosal tissues against infection. However, the systematic identification of galectin gene family is still lacking in teleost. Here, we characterized the galectin gene family and investigated their expression profiles post bacterial challenge in turbot (Scophthalmus maximus L.). In this study, a total of 13 galectin genes were characterized in turbot, phylogenetic analyses revealed their strong relationships to half smooth tongue sole and puffer fish, and syntenic analyses confirmed the orthology suggested by the phylogenetic analysis. In addition, the copy number of galectin genes is similar across a broad spectrum of species from fish to amphibians, birds, and mammals, ranging from 8 to 16 genes. Furthermore, the galectin genes were widely expressed in all the examined turbot tissues, and most of the galectin genes were strongly expressed in mucosal tissues (skin, gill and intestine). Moreover, majority of the galectin genes were significantly regulated after Vibrio anguillarum infection in the intestine, gill and skin, suggesting that galectins were involved in the mucosal immune response to V. anguillarum infection in turbot. In addition, subcellular localization analysis showed lgals3a was distributed in the cytoplasm and nucleus. However, the knowledge of galectins are still limited in teleost species, further studies should be carried out to better characterize its detailed roles in teleost mucosal immunity.

Genome assembly and annotation at the chromosomal level of first Pleuronectidae: Verasper variegatus provides a basis for phylogenetic study of Pleuronectiformes

High quality genome is of great significance for the mining of biological information resources of species. Up to now, the genomic information of several important economic flatfishes has been well explained. All these fishes are eyes on left side-type, and no high-quality genome of eyes on right side-type species has been reported. In this study, we applied a combined strategy involving stLFR and Hi-C technologies to generate sequencing data for constructing the chromosomal genome of Verasper variegates, which belongs to Pleuronectidae with characteristic of eyes on right side. The size of genome of V. variegatus is 556 Mb. More than 97.2% of BUSCO genes were detected, and N50 lengths of the contigs and scaffolds reached 79.8 Kb and 23.8 Mb, respectively, demonstrating the outstanding completeness and sequence continuity of the genome. A total of 22,199 protein-coding genes were predicted in the assembled genome, and more than 95% of those genes could be functionally annotated. Meanwhile, the genomic collinearity, gene family and phylogenetic analyses of similar species in Pleuronectiformes were also investigated and portrayed for metamorphosis and benthic adaptation. Sex related genes mapping has also been achieved at the chromosome level. This study is the first chromosomal level genome of a Pleuronectidae fish (V. variegatus). The chromosomal genome assembly constructed in this work will not only be valuable for conservation and aquaculture studies of the V. variegatus but will also be of general interest in the phylogenetic and taxonomic studies of Pleuronectiformes.

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.

The hemoglobin Gly16β1Asp polymorphism in turbot (Scophthalmus maximus) is differentially distributed across European populations

Turbot is an important flatfish widely distributed along the European coasts, whose fishery is centered in the North Sea. The commercial value of the species has boosted a successful aquaculture sector in Europe and China. Body growth is the main target of turbot breeding programs and is also a key trait related to local adaptation to temperature and salinity. Differences in growth rate and optimal growth temperature in turbot have been shown to be associated with a hemoglobin polymorphism reported more than 50 years ago. Here, we identified a Gly16Asp amino acid substitution in the β1 globin subunit by searching for genetic variation in the five functional globin genes within the whole annotated turbot genome. We predicted increased stability of the turbot hemoglobin by the replacement of the conserved Gly with the negative charged Asp residue that is consistent with the higher rate of αβ dimer assembly in the human J-Baltimore Gly16β->Asp mutant than in normal HbA. The turbot Hbβ1-Gly16 variant dominated in the northern populations examined, particularly in the Baltic Sea, while the Asp allele showed elevated frequencies in southern populations and was the prevalent variant in the Adriatic Sea. Body weight did not associate with the Hbβ1 genotypes at farming conditions (i.e., high oxygen levels, feeding ad libitum) after analyzing 90 fish with high growth dispersal from nine turbot families. Nevertheless, all data at hand suggest that the turbot hemoglobin polymorphism has an adaptive significance in the variable wild conditions regarding temperature and oxygen availability.

An improved de novo assembling and polishing of Solea senegalensis transcriptome shed light on retinoic acid signalling in larvae

Senegalese sole is an economically important flatfish species in aquaculture and an attractive model to decipher the molecular mechanisms governing the severe transformations occurring during metamorphosis, where retinoic acid seems to play a key role in tissue remodeling. In this study, a robust sole transcriptome was envisaged by reducing the number of assembled libraries (27 out of 111 available), fine-tuning a new automated and reproducible set of workflows for de novo assembling based on several assemblers, and removing low confidence transcripts after mapping onto a sole female genome draft. From a total of 96 resulting assemblies, two "raw" transcriptomes, one containing only Illumina reads and another with Illumina and GS-FLX reads, were selected to provide SOLSEv5.0, the most informative transcriptome with low redundancy and devoid of most single-exon transcripts. It included both Illumina and GS-FLX reads and consisted of 51,348 transcripts of which 22,684 code for 17,429 different proteins described in databases, where 9527 were predicted as complete proteins. SOLSEv5.0 was used as reference for the study of retinoic acid (RA) signalling in sole larvae using drug treatments (DEAB, a RA synthesis blocker, and TTNPB, a RA-receptor agonist) for 24 and 48 h. Differential expression and functional interpretation were facilitated by an updated version of DEGenes Hunter. Acute exposure of both drugs triggered an intense, specific and transient response at 24 h but with hardly observable differences after 48 h at least in the DEAB treatments. Activation of RA signalling by TTNPB specifically increased the expression of genes in pathways related to RA degradation, retinol storage, carotenoid metabolism, homeostatic response and visual cycle, and also modified the expression of transcripts related to morphogenesis and collagen fibril organisation. In contrast, DEAB mainly decreased genes related to retinal production, impairing phototransduction signalling in the retina. A total of 755 transcripts mainly related to lipid metabolism, lipid transport and lipid homeostasis were altered in response to both treatments, indicating non-specific drug responses associated with intestinal absorption. These results indicate that a new assembling and transcript sieving were both necessary to provide a reliable transcriptome to identify the many aspects of RA action during sole development that are of relevance for sole aquaculture.

Interactions between the Parasite Philasterides dicentrarchi and the Immune System of the Turbot Scophthalmus maximus. A Transcriptomic Analysis

The present study analyses the interactions between Philasterides dicentrarchi (a ciliate parasite that causes high mortalities in cultured flatfish) and the peritoneal cells of the turbot Scophthalmus maximus during an experimental infection. The transcriptomic response was evaluated in the parasites and in the fish peritoneal cells, at 1, 2 and 4 h post-infection (hpi) in turbot injected intraperitoneally (ip) with 10⁷ ciliates and at 12 and 48 hpi in turbot injected ip with 10⁵ ciliates. Numerous genes were differentially expressed (DE) in P. dicentrarchi, relative to their expression in control ciliates (0 hpi): 407 (369 were up-regulated) at 1 hpi, 769 (415 were up-regulated) at 2 hpi and 507 (119 were up-regulated) at 4 hpi. Gene ontology (GO) analysis of the DE genes showed that the most representative categories of biological processes affected at 1, 2 and 4 hpi were biosynthetic processes, catabolic processes, biogenesis, proteolysis and transmembrane transport. Twelve genes of the ABC transporter family and eight genes of the leishmanolysin family were DE at 1, 2 and 4 hpi. Most of these genes were strongly up-regulated (UR), suggesting that they are involved in P. dicentrarchi infection. A third group of UR genes included several genes related to ribosome biogenesis, DNA transcription and RNA translation. However, expression of tubulins and tubulin associated proteins, such as kinesins or dyneins, which play key roles in ciliate division and movement, was down-regulated (DR). Similarly, genes that coded for lysosomal proteins or that participate in the cell cycle mitotic control, glycolysis, the Krebs cycle and/or in the electron transport chain were also DR. The transcriptomic analysis also revealed that in contrast to many parasites, which passively evade the host immune system, P. dicentrarchi strongly stimulated turbot peritoneal cells. Many genes related to inflammation were DE in peritoneal cells at 1, 2 and 4 hpi. However, the response was much lower at 12 hpi and almost disappeared completely at 48 hpi in fish that were able to kill P. dicentrarchi during the first few hpi. The genes that were DE at 1, 2 and 4 hpi were mainly related to the apoptotic process, the immune response, the Fc-epsilon receptor signalling pathway, the innate immune response, cell adhesion, cell surface receptors, the NF-kappaB signalling pathway and the MAPK cascade. Expression of toll-like receptors 2, 5 and 13 and of several components of NF-κB, MAPK and JAK/STAT signalling pathways was UR in the turbot peritoneal cells. Genes expressing chemokines and chemokine receptors, genes involved in prostaglandin and leukotriene synthesis, prostaglandins, leukotriene receptors, proinflammatory cytokines and genes involved in apoptosis were strongly UR during the first four hours of infection. However, expression of anti-inflammatory cytokines such as Il-10 and lipoxygenases with anti-inflammatory activity (i.e., arachidonate 15-lipoxygenase) were only UR at 12 and/or 48 hpi, indicating an anti-inflammatory state in these groups of fish. In conclusion, the present study shows the regulation of several genes in P. dicentrarchi during the early stages of infection, some of which probably play important roles in this process. The infection induced a potent acute inflammatory response, and many inflammatory genes were regulated in peritoneal cells, showing that the turbot uses all the protective mechanisms it has available to prevent the entry of the parasite.

Genetic analysis redraws the management boundaries for the European sprat

Sustainable fisheries management requires detailed knowledge of population genetic structure. The European sprat is an important commercial fish distributed from Morocco to the Arctic circle, Baltic, Mediterranean, and Black seas. Prior to 2018, annual catch advice on sprat from the International Council for the Exploration of the Sea (ICES) was based on five putative stocks: (a) North Sea, (b) Kattegat-Skagerrak and Norwegian fjords, (c) Baltic Sea, (d) West of Scotland-southern Celtic Seas, and (e) English Channel. However, there were concerns that the sprat advice on stock size estimates management plan inadequately reflected the underlying biological units. Here, we used ddRAD sequencing to develop 91 SNPs that were thereafter used to genotype approximately 2,500 fish from 40 locations. Three highly distinct and relatively homogenous genetic groups were identified: (a) Norwegian fjords; (b) Northeast Atlantic including the North Sea, Kattegat-Skagerrak, Celtic Sea, and Bay of Biscay; and (c) Baltic Sea. Evidence of genetic admixture and possibly physical mixing was detected in samples collected from the transition zone between the North and Baltic seas, but not between any of the other groups. These results have already been implemented by ICES with the decision to merge the North Sea and the Kattegat-Skagerrak sprat to be assessed as a single unit, thus demonstrating that genetic data can be rapidly absorbed to align harvest regimes and biological units.

The Teleost Thymus in Health and Disease: New Insights from Transcriptomic and Histopathological Analyses of Turbot, Scophthalmus maximus

The thymus is a primary lymphoid organ that plays a pivotal role in the adaptive immune system. The immunobiology of the thymus in fish is considered to be similar to that of mammals, but it is actually poorly characterized in several cultured teleost species. In particular, while investigations in human and veterinary medicine have highlighted that the thymus can be affected by different pathological conditions, little is known about its response during disease in fish. To better understand the role of the thymus under physiological and pathological conditions, we conducted a study in turbot (Scophthalmus maximus), a commercially valuable flatfish species, combining transcriptomic and histopathological analyses. The myxozoan parasite Enteromyxum scophthalmi, which represents a major challenge to turbot production, was used as a model of infection. The thymus tissues of healthy fish showed overrepresented functions related to its immunological role in T-cell development and maturation. Large differences were observed between the transcriptomes of control and severely infected fish. Evidence of inflammatory response, apoptosis modulation, and declined thymic function associated with loss of cellularity was revealed by both genomic and morphopathological analyses. This study presents the first description of the turbot thymus transcriptome and provides novel insights into the role of this organ in teleosts’ immune responses.

Characterization of a novel lncRNA (SETD3-OT) in turbot (Scophthalmus maximus L.)

LncRNAs have been demonstrated to play pivotal roles in virous biological processes, especially the gene expression regulation, including transcriptional regulation, posttranscriptional control and epigenetic processes. However, most of the current studies of lncRNAs are still limited in mammalian species, the investigations of functional roles of lncRNAs in teleost species are still lacking. In current study, we identified a novel lncRNA (SETD3-OT) in turbot, with 2,504 bp full-length obtained by 5' and 3' RACE, located in turbot chromosome 17, ranged from 20,933,835 to 20,936,302 bp. In addition, 8 neighboring genes of SETD3-OT were identified within 100 kbp in genome location. From the annotation of the neighboring adjacent genes, SETD3-OT might involve in regulation of cell apoptosis and cycle, the immune cell development, and the immune response against infection, and its expression pattern is similar to majority of the neighboring genes following Aeromonas salmonicida challenge. Intriguingly, SETD3-OT showed significant high expression levels in mucosal surfaces (intestine, gill and skin), and was dramatically down-regulated in these mucosal tissues following Vibrio anguillarum challenge, especially in gill and skin. In addition, SETD3-OT was distributed in nucleus, it might regulate the neighboring genes in cis or in trans. Taken together, our results provide insights for lncRNA in fish innate immunity, further studies should be conduct to explore the detailed molecular mechanism of the gene regulation between SETD3-OT and its neighboring genes.

Genome-wide identification, expression signature and immune functional analysis of two cathepsin S (CTSS) genes in turbot (Scophthalmus maximus L.)

Cathepsins, a superfamily of hydrolytic enzymes produced and enclosed within lysosomes, play multiple roles at physiological and pathological states. Cathepsin S is a lysosomal cysteine endopeptidase of the papain family, and exerts critical roles in the regulation of MHC class II immune responses. In the present study, we captured two Cathepsin S genes in turbot (SmCTSS1 and SmCTSS2.1), characterized their expression patterns following V. anguillarum and S. iniae infections, and explored their binding ability and agglutination capability. Firstly, the SmCTSS1 contained a 990 bp ORF encoding 329 amino acids, while SmCTSS2.1 contained a 1,014 bp ORF encoding 337 amino acids. The phylogenetic analysis revealed that both genes showed the closest relationship to their counterparts of Japanese flounder (Paralichthys olivaceus). In addition, both genes were ubiquitously expressed in all examined healthy tissues, with the highest expression level observed in spleen and intestine, respectively, while the lowest expression level both observed in liver. Both SmCTSS1 and SmCTSS2.1 were significantly differentially expressed, and exhibited general down-regulations at most time points in skin and intestine after two bacterial infections. Finally, both rSmCTSS1 and rSmCTSS2.1 showed significant binding ability to three examined microbial ligands (LPS, PGN and LTA), and strong agglutination effect to different bacteria (E. tarda, S. agalactiae, S. aureus and V. anguillarum). Collectively, this study provided valuable data for understanding the roles of CTSS in the host defense against bacterial infections in turbot, and indicated the potential vital roles of CTSS in innate immune responses of teleost species.