Multi-omics for studying and understanding polar life.

Polar ecosystems are experiencing amongst the most rapid rates of regional warming on Earth. Here, we discuss 'omics' approaches to investigate polar biodiversity, including the current state of the art, future perspectives and recommendations. We propose a community road map to generate and more fully exploit multi-omics data from polar organisms. These data are needed for the comprehensive evaluation of polar biodiversity and to reveal how life evolved and adapted to permanently cold environments with extreme seasonality. We argue that concerted action is required to mitigate the impact of warming on polar ecosystems via conservation efforts, to sustainably manage these unique habitats and their ecosystem services, and for the sustainable bioprospecting of novel genes and compounds for societal gain.

A Weight of Evidence approach to support the assessment of the quality of Manila clam farming sites in a coastal lagoon.

Aquaculture productivity in coastal lagoons is endangered by a complex interplay of anthropogenic and environmental factors, amplified by the effects of climate change in these sensitive areas. To reach a more comprehensive assessment of farming sites quality, a quantitative Weight of Evidence approach (QWoE) is applied for the first time to data collected at four Manila clam (R. philippinarum) farming sites in the Venice lagoon (Italy). This included sediment quality, chemical bioaccumulation, and biological responses. Results revealed a greater hazard for sites closer to the open sea. In these areas, the combination of sediment characteristics and a higher frequency of salinity and temperature stress could explain the alterations measured at a transcriptional and biomarker level. The findings demonstrate that a QWoE approach that integrates multiple sources of evidence should also include physicochemical conditions in order to better understand the impacts of human activities and other stressors on clam aquaculture productivity.

Effects of environmental concentrations of the fragrance amyl salicylate on the mediterranean mussel Mytilus galloprovincialis.

Amyl salicylate (AS) is a fragrance massively used as a personal care product and following the discharged in wastewaters may end up in the aquatic environment representing a potential threat for the ecosystem and living organisms. AS was recently detected in water of the Venice Lagoon, a vulnerable area continuously subjected to the income of anthropogenic chemicals. The lagoon is a relevant area for mollusc farming, including the Mediterranean mussels (Mytilus galloprovincialis) having an important economic and ecological role. Despite high levels of AS occurred in water of the Lagoon of Venice, no studies investigated the possible consequences of AS exposures on species inhabiting this ecosystem to date. For the first time, we applied a multidisciplinary approach to investigate the potential effects of the fragrance AS on Mediterranean mussels. To reach such a goal, bioaccumulation, cellular, biochemical, and molecular analyses (RNA-seq and microbiota characterization) were measured in mussels treated for 7 and 14 days with different AS Venice lagoon environmental levels (0.1 and 0.5 µg L-1). Despite chemical investigations suggested low AS bioaccumulation capability, cellular and molecular analyses highlighted the disruption of several key cellular processes after the prolonged exposures to the high AS concentration. Among them, potential immunotoxicity and changes in transcriptional regulation of pathways involved in energy metabolism, stress response, apoptosis and cell death regulations have been observed. Conversely, exposure to the low AS concentration demonstrated weak transcriptional changes and transient increased representation of opportunistic pathogens, as Arcobacter genus and Vibrio aestuarianus. Summarizing, this study provides the first overview on the effects of AS on one of the most widely farmed mollusk species.

The effects of glyphosate and AMPA on the mediterranean mussel Mytilus galloprovincialis and its microbiota.

Glyphosate, the most widely used herbicide worldwide, targets the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) enzyme in the shikimate pathway found in plants and some microorganisms. While the potential for glyphosate to induce a broad range of biological effects in exposed organisms has been demonstrated, the global molecular mechanisms of toxicity and potential effects in bacterial symbionts remain unclear, in particular for ecologically important marine species such as bivalve molluscs. Here, the effects of glyphosate (GLY), its degradation product aminomethylphosphonic acid (AMPA), and a mixture of both (MIX) on the mussel M. galloprovincialis were assessed in a controlled experiment. For the first time, next generation sequencing (RNA-seq and 16S rRNA amplicon sequencing) was used to evaluate such effects at the molecular level in both the host and its respective microbiota. The results suggest that the variable capacity of bacterial species to proliferate in the presence of these compounds and the impairment of host physiological homeostasis due to AMPA and GLY toxicity may cause significant perturbations to the digestive gland microbiota, as well as elicit the spread of potential opportunistic pathogens such as Vibrio spp.. The consequent host-immune system activation identified at the molecular and cellular level could be aimed at controlling changes occurring in the composition of symbiotic microbial communities. Overall, our data raise further concerns about the potential adverse effects of glyphosate and AMPA in marine species, suggesting that both the effects of direct toxicity and the ensuing changes occurring in the host-microbial community must be taken into consideration to determine the overall ecotoxicological hazard of these compounds.

Ecotoxicological effects of the herbicide glyphosate in non-target aquatic species: Transcriptional responses in the mussel Mytilus galloprovincialis.

Glyphosate has been the most widely used herbicide worldwide over the last three decades, raising increasing concerns for its potential impacts on environmental and human health. Recent studies revealed that glyphosate occurs in soil, surface water, and groundwater, and residues are found at all levels of the food chain, such as drinking water, plants, animals, and even in humans. While research has demonstrated that glyphosate can induce a broad range of biological effects in exposed organisms, the global molecular mechanisms of action still need to be elucidated, in particular for marine species. In this study, we characterized for the first time the molecular mechanisms of action of glyphosate in a marine bivalve species after exposure to environmentally realistic concentrations. To reach such a goal, Mediterranean mussels Mytilus galloprovincialis, an ecologically and economically relevant species, were exposed for 21 days to 10, 100, and 1000 µg/L and digestive gland transcriptional profiles were investigated through RNA-seq. Differential expression analysis identified a total of 111, 124, and 211 differentially regulated transcripts at glyphosate concentrations of 10, 100, and 1000 µg/L, respectively. Five genes were found consistently differentially expressed at all investigated concentrations, including SERP2, which plays a role in the protection of unfolded target proteins against degradation, the antiapoptotic protein GIMAP5, and MTMR14, which is involved in macroautophagy. Functional analysis of differentially expressed genes reveals the disruption of several key biological processes, such as energy metabolism and Ca2+ homeostasis, cell signalling, and endoplasmic reticulum stress response. Together, the results obtained suggest that the presence of glyphosate in the marine ecosystem should raise particular concern because of its significant effects even at the lowest concentration.

Development of a real-time PCR assay for rapid detection and quantification of Photobacterium damselae subsp. piscicida in fish tissues.

The availability of a rapid and accurate method for the diagnosis of Photobacterium damselae subsp. piscicida (Phdp), able to discriminate its strictly correlated subsp. damselae (Phdd), formally known as Vibrio damsela, is essential for managing fish pasteurellosis outbreaks in farmed fish. A single-step, high-sensitivity real-time PCR assay for simultaneous detection and quantification of P. damselae was designed targeting partial of the sequence of the bamB gene and tested for specificity and sensitivity on laboratory-generated samples as well as on experimentally infected seabream tissue samples. With a limit of detection (LOD) of one copy in pure bacterial DNA, the sensitivity was higher than all methods previously reported. Validation in target and non-target bacterial species proved the assay was able to discriminate Phdd-Phdp subspecies from diverse hosts/geographical origins and between non-target species. In addition, two SNPs in the target amplicon region determine two distinctive qPCR dissociation curves distinguishing between Phdp-Phdd. This is the first time that a molecular method for P. damselae diagnosis combines detection, quantification and subspecies identification in one step. The assay holds the potential to improve the knowledge of infection dynamics and the development of better strategies to control an important fish disease.

Microbiota and environmental stress: how pollution affects microbial communities in Manila clams.

Given the crucial role of microbiota in host development, health, and environmental interactions, genomic analyses focusing on host-microbiota interactions should certainly be considered in the investigation of the adaptive mechanisms to environmental stress. Recently, several studies suggested that microbiota associated to digestive tract is a key, although still not fully understood, player that must be considered to assess the toxicity of environmental contaminants. Bacteria-dependent metabolism of xenobiotics may indeed modulate the host toxicity. Conversely, environmental variables (including pollution) may alter the microbial community and/or its metabolic activity leading to host physiological alterations that may contribute to their toxicity. Here, 16s rRNA gene amplicon sequencing has been applied to characterize the hepatopancreas microbiota composition of the Manila clam, Ruditapes philippinarum. The animals were collected in the Venice lagoon area, which is subject to different anthropogenic pressures, mainly represented by the industrial activities of Porto Marghera (PM). Seasonal and geographic differences in clam microbiotas were explored and linked to host response to chemical stress identified in a previous study at the transcriptome level, establishing potential interactions among hosts, microbes, and environmental parameters. The obtained results showed the recurrent presence of putatively detoxifying bacterial taxa in PM clams during winter and over-representation of several metabolic pathways involved in xenobiotic degradation, which suggested the potential for host-microbial synergistic detoxifying actions. Strong interaction between seasonal and chemically-induced responses was also observed, which partially obscured such potentially synergistic actions. Seasonal variables and exposure to toxicants were therefore shown to interact and substantially affect clam microbiota, which appeared to mirror host response to environmental variation. It is clear that understanding how animals respond to chemical stress cannot ignore a key component of such response, the microbiota.

Transcriptomic resources for environmental risk assessment: a case study in the Venice lagoon.

The development of new resources to evaluate the environmental status is becoming increasingly important representing a key challenge for ocean and coastal management. Recently, the employment of transcriptomics in aquatic toxicology has led to increasing initiatives proposing to integrate eco-toxicogenomics in the evaluation of marine ecosystem health. However, several technical issues need to be addressed before introducing genomics as a reliable tool in regulatory ecotoxicology. The Venice lagoon constitutes an excellent case, in which the assessment of environmental risks derived from the nearby industrial activities represents a crucial task. In this context, the potential role of genomics to assist environmental monitoring was investigated through the definition of reliable gene expression markers associated to chemical contamination in Manila clams, and their subsequent employment for the classification of Venice lagoon areas. Overall, the present study addresses key issues to evaluate the future outlooks of genomics in the environmental monitoring and risk assessment.

Phylogeography and demographic history of two widespread Indo-Pacific mudskippers (Gobiidae: Periophthalmus).

This study provides a first description of the phylogeographic patterns and evolutionary history of two species of the mudskipper genus Periophthalmus. These amphibious gobies are distributed throughout the whole Indo-Pacific region and Atlantic coast of Africa, in peritidal habitats of soft-bottom coastal ecosystems. Three sequence datasets of two widely distributed species, Periophthalmus argentilineatus and P. kalolo, were obtained by amplifying and sequencing two mtDNA markers (D-loop and 16S rDNA) and the nDNA rag1 region. The three datasets were then used to perform phylogeographic, demographic and population genetic analyses. Our results indicate that tectonic events and past climatic oscillations strongly contributed to shape present genetic differentiation, phylogeographic and demographic patterns. We found support for the monophyly of P. kalolo, and only shallow genetic differentiation between East-African and Indo-Malayan populations of this species. However, our collections of the morphospecies P. argentilineatus include three molecularly distinct lineages, one of them more closely related to P. kalolo. The presence of Miocenic timings for the most recent common ancestors of some of these morphologically similar clades, suggests the presence of strong stabilising selection in mudskippers' habitats. At population level, demographic analyses and palaeoecological records of mangrove ecosystems suggest that Pleistocene bottlenecks and expansion plus secondary contact events of the studied species were associated with recurrent sea transgressions during interglacials, and sea regressions or stable regimes during glacials, respectively.

Evidence for past and present hybridization in three Antarctic icefish species provides new perspectives on an evolutionary radiation.

Determining the timing, extent and underlying causes of interspecific gene exchange during or following speciation is central to understanding species' evolution. Antarctic notothenioid fish, thanks to the acquisition of antifreeze glycoproteins during Oligocene transition to polar conditions, experienced a spectacular radiation to >100 species during Late Miocene cooling events. The impact of recent glacial cycles on this group is poorly known, but alternating warming and cooling periods may have affected species' distributions, promoted ecological divergence into recurrently opening niches and/or possibly brought allopatric species into contact. Using microsatellite markers and statistical methods including Approximate Bayesian Computation, we investigated genetic differentiation, hybridization and the possible influence of the last glaciation/deglaciation events in three icefish species of the genus Chionodraco. Our results provide strong evidence of contemporary and past introgression by showing that: (i) a substantial fraction of contemporary individuals in each species has mixed ancestry, (ii) evolutionary scenarios excluding hybridization or including it only in ancient times have small or zero posterior probabilities, (iii) the data support a scenario of interspecific gene flow associated with the two most recent interglacial periods. Glacial cycles might therefore have had a profound impact on the genetic composition of Antarctic fauna, as newly available shelf areas during the warmer intervals might have favoured secondary contacts and hybridization between diversified groups. If our findings are confirmed in other notothenioids, they offer new perspectives for understanding evolutionary dynamics of Antarctic fish and suggest a need for new predictions on the effects of global warming in this group.

Detecting genome-wide gene transcription profiles associated with high pollution burden in the critically endangered European eel.

The European eel illustrates an example of a critically endangered fish species strongly affected by human stressors throughout its life cycle, in which pollution is considered to be one of the factors responsible for the decline of the stock. The objective of our study was to better understand the transcriptional response of European eels chronically exposed to pollutants in their natural environment. A total of 42 pre-migrating (silver) female eels from lowly, highly and extremely polluted environments in Belgium and, for comparative purposes, a lowly polluted habitat in Italy were measured for polychlorinated biphenyls (PCBs), organochlorine pesticides (OCPs) and brominated flame retardants (BFRs). Multipollutant level of bioaccumulation was linked to their genome-wide gene transcription using an eel-specific array of 14,913 annotated cDNAs. Shared responses to pollutant exposure were observed when comparing the highly polluted site in Belgium with the relatively clean sites in Belgium and Italy. First, an altered pattern of transcription of genes was associated with detoxification, with a novel European eel CYP3A gene and gluthatione S-transferase transcriptionally up-regulated. Second, an altered pattern of transcription of genes associated with the oxidative phosphorylation pathway, with the following genes involved in the generation of ATP being transcriptionally down-regulated in individuals from the highly polluted site: NADH dehydrogenase, succinate dehydrogenase, ubiquinol-cytochrome c reductase, cytochrome c oxidase and ATP synthase. Although we did not measure metabolism directly, seeing that the transcription level of many genes encoding enzymes involved in the mitochondrial respiratory chain and oxidative phosphorylation were down-regulated in the highly polluted site suggests that pollutants may have a significant effect on energy metabolism in these fish.

Quantitative trait loci for resistance to fish pasteurellosis in gilthead sea bream (Sparus aurata).

Fish pasteurellosis is a bacterial disease causing important losses in farmed fish, including gilthead sea bream, a teleost fish of great relevance in marine aquaculture. We report in this study a QTL analysis for resistance to fish pasteurellosis in this species. An experimental population of 500 offspring originating from eight sires and six dams in a single mass-spawning event was subjected to a disease challenge with Photobacterium damselae subsp. piscicida (Phdp), the causative agent of fish pasteurellosis. A total of 151 microsatellite loci were genotyped in the experimental population, and half-sib regression QTL analysis was carried out on two continuous traits, body length at time of death and survival, and for two binary traits, survival at day 7 and survival at day 15, when the highest peaks of mortality were observed. Two significant QTLs were detected for disease resistance. The first one was located on linkage group LG3 affecting late survival (survival at day 15). The second one, for overall survival, was located on LG21, which allowed us to highlight a potential marker (Id13) linked to disease resistance. A significant QTL was also found for body length at death on LG6 explaining 5-8% of the phenotypic variation.

Proposed new nomenclature for Bos taurus cytochromes P450 involved in xenobiotic drug metabolism.

The cytochrome P450 (CYP) superfamily of drug metabolizing enzymes (DMEs) plays a central role in the oxidative metabolism of xenobiotics to which living organisms are exposed. In Bos taurus (cattle), a definitive nomenclature for CYP proteins is still lacking, and to unambiguously settle cattle nomenclature a phylogenetic analysis of proteins belonging to CYP 1-4 families was performed. Sequences collected from GenBank and Dr Nelson's P450 homepage databases were analyzed according to the maximum likelihood method. Phylogenetic outputs showed that CYPs sharing the same name and collected from different species did not form, in several instances, monophyletic groups. Some cattle CYPs did not group with their supposed human orthologous counterparts, thus requiring a new nomenclature. Name changes mostly mirrored the orthologous counterparts established for other species, and new names were created when no clear orthologous sequences were identified. The new nomenclature will allow a more appropriate investigation of biochemical and molecular mechanisms involved in the expression and regulation of these DMEs.

Expression profiling of skeletal muscle in young bulls treated with steroidal growth promoters.

Dexamethasone (Dex), alone or in association with estrogens, is often illegally administered per os at very low dosage as a growth promoter in beef cattle, with effects that are opposite to the muscle wasting and atrophy induced by repeated administration at therapeutic dosages. In vitro and in vivo studies have investigated the catabolic effects of Dex at therapeutic doses on skeletal muscle, demonstrating an increase in the expression of GDF8 (myostatin) gene, a well-known negative regulator of skeletal muscle mass, in a dose-dependent way. This suggested a direct role of myostatin in Dex-induced muscle wasting. In the present study, an oligonucleotide microarray platform was used to compare expression profiles of beef cattle muscle in animals treated with either Dex or Dex plus 17-beta estradiol (Estr) administered at subtherapeutic dosage, against untreated controls. Data analysis demonstrates that the expression profiles were strongly affected by Dex treatment with hundreds of genes upregulated with relevant fold-change, whereas seven genes were downregulated including the myostatin gene. On the contrary, the number of differentially regulated genes was lower in response to the addition of Estr to the Dex treatment. Differentially regulated genes were analyzed to describe the effects of these treatments on muscle physiology, highlighting the importance of specific pathways (e.g., Wnt or cytokine signaling) and cellular processes (e.g., cell shape and motility). Finally, the observed differences in the expression profile will allow the development of indirect bio-markers to detect illegal Dex treatments in beef cattle using quantitative RT-PCR.

Demographic history and population structure of the Antarctic silverfish Pleuragramma antarcticum.

The Antarctic silverfish Pleuragramma antarcticum (Nototheniidae) is the most abundant pelagic fish inhabiting Antarctic waters. In this study we investigated, through partial sequencing of the D-loop mitochondrial region, samples collected at four different locations in the Southern Ocean, three in the Atlantic and one in the Pacific sector. Sampling was replicated in two different years at two locations. Sequence analysis showed a remarkably high polymorphism, with 110 haplotypes over the 256 investigated specimens, and about 80% of haplotypes occurring only once. Neutrality tests indicated that all samples were not at mutation-drift equilibrium, and suggested a past population expansion. This result was supported by the presence of a star-like topology in the D-loop gene tree, and by results of mismatch distribution. The start of the expansion was dated, using a specifically calibrated clock, between 111 and 126 thousand years ago. This value corresponds to the start of the cooling period that led to the last glaciation peak, and is in close agreement with a recently suggested range expansion for pelagic Antarctic ecosystems. Analysis of molecular variation indicated a small, though highly significant, value of differentiation between samples. This result, together with the lack of association between clades and geographical locations, indicates a weak population structure for the species.

Mitochondrial DNA reveals a mosaic pattern of phylogeographical structure in Atlantic and Mediterranean populations of anchovy (Engraulis encrasicolus).

This study extends the geographic coverage of a previous study of mitochondrial DNA restriction fragment length polymorphism in European anchovy. Both studies together include 24 samples representing 17 localities extending from the Black Sea, through the Mediterranean Sea to the eastern Atlantic as far south as Dakar, Senegal. Eighty-eight haplotypes define two clades (A and B) separated by 3.2% sequence divergence. Clade A has a star-like genealogy indicative of a recent population expansion. Clade B has a more complex genealogy, consisting of several haplotypes at intermediate frequencies. The distributions of these clades consist of a mosaic with abrupt changes between some areas and gradients between other areas. Clade A predominates the Black and Aegean seas, but is present throughout the Mediterranean. Unexpectedly, new data show that clade A is also at a high frequency in the Atlantic, from Portugal to at least Senegal. Overall, the level of genetic differentiation among populations is high (F(ST)=0.148, p<0.0001), with the greatest differences between basins. AMOVA reveals four main geographical groups: Atlantic, central Mediterranean, Aegean Sea, and Black Sea. Mismatch distribution clearly indicates historical bottleneck and population expansion for clade A, while for clade B such evidence is equivocal. This difference may reflect a range expansion for both clades, but with higher gene flow (Nm values) between demes for clade A. Both contemporary and historical processes are important in shaping the complex genetic population structure of European anchovy.

Intronic breakpoint definition and transcription analysis in DMD/BMD patients with deletion/duplication at the 5' mutation hot spot of the dystrophin gene.

Dystrophin mutations occurring at the 5' end of the gene frequently behave as exceptions to the "frame rule," their clinical severity being variable and often not related to the perturbation of the translation reading frame. The molecular mechanisms underlying the phenotypic variability of 5' dystrophin mutations have not been fully clarified. We have characterized the genomic breakpoints within introns 2, 6 and 7 and identified the splicing profiles in a cohort of DMD/BMD patients with deletion of dystrophin exons 3-7, 3-6 and duplication of exons 2-4. Our findings indicate that the occurrence of intronic cryptic promoter as well as corrective splicing events are unlikely to play a role in exons 3-7 deleted patients phenotypic variability. Our data suggest that re-initiation of translation could represent a major mechanism responsible for the production of a residual dystrophin in some patients with exons 3-7 deletion. Furthermore, we observed that the out-of-frame exon 2a is almost constantly spliced into a proportion of the dystrophin transcripts in the analysed patients. In the exons 2-4 duplicated DMD patient, producing both in-frame and out-of-frame transcripts, this splicing behaviour might represent a critical factor contributing to the severe phenotype. In conclusion, we suggest that multiple mechanisms may have a role in modulating the outcome of 5' dystrophin mutations, including recoding mechanisms and unusual splicing choices.

Expression analysis of androgen-responsive genes in the prostate of veal calves treated with anabolic hormones.

In order to identify indirect molecular biomarkers of anabolic treatments in veal calves, an animal experiment was performed using two combinations of growth promoters (consisting of boldenone undecylenate and estradiol benzoate, and of testosterone enantate and estradiol benzoate). We selected a set of 12 genes that are known to be androgen responsive in other mammalian species. The expression profile of this set of genes was analysed on prostate samples of veal calves using a real-time RT-PCR approach. For each selected gene the corresponding bovine sequence was obtained and a gene specific real-time assay was optimised and validated. The amplification was shown to be highly specific, linear and efficient. High reproducibility (<1%) and low-test variability (<2.5%) were also been achieved. Messenger RNA levels were quantified in prostate samples, non-parametric analysis of variance showed significant up-regulation of three genes (MAF, ESR1 and AR) and significant down-regulation of four genes (HMGCS1, HPGD, DBI, and LIM) in treated samples when compared with untreated controls. To assess the possibility of identifying hormone-treated animals by molecular means we performed a discriminant analysis that was effective in classifying treated and non-treated samples with an accuracy of 93%. Our results indicate that identification of treatment with steroid hormones in veal calves by means of gene expression analysis is a feasible approach and could be improved increasing both the number of genes and the number of controls analysed.

Direct identification of Photobacterium damselae subspecies piscicida by PCR-RFLP analysis.

Fish pasteurellosis is an infectious disease that affects several teleost species living in temperate marine waters. The pathogen responsible, Photobacterium damselae subspecies piscicida, shows high genetic similarity with P. damselae subsp. damselae, making subspecies discrimination extremely laborious. Here we report for the first time a PCR-RFLP method for the identification of P. damselae subsp. piscicida without prior isolation in pure culture. Genomic sequence information was obtained through cloning and sequencing of RAPD products. Two P. damselae-specific primer pairs were developed and tested on 17 strains of P. damselae subsp. piscicida, 10 strains of P. damselae subsp. damselae, and 6 closely related control species. High sensitivity was achieved in PCR amplification on serially diluted samples (<180 fg of pure bacterial DNA or <10 fg, depending on the amplified fragment). Restriction analysis of PCR products showed a unique digestion profile for all P. damselae subsp. piscicida strains. The same PCR-RFLP method was implemented on total DNA samples extracted from experimentally infected sea bream and sea bass. Positive results were obtained on fish with clear signs of the disease as well as on challenged, but asymptomatic, fish. The method presented here might provide a useful tool for both prevention and rapid diagnosis of fish pasteurellosis.