Extensive Loss and Gain of Conserved Noncoding Elements During Early Teleost Evolution.

Conserved noncoding elements in vertebrates are enriched around transcription factor loci associated with development. However, loss and rapid divergence of conserved noncoding elements has been reported in teleost fish, albeit taking only few genomes into consideration. Taking advantage of the recent increase in high-quality teleost genomes, we focus on studying the evolution of teleost conserved noncoding elements, carrying out targeted genomic alignments and comparisons within the teleost phylogeny to detect conserved noncoding elements and reconstruct the ancestral teleost conserved noncoding elements repertoire. This teleost-centric approach confirms previous observations of extensive vertebrate conserved noncoding elements loss early in teleost evolution, but also reveals massive conserved noncoding elements gain in the teleost stem-group over 300 million years ago. Using synteny-based association to link conserved noncoding elements to their putatively regulated target genes, we show the most teleost gained conserved noncoding elements are found in the vicinity of orthologous loci involved in transcriptional regulation and embryonic development that are also associated with conserved noncoding elements in other vertebrates. Moreover, teleost and vertebrate conserved noncoding elements share a highly similar motif and transcription factor binding site vocabulary. We suggest that early teleost conserved noncoding element gains reflect a restructuring of the ancestral conserved noncoding element repertoire through both extreme divergence and de novo emergence. Finally, we support newly identified pan-teleost conserved noncoding elements have potential for accurate resolution of teleost phylogenetic placements in par with coding sequences, unlike ancestral only elements shared with spotted gar. This work provides new insight into conserved noncoding element evolution with great value for follow-up work on phylogenomics, comparative genomics, and the study of gene regulation evolution in teleosts.

Sardines at a junction: Seascape genomics reveals ecological and oceanographic drivers of variation in the NW Mediterranean Sea.

By evaluating genetic variation across the entire genome, one can address existing questions in a novel way while raising new ones. The latter includes how different local environments influence adaptive and neutral genomic variation within and among populations, providing insights into local adaptation of natural populations and their responses to global change. Here, under a seascape genomic approach, ddRAD data of 4609 single nucleotide polymorphisms (SNPs) from 398 sardines (Sardina pilchardus) collected in 11 Mediterranean and one Atlantic site were generated. These were used along with oceanographic and ecological information to detect signals of adaptive divergence with gene flow across environmental gradients. The studied sardines constitute two clusters (FST  = 0.07), a pattern attributed to outlier loci, highlighting putative local adaptation. The trend in the number of days with sea surface temperature above 19°C, a critical threshold for successful sardine spawning, was crucial at all levels of population structuring with implications on the species' key biological processes. Outliers link candidate SNPs to the region's environmental heterogeneity. Our findings provide evidence for a dynamic equilibrium in which population structure is maintained by physical and ecological factors under the opposing influences of migration and selection. This dynamic in a natural system warrants continuous monitoring under a seascape genomic approach that might benefit from a temporal and more detailed spatial dimension. Our results may contribute to complementary studies aimed at providing deeper insights into the mechanistic processes underlying population structuring. Those are key to understanding and predicting future changes and responses of this highly exploited species in the face of climate change.

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.

SnakeCube: containerized and automated pipeline for de novo genome assembly in HPC environments.

OBJECTIVE: The rapid progress in sequencing technology and related bioinformatics tools aims at disentangling diversity and conservation issues through genome analyses. The foremost challenges of the field involve coping with questions emerging from the swift development and application of new algorithms, as well as the establishment of standardized analysis approaches that promote transparency and transferability in research. RESULTS: Here, we present SnakeCube, an automated and containerized whole de novo genome assembly pipeline that runs within isolated, secured environments and scales for use in High Performance Computing (HPC) domains. SnakeCube was optimized for its performance and tested for its effectiveness with various inputs, highlighting its safe and robust universal use in the field.

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.

0s and 1s in marine molecular research: a regional HPC perspective.

High-performance computing (HPC) systems have become indispensable for modern marine research, providing support to an increasing number and diversity of users. Pairing with the impetus offered by high-throughput methods to key areas such as non-model organism studies, their operation continuously evolves to meet the corresponding computational challenges. Here, we present a Tier 2 (regional) HPC facility, operating for over a decade at the Institute of Marine Biology, Biotechnology, and Aquaculture of the Hellenic Centre for Marine Research in Greece. Strategic choices made in design and upgrades aimed to strike a balance between depth (the need for a few high-memory nodes) and breadth (a number of slimmer nodes), as dictated by the idiosyncrasy of the supported research. Qualitative computational requirement analysis of the latter revealed the diversity of marine fields, methods, and approaches adopted to translate data into knowledge. In addition, hardware and software architectures, usage statistics, policy, and user management aspects of the facility are presented. Drawing upon the last decade's experience from the different levels of operation of the Institute of Marine Biology, Biotechnology, and Aquaculture HPC facility, a number of lessons are presented; these have contributed to the facility's future directions in light of emerging distribution technologies (e.g., containers) and Research Infrastructure evolution. In combination with detailed knowledge of the facility usage and its upcoming upgrade, future collaborations in marine research and beyond are envisioned.

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.

An evolutionary perspective on marine invasions.

Species distributions are rapidly changing as human globalization increasingly moves organisms to novel environments. In marine systems, species introductions are the result of a number of anthropogenic mechanisms, notably shipping, aquaculture/mariculture, the pet and bait trades, and the creation of canals. Marine invasions are a global threat to human and non-human populations alike and are often listed as one of the top conservation concerns worldwide, having ecological, evolutionary, and social ramifications. Evolutionary investigations of marine invasions can provide crucial insight into an introduced species' potential impacts in its new range, including: physiological adaptation and behavioral changes to exploit new environments; changes in resident populations, community interactions, and ecosystems; and severe reductions in genetic diversity that may limit evolutionary potential in the introduced range. This special issue focuses on current research advances in the evolutionary biology of marine invasions and can be broadly classified into a few major avenues of research: the evolutionary history of invasive populations, post-invasion reproductive changes, and the role of evolution in parasite introductions. Together, they demonstrate the value of investigating marine invasions from an evolutionary perspective, with benefits to both fundamental and applied evolutionary biology at local and broad scales.

A de novo transcriptome assembly for the bath sponge Spongia officinalis, adjusting for microsymbionts.

OBJECTIVES: We report a transcriptome acquisition for the bath sponge Spongia officinalis, a non-model marine organism that hosts rich symbiotic microbial communities. To this end, a pipeline was developed to efficiently separate between bacterial expressed genes from those of eukaryotic origin. The transcriptome was produced to support the assessment of gene expression and, thus, the response of the sponge, to elevated temperatures, replicating conditions currently occurring in its native habitat. DATA DESCRIPTION: We describe the assembled transcriptome along with the bioinformatic pipeline used to discriminate between signals of metazoan and prokaryotic origin. The pipeline involves standard read pre-processing steps and incorporates extra analyses to identify and filter prokaryotic reads out of the analysis. The proposed pipeline can be followed to overcome the technical RNASeq problems characteristic for symbiont-rich metazoan organisms with low or non-existent tissue differentiation, such as sponges and cnidarians. At the same time, it can be valuable towards the development of approaches for parallel transcriptomic studies of symbiotic communities and the host.

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.

Lateralized Feeding Behavior is Associated with Asymmetrical Neuroanatomy and Lateralized Gene Expressions in the Brain in Scale-Eating Cichlid Fish.

Lateralized behavior ("handedness") is unusual, but consistently found across diverse animal lineages, including humans. It is thought to reflect brain anatomical and/or functional asymmetries, but its neuro-molecular mechanisms remain largely unknown. Lake Tanganyika scale-eating cichlid fish, Perissodus microlepis show pronounced asymmetry in their jaw morphology as well as handedness in feeding behavior-biting scales preferentially only from one or the other side of their victims. This makes them an ideal model in which to investigate potential laterality in neuroanatomy and transcription in the brain in relation to behavioral handedness. After determining behavioral handedness in P. microlepis (preferred attack side), we estimated the volume of the hemispheres of brain regions and captured their gene expression profiles. Our analyses revealed that the degree of behavioral handedness is mirrored at the level of neuroanatomical asymmetry, particularly in the tectum opticum. Transcriptome analyses showed that different brain regions (tectum opticum, telencephalon, hypothalamus, and cerebellum) display distinct expression patterns, potentially reflecting their developmental interrelationships. For numerous genes in each brain region, their extent of expression differences between hemispheres was found to be correlated with the degree of behavioral lateralization. Interestingly, the tectum opticum and telencephalon showed divergent biases on the direction of up- or down-regulation of the laterality candidate genes (e.g., grm2) in the hemispheres, highlighting the connection of handedness with gene expression profiles and the different roles of these brain regions. Hence, handedness in predation behavior may be caused by asymmetric size of brain hemispheres and also by lateralized gene expressions in the brain.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons.

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

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.

Transcriptomics of two evolutionary novelties: how to make a sperm-transfer organ out of an anal fin and a sexually selected "sword" out of a caudal fin.

Swords are exaggerated male ornaments of swordtail fishes that have been of great interest to evolutionary biologists ever since Darwin described them in the Descent of Man (1871). They are a novel sexually selected trait derived from modified ventral caudal fin rays and are only found in the genus Xiphophorus. Another phylogenetically more widespread and older male trait is the gonopodium, an intromittent organ found in all poeciliid fishes, that is derived from a modified anal fin. Despite many evolutionary and behavioral studies on both traits, little is known so far about the molecular mechanisms underlying their development. By investigating transcriptomic changes (utilizing a RNA-Seq approach) in response to testosterone treatment in the swordtail fish, Xiphophorus hellerii, we aimed to better understand the architecture of the gene regulatory networks underpinning the development of these two evolutionary novelties. Large numbers of genes with tissue-specific expression patterns were identified. Among the "sword genes" those involved in embryonic organ development, sexual character development and coloration were highly expressed, while in the gonopodium rather more morphogenesis-related genes were found. Interestingly, many genes and genetic pathways are shared between both developing novel traits derived from median fins: the sword and the gonopodium. Our analyses show that a larger set of gene networks was co-opted during the development and evolution of the "older" gonopodium than in the "younger," and morphologically less complex trait, the sword. We provide a catalog of candidate genes for future efforts to dissect the development of those sexually selected exaggerated male traits in swordtails.

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.