Different approaches to identify markers for meat quality and production traits in pigs: analysis of the transcriptome in post mortem skeletal muscle and investigation of candidate genes

Due to the growing attention of consumers towards their food, improvement of quality of animal products has become one of the main focus of research. To this aim, the application of modern molecular genetics approaches has been proved extremely useful and effective. This innovative drive includes all livestock species productions, including pork. The Italian pig breeding industry is unique because needs heavy pigs slaughtered at about 160 kg for the production of high quality processed products. For this reason, it requires precise meat quality and carcass characteristics. Two aspects have been considered in this thesis: the application of the transcriptome analysis in post mortem pig muscles as a possible method to evaluate meat quality parameters related to the pre mortem status of the animals, including health, nutrition, welfare, and with potential applications for product traceability (chapters 3 and 4); the study of candidate genes for obesity related traits in order to identify markers associated with fatness in pigs that could be applied to improve carcass quality (chapters 5, 6, and 7). Chapter three addresses the first issue from a methodological point of view. When we considered this issue, it was not obvious that post mortem skeletal muscle could be useful for transcriptomic analysis. Therefore we demonstrated that the quality of RNA extracted from skeletal muscle of pigs sampled at different post mortem intervals (20 minutes, 2 hours, 6 hours, and 24 hours) is good for downstream applications. Degradation occurred starting from 48 h post mortem even if at this time it is still possible to use some RNA products. In the fourth chapter, in order to demonstrate the potential use of RNA obtained up to 24 hours post mortem, we present the results of RNA analysis with the Affymetrix microarray platform that made it possible to assess the level of expression of more of 24000 mRNAs. We did not identify any significant differences between the different post mortem times suggesting that this technique could be applied to retrieve information coming from the transcriptome of skeletal muscle samples not collected just after slaughtering. This study represents the first contribution of this kind applied to pork. In the fifth chapter, we investigated as candidate for fat deposition the TBC1D1 [TBC1 (tre-2/USP6, BUB2, cdc16) gene. This gene is involved in mechanisms regulating energy homeostasis in skeletal muscle and is associated with predisposition to obesity in humans. By resequencing a fragment of the TBC1D1 gene we identified three synonymous mutations localized in exon 2 (g.40A>G, g.151C>T, and g.172T>C) and 2 polymorphisms localized in intron 2 (g.219G>A and g.252G>A). One of these polymorphisms (g.219G>A) was genotyped by high resolution melting (HRM) analysis and PCR-RFLP. Moreover, this gene sequence was mapped by radiation hybrid analysis on porcine chromosome 8. The association study was conducted in 756 performance tested pigs of Italian Large White and Italian Duroc breeds. Significant results were obtained for lean meat content, back fat thickness, visible intermuscular fat and ham weight. In chapter six, a second candidate gene (tribbles homolog 3, TRIB3) is analyzed in a study of association with carcass and meat quality traits. The TRIB3 gene is involved in energy metabolism of skeletal muscle and plays a role as suppressor of adipocyte differentiation. We identified two polymorphisms in the first coding exon of the porcine TRIB3 gene, one is a synonymous SNP (c.132T> C), a second is a missense mutation (c.146C> T, p.P49L). The two polymorphisms appear to be in complete linkage disequilibrium between and within breeds. The in silico analysis of the p.P49L substitution suggests that it might have a functional effect. The association study in about 650 pigs indicates that this marker is associated with back fat thickness in Italian Large White and Italian Duroc breeds in two different experimental designs. This polymorphisms is also associated with lactate content of muscle semimembranosus in Italian Large White pigs. Expression analysis indicated that this gene is transcribed in skeletal muscle and adipose tissue as well as in other tissues. In the seventh chapter, we reported the genotyping results for of 677 SNPs in extreme divergent groups of pigs chosen according to the extreme estimated breeding values for back fat thickness. SNPs were identified by resequencing, literature mining and in silico database mining. analysis, data reported in the literature of 60 candidates genes for obesity. Genotyping was carried out using the GoldenGate (Illumina) platform. Of the analyzed SNPs more that 300 were polymorphic in the genotyped population and had minor allele frequency (MAF) >0.05. Of these SNPs, 65 were associated (P<0.10) with back fat thickness. One of the most significant gene marker was the same TBC1D1 SNPs reported in chapter 5, confirming the role of this gene in fat deposition in pig. These results could be important to better define the pig as a model for human obesity other than for marker assisted selection to improve carcass characteristics.

Functional relationship of colorectal cancer-associated glycans with malignant phenotype and transcriptome changes

Glycosyltransferases ST6GAL1 and B4GALNT2 (and their cognate antigens Sia6LacNAc and Sda, respectively) are associated with colorectal cancer (CRC) but it is not fully clear their biological and clinical significance. We explored the clinical relevance of both glycosyltransferases by interrogating The Cancer Genome Atlas (TCGA) database while the phenotypic/transcriptomic effects of ST6GAL1/B4GALNT2 overexpression were studied in genetically modified CRC cell lines. Transcriptomic data from CRC patients in TCGA database suggested a moderate impact of ST6GAL1 on CRC progression, although it was not possible to define a clear role for this glycosyltransferase. Transcriptomic analysis of ST6GAL1-transduced cell lines revealed a much deeper effect of ST6GAL1 on gene expression in SW948 than in SW48. The overexpression of ST6GAL1 induced opposite effects on soft agar growth and wound healing in both cell lines. These results indicate that the impact of a cancer-associated glycosyltransferase change on phenotype/transcriptome can be extremely variable, depending on the molecular context of the tumor cell. On the contrary, transcriptomic analysis of B4GALNT2-modified cell lines together with TCGA database survey demonstrated a strong impact of B4GALNT2 on the transcriptional activity of CRC cells, in particular its association with a better prognosis. We suggest an anti-tumoral role of B4GALNT2 in CRC. We also investigated the glycan changes related to ST6GAL1/B4GALNT2 expression in a small cohort of tissues/plasma as well as the N-glycomic profile of CRC, normal and polyp tissues. We found an increase of ST6GAL1 activity in CRC and inflammatory bowel disease plasma samples comparing with plasma from healthy donors. A different Sda protein carrier pattern was observed between healthy donors and CRC plasma samples. β-arrestin 1 is a possible candidate as Sda carrier protein in plasma samples although future validation studies are needed. The alterations found in the N-glycan pattern highlight the importance of N-glycome as a molecular signature in cancer.

Evaluation of the transcriptomic response of an Alzheimer's disease murine model induced at different ages: searching for predictors

Alzheimer’s disease (AD) is a chronic, progressive neurodegenerative disease, characterized by the impairment of mnesic and cognitive functions, that represents the most frequent type of dementia in older people worldwide. Aging is the most important risk factor for the sporadic form of the pathology and it is associated to the progressive impairment of the proteostasis network. The endoplasmic reticulum (ER), the main cellular actor involved in proteostasis, appears significantly compromised in AD due to the accumulation of β-amyloid (Aβ) protein and phosphorylated-tau protein. Increasing proteins misfolding activates a specific cellular response known as Unfolded Protein response (UPR) which orchestrates the recovery of ER function. The aim of the present study was to investigate the role of UPR and aging process in a murine model of AD induced by intracerebroventricular (i.c.v.) injection of Aβ1-42 oligomers at 3 or 18 months. The oligomers injection in aged animals caused the increased of memory impairment, oxidative stress, and the depletion of glutathione reserve. Furthermore, the RNA-sequencing analysis was performed and the bioinformatic analysis showed the enrichment of several pathways involved in neurodegeneration and protein regulations. The following analysis highlighted the significant dysregulation of the three branches of the UPR, the protein kinase RNA-like ER kinase (PERK), inositol-requiring protein 1α (IRE1α) and activating transcription factor 6 (ATF-6). In turn, ER stress affected the PI3K/Akt/Gsk3β and MAPK/ERK pathways, highlighting Mapkapk5 as a potential marker of the neurodegenerative process, which regulation could lead to the definition of new pharmacological and neuroprotective strategies to counteract AD.