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

Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds

In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.

Identification of core genes affecting IMF deposition in bovine

Intramuscular fat (IMF) content is an important economic factor in beef production. However, knowledge on the key factors controlling bovine IMF is limited. In this study, using weighted gene co-expression network analysis (WGCNA), nine modules were identified and the number of transcripts in these modules ranged from 36 to 3191. Two modules were found to be significantly associated with fat deposition and three genes (TCAP, MYH7, and TNNC1) were further identified by Protein-protein interaction (PPI), which may be the hub genes regulating bovine IMF deposition. In addition, considering the genetic variation, the PCK1 gene was found by functional enrichment analysis of overlapping genes, which was previously reported to be involved in IMF deposition. We noted that the core promoter region of buffalo PCK1 binds to transcription factors involved in lipid metabolism while cattle PCK1 binds transcription factors involved in muscle development. The results suggest that PCK1 participated in IMF deposition of buffalo and cattle in different ways. In summary, gene expression networks and new candidate genes associated with IMF deposition identified in this study. This would lay the foundation for further research into the molecular regulatory mechanisms underlying bovine IMF deposition.

Integrative analysis of microRNAs and mRNAs reveals the regulatory networks of triterpenoid saponin metabolism in Soapberry (Sapindus mukorossi Gaertn.)

Triterpenoid saponin are important secondary metabolites and bioactive constituents of soapberry (Sapindus mukorossi Gaertn.) and are widely used in medicine and toiletry products. However, little is known about the roles of miRNAs in the regulation of triterpenoid saponin biosynthesis in soapberry. In this study, a total of 3036 miRNAs were identified, of which 1372 miRNAs were differentially expressed at different stages of pericarp development. Important KEGG pathways, such as terpenoid backbone biosynthesis, sesquiterpenoid and triterpenoid biosynthesis, and basal transcription factors were highlighted, as well the roles of some key miRNAs, such as ath-miR5021, han-miR3630-3p, and ppe-miR858, which may play important roles in regulating triterpenoid saponin biosynthesis. In addition, 58 miRNAs might participate in saponin biosynthesis pathways by predicting the targets of those miRNAs to 53 saponin biosynthesis structural genes. And 75 miRNAs were identified to potentially play vital role in saponin accumulation by targeting transcript factor genes, bHLH, bZIP, ERF, MYB, and WRKY, respectively, which are candidate regulatory genes in the pathway of saponin biosynthesis. The results of weighted gene coexpression network analysis (WGCNA) suggested that two saponin-specific miRNA modules and 10 hub miRNAs may participate in saponin biosynthesis. Furthermore, multiple miRNA-mRNA regulatory networks potentially involved in saponin biosynthesis were generated, e.g., ath-miR5021-SmIDI2/SmGPS5/SmbAS1/SmCYP71D-3/SmUGT74G-2, han-miR3630-3p-SmCYP71A-14/SmbHLH54/SmMYB135/SmWRKY32, and ppe-miR858-SmMYB5/SmMYB32. qRT-PCR analysis validated the expression patterns of nine miRNAs and 12 corresponding target genes. This study represents the first comprehensive analysis of miRNAs in soapberry and lays the foundation for further understanding of miRNA-based regulation in triterpenoid saponin biosynthesis.

Identification of Key Genes and Biological Pathways Associated with Skeletal Muscle Maturation and Hypertrophy in Bos taurus, Ovis aries, and Sus scrofa

The aim of the current study was to identify the major genes and pathways involved in the process of hypertrophy and skeletal muscle maturation that is common for Bos taurus, Ovis aries, and Sus scrofa species. Gene expression profiles related to Bos taurus, Ovis aries, and Sus scrofa muscle, with accession numbers GSE44030, GSE23563, and GSE38518, respectively, were downloaded from the GEO database. Differentially expressed genes (DEGs) were screened out using the Limma package of R software. Genes with Fold Change > 2 and an adjusted p-value < 0.05 were identified as significantly different between two treatments in each species. Subsequently, gene ontology and pathway enrichment analyses were performed. Moreover, hub genes were detected by creating a protein−protein interaction network (PPI). The results of the analysis in Bos taurus showed that in the period of 280 dpc−3-months old, a total of 1839 genes showed a significant difference. In Ovis aries, however, during the period of 135dpc−2-months old, a total of 486 genes were significantly different. Additionally, in the 91 dpc−adult period, a total of 2949 genes were significantly different in Sus scrofa. The results of the KEGG pathway enrichment analysis and GO function annotation in each species separately revealed that in Bos taurus, DEGs were mainly enriched through skeletal muscle fiber development and skeletal muscle contraction, and the positive regulation of fibroblast proliferation, positive regulation of skeletal muscle fiber development, PPAR signaling pathway, and HIF-1 signaling pathway. In Ovis aries, DEGs were mainly enriched through regulating cell growth, skeletal muscle fiber development, the positive regulation of fibroblast proliferation, skeletal muscle cell differentiation, and the PI3K-Akt signaling, HIF-1 signaling, and Rap1 signaling pathways. In Sus scrofa, DEGs were mainly enriched through regulating striated muscle tissue development, the negative regulation of fibroblast proliferation and myoblast differentiation, and the HIF-1 signaling, AMPK signaling, and PI3K-Akt signaling pathways. Using a Venn diagram, 36 common DEGs were identified between Bos taurus, Ovis aries, and Sus scrofa. A biological pathways analysis of 36 common DEGs in Bos taurus, Ovis aries, and Sus scrofa allowed for the identification of common pathways/biological processes, such as myoblast differentiation, the regulation of muscle cell differentiation, and positive regulation of skeletal muscle fiber development, that orchestrated the development and maturation of skeletal muscle. As a result, hub genes were identified, including PPARGC1A, MYOD1, EPAS1, IGF2, CXCR4, and APOA1, in all examined species. This study provided a better understanding of the relationships between genes and their biological pathways in the skeletal muscle maturation process.

Effects of hypoxia and reoxygenation on mitochondrial functions and transcriptional profiles of isolated brain and muscle porcine cells

Oxygen fluctuations might occur in mammalian tissues under physiological (e.g. at high altitudes) or pathological (e.g. ischemia-reperfusion) conditions. Mitochondria are the key target and potential amplifiers of hypoxia-reoxygenation (H-R) stress. Understanding the mitochondrial responses to H-R stress is important for identifying adaptive mechanisms and potential therapeutic solutions for pathologies associated with oxygen fluctuations. We explored metabolic response to H-R stress in two tissue types (muscle and brain) with different degrees of hypoxia tolerance in a domestic pig Sus scrofa focusing on the cellular responses independent of the systemic regulatory mechanisms. Isolated cells from the skeletal muscle (masseter) and brain (thalamus) were exposed to acute short-term (15 min) hypoxia followed by reoxygenation. The mitochondrial oxygen consumption, reactive oxygen species (ROS) production rates and transcriptional profiles of hypoxia-responsive mRNA and miRNA were determined. Mitochondria of the porcine brain cells showed a decrease in the resting respiration and ATP synthesis capacity whereas the mitochondria from the muscle cells showed robust respiration and less susceptibility to H-R stress. ROS production was not affected by the short-term H-R stress in the brain or muscle cells. Transcriptionally, prolyl hydroxylase domain protein EGLN3 was upregulated during hypoxia and suppressed during reoxygenation in porcine muscle cells. The decline in EGLN3 mRNA during reoxygenation was accompanied by an upregulation of hypoxia-inducible factor subunit α (HIF1A) transcripts in the muscle cells. However, in the brain cells, HIF1A mRNA levels were suppressed during reoxygenation. Other functionally important transcripts and miRNAs involved in antioxidant response, apoptosis, inflammation, and substrate oxidation were also differentially expressed between the muscle and brain cells. Suppression of miRNA levels during acute intermittent hypoxia was stronger in the brain cells affecting ~ 55% of all studied miRNA transcripts than in the muscle cells (~ 25% of miRNA) signifying transcriptional derepression of the respective mRNA targets. Our study provides insights into the potential molecular and physiological mechanisms contributing to different hypoxia sensitivity of the studied tissues and can serve as a starting point to better understand the biological processes associated with hypoxia stress, e.g. during ischemia and reperfusion.

miRNA-mRNA associations with inosine monophosphate specific deposition in the muscle of Jingyuan chicken

1. Inosine monophosphate (IMP), is an essential component for meat flavour and microRNAs (miRNAs) play a vital role in its post-transcriptional regulation. However, the mechanism of how miRNA expression affects muscle-specific IMP deposition is unclear.2. The following study performed transcriptome sequencing and bioinformatics analysis of breast and leg muscle, which have significantly different IMP content in Jingyuan chicken. The differential miRNA-mRNAs were screened out and correlation analysis with IMP content was performed.3. A total of 39 differentially expressed miRNAs (DE miRNAs) and 666 differentially expressed mRNAs (DE mRNAs) were identified between breast muscles and leg muscles. Using miRNA-mRNA integrated analysis, 29 miRNA-target gene pairs were obtained, composed of 13 DE miRNAs and 28 DE mRNAs. Next, purine metabolism, glycolysis/gluconeogenesis, pyruvate metabolism and the biosynthesis of amino acid pathways as necessary for muscle IMP-specific deposition were identified. The differentially expressed gene PKM2, which was significantly enriched in all four pathways, is involved in IMP anabolism in the form of energy metabolism and enzyme activity regulation. The correlation analysis suggested that the gga-miR-107-3p-KLHDC2 negative interaction may be a key regulator in IMP deposition.4. This study explores the functional mechanism of IMP-specific deposition in Jingyuan chicken muscles at the miRNA and mRNA levels and highlights multiple candidate miRNAs and mRNAs for molecular-assisted breeding.

Weighted Co-Expression Network Analysis Identifies RNF181 as a Causal Gene of Coronary Artery Disease

Background: Coronary artery disease (CAD) exerts a global challenge to public health. Genetic heritability is one of the most vital contributing factors in the pathophysiology of CAD. Co-expression network analysis is an applicable and robust method for the interpretation of biological interaction from microarray data. Previous CAD studies have focused on peripheral blood samples since the processes of CAD may vary from tissue to blood. It is therefore necessary to find biomarkers for CAD in heart tissues; their association also requires further illustration.

Materials and Methods: To filter for causal genes, an analysis of microarray expression profiles, GSE12504 and GSE22253, was performed with weighted gene co-expression network analysis (WGCNA). Co-expression modules were constructed after batch effect removal and data normalization. The results showed that 7 co-expression modules with 8,525 genes and 1,210 differentially expressed genes (DEGs) were identified. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted. Four major pathways in CAD tissue and hub genes were addressed in the Hybrid Mouse Diversity Panel (HMDP) and Human Protein Atlas (HPA), and isoproterenol (ISO)/doxycycline (DOX)-induced heart toxicity models were used to validate the hub genes. Lastly, the hub genes and risk variants were verified in the CAD cohort and in genome-wide association studies (GWAS).

Results: The results showed that RNF181 and eight other hub genes are perturbed during CAD in heart tissues. Additionally, the expression of RNF181 was validated using RT-PCR and immunohistochemistry (IHC) staining in two cardiotoxicity mouse models. The association was further verified in the CAD patient cohort and in GWAS.

Conclusion: Our findings illustrated for the first time that the E3 ubiquitination ligase protein RNF181 may serve as a potential biomarker in CAD, but further in vivo validation is warranted.

The intra-class heterogeneity of immunophenotyping and immune landscape in oesophageal cancer and clinical implications

BACKGROUND

The response rate and survival benefit of immunotherapy vary among patients, implying specific immune status of an individual could be associated with the effect of immunotherapy. However, in-depth studies of immune subtypes (ISs), immune landscape and tumour microenvironment of oesophageal cancer (ESCA) and their clinical implications are less reported.

METHODS

We first accessed data from publicly available databases and preprocessed it based on a standard protocol. Then, ISs were identified by unsupervised learning. Thereafter, the association of these ISs and tumour mutation burden (TMB), biomarkers of chemotherapy-induced immune response, tumour markers were also assessed. In addition, the immune characteristics, immune landscape, co-expression network of immune genes, and clinical implications were visualized and analysed.

RESULTS

We identified three immunoclusters based on immune-associated genes with intra-class heterogeneity and prognostic value. Cluster-specific associations with TMB, markers of chemotherapy-induced immune response, and tumour markers were revealed. A 4-gene signature (risk score= -0.16514291×BHLHE22-0.03964046×MXRA8-0.15242778×SLIT2-0.05553572×SPON1) based on co-expressed genes in the immunoclusters was developed and externally validated.

CONCLUSIONS

In summary, we identified clinically relevant immunoclusters in both adenocarcinoma and squamous cell carcinoma of oesophagus, revealing the necessity of assessing the complexity and diversity of immune microenvironment for cancer immunotherapy.

Genetic regulation and heritability of miRNA and mRNA expression link to phosphorus utilization and gut microbiome

Improved utilization of phytates and mineral phosphorus (P) in monogastric animals contributes significantly to preserving the finite resource of mineral P and mitigating environmental pollution. In order to identify pathways and to prioritize candidate genes related to P utilization (PU), the genomic heritability of 77 and 80 trait-dependent expressed miRNAs and mRNAs in 482 Japanese quail were estimated and eQTL (expression quantitative trait loci) were detected. In total, 104 miR-eQTL (microRNA expression quantitative traits loci) were associated with SNP markers (false discovery rate less than 10%) including 41 eQTL of eight miRNAs. Similarly, 944 mRNA-eQTL were identified at the 5% False discovery rate threshold, with 573 being cis-eQTL of 36 mRNAs. High heritabilities of miRNA and mRNA expression coincide with highly significant eQTL. Integration of phenotypic data with transcriptome and microbiome data of the same animals revealed genetic regulated mRNA and miRNA transcripts (SMAD3, CAV1, ENNPP6, ATP2B4, miR-148a-3p, miR-146b-5p, miR-16-5p, miR-194, miR-215-5p, miR-199-3p, miR-1388a-3p) and microbes (Candidatus Arthromitus, Enterococcus) that are associated with PU. The results reveal novel insights into the role of mRNAs and miRNAs in host gut tissue functions, which are involved in PU and other related traits, in terms of the genetic regulation and inheritance of their expression and in association with microbiota components.

Functional contexts of adipose and gluteal muscle tissue gene co-expression networks in the domestic horse

A gene's response to an environment is tightly bound to the underlying genetic variation present in an individual's genome and varies greatly depending on the tissue it is being expressed in. Gene co-expression networks provide a mechanism to understand and interpret the collective transcriptional responses of genes. Here, we use the Camoco co-expression network framework to characterize the transcriptional landscape of adipose and gluteal muscle tissue in 83 domestic horses (Equus caballus) representing 5 different breeds. In each tissue, gene expression profiles, capturing transcriptional response due to variation across individuals, were used to build two separate, tissue-focused, genotypically-diverse gene co-expression networks. The aim of our study was to identify significantly co-expressed clusters of genes in each tissue, then compare the clusters across networks to quantify the extent that clusters were found in both networks as well as to identify clusters found in a single network. The known and unknown functions for each network were quantified using complementary, supervised and unsupervised approaches. First, supervised ontological enrichment was utilized to quantify biological functions represented by each network. Curated ontologies (GO and KEGG) were used to measure the known biological functions present in each tissue. Overall, a large percentage of terms (40.3% of GO and 41% of KEGG) were co-expressed in at least one tissue. Many terms were co-expressed in both tissues, however a small proportion of terms exhibited single tissue co-expression suggesting functional differentiation based on curated, functional annotation. To complement this, an unsupervised approach not relying on ontologies was employed. Strongly co-expressed sets of genes defined by Markov clustering identified sets of unannotated genes showing similar patterns of co-expression within a tissue. We compared gene sets across tissues and identified clusters of genes the either segregate in co-expression by tissue or exhibit high levels of co-expression in both tissues. Clusters were also integrated with GO and KEGG ontologies to identify gene sets containing previously curated annotations versus unannotated gene sets indicating potentially novel biological function. Coupling together these transcriptional datasets, we mapped the transcriptional landscape of muscle and adipose setting up a generalizable framework for interpreting gene function for additional tissues in the horse and other species.

Mitonuclear Interactions in the Maintenance of Mitochondrial Integrity

In eukaryotic cells, mitochondria originated in an α-proteobacterial endosymbiont. Although these organelles harbor their own genome, the large majority of genes, originally encoded in the endosymbiont, were either lost or transferred to the nucleus. As a consequence, mitochondria have become semi-autonomous and most of their processes require the import of nuclear-encoded components to be functional. Therefore, the mitochondrial-specific translation has evolved to be coordinated by mitonuclear interactions to respond to the energetic demands of the cell, acquiring unique and mosaic features. However, mitochondrial-DNA-encoded genes are essential for the assembly of the respiratory chain complexes. Impaired mitochondrial function due to oxidative damage and mutations has been associated with numerous human pathologies, the aging process, and cancer. In this review, we highlight the unique features of mitochondrial protein synthesis and provide a comprehensive insight into the mitonuclear crosstalk and its co-evolution, as well as the vulnerabilities of the animal mitochondrial genome.

Identification of Four Potential Biomarkers Associated With Coronary Artery Disease in Non-diabetic Patients by Gene Co-expression Network Analysis

Background

Coronary artery disease (CAD) is a type of cardiovascular disease that greatly hurts the health of human beings. Diabetic status is one of the largest clinical factors affecting CAD-associated gene expression changes. Most of the studies focus on diabetic patients, whereas few have been done for non-diabetic patients. Since the pathophysiological processes may vary among these patients, we cannot simply follow the standard based on the data from diabetic patients. Therefore, the prognostic and predictive diagnostic biomarkers for CAD in non-diabetic patient need to be fully recognized.

Materials and Methods

To screen out candidate genes associated with CAD in non-diabetic patients, weighted gene co-expression network analysis (WGCNA) was constructed to conduct an analysis of microarray expression profiling in patients with CAD. First, the microarray data GSE20680 and GSE20681 were downloaded from NCBI. We constructed co-expression modules via WGCNA after excluding the diabetic patients. As a result, 18 co-expression modules were screened out, including 1,225 differentially expressed genes (DEGs) that were obtained from 152 patients (luminal stenosis ≥50% in at least one major vessel) and 170 patients (stenosis of <50%). Subsequently, a Pearson's correlation analysis was conducted between the modules and clinical traits. Then, a functional enrichment analysis was conducted, and we used gene network analysis to reveal hub genes. Last, we validated the hub genes with peripheral blood samples in an independent patient cohort using RT-qPCR.

Results

The results showed that the midnight blue module and the yellow module played vital roles in the pathogenesis of CAD in non-diabetic patients. Additionally, CD40, F11R, TNRC18, and calcium/calmodulin-dependent protein kinase type II gamma (CAMK2G) were screened out and validated using enzyme-linked immunosorbent assay (ELISA) in an independent patient cohort and immunohistochemical (IHC) staining in an atherosclerosis mouse model.

Conclusion

Our findings demonstrate that hub genes, CD40, F11R, TNRC18, and CAMK2G, are surrogate diagnostic biomarkers and/or therapeutic targets for CAD in non-diabetic patients and require deeper validation.

Weighted gene correlation network analysis reveals novel biomarkers associated with mesenchymal stromal cell differentiation in early phase

Osteoporosis is a major public health problem that is associated with high morbidity and mortality, and its prevalence is increasing as the world’s population ages. Therefore, understanding the molecular basis of the disease is becoming a high priority. In this regard, studies have shown that an imbalance in adipogenic and osteogenic differentiation of bone marrow mesenchymal stem cells (MSCs) is associated with osteoporosis. In this study, we conducted a Weighted Gene Co-Expression Network Analysis to identify gene modules associated with the differentiation of bone marrow MSCs. Gene Ontology and Kyoto Encyclopedia of Genes and Genome enrichment analysis showed that the most significant module, the brown module, was enriched with genes involved in cell cycle regulation, which is in line with the initial results published using these data. In addition, the Cytoscape platform was used to identify important hub genes and lncRNAs correlated with the gene modules. Furthermore, differential gene expression analysis identified 157 and 40 genes that were upregulated and downregulated, respectively, after 3 h of MSCs differentiation. Interestingly, regulatory network analysis, and comparison of the differentially expressed genes with those in the brown module identified potential novel biomarker genes, including two transcription factors (ZNF740, FOS) and two hub genes (FOXQ1, SGK1), which were further validated for differential expression in another data set of differentiation of MSCs. Finally, Gene Set Enrichment Analysis suggested that the two most important candidate hub genes are involved in regulatory pathways, such as the JAK-STAT and RAS signaling pathways. In summary, we have revealed new molecular mechanisms of MSCs differentiation and identified novel genes that could be used as potential therapeutic targets for the treatment of osteoporosis.

The long noncoding RNA MyHC IIA/X-AS contributes to skeletal muscle myogenesis and maintains the fast fiber phenotype

Mammalian skeletal muscles comprise different types of muscle fibers, and this muscle fiber heterogeneity is generally characterized by the expression of myosin heavy chain (MyHC) isoforms. A switch in MyHC expression leads to muscle fiber-type transition under various physiological and pathological conditions, but the underlying regulator coordinating the switch of MyHC expression remains largely unknown. Experiments reported in this study revealed the presence of a skeletal muscle-specific antisense transcript generated from the intergenic region between porcine MyHC IIa and IIx and is referred to here as MyHC IIA/X-AS. We found that MyHC IIA/X-AS is identified as a long noncoding RNA (lncRNA) that is strictly expressed in skeletal muscles and is predominantly distributed in the cytoplasm. Genetic analysis disclosed that MyHC IIA/X-AS stimulates cell cycle exit of skeletal satellite cells and their fusion into myotubes. Moreover, we observed that MyHC IIA/X-AS is more enriched in fast-twitch muscle and represses slow-type gene expression and thereby maintains the fast phenotype. Furthermore, we found that MyHC IIA/X-AS acts as a competing endogenous RNA that sponges microRNA-130b (miR-130b) and thereby maintains MyHC IIx expression and the fast fiber type. We also noted that miR-130b was proved to down-regulate MyHC IIx by directly targeting its 3'-UTR. Together, the results of our study uncovered a novel pathway, which revealed that lncRNA derived from the skeletal MyHC cluster could modulate local MyHC expression in trans, highlighting the role of lncRNAs in muscle fiber-type switching.

Detection of Co-expressed Pathway Modules Associated With Mineral Concentration and Meat Quality in Nelore Cattle

Meat quality is a complex trait that is influenced by genetic and environmental factors, which includes mineral concentration. However, the association between mineral concentration and meat quality, and the specific molecular pathways underlying this association, are not well explored. We therefore analyzed gene expression as measured with RNA-seq in Longissimus thoracis muscle of 194 Nelore steers for association with three meat quality traits (intramuscular fat, meat pH, and tenderness) and the concentration of 13 minerals (Ca, Cr, Co, Cu, Fe, K, Mg, Mn, Na, P, S, Se, and Zn). We identified seven sets of co-expressed genes (modules) associated with at least two traits, which indicates that common pathways influence these traits. From pathway analysis of module hub genes, we further found an over-representation for energy and protein metabolism (AMPK and mTOR signaling pathways) in addition to muscle growth, and protein turnover pathways. Among the identified hub genes FASN, ELOV5, and PDE3B are involved with lipid metabolism and were affected by previously identified eQTLs associated to fat deposition. The reported hub genes and over-represented pathways provide evidence of interplay among gene expression, mineral concentration, and meat quality traits. Future studies investigating the effect of different levels of mineral supplementation in the gene expression and meat quality traits could help us to elucidate the regulatory mechanism by which the genes/pathways are affected.

Gene co-expression networks associated with carcass traits reveal new pathways for muscle and fat deposition in Nelore cattle

BACKGROUND

Positively correlated with carcass weight and animal growth, the ribeye area (REA) and the backfat thickness (BFT) are economic important carcass traits, which impact directly on producer's payment. The selection of these traits has not been satisfactory since they are expressed later in the animal's life and multigene regulated. So, next-generation technologies have been applied in this area to improve animal's selection and better understand the molecular mechanisms involved in the development of these traits. Correlation network analysis, performed by tools like WGCNA (Weighted Correlation Network Analysis), has been used to explore gene-gene interactions and gene-phenotype correlations. Thus, this study aimed to identify putative candidate genes and metabolic pathways that regulate REA and BFT by constructing a gene co-expression network using WGCNA and RNA sequencing data, to better understand genetic and molecular variations behind these complex traits in Nelore cattle.

RESULTS

The gene co-expression network analysis, using WGCNA, were built using RNA-sequencing data normalized by transcript per million (TPM) from 43 Nelore steers. Forty-six gene clusters were constructed, between them, three were positively correlated (p-value< 0.1) to the BFT (Green Yellow, Ivory, and Light Yellow modules) and, one cluster was negatively correlated (p-value< 0.1) with REA (Salmon module). The enrichment analysis performed by DAVID and WebGestalt (FDR 5%) identified eight Gene Ontology (GO) terms and three KEGG pathways in the Green Yellow module, mostly associated with immune response and inflammatory mechanisms. The enrichment of the Salmon module demonstrated 19 GO terms and 21 KEGG pathways, related to muscle energy metabolism, lipid metabolism, muscle degradation, and oxidative stress diseases. The Ivory and Light yellow modules have not shown significant results in the enrichment analysis.

CONCLUSION

With this study, we verified that inflammation and immune response pathways modulate the BFT trait. Energy and lipid metabolism pathways, highlighting fatty acid metabolism, were the central pathways associated with REA. Some genes, as RSAD2, EIF2AK2, ACAT1, and ACSL1 were considered as putative candidate related to these traits. Altogether these results allow us to a better comprehension of the molecular mechanisms that lead to muscle and fat deposition in bovine.

An integrative transcriptome analysis indicates regulatory mRNA-miRNA networks for residual feed intake in Nelore cattle

Residual Feed Intake (RFI) is an economically relevant trait in beef cattle. Among the molecular regulatory mechanisms, microRNAs (miRNAs) are an important dimension in post-transcriptional regulation and have been associated with different biological pathways. Here, we performed differential miRNAs expression and weighted gene co-expression network analyses (WGCNA) to better understand the complex interactions between miRNAs and mRNAs expressed in bovine skeletal muscle and liver. MiRNA and mRNA expression data were obtained from Nelore steers that were genetically divergent for RFI (N = 10 [low RFI or feed efficient]; N = 10 [high RFI or feed inefficient]). Differentially expressed and hub miRNAs such as bta-miR-486, bta-miR-7, bta-miR15a, bta-miR-21, bta-miR 29, bta- miR-30b, bta-miR-106b, bta-miR-199a-3p, bta-miR-204, and bta-miR 296 may have a potential role in variation of RFI. Functional enrichment analysis of differentially expressed (DE) miRNA's target genes and miRNA-mRNA correlated modules revealed that insulin, lipid, immune system, oxidative stress and muscle development signaling pathways might potentially be involved in RFI in this population. Our study identified DE miRNAs, miRNA - mRNA regulatory networks and hub miRNAs related to RFI. These findings suggest a possible role of miRNAs in regulation of RFI, providing new insights into the potential molecular mechanisms that control feed efficiency in Nelore cattle.

Gene Co-expression Analysis Indicates Potential Pathways and Regulators of Beef Tenderness in Nellore Cattle

Beef tenderness, a complex trait affected by many factors, is economically important to beef quality, industry, and consumer’s palatability. In this study, RNA-Seq was used in network analysis to better understand the biological processes that lead to differences in beef tenderness. Skeletal muscle transcriptional profiles from 24 Nellore steers, selected by extreme estimated breeding values (EBVs) for shear force after 14 days of aging, were analyzed and 22 differentially expressed transcripts were identified. Among these were genes encoding ribosomal proteins, glutathione transporter ATP-binding cassette, sub-family C (CFTR/MRP), member 4 (ABCC4), and synaptotagmin IV (SYT4). Complementary co-expression analyses using Partial Correlation with Information Theory (PCIT), Phenotypic Impact Factor (PIF) and the Regulatory Impact Factor (RIF) methods identified candidate regulators and related pathways. The PCIT analysis identified ubiquitin specific peptidase 2 (USP2), growth factor receptor-bound protein 10 (GBR10), anoctamin 1 (ANO1), and transmembrane BAX inhibitor motif containing 4 (TMBIM4) as the most differentially hubbed (DH) transcripts. The transcripts that had a significant correlation with USP2, GBR10, ANO1, and TMBIM4 enriched for proteasome KEGG pathway. RIF analysis identified microRNAs as candidate regulators of variation in tenderness, including bta-mir-133a-2 and bta-mir-22. Both microRNAs have target genes present in the calcium signaling pathway and apoptosis. PIF analysis identified myoglobin (MB), enolase 3 (ENO3), and carbonic anhydrase 3 (CA3) as potentially having fundamental roles in tenderness. Pathways identified in our study impacted in beef tenderness included: calcium signaling, apoptosis, and proteolysis. These findings underscore some of the complex molecular mechanisms that control beef tenderness in Nellore cattle.

Co-Expression Network Analysis Identifies miRNA⁻mRNA Networks Potentially Regulating Milk Traits and Blood Metabolites

MicroRNAs (miRNA) regulate mRNA networks to coordinate cellular functions. In this study, we constructed gene co-expression networks to detect miRNA modules (clusters of miRNAs with similar expression patterns) and miRNA–mRNA pairs associated with blood (triacylglyceride and nonesterified fatty acids) and milk (milk yield, fat, protein, and lactose) components and milk fatty acid traits following dietary supplementation of cows’ diets with 5% linseed oil (LSO) (n = 6 cows) or 5% safflower oil (SFO) (n = 6 cows) for 28 days. Using miRNA transcriptome data from mammary tissues of cows for co-expression network analysis, we identified three consensus modules: blue, brown, and turquoise, composed of 70, 34, and 86 miRNA members, respectively. The hub miRNAs (miRNAs with the most connections with other miRNAs) were miR-30d, miR-484 and miR-16b for blue, brown, and turquoise modules, respectively. Cell cycle arrest, and p53 signaling and transforming growth factor–beta (TGF-β) signaling pathways were the common gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched for target genes of the three modules. Protein percent (p = 0.03) correlated with the turquoise module in LSO treatment while protein yield (p = 0.003) and milk yield (p = 7 × 10⁻⁰⁴) correlated with the turquoise model, protein and milk yields and lactose percent (p < 0.05) correlated with the blue module and fat percent (p = 0.04) correlated with the brown module in SFO treatment. Several fatty acids correlated (p < 0.05) with the blue (CLA:9,11) and brown (C4:0, C12:0, C22:0, C18:1n9c and CLA:10,12) modules in LSO treatment and with the turquoise (C14:0, C18:3n3 and CLA:9,11), blue (C14:0 and C23:0) and brown (C6:0, C16:0, C22:0, C22:6n3 and CLA:10,12) modules in SFO treatment. Correlation of miRNA and mRNA data from the same animals identified the following miRNA–mRNA pairs: miR-183/RHBDD2 (p = 0.003), miR-484/EIF1AD (p = 0.011) and miR-130a/SBSPON (p = 0.004) with lowest p-values for the blue, brown, and turquoise modules, respectively. Milk yield, protein yield, and protein percentage correlated (p < 0.05) with 28, 31 and 5 miRNA–mRNA pairs, respectively. Our results suggest that, the blue, brown, and turquoise modules miRNAs, hub miRNAs, miRNA–mRNA networks, cell cycle arrest GO term, p53 signaling and TGF-β signaling pathways have considerable influence on milk and blood phenotypes following dietary supplementation of dairy cows’ diets with 5% LSO or 5% SFO.

miRNAs regulate acute transcriptional changes in broiler embryos in response to modification of incubation temperature

MicroRNAs are post-transcriptional regulators that play critical roles in diverse biological processes. We hypothesize that miRNAs may be involved in regulating transcriptome responses to changes in embryonic incubation temperature in chickens affecting differentiation and proliferation processes during tissue development. Therefore, we conducted comparative transcriptome profiling of miRNAs to examine altered expression in breast and hind muscle of embryos and day 35 chickens experiencing high (38.8 °C), control (37.8 °C), or low (36.8 °C) embryonic incubation temperature during embryonic day (ED) 7–10 or ED10–13. The results revealed differential expression of miRNAs due to modification of embryonic incubation temperature in a muscle type-specific and a developmental stage-specific manner. The immediate effects of thermal change observed in embryos were substantial compared to the subtle long-term effects in chickens at day 35 post-hatch. Upregulation of miR-133 in breast muscle and downregulation of miR-199a-5p, miR-1915, and miR-638 in hind muscle post ED7–10 high-temperature treatment are functionally associated with myogenesis and body size. ED10–13 low-temperature treatment led to downregulation of let-7, miR-93, and miR-130c that are related to proliferation and differentiation. The results provide insight into the dynamics of miRNA expression at variable embryonic incubation temperatures during developmental processes and indicate a major regulatory role of miRNAs in acute responses to modified environmental conditions that affect remodelling of cells and tissues.

Genetic architecture and regulatory impact on hepatic microRNA expression linked to immune and metabolic traits

Regulation of microRNA (miRNA) expression contributes to a wide range of target gene expression and phenotypes. The miRNA expression in the liver, the central metabolic organ, was examined in 209 pigs, and integrated with haematological and clinical biomarkers of metabolic and overall health, mRNA-target expression levels and single-nucleotide polymorphism (SNP) genotypes. The expression levels of 426 miRNA species correlated with plasma haematological or biochemical traits (r² = |0.19–0.45|, false discovery rate < 5%). Pairs of these miRNAs and their predicted target mRNAs showing expressing levels associated with the identical traits were examined to understand how immune and metabolic traits are affected by miRNA-mediated regulatory networks derived by mapping miRNA abundance as an expression quantitative trait. In total, 221 miRNA-expression-QTL correspond to 164 SNPs and 108 miRNAs, including miR-34a, miR-30e, miR-148-3p, miR-204, miR-181-5p, miR-143-5p and let-7 g that also correlate with the biomarkers. Sixty-one SNPs were simultaneously associated with 29 miRNA and 41 mRNA species. The expression levels of 13 out of 29 miRNA were correlated with one of the biochemical or haematological traits. For example, the expression levels of miR-34a were correlated with serum phosphorus and cholesterin levels; miR-204, miR-15a and miR-16b were correlated with triglyceride. For haematological traits, the expression levels of miR-652 and miR-204 were correlated with the mean corpuscular haemoglobin concentration, and the expression of miR-143 was correlated with plateletcrit. Pleiotropic association analyses revealed genetic links between mRNA and miRNA on SSC6 for miR-34a, SSC9 for miR-708 and SSC14 for miR-652. Our analysis of miRNA and mRNA transcript profiles, their correlation with clinically important plasma parameters of hepatic functions as well as information on their genetic regulation provide novel regulatory networks and potential new biomarkers for immune and metabolic traits.

Integrative analysis of microRNAs and mRNAs revealed regulation of composition and metabolism in Nelore cattle

BACKGROUND

The amount of intramuscular fat can influence the sensory characteristics and nutritional value of beef, thus the selection of animals with adequate fat deposition is important to the consumer. There is growing knowledge about the genes and pathways that control the biological processes involved in fat deposition in muscle. MicroRNAs (miRNAs) belong to a well-conserved class of non-coding small RNAs that modulate gene expression across a range of biological functions in animal development and physiology. The aim of this study was to identify differentially expressed (DE) miRNAs, regulatory candidate genes and co-expression networks related to intramuscular fat (IMF) deposition. To achieve this, we used mRNA and miRNA expression data from the Longissimus dorsi muscle of 30 Nelore steers with high (H) and low (L) genomic estimated breeding values (GEBV) for IMF deposition.

RESULTS

Differential miRNA expression analysis between animals with extreme GEBV values for IMF identified six DE miRNAs (FDR 10%). Functional annotation of the target genes for these microRNAs indicated that the PPARs signaling pathway is involved with IMF deposition. Candidate regulatory genes such as SDHAF4, FBXO17, ALDOA and PKM were identified by partial correlation with information theory (PCIT), phenotypic impact factor (PIF) and regulatory impact factor (RIF) co-expression approaches from integrated miRNA-mRNA expression data. Two DE miRNAs (FDR 10%), bta-miR-143 and bta-miR-146b, which were upregulated in the Low IMF group, were correlated with regulatory candidate genes, which were functionally enriched for fatty acid oxidation GO terms. Co-expression patterns obtained by weighted correlation network analysis (WGCNA), which showed possible interaction and regulation between mRNAs and miRNAs, identified several modules related to immune system function, protein metabolism, energy metabolism and glucose catabolism according to in silico analysis performed herein.

CONCLUSION

In this study, several genes and miRNAs were identified as candidate regulators of IMF by analyzing DE miRNAs using two different miRNA-mRNA co-expression network methods. This study contributes to the understanding of potential regulatory mechanisms of gene signaling networks involved in fat deposition processes measured in muscle. Glucose metabolism and inflammation processes were the main pathways found in silico to influence intramuscular fat deposition in beef cattle in the integrative mRNA-miRNA co-expression analysis.

Mitochondrial-nuclear crosstalk, haplotype and copy number variation distinct in muscle fiber type, mitochondrial respiratory and metabolic enzyme activities

Genes expressed in mitochondria work in concert with those expressed in the nucleus to mediate oxidative phosphorylation (OXPHOS), a process that is relevant for muscle metabolism and meat quality. Mitochondrial genome activity can be efficiently studied and compared in Duroc and Pietrain pigs, which harbor different mitochondrial haplotypes and distinct muscle fiber types, mitochondrial respiratory activities, and fat content. Pietrain pigs homozygous-positive for malignant hyperthermia susceptibility (PiPP) carried only haplotype 8 and showed the lowest absolute mtDNA copy number accompanied by a decrease transcript abundance of mitochondrial-encoded subunits ND1, ND6, and ATP6 and nuclear-encoded subunits NDUFA11 and NDUFB8. In contrast, we found that haplotype 4 of Duroc pigs had significantly higher mitochondrial DNA (mtDNA) copy numbers and an increase transcript abundance of mitochondrial-encoded subunits ND1, ND6, and ATP6. These results suggest that the variation in mitochondrial and nuclear genetic background among these animals has an effect on mitochondrial content and OXPHOS system subunit expression. We observed the co-expression pattern of mitochondrial and nuclear encoded OXPHOS subunits suggesting that the mitochondrial-nuclear crosstalk functionally involves in muscle metabolism. The findings provide valuable information for understanding muscle biology processes and energy metabolism, and may direct use for breeding strategies to improve meat quality and animal health.

Transcriptomic profiling in muscle and adipose tissue identifies genes related to growth and lipid deposition

Growth performance and meat quality are important traits for the pig industry and consumers. Adipose tissue is the main site at which fat storage and fatty acid synthesis occur. Therefore, we combined high-throughput transcriptomic sequencing in adipose and muscle tissues with the quantification of corresponding phenotypic features using seven Chinese indigenous pig breeds and one Western commercial breed (Yorkshire). We obtained data on 101 phenotypic traits, from which principal component analysis distinguished two groups: one associated with the Chinese breeds and one with Yorkshire. The numbers of differentially expressed genes between all Chinese breeds and Yorkshire were shown to be 673 and 1056 in adipose and muscle tissues, respectively. Functional enrichment analysis revealed that these genes are associated with biological functions and canonical pathways related to oxidoreductase activity, immune response, and metabolic process. Weighted gene coexpression network analysis found more coexpression modules significantly correlated with the measured phenotypic traits in adipose than in muscle, indicating that adipose regulates meat and carcass quality. Using the combination of differential expression, QTL information, gene significance, and module hub genes, we identified a large number of candidate genes potentially related to economically important traits in pig, which should help us improve meat production and quality.

Comparative muscle transcriptome associated with carcass traits of Nellore cattle

Background

Commercial cuts yield is an important trait for beef production, which affects the final value of the products, but its direct determination is a challenging procedure to be implemented in practice. The measurement of ribeye area (REA) and backfat thickness (BFT) can be used as indirect measures of meat yield. REA and BFT are important traits studied in beef cattle due to their strong implication in technological (carcass yield) and nutritional characteristics of meat products, like the degree of muscularity and total body fat. Thus, the aim of this work was to study the Longissimus dorsi muscle transcriptome of Nellore cattle, associated with REA and BFT, to find differentially expressed (DE) genes, metabolic pathways, and biological processes that may regulate these traits.

Results

By comparing the gene expression level between groups with extreme genomic estimated breeding values (GEBV), 101 DE genes for REA and 18 for BFT (false discovery rate, FDR 10%) were identified. Functional enrichment analysis for REA identified two KEGG pathways, MAPK (Mitogen-Activated Protein Kinase) signaling pathway and endocytosis pathway, and three biological processes, response to endoplasmic reticulum stress, cellular protein modification process, and macromolecule modification. The MAPK pathway is responsible for fundamental cellular processes, such as growth, differentiation, and hypertrophy. For BFT, 18 biological processes were found to be altered and grouped into 8 clusters of semantically similar terms. The DE genes identified in the biological processes for BFT were ACHE, SRD5A1, RSAD2 and RSPO3. RSAD2 has been previously shown to be associated with lipid droplet content and lipid biosynthesis.

Conclusion

In this study, we identified genes, metabolic pathways, and biological processes, involved in differentiation, proliferation, protein turnover, hypertrophy, as well as adipogenesis and lipid biosynthesis related to REA and BFT. These results enlighten some of the molecular processes involved in muscle and fat deposition, which are economically important carcass traits for beef production.

Electronic supplementary material

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

Differences in Circulating microRNAs between Grazing and Grain-Fed Wagyu Cattle Are Associated with Altered Expression of Intramuscular microRNA, the Potential Target PTEN, and Lipogenic Genes

We aimed to understand the roles of miRNAs in the muscle tissue maturation and those of circulating microRNAs (c-miRNAs) in beef production of Japanese Black (JB) cattle (Wagyu), a breed with genetically background of superior intermuscular fat depot, by comparing different feeding conditions (indoor grain-feeding vs. grazing on pasture). The cattle at 18 months old were assigned to pasture feeding or conventional indoor grain feeding conditions for 5 months. Microarray analysis of c-miRNAs from the plasma extracellular vesicles led to the detection of a total of 202 bovine miRNAs in the plasma, including 15 miRNAs that differed between the feeding conditions. Validation of the microarray results by qPCR showed that the circulating miR-10b level in the grazing cattle was upregulated compared to that of the grain-fed cattle. In contrast, the levels of miR-17-5p, miR-19b, miR-29b, miR-30b-5p, miR-98, miR-142-5p, miR-301a, miR-374b, miR-425-5p, and miR-652 were lower in the grazing cattle than in the grain-fed cattle. Bioinformatic analysis indicated that the predicted target genes of those c-miRNAs were enriched in gene ontology terms associated with blood vessel morphogenesis, plasma membrane, focal adhesion, endocytosis, collagen, ECM-receptor interaction, and phosphorylation. In the grazing cattle, the elevation of miR-10b expression in the plasma was coincident with its elevation in the longissimus lumborum (LL) muscle. Expression of bovine-specific miR-2478, the most plasma-enriched miRNA, tended to be also upregulated in the muscle but not in the plasma. Furthermore, grazing caused the downregulated mRNA expression of predicted miR-10b and/or miR-2478 target genes, such as DNAJB2, PTEN, and SCD1. Thus, the feeding system used for JB cattle affected the c-miRNAs that could be indicators of grain feeding. Among these, miR-10b expression was especially associated with feeding-induced changes and with the expression of the potential target genes responsible for glucose homeostasis and intramuscular fat depot in the LL muscle of JB cattle.

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.

Multi-omic data integration and analysis using systems genomics approaches: methods and applications in animal production, health and welfare

In the past years, there has been a remarkable development of high-throughput omics (HTO) technologies such as genomics, epigenomics, transcriptomics, proteomics and metabolomics across all facets of biology. This has spearheaded the progress of the systems biology era, including applications on animal production and health traits. However, notwithstanding these new HTO technologies, there remains an emerging challenge in data analysis. On the one hand, different HTO technologies judged on their own merit are appropriate for the identification of disease-causing genes, biomarkers for prevention and drug targets for the treatment of diseases and for individualized genomic predictions of performance or disease risks. On the other hand, integration of multi-omic data and joint modelling and analyses are very powerful and accurate to understand the systems biology of healthy and sustainable production of animals. We present an overview of current and emerging HTO technologies each with a focus on their applications in animal and veterinary sciences before introducing an integrative systems genomics framework for analysing and integrating multi-omic data towards improved animal production, health and welfare. We conclude that there are big challenges in multi-omic data integration, modelling and systems-level analyses, particularly with the fast emerging HTO technologies. We highlight existing and emerging systems genomics approaches and discuss how they contribute to our understanding of the biology of complex traits or diseases and holistic improvement of production performance, disease resistance and welfare.

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.

Integrated analysis of miRNA and mRNA paired expression profiling of prenatal skeletal muscle development in three genotype pigs

MicroRNAs (miRNAs) play a vital role in muscle development by binding to messenger RNAs (mRNAs). Based on prenatal skeletal muscle at 33, 65 and 90 days post-coitus (dpc) from Landrace, Tongcheng and Wuzhishan pigs, we carried out integrated analysis of miRNA and mRNA expression profiling. We identified 33, 18 and 67 differentially expressed miRNAs and 290, 91 and 502 mRNA targets in Landrace, Tongcheng and Wuzhishan pigs, respectively. Subsequently, 12 mRNAs and 3 miRNAs differentially expressed were validated using quantitative real-time PCR (qPCR), and 5 predicted miRNA targets were confirmed via dual luciferase reporter or western blot assays. We identified a set of miRNAs and mRNA genes differentially expressed in muscle development. Gene ontology (GO) enrichment analysis suggests that the miRNA targets are primarily involved in muscle contraction, muscle development and negative regulation of cell proliferation. Our data indicated that more mRNAs are regulated by miRNAs at earlier stages than at later stages of muscle development. Landrace and Tongcheng pigs also had longer phases of myoblast proliferation than Wuzhishan pigs. This study will be helpful to further explore miRNA-mRNA interactions in myogenesis and aid to uncover the molecular mechanisms of muscle development and phenotype variance in pigs.

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.

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.

Wnt antagonist, secreted frizzled-related protein 1, is involved in prenatal skeletal muscle development and is a target of miRNA-1/206 in pigs

Background

The Wnt signaling pathway is involved in the control of cell proliferation and differentiation during skeletal muscle development. Secreted frizzled-related proteins (SFRPs), such as SFRP1, function as inhibitors of Wnt signaling. MicroRNA-1/206(miRNA-1/206) is specifically expressed in skeletal muscle and play a critical role in myogenesis. The miRNA-mRNA profiles and bioinformatics study suggested that the SFRP1 gene was potentially regulated by miRNA-1/206 during porcine skeletal muscle development.

Methods

To understand the function of SFRP1 and miRNA-1/206 in swine myogenesis, we first predicted the targets of miRNA-1/206 with the TargetScan and PicTar programs, and analyzed the molecular characterization of the porcine SFRP1 gene. We performed a temporal-spatial expression analysis of SFRP1 mRNA and miRNA-206 in Tongcheng pigs (a Chinese indigenous breed) by quantitative real-time polymerase chain reaction, and conducted the co-expression analyses of SFRP1 and miRNA-1/206. Subsequently, the interaction between SFRP1 and miRNA-1/206 was validated via dual luciferase and Western blot assays.

Results

The bioinformatics analysis predicted SFRP1 to be a target of miRNA-1/206. The expression level of the SFRP1 was highly varied across numerous pig tissues and it was down-regulated during porcine skeletal muscle development. The expression level of the SFRP1 was significantly higher in the embryonic skeletal compared with postnatal skeletal muscle, whereas miR-206 showed the inverse pattern of expression. A significant negative correlation was observed between the expression of miR-1/206 and SFRP1 during porcine skeletal muscle development (p <0.05). Dual luciferase assay and Western-blot results demonstrated that SFRP1 was a target of miR-1/206 in porcine iliac endothelial cells.

Conclusions

Our results indicate that the SFRP1 gene is regulated by miR-1/206 and potentially affects skeletal muscle development. These findings increase understanding of the biological functions and the regulation of the SFRP1 gene in mammals.

Mito-nuclear co-evolution: the positive and negative sides of functional ancient mutations

Most cell functions are carried out by interacting factors, thus underlying the functional importance of genetic interactions between genes, termed epistasis. Epistasis could be under strong selective pressures especially in conditions where the mutation rate of one of the interacting partners notably differs from the other. Accordingly, the order of magnitude higher mitochondrial DNA (mtDNA) mutation rate as compared to the nuclear DNA (nDNA) of all tested animals, should influence systems involving mitochondrial-nuclear (mito-nuclear) interactions. Such is the case of the energy producing oxidative phosphorylation (OXPHOS) and mitochondrial translational machineries which are comprised of factors encoded by both the mtDNA and the nDNA. Additionally, the mitochondrial RNA transcription and mtDNA replication systems are operated by nDNA-encoded proteins that bind mtDNA regulatory elements. As these systems are central to cell life there is strong selection toward mito-nuclear co-evolution to maintain their function. However, it is unclear whether (A) mito-nuclear co-evolution befalls only to retain mitochondrial functions during evolution or, also, (B) serves as an adaptive tool to adjust for the evolving energetic demands as species' complexity increases. As the first step to answer these questions we discuss evidence of both negative and adaptive (positive) selection acting on the mtDNA and nDNA-encoded genes and the effect of both types of selection on mito-nuclear interacting factors. Emphasis is given to the crucial role of recurrent ancient (nodal) mutations in such selective events. We apply this point-of-view to the three available types of mito-nuclear co-evolution: protein-protein (within the OXPHOS system), protein-RNA (mainly within the mitochondrial ribosome), and protein-DNA (at the mitochondrial replication and transcription machineries).

Metabolic and co-expression network-based analyses associated with nitrate response in rice

BACKGROUND

Understanding gene expression and metabolic re-programming that occur in response to limiting nitrogen (N) conditions in crop plants is crucial for the ongoing progress towards the development of varieties with improved nitrogen use efficiency (NUE). To unravel new details on the molecular and metabolic responses to N availability in a major food crop, we conducted analyses on a weighted gene co-expression network and metabolic profile data obtained from leaves and roots of rice plants adapted to sufficient and limiting N as well as after shifting them to limiting (reduction) and sufficient (induction) N conditions.

RESULTS

A gene co-expression network representing clusters of rice genes with similar expression patterns across four nitrogen conditions and two tissue types was generated. The resulting 18 clusters were analyzed for enrichment of significant gene ontology (GO) terms. Four clusters exhibited significant correlation with limiting and reducing nitrate treatments. Among the identified enriched GO terms, those related to nucleoside/nucleotide, purine and ATP binding, defense response, sugar/carbohydrate binding, protein kinase activities, cell-death and cell wall enzymatic activity are enriched. Although a subset of functional categories are more broadly associated with the response of rice organs to limiting N and N reduction, our analyses suggest that N reduction elicits a response distinguishable from that to adaptation to limiting N, particularly in leaves. This observation is further supported by metabolic profiling which shows that several compounds in leaves change proportionally to the nitrate level (i.e. higher in sufficient N vs. limiting N) and respond with even higher levels when the nitrate level is reduced. Notably, these compounds are directly involved in N assimilation, transport, and storage (glutamine, asparagine, glutamate and allantoin) and extend to most amino acids. Based on these data, we hypothesize that plants respond by rapidly mobilizing stored vacuolar nitrate when N deficit is perceived, and that the response likely involves phosphorylation signal cascades and transcriptional regulation.

CONCLUSIONS

The co-expression network analysis and metabolic profiling performed in rice pinpoint the relevance of signal transduction components and regulation of N mobilization in response to limiting N conditions and deepen our understanding of N responses and N use in crops.

A structural approach to understanding the interactions between colour, water-holding capacity and tenderness

The colour, water-holding capacity (WHC) and tenderness of meat are primary determinants of visual and sensory appeal. Although there are many factors which influence these quality traits, the end-results of their influence is often through key changes to the structure of muscle proteins and their spatial arrangement. Water acts as a plasticiser of muscle proteins and water is lost from the myofibrillar lattice structure as a result of protein denaturation and consequent reductions in the muscle fibre volume with increasing cooking temperature. Changes in the myofilament lattice arrangement also impact the light scattering properties and the perceived paleness of the meat. Causes of variation in the quality traits of raw meat do not generally correspond to variations in cooked meat and the differences observed between the raw muscle and cooked or further processed meat are discussed. The review will also identify the gaps in our knowledge and where further investigation would beneficial.

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.