UBE3B and ZRANB1 polymorphisms and transcript abundance are associated with water holding capacity of porcine M. longissimus dorsi

Haplotype analysis incorporating ancestral origins identified novel genetic loci associated with chicken body weight using an advanced intercross line

BACKGROUND

The genome-wide association study (GWAS) is a powerful method for mapping quantitative trait loci (QTL). However, standard GWAS can detect only QTL that segregate in the mapping population. Crossing populations with different characteristics increases genetic variability but F2 or back-crosses lack mapping resolution due to the limited number of recombination events. This drawback can be overcome with advanced intercross line (AIL) populations, which increase the number recombination events and provide a more accurate mapping resolution. Recent studies in humans have revealed ancestry-dependent genetic architecture and shown the effectiveness of admixture mapping in admixed populations.

RESULTS

Through the incorporation of line-of-origin effects and GWAS on an F9 AIL population, we identified genes that affect body weight at eight weeks of age (BW8) in chickens. The proposed ancestral-haplotype-based GWAS (testing only the origin regardless of the alleles) revealed three new QTLs on GGA12, GGA15, and GGA20. By using the concepts of ancestral homozygotes (individuals that carry two haplotypes of the same origin) and ancestral heterozygotes (carrying one haplotype of each origin), we identified 632 loci that exhibited high-parent (the heterozygote is better than both parents) and mid-parent (the heterozygote is better than the median of the parents) dominance across 12 chromosomes. Out of the 199 genes associated with BW8, EYA1, PDE1C, and MYC were identified as the best candidate genes for further validation.

CONCLUSIONS

In addition to the candidate genes reported in this study, our research demonstrates the effectiveness of incorporating ancestral information in population genetic analyses, which can be broadly applicable for genetic mapping in populations generated by ancestors with distinct phenotypes and genetic backgrounds. Our methods can benefit both geneticists and biologists interested in the genetic determinism of complex traits.

Effect of Fatty Acids on Backfat Quality in Beijing Black Pigs

The quality of pig backfat affects both pork quality and consumer preferences. Fatty acids (FAs) are crucial in determining the backfat quality. This study assessed the effect of FAs on the backfat quality and identified candidate genes associated with these FAs. The differential fatty acids (DFAs) were compared in pigs with varying backfat firmness and four DFAs-caproic acid, stearic acid, linoleic acid and alpha-linolenic acid-were selected based on T-tests (p < 0.05), fold changes (FC > 2 or FC < 0.5), and variable importance (VIP > 1). Genome-wide association studies on the DFAs and linoleic acid/alpha-linolenic acid ratios in 413 Beijing Black pigs identified 22 single-nucleotide polymorphisms significantly associated with one or more traits. The genes PLPP3, MGLL, CYP27A1 and UBE3C were identified as candidates associated with these traits influencing the backfat quality. These findings enhance our understanding of the backfat quality in Beijing Black pigs and provide a basis for further research.

Phase-specific outcmes of arginine or branched-chain amino acids supplementation in low crude protein diets on performance, nutrient digestibility, and expression of tissue protein synthesis and degradation in broiler chickens infected with mixed Eimeria spp

A 35-d study investigated the impact of dietary supplementation with Arginine (Arg) or branched-chain amino acids (BCAA) of broilers receiving low-protein diets whilst infected with mixed Eimeria species. All birds were given the same starter (d0-10) and finisher (d28-35) diets. The 4 grower diets used were a positive control (PC) with adequate protein (18.5%), a low protein diet (NC;16.5% CP), or the NC supplemented with Arg or BCAA. Supplemental AA was added at 50% above the recommended levels. The treatments were in a 4 × 2 factorial arrangement, with 4 diets, with or without Eimeria inoculation on d14. Birds and feed were weighed after inoculation in phases: prepatent (d14-17), acute (d18-21), recovery (d22-28), and compensatory (d29-35). Ileal digesta, jejunum, and breast tissue were collected on d21, 28, and 35. There was no diet × Eimeria inoculation on growth performance at any phase. Infected birds weighed less and consumed less feed (P < 0.05) in all phases. In the prepatent and acute phases, birds on the Arg diets had higher weight gain (P < 0.05) and lower FCR, similar to PC, when compared to NC and BCAA-fed ones. Infection reduced AA digestibility on d21 and 28 (Met and Cys). However, birds that received supplemental AA had higher digestibility (P < 0.05) of their respective supplemented AA on d 21 only. Infected birds had lower (P < 0.05) BO + AT and higher PEPT1 expression on d21. There was a diet × Eimeria interaction (P = 0.004) on gene expression at d28; 4EBP1 genes were significantly downwardly expressed (P < 0.05) in birds fed Arg diet, irrespective of infection. Infected birds exhibited an upward expression (P < 0.05) of Eef2 on d21 and d28 but experienced a downward expression on d35. Supplemental Arg and BCAA had variable effects on growth performance, apparent ileal AA digestibility, and genes of protein synthesis and degradation, but the effect of Arg on promoting weight gain, irrespective of the Eimeria challenge, was more consistent.

Integration of transcriptome and machine learning to identify the potential key genes and regulatory networks affecting drip loss in pork

Low level of drip loss (DL) is an important quality characteristic of meat with high economic value. However, the key genes and regulatory networks contributing to DL in pork remain largely unknown. To accurately identify the key genes affecting DL in muscles postmortem, 12 Duroc × (Landrace × Yorkshire) pigs with extremely high (n = 6, H group) and low (n = 6, L group) DL at both 24 and 48 h postmortem were selected for transcriptome sequencing. The analysis of differentially expressed genes and weighted gene co-expression network analysis (WGCNA) were performed to find the overlapping genes using the transcriptome data, and functional enrichment and protein-protein interaction (PPI) network analysis were conducted using the overlapping genes. Moreover, we used machine learning to identify the key genes and regulatory networks related to DL based on the interactive genes of the PPI network. Finally, nine potential key genes (IRS1, ESR1, HSPA6, INSR, SPOP, MSTN, LGALS4, MYLK2, and FRMD4B) mainly associated with the MAPK signaling pathway, the insulin signaling pathway, and the calcium signaling pathway were identified, and a single-gene set enrichment analysis (GSEA) was performed to further annotate the functions of these potential key genes. The GSEA results showed that these genes are mainly related to ubiquitin-mediated proteolysis and oxidative reactions. Taken together, our results indicate that the potential key genes influencing DL are mainly related to insulin signaling mediated differences in glycolysis and ubiquitin-mediated changes in muscle structure and improve the understanding of gene expression and regulation related to DL and contribute to future molecular breeding for improving pork quality.

RNA-Seq-based discovery of genetic variants and allele-specific expression of two layer lines and broiler chicken

Recent advances in the selective breeding of broilers and layers have made poultry production one of the fastest-growing industries. In this study, a transcriptome variant calling approach from RNA-seq data was used to determine population diversity between broilers and layers. In total, 200 individuals were analyzed from three different chicken populations (Lohmann Brown (LB), n = 90), Lohmann Selected Leghorn (LSL, n = 89), and Broiler (BR, n = 21). The raw RNA-sequencing reads were pre-processed, quality control checked, mapped to the reference genome, and made compatible with Genome Analysis ToolKit for variant detection. Subsequently, pairwise fixation index (F ST) analysis was performed between broilers and layers. Numerous candidate genes were identified, that were associated with growth, development, metabolism, immunity, and other economically significant traits. Finally, allele-specific expression (ASE) analysis was performed in the gut mucosa of LB and LSL strains at 10, 16, 24, 30, and 60 weeks of age. At different ages, the two-layer strains showed significantly different allele-specific expressions in the gut mucosa, and changes in allelic imbalance were observed across the entire lifespan. Most ASE genes are involved in energy metabolism, including sirtuin signaling pathways, oxidative phosphorylation, and mitochondrial dysfunction. A high number of ASE genes were found during the peak of laying, which were particularly enriched in cholesterol biosynthesis. These findings indicate that genetic architecture as well as biological processes driving particular demands relate to metabolic and nutritional requirements during the laying period shape allelic heterogeneity. These processes are considerably affected by breeding and management, whereby elucidating allele-specific gene regulation is an essential step towards deciphering the genotype to phenotype map or functional diversity between the chicken populations. Additionally, we observed that several genes showing significant allelic imbalance also colocalized with the top 1% of genes identified by the FST approach, suggesting a fixation of genes in cis-regulatory elements.

Whole-genome resequencing reveals genetic diversity, differentiation, and selection signatures of yak breeds/populations in Qinghai, China

The Qinghai Province of China is located in the northeast region of the Qinghai-Tibetan Plateau (QTP) and carries abundant yak genetic resources. Previous investigations of archaeological records, mitochondrial DNA, and Y chromosomal markers have suggested that Qinghai was the major center of yak domestication. In the present study, we examined the genomic diversity, differentiation, and selection signatures of 113 Qinghai yak, including 42 newly sequenced Qinghai yak and 71 publicly available individuals, from nine yak breeds/populations (wild, Datong, Huanhu, Xueduo, Yushu, Qilian, Geermu, Tongde, and Huzhu white) using high-depth whole-genome resequencing data. We observed that most of Qinghai yak breeds/populations have abundant genomic diversity based on four genomic parameters (nucleotide diversity, inbreeding coefficients, linkage disequilibrium decay, and runs of homozygosity). Population genetic structure analysis showed that Qinghai yak have two lineages with two ancestral origins and that nine yak breeds/populations are clustered into three distinct groups of wild yak, Geermu yak, and seven other domestic yak breeds/populations. F _ST values showed moderate genetic differentiation between wild yak, Geermu yak, and the other Qinghai yak breeds/populations. Positive selection signals were detected in candidate genes associated with disease resistance (CDK2AP2, PLEC, and CYB5B), heat stress (NFAT5, HSF1, and SLC25A48), pigmentation (MCAM, RNF26, and BOP1), vision (C1QTNF5, MFRP, and TAX1BP3), milk quality (OPLAH and GRINA), neurodevelopment (SUSD4, INSYN1, and PPP1CA), and meat quality (ZRANB1), using the integrated PI, composite likelihood ratio (CLR), and F _ST methods. These findings offer new insights into the genetic mechanisms underlying target traits in yak and provide important information for understanding the genomic characteristics of yak breeds/populations in Qinghai.

Effects of nicotinamide and sodium butyrate on meat quality and muscle ubiquitination degradation genes in broilers reared at a high stocking density

This study was conducted to investigate the effects of supplementation with nicotinamide (NAM) and sodium butyrate (BA) on meat quality and expression of muscle development genes in broilers reared at a high stocking density. A total of 567, 21-day-old AA broilers were randomly assigned to 5 treatment groups and 2 control groups, with 7 replicates of each group. The control groups included a low stocking density (LD; 12.9 birds/m²) and were fed a basal diet. The treatment groups were kept at a high stocking density (HD; 18.6 birds/m²) and received either a low dose of NAM (50 mg/kg; treatment LN), a high dose of NAM (100 mg/kg; treatment HN), a low dose of BA (500 mg/kg; treatment LB), a high dose of BA (1,000 mg/kg; treatment HB), or a compound supplement (50 mg/kg NAM+500 mg/kg BA; treatment COMB); broilers were reared till 42 D of age. The control groups were kept at HD or at LD (12.9 birds/m²) and were fed a basal diet. The heterophil-to-lymphocyte ratio was significantly higher in the HD control group than that in the LD group; this ratio was significantly lower in treatments LN, HN, HB, and COMB than that in the HD control group. The lightness of breast muscles at 45 min and 24 h after slaughter was significantly higher in the HD group than that in the LD group, and the HD group showed the highest drip loss at 24 h and 48 h. Lightness and drip loss were lower in the HN, LB, and COMB treatments than those in the HD group. HD rearing significantly reduced gene expression of myogenic regulatory factor 5 (MYF5) while significantly increased expression of the protein ubiquitin degradation genes FBXO9, FBXO22, and FBXO32. All treatments significantly reduced FBXO9 and FBXO32 expression. Our results suggest dietary supplementation with NAM and BA can improve meat quality of broilers under high stocking density by upregulating the expression of myogenic genes, and inhibiting protein ubiquitination.

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

Background

MicroRNAs (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.

Results

A 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.

Conclusions

We 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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2850-8) contains supplementary material, which is available to authorized users.

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.

Correlated mRNAs and miRNAs from co-expression and regulatory networks affect porcine muscle and finally meat properties

Background

Physiological processes aiding the conversion of muscle to meat involve many genes associated with muscle structure and metabolic processes. MicroRNAs regulate networks of genes to orchestrate cellular functions, in turn regulating phenotypes.

Results

We applied weighted gene co-expression network analysis to identify co-expression modules that correlated to meat quality phenotypes and were highly enriched for genes involved in glucose metabolism, response to wounding, mitochondrial ribosome, mitochondrion, and extracellular matrix. Negative correlation of miRNA with mRNA and target prediction were used to select transcripts out of the modules of trait-associated mRNAs to further identify those genes that are correlated with post mortem traits.

Conclusions

Porcine muscle co-expression transcript networks that correlated to post mortem traits were identified. The integration of miRNA and mRNA expression analyses, as well as network analysis, enabled us to interpret the differentially-regulated genes from a systems perspective. Linking co-expression networks of transcripts and hierarchically organized pairs of miRNAs and mRNAs to meat properties yields new insight into several biological pathways underlying phenotype differences. These pathways may also be diagnostic for many myopathies, which are accompanied by deficient nutrient and oxygen supply of muscle fibers.