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.

Structural and functional analysis of centromeric chromatin

Animal neocentromeres are defined as ectopic centromeres that have formed in non-centromeric locations and avoid some of the features, like the DNA satellite sequence, that normally characterize canonical centromeres. Despite this, they are stable functional centromeres inherited through generations. The only existence of neocentromeres provide convincing evidence that centromere specification is determined by epigenetic rather than sequence-specific mechanisms. For all this reasons, we used them as simplified models to investigate the molecular mechanisms that underlay the formation and the maintenance of functional centromeres. We collected human cell lines carrying neocentromeres in different positions. To investigate the region involved in the process at the DNA sequence level we applied a recent technology that integrates Chromatin Immuno-Precipitation and DNA microarrays (ChIP-on-chip) using rabbit polyclonal antibodies directed against CENP-A or CENP-C human centromeric proteins. These DNA binding-proteins are required for kinetochore function and are exclusively targeted to functional centromeres. Thus, the immunoprecipitation of DNA bound by these proteins allows the isolation of centromeric sequences, including those of the neocentromeres. Neocentromeres arise even in protein-coding genes region. We further analyzed if the increased scaffold attachment sites and the corresponding tighter chromatin of the region involved in the neocentromerization process still were permissive or not to transcription of within encoded genes. Centromere repositioning is a phenomenon in which a neocentromere arisen without altering the gene order, followed by the inactivation of the canonical centromere, becomes fixed in population. It is a process of chromosome rearrangement fundamental in evolution, at the bases of speciation. The repeat-free region where the neocentromere initially forms, progressively acquires extended arrays of satellite tandem repeats that may contribute to its functional stability. In this view our attention focalized to the repositioned horse ECA11 centromere. ChIP-on-chip analysis was used to define the region involved and SNPs studies, mapping within the region involved into neocentromerization, were carried on. We have been able to describe the structural polymorphism of the chromosome 11 centromeric domain of Caballus population. That polymorphism was seen even between homologues chromosome of the same cells. That discovery was the first described ever. Genomic plasticity had a fundamental role in evolution. Centromeres are not static packaged region of genomes. The key question that fascinates biologists is to understand how that centromere plasticity could be combined to the stability and maintenance of centromeric function. Starting from the epigenetic point of view that underlies centromere formation, we decided to analyze the RNA content of centromeric chromatin. RNA, as well as secondary chemically modifications that involve both histones and DNA, represents a good candidate to guide somehow the centromere formation and maintenance. Many observations suggest that transcription of centromeric DNA or of other non-coding RNAs could affect centromere formation. To date has been no thorough investigation addressing the identity of the chromatin-associated RNAs (CARs) on a global scale. This prompted us to develop techniques to identify CARs in a genome-wide approach using high-throughput genomic platforms. The future goal of this study will be to focalize the attention on what strictly happens specifically inside centromere chromatin.

Molecular strategies for genetic diversity analysis and development of markers linked to resistance traits in apple

The goal of many plant scientists’ research is to explain natural phenotypic variation in term of simple changes in DNA sequence. DNA-based molecular markers are extensively used for the construction of genome-wide molecular maps and to perform genetic analysis for simple and complex traits. The PhD thesis was divided into two main research lines according to the different approaches adopted. The first research line is to analyze the genetic diversity in an Italian apple germplasm collection for the identification of markers tightly linked to targeted genes by an association genetic method. This made it possible to identify synomym and homonym accessions and triploids. The fruit red skin color trait has been used to test the reliability of the genetic approaches in this species. The second line is related to the development of molecular markers closely linked to the Rvi13 and Rvi5 scab resistance genes, previously mapped on apple’s chromosome 10 and 17 respectively by using the traditional linkage mapping method. Both region have been fine-mapped with various type of markers that could be used for marker-assisted selection in future breeding programs and to isolate the two resistance genes.

Holocentric chromosomes: convergent evolution, meiotic adaptations, and genomic analysis

In most eukaryotes, the kinetochore protein complex assembles at a single locus termed the centromere to attach chromosomes to spindle microtubules. Holocentric chromosomes have the unusual property of attaching to spindle microtubules along their entire length. Our mechanistic understanding of holocentric chromosome function is derived largely from studies in the nematode Caenorhabditis elegans, but holocentric chromosomes are found over a broad range of animal and plant species. In this review, we describe how holocentricity may be identified through cytological and molecular methods. By surveying the diversity of organisms with holocentric chromosomes, we estimate that the trait has arisen at least 13 independent times (four times in plants and at least nine times in animals). Holocentric chromosomes have inherent problems in meiosis because bivalents can attach to spindles in a random fashion. Interestingly, there are several solutions that have evolved to allow accurate meiotic segregation of holocentric chromosomes. Lastly, we describe how extensive genome sequencing and experiments in nonmodel organisms may allow holocentric chromosomes to shed light on general principles of chromosome segregation.

Bipolar disorder susceptibility region on chromosome 3q29 not confirmed in a case-control association study

We identified a bipolar disorder (BPD) susceptibility region on chromosome 3q29 in a genome-wide linkage scan (Bailer et al. 2002 (Biol Psychiatry 52: 40), NPL-score 4.09) and follow-up linkage analysis (Schosser et al. 2004 (J Psychiatr Res 38(3): 357), NPL-scores >3 with five markers). These findings were supported by further fine-mapping of this region (Schosser et al. 2007 (Eur Neuropsychopharmacol 17(6-7): 501)), finding NPL-scores >3.9 with SNPs (single nucleotide polymorphisms) spanning a region of 3.46 Mbp in BPD families. Since genetic association studies are more powerful than linkage studies for detecting susceptibility genes of small effect size, we aimed to replicate these findings in an independent case-control sample collected in London (UK) and Vienna (Austria).

Genetic analysis of syndactyly in German Holstein cattle

Congenital syndactyly with a variable number of affected feet was observed in eight black and white German Holstein calves. Analysis of the pedigree data revealed that all affected individuals could be traced back to a single founder. The pedigree was consistent with monogenic autosomal recessive inheritance and variable expressivity. Bovine syndactyly or "mulefoot" has been previously shown to map on the telomeric end of bovine chromosome 15 and we performed PCR genotyping of microsatellite markers spanning 27 cM of this chromosomal region to test the new cases for genetic linkage with the phenotype. The haplotype segregation confirmed the suggested inheritance pattern of the mulefoot mutation in this family and markers RM004, BM848 and BMS820 showed significant linkage to the phenotype. The results confirmed the chromosomal location of the mulefoot gene in this pedigree. Furthermore the study demonstrated that although marker testing has been available for nearly a decade the use of mulefoot carriers in cattle breeding remains uncontrolled. The presented family provides a resource for positional cloning of the causative mutation.

Linkage mapping of ovine microphthalmia to chromosome 23, the sheep orthologue of human chromosome 18

PURPOSE: To characterize the phenotype and map the locus responsible for autosomal recessive inherited ovine microphthalmia (OMO) in sheep. METHODS: Microphthalmia-affected lambs and their available relatives were collected in a field, and experimental matings were performed to obtain affected and normal lambs for detailed necropsy and histologic examinations. The matings resulted in 18 sheep families with 48 cases of microphthalmia. A comparative candidate gene approach was used to map the disease locus within the sheep genome. Initially, 27 loci responsible for the microphthalmia-anophthalmia phenotypes in humans or mice were selected to test for comparative linkage. Fifty flanking markers that were predicted from comparative genomic analysis to be closely linked to these genes were tested for linkage to the disease locus. After observation of statistical evidence for linkage, a confirmatory fine mapping strategy was applied by further genotyping of 43 microsatellites. RESULTS: The clinical and pathologic examinations showed slightly variable expressivity of isolated bilateral microphthalmia. The anterior eye chamber was small or absent, and a white mass admixed with cystic spaces extended from the papilla to the anterior eye chamber, while no recognizable vitreous body or lens was found within the affected eyes. Significant linkage to a single candidate region was identified at sheep chromosome 23. Fine mapping and haplotype analysis assigned the candidate region to a critical interval of 12.4 cM. This ovine chromosome segment encompasses an ancestral chromosomal breakpoint corresponding to two orthologue segments of human chromosomes 18, short and long arms. For the examined animals, we excluded the complete coding region and adjacent intronic regions of ovine TGIF1 to harbor disease-causing mutations. CONCLUSIONS: This is the first genetic localization for hereditary ovine isolated microphthalmia. It seems unlikely that a mutation in the TGIF1 gene is responsible for this disorder. The studied sheep represent a valuable large animal model for similar human ocular phenotypes.

A radiation hybrid map of sheep chromosome 23 based on ovine BAC-end sequences

More than 375,000 BAC-end sequences (BES) of the CHORI-243 ovine BAC library have been deposited in public databases. blastn searches with these BES against HSA18 revealed 1806 unique and significant hits. We used blastn-anchored BES for an in silico prediction of gene content and chromosome assignment of comparatively mapped ovine BAC clones. Ovine BES were selected at approximately 1.3-Mb intervals of HSA18 and incorporated into a human-sheep comparative map. An ovine 5000-rad whole-genome radiation hybrid panel (USUoRH5000) was typed with 70 markers, all of which mapped to OAR23. The resulting OAR23 RH map included 43 markers derived from BES with high and unique BLAST hits to the sequence of the orthologous HSA18, nine EST-derived markers, 16 microsatellite markers taken from the ovine linkage map and two bovine microsatellite markers. Six new microsatellite markers derived from the 43 mapped BES and the two bovine microsatellite markers were linkage-mapped using the International Mapping Flock (IMF). Thirteen additional microsatellite markers were derived from other ovine BES with high and unique BLAST hits to the sequence of the orthologous HSA18 and also positioned on the ovine linkage map but not incorporated into the OAR23 RH map. This resulted in 24 markers in common and in the same order between the RH and linkage maps. Eight of the BES-derived markers were mapped using fluorescent in situ hybridization (FISH), to thereby align the RH and cytogenetic maps. Comparison of the ovine chromosome 23 RH map with the HSA18 map identified and localized three major breakpoints between HSA18 and OAR23. The positions of these breakpoints were equivalent to those previously shown for syntenic BTA24 and HSA18. This study presents evidence for the usefulness of ovine BES when constructing a high-resolution comprehensive map for a single sheep chromosome. The comparative analysis confirms and refines knowledge about chromosomal conservation and rearrangements between sheep, cattle and human. The constructed RH map demonstrates the resolution and utility of the newly constructed ovine RH panel.

The locus for bovine dilated cardiomyopathy maps to chromosome 18

Bovine dilated cardiomyopathy (BDCMP) is a severe and terminal disease of the heart muscle observed in Holstein-Friesian cattle over the last 30 years. There is strong evidence for an autosomal recessive mode of inheritance for BDCMP. The objective of this study was to genetically map BDCMP, with the ultimate goal of identifying the causative mutation. A whole-genome scan using 199 microsatellite markers and one SNP revealed an assignment of BDCMP to BTA18. Fine-mapping on BTA18 refined the candidate region to the MSBDCMP06-BMS2785 interval. The interval containing the BDCMP locus was confirmed by multipoint linkage analysis using the software loki. The interval is about 6.7 Mb on the bovine genome sequence (Btau 3.1). The corresponding region of HSA19 is very gene-rich and contains roughly 200 genes. Although telomeric of the marker interval, TNNI3 is a possible positional and a functional candidate for BDCMP given its involvement in a human form of dilated cardiomyopathy. Sequence analysis of TNNI3 in cattle revealed no mutation in the coding sequence, but there was a G-to-A transition in intron 6 (AJ842179:c.378+315G>A). The analysis of this SNP using the study's BDCMP pedigree did not conclusively exclude TNNI3 as a candidate gene for BDCMP. Considering the high density of genes on the homologous region of HSA19, further refinement of the interval on BTA18 containing the BDCMP locus is needed.

A region on equine chromosome 13 is linked to recurrent airway obstruction in horses

REASONS FOR STUDY: Equine recurrent airway obstruction (RAO) is probably dependent on a complex interaction of genetic and environmental factors and shares many characteristic features with human asthma. Interleukin 4 receptor a chain (IL4RA) is a candidate gene because of its role in the development of human asthma, confirmation of this association is therefore required. METHODS: The equine BAC clone containing the IL4RA gene was localised to ECA13q13 by the FISH method. Microsatellite markers in this region were investigated for possible association and linkage with RAO in 2 large Warmblood halfsib families. Based on a history of clinical signs (coughing, nasal discharge, abnormal breathing and poor performance), horses were classified in a horse owner assessed respiratory signs index (HOARSI 1-4: from healthy, mild, moderate to severe signs). Four microsatellite markers (AHT133, LEX041, VHL47, ASB037) were analysed in the offspring of Sire 1 (48 unaffected HOARSI 1 vs. 59 affected HOARSI 2-4) and Sire 2 (35 HOARSI 1 vs. 50 HOARSI 2-4), age 07 years. RESULTS: For both sires haplotypes could be established in the order AHT133-LEXO47-VHL47-ASB37. The distances in this order were estimated to be 2.9, 0.9 and 2.3 centiMorgans, respectively. Haplotype association with mild to severe clinical signs of chronic lower airway disease (HOARSI 2-4) was significant in the offspring of Sire 1 (P = 0.026) but not significant for the offspring of Sire 2 (P = 0.32). Linkage analysis showed the ECA13q13 region containing IL4RA to be linked to equine chronic lower airway disease in one family (P<0.01), but not in the second family. CONCLUSIONS: This supports a genetic background for equine RAO and indicates that IL4RA is a candidate gene with possible locus heterogeneity for this disease. POTENTIAL RELEVANCE: Identification of major genes for RAO may provide a basis for breeding and individual prevention for this important disease.

Analysis and mapping of CACNB4, CHRNA1, KCNJ3, SCN2A and SPG4, physiological candidate genes for porcine congenital progressive ataxia and spastic paresis

The cause of porcine congenital progressive ataxia and spastic paresis (CPA) is unknown. This severe neuropathy manifests shortly after birth and is lethal. The disease is inherited as a single autosomal recessive allele, designated cpa. In a previous study, we demonstrated close linkage of cpa to microsatellite SW902 on porcine chromosome 3 (SSC3), which corresponds syntenically to human chromosome 2. This latter chromosome contains ion channel genes (Ca(2+), K(+) and Na(+)), a cholinergic receptor gene and the spastin (SPG4) gene, which cause human epilepsy and ataxia when mutated. We mapped porcine CACNB4, KCNJ3, SCN2A and CHRNA1 to SSC15 and SPG4 to SSC3 with the INRA-Minnesota porcine radiation hybrid panel (IMpRH) and we sequenced the entire open reading frames of CACNB4 and SPG4 without finding any differences between healthy and affected piglets. An anti-epileptic drug treatment with ethosuximide did not change the severity of the disease, and pigs with CPA did not exhibit the corticospinal tract axonal degeneration found in humans suffering from hereditary spastic paraplegia, which is associated with mutations in SPG4. For all these reasons, the hypothesis that CACNB4, CHRNA1, KCNJ3, SCN2A or SPG4 are identical with the CPA gene was rejected.

A de novo 1.1-1.6 Mb subtelomeric deletion of chromosome 20q13.33 in a patient with learning difficulties but without obvious dysmorphic features

We report on a de novo submicroscopic deletion of 20q13.33 identified by subtelomeric fluorescence in situ hybridization (FISH) in a 4-year-old girl with learning difficulties, hyperlaxity and strabismus, but without obvious dysmorphic features. Further investigations by array-based comparative genomic hybridization (array-CGH) and FISH analysis allowed us to delineate the smallest reported subterminal deletion of chromosome 20q, spanning a 1.1-1.6 Mb with a breakpoint localized between BAC RP5-887L7 and RP11-261N11. The genes CHRNA4 and KCNQ2 implicated in autosomal dominant epilepsy are included in the deletion interval. Subterminal 20q deletions as found in the present patient have, to our knowledge, only been reported in three patients. We review the clinical and behavioral phenotype of such "pure" subterminal 20q deletions.

Losses of chromosome arms 4q, 8p, 13q and gain of 8q are correlated with increasing chromosomal instability in hepatocellular carcinoma

OBJECTIVE: Chromosomal instability is a key feature in hepatocellular carcinoma (HCC). Array comparative genomic hybridization (aCGH) revealed recurring structural aberrations, whereas fluorescence in situ hybridization (FISH) indicated an increasing number of numerical aberrations in dedifferentiating HCC. Therefore, we examined whether there was a correlation between structural and numerical aberrations of chromosomal instability in HCC. METHODS AND RESULTS: 27 HCC (5 well, 10 moderately, 12 lower differentiated) already cytogenetically characterized by aCGH were analyzed. FISH analysis using probes for chromosomes 1, 3, 7, 8 and 17 revealed 1.46-4.24 signals/nucleus, which correlated with the histological grade (well vs. moderately,p < 0.0003; moderately vs. lower, p < 0.004). The number of chromosomes to each other was stable with exceptions only seen for chromosome 8. Loss of 4q and 13q, respectively, were correlated with the number of aberrations detected by aCGH (p < 0.001, p < 0.005; Mann-Whitney test). Loss of 4q and gain of 8q were correlated with an increasing number of numerical aberrations detected by FISH (p < 0.020, p < 0.031). Loss of 8p was correlated with the number of structural imbalances seen in aCGH (p < 0.048), but not with the number of numerical changes seen in FISH. CONCLUSION: We found that losses of 4q, 8p and 13q were closely correlated with an increasing number of aberrations detected by aCGH, whereas a loss of 4q and a gain of 8q were also observed in the context of polyploidization, the cytogenetic correlate of morphological dedifferentiation.

Genome sequence, comparative analysis, and population genetics of the domestic horse

We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

Phylogeny of horse chromosome 5q in the genus Equus and centromere repositioning

Horses, asses and zebras belong to the genus Equus and are the only extant species of the family Equidae in the order Perissodactyla. In a previous work we demonstrated that a key factor in the rapid karyotypic evolution of this genus was evolutionary centromere repositioning, that is, the shift of the centromeric function to a new position without alteration of the order of markers along the chromosome. In search of previously undiscovered evolutionarily new centromeres, we traced the phylogeny of horse chromosome 5, analyzing the order of BAC markers, derived from a horse genomic library, in 7 Equus species (E. caballus, E. hemionus onager, E. kiang, E. asinus, E. grevyi, E. burchelli and E. zebra hartmannae). This analysis showed that repositioned centromeres are present in E. asinus (domestic donkey, EAS) chromosome 16 and in E. burchelli (Burchell's zebra, EBU) chromosome 17, confirming that centromere repositioning is a strikingly frequent phenomenon in this genus. The observation that the neocentromeres in EAS16 and EBU17 are in the same chromosomal position suggests that they may derive from the same event and therefore, E. asinus and E. burchelli may be more closely related than previously proposed; alternatively, 2 centromere repositioning events, involving the same chromosomal region, may have occurred independently in different lineages, pointing to the possible existence of hot spots for neocentromere formation. Our comparative analysis also showed that, while E. caballus chromosome 5 seems to represent the ancestral configuration, centric fission followed by independent fusion events gave rise to 3 different submetacentric chromosomes in other Equus lineages.