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.

Microbial community differentiation between active and inactive sulfide chimneys of the Kolumbo submarine volcano, Hellenic Volcanic Arc.

Over the last decades, there has been growing interest about the ecological role of hydrothermal sulfide chimneys, their microbial diversity and associated biotechnological potential. Here, we performed dual-index Illumina sequencing of bacterial and archaeal communities on active and inactive sulfide chimneys collected from the Kolumbo hydrothermal field, situated on a geodynamic convergent setting. A total of 15,701 OTUs (operational taxonomic units) were assigned to 56 bacterial and 3 archaeal phyla, 133 bacterial and 16 archaeal classes. Active chimney communities were dominated by OTUs related to thermophilic members of Epsilonproteobacteria, Aquificae and Deltaproteobacteria. Inactive chimney communities were dominated by an OTU closely related to the archaeon Nitrosopumilus sp., and by members of Gammaproteobacteria, Deltaproteobacteria, Planctomycetes and Bacteroidetes. These lineages are closely related to phylotypes typically involved in iron, sulfur, nitrogen, hydrogen and methane cycling. Overall, the inactive sulfide chimneys presented highly diverse and uniform microbial communities, in contrast to the active chimney communities, which were dominated by chemolithoautotrophic and thermophilic lineages. This study represents one of the most comprehensive investigations of microbial diversity in submarine chimneys and elucidates how the dissipation of hydrothermal activity affects the structure of microbial consortia in these extreme ecological niches.

Correction to: Microbial community differentiation between active and inactive sulfide chimneys of the Kolumbo submarine volcano, Hellenic Volcanic Arc.

In the original publication there is a mistake in the supplementary material. The correct supplementary material is provided in this correction article.

Metagenomic investigation of the geologically unique Hellenic Volcanic Arc reveals a distinctive ecosystem with unexpected physiology.

Hydrothermal vents represent a deep, hot, aphotic biosphere where chemosynthetic primary producers, fuelled by chemicals from Earth's subsurface, form the basis of life. In this study, we examined microbial mats from two distinct volcanic sites within the Hellenic Volcanic Arc (HVA). The HVA is geologically and ecologically unique, with reported emissions of CO2 -saturated fluids at temperatures up to 220°C and a notable absence of macrofauna. Metagenomic data reveals highly complex prokaryotic communities composed of chemolithoautotrophs, some methanotrophs, and to our surprise, heterotrophs capable of anaerobic degradation of aromatic hydrocarbons. Our data suggest that aromatic hydrocarbons may indeed be a significant source of carbon in these sites, and instigate additional research into the nature and origin of these compounds in the HVA. Novel physiology was assigned to several uncultured prokaryotic lineages; most notably, a SAR406 representative is attributed with a role in anaerobic hydrocarbon degradation. This dataset, the largest to date from submarine volcanic ecosystems, constitutes a significant resource of novel genes and pathways with potential biotechnological applications.

The bioinformatics landscape in environmental omics: Lessons from a national ELIXIR survey.

As a research infrastructure with a mission to provide services for bioinformatics, ELIXIR aims to identify and inform its target audiences. Here, we present a survey on a community of researchers studying the environment with omics approaches in Greece, one of the youngest member countries of ELIXIR. Personal interviews followed by quantitative and qualitative analysis were employed to document interactions and practices of the community and to perform a gap analysis for the transition toward multiomics and systems biology. Environmental omics in Greece mostly concerns production of data, in large majority on microbes and non-model organisms. Our survey highlighted (1) the popularity and suitability of targeted hands-on training events; (2) data quality and management issues as important elements for the transition to multiomics, and (3) lack of knowledge and misconceptions regarding interoperability, metadata standards, and pre-registration. The publicly available collected answers represent a valuable resource in view of future strategic planning.

Microbial diversity in four Mediterranean irciniid sponges.

This paper describes a dataset of microbial communities from four different sponge species: Irciniaoros (Schmidt, 1864), Irciniavariabilis (Schmidt, 1862), Sarcotragusspinosulus Schmidt, 1862 and Sarcotragusfasciculatus (Pallas, 1766). The examined sponges all belong to Demospongiae (Class); Keratosa (Subclass); Dictyoceratida (Order); Irciniidae (Family). Samples were collected by scuba diving at depths between 6-14 m from two sampling sites of rocky formations at the northern coast of Crete (Cretan Sea, eastern Mediterranean) and were subjected to metabarcoding for the V5-V6 region of the 16S rRNA gene.

Assessing the genetic landscape of a contact zone: the case of European hare in northeastern Greece.

The European hare populations of the Balkan Peninsula comprise two divergent phylogenetic lineages with discrete geographical distribution slightly overlapping in the area of northeastern Greece and Bulgaria. Here we elucidate their contact zone, by defining the spatial distributional pattern of the two highly divergent groups, detecting individuals of hybrid origin, and identifying genetic barriers present in the area of their co-existence. Specimens from northeastern Greece were assayed for lineage assignment and population genetic inference based on a 511 bp fragment of mitochondrial DNA control region and allelic data from 10 microsatellite loci. Bayesian analyses on original and simulated genotypes were performed allowing for the contact zone delineation. Our results indicate high genetic diversity in both nuclear and mitochondrial DNA, strong population structure and non random spatial distribution of the differentiated gene pools. The information provided by the two types of molecular markers yielded consistent results. This study comprises a fine scale analysis of the contact zone between the two evolutionary lineages of European brown hares in northeastern Greece. Specific questions on the spatial patterns where addressed for the first time. Furthermore, hypotheses regarding the presence of hybrids were also tested. As a result, interpretive power to the diversity patterns observed today in the Balkans was added and previously overlooked aspects of the species biology were highlighted.

Taxonomic diversity of microbial communities in sub-seafloor hydrothermal sediments of the active Santorini-Kolumbo volcanic field.

Introduction: Active hydrothermal vents of volcanic origin provide a remarkable manifestation of life on Earth under extreme conditions, which may have consequences for our understanding of habitability on other terrestrial bodies as well. Methods: Here, we performed for the first time Illumina sequencing of bacterial and archaeal communities on sub-seafloor samples collected from the Santorini-Kolumbo volcanic field. A total of 19 (3-m long) gravity corers were collected and processed for microbial community analysis. Results: From a total of 6,46,671 produced V4 sequences for all samples, a total of 10,496 different Operational Taxonomic Units (OTUs) were identified that were assigned to 40 bacterial and 9 archaeal phyla and 14 candidate divisions. On average, the most abundant phyla in all samples were Chloroflexi (Chloroflexota) (24.62%), followed by Proteobacteria (Pseudomonadota) (11.29%), Firmicutes (Bacillota) (10.73%), Crenarchaeota (Thermoproteota) (8.55%), and Acidobacteria (Acidobacteriota) (8.07%). At the genus level, a total of 286 known genera and candidate genera were mostly dominated by members of Bacillus, Thermoflexus, Desulfatiglans, Pseudoalteromonas, and Pseudomonas. Discussion: In most of the stations, the Chao1 values at the deeper layers were comparable to the surface sediment samples denoting the high diversity in the subsurface of these ecosystems. Heatmap analysis based on the 100 most abundant OTUs, grouped the sampling stations according to their geographical location, placing together the two hottest stations (up to 99°C). This result indicates that this specific area within the active Kolumbo crater create a distinct niche, where microorganisms with adaptation strategies to withstand heat stresses can thrive, such as the endospore-forming Firmicutes.

A Comparative BAC map for the gilthead sea bream (Sparus aurata L.).

This study presents the first comparative BAC map of the gilthead sea bream (Sparus aurata), a highly valuated marine aquaculture fish species in the Mediterranean. High-throughput end sequencing of a BAC library yielded 92,468 reads (60.6 Mbp). Comparative mapping was achieved by anchoring BAC end sequences to the three-spined stickleback (Gasterosteus aculeatus) genome. BACs that were consistently ordered along the stickleback chromosomes accounted for 14,265 clones. A fraction of 5,249 BACs constituted a minimal tiling path that covers 73.5% of the stickleback chromosomes and 70.2% of the genes that have been annotated. The N50 size of 1,485 "BACtigs" consisting of redundant BACs is 337,253 bp. The largest BACtig covers 2.15 Mbp in the stickleback genome. According to the insert size distribution of mapped BACs the sea bream genome is 1.71-fold larger than the stickleback genome. These results represent a valuable tool to researchers in the field and may support future projects to elucidate the whole sea bream genome.

Functional and transcriptomic investigation of laccase activity in the presence of PCB29 identifies two novel enzymes and the multicopper oxidase repertoire of a marine-derived fungus.

Polychlorinated biphenyls (PCBs) are persistent organic pollutants (POPs), that can be detected in a variety of environments including the human body, adversely affecting global health. Bioremediation is an emerging field for the detoxification and removal of environmental pollutants, with novel biocatalysts appropriate for this task being in high demand. In this study, a biobank of novel fungal strains isolated as symbionts of marine invertebrates was screened for their ability to remove 2,4,5-trichlorobiphenyl (PCB29). The most efficient strains were studied further for their ability to express laccase activity, the most commonly associated extracellular activity involved in the removal of aromatic pollutants and encoded in fungi by the enzymatic class of multicopper oxidases (MCOs). The strain expressing the highest laccase activity, Cladosporium sp. TM138-S3, was cultivated in the presence of copper ions in a 12 L bioreactor and two enzymes exhibiting laccase activity were isolated from the culture broth through ion-exchange chromatography. The two enzymes, Lac1 and Lac2, were biochemically characterized and showed similar characteristics, although an improved ability to remove PCB29 (up to 71.2%) was observed for Lac2 in the presence of mediators. In parallel, we performed RNAseq of the strain growing in presence and absence of PCB29 and reconstructed its transcriptome assembly. Functional annotation allowed identifying the MCO repertoire of the fungus, consisting of 13 enzymes. Phylogenetic analysis of Ascomycete MCOs further allowed classifying these enzymes, revealing the diversity of laccase activities in Cladosporium sp. TM138-S3.

Genomic adaptation of Pseudomonas strains to acidity and antibiotics in hydrothermal vents at Kolumbo submarine volcano, Greece.

Although the rise of antibiotic and multidrug resistant bacteria is one of the biggest current threats to human health, our understanding of the mechanisms involved in antibiotic resistance selection remains scarce. We performed whole genome sequencing of 21 Pseudomonas strains, previously isolated from an active submarine volcano of Greece, the Kolumbo volcano. Our goal was to identify the genetic basis of the enhanced co-tolerance to antibiotics and acidity of these Pseudomonas strains. Pangenome analysis identified 10,908 Gene Clusters (GCs). It revealed that the numbers of phage-related GCs and sigma factors, which both provide the mechanisms of adaptation to environmental stressors, were much higher in the high tolerant Pseudomonas strains compared to the rest ones. All identified GCs of these strains were associated with antimicrobial and multidrug resistance. The present study provides strong evidence that the CO2-rich seawater of the volcano associated with low pH might be a reservoir of microorganisms carrying multidrug efflux-mediated systems and pumps. We, therefore, suggest further studies of other extreme environments (or ecosystems) and their associated physicochemical parameters (or factors) in the rise of antibiotic resistance.

Micro-CT for Biological and Biomedical Studies: A Comparison of Imaging Techniques.

Several imaging techniques are used in biological and biomedical studies. Micro-computed tomography (micro-CT) is a non-destructive imaging technique that allows the rapid digitisation of internal and external structures of a sample in three dimensions and with great resolution. In this review, the strengths and weaknesses of some common imaging techniques applied in biological and biomedical fields, such as optical microscopy, confocal laser scanning microscopy, and scanning electron microscopy, are presented and compared with the micro-CT technique through five use cases. Finally, the ability of micro-CT to create non-destructively 3D anatomical and morphological data in sub-micron resolution and the necessity to develop complementary methods with other imaging techniques, in order to overcome limitations caused by each technique, is emphasised.

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.

Evolution of a multifunctional gene: The warm temperature acclimation protein Wap65 in the European seabass Dicentrarchus labrax.

The warm temperature acclimation protein Wap65 has been shown to be involved in temperature acclimation, in immune response as well as in development. In teleosts, two types of Wap65 proteins, Wap65-1 and Wap65-2 are found, both acting as a multifunctional agent in several biological processes. In the present study we identified both transcripts Wap65-1 and Wap65-2 for the European seabass (Dicentrarchus labrax), examined their evolutionary rate and performed selection tests. The two paralogues were shown to be under moderate positive selection indicating their evolutionary adaptation. This functional diversification was further explored through expression studies. Both transcripts were differentially expressed during development as well as in various tissues and pathogen challenges, showing that Wap65-1 and Wap65-2 have evolved diverse functions. These results direct to the hypothesis that Wap65 proteins may, similarly to heat-shock proteins, have a general role in cell physiology.

Profiling of infection specific mRNA transcripts of the European seabass Dicentrarchus labrax.

BACKGROUND: The European seabass (Dicentrarchus labrax), one of the most extensively cultured species in European aquaculture productions, is, along with the gilthead sea bream (Sparus aurata), a prospective model species for the Perciformes which includes several other commercially important species. Massive mortalities may be caused by bacterial or viral infections in intensive aquaculture production. Revealing transcripts involved in immune response and studying their relative expression enhances the understanding of the immune response mechanism and consequently also the creation of vaccines. The analysis of expressed sequence tags (EST) is an efficient and easy approach for gene discovery, comparative genomics and for examining gene expression in specific tissues in a qualitative and quantitative way. RESULTS: Here we describe the construction, analysis and comparison of a total of ten cDNA libraries, six from different tissues infected with V. anguillarum (liver, spleen, head kidney, gill, peritoneal exudates and intestine) and four cDNA libraries from different tissues infected with Nodavirus (liver, spleen, head kidney and brain). In total 9605 sequences representing 3075 (32%) unique sequences (set of sequences obtained after clustering) were obtained and analysed. Among the sequences several immune-related proteins were identified for the first time in the order of Perciformes as well as in Teleostei. CONCLUSION: The present study provides new information to the Gene Index of seabass. It gives a unigene set that will make a significant contribution to functional genomic studies and to studies of differential gene expression in relation to the immune system. In addition some of the potentially interesting genes identified by in silico analysis and confirmed by real-time PCR are putative biomarkers for bacterial and viral infections in fish.

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.

Microbial strains isolated from CO2-venting Kolumbo submarine volcano show enhanced co-tolerance to acidity and antibiotics.

As ocean acidification intensifies, there is growing global concern about the impacts that future pH levels are likely to have on marine life and ecosystems. By analogy, a steep decrease of seawater pH with depth is encountered inside the Kolumbo submarine volcano (northeast Santorini) as a result of natural CO2 venting, making this system ideal for ocean acidification research. Here, we investigated whether the increase of acidity towards deeper layers of Kolumbo crater had any effect on relevant phenotypic traits of bacterial isolates. A total of 31 Pseudomonas strains were isolated from both surface- (SSL) and deep-seawater layers (DSL), with the latter presenting a significantly higher acid tolerance. In particular, the DSL strains were able to cope with H+ levels that were 18 times higher. Similarly, the DSL isolates exhibited a significantly higher tolerance than SSL strains against six commonly used antibiotics and As(III). More importantly, a significant positive correlation was revealed between antibiotics and acid tolerance across the entire set of SSL and DSL isolates. Our findings imply that Pseudomonas species with higher resilience to antibiotics could be favored by the prospect of acidifying oceans. Further studies are required to determine if this feature is universal across marine bacteria and to assess potential ecological impacts.

Absence of spatial genetic structure in common dentex (Dentex dentex Linnaeus, 1758) in the Mediterranean Sea as evidenced by nuclear and mitochondrial molecular markers.

The common dentex, Dentex dentex, is a fish species which inhabits marine environments in the Mediterranean and Northeast Atlantic regions. This is an important species from an ecological, economic and conservation perspective, however critical information on its population genetic structure is lacking. Most samples were obtained from the Mediterranean Sea (17 sites) with an emphasis around Corsica (5 sites), plus one Atlantic Ocean site. This provided an opportunity to examine genetic structuring at local and broader scales to provide science based data for the management of fishing stocks in the region. Two mitochondrial regions were examined (D-loop and COI) along with eight microsatellite loci. The COI data was combined with publicly available sequences and demonstrated past misidentification of common dentex. All markers indicated the absence of population genetic structure from the Bay of Biscay to the eastern Mediterranean Sea. Bayesian approaches, as well as the statistical tests performed on the allelic frequencies from microsatellite loci, indicated low differentiation between samples; there was only a slight (p = 0.05) indication of isolation by distance. Common dentex is a marine fish species with a unique panmictic population in the Mediterranean and likely in the Atlantic Ocean as well.

A gene-based radiation hybrid map of the gilthead sea bream Sparus aurata refines and exploits conserved synteny with Tetraodon nigroviridis.

BACKGROUND: Comparative teleost studies are of great interest since they are important in aquaculture and in evolutionary issues. Comparing genomes of fully sequenced model fish species with those of farmed fish species through comparative mapping offers shortcuts for quantitative trait loci (QTL) detections and for studying genome evolution through the identification of regions of conserved synteny in teleosts. Here a comparative mapping study is presented by radiation hybrid (RH) mapping genes of the gilthead sea bream Sparus aurata, a non-model teleost fish of commercial and evolutionary interest, as it represents the worldwide distributed species-rich family of Sparidae. RESULTS: An additional 74 microsatellite markers and 428 gene-based markers appropriate for comparative mapping studies were mapped on the existing RH map of Sparus aurata. The anchoring of the RH map to the genetic linkage map resulted in 24 groups matching the karyotype of Sparus aurata. Homologous sequences to Tetraodon were identified for 301 of the gene-based markers positioned on the RH map of Sparus aurata. Comparison between Sparus aurata RH groups and Tetraodon chromosomes (karyotype of Tetraodon consists of 21 chromosomes) in this study reveals an unambiguous one-to-one relationship suggesting that three Tetraodon chromosomes correspond to six Sparus aurata radiation hybrid groups. The exploitation of this conserved synteny relationship is furthermore demonstrated by in silico mapping of gilthead sea bream expressed sequence tags (EST) that give a significant similarity hit to Tetraodon. CONCLUSION: The addition of primarily gene-based markers increased substantially the density of the existing RH map and facilitated comparative analysis. The anchoring of this gene-based radiation hybrid map to the genome maps of model species broadened the pool of candidate genes that mainly control growth, disease resistance, sex determination and reversal, reproduction as well as environmental tolerance in this species, all traits of great importance for QTL mapping and marker assisted selection. Furthermore this comparative mapping approach will facilitate to give insights into chromosome evolution and into the genetic make up of the gilthead sea bream.