Gene Expression Profiles of the Immuno-Transcriptome in Equine Asthma.

BACKGROUND: Mild equine asthma (MEA) and severe equine asthma (SEA) are two of the most frequent equine airway inflammatory diseases, but knowledge about their pathogenesis is limited. The goal of this study was to investigate gene expression differences in the respiratory tract of MEA- and SEA-affected horses and their relationship with clinical signs. METHODS: Clinical examination and endoscopy were performed in 8 SEA- and 10 MEA-affected horses and 7 healthy controls. Cytological and microbiological analyses of bronchoalveolar lavage (BAL) fluid were performed. Gene expression profiling of BAL fluid was performed by means of a custom oligo-DNA microarray. RESULTS: In both MEA and SEA, genes involved in the genesis, length, and motility of respiratory epithelium cilia were downregulated. In MEA, a significant overexpression for genes encoding inflammatory mediators was observed. In SEA, transcripts involved in bronchoconstriction, apoptosis, and hypoxia pathways were significantly upregulated, while genes involved in the formation of the protective muco-protein film were underexpressed. The SEA group also showed enrichment of gene networks activated during human asthma. CONCLUSIONS: The present study provides new insight into equine asthma pathogenesis, representing the first step in transcriptomic analysis to improve diagnostic and therapeutic approaches for this respiratory disease.

Long-lasting effects of chronic exposure to chemical pollution on the hologenome of the Manila clam.

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.

fshr: a fish sex-determining locus shows variable incomplete penetrance across flathead grey mullet populations.

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.

Potential for Genetic Improvement of Resistance to Perkinsus olseni in the Manila Clam, Ruditapes philippinarum, Using DNA Parentage Assignment and Mass Spawning.

The Manila clam Ruditapes philippinarum, a major cultured shellfish species, is threatened by infection with the microparasite Perkinsus olseni, whose prevalence increases with high water temperatures. Under the current trend of climate change, the already severe effects of this parasitic infection might rapidly increase the frequency of mass mortality events. Treating infectious diseases in bivalves is notoriously problematic, therefore selective breeding for resistance represents a key strategy for mitigating the negative impact of pathogens. A crucial step in initiating selective breeding is the estimation of genetic parameters for traits of interest, which relies on the ability to record parentage and accurate phenotypes in a large number of individuals. Here, to estimate the heritability of resistance against P. olseni, a field experiment mirroring conditions in industrial clam production was set up, a genomic tool was developed for parentage assignment, and parasite load was determined through quantitative PCR. A mixed-family cohort of potentially 1,479 clam families was produced in a hatchery by mass spawning of 53 dams and 57 sires. The progenies were seeded in a commercial clam production area in the Venice lagoon, Italy, where high prevalence of P. olseni had previously been reported. Growth and parasite load were monitored every month and, after 1 year, more than 1,000 individuals were collected for DNA samples and phenotype recording. A pooled sequencing approach was carried out using DNA samples from the hatchery broodstock and from a Venice lagoon clam population, providing candidate markers used to develop a 245-SNP panel. Parentage assignment for 246 F1 individuals showed sire and dam representation were high (75 and 85%, respectively), indicating a very limited risk of inbreeding. Moderate heritability (0.23 ± 0.11-0.35 ± 0.13) was estimated for growth traits (shell length, shell weight, total weight), while parasite load showed high heritability, estimated at 0.51 ± 0.20. No significant genetic correlations were found between growth-associated traits and parasite load. Overall, the preliminary results provided by this study show high potential for selecting clams resistant to parasite load. Breeding for resistance may help limit the negative effects of climate change on clam production, as the prevalence of the parasite is predicted to increase under a future scenario of higher temperatures. Finally, the limited genetic correlation between resistance and growth suggests that breeding programs could incorporate dual selection without negative interactions.

Integrated analysis of transcriptome, methylome and copy number aberrations data of marginal zone lymphoma and follicular lymphoma in dog.

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.

Cyclin D1 immunohistochemical expression and somatic mutations in canine oral melanoma.

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.

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.

Long Non-Coding RNAs as Molecular Signatures for Canine B-Cell Lymphoma Characterization.

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.

Mutational landscape of canine B-cell lymphoma profiled at single nucleotide resolution by RNA-seq.

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.

Tracing seafood at high spatial resolution using NGS-generated data and machine learning: Comparing microbiome versus SNPs.

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.

Longitudinal transcriptomic and genetic landscape of radiotherapy response in canine melanoma.

Canine malignant melanoma (MM) is a highly aggressive tumour 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 MM at diagnosis and at disease progression to identify possible genetic mechanisms associated with therapy failure. According to previous studies, a genetic similarity between canine MM and its human counterpart was observed. Several somatic mutations were functionally related to MAPK, PI3K/AKT and p53 signalling 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 tumour. 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.

Array Comparative Genomic Hybridization Analysis Reveals Significantly Enriched Pathways in Canine Oral Melanoma.

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.

Validation of epigenetic mechanisms regulating gene expression in canine B-cell lymphoma: An in vitro and in vivo approach.

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.

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.

DNA methylation profiling reveals common signatures of tumorigenesis and defines epigenetic prognostic subtypes of canine Diffuse Large B-cell Lymphoma.

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.

Transcriptome analysis of the brain of the sea bream (Sparus aurata) after exposure to human pharmaceuticals at realistic environmental concentrations.

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.

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.

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.

Ontogenetic onset of immune-relevant genes in the common sole (Solea solea).

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.