Serena Ferraresso

Development and validation of a gene expression oligo microarray for the gilthead sea bream (Sparus aurata)

BackgroundAquaculture represents the most sustainable alternative of seafood supply to substitute for the declining marine fisheries, but severe production bottlenecks remain to be solved. The application of genomic technologies offers much promise to rapidly increase our knowledge on biological processes in farmed species and overcome such bottlenecks. Here we present an integrated platform for mRNA expression profiling in the gilthead sea bream (Sparus aurata), a marine teleost of great importance for aquaculture.ResultsA public data base was constructed, consisting of 19,734 unique clusters (3,563 contigs and 16,171 singletons). Functional annotation was obtained for 8,021 clusters. Over 4,000 sequences were also associated with a GO entry. Two 60mer probes were designed for each gene and in-situ synthesized on glass slides using Agilent SurePrint™ technology. Platform reproducibility and accuracy were assessed on two early stages of sea bream development (one-day and four days old larvae). Correlation between technical replicates was always > 0.99, with strong positive correlation between paired probes. A two class SAM test identified 1,050 differentially expressed genes between the two developmental stages. Functional analysis suggested that down-regulated transcripts (407) in older larvae are mostly essential/housekeeping genes, whereas tissue-specific genes are up-regulated in parallel with the formation of key organs (eye, digestive system). Cross-validation of microarray data was carried out using quantitative qRT-PCR on 11 target genes, selected to reflect the whole range of fold-change and both up-regulated and down-regulated genes. A statistically significant positive correlation was obtained comparing expression levels for each target gene across all biological replicates. Good concordance between qRT-PCR and microarray data was observed between 2- and 7-fold change, while fold-change compression in the microarray was present for differences greater than 10-fold in the qRT-PCR.ConclusionA highly reliable oligo-microarray platform was developed and validated for the gilthead sea bream despite the presently limited knowledge of the species transcriptome. Because of the flexible design this array will be able to accommodate additional probes as soon as novel unique transcripts are available.

Effects of the total replacement of fish-based diet with plant-based diet on the hepatic transcriptome of two European sea bass (Dicentrarchus labrax) half-sibfamilies showing different growth rates with the plant-based diet.

BackgroundEfforts towards utilisation of diets without fish meal (FM) or fish oil (FO) in finfish aquaculture have been being made for more than two decades. Metabolic responses to substitution of fishery products have been shown to impact growth performance and immune system of fish as well as their subsequent nutritional value, particularly in marine fish species, which exhibit low capacity for biosynthesis of long-chain poly-unsaturated fatty acids (LC-PUFA). The main objective of the present study was to analyse the effects of a plant-based diet on the hepatic transcriptome of European sea bass (Dicentrarchus labrax).ResultsWe report the first results obtained using a transcriptomic approach on the liver of two half-sibfamilies of the European sea bass that exhibit similar growth rates when fed a fish-based diet (FD), but significantly different growth rates when fed an all-plant diet (VD). Overall gene expression was analysed using oligo DNA microarrays (GPL9663). Statistical analysis identified 582 unique annotated genes differentially expressed between groups of fish fed the two diets, 199 genes regulated by genetic factors, and 72 genes that exhibited diet-family interactions. The expression of several genes involved in the LC-PUFA and cholesterol biosynthetic pathways was found to be up-regulated in fish fed VD, suggesting a stimulation of the lipogenic pathways. No significant diet-family interaction for the regulation of LC-PUFA biosynthesis pathways could be detected by microarray analysis. This result was in agreement with LC-PUFA profiles, which were found to be similar in the flesh of the two half-sibfamilies. In addition, the combination of our transcriptomic data with an analysis of plasmatic immune parameters revealed a stimulation of complement activity associated with an immunodeficiency in the fish fed VD, and different inflammatory status between the two half-sibfamilies. Biological processes related to protein catabolism, amino acid transaminations, RNA splicing and blood coagulation were also found to be regulated by diet, while the expression of genes involved in protein and ATP synthesis differed between the half-sibfamilies.ConclusionsOverall, the combined gene expression, compositional and biochemical studies demonstrated a large panel of metabolic and physiological effects induced by total substitution of both FM and FO in the diets of European sea bass and revealed physiological characteristics associated with the two half-sibfamilies.

Skin healing and scale regeneration in fed and unfed sea bream, Sparus auratus

BackgroundFish scales are an important reservoir of calcium and phosphorus and together with the skin function as an integrated barrier against environmental changes and external aggressors. Histological studies have revealed that the skin and scales regenerate rapidly in fish when they are lost or damaged. In the present manuscript the histological and molecular changes underlying skin and scale regeneration in fed and fasted sea bream (Sparus auratus) were studied using a microarray 3 and 7 days after scale removal to provide a comprehensive molecular understanding of the early stages of these processes.ResultsHistological analysis of skin/scales revealed 3 days after scale removal re-epithelisation and formation of the scale pocket had occurred and 53 and 109 genes showed significant up or down-regulation, respectively. Genes significantly up-regulated were involved in cell cycle regulation, cell proliferation and adhesion, immune response and antioxidant activities. 7 days after scale removal a thin regenerated scale was visible and only minor changes in gene expression occurred. In animals that were fasted to deplete mineral availability the expression profiles centred on maintaining energy homeostasis. The utilisation of fasting as a treatment emphasised the competing whole animal physiological requirements with regard to barrier repair, infection control and energy homeostasis.ConclusionsThe identification of numerous genes involved in the mitotic checkpoint and cell proliferation indicate that the experimental procedure may be useful for understanding cell proliferation and control in vertebrates within the context of the whole animal physiology. In response to skin damage genes of immune surveillance were up-regulated along with others involved in tissue regeneration required to rapidly re-establish barrier function. Additionally, candidate fish genes were identified that may be involved in cytoskeletal re-modelling, mineralization and stem cells, which are of potential use in aquaculture and fish husbandry, as they may impact on the ability of the fish to produce structural proteins, such as muscle, efficiently.

Development of an oligo DNA microarray for the European sea bass and its application to expression profiling of jaw deformity

BackgroundThe European sea bass (Dicentrarchus labrax) is a marine fish of great importance for fisheries and aquaculture. Functional genomics offers the possibility to discover the molecular mechanisms underlying productive traits in farmed fish, and a step towards the application of marker assisted selection methods in this species. To this end, we report here on the development of an oligo DNA microarray for D. labrax.ResultsA database consisting of 19,048 unique transcripts was constructed, of which 12,008 (63%) could be annotated by similarity and 4,692 received a GO functional annotation. Two non-overlapping 60mer probes were designed for each unique transcript and in-situ synthesized on glass slides using Agilent SurePrint™ technology. Probe design was positively completed for 19,035 target clusters; the oligo microarray was then applied to profile gene expression in mandibles and whole-heads of fish affected by prognathism, a skeletal malformation that strongly affects sea bass production. Statistical analysis identified 242 transcripts that are significantly down-regulated in deformed individuals compared to normal fish, with a significant enrichment in genes related to nervous system development and functioning. A set of genes spanning a wide dynamic range in gene expression level were selected for quantitative RT-PCR validation. Fold change correlation between microarray and qPCR data was always significant.ConclusionsThe microarray platform developed for the European sea bass has a high level of flexibility, reliability, and reproducibility. Despite the well known limitations in achieving a proper functional annotation in non-model species, sufficient information was obtained to identify biological processes that are significantly enriched among differentially expressed genes. New insights were obtained on putative mechanisms involved on mandibular prognathism, suggesting that bone/nervous system development might play a role in this phenomenon.

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

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

Exploring the larval transcriptome of the common sole ( Solea solea L.)

BackgroundThe common sole (Solea solea) is a promising candidate for European aquaculture; however, the limited knowledge of the physiological mechanisms underlying larval development in this species has hampered the establishment of successful flatfish aquaculture. Although the fact that genomic tools and resources are available for some flatfish species, common sole genomics remains a mostly unexplored field. Here, we report, for the first time, the sequencing and characterisation of the transcriptome of S. solea and its application for the study of molecular mechanisms underlying physiological and morphological changes during larval-to-juvenile transition.ResultsThe S. solea transcriptome was generated from whole larvae and adult tissues using the Roche 454 platform. The assembly process produced a set of 22,223 Isotigs with an average size of 726 nt, 29 contigs and a total of 203,692 singletons. Of the assembled sequences, 75.2% were annotated with at least one known transcript/protein; these transcripts were then used to develop a custom oligo-DNA microarray. A total of 14,674 oligonucleotide probes (60 nt), representing 12,836 transcripts, were in situ synthesised onto the array using Agilent non-contact ink-jet technology. The microarray platform was used to investigate the gene expression profiles of sole larvae from hatching to the juvenile form. Genes involved in the ontogenesis of the visual system are up-regulated during the early stages of larval development, while muscle development and anaerobic energy pathways increase in expression over time. The gene expression profiles of key transcripts of the thyroid hormones (TH) cascade and the temporal regulation of the GH/IGF1 (growth hormone/insulin-like growth factor I) system suggest a pivotal role of these pathways in fish growth and initiation of metamorphosis. Pre-metamorphic larvae display a distinctive transcriptomic landscape compared to previous and later stages. Our findings highlighted the up-regulation of gene pathways involved in the development of the gastrointestinal system as well as biological processes related to folic acid and retinol metabolism. Additional evidence led to the formation of the hypothesis that molecular mechanisms of cell motility and ECM adhesion may play a role in tissue rearrangement during common sole metamorphosis.ConclusionsNext-generation sequencing provided a good representation of the sole transcriptome, and the combination of different approaches led to the annotation of a high number of transcripts. The construction of a microarray platform for the characterisation of the larval sole transcriptome permitted the definition of the main processes involved in organogenesis and larval growth.

Genomic Prediction of Resistance to Pasteurellosis in Gilthead Sea Bream (Sparus aurata) Using 2b-RAD Sequencing.

Gilthead sea bream (Sparus aurata) is a species of paramount importance to the Mediterranean aquaculture industry, with an annual production exceeding 140,000 metric tons. Pasteurellosis due to the Gram-negative bacterium Photobacterium damselae subsp. piscicida (Phdp) causes significant mortality, especially during larval and juvenile stages, and poses a serious threat to bream production. Selective breeding for improved resistance to pasteurellosis is a promising avenue for disease control, and the use of genetic markers to predict breeding values can improve the accuracy of selection, and allow accurate calculation of estimated breeding values of nonchallenged animals. In the current study, a population of 825 sea bream juveniles, originating from a factorial cross between 67 broodfish (32 sires, 35 dams), were challenged by 30 min immersion with 1 × 105 CFU virulent Phdp. Mortalities and survivors were recorded and sampled for genotyping by sequencing. The restriction-site associated DNA sequencing approach, 2b-RAD, was used to generate genome-wide single nucleotide polymorphism (SNP) genotypes for all samples. A high-density linkage map containing 12,085 SNPs grouped into 24 linkage groups (consistent with the karyotype) was constructed. The heritability of surviving days (censored data) was 0.22 (95% highest density interval: 0.11–0.36) and 0.28 (95% highest density interval: 0.17–0.4) using the pedigree and the genomic relationship matrix respectively. A genome-wide association study did not reveal individual SNPs significantly associated with resistance at a genome-wide significance level. Genomic prediction approaches were tested to investigate the potential of the SNPs obtained by 2b-RAD for estimating breeding values for resistance. The accuracy of the genomic prediction models (r = 0.38–0.46) outperformed the traditional BLUP approach based on pedigree records (r = 0.30). Overall results suggest that major quantitative trait loci affecting resistance to pasteurellosis were not present in this population, but highlight the effectiveness of 2b-RAD genotyping by sequencing for genomic selection in a mass spawning fish species.

Exploring the effects of seasonality and chemical pollution on the hepatopancreas transcriptome of the Manila clam

The assessment of marine environmental health is a complex but fundamental task both for ecosystem conservation and food safety related to the human consumption of marine products. Manila clams inhabiting the Venice Lagoon constitute an excellent case study for evaluating the effects of complex mixtures of industrial and urban effluents on aquatic organisms. Clams were collected in different seasons at four locations within the Venice Lagoon. The sampling sites were characterized by a range of pollutant concentrations and included Porto Marghera, a highly polluted industrial area where clam harvesting for human consumption is strictly forbidden. Pooled soft tissues were subjected to mass spectroscopy analysis to measure the concentrations of PCDDs/PCDFs/PCBs‐DL, PCBs, PBDEs, HCB and PAHs, and pooled digestive gland samples were used for gene expression profiling. While seasonal variation was found to be responsible for the largest proportion of transcriptional changes, significance analysis of microarrays quantitative correlation analysis identified 162 transcripts that were correlated with at least one class of chemicals measured in the samples from the four different sampling sites. Prediction Analysis of Microarrays (PAM) identified a minimal set of seven genes that correctly assigned samples collected in the restricted polluted area (Porto Marghera), independent of the season in which they were collected. An integrated approach combining transcriptomics and chemical analyses of the Manila clam provided a global picture of how Manila clams respond to complex mixtures of xenobiotics and their interplay with other biotic and abiotic factors. We were also able to identify gene expression signatures for different classes of chemicals and a set of robust biomarkers of exposure to these chemicals.

Epigenetic Silencing of TFPI-2 in Canine Diffuse Large B-Cell Lymphoma

Epigenetic modifications are important early events during carcinogenesis. In particular, hypermethylation of CpG islands in the promoter region of tumor suppressor genes is a well-known mechanism of gene silencing that contributes to cancer development and progression. Tissue factor pathway inhibitor 2 (TFPI-2) is a tumor suppressor involved in invasiveness inhibition. Although TFPI-2 transcriptional silencing, through promoter hypermethylation, has been widely reported in several human malignancies, it has never been explored in lymphoma. In the present study TFPI-2 methylation and gene expression have been investigated in canine Diffuse Large B-cell lymphomas (cDLBCL). The methylation level of 23 CpGs located within the TFPI-2 promoter was investigated by bisulfite-specific PCR and next generation amplicon deep sequencing (GS Junior 454, Roche) in 22 cDLBCLs and 9 controls. For the same specimens, TFPI-2 gene expression was assessed by means of Real-time RT-PCR. Sequence analysis clearly demonstrated that TFPI2 is frequently hypermethylated in cDLBCL. Hypermethylation of the TFPI-2 promoter was found in 77% of DLBCLs (17 out of 22) and in one normal lymph node. Globally, dogs with DLBCL showed a mean methylation level significantly increased compared to controls (p<0.01) and analysis of hypermethylation by site identified 19 loci out of 23 (82%) with mean methylation levels from 2- to 120-fold higher in cDLBCL. Gene expression analysis confirmed a significant down-regulation of TFPI-2 (p<0.05) in DLBCLs compared with normal lymph nodes, suggesting that TFPI-2 hypermethylation negatively regulates its transcription. In addition, a significant positive correlation (p<0.01) was found between TFPI-2 methylation levels and age providing the first indication of age-associated epigenetic modifications in canine DLBCL. To conclude, our findings demonstrated that epigenetic dysregulation of TFPI-2, leading to its reduced expression, is frequently detected in canine DLBCL. In the next future, the aberrant TFPI-2 promoter hypermethylation may be considered in association with prognosis and therapy.

Real time RT-PCR analysis of inflammatory mediator expression in recurrent airway obstruction-affected horses

The goal of the present study was to investigate mRNA expression levels of several cytokines and inflammatory mediators in broncho-alveolar lavage (BAL) fluid and respiratory epithelium in recurrent airway obstruction (RAO)-affected horses. RAO, also called heaves, is a common, performance-limiting, equine respiratory disease with clinical signs and pathophysiological similarities to human asthma, and characterized by bronchospasm, neutrophilic infiltration and increased mucus in the airways. Six RAO-affected horses were examined twice within 15 days and seven clinically healthy horses were examined for comparison. Quantitative real-time RT-PCR was used to assess mRNA expression of the inflammatory mediators IL-1β, IL-6, IL-8, IL-13, IL-17, TNFα, INFγ, TGFβ1, NFκ-β and TRL4 in bronchial biopsies and in BAL fluid. Gene expression levels were then compared with clinical signs, endoscopic examination, complete blood cell count, cytology of BAL fluid, histological examination of bronchial tissue and bacteriological and mycological examinations. Expression of IL1β, IL8, TLR4, TNFα, TGFβ1 and NFkβ transcripts was significantly up-regulated in RAO-affected compared to healthy horses. A similar trend, albeit not significant, was showed for IL17 and INFγ. A highly significant correlation was observed among IL-1β, IL8, TGFβ1, NFkβ, TRL4, and INFγ expression patterns as well as between expression levels of these genes and clinical parameters. In the present study, the comparison between clinically healthy and RAO-affected horses gave new insights on the cytokine expression in equine health and disease status. The identification of cytokines implicated in the pathogenesis of RAO may contribute to the diagnosis and treatment of this disease.