Molecular characterization and methylation study of matrix gla protein in articular cartilage from pig with osteochondrosis

Comparative label-free proteomic analysis of equine osteochondrotic chondrocytes

Osteochondrosis is a developmental orthopedic disease affecting growing cartilage in young horses. In this study we compared the proteomes of equine chondrocytes obtained from healthy and osteochondrotic cartilage using a label-free mass spectrometry approach. Quantitative changes of some proteins selected for their involvement in different functional pathways highlighted by the bioinformatics analysis, were validated by western blotting, while biochemical alterations of extracellular matrix were confirmed via Raman spectroscopy analysis. In total 1637 proteins were identified, of which 59 were differentially abundant. Overall, the results highlighted differentially represented proteins involved in metabolic and functional pathways that may be related to the failure of the endochondral ossification process occurring in osteochondrosis. In particular, we identified proteins involved in extracellular matrix degradation and organization, vitamin metabolism, osteoblast differentiation, apoptosis, protein folding and localization, signalling and gene expression modulation and lysosomal activities. These results provide valuable new insights to elucidate the underlying molecular mechanisms associated with the development and progression of osteochondrosis. SIGNIFICANCE: Osteochondrosis is a common articular disorder in young horses mainly due to defects in endochondral ossification. The pathogenesis of osteochondrosis is still poorly understood and only a limited number of proteomic studies have been conducted. This study provides a comprehensive characterization of proteomic alterations occurring in equine osteochondrotic chondrocytes, the only resident cell type that modulates differentiation and maturation of articular cartilage. The results evidenced alterations in abundance of proteins involved in functional and metabolic pathways and in extracellular matrix remodelling. These findings could help clarify some molecular aspects of osteochondrosis and open new fields of research for elucidating the pathogenesis of this disease.

Association of matrix Gla protein polymorphism and knee osteoarthritis in a chinese population

Several studies have explored the association between matrix Gla protein (MGP) gene polymorphism and knee osteoarthritis (OA) risk; however, they obtained conflicting findings. The present study aims to explore the association of MGP gene polymorphism and OA risk in a Chinese Han population. A total of 256 patients with radiographic knee OA and 327 control subjects were recruited in this case-control study. The genotypes of MGP gene rs1800802 polymorphism was determined by standard PCR and restriction fragment length polymorphism (PCR-RLFP). In this case-control study, we observed that MGP gene rs1800802 polymorphism increased the risk of knee OA. Subgroup analyses also found that rs1800802 polymorphism was related to the elevated risk for knee OA among the female, smoker, drinker, and body mass index (BMI) ≥25 kg/m2 groups. In conclusion, this study shows that MGP gene rs1800802 polymorphism is associated with increased risk for knee OA in Chinese Han population and the rs1800802 polymorphism may be a diagnostic marker of radiographic knee OA.

Balancing selection at a premature stop mutation in the myostatin gene underlies a recessive leg weakness syndrome in pigs

Balancing selection provides a plausible explanation for the maintenance of deleterious alleles at moderate frequency in livestock, including lethal recessives exhibiting heterozygous advantage in carriers. In the current study, a leg weakness syndrome causing mortality of piglets in a commercial line showed monogenic recessive inheritance, and a region on chromosome 15 associated with the syndrome was identified by homozygosity mapping. Whole genome resequencing of cases and controls identified a mutation causing a premature stop codon within exon 3 of the porcine Myostatin (MSTN) gene, similar to those causing a double-muscling phenotype observed in several mammalian species. The MSTN mutation was in Hardy-Weinberg equilibrium in the population at birth, but significantly distorted amongst animals still in the herd at 110 kg, due to an absence of homozygous mutant genotypes. In heterozygous form, the MSTN mutation was associated with a major increase in muscle depth and decrease in fat depth, suggesting that the deleterious allele was maintained at moderate frequency due to heterozygous advantage (allele frequency, q = 0.22). Knockout of the porcine MSTN by gene editing has previously been linked to problems of low piglet survival and lameness. This MSTN mutation is an example of putative balancing selection in livestock, providing a plausible explanation for the lack of disrupting MSTN mutations in pigs despite many generations of selection for lean growth.

Transcriptome signature of liver tissue with divergent mutton odour and flavour using RNA deep sequencing

Mutton consumption is less popular in many Asian countries including Indonesia, whose consumers often complain about the unpleasant flavour and odour of the meat. The main causes of mutton odour are the two compounds of branched chain fatty acid (BCFA): methylnonanoic (MNA), phenol, 3-methyl (MP), 4-methylnonanoic (MNA) and 4-ethyloctanoic (EOA) present in all the adipose tissue; and the 3-methylindole (MI) or skatole and indole, which are originated from pastoral diets. It is crucial to understand the genetic mechanism of mutton odour and flavour (MOF) to select sheep for lower BCFA and indole thus reduce the unpleasant flavour of meat. The aim of the present study was to investigate transcriptome profiling in liver tissue with divergent MOF using RNA deep sequencing. Liver tissues from higher (n = 3) and lower (n = 3) MOF sheep were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample ranged from 21.37 to 25.37 million. Approximately 103 genes were differentially expressed (DEGs) with significance level of p-adjusted value <0.05. Among them, 60 genes were up-regulated, and 43 were down-regulated (p < 0.01, FC > 1.5) in higher MOF group. Differentially regulated genes in high MOF liver samples were enriched in biological processes such as cellular response to chemical stimulus and endogenous stimulus; cellular components such as such as basement membrane and extracellular matrix; and molecular functions such as haeme binding and oxidoreductase activity. Among the DEGs, metabolic phase I related genes belonging to the cytochrome P450 CYP2A6 were dominantly expressed. Additionally, phase II conjugation genes including UDP glucuronosyltransferases UGT2B18, sulfotransferase SULT1C1, and glutathione S-transferase GSTM1 were identified. The dominant candidate genes for SOF could be cytochrome P450, sodium-channel protein, transmembrane protein, glutathione transferase, UDP glucuronosyltransferases and sulfotransferase. Pathway analysis identified steroid hormone biosynthesis and chemical carcinogenesis by cytochrome P450 pathways which may play important roles in MOF-related molecules metabolism. This work highlighted potential genes and gene-networks that may affect meat off flavour and odour in sheep.

Polymorphisms in ten candidate genes are associated with conformational and locomotive traits in Spanish Purebred horses

The Spanish Purebred horses, also known as Andalusian horses, compete to the highest standards in international dressage events. Gait and conformation could be used as early selection criteria to detect young horses with promising dressage ability. Although the genetic background of equine size variation has been recently uncovered, the genetic basis of horse conformational and locomotive traits is not known, hampered by the complex genetic architecture underlying quantitative traits and the lack of phenotypic data. The aim of this study was to validate the loci associated with size in 144 Spanish Purebred horses, and to seek novel associations between loci previously associated with the development of osteochondrosis (OC) lesions and 20 conformational and locomotive traits. Ten loci were associated with different conformational and locomotive traits (LCORL/NCAPG, HMGA2, USP31, MECR, COL24A1, MGP, FAM184B, PTH1R, KLF3 and SGK1), and the LCORL/NCAPG association with size in the Spanish Purebred horse was validated. Except for HMGA2, all polymorphisms seem to influence both the prevalence of OC lesions and morphological characters, supporting the link between conformation and OC. Also, the implication of most genes in either immune and inflammatory responses and cellular growth, or ossification processes, reinforces the role that these mechanisms have in the aetiology of OC, as well as their reflection on the general conformation of the individual. These polymorphisms could be used in marker-assisted selection (MAS) programmes to improve desirable conformational traits, but taking into account their possible detrimental effect on OC prevalence.

An Update on the Pathogenesis of Osteochondrosis

Osteochondrosis is defined as a focal disturbance in endochondral ossification. The cartilage superficial to an osteochondrosis lesion can fracture, giving rise to fragments in joints known as osteochondrosis dissecans (OCD). In pigs and horses, it has been confirmed that the disturbance in ossification is the result of failure of the blood supply to epiphyseal growth cartilage and associated ischemic chondronecrosis. The earliest lesion following vascular failure is an area of ischemic chondronecrosis at an intermediate depth of the growth cartilage (osteochondrosis latens) that is detectable ex vivo, indirectly using contrast-enhanced micro- and conventional computed tomography (CT) or directly using adiabatic T1ρ magnetic resonance imaging. More chronic lesions of ischemic chondronecrosis within the ossification front (osteochondrosis manifesta) are detectable by the same techniques and have also been followed longitudinally in pigs using plain CT. The results confirm that lesions sometimes undergo spontaneous resolution, and in combination, CT and histology observations indicate that this occurs by filling of radiolucent defects with bone from separate centers of endochondral ossification that form superficial to lesions and by phagocytosis and intramembranous ossification of granulation tissue that forms deep to lesions. Research is currently aimed at discovering the cause of the vascular failure in osteochondrosis, and studies of spontaneous lesions suggest that failure is associated with the process of incorporating blood vessels into the advancing ossification front during growth. Experimental studies also show that bacteremia can lead to vascular occlusion. Future challenges are to differentiate between causes of vascular failure and to discover the nature of the heritable predisposition for osteochondrosis.

Genome-wide association analysis of osteochondrosis of the tibiotarsal joint in Norwegian Standardbred trotters

Osteochondrosis (OC), a disturbance in the process of endochondral ossification, is by far the most important equine developmental orthopaedic disease and is also common in other domestic animals and humans. The purpose of this study was to identify quantitative trait loci (QTL) associated with osteochondrosis dissecans (OCD) at the intermediate ridge of the distal tibia in Norwegian Standardbred (SB) using the Illumina Equine SNP50 BeadChip whole-genome single-nucleotide polymorphism (SNP) assay. Radiographic data and blood samples were obtained from 464 SB yearlings. Based on the radiographic examination, 162 horses were selected for genotyping; 80 of these were cases with an OCD at the intermediate ridge of the distal tibia, and 82 were controls without any developmental lesions in the joints examined. Genotyped horses descended from 22 sires, and the number of horses in each half-sib group ranged from 3 to 14. The population structure necessitated statistical correction for stratification. When conducting a case-control genome-wide association study (GWAS), mixed-model analyses displayed regions on chromosomes (Equus callabus chromosome - ECA) 5, 10, 27 and 28 that showed moderate evidence of association (P ≤ 5 × 10(-5); this P-value is uncorrected i.e. not adjusted for multiple comparisons) with OCD in the tibiotarsal joint. Two SNPs on ECA10 represent the most significant hits (uncorrected P=1.19 × 10(-5) in the mixed-model). In the basic association (chi-square) test, these SNPs achieved statistical significance with the Bonferroni correction (P=0.038) and were close in the permuted logistic regression test (P=0.054). Putative QTL on ECA 5, 10, 27 and 28 represent interesting areas for future research, validation studies and fine mapping of candidate regions. Results presented here represent the first GWAS of OC in horses using the recently released Illumina Equine SNP50 BeadChip.

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0068-2) contains supplementary material, which is available to authorized users.

Fibrocalcific aortic valve disease: opportunity to understand disease mechanisms using mouse models

Studies in vitro and in vivo continue to identify complex-regulated mechanisms leading to overt fibrocalcific aortic valve disease (FCAVD). Assessment of the functional impact of those processes requires careful studies of models of FCAVD in vivo. Although the genetic basis for FCAVD is unknown for most patients with FCAVD, several disease-associated genes have been identified in humans and mice. Some gene products which regulate valve development in utero also protect against fibrocalcific disease during postnatal aging. Valve calcification can occur via processes that resemble bone formation. But valve calcification can also occur by nonosteogenic mechanisms, such as formation of calcific apoptotic nodules. Anticalcific interventions might preferentially target either osteogenic or nonosteogenic calcification. Although FCAVD and atherosclerosis share several risk factors and mechanisms, there are fundamental differences between arteries and the aortic valve, with respect to disease mechanisms and responses to therapeutic interventions. Both innate and acquired immunity are likely to contribute to FCAVD. Angiogenesis is a feature of inflammation, but may also contribute independently to progression of FCAVD, possibly by actions of pericytes that are associated with new blood vessels. Several therapeutic interventions seem to be effective in attenuating the development of FCAVD in mice. Therapies which are effective early in the course of FCAVD, however, are not necessarily effective in established disease.

Gene expression profiling of articular cartilage reveals functional pathways and networks of candidate genes for osteochondrosis in pigs

Osteochondrosis (OC) is a joint disorder that frequently causes leg weakness in growing pigs, resulting in welfare problems and economic losses. We aimed to detect molecular pathways relevant to the emergence of the disease and to identify candidate genes for the liability to the disorder. Therefore, we compared microarray-based expression patterns of articular cartilage with (n=11) and without (n=11) histologically diagnosed OC lesions obtained from discordant sib-pairs. A total of 1,564 genes were found with different transcript abundance [differentially expressed (DE) genes] at q≤0.05. To further identify candidate genes, we integrated data from quantitative trait loci (QTL) and genome-wide association (GWA) studies with the expression analysis. We detected 317 DE genes within the QTL confidence intervals, of which 26 DE genes also overlapped GWA regions. Ingenuity Pathway Analysis suggests a pathogenic role of immune response, angiogenesis, and synthesis of extracellular matrix pathways for OC. These processes could facilitate the emergence of defects. But they may also promote the degradation of articular cartilage and the worsening of the disease. A functional network was derived that comprised genes with functional and positional clues of their role in bone and cartilage metabolisms and development, including extracellular matrix genes (e.g., LOX, OGN, and ASPN), angiogenesis genes (e.g., ANGPTL4 and PDGFA), and immune response genes (e.g., ICAM1, AZGP1, C1QB, C1QC, PDE4B, and CDA). The study identified molecular processes linked to OC and several genes with positional, genetic-statistical, and functional evidence for their role in the emergence of articular cartilage lesions and the liability to OC.

Association and expression analysis of porcine HNF1A gene related to meat and carcass quality traits

The aim of this study was to investigate the association and expression of HNF1A gene as a candidate gene for meat and carcass quality traits in pigs. Statistical analysis revealed that the g.8260 A>G polymorphism significantly associated with pH 24(H), meat percentage and muscle area in the F₂ Duroc × Pietrain (DuPi, n=313) and with pH 24(L), fat area and backfat thickness in the Pietrain (Pi, n=110) population. HNF1A mRNA and protein expressions were higher (p<0.05) in animals with the low post-mortem muscle pH 24(L). The promoter methylation profiling suggested that methylation was not involved on HNF1A expression regulation (p>0.05) in animal with divergent muscle pH. In conclusion, polymorphism in porcine HNF1A gene could be used as a candidate marker to improve the meat and carcass quality traits, with the consideration of breed-specific effect.

KRT8, FAF1 and PTH1R gene polymorphisms are associated with leg weakness traits in pigs

The liability to lesions of dysfunctions of bone and joints in pigs, summarized as leg weakness and mostly expressed as osteochondrosis, is an animal welfare and economic issue in pig production. The objective of this study was to identify polymorphisms in the functional and positional candidate genes keratin 8 (KRT8), Fas-associated factor 1 (FAF1) and parathyroid hormone type I receptor (PTH1R) and to evaluate their association with leg weakness traits. Therefore, osteochondrosis lesions were scored in animals of a Duroc × Pietrain F2 population (DuPi; n = 310) and commercial herds of the breed Large White (n = 298). In addition, bone mineralization traits were observed in DuPi population. SNPs were identified in genes KRT8 (g.8,039G > A), FAF1 (g.380,914T > C) and PTH1R (c.1,672C > T). KRT8 showed significant association with bone mineral density and content (P ≤ 0.05). FAF1 was association with OC lesions score of all joints inspected (P ≤ 0.05). PTH1R showed significant dominance effects on OC lesion scores of the distal femur articular cartilage (P = 0.01) and epiphysis of the distal ulna (P = 0.05) as well as sums of scores of all joints (OCsum, P = 0.04) and assignment to groups of either severely or gently affected animals (OCcat, P = 0.01). This study reveals clear genetic-statistical evidence for a link of KRT8, FAF1 and PTH1R with some of leg weakness related traits in pigs.

Evaluation of suitable reference genes for gene expression studies in porcine alveolar macrophages in response to LPS and LTA

Background

To obtain reliable quantitative real-time PCR data, normalization relative to stable housekeeping genes (HKGs) is required. However, in practice, expression levels of 'typical' housekeeping genes have been found to vary between tissues and under different experimental conditions. To date, validation studies of reference genes in pigs are relatively rare and have never been performed in porcine alveolar macrophages (AMs). In this study, expression stability of putative housekeeping genes were identified in the porcine AMs in response to the stimulation with two pathogen-associated molecular patterns (PAMPs) lipopolysaccharide (LPS) and lipoteichoic acid (LTA). Three different algorithms (geNorm, Normfinder and BestKeeper) were applied to assess the stability of HKGs.

Results

The mRNA expression stability of nine commonly used reference genes (B2M, BLM, GAPDH, HPRT1, PPIA, RPL4, SDHA, TBP and YWHAZ) was determined by qRT-PCR in AMs that were stimulated by LPS and LTA in vitro. mRNA expression levels of all genes were found to be affected by the type of stimulation and duration of the stimulation (P < 0.0001). geNorm software revealed that SDHA, B2M and RPL4 showed a high expression stability in the irrespective to the stimulation group, while SDHA, YWHAZ and RPL4 showed high stability in non-stimulated control group. In all cases, GAPDH showed the least stability in geNorm. NormFinder revealed that SDHA was the most stable gene in all the groups. Moreover, geNorm software suggested that the geometric mean of the three most stable genes would be the suitable combination for accurate normalization of gene expression study.

Conclusions

There was discrepancy in the ranking order of reference genes obtained by different analysing algorithms. In conclusion, the geometric mean of the SDHA, YWHAZ and RPL4 seemed to be the most appropriate combination of HKGs for accurate normalization of gene expression data in porcine AMs without knowing the type of bacterial pathogenic status of the animals.

Age-related changes in relative expression stability of commonly used housekeeping genes in selected porcine tissues

Background

Gene expression analysis using real-time RT-PCR (qRT-PCR) is increasingly important in biological research due to the high-throughput and accuracy of qRT-PCR. For accurate and reliable gene expression analysis, normalization of gene expression data against housekeeping genes or internal control genes is required. The stability of reference genes has a tremendous effect on the results of relative quantification of gene expression by qRT-PCR. The expression stability of reference genes could vary according to tissues, age of individuals and experimental conditions. In the pig however, very little information is available on the expression stability of reference genes. The aim of this research was therefore to develop a new set of reference genes which can be used for normalization of mRNA expression data of genes expressed in varieties of porcine tissues at different ages.

Results

The mRNA expression stability of nine commonly used reference genes (B2M, BLM, GAPDH, HPRT1, PPIA, RPL4, SDHA, TBP and YWHAZ) was determined in varieties of tissues collected from newborn, young and adult pigs. geNorm, NormFinder and BestKeeper software were used to rank the genes according to their stability. geNorm software revealed that RPL4, PPIA and YWHAZ showed high stability in newborn and adult pigs, while B2M, YWHAZ and SDHA showed high stability in young pigs. In all cases, GAPDH showed the least stability in geNorm. NormFinder revealed that TBP was the most stable gene in newborn and young pigs, while PPIA was most stable in adult pigs. Moreover, geNorm software suggested that the geometric mean of three most stable gene would be the suitable combination for accurate normalization of gene expression study.

Conclusions

Although, there was discrepancy in the ranking order of reference genes obtained by different analysing software methods, the geometric mean of the RPL4, PPIA and YWHAZ seems to be the most appropriate combination of housekeeping genes for accurate normalization of gene expression data in different porcine tissues at different ages.

Association and expression study of MMP3, TGFβ1 and COL10A1 as candidate genes for leg weakness-related traits in pigs

The present study was aimed to determine the association between metalloproteinase 3 (MMP3), transforming growth factor beta 1 (TGFβ1) and collagen type X alpha I (COL10A1) gene polymorphisms with traits related to leg weakness in pigs. Three hundred Duroc × Pietrain cross breds (DuPi) and 299 pigs of a commercial population (CP) were used for the experiment. DuPi animals were examined for 10 different traits describing leg and feet structure, osteochondrosis (OC) scores and bone density status. Data of OC score at condylus medialis humeri, condylus medialis femoris and distal epiphysis ulna regions of CP were used for association analysis. Significant association (P < 0.05) was found for MMP3 SNP (g.158 C>T) with OC at head of femur and bone mineral density in the DuPi population. Association (P < 0.05) was found between SNP of TGFβ1 (g.180 G>A) with rear leg score and the principle component denoting both OC and feet and leg scores in the DuPi population. No association was found between COL10A1 (g.72 C>T) and leg weakness related traits. The associations of SNPs with OC traits could not be confirmed in the commercial population. Expression analysis of the three candidate genes was performed to compare between healthy and OC. TGFβ1 was found to be highly expressed (P < 0.05) in the OC compared to healthy cartilages, but no significant different expressions were observed for MMP3 and COL10A1 genes. The present finding suggested that TGFβ1 and MMP3 genes variants have an effect on some of the leg weakness related traits.

Quantitative trait loci analysis for leg weakness-related traits in a Duroc × Pietrain crossbred population

Background

Leg weakness issues are a great concern for the pig breeding industry, especially with regard to animal welfare. Traits associated with leg weakness are partly influenced by the genetic background of the animals but the genetic basis of these traits is not yet fully understood. The aim of this study was to identify quantitative trait loci (QTL) affecting leg weakness in pigs.

Methods

Three hundred and ten F₂ pigs from a Duroc × Pietrain resource population were genotyped using 82 genetic markers. Front and rear legs and feet scores were based on the standard scoring system. Osteochondrosis lesions were examined histologically at the head and the condylus medialis of the left femur and humerus. Bone mineral density, bone mineral content and bone mineral area were measured in the whole ulna and radius bones using dual energy X-ray absorptiometry. A line-cross model was applied to determine QTL regions associated with leg weakness using the QTL Express software.

Results

Eleven QTL affecting leg weakness were identified on eight autosomes. All QTL reached the 5% chromosome-wide significance level. Three QTL were associated with osteochondrosis on the humerus end, two with the fore feet score and two with the rear leg score. QTL on SSC2 and SSC3 influencing bone mineral content and bone mineral density, respectively, reached the 5% genome-wide significance level.

Conclusions

Our results confirm previous studies and provide information on new QTL associated with leg weakness in pigs. These results contribute towards a better understanding of the genetic background of leg weakness in pigs.