Establishment and characterization of cytochrome P450 1A1 CRISPR/Cas9 Knockout Bovine Foetal Hepatocyte Cell Line (BFH12)

The cytochrome P450 1A (CYP1A) subfamily of xenobiotic metabolizing enzymes (XMEs) consists of two different isoforms, namely CYP1A1 and CYP1A2, which are highly conserved among species. These two isoenzymes are involved in the biotransformation of many endogenous compounds as well as in the bioactivation of several xenobiotics into carcinogenic derivatives, thereby increasing the risk of tumour development. Cattle (Bos taurus) are one of the most important food-producing animal species, being a significant source of nutrition worldwide. Despite daily exposure to xenobiotics, data on the contribution of CYP1A to bovine hepatic metabolism are still scarce. The CRISPR/Cas9-mediated knockout (KO) is a useful method for generating in vivo and in vitro models for studying xenobiotic biotransformations. In this study, we applied the ribonucleoprotein (RNP)-complex approach to successfully obtain the KO of CYP1A1 in a bovine foetal hepatocyte cell line (BFH12). After clonal expansion and selection, CYP1A1 excision was confirmed at the DNA, mRNA and protein level. Therefore, RNA-seq analysis revealed significant transcriptomic changes associated with cell cycle regulation, proliferation, and detoxification processes as well as on iron, lipid and mitochondrial homeostasis. Altogether, this study successfully generates a new bovine CYP1A1 KO in vitro model, representing a valuable resource for xenobiotic metabolism studies in this important farm animal species.

Discovering the Protective Effects of Quercetin on Aflatoxin B1-Induced Toxicity in Bovine Foetal Hepatocyte-Derived Cells (BFH12)

Aflatoxin B1 (AFB1) induces lipid peroxidation and mortality in bovine foetal hepatocyte-derived cells (BFH12), with underlying transcriptional perturbations associated mainly with cancer, cellular damage, inflammation, bioactivation, and detoxification pathways. In this cell line, curcumin and resveratrol have proven to be effective in mitigating AFB1-induced toxicity. In this paper, we preliminarily assessed the potential anti-AFB1 activity of a natural polyphenol, quercetin (QUE), in BFH12 cells. To this end, we primarily measured QUE cytotoxicity using a WST-1 reagent. Then, we pre-treated the cells with QUE and exposed them to AFB1. The protective role of QUE was evaluated by measuring cytotoxicity, transcriptional changes (RNA-sequencing), lipid peroxidation (malondialdehyde production), and targeted post-transcriptional modifications (NQO1 and CYP3A enzymatic activity). The results demonstrated that QUE, like curcumin and resveratrol, reduced AFB1-induced cytotoxicity and lipid peroxidation and caused larger transcriptional variations than AFB1 alone. Most of the differentially expressed genes were involved in lipid homeostasis, inflammatory and immune processes, and carcinogenesis. As for enzymatic activities, QUE significantly reverted CYP3A variations induced by AFB1, but not those of NQO1. This study provides new knowledge about key molecular mechanisms involved in QUE-mediated protection against AFB1 toxicity and encourages in vivo studies to assess QUE's bioavailability and beneficial effects on aflatoxicosis.

An In Vivo Whole-Transcriptomic Approach to Assess Developmental and Reproductive Impairments Caused by Flumequine in Daphnia magna

Among veterinary antibiotics, flumequine (FLU) is still widely used in aquaculture due to its efficacy and cost-effectiveness. Although it was synthesized more than 50 years ago, a complete toxicological framework of possible side effects on non-target species is still far from being achieved. The aim of this research was to investigate the FLU molecular mechanisms in Daphnia magna, a planktonic crustacean recognized as a model species for ecotoxicological studies. Two different FLU concentrations (2.0 mg L^-1 and 0.2 mg L^-1) were assayed in general accordance with OECD Guideline 211, with some proper adaptations. Exposure to FLU (2.0 mg L^-1) caused alteration of phenotypic traits, with a significant reduction in survival rate, body growth, and reproduction. The lower concentration (0.2 mg L^-1) did not affect phenotypic traits but modulated gene expression, an effect which was even more evident under the higher exposure level. Indeed, in daphnids exposed to 2.0 mg L^-1 FLU, several genes related with growth, development, structural components, and antioxidant response were significantly modulated. To the best of our knowledge, this is the first work showing the impact of FLU on the transcriptome of D. magna.

Effect of the dietary supplementation with extracts of chestnut wood and grape pomace on performance and jejunum response in female and male broiler chickens at different ages

BACKGROUND

Recently, interest in the use of herbs and phytogenic compounds has grown because of their potential role in the production and health of livestock animals. Among these compounds, several tannins have been tested in poultry, but those from chestnut wood and grape-industry byproducts have attracted remarkable interest. Thus, the present study aimed to gain further insights into the mechanisms involved in the response to the dietary supplementation with extracts of chestnut wood or grape pomace. To this purpose, 864 broiler chickens were fed a control diet (C) or the same diet supplemented 0.2% chestnut wood (CN) extract or 0.2% grape pomace (GP) extract from hatching until commercial slaughtering (at 45 days of age) to assess their effects on performance, meat quality, jejunum immune response and whole-transcriptome profiling in both sexes at different ages (15 and 35 d).

RESULTS

Final live weight and daily weight gain significantly increased (P < 0.01) in chickens fed GP diets compared to CN and C diets. The villi height was lower in chickens fed the CN diet than in those fed the C diet (P < 0.001); moreover, a lower density of CD45⁺ cells was observed in chickens fed the CN diet (P < 0.05) compared to those fed the C and GP diets. Genes involved in either pro- or anti-inflammatory response pathways, and antimicrobial and antioxidant responses were affected by GP and CN diets. There was no effect of the dietary treatment on meat quality. Regarding sex, in addition to a lower growth performance, females showed a lower occurrence of wooden breast (16.7% vs. 55.6%; P < 0.001) and a higher occurrence of spaghetti meat (48.6% vs. 4.17%; P < 0.001) in pectoralis major muscles after slaughtering than those in males. Based on the results of whole-transcriptome profiling, a significant activation of some molecular pathways related to immunity was observed in males compared with those of females.

CONCLUSIONS

The GP supplementation improved chicken performance and promoted immune responses in the intestinal mucosa; moreover, age and sex were associated with the most relevant transcriptional changes.

Deepening the Whole Transcriptomics of Bovine Liver Cells Exposed to AFB1: A Spotlight on Toll-like Receptor 2

Aflatoxin B1 (AFB1) is a food contaminant metabolized mostly in the liver and leading to hepatic damage. Livestock species are differently susceptible to AFB1, but the underlying mechanisms of toxicity have not yet been fully investigated, especially in ruminants. Thus, the aim of the present study was to better characterize the molecular mechanism by which AFB1 exerts hepatotoxicity in cattle. The bovine fetal hepatocyte cell line (BFH12) was exposed for 48 h to three different AFB1 concentrations (0.9 µM, 1.8 µM and 3.6 µM). Whole-transcriptomic changes were measured by RNA-seq analysis, showing significant differences in the expression of genes mainly involved in inflammatory response, oxidative stress, drug metabolism, apoptosis and cancer. As a confirmatory step, post-translational investigations on genes of interest were implemented. Cell death associated with necrosis rather than apoptosis events was noted. As far as the toxicity mechanism is concerned, a molecular pathway linking inflammatory response and oxidative stress was postulated. Toll-Like Receptor 2 (TLR2) activation, consequent to AFB1 exposure, triggers an intracellular signaling cascade involving a kinase (p38β MAPK), which in turn allows the nuclear translocation of the activator protein-1 (AP-1) and NF-κB, finally leading to the release of pro-inflammatory cytokines. Furthermore, a p38β MAPK negative role in cytoprotective genes regulation was postulated. Overall, our investigations improved the actual knowledge on the molecular effects of this worldwide relevant natural toxin in cattle.

Curcumin Mitigates AFB1-Induced Hepatic Toxicity by Triggering Cattle Antioxidant and Anti-inflammatory Pathways: A Whole Transcriptomic In Vitro Study

Aflatoxin B1 (AFB1) toxicity in livestock and human beings is a major economic and health concern. Natural polyphenolic substances with antioxidant properties have proven to be effective in ameliorating AFB1-induced toxicity. Here we assessed the potential anti-AFB1 activity of curcumin (pure curcumin, C, and curcumin from Curcuma longa, CL) in a bovine fetal hepatocyte-derived cell line (BFH12). First, we measured viability of cells exposed to AFB1 in presence or absence of curcumin treatment. Then, we explored all the transcriptional changes occurring in AFB1-exposed cells cotreated with curcumin. Results demonstrated that curcumin is effective in reducing AFB1-induced toxicity, decreasing cells mortality by approximately 30%. C and CL induced similar transcriptional changes in BFH12 exposed to AFB1, yet C treatment resulted in a larger number of significant genes compared to CL. The mitigating effects of curcuminoids towards AFB1 toxicity were mainly related to molecular pathways associated with antioxidant and anti-inflammatory response, cancer, and drug metabolism. Investigating mRNA changes induced by curcumin in cattle BFH12 cells exposed to AFB1 will help us to better characterize possible tools to reduce its consequences in this susceptible and economically important food-producing species.

Midazolam oxidation in cattle liver microsomes: The role of cytochrome P450 3A

In humans, the cytochrome P450 3A (CYP3A) subfamily is involved in midazolam (MDZ) biotransformation into 1'- and 4-hydroxy metabolites, and the former serves as a probe for CYP3A catalytic activity. In veterinary species is still crucial to identify enzyme- and species-specific CYP substrates; thus, the aim of this study was to characterize MDZ oxidation in cattle liver. A HPLC-UV method was used to measure 1'- and 4-hydroxy MDZ (1'- and 4-OHMDZ, respectively) formation in cattle liver microsomes and assess the role of CYP3A by an immunoinhibition study. Moreover, MDZ hydroxylation was evaluated in 300 cattle liver samples and results were correlated with testosterone hydroxylation. Formation of both metabolites conformed to a single-enzyme Michaelis-Menten kinetics. Values of V_max and K_m were 0.67 nmol/min/mg protein and 6.16 μM for 4-OHMDZ, and 0.06 nmol/min/mg protein and 10.08 μM for 1'-OHMDZ. An anti-rat CYP3A1 polyclonal antibody inhibited up to 50% and 94% 1'- and 4-OHMDZ formation, respectively. MDZ oxidation in liver microsomes was poorly correlated with testosterone hydroxylation. In conclusion, cattle metabolized MDZ to 1'-OHMDZ and 4-OHMDZ. The immunoinhibition results indicated a major contribution of CYP3As to 4-OHMDZ formation and the involvement of other CYPs in 1'-OHMDZ production, paving the way for further investigations.

Effect of dietary supplementation with yeast cell wall extracts on performance and gut response in broiler chickens

Background

The dietary supplementation of yeast cell wall extracts (YCW) has been found to reduce pathogenic bacteria load, promote immunoglobulin production, prevent diseases by pro-inflammatory responses, and alter gut microbiota composition. This study evaluated growth and slaughter results, health, gut morphology, immune status and gut transcriptome of 576 male chickens fed two diets, i.e. C (control) or Y (with 250-500 g/t of YCW fractions according to the growth period). At 21 and 42 d the jejunum of 12 chickens per diet were sampled and stained with hematoxylin/eosin for morphometric evaluation, with Alcian-PAS for goblet cells, and antibodies against CD3+ intraepithelial T-cells and CD45+ intraepithelial leukocytes. The jejunum sampled at 42 d were also used for whole-transcriptome profiling.

Results

Dietary YCW supplementation did not affect final live weight, whereas it decreased feed intake (114 to 111 g/d; P ≤ 0.10) and improved feed conversion (1.74 to 1.70; P ≤ 0.01). Regarding the gut, YCW supplementation tended to increase villi height (P = 0.07); it also increased the number of goblet cells and reduced the density of CD45+ cells compared to diet C (P < 0.001). In the gut transcriptome, four genes were expressed more in broilers fed diet Y compared to diet C, i.e. cytochrome P450, family 2, subfamily C, polypeptide 23b (CYP2C23B), tetratricopeptide repeat domain 9 (TTC9), basic helix-loop-helix family member e41 (BHLHE41), and the metalloreductase STEAP4. Only one gene set (HES_PATHWAY) was significantly enriched among the transcripts more expressed in broilers fed diet Y. However, a total of 41 gene sets were significantly over-represented among genes up-regulated in control broilers. Notably, several enriched gene sets are implicated in immune functions and related to NF-κB signaling, apoptosis, and interferon signals.

Conclusions

The dietary YCW supplementation improved broiler growth performance, increased gut glycoconjugate secretion and reduced the inflammatory status together with differences in the gut transcriptome, which can be considered useful to improve animal welfare and health under the challenging conditions of intensive rearing systems in broiler chickens.

Molecular insights into post-mating immune response in a fish with internal fertilization

The tight connection between immunity and reproduction has been studied for decades. However, basic knowledge at the molecular level of the effect of mating on immune function is still lacking in many taxa. Determining whether and how the immune system is engaged after mating is a crucial step in understanding post-mating mechanisms of reproduction and sexual selection. Here, we study the transcriptional changes in immunity-related genes caused by the ejaculate in the female reproductive tract using a model species for sexual selection studies, the guppy Poecilia reticulata. To study changes triggered by the ejaculate only, rather than caused by mating, we used artificial inseminations to transfer ejaculate into females. We then compared gene expression in the reproductive tract (gonoduct and ovary) of females artificially inseminated either with ejaculate or with a control solution, after 1 hr and after 6 hr. Overall, contact with ejaculate caused short-term changes in the expression of immune-related genes in the female reproductive tract, with a complex pattern of up- and down-regulation of immune-related pathways, but with clear indication of a marked down-regulation of the immune system shortly after ejaculate contact. This suggests a link between immune function and processes occurring between mating and fertilization in this species.

High temperature induces transcriptomic changes in Crassostrea gigas that hinders progress of Ostreid herpesvirus (OsHV-1) and promotes survival

Among all the environmental factors, seawater temperature plays a decisive role in triggering marine diseases. Like fever in vertebrates, high seawater temperature could modulate the host response to the pathogens in ectothermic animals. In France, massive mortality of Pacific oysters Crassostrea gigas caused by the ostreid herpesvirus 1 (OsHV-1) is markedly reduced when temperatures exceed 24°C in the field. In the present study we assess how high temperature influences the host response to the pathogen by comparing transcriptomes (RNA-sequencing) during the course of experimental infection at 21°C (reference) and 29°C. We show that high temperature induced host physiological processes that are unfavorable to the viral infection. Temperature influenced the expression of transcripts related to the immune process and increased the transcription of genes related to apoptotic process, synaptic signaling, and protein processes at 29°C. Concomitantly, the expression of genes associated to catabolism, metabolites transport, macromolecules synthesis and cell growth remained low since the first stage of infection at 29°C. Moreover, viral entry into the host might have been limited at 29°C by changes in extracellular matrix composition and protein abundance. Overall, these results provide new insights into how environmental factors modulate the host-pathogen interactions.

Draft genome assembly and transcriptome data of the icefish Chionodraco myersi reveal the key role of mitochondria for a life without hemoglobin at subzero temperatures

Antarctic fish belonging to Notothenioidei represent an extraordinary example of radiation in the cold. In addition to the absence of hemoglobin, icefish show a number of other striking peculiarities including large-diameter blood vessels, high vascular densities, mitochondria-rich muscle cells, and unusual mitochondrial architecture. In order to investigate the bases of icefish adaptation to the extreme Southern Ocean conditions we sequenced the complete genome of the icefish Chionodraco myersi. Comparative analyses of the icefish genome with those of other teleost species, including two additional white-blooded and five red-blooded notothenioids, provided a new perspective on the evolutionary loss of globin genes. Muscle transcriptome comparative analyses against red-blooded notothenioids as well as temperate fish revealed the peculiar regulation of genes involved in mitochondrial function in icefish. Gene duplication and promoter sequence divergence were identified as genome-wide patterns that likely contributed to the broad transcriptional program underlying the unique features of icefish mitochondria.

Significance of the goby Zosterisessor ophiocephalus as a sentinel species for Venice Lagoon contamination: Combining biomarker responses and bioaccumulation

The Venice Lagoon is an interesting example of an ecosystem suffering for a considerable anthropogenic impact, resulting in high concentrations of persistent organic pollutants (POPs) in lagoon sediments and seafood. In this context, biomonitoring is a crucially important task. The present study aimed at evaluating the validity of a multiple biomarker approach in a benthic fish species. A total of 567 Zosterisessor ophiocephalus (Gobiidae) fish were collected in spring and autumn from three areas of Venice Lagoon (Porto Marghera, Val di Brenta, and Cà Roman) showing high, intermediate and low amounts of POPs, respectively. Aryl hydrocarbon receptor (AHR) and cytochrome P450 1A (CYP1A) mRNA levels, CYP1A protein amount and ethoxyresorufin O-deethylase activity (EROD) were measured in pooled liver and gills (mRNA levels only). Such biological data were then compared with polychlorinated biphenyls (PCBs) residues, measured in grass goby muscle by gas chromatography. Aryl hydrocarbon receptor and CYP1A mRNAs, protein and EROD were upregulated in accordance with PCB amounts measured in Z. ophiocephalus muscles. In fact, the highest AHR and CYP1A induction was observed in fish sampled in close proximity of the industrial area of Porto Marghera. Overall, the present study confirm the grass goby as a reliable sentinel species for Venice Lagoon, and AHR/CYP1A/EROD as a sensitive set of biomarkers of exposure for AHR ligands.

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.

Revealing Mytilus galloprovincialis transcriptomic profiles during ontogeny

Mediterranean mussels are a worldwide spread bivalve species with extraordinary biological success. One of the reasons of this success could be the reproduction strategy of bivalves, characterized by the presence of trochophore larvae. Larval development in bivalves has been a topic of raising interest in the scientific community but it deserves much more attention. The principal objective of this work was to study the transcriptomic profile of the ontogeny of Mytilus galloprovincialis analyzing the gene expression in different developmental stages, from oocytes to juveniles. For this purpose, after conducting a 454 sequencing of the transcriptomes of mussel hemocytes, adult tissues and larvae, a new DNA microarray was designed and developed. The studied developmental stages: unfertilized oocytes, veliger, pediveliger, settled larvae and juveniles, showed very different transcriptomic profiles and clustered in groups defining their characteristic gene expression along ontogeny. Our results show that oocytes present a distinct and characteristic transcriptome. After metamorphosis, both settled larvae and juveniles showed a very similar transcriptome, with no enriched GO terms found between these two stages. This suggests: 1.- the progressive loss of RNA of maternal origin through larval development and 2.- the stabilization of the gene expression after settlement. On the other hand during metamorphosis a specific profile of differentially expressed genes was found. These genes were related to processes such as differentiation and biosynthesis. Processes related to the immune response were strongly down regulated. These suggest a development commitment at the expense of other non-essential functions, which are temporary set aside. Immune genes such as antimicrobial peptides suffer a decreased expression during metamorphosis. In fact, we found that the oocytes which express a higher quantity of genes such as myticins are more likely to reach success of the offspring, compared to oocytes poor in such mRNAs, whose progeny died before reaching metamorphosis.

The Voltage-Dependent Anion Channel (VDAC) of Pacific Oysters Crassostrea gigas Is Upaccumulated During Infection by the Ostreid Herpesvirus-1 (OsHV-1): an Indicator of the Warburg Effect

Voltage-dependent anion channel (VDAC) is a key mitochondrial protein. VDAC drives cellular energy metabolism by controlling the influx and efflux of metabolites and ions through the mitochondrial membrane, playing a role in its permeabilization. This protein exerts a pivotal role during the white spot syndrome virus (WSSV) infection in shrimp, through its involvement in a particular metabolism that plays in favor of the virus, the Warburg effect. The Warburg effect corresponds to an atypical metabolic shift toward an aerobic glycolysis that provides energy for rapid cell division and resistance to apoptosis. In the Pacific oyster Crassostrea gigas, the Warburg effect occurs during infection by Ostreid herpesvirus (OsHV-1). At present, the role of VDAC in the Warburg effect, OsHV-1 infection and apoptosis is unknown. Here, we developed a specific antibody directed against C. gigas VDAC. This tool allowed us to quantify the tissue-specific expression of VDAC, to detect VDAC oligomers, and to follow the amount of VDAC in oysters deployed in the field. We showed that oysters sensitive to a mortality event in the field presented an accumulation of VDAC. Finally, we propose to use VDAC quantification as a tool to measure the oyster susceptibility to OsHV-1 depending on its environment.

Long-lasting antiviral innate immune priming in the Lophotrochozoan Pacific oyster, Crassostrea gigas

In the last decade, a paradigm shift has emerged in comparative immunology. Invertebrates can no longer be considered to be devoid of specific recognition and immune memory. However, we still lack a comprehensive view of these phenomena and their molecular mechanisms across phyla, especially in terms of duration, specificity, and efficiency in a natural context. In this study, we focused on a Lophotrochozoan/virus interaction, as antiviral priming is mostly overlooked in molluscs. Juvenile Crassostrea gigas oysters experience reoccurring mass mortalities events from Ostreid herpes virus 1 with no existing therapeutic treatment. Our results showed that various nucleic acid injections can prime oysters to trigger an antiviral state ultimately protecting them against a subsequent viral infection. Focusing on poly(I:C) as elicitor, we evidenced that it protected from an environmental infection, by mitigating viral replication. That protection seemed to induce a specific antiviral response as poly(I:C) fails to protect against a pathogenic bacteria. Finally, we showed that this phenomenon was long-lasting, persisting for at least 5 months thus suggesting for the first time the existence of innate immune memory in this invertebrate species. This study strengthens the emerging hypotheses about the broad conservation of innate immune priming and memory mechanisms in Lophotrochozoans.

The highly rearranged mitochondrial genomes of the crabs Maja crispata and Maja squinado (Majidae) and gene order evolution in Brachyura

We sequenced the mitochondrial genomes of the spider crabs Maja crispata and Maja squinado (Majidae, Brachyura). Both genomes contain the whole set of 37 genes characteristic of Bilaterian genomes, encoded on both α- and β-strands. Both species exhibit the same gene order, which is unique among known animal genomes. In particular, all the genes located on the β-strand form a single block. This gene order was analysed together with the other nine gene orders known for the Brachyura. Our study confirms that the most widespread gene order (BraGO) represents the plesiomorphic condition for Brachyura and was established at the onset of this clade. All other gene orders are the result of transformational pathways originating from BraGO. The different gene orders exhibit variable levels of genes rearrangements, which involve only tRNAs or all types of genes. Local homoplastic arrangements were identified, while complete gene orders remain unique and represent signatures that can have a diagnostic value. Brachyura appear to be a hot-spot of gene order diversity within the phylum Arthropoda. Our analysis, allowed to track, for the first time, the fully evolutionary pathways producing the Brachyuran gene orders. This goal was achieved by coupling sophisticated bioinformatic tools with phylogenetic analysis.

Transcriptomic features of Pecten maximus oocyte quality and maturation

The king scallop Pecten maximus is a high valuable species of great interest in Europe for both fishery and aquaculture. Notably, there has been an increased investment to produce seed for enhancement programmes of wild scallop populations. However, hatchery production is a relatively new industry and it is still underdeveloped. Major hurdles are spawning control and gamete quality. In the present study, a total of 14 scallops were sampled in the bay of Brest (Brittany, France) to compare transcriptomic profiles of mature oocytes collected by spawning induction or by stripping. To reach such a goal, a microarray analysis was performed by using a custom 8x60K oligonucleotide microarray representing 45,488 unique scallop contigs. First we identified genes that were differentially expressed depending on oocyte quality, estimated as the potential to produce D-larvae. Secondly, we investigated the transcriptional features of both stripped and spawned oocytes. Genes coding for proteins involved in cytoskeletal dynamics, serine/threonine kinases signalling pathway, mRNA processing, response to DNA damage, apoptosis and cell-cycle appeared to be of crucial importance for both oocyte maturation and developmental competence. This study allowed us to dramatically increase the knowledge about transcriptional features of oocyte quality and maturation, as well as to propose for the first time putative molecular markers to solve a major bottleneck in scallop aquaculture.

Long dsRNAs promote an anti-viral response in Pacific oyster hampering ostreid herpesvirus 1 replication

Double stranded RNA-mediated genetic interference (RNAi) is a widely used reverse genetic tool for determining the loss-of-function phenotype of a gene. Here, the possible induction of an immune response by long dsRNA was tested in a marine bivalve, i.e. Crassostrea gigas, as well as the specific role of the subunit 2 of the nuclear factor κB inhibitor (IκB2). This gene is a candidate of particular interest for functional investigations in the context of massive mortality oyster events as Cg-IκB2 mRNA levels exhibited significant variation depending on the amount of ostreid herpesvirus 1 (OsHV-1) DNA detected. In the present study, dsRNAs targeting Cg-IκB2 and Green Fluorescence Protein genes were injected in vivo into oysters before being challenged by OsHV-1. Survival appeared close to 100% in both dsRNA injected conditions associated with a low detection of viral DNA and a low expression of a panel of 39 OsHV-1 genes as compared to infected control. Long dsRNA molecules, both Cg-IκB2- and GFP-dsRNA, may have induced an anti-viral state controlling the OsHV-1 replication and precluding the understanding of the Cg-IκB2 specific role. Immune-related genes including Cg-IκB1, Cg-Rel1, Cg-IFI44, Cg-PKR, and Cg-IAP appeared activated in dsRNA-injected condition potentially hampering viral replication and thus conferring a better resistance to OsHV-1 infection. We revealed that long dsRNA-mediated genetic interference triggered an anti-viral state in the oyster, emphasizing the need of new reverse genetics tools for assessing immune gene function and avoiding off-target effects in bivalves.

An integrated genomic approach for the study of mandibular prognathism in the European seabass (Dicentrarchus labrax)

Skeletal anomalies in farmed fish are a relevant issue affecting animal welfare and health and causing significant economic losses. Here, a high-density genetic map of European seabass for QTL mapping of jaw deformity was constructed and a genome-wide association study (GWAS) was carried out on a total of 298 juveniles, 148 of which belonged to four full-sib families. Out of 298 fish, 107 were affected by mandibular prognathism (MP). Three significant QTLs and two candidate SNPs associated with MP were identified. The two GWAS candidate markers were located on ChrX and Chr17, both in close proximity with the peaks of the two most significant QTLs. Notably, the SNP marker on Chr17 was positioned within the Sobp gene coding region, which plays a pivotal role in craniofacial development. The analysis of differentially expressed genes in jaw-deformed animals highlighted the “nervous system development” as a crucial pathway in MP. In particular, Zic2, a key gene for craniofacial morphogenesis in model species, was significantly down-regulated in MP-affected animals. Gene expression data revealed also a significant down-regulation of Sobp in deformed larvae. Our analyses, integrating transcriptomic and GWA methods, provide evidence for putative mechanisms underlying seabass jaw deformity.

Transcriptomic Changes in Liver of Young Bulls Caused by Diets Low in Mineral and Protein Contents and Supplemented with n-3 Fatty Acids and Conjugated Linoleic Acid

The aim of the present study was to identify transcriptional modifications and regulatory networks accounting for physiological and metabolic responses to specific nutrients in the liver of young Belgian Blue × Holstein bulls using RNA-sequencing. A larger trial has been carried out in which animals were fed with different diets: 1] a conventional diet; 2] a low-protein/low-mineral diet (low-impact diet) and 3] a diet enriched in n-3 fatty acids (FAs), conjugated linoleic acid (CLA) and vitamin E (nutraceutical diet). The initial hypothesis was that the administration of low-impact and nutraceutical diets might influence the transcriptional profiles in bovine liver and the resultant nutrient fluxes, which are essential for optimal liver function and nutrient interconversion. Results showed that the nutraceutical diet significantly reduced subcutaneous fat covering in vivo and liver pH. Dietary treatments did not affect overall liver fat content, but significantly modified the liver profile of 33 FA traits (out of the total 89 identified by gas-chromatography). In bulls fed nutraceutical diet, the percentage of n-3 and CLA FAs increased around 2.5-fold compared with the other diets, whereas the ratio of n6/n3 decreased 2.5-fold. Liver transcriptomic analyses revealed a total of 198 differentially expressed genes (DEGs) when comparing low-impact, nutraceutical and conventional diets, with the nutraceutical diet showing the greatest effects on liver transcriptome. Functional analyses using ClueGo and Ingenuity Pathway Analysis evidenced that DEGs in bovine liver were variously involved in energy reserve metabolic process, glutathione metabolism, and carbohydrate and lipid metabolism. Modifications in feeding strategies affected key transcription factors regulating the expression of several genes involved in fatty acid metabolism, e.g. insulin-induced gene 1, insulin receptor substrate 2, and RAR-related orphan receptor C. This study provides noteworthy insights into the molecular changes occurring as a result of nutrient variation in diets (aimed at reducing the environmental impact and improving human health) and broadens our understanding of the relationship between nutrients variation and phenotypic effects.

Can ecological history influence response to pollutants? Transcriptomic analysis of Manila clam collected in different Venice lagoon areas and exposed to heavy metal

Chronic exposure to environmental pollutants can exert strong selective pressures on natural populations, favoring the transmission over generations of traits that enable individuals to survive and thrive in highly impacted environments. The lagoon of Venice is an ecosystem subject to heavy anthropogenic impact, mainly due to the industrial activities of Porto Marghera (PM), which led to a severe chemical contamination of soil, groundwater, and sediments. Gene expression analysis on wild Manila clams collected in different Venice lagoon areas enabled to identify differences in gene expression profiles between clams collected in PM and those sampled in clean areas, and the definition of molecular signatures of chemical stress. However, it remains largely unexplored to which extent modifications of gene expression patterns persists after removing the source of contamination. It is also relatively unknown whether chronic exposure to xenobiotics affects the response to other chemical pollutants. To start exploring such issues, in the present study a common-garden experiment was coupled with transcriptomic analysis, to compare gene expression profiles of PM clams with those of clams collected in the less impacted area of Chioggia (CH) during a period under the same control conditions. Part of the two experimental groups were also exposed to copper for seven days to assess whether different "ecological history" does influence response to such pollutant. The results obtained suggest that the chronic exposure to chemical pollution generated a response at the transcriptional level that persists after removal for the contaminated site. These transcriptional changes are centered on key biological processes, such as defense against either oxidative stress or tissue/protein damage, and detoxification, suggesting an adaptive strategy for surviving in the deeply impacted environment of Porto Marghera. On the other hand, CH clams appeared to respond more effectively to copper exposure than PM animals, proposing that chronic exposure to chemical toxicants either lowers the sensitivity to additional toxicants or blunts the capacity to respond to novel chemical challenges in PM clams.

Extending RAD tag analysis to microbial ecology: a comparison between MultiLocus Sequence Typing and 2b-RAD to investigate Listeria monocytogenes genetic structure

The advent of next-generation sequencing (NGS) has dramatically changed bacterial typing technologies, increasing our ability to differentiate bacterial isolates. Despite it is now possible to sequence a bacterial genome in a few days and at reasonable costs, most genetic analyses do not require whole-genome sequencing, which also remains impractical for large population samples due to the cost of individual library preparation and bioinformatics. More traditional sequencing approaches, however, such as MultiLocus Sequence Typing (mlst) are quite laborious and time-consuming, especially for large-scale analyses. In this study, a genotyping approach based on restriction site-associated (RAD) tag sequencing, 2b-RAD, was applied to characterize Listeria monocytogenes strains. To verify the feasibility of the method, an in silico analysis was performed on 30 available complete genomes. For the same set of strains, in silico mlst analysis was conducted as well. Subsequently, 2b-RAD and mlst analyses were experimentally carried out on 58 isolates collected from food samples or food-processing sites. The obtained results demonstrate that 2b-RAD predicts mlst types and often provides more detailed information on population structure than mlst. Moreover, the majority of variants differentiating identical sequence type isolates mapped against accessory fragments, thus providing additional information to characterize strains. Although mlst still represents a reliable typing method, large-scale studies on molecular epidemiology and public health, as well as bacterial phylogenetics, population genetics and biosafety could benefit of a low cost and fast turnaround time approach such as the 2b-RAD analysis proposed here.

Pollutants bioavailability and toxicological risk from microplastics to marine mussels

Microplastics represent a growing environmental concern for the oceans due to their potential of adsorbing chemical pollutants, thus representing a still unexplored source of exposure for aquatic organisms. In this study polyethylene (PE) and polystyrene (PS) microplastics were shown to adsorb pyrene with a time and dose-dependent relationship. Results also indicated a marked capability of contaminated microplastics to transfer this model PAH to exposed mussels Mytilus galloprovincialis; tissue localization of microplastics occurred in haemolymph, gills and especially digestive tissues where a marked accumulation of pyrene was also observed. Cellular effects included alterations of immunological responses, lysosomal compartment, peroxisomal proliferation, antioxidant system, neurotoxic effects, onset of genotoxicity; changes in gene expression profile was also demonstrated through a new DNA microarray platform. The study provided the evidence that microplastics adsorb PAHs, emphasizing an elevated bioavailability of these chemicals after the ingestion, and the toxicological implications due to responsiveness of several molecular and cellular pathways to microplastics.

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.

Insights into molecular features of Venerupis decussata oocytes: a microarray-based study

The production of Venerupis decussata relies on wild seed collection, which has been recently compromised due to recruitment failure and severe mortalities. To address this issue and provide an alternative source of seed, artificial spawning and larval rearing programs were developed. However, hatchery-based seed production is a relatively new industry and it is still underdeveloped. A major hurdle in the European clam seed production is the control of spawning and reproduction, which is further hindered by the impossibility of obtaining fertile gametes by gonadal "stripping", as meiosis re-initiation is constrained to a maturation process along the genital ducts. In the present study, oocytes were collected from 15 females and microarray analyses was performed to investigate gene expression profiles characterizing released and stripped ovarian oocytes. A total of 198 differentially expressed transcripts between stripped and spawned oocytes were detected. Functional analysis carried out on these transcripts highlighted the importance of a few biological processes, which are most probably implicated in the control of oocyte competence. Significant differences were observed for transcripts encoding proteins involved in meiosis progression (e.g. dual specificity phosphatase CDC25), WNT signalling (e.g. frizzled class receptor 8, wingless-type MMTV integration site family member 4), steroid synthesis (e.g. progestin and adipoQ receptor family member 3, cytochrome P450-C17), mRNA processing (e.g. zinc finger protein XlCOF28), calcium regulation (e.g. regucalcin, calmodulin) and ceramide metabolism (ceramidase B, sphingomyelinase). This study provides new information on transcriptional profiles putatively associated with ovarian egg infertility, and suggests potential mechanisms regulating early oocyte development in clams. Genes which were differentially expressed between stripped and spawned oocytes might have a pivotal role during maturation process in the gonadal duct and could be interesting targets for further functional studies aiming to make ovarian oocytes fertilizable.

A microarray-based analysis of gametogenesis in two Portuguese populations of the European clam Ruditapes decussatus

The European clam, Ruditapes decussatus is a species with a high commercial importance in Portugal and other Southern European countries. Its production is almost exclusively based on natural recruitment, which is subject to high annual fluctuations. Increased knowledge of the natural reproductive cycle of R. decussatus and its molecular mechanisms would be particularly important in providing new highly valuable genomic information for better understanding the regulation of reproduction in this economically important aquaculture species. In this study, the transcriptomic bases of R. decussatus reproduction have been analysed using a custom oligonucleotide microarray representing 51,678 assembled contigs. Microarray analyses were performed in four gonadal maturation stages from two different Portuguese wild populations, characterized by different responses to spawning induction when used as progenitors in hatchery. A comparison between the two populations elucidated a specific pathway involved in the recognition signals and binding between the oocyte and components of the sperm plasma membrane. We suggest that this pathway can explain part of the differences in terms of spawning induction success between the two populations. In addition, sexes and reproductive stages were compared and a correlation between mRNA levels and gonadal area was investigated. The lists of differentially expressed genes revealed that sex explains most of the variance in gonadal gene expression. Additionally, genes like Foxl2, vitellogenin, condensing 2, mitotic apparatus protein p62, Cep57, sperm associated antigens 6, 16 and 17, motile sperm domain containing protein 2, sperm surface protein Sp17, sperm flagellar proteins 1 and 2 and dpy-30, were identified as being correlated with the gonad area and therefore supposedly with the number and/or the size of the gametes produced.

Deep transcriptome sequencing of Pecten maximus hemocytes: a genomic resource for bivalve immunology

Pecten maximus, the king scallop, is a bivalve species with important commercial value for both fisheries and aquaculture, traditionally consumed in several European countries. Major problems in larval rearing, however, still limit hatchery-based seed production. High mortalities during early larval stages, likely related to bacterial pathogens, represent the most relevant bottleneck. To address this issue, understanding host defense mechanisms against microbes is extremely important. In this study next-generation RNA-sequencing was carried on scallop hemocytes. To enrich for immune-related transcripts, cDNA libraries from hemocytes challenged in vivo with inactivated-Vibrio anguillarum and in vitro with pathogen-associated molecular patterns, as well as unchallenged controls, were sequenced yielding 216,444,674 sequence reads. De novo assembly of the scallop hemocyte transcriptome consisted of 73,732 contigs (31% annotated). A total of 934 contigs encoded proteins with a known immune function, grouped into several functional categories. Particular attention was reserved to Toll-like receptors (TLRs), a family of pattern recognition receptors (PRRs) involved in non-self recognition. Through mining the scallop hemocyte transcriptome, at least four TLRs could be identified. The organization of canonical TLR domains demonstrated that single cysteine cluster and multiple cysteine cluster TLRs co-exist in this species. In addition, preliminary data concerning their mRNA level following bacterial challenge suggested that different members of this family could exhibit opposite responses to pathogenic stimuli. Finally, a global analysis of differential expression comparing gene-expression levels in in vitro and in vivo stimulated hemocytes against controls provided evidence on a large set of transcripts involved in the great scallop immune response.

Gene transcription and biomarker responses in the clam Ruditapes philippinarum after exposure to ibuprofen

Pharmaceuticals are a class of emerging environmental contaminants that continuously enter aquatic environments. Presently, little information is available about the effects of these substances on non-target organisms, such as bivalves. We investigated the effects of ibuprofen (IBU) on the clam Ruditapes philippinarum. Clams were exposed for 1, 3, 5 and 7 days to 0, 100 and 1000 μgIBU/L, and established biomarker responses (haemolymph lysozyme, gill acetylcholinesterase and digestive gland superoxide dismutase activities) as well as digestive gland transcriptome were evaluated. A two-way ANOVA revealed significant effects of both "IBU concentration" and "exposure duration" on biomarker responses. Overall, the enzyme activities were generally lower in IBU-exposed clams than in controls. Although limited knowledge of the mollusc transcriptome makes it difficult to interpret the effects of IBU on clams, the gene transcription analysis using DNA microarrays enabled the identification of the putative molecular mode of action of the IBU. The functional analysis of differentially transcribed genes suggests that IBU can interfere with various signalling pathways in clams, such as arachidonic acid metabolism, apoptosis, peroxisomal proliferator-activated receptors, and nuclear factor-kappa B. In addition, several genes involved in the metabolism of xenobiotics (e.g., glutathione S-transferase, sulfotransferase, cytochrome P450) were also found to be significantly affected by IBU exposure. In summary, the integrated approach of gene transcription analysis and biomarker responses facilitated the elucidation of the putative mechanisms of action of IBU in non-target species.