Search

Chronic exposure to pollutants affects natural populations, creating specific molecular and biochemical signatures. In the present study, we tested the hypothesis that chronic exposure to pollutants might have substantial effects on the Manila clam hologenome long after removal from contaminated sites. To reach this goal, a highly integrative approach was implemented, combining transcriptome, genetic and microbiota analyses with the evaluation of biochemical and histological profiles of the edible Manila clam Ruditapes philippinarum, as it was transplanted for 6 months from the polluted area of Porto Marghera (PM) to the clean area of Chioggia (Venice lagoon, Italy). One month post‐transplantation, PM clams showed several modifications to its resident microbiota, including an overrepresentation of the opportunistic pathogen Arcobacter spp. This may be related to the upregulation of several immune genes in the PM clams, potentially representing a host response to the increased abundance of deleterious bacteria. Six months after transplantation, PM clams demonstrated a lower ability to respond to environmental/physiological stressors related to the summer season, and the hepatopancreas‐associated microbiota still showed different compositions among PM and CH clams. This study confirms that different stressors have predictable effects in clams at different biological levels and demonstrates that chronic exposure to pollutants leads to long‐lasting effects on the animal hologenome. In addition, no genetic differentiation between samples from the two areas was detected, confirming that PM and CH clams belong to a single population. Overall, the obtained responses were largely reversible and potentially related to phenotypic plasticity rather than genetic adaptation. The results here presented will be functional for the assessment of the environmental risk imposed by chemicals on an economically important bivalve species.

The gilthead sea bream (Sparus aurata) is a marine fish of great importance for Mediterranean aquaculture. This species has long been considered resistant to Nervous Necrosis Virus (NNV), an RNA virus that causes massive mortalities in several farmed fish animals. However, the recent appearance of RGNNV/SJNNV reassortant strains started to pose a serious threat to sea bream hatcheries, as it is able to infect larvae and juveniles of this species. While host response to NNV has been extensively studied in adult fish, little attention has been devoted to early life history stages, which are generally the most sensitive ones. Here we report for the first time a time-course RNA-seq analysis on 21-day old fish gilthead sea bream larvae experimentally infected with a RGNNV/SJNNV strain. NNV-infected and mock-infected samples were collected at four time points (6 h, 12 h, 24 h, and 48 h post infection). Four biological replicates, each consisting of five pooled larvae, were analysed for each time point and group. A large set of genes were found to be significantly regulated, especially at early time points (6 h and 12 h), with several heat shock protein encoding transcripts being up-regulated (e.g. hspa5, dnaj4, hspa9, hsc70), while many immune genes were down-regulated (e.g. myd88 and irf5 at T06, pik3r1, stat3, jak1, il12b and il6st at T12). A gene set enrichment analysis (GSEA) identified several altered pathways/processes. For instance, the formation of peroxisomes, which are important anti-viral components as well as essential for nervous system homeostasis, and the autophagy pathway were down-regulated at 6 h and 24 h post infection (hpi). Finally, two custom “reactomes” (i.e. significant gene sets observed in other studies) were defined and used. The first reactome integrated the transcriptomic response to NNV in different fish species, while the second one included all genes found to be stimulated either by interferon (IFN) or by IFN and Chikungunya virus in zebrafish. Genes in both reactomes showed predominant up-regulation at 6hpi and 12hpi and a general down-regulation at 24hpi. Such evidence suggest a certain degree of similarity between the response of sea bream and that of other fish species to NNV, while the observed down-regulation of IFN- and viral-stimulated pathways argues for a possible interference of NNV against the host response.

Summary Whole-genome sequencing data were produced from a single flathead grey mullet female and assembled into a draft genome sequence, whereas publicly available sequence data were used to obtain a male draft sequence. Two pools, each consisting of 60 unrelated individuals, respectively, of male and female fish were analyzed using Pool-Sequencing. Mapping and analysis of Pool-Seq data against the draft genome(s) revealed >30 loci potentially associated with sex, the most promising locus of which, encoding the follicle-stimulating hormone receptor (fshr) and harboring two missense variants, was genotyped on 245 fish from four Mediterranean populations. Genotype data showed that fshr represents a previously unknown sex-determining locus, although the incomplete association pattern between fshr genotype and sex-phenotype, the variability of such pattern across different populations, and the presence of other candidate loci reveal that a greater complexity underlies sex determination in the flathead grey mullet., Graphical abstract, Highlights • Three SNPs with sex-specific allele frequencies found in the grey mullet fshr gene • Two missense variants at the fshr locus are significantly associated with male sex • fshr male genotypes show divergent distribution across geographical populations • Incomplete penetrance of fshr variants suggests that additional factors exist, Genetics; Evolutionary Biology; Genomics

Canine oral melanoma (COM) is the most frequent tumour with oral localization in dogs. Copy number gains and amplifications of CCND1, a gene coding for Cyclin D1, are the most frequent chromosomal aberrations described in human non-UV induced melanomas. Twenty-eight cases of COM were retrieved from paraffin-blocks archives. A total of 4 μm thick sections were immunostained with an antibody against human Cyclin D1 and Ki-67. Cyclin D1 and Ki-67 expressions were scored through two counting methods. DNA was extracted from 20 μm thick sections of formalin-fixed paraffin-embedded blocks. Pathological and surrounding healthy tissue was extracted independently. Cyclin D1 immunolabelling was detected in 69% (18/26) while Ki-67 was present in 88.5% (23/26) of cases. Statistical analysis revealed correlation between two counting methods for Cyclin D1 (r = 0.54; P = .004) and Ki-67 (r = 0.56; P = .003). The correlation found between Ki-67 and Cyclin D1 indexes in 16/26 cases labelled by both antibodies (r = 0.7947; P = .0002) suggests a possible use of Cyclin D1 index as prognostic marker. Polymerase chain reaction analysis on CCND1 coding sequence revealed the presence of nine somatic mutations in seven samples producing synonymous, missense and stop codons. Since none of the single-nucleotide polymorphisms was found to be recurrent, it is suggested that overexpression of Cyclin D1 may be the consequence of alterations of CCND1 upstream regions or other genetic aberrations not detectable with the methodology used in this study. Future studies are needed to verify the potential use of Cyclin D1 index as prognostic indicator and to highlight the molecular events responsible for Cyclin D1 overexpression in COMs.

Whole-genome sequence data were produced from a single flathead grey mullet female and assembled into a draft genome sequence, whereas publicly available sequence data were used to obtain a male draft sequence. Two pools, each consisting of 60 unrelated individuals respectively of male and female fish were analysed using Pool-Sequencing. Mapping and analysis of Pool-Seq data against the draft genome(s) revealed >30 loci potentially associated with sex. The most promising of these loci encoding the follicle stimulating hormone receptor (FSHR) and harbouring two missense variants, was genotyped on 245 fish from four Mediterranean populations. Genotype data showed that FSHR represents a previously unknown sex-determining locus. However, the incomplete association pattern between FSHR genotype and sex-phenotype, the variability of such pattern across different natural populations and the presence of other candidate loci suggest a greater complexity underling sex determination in the flathead grey mullet.

Human Mucosal Melanoma (hMM) is an aggressive neoplasm of neuroectodermal origin with distinctive features from the more common cutaneous form of malignant melanoma (cMM). At the molecular level, hMMs are characterized by large chromosomal aberrations rather than single-nucleotide mutations, typically observed in cMM. Given the scarcity of available cases, there have been many attempts to establish a reliable animal model. In pet dogs, Canine Oral Melanoma (COM) is the most common malignant tumor of the oral cavity, sharing clinical and histological aspects with hMM. To improve the knowledge about COM's genomic DNA alterations, in the present work, formalin-fixed, paraffin-embedded (FFPE) samples of COM from different European archives were collected to set up an array Comparative Genomic Hybridization (aCGH) analysis to estimate recurrent Copy Number Aberrations (CNAs). DNA was extracted in parallel from tumor and healthy fractions and 19 specimens were successfully submitted to labeling and competitive hybridization. Data were statistically analyzed through GISTIC2.0 and a pathway-enrichment analysis was performed with ClueGO. Recurrent gained regions were detected, affecting chromosomes CFA 10, 13 and 30, while lost regions involved chromosomes CFA 10, 11, 22, and 30. In particular, CFA 13 showed a whole-chromosome gain in 37% of the samples, while CFA 22 showed a whole-chromosome loss in 25%. A distinctive sigmoidal trend was observed in CFA 10 and 30 in 25 and 30% of the samples, respectively. Comparative analysis revealed that COM and hMM share common chromosomal changes in 32 regions. MAPK- and PI3K-related genes were the most frequently involved, while pathway analysis revealed statistically significant perturbation of cancer-related biological processes such as immune response, drug metabolism, melanocytes homeostasis, and neo-angiogenesis. The latter is a new evidence of a significant involvement of neovascularization-related pathways in COMs and can provide the rationale for future application in anti-cancer targeted therapies.

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.

Background: Diffuse large B-cell lymphoma (DLBCL), marginal zone lymphoma (MZL) and follicular lymphoma (FL) are the most common B-cell lymphomas (BCL) in dogs. Recent investigations have demonstrated overlaps of these histotypes with the human counterparts, including clinical presentation, biologic behavior, tumor genetics, and treatment response. The molecular mechanisms that underlie canine BCL are still unknown and new studies to improve diagnosis, therapy, and the utilization of canine species as spontaneous animal tumor models are undeniably needed. Recent work using human DLBCL transcriptomes has suggested that long non-coding RNAs (lncRNAs) play a key role in lymphoma pathogenesis and pinpointed a restricted number of lncRNAs as potential targets for further studies. Results: To expand the knowledge of non-coding molecules involved in canine BCL, we used transcriptomes obtained from a cohort of 62 dogs with newly-diagnosed multicentric DLBCL, MZL and FL that had undergone complete staging work-up and were treated with chemotherapy or chemo-immunotherapy. We developed a customized R pipeline performing a transcriptome assembly by multiple algorithms to uncover novel lncRNAs, and delineate genome-wide expression of unannotated and annotated lncRNAs. Our pipeline also included a new package for high performance system biology analysis, which detects high-scoring network biological neighborhoods to identify functional modules. Moreover, our customized pipeline quantified the expression of novel and annotated lncRNAs, allowing us to subtype DLBCLs into two main groups. The DLBCL subtypes showed statistically different survivals, indicating the potential use of lncRNAs as prognostic biomarkers in future studies. Conclusions: In this manuscript, we describe the methodology used to identify lncRNAs that differentiate B-cell lymphoma subtypes and we interpreted the biological and clinical values of the results. We inferred the potential functions of lncRNAs to obtain a comprehensive and integrative insight that highlights their impact in this neoplasm.

Marginal zone lymphoma (MZL) and follicular lymphoma (FL) are classified as indolent B-cell lymphomas in dogs. Aside from the clinical and histopathological similarities with the human counterpart, the molecular pathogenesis remains unclear. We integrated transcriptome, genome-wide DNA methylation and copy number aberration analysis to provide insights on the pathogenesis of canine MZL (n = 5) and FL (n = 7), also comparing them with diffuse large B-cell lymphoma (DLBCL). Transcriptome profiling highlighted the presence of similar biological processes affecting both histotypes, including BCR and TLR signalling pathways. However, FLs showed an enrichment of E2F targets, whereas MZLs were characterized by MYC-driven transcriptional activation signatures. FLs showed a distinctive loss on chr1 containing CEACAM23 and 24, conversely MZLs presented multiple recurrent gains on chr13, where MYC is located. The distribution of methylation peaks was similar between the two histotypes. Integrating data from the three omics, FLs resulted clearly separated from MZLs and DLBCL dataset. MZLs showed the enrichment of FoxM1 network and TLR associated TICAM1-dependent IRFs activation pathway. However, no specific signatures differentiated MZLs from DLBCLs. In conclusion, our study presents the first comprehensive analysis of molecular and epigenetic pathogenesis of canine FL and MZL.

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., Bargelloni et al. assemble the genome of the icefish Chionodraco myersi and use comparative analysis to examine the basis of adaptation to extreme Southern Ocean conditions. They show evidence for the evolutionary loss of globin genes and the peculiar regulation of genes involved in mitochondrial function.

Background. Diffuse large B-cell lymphoma (DLBCL) is the commonest lymphoma in both humans and dogs. Canine DLBCL (cDLBCL) is considered an ideal comparative model for drug development, but a complete genomic characterization of this tumor is still lacking. In this study, we report an integrated analysis to comprehensively define the molecular mechanisms of cDLBCL and possible associations with clinical outcome. Methods. Fifty cDLBCLs were analyzed by RNA-Seq, methyl-CpG-binding sequencing and array comparative genomic hybridization. Normal B-cells derived from lymph nodes of 11 healthy dogs were used as controls.Additionally, immunohistochemistry, in vitroand in vivoexperiments were performed as validation analyses. Results.Compared to normal B-cells, cDLBCL showed a marked up-regulation of genes involved in the PI3K/mTOR and NF-κB pathways, including several TLRs in association with MYD88, indicating mechanisms similar to the human activated B cell-like subtype DLBCL. Both RNA-Seq and methylation sequencing led to the identification of two groups of cDLBCLs bearing different clinical outcome. The two groups did not overlap with the human germinal center B-cell (GCB) and the activated B-cell-like (ABC) DLBCL subtypes or the human DLBCL consensus clusters. The dogs with the poorest outcome presented a signature largely defined by markers of T-cell-mediated immune responses, with a high expression of PDL-1, PD-1 and CTLA-4, also validated in an independent cohort of cDLBCL by immunohistochemistry. These data provide a strong rationale for the use of cDLBCL to study immune checkpoint modulators. The observed high expression of PI3K/mTOR pathway genes was confirmed and validated achieving a clear anti-tumor activity with the use of the PI3K-delta inhibitor idelalisib and of the novel dual PI3K/mTOR inhibitor bimiralisib in the cDLBCL cell line CLBL-1. The cDLBCLs showed an up-regulation of MYC and of its targets, sustained by recurrent gains in the chromosome 13, where the oncogene is located, in approximately half of the cases. Thus, we have exposed the cDLBCL cell line CLBL-1 to the BET inhibitor birabresib (OTX015) and to the BRD4 degrader MZ1. Both compounds caused a significant reduction in the proliferation of tumor cells, and this effect was stronger especially with the second compound. Exposure to MZ1 determined an important downregulation of MYC and also of LIN28B, the most overexpressed transcript in cDLBCL when compared to controls. While LIN28B does not seem to be a relevant gene for human DLBCL, its overexpression causes murine T-cell lymphomas (Beachy et al, Blood 2011), and there is a direct association of MYC with LIN28B promoter resulting in transcriptional transactivation (Chang et al, PNAS 2009). Here, LIN28B genetic silencing in the CLBL-1 lead to a reduction in cell growth, opening new therapeutic target perspectives in canine lymphoma. Conclusions. We have reported the first large next generation sequencing study investigating the cDLBCL transcriptome, methylome and the genome-wide CNVs. We identified deregulated pathways and individual transcripts providing therapeutic targets, including an immune-related signature affecting the outcome of a subgroup of cDLBCL. Our data sustain the use of cDLBCL as comparative models for human DLBCL but also highlight differences that must be kept in consideration. Disclosures Hillmann: PIQUR Therapeutics AG: Employment. Wymann:PIQUR Therapeutics AG: Employment, Equity Ownership, Patents & Royalties.

Canine malignant melanoma is a highly aggressive tumor with a low survival rate and represents an ideal spontaneous model for the human counterpart. Considerable progress has been recently obtained, but the therapeutic success for canine melanoma is still challenging. Little is known about the mechanisms beyond pathogenesis and melanoma development, and the molecular response to radiotherapy has never been explored before. A faster and deeper understanding of cancer mutational processes and developing mechanisms are now possible through next generation sequencing technologies. In this study, we matched whole exome and transcriptome sequencing in four dogs affected by malignant melanoma at diagnosis and at disease progression to identify possible genetic mechanisms associated with therapy failure. According to previous studies, a genetic similarity between canine malignant melanoma and its human counterpart was observed. Several somatic mutations were functionally related to MAPK, PI3K/AKT and p53 signaling pathways, but located in genes other than BRAF, RAS and KIT. At disease progression, several mutations were related to therapy effects. Natural killer cell-mediated cytotoxicity and several immune-system-related pathways resulted activated opening a new scenario on the microenvironment in this tumor. In conclusion, this study suggests a potential role of the immune system associated to radiotherapy in canine melanoma, but a larger sample size associated with functional studies are needed. This article is protected by copyright. All rights reserved.

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.

WOS:000442900300025; International audience; Among reared bivalves, some "novel species", such as the great scallop, Pecten maximus, have experienced more difficulty with routine reproduction due to their high sensitivity to biological, chemical, and physical stress during stages of early development. Working with high larval densities requires the use of aeration systems to provide optimal larval suspension and feed distribution. The high susceptibility of the great scallop to aeration in small-volume systems may impose an important limitation in hatchery-based practices. The present study aimed to investigate the processes impacted by aeration in P. maximus veliger larvae exposed to continuous aeration in small-volume tanks (10 L). Aeration appeared as major stressor that was responsible for early mortality among exposed animals (at 96 h after aeration started, haa) when exposure started 13 days after fertilization. Exposed larvae and controls were collected at 12, 24 and 72 haa, and a total of 18 cDNA libraries, each representing a pool of approximately 7,500 larvae, were sequenced, obtaining 358,817,016 raw reads. RNA-seq data were first used to build a de novo transcriptome assembly, and differential transcript abundance was assessed in exposed and control groups; thus, the molecular mechanisms involved in high sensitivity to aeration were deciphered. More than 2,000 transcripts were differentially expressed between exposed and control larvae across the entire time series (logFC \textgreater 1, FDR \textless 5%). Functional analysis revealed that transcriptional changes in larvae exposed to aeration mainly involved the genes that regulate digestive activity and energy metabolism, immune defense, inflammation, apoptosis, larval growth, and development. The results of this study demonstrate that, overall, aeration affects the feeding capacity and energy metabolism of larvae, with expected consequences on the animal's fitness, including its swimming efficiency. Aeration also triggered immune responses and apoptosis, which then increased through opportunistic infections. Notably, infections may be a consequence of a bacterial bloom triggered by the first mortality events that occurred in the culture. This study provides insights into the interactions between environmental variables and great scallop larvae physiology, and the results may contribute to the development of strategies for improving larval rearing practices and ensuring long-term sustainability of P. maximus aquaculture.

Developing reliable tools to trace food origin represents a major goal for producers and control authorities. Here, we test the hypothesis whether NGS-generated data could provide a reliable tool to ensure seafood traceability. As a test case, we used the Manila clam Ruditapes philippinarum, a bivalve mollusk of high commercial interest with worldwide distribution, collected in the Venice lagoon sites subjected to prohibition of clam harvesting because of chemical contamination as well as in authorized clam harvesting areas. The results obtained demonstrated that the geographic origin of Manila clam may be more accurately determined basing on microbiome data than single nucleotide polymorphisms. In particular, combining microbiome data with machine-learning techniques, we provide the experimental evidence that it is possible to trace the clam place of origin at high spatial resolution. Considering its low cost and portability, NGS-analysis of microbiome data might represent a cost-effective, high-resolution tool for reliable food traceability.

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

Despite canine B-cell Lymphoma (BCL) representing the most common haematological tumour, epigenetic events driving development and progression are scarcely known. Recently, canine Diffuse Large BCL (DLBCL) DNA methylome by genome-wide CpG microarray has identified genes and pathways associated to pathogenesis. To validate data previously obtained by array analysis, the CLBL-1 cell line was used and the HOXD10, FGFR2, ITIH5 and RASAL3 genes were selected. CLBL-1 cells were treated with two hypomethylating drugs (HDs; IC50, 50% inhibitory concentration), i.e. azacytidine and decitabine (DEC), either alone or in combination with three histone deacetylase inhibitors (HDACis; IC20), i.e. valproic acid, trichostatin and vorinostat. Following the incubation with both HDs, an overall decrease of promoter methylation was highlighted, thus confirming target genes hypermethylation. The highest mRNA restoration was observed following the exposure to HDs combined with HDACis, and mostly with valproic acid. Contrasting results were only obtained for RASAL3. An in vivo confirmation was finally attempted treating Nod-Scid mice engrafted with CLBL-1 cells with DEC. Although DEC did not arrest tumour growth, target genes promoter methylation was significantly reduced in DEC-treated mice vs controls. Overall, this work demonstrates that CLBL-1 cell line represents a reliable in vitro model to validate the methylation-dependent silencing of key genes for BCL; moreover, it may be useful for xenograft models in mice, despite its aggressive behaviour. In future, functional studies will be performed to deepen the role of selected genes on BCL pathogenesis and progression, and their methylation-dependent mechanism of regulation.

International audience; 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 mu 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 mu 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.

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. European Union KBBE.2013.1.2-10 European Community 311920 Fondazione Cassa di Risparmio Padova e Rovigo FCT - Foundation for Science and Technology research grant SPARCOMP under the Call ARISTEIA I of the National Strategic Reference Framework - by the EU 36 Hellenic Republic through the European Social Fund info:eu-repo/semantics/publishedVersion

The genomic landscape in human B-cell lymphoma has revealed several somatic mutations and potentially relevant germline alterations affecting therapy and prognosis. Also, mutations originally described as somatic aberrations have been shown to confer cancer predisposition when occurring in the germline. The relevance of mutations in canine B-cell lymphoma is scarcely known and gene expression profiling has shown similar molecular signatures among different B-cell histotypes, suggesting other biological mechanisms underlining differences. Here, we present a highly accurate approach to identify single nucleotide variants (SNVs) in RNA-seq data obtained from 62 completely staged canine B-cell lymphomas and 11 normal B-cells used as controls. A customized variant discovery pipeline was applied and SNVs were found in tumors and differentiated for histotype. A number of known and not previously identified SNVs were significantly associated to MAPK signaling pathway, negative regulation of apoptotic process and cell death, B-cell activation, NF-kB and JAK-STAT signaling. Interestingly, no significant genetic fingerprints were found separating diffuse large B-cell lymphoma from indolent lymphomas suggesting that differences of genetic landscape are not the pivotal causative factor of indolent behavior. We also detected several variants in expressed regions of canine B-cell lymphoma and identified SNVs having a direct impact on genes. Using this brand-new approach the consequence of a gene variant is directly associated to expression. Further investigations are in progress to deeply elucidate the mechanisms by which altered genes pathways may drive lymphomagenesis and a higher number of cases is also demanded to confirm this evidence.

Human pharmaceuticals such as Acetaminophen, Atenolol and Carbamazepine are pseudo persistent aquatic pollutants with yet unknown sub-lethal effects at environmentally relevant concentrations. Gilthead seabream (Sparus aurata) were exposed to Acetaminophen: 31.90 ± 11.07 μg L−1; Atenolol: 0.95 ± 0.38 μg L−1 and Carbamazepine: 6.95 ± 0.13 μg L−1 in a 28 day flow through experiment to (1) determine whether exposure to low concentrations in the μg·L−1 range of the pharmaceuticals alters the brain transcriptome and, (2) identify different expression profiles and treatment specific modes of action and pathways. Despite low exposure concentrations, 411, 7 and 612 differently expressed transcripts were identified in the individual treatments with Acetaminophen, Atenolol and Carbamazepine, respectively. Functional analyses of differentially expressed genes revealed a significant over representation of several biological processes, cellular compartment features and molecular functions for both Acetaminophen and Carbamazepine treatments. Overall, the results obtained in seabream brain suggest similar physiological responses to those observed in humans also at environmental concentrations, as well as the existence of treatment specific processes that may be useful for the development of biomarkers of contamination., This study was funded by a Marie Curie European Reintegration Grant (Proposal N° 239325-ERA4PHARM, FP7-PEOPLE-ERG-2008), the project SCARCE, funded by the Spanish Ministry for Sciences and Innovation (Consolider-Ingenio 2010 CSD2009-00065) and the project I2TEP for Cross Border Research and Technology Transfer between Spain and Portugal funded by the European Commission (0432_I2TEP_5_E Interreg EU). M.H. is supported by a Ramón y Cajal contract (RYC-2012-12217) from the Spanish Ministry of Economy and Competitiveness (MINECO).

Epigenetic deregulation is a hallmark of cancer characterized by frequent acquisition of new DNA methylation in CpG islands. To gain insight into the methylation changes of canine DLBCL, we investigated the DNA methylome in primary DLBCLs in comparison with control lymph nodes by genome-wide CpG microarray. We identified 1,194 target loci showing different methylation levels in tumors compared with controls. The hypermethylated CpG loci included promoter, 5′-UTRs, upstream and exonic regions. Interestingly, targets of polycomb repressive complex in stem cells were mostly affected suggesting that DLBCL shares a stem cell-like epigenetic pattern. Functional analysis highlighted biological processes strongly related to embryonic development, tissue morphogenesis and cellular differentiation, including HOX, BMP and WNT. In addition, the analysis of epigenetic patterns and genome-wide methylation variability identified cDLBCL subgroups. Some of these epigenetic subtypes showed a concordance with the clinical outcome supporting the hypothesis that the accumulation of aberrant epigenetic changes results in a more aggressive behavior of the tumor. Collectively, our results suggest an important role of DNA methylation in DLBCL where aberrancies in transcription factors were frequently observed, suggesting an involvement during tumorigenesis. These findings warrant further investigation to improve cDLBCL prognostic classification and provide new insights on tumor aggressiveness.

Phospholipidosis (PLD) is characterized by an intracellular accumulation of phospholipids in lysosomes and concurrent development of concentric lamellar bodies. It is induced in humans and in animals by drugs with a cationic amphiphilic structure. The purpose of the present study was to identify a set of molecular biomarkers of PLD in rat blood and heart, hypothetically applicable in preclinical screens within the drug development process. A toxicological study was set up in rats orally treated up to 11 days with 300 mg kg–1 per day–1 amiodarone (AMD). Light and transmission electron microscopy investigations were performed to confirm the presence of lamellar bodies indicative of phospholipid accumulation. The effects of AMD upon the transcriptome of these tissues were estimated using DNA microarray technology. Microarray data analysis showed that a total of 545 and 8218 genes were modulated by AMD treatment in heart and blood, respectively. Some genes implicated in the phospholipid accumulation in cells, such as phospholipase A2, showed similar alterations of gene expression. After transcriptome criteria of analysis and target selection, including also the involvement in the onset of PLD, 7 genes (Pla2g2a, Pla2g7, Gal, Il1b, Cebpb, Fcgr2b, Acer 2) were selected as candidate biomarkers of PLD in heart and blood tissues, and their potential usefulness as a sensitive screening test was screened and confirmed by quantitative Real-Time PCR analysis. Collectively, these data underscore the importance of transcriptional profiling in drug discovery and development, and suggest blood as a surrogate tissue for possible phospholipid accumulation in cardiomyocytes. Copyright © 2014 John Wiley & Sons, Ltd.

Posttransplant lymphoproliferative disorders (PTLDs) are a heterogeneous group of lymphoid proliferations that occur in the setting of depressed T-cell function due to immunosuppressive therapy used following solid organ transplantation, hematopoietic stem cell transplantation, and also xenotransplantation. In the present study, 28 immunosuppressed parkinsonian Macaca fascicularis were intracerebrally injected with wild-type or CTLA4-Ig transgenic porcine xenografts to identify a suitable strategy to enable long-term cell survival, maturation, and differentiation. Nine of 28 (32%) immunosuppressed primates developed masses compatible with PTLD, located mainly in the gastrointestinal tract and/or nasal cavity. The masses were classified as monomorphic PTLD according to the World Health Organization classification. Immunohistochemistry and polymerase chain reaction (PCR) analyses revealed that the PTLDs were associated with macaca lymphocryptovirus as confirmed by double-labeling immunohistochemistry for CD20 and Epstein-Barr nuclear antigen 2 (EBNA-2), where the viral protein was located within the CD20+ neoplastic B cells. In sera from 3 distinct phases of the experimental life of the primates, testing by quantitative PCR revealed a progression of the viral load that paralleled the PTLD progression and no evidence of zoonotic transmission of porcine lymphotropic herpesvirus through xenoneuronal grafts. These data suggest that monitoring the variation of macaca lymphocryptovirus DNA in primates could be used as a possible early diagnostic tool for PTLD progression, allowing preemptive treatment such as immunosuppression therapy reduction.

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.

In the present study we address the investigation of TLR2 evolutionary selection in two Antarctic teleosts, Trematomus bernacchii (Nototheniidae) and Chionodraco hamatus (Channichthyidae). The nucleotide sequence of TLR2 has been determined in both species, encoding 20 leucine-rich repeats (LRRs) in the extracellular region and a classical Toll/IL-1R (TIR) domain in the intracellular region. High expression level of T. bernacchii TLR2 was found in spleen and skin. Using different methods we identified six codons that underwent Darwinian selection while 20 were found to be negatively selected. Molecular models of C. hamatus and T. bernacchii TLR2 ectodomain as well as of the TIR domain were built by Homology Modeling. Molecular Dynamics simulations were performed in water for 15 ns. The sites under positive selection were residing on the convex side of the solenoid, four out of six were in a 35-residue-long region including the central/N-terminal domain boundary: two in the external loop of LRR11 and the other two in the LRR12 loop. This region has been demonstrated to be the functional site of ligand interaction in human TLR2 structure. Antarctic TLR2 models showed more flexibility than TLR2 from the temperate species Gasterosteus aculeatus. These results suggest that the selective pressure has shaped TLR2 molecule in such a way that increased its activity under the peculiar Antarctic environmental conditions.

In the last decade, a new gene family encoding non-rearranging receptors, called novel immune-type receptors (NITRs), has been discovered in teleost fish. NITRs belong to the immunoglobulin superfamily and represent an extraordinarily divergent and rapidly evolving gene complex. Genomic analysis of a region spanning 270 kb led to the discovery of a NITR gene cluster in the European sea bass (Dicentrarchus labrax). In total, 27 NITR genes and three putative pseudogenes, organised in a tandemly arrayed cluster, were identified. Sea bass NITR genes maintain the three major genomic organisations that appear to be essentially conserved among fish species along with new features presumably involving processes of intron loss, exon deletion and acquisition of new exons. Comparative and evolutionary analyses suggest that these receptors have evolved following a “birth-and-death” model of gene evolution in which duplication events together with lineage-specific gain and loss of individual members contributed to the rapid diversification of individual gene families. In this study, we demonstrate that species-specific gene expansions provide the raw material for diversifying, positive Darwinian selection favouring the evolution of a highly diverse array of molecules. Electronic supplementary material The online version of this article (doi:10.1007/s00251-009-0398-3) contains supplementary material, which is available to authorized users.

Fish are free-living organisms since initial stages of development and are exposed to numerous pathogens before their lymphoid organs have matured and adaptive immunity has developed. Susceptibility to diseases and juvenile mortality represent key critical factors for aquaculture. In this context, the characterization of the appearance kinetics of the immune system key members will be useful in understanding the ability of a particular species in generating immune protection against invading pathogens at different developmental stages. The present study characterized, for the first time, the transcriptional onset of un-explored relevant genes of both innate and adaptive immune system during the Solea solea ontogenesis. Gene expression profiles of immune relevant genes was investigated, by means of DNA microarray, in ten developmental stages, from hatching (1 day post-hatching, dph) to accomplishment of the juvenile form (33 dph). The obtained results revealed that transcripts encoding relevant members of innate immune repertoire, such as lysozyme, AMPs (hepcidin, β-defensin), PPRs and complement components are generally characterized by high expression levels at first stages (i.e. hatch and first feeding) indicating protection from environmental pathogens even at early development. Transcription of adaptive immune genes (i.e. Class I and class II MHC, TCRs) differs from that of the innate immune system. Their onset coincides with metamorphosis and larvae-to-juvenile transition, and likely overlaps with the appearance and maturation of the main lymphoid organs. Finally, data collected suggest that at the end of metamorphosis S. solea cell-mediated immune system hasn't still undergone full maturation.

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.

International audience; Coastal systems could be affected by hypoxic events brought about by global change. These areas are essential nursery habitats for several fish species including the common sole (Solea solea L.). Tolerance of fish to hypoxia depends on species and also on their physiological condition and nutritional status. Indeed, high dietary lipid content has been recently shown to negatively impact the resistance of sole to a severe hypoxic challenge. In order to study the molecular mechanisms involved in the early response to hypoxic stress, the present work examined the hepatic transcriptome in common sole fed diets with low and high lipid content, exposed to severe hypoxia. The activity of AMP-activated protein kinase (AMPK) was also investigated through the quantification of threonine-172 phosphorylation in the alpha subunit. The results show that hypoxia consistently regulates several actors involved in energy metabolism pathways and particularly AMPKα, as well as some involved in cell growth and maintenance or unfolded protein response. Our findings reveal that (1) the expression of genes involved in biological processes with high energy cost or implicated in aerobic ATP synthesis was down-regulated by hypoxia, contrary to genes involved in neoglucogenesis or in angiogenesis, (2) the consumption of high lipid induced regulation of metabolic pathways going against this energy saving, and (3) this control was fine-tuned by the regulation of several transcriptomic factors. These results provide insight into the biological processes involved in the hepatic response to hypoxic stress and underline the negative impact of high lipid consumption on the tolerance of common sole to hypoxia.

Hypoxia episodes have been generally studied over the last decade for their lethal consequences in coastal areas, which are nurseries for several fish species. The possible long-lasting effects of non-lethal hypoxia exposures at larval stage are, however, poorly documented. We investigated the long-lasting phenotypic impact of an experimental exposure to moderate hypoxia (40 % air saturation) between 30 and 38 days post-hatching in European sea bass (Dicentrarchus labrax). At the juvenile stage, our data indicate that this early hypoxia exposure impacted on juvenile growth rate but not on hypoxia tolerance. Growth depression was associated with a regulation of metabolism as suggested by a decrease in blood glycaemia, a decreased hepatic expression of the glucose transporter glut2 and of acaca involved in fatty acid synthesis as well as an increased hepatic expression of glycolytic ldh-a, which is a target gene of hypoxia inducible factors (HIFs). Interestingly, while fish were not exposed to hypoxia, the cellular oxygen sensor egln3 gene, which is also a target of HIFs, exhibited higher expression level in the hepatic tissue of juveniles that have experienced hypoxia at larval stage. Even if further studies are needed to evaluate the functional significance of such effects and to decipher the underlying molecular mechanisms, the present data evidence that early life exposure to a hypoxia event has long-lasting impact on sea bass juvenile physiology.

Background Xenotransplantation is a potential answer to the current organ shortage, but the risk of infections related to overimmunosuppression is an important parameter that may predict the recipient's long-term survival. Cytomegalovirus (CMV) in xenotransplanted and immunosuppressed primates is a well-known cause of disease particularly affecting the gastrointestinal (GI) tract and a zoonotic concern. Methods Post-mortem sera and tissues from 45 immunosuppressed and xenografted Macaca fascicularis were evaluated for CMV using antisera specific for the immediate early 1 (IE1), anti-RhCMV, and QPCR for virus. Results Serological analysis showed 100% positivity for the presence of CMV antibodies following the application of a specific test designed for RhCMV. Five of 45 primates showed typical lesions of CMV infection in the GI tract, including neutrophilic enteritis and inclusion bodies. Molecular analysis confirmed the presence of recipient's CMV in the tissues with CMV histopathology. Porcine CMV from the donor animals was not found in any of the CMV-specific IHC-positive recipients. Conclusion The presence of active CMV infection in animals intended for xenograft experiments can lead to severe gastrointestinal lesions that could impact the overall aims of the study. In such cases, the animals should be investigated using appropriate (non-human primate-specific) diagnostic tools prior to use and treated aggressively with state-of-the-art antiviral therapy.

Background Water temperature greatly influences the physiology and behaviour of teleost fish as other aquatic organisms. While fish are able to cope with seasonal temperature variations, thermal excursions outside their normal thermal range might exceed their ability to respond leading to severe diseases and death. Profound differences exist in thermal tolerance across fish species living in the same geographical areas, promoting for investigating the molecular mechanisms involved in susceptibility and resistance to low and high temperatures toward a better understanding of adaptation to environmental challenges. The gilthead sea bream, Sparus aurata, is particularly sensitive to cold and the prolonged exposure to low temperatures may lead to the "winter disease", a metabolic disorder that significantly affects the aquaculture productions along the Northern Mediterranean coasts during winter-spring season. While sea bream susceptibility to low temperatures has been extensively investigated, the cascade of molecular events under such stressful condition is not fully elucidated. Results In the present study two groups of wild sea bream were exposed for 21 days to two temperature regimes: 16 ± 0.3°C (control group) and 6.8 ± 0.3°C (cold-exposed group) and DNA microarray analysis of liver transcriptome was carried out at different time points during cold exposure. A large set of genes was found to be differentially expressed upon cold-exposure with increasingly relevant effects being observed after three weeks at low temperature. All major known responses to cold (i.e. anti-oxidant response, increased mitochondrial function, membrane compositional changes) were found to be conserved in the gilthead sea bream, while, evidence for a key role of unfolded protein response (UPR) to endoplasmic reticulum (ER) stress, during short- and long-term exposure to cold is reported here for the first time. Conclusions Transcriptome data suggest a scenario where oxidative stress, altered lipid metabolism, ATP depletion and protein denaturation converge to induce ER stress. The resulting UPR activation further promotes conditions for cell damage, and the inability to resolve ER stress leads to severe liver dysfunction and potentially to death. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-765) contains supplementary material, which is available to authorized users.

Platelet-derived growth factors (PDGFs) belong to a family of polypeptide growth factors that signal through cell surface tyrosine kinase receptors to stimulate growth, proliferation and differentiation. Platelet-derived growth factor receptors (PDGFRs) are also considered important targets for specific kinase inhibitors in the treatment of several human tumours. The aim of this study was to investigate the role of PDGF-A, PDGFB, PDGFR-alpha and PDGFR-beta in canine lymphoma by determining gene and protein expression in lymph nodes of dogs with diffuse large B-cell lymphoma (DLBCL), peripheral T-cell lymphoma (PTCL), T-Iymphoblastic lymphoma (T-LBL) and in healthy control dogs. One lymph node was also studied at the end of therapy in a subset of dogs in remission for DLBCL. In controls, PDGF-A, PDGFR-alpha and PDGFR-beta mRNA levels were significantly higher than in DLBCLs, PTCLs and T-LBLs. However, PDGFR-alpha and PDGFR-beta were minimally expressed by lymphocytes and plasma cells in normal lymph nodes as determined by immunohistochemistry, while neoplastic B and T cells showed the highest score (P

Canine Diffuse Large B-cell Lymphoma (cDLBCL) is an aggressive cancer with variable clinical response. Despite recent attempts by gene expression profiling to identify the dog as a potential animal model for human DLBCL, this tumor remains biologically heterogeneous with no prognostic biomarkers to predict prognosis. The aim of this work was to identify copy number aberrations (CNAs) by high-resolution array comparative genomic hybridization (aCGH) in 12 dogs with newly diagnosed DLBCL. In a subset of these dogs, the genetic profiles at the end of therapy and at relapse were also assessed. In primary DLBCLs, 90 different genomic imbalances were counted, consisting of 46 gains and 44 losses. Two gains in chr13 were significantly correlated with clinical stage. In addition, specific regions of gains and losses were significantly associated to duration of remission. In primary DLBCLs, individual variability was found, however 14 recurrent CNAs (>30%) were identified. Losses involving IGK, IGL and IGH were always found, and gains along the length of chr13 and chr31 were often observed (>41%). In these segments, MYC, LDHB, HSF1, KIT and PDGFR alpha are annotated. At the end of therapy, dogs in remission showed four new CNAs, whereas three new CNAs were observed in dogs at relapse compared with the previous profiles. One ex novo CNA, involving TCR, was present in dogs in remission after therapy, possibly induced by the autologous vaccine. Overall, aCGH identified small CNAs associated with outcome, which, along with future expression studies, may reveal target genes relevant to cDLBCL.