Association and expression quantitative trait loci (eQTL) analysis of porcine AMBP, GC and PPP1R3B genes with meat quality traits

Association and expression study of SEPW1 and JAML as preliminary candidate genes related to lamb odor and flavor

The current study aimed to identify the SEPW1 and JAML genes in lamb as candidate genes related to lamb odor and flavor. The polymorphism study showed that the SEPW1 gene was polymorphic at the BanI restriction site with three genotypes (AA, AG, and GG), whereas the JAML gene was monomorphic at HhaI with genotype (GG). The association of SEPW1 between genotype and lamb odor and flavor (BCFAs and skatole) was analyzed using GLM (General Linear Model). MNA (4-methylnonanoic) was significantly associated (p < 0.05) with lamb odor and flavor. AA genotype has a lower level of MNA than AG and GG, while MOA (4-methyloctanoic), EOA(4-ethyloctanoic), MI (3-methylindole) and MP (3-methylphenol) was not significantly associated with lamb odor and flavor (p > 0.05). Furthermore, to analyze the mRNA expression of SEPW1 in liver tissues, the lambs were divided into three groups based on the genotypes AA, AG, and GG, however, mRNA expression was not differentially expressed between AA, AG, and GG (p > 0.05). These results will enhance the understanding of the functions of SEPW1 gene relation to odor and flavor traits and will shed light on the polymorphism of SEPW1 gene in lamb as a candidate gene for reducing MNA in lamb.

A look under the hood of genomic-estimated breed compositions for brangus cattle: What have we learned?

The Brangus cattle were developed to utilize the superior traits of Angus and Brahman cattle. Their genetic compositions are expected to be stabilized at 3/8 Brahman and 5/8 Angus. Previous studies have shown more than expected Angus lineage with Brangus cattle, and the reasons are yet to be investigated. In this study, we revisited the breed compositions for 3,605 Brangus cattle from three perspectives: genome-wise (GBC), per chromosomes (CBC), and per chromosome segments (SBC). The former (GBC) depicted an overall picture of the “mosaic” genome of the Brangus attributable to their ancestors, whereas the latter two criteria (CBC and SBC) corresponded to local ancestral contributions. The average GBC for the 3,605 Brangus cattle were 70.2% Angus and 29.8% Brahman. The K-means clustering supported the postulation of the mixture of 1/2 Ultrablack (UB) animals in Brangus. For the non-UB Brangus animals, the average GBC were estimated to be 67.4% Angus and 32.6% Brahman. The 95% confidence intervals of their overall GBC were 60.4%–73.5% Angus and 26.5%–39.6% Brahman. Possibly, genetic selection and drifting have resulted in an approximately 5% average deviation toward Angus lineage. The estimated ancestral contributions by chromosomes were heavily distributed toward Angus, with 27 chromosomes having an average Angus CBC greater than 62.5% but only two chromosomes (5 and 20) having Brahman CBC greater than 37.5%. The chromosomal regions with high Angus breed proportions were prevalent, tending to form larger blocks on most chromosomes. In contrast, chromosome segments with high Brahman breed proportion were relatively few and isolated, presenting only on seven chromosomes. Hence, genomic hitchhiking effects were strong where Angus favorable alleles resided but weak where Brahman favorable alleles were present. The functions of genes identified in the chromosomal regions with high (≥75%) Angus compositions were diverse yet may were related to growth and body development. In contrast, the genes identified in the regions with high (≥37.5%) Brahman compositions were primarily responsible for disease resistance. In conclusion, we have addressed the questions concerning the Brangus genetic make-ups. The results can help form a dynamic picture of the Brangus breed formation and the genomic reshaping.

Hepatic Transcriptome Analysis Reveals Genes, Polymorphisms, and Molecules Related to Lamb Tenderness

Tenderness is a key meat quality trait that determines the public acceptance of lamb consumption, so genetic improvement toward lamb with higher tenderness is pivotal for a sustainable sheep industry. However, unravelling the genomics controlling the tenderness is the first step. Therefore, this study aimed to identify the transcriptome signatures and polymorphisms related to divergent lamb tenderness using RNA deep sequencing. Since the molecules and enzymes that control muscle growth and tenderness are metabolized and synthesized in the liver, hepatic tissues of ten sheep with divergent phenotypes: five high- and five low-lamb tenderness samples were applied for deep sequencing. Sequence analysis identified the number of reads ranged from 21.37 to 25.37 million bases with a mean value of 22.90 million bases. In total, 328 genes are detected as differentially expressed (DEGs) including 110 and 218 genes that were up- and down-regulated, respectively. Pathway analysis showed steroid hormone biosynthesis as the dominant pathway behind the lamb tenderness. Gene expression analysis identified the top high (such as TP53INP1, CYP2E1, HSD17B13, ADH1C, and LPIN1) and low (such as ANGPTL2, IGFBP7, FABP5, OLFML3, and THOC5) expressed candidate genes. Polymorphism and association analysis revealed that mutation in OLFML3, ANGPTL2, and THOC5 genes could be potential candidate markers for tenderness in sheep. The genes and pathways identified in this study cause variation in tenderness, thus could be potential genetic markers to improve meat quality in sheep. However, further validation is needed to confirm the effect of these markers in different sheep populations so that these could be used in a selection program for lamb with high tenderness.

Transcriptome-based analysis of early post-mortem formation of pale, soft, and exudative (PSE) pork

Pale, soft, and exudative (PSE) meat can cause consumer dissatisfaction and economic losses. This study determined meat quality, glycolytic enzyme activity, and differential gene expression in the longissimus lumborum (LL) and semimembranosus (SM) of normal and PSE pork carcasses. The SM did not result in PSE meat. Hexokinase, lactate dehydrogenase, and pyruvate kinase activities were lower in the SM of PSE carcasses than in the normal carcasses. Functional enrichment analysis revealed that immune, inflammatory, and muscle fibre genes were significantly enriched in PSE pork. More specifically, PPP1R3G and MSS51 may be key genes regulating pork quality in the SM. Meanwhile, the differential expression of PLVAB, ADIPOQ, LEP, MYH4, MYH7, MYL3, MYL6B, FOS, ATF3, and HSPA6 may induce PSE formation in the LL. These results may provide insights into PSE pork formation mechanisms and reveal candidate genes for improving meat quality after validation.

Association study and expression analysis of olfactomedin like 3 gene related to meat quality, carcass characteristics, retail meat cut, and fatty acid composition in sheep

Objective

The objective of this study was to identify polymorphism in olfactomedin like 3 (OLFML3) gene, and association analysis with meat quality, carcass characteristics, retail meat cut, and fatty acid composition in sheep, and expression quantification of OLFML3 gene in phenotypically divergent sheep.

Methods

A total of 328 rams at the age of 10 to 12 months with an average body weight of 26.13 kg were used. A novel polymorphism was identified using high-throughput sequencing in sheep and genotyping of OLFML3 polymorphism was performed using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Among 328 rams, 100 rams representing various sheep genotypes were used for association study and proc general linear model was used to analyse association between genotypes and phenotypic traits. Quantitative real-time polymerase chain reaction (qRT-PCR) was used for the expression analysis of OLFML3 mRNA in phenotypically divergent sheep population.

Results

The findings revealed a novel polymorphism in the OLFML3 gene (g.90317673 C>T). The OLFML3 gene revealed three genotypes: CC, CT, and TT. The single nucleotide polymorphism (SNP) was found to be significantly (p<0.05) associated with meat quality traits such as tenderness and cooking loss; carcass characteristics such as carcass length; retail meat cut such as pelvic fat in leg, intramuscular fat in loin and tenderloin, muscle in flank and shank; fatty acids composition such as tridecanoic acid (C13:0), palmitoleic acid (C16:1), heptadecanoic acid (C17:0), ginkgolic acid (C17:1), linolenic acid (C18:3n3), arachidic acid (C20:0), eicosenoic acid (C20:1), arachidonic acid (C20:4n6), heneicosylic acid (C21:0), and nervonic acid (C24:1). The TT genotype was associated with higher level of meat quality, carcass characteristics, retail meat cut, and some fatty acids composition. However, the mRNA expression analysis was not different among genotypes.

Conclusion

The OLFML3 gene could be a potential putative candidate for selecting higher quality sheep meat, carcass characteristics, retail meat cuts, and fatty acid composition in sheep.

Identification of Genes Related to Squab Muscle Growth and Lipid Metabolism from Transcriptome Profiles of Breast Muscle and Liver in Domestic Pigeon (Columba livia)

Simple Summary

Domestic pigeon is an important small poultry species raised for high-quality meat production. However, the relevant gene associated with meat growth and lipid metabolism during the period from dehulling to marketing are not known. Therefore, we aim to identify genes related to squab muscle growth and lip metabolism from transcriptome profiles of breast muscle and liver in domestic pigeon. In this study, we totally found that 4465 differentially expressed genes (DEGs) identified in the breast muscle and liver libraries, which include 2585 genes that were up-regulated and 2122 genes that were down-regulated. Most genes are involved in cell proliferation and differentiation, lipid metabolism and energy metabolism according to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs. We also detected 16 DEGs to verify data from RNA-Seq and qPCR, which were consistent in quantitatively estimating the transcription levels of the tested transcripts by qPCR analysis. The results of this study will lay the foundation for understanding the mechanisms of muscle growth and lipid metabolism in domestic pigeons.

Abstract

The improvements in muscle growth rate and meat quality are the major breeding aims in pigeon industry. Liver and muscle are recognized as important sites for fatty acid metabolism; understanding the role of specific transcripts in the breast muscle and liver might lead to the elucidation of interrelated biological processes. In this study, RNA-Sequencing (RNA-Seq) was applied to compare the transcriptomes of breast muscle and liver tissues among pigeons at five developmental periods (0, 1, 2, 3, 4 weeks post-hatching) to identify candidate genes related to muscle growth and lipid metabolism. There were 3142 differentially expressed genes (DEGs) identified in the breast muscle libraries; 1794 genes were up-regulated while 1531 genes were down-regulated. A total of 1323 DEGs were acquired from the liver libraries, with 791 up-regulated genes and 591 down-regulated genes. By pathway enrichment analysis, a set of significantly enriched pathways were identified for the DEGs, which are potentially involved in cell proliferation and differentiation, lipid metabolism and energy metabolism in pigeon breast muscle and liver. Our results are consistent with previous partial reports from domestic animals and poultry and provide some unidentified genes involved in muscle growth and lipid metabolism. The reliability of the sequencing data was verified through qPCR analysis of 16 genes from eight comparison groups (two genes per group). The findings from this study could contribute to future investigations of muscle growth and lipid metabolism mechanisms and establish molecular approaches to improve muscle growth rate and meat quality in domestic pigeon breeding.

Hepatic transcriptome analysis identifies genes, polymorphisms and pathways involved in the fatty acids metabolism in sheep

Fatty acids (FA) in ruminants, especially unsaturated FA (USFA) have important impact in meat quality, nutritional value, and flavour quality of meat, and on consumer’s health. Identification of the genetic factors controlling the FA composition and metabolism is pivotal to select sheep that produce higher USFA and lower saturated (SFA) for the benefit of sheep industry and consumers. Therefore, this study was aimed to investigate the transcriptome profiling in the liver tissues collected from sheep with divergent USFA content in longissimus muscle using RNA deep-sequencing. From sheep (n = 100) population, liver tissues with higher (n = 3) and lower (n = 3) USFA content were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample were ranged from 21.28 to 28.51 million with a median of 23.90 million. Approximately, 198 genes were differentially regulated with significance level of p-adjusted value <0.05. Among them, 100 genes were up-regulated, and 98 were down-regulated (p<0.01, FC>1.5) in the higher USFA group. A large proportion of key genes involved in FA biosynthesis, adipogenesis, fat deposition, and lipid metabolism were identified, such as APOA5, SLC25A30, GFPT1, LEPR, TGFBR2, FABP7, GSTCD, and CYP17A. Pathway analysis revealed that glycosaminoglycan biosynthesis- keratan sulfate, adipokine signaling, galactose metabolism, endocrine and other factors-regulating calcium metabolism, mineral metabolism, and PPAR signaling pathway were playing important regulatory roles in FA metabolism. Importantly, polymorphism and association analyses showed that mutation in APOA5, CFHR5, TGFBR2 and LEPR genes could be potential markers for the FA composition in sheep. These polymorphisms and transcriptome networks controlling the FA variation could be used as genetic markers for FA composition-related traits improvement. However, functional validation is required to confirm the effect of these SNPs in other sheep population in order to incorporate them in the sheep breeding program.

Genetic Diversity and Identification of Homozygosity-Rich Genomic Regions in Seven Italian Heritage Turkey (Meleagris gallopavo) Breeds

Italian autochthonous turkey breeds are an important reservoir of genetic biodiversity that should be maintained with an in vivo approach. The aim of this study, part of the TuBAvI national project on biodiversity, was to use run of homozygosity (ROH), together with others statistical approaches (e.g., Wright's F-statistics, principal component analysis, ADMIXTURE analysis), to investigate the genomic diversity in several heritage turkey breeds. We performed a genome-wide characterization of ROH-rich regions in seven autochthonous turkey breeds, i.e., Brianzolo (Brzl), Bronzato Comune Italiano (BrCI), Bronzato dei Colli Euganei (CoEu), Parma e Piacenza (PrPc), Nero d'Italia (NeIt), Ermellinato di Rovigo (ErRo) and Romagnolo (Roma). ROHs were detected based on a 650K SNP genotyping. ROH_islands were identified as homozygous ROH regions shared by at least 75% of birds (within breed). Annotation of genes was performed with DAVID. The admixture analyses revealed that six breeds are unique populations while the Roma breed consists in an admixture of founder populations. Effective population size estimated on genomic data shows a numeric contraction. ROH_islands harbour genes that may be interesting for target selection in commercial populations also. Among them the PTGS2 and PLA2G4A genes on chr10 were related to reproduction efficiency. This is the first study mapping genetic variation in autochthonous turkey populations. Breeds were genetically different among them, with the Roma breed proving to be a mixture of the other breeds. The ROH_islands identified harboured genes peculiar to the selection that occurred in heritage breeds. Finally, this study releases previously undisclosed information on existing genetic variation in the turkey species.

Genetic Architecture of Carcass and Meat Quality Traits in Montana Tropical® Composite Beef Cattle

The Montana Tropical^® Composite is a recently developed beef cattle population that is rapidly expanding in Brazil and other tropical countries. This is mainly due to its improved meat quality and adaptation to tropical climate conditions compared to Zebu and Taurine cattle breeds, respectively. This study aimed to investigate the genetic architecture of ultrasound-based carcass and meat quality traits in Montana Tropical^® Composite beef cattle. Therefore, we estimated variance components and genetic parameters and performed genome-wide association studies using the weighted single-step Genomic Best Linear Unbiased Prediction (GBLUP) approach. A pedigree dataset containing 28,480 animals was used, in which 1,436 were genotyped using a moderate-density Single Nucleotide Polymorphism panel (30K; 30,105 SNPs). A total of 9,358, 5,768, 7,996, and 1,972 phenotypic records for the traits Longissimus muscle area (LMA), backfat thickness (BFT), rump fat thickness (RFT), and for marbling score (MARB), respectively, were used for the analyses. Moderate to high heritability estimates were obtained and ranged from 0.16 ± 0.03 (RFT) to 0.33 ± 0.05 (MARB). A high genetic correlation was observed between BFT and RFT (0.97 ± 0.02), suggesting that a similar set of genes affects both traits. The most relevant genomic regions associated with LMA, BFT, RFT, and MARB were found on BTA10 (5.4–5.8 Mb), BTA27 (25.2–25.5 Mb), BTA18 (60.6–61.0 Mb), and BTA21 (14.8–15.4 Mb). Two overlapping genomic regions were identified for RFT and MARB (BTA13:47.9–48.1 Mb) and for BFT and RFT (BTA13:61.5–62.3 Mb). Candidate genes identified in this study, including PLAG1, LYN, WWOX, and PLAGL2, were previously reported to be associated with growth, stature, skeletal muscle growth, fat thickness, and fatty acid composition. Our results indicate that ultrasound-based carcass and meat quality traits in the Montana Tropical^® Composite beef cattle are heritable, and therefore, can be improved through selective breeding. In addition, various novel and already known genomic regions related to these traits were identified, which contribute to a better understanding of the underlying genetic background of LMA, BFT, RFT, and MARB in the Montana Tropical Composite population.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles

Background

Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles.

Results

A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110–114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq.

Conclusions

We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5557-9) contains supplementary material, which is available to authorized users.

RNA-Seq transcriptome analysis of breast muscle in Pekin ducks supplemented with the dietary probiotic Clostridium butyricum

BACKGROUND

Increased attention is being paid to breast muscle yield and meat quality in the duck breeding industry. Our previous report has demonstrated that dietary Clostridium butyricum (C. butyricum) can improve meat quality of Pekin ducks. However, the potential biological processes and molecular mechanisms that are modulated by dietary C. butyricum in the breast muscle of Pekin ducks remain unknown.

RESULTS

Supplementation with C. butyricum increased growth performance and meat yield. Therefore, we utilized de novo assembly methods to analyze the RNA-Seq transcriptome profiles in breast muscle to explore the differentially expressed genes between C. butyricum-treated and control Pekin ducks. A total of 1119 differentially expressed candidate genes were found of which 403 genes were significantly up-regulated and 716 genes were significantly down-regulated significantly. qRT-PCR analysis was used to confirm the accuracy of the of RNA-Seq results. GO annotations revealed potential genes, processes and pathways that may participate in meat quality and muscle development. KEGG pathway analysis showed that the differentially expressed genes participated in numerous pathways related to muscle development, including ECM-receptor interaction, the MAPK signaling pathway and the TNF signaling pathway.

CONCLUSIONS

This study suggests that long-time dietary supplementation with C. butyricum can modulate muscle development and meat quality via altering the expression patterns of genes involved in crucial metabolic pathways. The findings presented here provide unique insights into the molecular mechanisms of muscle development in Pekin ducks in response to dietary C. butyricum.

Deciphering the genetic regulation of peripheral blood transcriptome in pigs through expression genome-wide association study and allele-specific expression analysis

BACKGROUND

Efforts to improve sustainability in livestock production systems have focused on two objectives: investigating the genetic control of immune function as it pertains to robustness and disease resistance, and finding predictive markers for use in breeding programs. In this context, the peripheral blood transcriptome represents an important source of biological information about an individual's health and immunological status, and has been proposed for use as an intermediate phenotype to measure immune capacity. The objective of this work was to study the genetic architecture of variation in gene expression in the blood of healthy young pigs using two approaches: an expression genome-wide association study (eGWAS) and allele-specific expression (ASE) analysis.

RESULTS

The blood transcriptomes of 60-day-old Large White pigs were analyzed by expression microarrays for eGWAS (242 animals) and by RNA-Seq for ASE analysis (38 animals). Using eGWAS, the expression levels of 1901 genes were found to be associated with expression quantitative trait loci (eQTLs). We recovered 2839 local and 1752 distant associations (Single Nucleotide Polymorphism or SNP located less or more than 1 Mb from expression probe, respectively). ASE analyses confirmed the extensive cis-regulation of gene transcription in blood, and revealed allelic imbalance in 2286 SNPs, which affected 763 genes. eQTLs and ASE-genes were widely distributed on all chromosomes. By analyzing mutually overlapping eGWAS results, we were able to describe putative regulatory networks, which were further refined using ASE data. At the functional level, genes with genetically controlled expression that were detected by eGWAS and/or ASE analyses were significantly enriched in biological processes related to RNA processing and immune function. Indeed, numerous distant and local regulatory relationships were detected within the major histocompatibility complex region on chromosome 7, revealing ASE for most class I and II genes.

CONCLUSIONS

This study represents, to the best of our knowledge, the first genome-wide map of the genetic control of gene expression in porcine peripheral blood. These results represent an interesting resource for the identification of genetic markers and blood biomarkers associated with variations in immunity traits in pigs, as well as any other complex traits for which blood is an appropriate surrogate tissue.

Expression of the SIRT2 gene and its relationship with body size traits in Qinchuan cattle (Bos taurus)

Silent information regulator 2 (SIRT2) is a member of the sirtuin family of class III NAD (nicotinamide adenine dinucleotide)-dependent protein deacetylases and may regulate senescence, metabolism and apoptosis. The aims of this study were to investigate whether the SIRT2 gene could be used as a candidate gene in the breeding of Qinchuan cattle. Real-time polymerase chain reaction (RT-PCR) results showed that among all types of tissue that were analyzed, the highest mRNA expression levels of the gene were found in subcutaneous fat. DNA sequencing of 468 individual Qinchuan cattle identified two novel, single nucleotide polymorphisms (g.19501 C > T and g.19518 C > T) in the 3' untranslated region (3'UTR) of the SIRT2 gene. The frequencies of SNP g.19501 C > T and g.19518 C > T were in Hardy-Weinberg disequilibrium in all the samples (chi-square test, χ² < χ_0.05²). An association analysis showed that the two loci were significantly correlated with some body size traits and the H₂H₂ (-CT-CT-) diplotypes performed better than other combinations. These results indicated that the variations in the SIRT2 gene and their corresponding genotypes may be considered as molecular markers for economic traits in cattle breeding.

Molecular characterization, expression profiles, and analysis of Qinchuan cattle SIRT1 gene association with meat quality and body measurement traits (Bos taurus)

Silent information regulator (SIRT1), was closely associated with senescence, metabolism, and apoptosis. The objective of this study was to investigate whether SIRT1 gene could be used as a candidate gene in the breeding process of Qinchuan cattle. Via sequencing technology conducted in 453 individuals of Qinchuan cattle, single nucleotide polymorphisms (G25764A, A25846G, and T25868C) with 5 haplotypes and 6 combined genotypes in 3' untranslated region of SIRT1 gene were identified. In addition, three loci were significantly associated with some of the body measurements and meat quality traits in Qinchuan cattle (P < 0.05), and the H2H2 (GG-AA-CC) diplotypes had better performance than other combinations in Qinchuan cattle. These results suggest that the SIRT1 gene could be used in marker assisted selection to improve the production traits of Qinchuan cattle.

A genetical genomics approach reveals new candidates and confirms known candidate genes for drip loss in a porcine resource population

In this study lean meat water-holding capacity (WHC) of a Duroc × Pietrain (DuPi) resource population with corresponding genotypes and transcriptomes was investigated using genetical genomics. WHC was characterized by drip loss measured in M. longissimus dorsi. The 60K Illumina SNP chips identified genotypes of 169 F2 DuPi animals. Whole-genome transcriptomes of muscle samples were available for 132 F2 animals using the Affymetrix 24K GeneChip® Porcine Genome Array. Performing genome-wide association studies of transcriptional profiles, which are correlated with phenotypes, allows elucidation of cis- and trans-regulation. Expression levels of 1,228 genes were significantly correlated with drip loss and were further analyzed for enrichment of functional annotation groups as defined by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. A hypergeometric gene set enrichment test was performed and revealed glycolysis/glyconeogenesis, pentose phosphate pathway, and pyruvate metabolism as the most promising pathways. For 267 selected transcripts, expression quantitative trait loci (eQTL) analysis was performed and revealed a total of 1,541 significant associations. Because of positional accordance of the gene underlying transcript and the eQTL location, it was possible to identify eight eQTL that can be assumed to be cis-regulated. Comparing the results of gene set enrichment and the eQTL detection tests, molecular networks and potential candidate genes, which seemed to play key roles in the expression of WHC, were detected. The α-1-microglobulin/bikunin precursor (AMBP) gene was assumed to be cis-regulated and was part of the glycolysis pathway. This approach supports the identification of trait-associated SNPs and the further biological understanding of complex traits.

Identification of the novel candidate genes and variants in boar liver tissues with divergent skatole levels using RNA deep sequencing

Boar taint is the unpleasant odour of meat derived from non-castrated male pigs, caused by the accumulation of androstenone and skatole in fat. Skatole is a tryptophan metabolite produced by intestinal bacteria in gut and catabolised in liver. Since boar taint affects consumer's preference, the aim of this study was to perform transcriptome profiling in liver of boars with divergent skatole levels in backfat by using RNA-Seq. The total number of reads produced for each liver sample ranged from 11.8 to 39.0 million. Approximately 448 genes were differentially regulated (p-adjusted <0.05). Among them, 383 genes were up-regulated in higher skatole group and 65 were down-regulated (p<0.01, FC>1.5). Differentially regulated genes in the high skatole liver samples were enriched in metabolic processes such as small molecule biochemistry, protein synthesis, lipid and amino acid metabolism. Pathway analysis identified the remodeling of epithelial adherens junction and TCA cycle as the most dominant pathways which may play important roles in skatole metabolism. Differential gene expression analysis identified candidate genes in ATP synthesis, cytochrome P450, keratin, phosphoglucomutase, isocitrate dehydrogenase and solute carrier family. Additionally, polymorphism and association analysis revealed that mutations in ATP5B, KRT8, PGM1, SLC22A7 and IDH1 genes could be potential markers for skatole levels in boars. Furthermore, expression analysis of exon usage of three genes (ATP5B, KRT8 and PGM1) revealed significant differential expression of exons of these genes in different skatole levels. These polymorphisms and exon expression differences may have impacts on the gene activity ultimately leading to skatole variation and could be used as genetic marker for boar taint related traits. However, further validation is required to confirm the effect of these genetic markers in other pig populations in order to be used in genomic selection against boar taint in pig breeding programs.

RNA deep sequencing reveals novel candidate genes and polymorphisms in boar testis and liver tissues with divergent androstenone levels

Boar taint is an unpleasant smell and taste of pork meat derived from some entire male pigs. The main causes of boar taint are the two compounds androstenone (5α-androst-16-en-3-one) and skatole (3-methylindole). It is crucial to understand the genetic mechanism of boar taint to select pigs for lower androstenone levels and thus reduce boar taint. The aim of the present study was to investigate transcriptome differences in boar testis and liver tissues with divergent androstenone levels using RNA deep sequencing (RNA-Seq). The total number of reads produced for each testis and liver sample ranged from 13,221,550 to 33,206,723 and 12,755,487 to 46,050,468, respectively. In testis samples 46 genes were differentially regulated whereas 25 genes showed differential expression in the liver. The fold change values ranged from −4.68 to 2.90 in testis samples and −2.86 to 3.89 in liver samples. Differentially regulated genes in high androstenone testis and liver samples were enriched in metabolic processes such as lipid metabolism, small molecule biochemistry and molecular transport. This study provides evidence for transcriptome profile and gene polymorphisms of boars with divergent androstenone level using RNA-Seq technology. Digital gene expression analysis identified candidate genes in flavin monooxygenease family, cytochrome P450 family and hydroxysteroid dehydrogenase family. Moreover, polymorphism and association analysis revealed mutation in IRG6, MX1, IFIT2, CYP7A1, FMO5 and KRT18 genes could be potential candidate markers for androstenone levels in boars. Further studies are required for proving the role of candidate genes to be used in genomic selection against boar taint in pig breeding programs.

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.