Transcriptome analysis of flathead grey mullet (Mugil cephalus) ovarian development induced by recombinant gonadotropin hormones.

Background: Treatment with recombinant gonadotropin hormones (rGths), follicle-stimulating hormone (rFsh) and luteinizing hormone (rLh), was shown to induce and complete vitellogenesis to finally obtain viable eggs and larvae in the flathead grey mullet (Mugil cephalus), a teleost arrested at early stages of gametogenesis in intensive captivity conditions. This study aimed to investigate the transcriptomic changes that occur in the ovary of females during the rGths-induced vitellogenesis. Methods: Ovarian samples were collected through biopsies from the same five females at four stages of ovarian development. RNASeq libraries were constructed for all stages studied, sequenced on an Illumina HiSeq4000, and a de novo transcriptome was constructed. Differentially expressed genes (DEGs) were identified between stages and the functional properties of DEGs were characterized by comparison with the gene ontology and Kyoto Encyclopedia. An enrichment analysis of molecular pathways was performed. Results: The de novo transcriptome comprised 287,089 transcripts after filtering. As vitellogenesis progressed, more genes were significantly upregulated than downregulated. The rFsh application induced ovarian development from previtellogenesis to early-to-mid-vitellogenesis with associated pathways enriched from upregulated DEGs related to ovarian steroidogenesis and reproductive development, cholesterol metabolism, ovarian growth and differentiation, lipid accumulation, and cell-to-cell adhesion pathways. The application of rFsh and rLh at early-to-mid-vitellogenesis induced the growth of oocytes to late-vitellogenesis and, with it, the enrichment of pathways from upregulated DEGs related to the production of energy, such as the lysosomes activity. The application of rLh at late-vitellogenesis induced the completion of vitellogenesis with the enrichment of pathways linked with the switch from vitellogenesis to oocyte maturation. Conclusion: The DEGs and enriched molecular pathways described during the induced vitellogenesis of flathead grey mullet with rGths were typical of natural oogenesis reported for other fish species. Present results add new knowledge to the rGths action to further raise the possibility of using rGths in species that present similar reproductive disorders in aquaculture, the aquarium industry as well as the conservation of endangered species.

Chromosome genome assembly for the meagre, Argyrosomus regius, reveals species adaptations and sciaenid sex-related locus evolution.

The meagre, Argyrosomus regius, has recently become a species of increasing economic interest for the Mediterranean aquaculture and there is ongoing work to boost production efficiency through selective breeding. Access to the complete genomic sequence will provide an essential resource for studying quantitative trait-associated loci and exploring the genetic diversity of different wild populations and aquaculture stocks in more detail. Here, we present the first complete genome for A. regius, produced through a combination of long and short read technologies and an efficient in-house developed pipeline for assembly and polishing. Scaffolding using previous linkage map data allowed us to reconstruct a chromosome level assembly with high completeness, complemented with gene annotation and repeat masking. The 696 Mb long assembly has an N50 = 27.87 Mb and an L50 = 12, with 92.85% of its length placed in 24 chromosomes. We use this new resource to study the evolution of the meagre genome and other Sciaenids, via a comparative analysis of 25 high-quality teleost genomes. Combining a rigorous investigation of gene duplications with base-wise conservation analysis, we identify candidate loci related to immune, fat metabolism and growth adaptations in the meagre. Following phylogenomic reconstruction, we show highly conserved synteny within Sciaenidae. In contrast, we report rapidly evolving syntenic rearrangements and gene copy changes in the sex-related dmrt1 neighbourhood in meagre and other members of the family. These novel genomic datasets and findings will add important new tools for aquaculture studies and greatly facilitate husbandry and breeding work in the species.

Genome Analysis of Lagocephalus sceleratus: Unraveling the Genomic Landscape of a Successful Invader.

The Tetraodontidae family encompasses several species which attract scientific interest in terms of their ecology and evolution. The silver-cheeked toadfish (Lagocephalus sceleratus) is a well-known "invasive sprinter" that has invaded and spread, in less than a decade, throughout the Eastern and part of the Western Mediterranean Sea from the Red Sea through the Suez Canal. In this study, we built and analysed the first near-chromosome level genome assembly of L. sceleratus and explored its evolutionary landscape. Through a phylogenomic analysis, we positioned L. sceleratus closer to T. nigroviridis, compared to other members of the family, while gene family evolution analysis revealed that genes associated with the immune response have experienced rapid expansion, providing a genetic basis for studying how L. sceleratus is able to achieve highly successful colonisation. Moreover, we found that voltage-gated sodium channel (NaV 1.4) mutations previously connected to tetrodotoxin resistance in other pufferfishes are not found in L. sceleratus, highlighting the complex evolution of this trait. The high-quality genome assembly built here is expected to set the ground for future studies on the species biology.

Community Parameters and Genome-Wide RAD-Seq Loci of Ceratothoa oestroides Imply Its Transfer between Farmed European Sea Bass and Wild Farm-Aggregating Fish.

Wild fish assemblages that aggregate within commercial marine aquaculture sites for feeding and shelter have been considered as a primary source of pathogenic parasites vectored to farmed fish maintained in net pens at an elevated density. In order to evaluate whether Ceratothoa oestroides (Isopoda, Cymothoidae), a generalist and pestilent isopod that is frequently found in Adriatic and Greek stocks of farmed European sea bass (Dicentrarchus labrax), transfers between wild and farmed fish, a RAD-Seq (restriction-site-associated DNA sequencing)-mediated genetic screening approach was employed. The double-digest RAD-Seq of 310 C. oestroides specimens collected from farmed European sea bass (138) and different wild farm-aggregating fish (172) identified 313 robust SNPs that evidenced a close genetic relatedness between the "wild" and "farmed" genotypes. ddRAD-Seq proved to be an effective method for detecting the discrete genetic structuring of C. oestroides and genotype intermixing between two populations. The parasite prevalence in the farmed sea bass was 1.02%, with a mean intensity of 2.0 and mean abundance of 0.02, while in the wild fish, the prevalence was 8.1%; the mean intensity, 1.81; and the mean abundance, 0.15. Such differences are likely a consequence of human interventions during the farmed fish's rearing cycle that, nevertheless, did not affect the transfer of C. oestroides.

Phylogenomics investigation of sparids (Teleostei: Spariformes) using high-quality proteomes highlights the importance of taxon sampling.

Sparidae (Teleostei: Spariformes) are a family of fish constituted by approximately 150 species with high popularity and commercial value, such as porgies and seabreams. Although the phylogeny of this family has been investigated multiple times, its position among other teleost groups remains ambiguous. Most studies have used a single or few genes to decipher the phylogenetic relationships of sparids. Here, we conducted a thorough phylogenomic analysis using five recently available Sparidae gene-sets and 26 high-quality, genome-predicted teleost proteomes. Our analysis suggested that Tetraodontiformes (puffer fish, sunfish) are the closest relatives to sparids than all other groups used. By analytically comparing this result to our own previous contradicting finding, we show that this discordance is not due to different orthology assignment algorithms; on the contrary, we prove that it is caused by the increased taxon sampling of the present study, outlining the great importance of this aspect in phylogenomic analyses in general.

Scanning of Genetic Variants and Genetic Mapping of Phenotypic Traits in Gilthead Sea Bream Through ddRAD Sequencing.

Gilthead sea bream (Sparus aurata) is a teleost of considerable economic importance in Southern European aquaculture. The aquaculture industry shows a growing interest in the application of genetic methods that can locate phenotype-genotype associations with high economic impact. Through selective breeding, the aquaculture industry can exploit this information to maximize the financial yield. Here, we present a Genome Wide Association Study (GWAS) of 112 samples belonging to seven different sea bream families collected from a Greek commercial aquaculture company. Through double digest Random Amplified DNA (ddRAD) Sequencing, we generated a per-sample genetic profile consisting of 2,258 high-quality Single Nucleotide Polymorphisms (SNPs). These profiles were tested for association with four phenotypes of major financial importance: Fat, Weight, Tag Weight, and the Length to Width ratio. We applied two methods of association analysis. The first is the typical single-SNP to phenotype test, and the second is a feature selection (FS) method through two novel algorithms that are employed for the first time in aquaculture genomics and produce groups with multiple SNPs associated to a phenotype. In total, we identified 9 single SNPs and 6 groups of SNPs associated with weight-related phenotypes (Weight and Tag Weight), 2 groups associated with Fat, and 16 groups associated with the Length to Width ratio. Six identified loci (Chr4:23265532, Chr6:12617755, Chr:8:11613979, Chr13:1098152, Chr15:3260819, and Chr22:14483563) were present in genes associated with growth in other teleosts or even mammals, such as semaphorin-3A and neurotrophin-3. These loci are strong candidates for future studies that will help us unveil the genetic mechanisms underlying growth and improve the sea bream aquaculture productivity by providing genomic anchors for selection programs.

Genomic analysis of Sparus aurata reveals the evolutionary dynamics of sex-biased genes in a sequential hermaphrodite fish.

Sexual dimorphism is a fascinating subject in evolutionary biology and mostly results from sex-biased expression of genes, which have been shown to evolve faster in gonochoristic species. We report here genome and sex-specific transcriptome sequencing of Sparus aurata, a sequential hermaphrodite fish. Evolutionary comparative analysis reveals that sex-biased genes in S. aurata are similar in number and function, but evolved following strikingly divergent patterns compared with gonochoristic species, showing overall slower rates because of stronger functional constraints. Fast evolution is observed only for highly ovary-biased genes due to female-specific patterns of selection that are related to the peculiar reproduction mode of S. aurata, first maturing as male, then as female. To our knowledge, these findings represent the first genome-wide analysis on sex-biased loci in a hermaphrodite vertebrate species, demonstrating how having two sexes in the same individual profoundly affects the fate of a large set of evolutionarily relevant genes.

The Gene Toolkit Implicated in Functional Sex in Sparidae Hermaphrodites: Inferences From Comparative Transcriptomics.

Sex-biased gene expression is the mode through which sex dimorphism arises from a nearly identical genome, especially in organisms without genetic sex determination. Teleost fishes show great variations in the way the sex phenotype forms. Among them, Sparidae, that might be considered as a model family displays a remarkable diversity of reproductive modes. In this study, we sequenced and analyzed the sex-biased transcriptome in gonads and brain (the tissues with the most profound role in sexual development and reproduction) of two sparids with different reproductive modes: the gonochoristic common dentex, Dentex dentex, and the protandrous hermaphrodite gilthead seabream, Sparus aurata. Through comparative analysis with other protogynous and rudimentary protandrous sparid transcriptomes already available, we put forward common male and female-specific genes and pathways that are probably implicated in sex-maintenance in this fish family. Our results contribute to the understanding of the complex processes behind the establishment of the functional sex, especially in hermaphrodite species and set the groundwork for future experiments by providing a gene toolkit that can improve efforts to control phenotypic sex in finfish in the ever-increasingly important field of aquaculture.

Exploring a Nonmodel Teleost Genome Through RAD Sequencing-Linkage Mapping in Common Pandora, Pagellus erythrinus and Comparative Genomic Analysis.

Common pandora (Pagellus erythrinus) is a benthopelagic marine fish belonging to the teleost family Sparidae, and a newly recruited species in Mediterranean aquaculture. The paucity of genetic information relating to sparids, despite their growing economic value for aquaculture, provides the impetus for exploring the genomics of this fish group. Genomic tool development, such as genetic linkage maps provision, lays the groundwork for linking genotype to phenotype, allowing fine-mapping of loci responsible for beneficial traits. In this study, we applied ddRAD methodology to identify polymorphic markers in a full-sib family of common pandora. Employing the Illumina MiSeq platform, we sampled and sequenced a size-selected genomic fraction of 99 individuals, which led to the identification of 920 polymorphic loci. Downstream mapping analysis resulted in the construction of 24 robust linkage groups, corresponding to the karyotype of the species. The common pandora linkage map showed varying degrees of conserved synteny with four other teleost genomes, namely the European seabass (Dicentrarchus labrax), Nile tilapia (Oreochromis niloticus), stickleback (Gasterosteus aculeatus), and medaka (Oryzias latipes), suggesting a conserved genomic evolution in Sparidae. Our work exploits the possibilities of genotyping by sequencing to gain novel insights into genome structure and evolution. Such information will boost the study of cultured species and will set the foundation for a deeper understanding of the complex evolutionary history of teleosts.

The sex-specific transcriptome of the hermaphrodite sparid sharpsnout seabream (Diplodus puntazzo).

BACKGROUND: Teleosts are characterized by a remarkable breadth of sexual mechanisms including various forms of hermaphroditism. Sparidae is a fish family exhibiting gonochorism or hermaphroditism even in closely related species. The sparid Diplodus puntazzo (sharpsnout seabream), exhibits rudimentary hermaphroditism characterized by intersexual immature gonads but single-sex mature ones. Apart from the intriguing reproductive biology, it is economically important with a continuously growing aquaculture in the Mediterranean Sea, but limited available genetic resources. Our aim was to characterize the expressed transcriptome of gonads and brains through RNA-Sequencing and explore the properties of genes that exhibit sex-biased expression profiles. RESULTS: Through RNA-Sequencing we obtained an assembled transcriptome of 82,331 loci. The expression analysis uncovered remarkable differences between male and female gonads, while male and female brains were almost identical. Focused search for known targets of sex determination and differentiation in vertebrates built the sex-specific expression profile of sharpsnout seabream. Finally, a thorough genetic marker discovery pipeline led to the retrieval of 85,189 SNPs and 29,076 microsatellites enriching the available genetic markers for this species. CONCLUSIONS: We obtained a nearly complete source of transcriptomic sequence as well as marker information for sharpsnout seabream, laying the ground for understanding the complex process of sex differentiation of this economically valuable species. The genes involved include known candidates from other vertebrate species, suggesting a conservation of the toolkit between gonochorists and hermaphrodites.