Nares Trakooljul

Discovery of Candidate Genes for Muscle Traits Based on GWAS Supported by eQTL-analysis

Biochemical and biophysical processes that take place in muscle under relaxed and stressed conditions depend on the abundance and activity of gene products of metabolic and structural pathways. In livestock at post-mortem, these muscle properties determine aspects of meat quality and are measurable. The conversion of muscle to meat mimics pathological processes associated with muscle ischemia, injury or damage in humans and it is an economic factor in pork production. Linkage, association, and expression analyses independently contributed to the identification of trait-associated molecular pathways and genes. We aim at providing multiple evidences for the role of specific genes in meat quality by integrating a genome-wide association study (GWAS) for meat quality traits and the detection of eQTL based on trait-correlated expressed genes and trait-associated markers. The GWAS revealed 51 and 200 SNPs significantly associated with meat quality in a crossbred Pietrain×(German Landrace×Large White) (Pi×(GL×LW)) and a purebred German Landrace (GL) population, respectively. Most significant SNPs in Pi×(GL×LW) were located on chromosomes (SSC) 4 and 6. The data of 47,836 eQTLs at a significance level of p<10-5 were used to scale down the number candidate genes located in these regions. These SNPs on SSC4 showed association with expression levels of ZNF704, IMPA1, and OXSR1; SSC6 SNPs were associated with expression of SIGLEC10 and PIH1D1. Most significant SNPs in GL were located on SSC6 and associated with expression levels of PIH1D1, SIGLEC10, TBCB, LOC100518735, KIF1B, LOC100514845, and two unknown genes. The abundance of transcripts of these genes in muscle, in turn, is significantly correlated with meat quality traits. We identified several genes with evidence for their candidacy for meat quality arising from the integrative approach of a genome-wide association study and eQTL analysis.

Pre- and post-natal muscle microRNA expression profiles of two pig breeds differing in muscularity

miRNAs regulate the expression of target genes in diverse cellular processes and hence play important roles in physiological processes including developmental timing, patterning, embryogenesis, organogenesis, cell lineage, myogenesis and growth control. A comparative expression analysis of miRNAs expressed in the longissimus dorsi muscle at two prenatal stages (63 and 91 days post-conception (dpc)), and one adult stage (180 days post-natum) in both German Landrace (DL) and Pietrain (Pi) pig breeds was performed using a custom-designed array. During the prenatal stages, miR-199 and the miR-17 families were significantly up-regulated at 63 dpc, whereas miR-1 and miR-133a were overexpressed at 91 dpc. The abundance of several miRNAs was increased in the adult stage compared to 91 dpc including miR-1, miR-133, miR-22(a/b) and miR-29a. Some miRNAs were breed-specific, such as miR-199 and the miR-17 families which were all up-regulated in Pi pigs, while miR-133, miR-181 and miR-214 were up-regulated in DL pigs. Several pathways related to muscle development were enriched with predicted targets for the differentially expressed miRNAs. The dynamic expression and breed-associated regulation of porcine muscle miRNAs suggests a functional role for miRNA-mediated gene regulation during muscle development and phenotypic variations of muscle traits.

Integrated Genome-wide association and hypothalamus eQTL studies indicate a link between the circadian rhythm-related gene PER1 and coping behavior

Animal personality and coping styles are basic concepts for evaluating animal welfare. Struggling response of piglets in so-called backtests early in life reflects their coping strategy. Behavioral reactions of piglets in backtests have a moderate heritability, but their genetic basis largely remains unknown. Here, latency, duration and frequency of struggling attempts during one-minute backtests were repeatedly recorded of piglets at days 5, 12, 19, and 26. A genome-wide association study for backtest traits revealed 465 significant SNPs (FDR ≤ 0.05) mostly located in QTL (quantitative trait locus) regions on chromosome 3, 5, 12 and 16. In order to capture genes in these regions, 37 transcripts with significant SNPs were selected for expressionQTL analysis in the hypothalamus. Eight genes (ASGR1, CPAMD8, CTC1, FBXO39, IL19, LOC100511790, RAD51B, UBOX5) had cis- and five (RANGRF, PER1, PDZRN3, SH2D4B, LONP2) had trans-expressionQTL. In particular, for PER1, with known physiological implications for maintenance of circadian rhythms, a role in coping behavior was evidenced by confirmed association in an independent population. For CTC1 a cis-expression QTL and the consistent relationship of gene polymorphism, mRNA expression level and backtest traits promoted its link to coping style. GWAS and eQTL analyses uncovered positional and functional gene candidates for coping behavior.

MicroRNAs Regulate Cellular ATP Levels by Targeting Mitochondrial Energy Metabolism Genes during C2C12 Myoblast Differentiation.

In our previous study, we identified an miRNA regulatory network involved in energy metabolism in porcine muscle. To better understand the involvement of miRNAs in cellular ATP production and energy metabolism, here we used C2C12 myoblasts, in which ATP levels increase during differentiation, to identify miRNAs modulating these processes. ATP level, miRNA and mRNA microarray expression profiles during C2C12 differentiation into myotubes were assessed. The results suggest 14 miRNAs (miR-423-3p, miR-17, miR-130b, miR-301a/b, miR-345, miR-15a, miR-16a, miR-128, miR-615, miR-1968, miR-1a/b, and miR-194) as cellular ATP regulators targeting genes involved in mitochondrial energy metabolism (Cox4i2, Cox6a2, Ndufb7, Ndufs4, Ndufs5, and Ndufv1) during C2C12 differentiation. Among these, miR-423-3p showed a high inverse correlation with increasing ATP levels. Besides having implications in promoting cell growth and cell cycle progression, its function in cellular ATP regulation is yet unknown. Therefore, miR-423-3p was selected and validated for the function together with its potential target, Cox6a2. Overexpression of miR-423-3p in C2C12 myogenic differentiation lead to decreased cellular ATP level and decreased expression of Cox6a2 compared to the negative control. These results suggest miR-423-3p as a novel regulator of ATP/energy metabolism by targeting Cox6a2.

Muscle Transcriptional Profile Based on Muscle Fiber, Mitochondrial Respiratory Activity, and Metabolic Enzymes

Skeletal muscle is a highly metabolically active tissue that both stores and consumes energy. Important biological pathways that affect energy metabolism and metabolic fiber type in muscle cells may be identified through transcriptomic profiling of the muscle, especially ante mortem. Here, gene expression was investigated in malignant hyperthermia syndrome (MHS)-negative Duroc and Pietrian (PiNN) pigs significantly differing for the muscle fiber types slow-twitch-oxidative fiber (STO) and fast-twitch-oxidative fiber (FTO) as well as mitochondrial activity (succinate-dependent state 3 respiration rate). Longissimus muscle samples were obtained 24 h before slaughter and profiled using cDNA microarrays. Differential gene expression between Duroc and PiNN muscle samples were associated with protein ubiquitination, stem cell pluripotency, amyloid processing, and 3-phosphoinositide biosynthesis and degradation pathways. In addition, weighted gene co-expression network analysis within both breeds identified several co-expression modules that were associated with the proportion of different fiber types, mitochondrial respiratory activity, and ATP metabolism. In particular, Duroc results revealed strong correlations between mitochondrion-associated co-expression modules and STO (r = 0.78), fast-twitch glycolytic fiber (r = -0.98), complex I (r=0.72) and COX activity (r = 0.86). Other pathways in the protein-kinase-activity enriched module were positively correlated with STO (r=0.93), while negatively correlated with FTO (r = -0.72). In contrast to PiNN, co-expression modules enriched in macromolecule catabolic process, actin cytoskeleton, and transcription activator activity were associated with fiber types, mitochondrial respiratory activity, and metabolic enzyme activities. Our results highlight the importance of mitochondria for the oxidative capacity of porcine muscle and for breed-dependent molecular pathways in muscle cell fibers.

Identification of Common Regulators of Genes in Co-Expression Networks Affecting Muscle and Meat Properties

Understanding the genetic contributions behind skeletal muscle composition and metabolism is of great interest in medicine and agriculture. Attempts to dissect these complex traits combine genome-wide genotyping, expression data analyses and network analyses. Weighted gene co-expression network analysis (WGCNA) groups genes into modules based on patterns of co-expression, which can be linked to phenotypes by correlation analysis of trait values and the module eigengenes, i.e. the first principal component of a given module. Network hub genes and regulators of the genes in the modules are likely to play an important role in the emergence of respective traits. In order to detect common regulators of genes in modules showing association with meat quality traits, we identified eQTL for each of these genes, including the highly connected hub genes. Additionally, the module eigengene values were used for association analyses in order to derive a joint eQTL for the respective module. Thereby major sites of orchestrated regulation of genes within trait-associated modules were detected as hotspots of eQTL of many genes of a module and of its eigengene. These sites harbor likely common regulators of genes in the modules. We exemplarily showed the consistent impact of candidate common regulators on the expression of members of respective modules by RNAi knockdown experiments. In fact, Cxcr7 was identified and validated as a regulator of genes in a module, which is involved in the function of defense response in muscle cells. Zfp36l2 was confirmed as a regulator of genes of a module related to cell death or apoptosis pathways. The integration of eQTL in module networks enabled to interpret the differentially-regulated genes from a systems perspective. By integrating genome-wide genomic and transcriptomic data, employing co-expression and eQTL analyses, the study revealed likely regulators that are involved in the fine-tuning and synchronization of genes with trait-associated expression.

Single- and Bayesian Multi-Marker Genome-Wide Association for Haematological Parameters in Pigs.

Haematological traits are important traits that show associations with immune and metabolic status, as well as diseases in humans and animals. Mapping genome regions that affect the blood cell traits can contribute to the identification of genomic features useable as biomarkers for immune, disease and metabolic status. A genome-wide association study (GWAS) was conducted using PorcineSNP60 BeadChips. Single-marker and Bayesian multi-marker approaches were integrated to identify genomic regions and corresponding genes overlapping for both methods. GWAS was performed for haematological traits of 591 German Landrace pig. Heritability estimates for haematological traits were medium to high. In total 252 single SNPs associated with 12 haematological traits were identified (NegLog10 of p-value > 5). The Bayesian multi-marker approach revealed 102 QTL regions across the genome, indicated by 1-Mb windows with contribution to additive genetic variance above 0.5%. The integration of both methods resulted in 24 overlapping QTL regions. This study identified overlapping QTL regions from single- and multi-marker approaches for haematological traits. Identifying candidate genes that affect blood cell traits provides the first step towards the understanding of the molecular basis of haematological phenotypes.

MicroRNA-mRNA regulatory networking fine-tunes the porcine muscle fiber type, muscular mitochondrial respiratory and metabolic enzyme activities

BackgroundMicroRNAs (miRNAs) are small non-coding RNAs that play critical roles in diverse biological processes via regulation of gene expression including in skeletal muscles. In the current study, miRNA expression profile was investigated in longissimus muscle biopsies of malignant hyperthermia syndrome-negative Duroc and Pietrain pigs with distinct muscle metabolic properties in order to explore the regulatory role of miRNAs related to mitochondrial respiratory activity and metabolic enzyme activity in skeletal muscle.ResultsA comparative analysis of the miRNA expression profile between Duroc and Pietrain pigs was performed, followed by integration with mRNA profiles based on their pairwise correlation and computational target prediction. The identified target genes were enriched in protein ubiquitination pathway, stem cell pluripotency and geranylgeranyl diphosphate biosynthesis, as well as skeletal and muscular system development. Next, we analyzed the correlation between individual miRNAs and phenotypical traits including muscle fiber type, mitochondrial respiratory activity, metabolic enzyme activity and adenosine phosphate concentrations, and constructed the regulatory miRNA-mRNA networks associated with energy metabolism. It is noteworthy that miR-25 targeting BMPR2 and IRS1, miR-363 targeting USP24, miR-28 targeting HECW2 and miR-210 targeting ATP5I, ME3, MTCH1 and CPT2 were highly associated with slow-twitch oxidative fibers, fast-twitch oxidative fibers, ADP and ATP concentration suggesting an essential role of the miRNA-mRNA regulatory networking in modulating the mitochondrial energy expenditure in the porcine muscle. In the identified miRNA-mRNA network, a tight relationship between mitochondrial and ubiquitin proteasome system at the level of gene expression was observed. It revealed a link between these two systems contributing to energy metabolism of skeletal muscle under physiological conditions.ConclusionsWe assembled miRNA-mRNA regulatory networks based on divergent muscle properties between different pig breeds and further with the correlation analysis of expressed genes and phenotypic measurements. These complex networks relate to muscle fiber type, metabolic enzyme activity and ATP production and may contribute to divergent muscle phenotypes by fine-tuning the expression of genes. Altogether, the results provide an insight into a regulatory role of miRNAs in muscular energy metabolisms and may have an implication on meat quality and production.

Sex-Specific Muscular Maturation Responses Following Prenatal Exposure to Methylation-Related Micronutrients in Pigs

Supplementation of micronutrients involved in DNA methylation, particularly during pregnancy, is recommended because of its impacts on human health, but further evidence is needed regarding the effects of over-supplementation and differences between sexes. Here, a porcine model was used to assess effects of maternal supplementation with one-carbon-cycle compounds during prenatal and postnatal stages on offspring muscle development. Sows received either a standard diet (CON) or a standard diet supplemented with folate, B6, B12, methionine, choline, and zinc (MET) throughout gestation. Myogenesis-, growth-, and nutrient utilization-related transcript expression was assessed using quantitative PCR. Organismal phenotype and gene expression effects differed significantly between males and females. Male MET-offspring showed increased fetal weight during late pregnancy but decreased live weight postnatally, with compensatory transcriptional responses comprising myogenic key drivers (Pax7, MyoD1, myogenin). In contrast, female weights were unaffected by diet, and mRNA abundances corresponded to a phenotype of cellular reorganization via FABP3, FABP4, SPP1 and Insulin-like Growth Factor-signaling. These findings in an animal model suggest that supplementation during pregnancy with methylation-related micronutrients can promote sex-specific myogenic maturation processes related to organismal growth and muscle metabolism. The usage of maternal dietary supplements should be more carefully considered regarding its ability to promote fetal and postnatal health.

Genetically regulated hepatic transcripts and pathways orchestrate haematological, biochemical and body composition traits.

The liver is the central metabolic organ and exhibits fundamental functions in haematological traits. Hepatic expression, haematological, plasma biochemical, and body composition traits were assessed in a porcine model (n = 297) to establish tissue-specific genetic variations that influence the function of immune-metabolism-correlated expression networks. At FDR (false discovery rate) <1%, more than 3,600 transcripts were jointly correlated (r = |0.22–0.48|) with the traits. Functional enrichment analysis demonstrated common links of metabolic and immune traits. To understand how immune and metabolic traits are affected via genetic regulation of gene expression, eQTLs were assessed. 20517 significant (FDR < 5%) eQTLs for 1401 transcripts were identified, among which 443 transcripts were associated with at least one of the examined traits and had cis-eQTL (such as ACO1 (6.52 × 10−7) and SOD1 (6.41 × 10−30). The present study establishes a comprehensive view of hepatic gene activity which links together metabolic and immune traits in a porcine model for medical research.