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.

Mitochondrial DNA rearrangements during human aging: a study on liver, muscle and adipose tissue

Aging is characterized by a chronic, low-grade inflammatory state called “inflammaging”. Mitochondria are the main source of reactive oxygen species (ROS), which trigger the production of pro-inflammatory molecules. We are interested in studying the age-related modifications of the mitochondrial DNA (mtDNA), which can be affected by the lifelong exposure to ROS and are responsible of mitochondrial dysfunction. Moreover, increasing evidences show that telomere shortening, naturally occurring with aging, is involved in mtDNA damage processes and thus in the pathogenesis of age-related disorders. Thus the primary aim of this thesis was the analysis of mtDNA copy number, deletion level and integrity in different-age human biopsies from liver, vastus lateralis skeletal muscle of healthy subjects and patients with limited mobility of lower limbs (LMLL), as well as adipose tissue. The telomere length and the expression of nuclear genes related to mitobiogenesis, fusion and fission, mitophagy, mitochondrial protein quality control system, hypoxia, production and protection from ROS were also evaluated. In liver the decrease in mtDNA integrity with age is accompanied with an increase in mtDNA copy number, suggesting the existence of a “compensatory mechanism” able to maintain the functionality of this organ. Different is the case of vastus lateralis muscle, where any “compensatory pathway” is activated and mtDNA integrity and copy number decrease with age, both in healthy subjects and in patients. Interestingly, mtDNA rearrangements do not incur in adipose tissue with advancing age. Finally, in all tissues a marked gender difference appears, suggesting that aging and also gender diversely affect mtDNA rearrangements and telomere length in the three human tissues considered, likely depending on their different metabolic needs and inflammatory status.

Azione antiossidante del sulforafane: analisi in cellule cardiache in coltura ed in un modello animale di esercizio fisico

Oxidative stress is considered to be of major relevance for a variety of pathological processes. Thus, it is valuable to identify compounds, which might act as antioxidants, i.e. compounds that antagonize the deleterious action of reactive oxygen species (ROS) on biomolecules. The mode of action of these compounds could be either to scavenge ROS directly or to trigger protective mechanisms inside the cell, thereby resulting in improved defense against ROS. Sulforaphane (SF) (1-isothiocyanato-(4R)-(methylsulfinyl)butane) is a naturally occurring cancer chemopreventive agent found as a precursor glucosinolate in Cruciferous vegetables like broccoli. Although SF is not a direct-acting antioxidant, there is substantial evidence that SF acts indirectly to increase the antioxidant capacity of animal cells and their abilities to cope with oxidative stress. Induction of phase 2 enzymes is one means by which SF enhances the cellular antioxidant capacity. Enzymes induced by SF include Glutathione S-transferases (GST) and NAD[P]H:quinone oxidoreductase (NQO1) which can function as protectors against oxidative stress. To protect themselves from oxidative stress, cells are equipped with reducing buffer systems including the GSH and thioredoxin (Trx) reductase. GSH is an important tripeptide thiol which in addition to being the substrate for GSTs maintains the cellular oxidation– reduction balance and protects cells against free radical species. Aim of the first part of this thesis was to investigate the ability of SF to induce the expression and the activity of different phase 2 and antioxidant enzymes (such as GST, GR, GPx, NQO1, TR, SOD, CAT) in an in vitro model of rat cardiomyocytes, and also to define if SF treatment supprts cells in counteracting oxidative stress induced by H2O2 It is well known that acute exhaustive exercise causes significant reactive oxygen species generation that results in oxidative stress, which can induce negative effects on health and well being. In fact, increased oxidative stress and biomarkers (e.g., protein carbonyls, MDA, and 8- hydroxyguanosine) as well as muscle damage biomarkers (e.g. plasmatic Creatine cinase and Lactate dehydrogenase) have been observed after supramaximal sprint exercises, exhaustive longdistance cycling or running as well as resistance-type exercises, both in trained and untrained humans. Markers of oxidative stress also increase in rodents following exhaustive exercise. Moreover, antioxidant enzyme activities and expressions of antioxidant enzymes are known to increase in response to exhaustive exercise in both animal and human tissues. Aim of this project was to evaluate the effect of SF supplementation in counteracting oxidative stress induced by physical activity through its ability to induce phase 2, and antioxidant enzymes in rat muscle. The results show that SF is a nutraceutical compound able to induce the activity of different phase 2 and antioxidant enzymes in both cardiac muscle and skeletal muscle. Thanks to its actions SF is becoming a promising molecule able to prevent cardiovascular damages induced by oxidative stress and muscle damages induced by acute exhaustive exercise.

Characterization of L-type Calcium Channel Binding-Site of a new class of Calcium modulators by a Multidisciplinary approach

Many potential diltiazem related L-VDCC blockers were developed using a multidisciplinary approach. This current study was to investigate and compare diltiazem with to the newly developed compounds by mouse Langendorff-perfused heart, Ca2+-transient and on recombinant L-VDCC. Twenty particular compounds were selected by the ligand-based virtual screening procedure (LBVS). From these compounds, five of them (5b, M2, M7, M8 and P1) showed a potent and selective inotropic activity on guinea-pig left atria driven 1 Hz. Further assays displayed an interesting negative inotropic effect of M2, M8, P1 and M7 on guinea pig isolated left papillary muscle driven at 1 Hz, a relevant vasorelaxant activity of 5b, M2, M7, M8 and P1 on K+-depolarized guinea-pig ileum longitudinal smooth muscle and a significant inhibition of contraction of 5b, M2, M8 and P1 on carbachol stimulated ileum longitudinal smooth muscle. Wild-type human heart and rabbit lung α1 subunits were expressed (combined with the regulatory α2δ and β3 subunits) in Xenopus Leavis oocytes using a two-electrode voltage clamp technique. Diltiazem is a benzothiazepine Ca2+ channel blocker used clinically for its antihypertensive and antiarrhythmic effects. Previous radioligand binding assays revealed a complex interaction with the benzothiazepine binding site for M2, M7 and M8. (Carosati E. et al. J. Med Chem. 2006, 49; 5206). In agreement with this findings, the relative order of increased rates of contraction and relaxation at lower concentrations s(≤10-6M) in unpaced hearts was M7>M2>M8>P1. Similar increases in Ca2+ transient were observed in cardiomyocytes. Diltiazem showed negative inotropic effects whereas 5b had no significant effect. Diltiazem blocks Ca2+current in a use-dependent manner and facilitates the channel by accelerating the inactivation and decelerating the recovery from inactivation. In contrast to diltiazem, the new analogs had no pronounced use-dependence. Application of 100 μM M8, M2 showed ~ 10% tonic block; in addition, M8, M2 and P1 shifted the steady state inactivation in hyperpolarized direction and the current inactivation time was significantly decreased compared with control (219.6 ± 11.5 ms, 226 ± 14.5 vs. 269 ± 12.9 vs. 199.28 ± 8.19 ms). Contrary to diltiazem, the recovery from the block by M8 and M2 was comparable to control. Only P1 showed a significantly decrease of the time for the recovery from inactivation. All of the compounds displayed the same sensitivity on the Ca2+ channel rabbit lung α1 except P1. Taken together, these findings suggest that M8, M2 and P1 might directly decrease the binding affinity or allow rapid dissociation from the benzothiazepine binding site.

Wide-scale population genomics of Atlantic bluefin tuna (Thunnus thynnus) inferred by novel high-throughput technology

My PhD project was focused on Atlantic bluefin tuna, Thunnus thynnus, a fishery resource overexploited in the last decades. For a better management of stocks, it was necessary to improve scientific knowledge of this species and to develop novel tools to avoid collapse of this important commercial resource. To do this, we used new high throughput sequencing technologies, as Next Generation Sequencing (NGS), and markers linked to expressed genes, as SNPs (Single Nucleotide Polymorphisms). In this work we applied a combined approach: transcriptomic resources were used to build cDNA libreries from mRNA isolated by muscle, and genomic resources allowed to create a reference backbone for this species lacking of reference genome. All cDNA reads, obtained from mRNA, were mapped against this genome and, employing several bioinformatics tools and different restricted parameters, we achieved a set of contigs to detect SNPs. Once a final panel of 384 SNPs was developed, following the selection criteria, it was genotyped in 960 individuals of Atlantic bluefin tuna, including all size/age classes, from larvae to adults, collected from the entire range of the species. The analysis of obtained data was aimed to evaluate the genetic diversity and the population structure of Thunnus thynnus. We detect a low but significant signal of genetic differentiation among spawning samples, that can suggest the presence of three genetically separate reproduction areas. The adult samples resulted instead genetically undifferentiated between them and from the spawning populations, indicating a presence of panmictic population of adult bluefin tuna in the Mediterranean Sea, without different meta populations.

Subject-specific musculoskeletal models of the lower limbs for the prediction of skeletal loads during motion

The determination of skeletal loading conditions in vivo and their relationship to the health of bone tissues, remain an open question. Computational modeling of the musculoskeletal system is the only practicable method providing a valuable approach to muscle and joint loading analyses, although crucial shortcomings limit the translation process of computational methods into the orthopedic and neurological practice. A growing attention focused on subject-specific modeling, particularly when pathological musculoskeletal conditions need to be studied. Nevertheless, subject-specific data cannot be always collected in the research and clinical practice, and there is a lack of efficient methods and frameworks for building models and incorporating them in simulations of motion. The overall aim of the present PhD thesis was to introduce improvements to the state-of-the-art musculoskeletal modeling for the prediction of physiological muscle and joint loads during motion. A threefold goal was articulated as follows: (i) develop state-of-the art subject-specific models and analyze skeletal load predictions; (ii) analyze the sensitivity of model predictions to relevant musculotendon model parameters and kinematic uncertainties; (iii) design an efficient software framework simplifying the effort-intensive phases of subject-specific modeling pre-processing. The first goal underlined the relevance of subject-specific musculoskeletal modeling to determine physiological skeletal loads during gait, corroborating the choice of full subject-specific modeling for the analyses of pathological conditions. The second goal characterized the sensitivity of skeletal load predictions to major musculotendon parameters and kinematic uncertainties, and robust probabilistic methods were applied for methodological and clinical purposes. The last goal created an efficient software framework for subject-specific modeling and simulation, which is practical, user friendly and effort effective. Future research development aims at the implementation of more accurate models describing lower-limb joint mechanics and musculotendon paths, and the assessment of an overall scenario of the crucial model parameters affecting the skeletal load predictions through probabilistic modeling.

Transcriptional dynamics of the human DMD locus

The human DMD locus encodes dystrophin protein. Absence or reduced levels of dystrophin (DMD or BMD phenotype, respectively) lead to progressive muscle wasting. Little is known about the complex coordination of dystrophin expression and its transcriptional regulation is a field of intense interest. In this work we found that DMD locus harbours multiple long non coding RNAs which orchestrate and control transcription of muscle dystrophin mRNA isoforms. These lncRNAs are tissue-specific and highly expressed during myogenesis, suggesting a possible role in tissue-specific expression of DMD gene isoforms. Their forced ectopic expression in human muscle and neuronal cells leads to a specific and negative regulation of endogenous dystrophin full lenght isoforms. An intriguing aspect regarding the transcription of the DMD locus is the gene size (2.4Mb). The mechanism that ensures the complete synthesis of the primary transcript and the coordinated splicing of 79 exons is still completely unknown. By ChIP-on-chip analyses, we discovered novel regions never been involved before in the transcription regulation of the DMD locus. Specifically, we observed enrichments for Pol II, P-Ser2, P-Ser5, Ac-H3 and 2Me-H3K4 in an intronic region of 3Kb (approximately 21Kb) downstream of the end of DMD exon 52 and in a region of 4Kb spanning the DMD exon 62. Interestingly, this latter region and the TSS of Dp71 are strongly marked by 3Me-H3K36, an histone modification associated with the regulation of splicing process. Furthermore, we also observed strong presence of open chromatin marks (Ac-H3 and 2Me-H3K4) around intron 34 and the exon 45 without presence of RNA pol II. We speculate that these two regions may exert an enhancer-like function on Dp427m promoter, although further investigations are necessary. Finally, we investigated the nuclear-cytoplasmic compartmentalization of the muscular dystrophin mRNA and, specifically, we verified whether the exon skipping therapy could influence its cellular distribution.