RNA sequencing analysis revealed differentially expressed genes and their functional annotation in porcine longissimus dorsi muscle affected by dietary lysine restriction

The objective of this study was to investigate the effects of dietary lysine restriction on the global gene expression profile of skeletal muscle in growing pigs. Twelve crossbred (Yorkshire × Landrace) barrows (initial BW 22.6 ± 2.04 kg) were randomly assigned to two dietary treatments (LDD: a lysine-deficient diet; LAD: a lysine-adequate diet) according to a completely randomized experiment design (n = 6). After feeding for 8 weeks, skeletal muscle was sampled from the longissimus dorsi of individual pigs. The muscle total RNA was isolated and cDNA libraries were prepared for RNA sequencing (RNA-Seq) analysis. The RNA-Seq data obtained was then analyzed using the CLC Genomics Workbench to identify differentially expressed genes (DEGs). A total of 80 genes (padj ≤ 0.05) were differentially expressed in the longissimus dorsi muscle of the pigs fed LDD vs. LAD, of which 46 genes were downregulated and 34 genes were upregulated. Gene Ontology (GO) analysis of the DEGs (padj ≤ 0.05) for functional annotation identified those GO terms that are mostly associated with the molecular functions of structural molecules and metabolic enzymes (e.g., oxidoreductase and endopeptidase), biological process of acute-phase response, and amino acid metabolism including synthesis and degradation in the extracellular matrix region. Collectively, the results of this study have provided some novel insight regarding the molecular mechanisms of muscle growth that are associated with dietary lysine supply.


Introduction
Lysine, an essential nutrient for swine, is typically the first limiting amino acid (AA) in those common grain-based swine diets.Studies have shown that dietary deficiency in lysine can lead to reduced growth performance in pigs due to decreased body protein synthesis and increased fat deposition (1,2).Therefore, in the swine industry, the pig diets are commonly supplied with the non-protein-bound free lysine (called crystalline lysine) to ensure that pigs receive enough lysine for their optimal growth (3).
At the molecular biological level, lysine has regulatory effects on gene expression, which can alter animal metabolic and signaling pathways associated with protein and lipid metabolism (4).For example, Wang et al. (2) reported that dietary lysine deficiency may activate the ubiquitination pathway to increase muscle protein degradation and up-regulate the expression of genes associated with lipid biosynthesis in finishing pigs.Additionally, Jin et al. (4,5) also found that dietary lysine deficiency inhibited the satellite cell proliferation and reduced protein synthesis via inhibiting the mTOR signaling pathway in the longissimus dorsi muscle of weaning pigs.However, the global gene expression profile or the transcriptional responses to the dietary lysine restriction or deficiency has not been reported in growing pigs.
The transcriptome of a specific tissue under a certain environmental, nutritional, and physiological condition comprises all types of RNAs that include coding and non-coding RNAs (6,7).In recent years, RNA sequencing (RNA-Seq) -using next-generation sequencing technology -has been employed as a powerful tool to perform a comprehensive analysis and quantification of all RNA species expressed in tissues or cells (7)(8)(9).The main objective of this study was to use RNA-Seq to investigate the effects of dietary lysine restriction on the global gene expression profile in the longissimus dorsi muscle of young growing pigs.The functions of those differentially expressed genes (DEGs) were annotated as well.

Animals and dietary treatments
Twelve crossbred (Yorkshire × Landrace) growing barrows (initial BW 22.6 ± 2.04 kg) were purchased from Prestage Farms of Mississippi (West Point, MS), and housed to an environment-controlled swine barn at the Leveck Animal Research Center of Mississippi State University.Pigs were randomly allocated into 12 individual feeding pens and were allowed to acclimatize to the barn environment for 1 week, while a commercial diet for growing pigs was fed (ad libitum).Thereafter, the pigs were assigned to 2 dietary treatments (n = 6) according to a completely randomized experimental design with pen or pig as experimental unit.
A corn-and soybean meal-based diet (a lysine-deficient diet; LDD) was formulated (Table 1) to meet or exceed NRC (10) recommended requirements for various nutrients including crude protein and various essential AAs but not lysine.A control diet (a lysine-adequate diet; LAD) was formulated by adding L-lysine monohydrochloride (the commonly used commercial form of crystalline lysine) to the LDD at a rate of 0.40%.No effort was made to maintain the constant ratios of other dietary essential AAs relative to lysine.The samples of diets were collected a few times during the feeding trial, mixed, subsampled, and submitted to the Essig Animal Nutrition Laboratory at Mississippi State University for proximate analysis to confirm the contents of major nutrients and energy.Also, the AA composition of the diets were analyzed at an analytical laboratory of Ajinomoto Heartland, Inc. (Chicago, IL).The analyzed nutrients and AA composition of the two experimental diets are presented in Table 2.

Animal trial and sample collection
The animal feeding trial lasted for 8 weeks, and the pigs had ad libitum access to the experimental diets and fresh water during the 8-week period.All the pigs, feeders, and waterers were checked 2 to 3 times daily (0600 to 2000 h).At the end of the trial, the pigs were slaughtered in the Meat Science and Muscle Biology Laboratory of Mississippi State University.Skeletal muscle samples (approximately 2 g/pig) were collected from the middle portion of longissimus dorsi (between the 10th and 12th ribs) of each pig, and immediately snap frozen in liquid nitrogen.The frozen muscle samples were then stored at −80°C freezer until the RNA-Seq analysis was started.All the animal-related experimental protocols (e.g., caring, handling, and

mRNA isolation, and cDNA library construction and sequencing
Total RNA was extracted from the frozen muscle samples using TRIzol Reagent (Invitrogen Corporation, Carlsbad, CA) according to the manufacturer's instructions.The RNA samples were further purified by DNase I (Ambion, Austin, TX) treatment to avoid any DNA contamination.Thereafter, the RNA integrity (RIN) of the purified RNA samples was assessed using Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA), and only the samples with RIN values ≥7 were used for cDNA library preparation.The concentrations and purity of the RNA samples were also determined using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE).
Following the quality control analyses, a paired-end (i.e., 250 bp) cDNA libraries were constructed for each sample by using the NEBNext ® Ultra™ Directional RNA Library Prep Kit (New England Biolabs, Inc., Ipswich, MA) following the manufacturer's protocols.To generate paired-end reads of 100 bp in FASTQ format, a multiplex sequencing of the cDNA libraries was carried out on a HiSeq 4000 sequence analyzer (llumina, San Diego, CA) at the Genome Center of DNA Technologies Core in the University of California, Davis, CA.

Bioinformatics analyses and data interpretation 2.4.1. Sequence mapping and differential gene expression
The raw RNA-Seq data were analyzed using the CLC Genomic Workbench (Version 22.0; Qiagen, Germantown, MD) program.The quality of the sequence reads was first checked by the quality control pipeline with the default program parameters, and the high-quality reads were mapped against the recent version of swine reference genome, Sus scrofa 11.1, housed in a National Center for Biotechnology Information (NCBI) genome database. 1Then, the gene expression data were normalized by calculating the fragments per kilobase per million mapped reads (FPKM).The DEG analytical tool of the CLC Genomic Workbench was then used to perform statistical analysis, for which a modified 'Exact Test, ' developed by Robinson and Smyth (11) and incorporated in the Edge R Bioconductor package by Robinson et al. (12), was implemented.The obtained p-values were corrected by Benjamini-Hochberg (B-H) multiple testing procedure.

Gene enrichment and functional analysis of differentially expressed genes
For Gene Ontology (GO) analysis, the web-based functional analysis tool, DAVID (Database for Annotation, Visualization, and Integrated Discovery) bioinformatics program, was used to analyze gene enrichment and functional annotation of the DEG.Briefly, a list of DEG names was uploaded and stored in the centralized list manager panel of DAVID.Since the program is a case insensitive tool for all the accessions (i.e., the IDs), the IDs of DEG were firstly converted to DAVID gene IDs; A total of 78 (80) DAVID IDs were used for further functional analysis while Sus scrofa genome annotation was selected as the background species.
Reports for the functional annotation chart were generated for Molecular Functions (MF), Biological Processes (BP), Cellular Components (CC), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment.The Fisher Exact Test was used with the DAVID gene IDs, and only the results that met specific criteria (such as a maximum probability of 0.1 and a minimum count of 2) were shown in the chart.The EASE Score threshold, which is a modified Fisher Exact value of p, was used for analyzing those enriched genes and pathways.

Mapping and identification of differentially expressed genes
As shown in Table 3, the RNA-Seq analysis generated approximately 822 million raw reads from the twelve longissimus dorsi tissue samples (i.e., 6 pigs in the LAD group and 6 pigs in the LDD group), and these raw sequence data have been submitted to the NCBI Sequence Read Archive (SRA) and are available under the BioProject accession number PRJNA975688.Following the quality control check, approximately 59 million of low quality and ambiguous sequence reads were removed and those resulted clean sequence reads were further used for analysis.The sequence alignment against the pig genome yielded mapping rates from 87.3 to 95.3% of uniquely aligned sequence reads in pairs.The sequence reads aligned in broken pairs and those unmapped were excluded from the further analysis.
The DEGs were identified using the mode of "against control group" in the CLC Genomic Workbench.This mode tests the differences in gene expression level using "Wald Tests" between the experimental group (i.e., the LDD group) and the control group (i.e., the LAD group).A total of 943 DEGs (p ≤ 0.05) were identified which include 101 upregulated genes (Log 2 Fold Change (FC) ≥ 1) and 337 downregulated genes (Log 2 FC ≤ −1) in the longissimus dorsi muscle of the pigs in the LDD group (Figure 1).Next, the p-values were adjusted (padj) using the Benjamini-Hochberg (B-H) method.A total of 80 genes (padj ≤ 0.05) were found to be differentially expressed due to dietary lysine restriction.The gene names, descriptions, Ensembl identification numbers, and fold changes are listed in Supplementary Table S1.

Functional annotation of the differentially expressed genes
To understand the biological functions related to those DEGs, a functional annotation analysis was conducted with the DAVID bioinformatics program.The gene accession conversion tool in DAVID was able to convert 78 out of 80 DEGs (97.5%) to DAVID gene IDs.The results of GO analysis of the DEGs demonstrate that the alteration of gene expression levels by dietary lysine restriction may affect the cellular BP, MF, and KEGG pathways.As shown in Supplementary Table S2, the results of gene functional annotation were categorized in MF, BP, CC, and KEGG pathway and described in GO types and terms.The enrichment analysis for MF shows that 60% (47 out of 78) genes are enriched in various MF related GO terms.There were 61 records of MF GO terms found, of which 78% are matched with swine population background.As shown in Figure 3, the major MF-related GO terms include the structural molecule activity, serine-type endopeptidase inhibitor activity, small molecule binding, heme binding, oxidoreductase activity, extracellular matrix (ECM) structural constituent, heparan sulfate proteoglycan binding, and cysteine-type endopeptidases inhibitor activity.
Significantly enriched in different BP-related GO terms were 57% of the DEGs (Figure 3), with the terms including the acute-phase response, cell adhesion, positive regulation of phagocytosis, negative regulation of endopeptidase activity, heterocycle metabolic process, heme transport, L-serine biosynthetic process, growth plate cartilage development, negative regulation of growth, regulation of peptidyltyrosine phosphorylation, protein polymerization, high-density lipoprotein particle remodeling, and glutamine metabolic process.
The CC-indicating genes (58% of the DEGs) were detected as highly enriched with the following GO terms (Figure 3): blood microparticle, extracellular space, extracellular region, extracellular exosome, extracellular matrix, keratin filament, intermediate filament, and fibrinogen complex.The KEGG pathway analysis of these 78 DAVID IDs revealed that several cell-signaling and metabolic pathways could be affected by dietary lysine restriction (Figure 3), of which the top 5 are responsible for complement and coagulation cascades, chemical carcinogenesis, biosynthesis of amino acids, metabolism of xenobiotics by cytochrome P450, and steroid hormone biosynthesis pathways.

Discussion
Lysine is typically the first limiting essential AA in grain-based swine diets.It has been known that dietary deficiency in lysine will  negatively affect not only pigs' muscle -the largest part of the body mass -protein deposition, but also the overall heath and growth performance (1,2).The growth performance data generated from this study have been reported previously, which showed that the average daily weight gain and the gain to feed ratio were both significantly decreased in the pigs fed LDD than fed LAD, although there was no difference in the average daily feed intake between the two groups of the pigs (3).The reduced weight gain observed in this study must be the consequence of pigs' poor utilization of dietary nutrients and, especially the proteinogenic AA, for body protein biosynthesis and muscle growth, since the skeletal muscle is the largest organ, as well as the largest protein pool, in the body (1, 3).To further explore the molecular mechanism responsible for the reduced growth (mainly the muscle growth) in response to the dietary lysine restriction, we conducted this global gene expression profiling to identify the DEGs that may have diverse molecular functions associated with the growth and development of the skeletal muscle.Through this study, we found that dietary lysine restriction downregulated some transporter-protein coding genes (Table 4), such as transthyretin (TTR) and albumin (ALB).Although the TTR and ALB genes are highly expressed in the liver according to the human gene annotation, Soprano et al. (13) and Wagatsuma et al. (14) reported that these genes were also expressed in animal skeletal muscle.The TTR gene product is a 55-kDa homotetraeric transporter protein for T3 and T4 (both are thyroid hormones) in the blood (15).According to Monk et al. ( 16), TTR deficient mice show a slower transition from suckling to weaning, slower growth, less muscle mass, and impaired bone growth.In addition, some other studies demonstrated that TTR plays a role in promoting the differentiation of muscle cells and the repair of muscle tissue after injury (15, 17, 18).The ALB gene product is another acute phase reactant and a non-specific protein transporter for numerous hormones existing in the blood.It has the highest binding affinity for T3 and T4 followed by TTR, and it determines the transportation time and distribution of thyroid hormones to target tissues.Hypoalbuminemia can be caused by poor nutrition and ongoing inflammation in animals (19).A crosssectional study reported a significant association between low serum albumin and low muscle strength (20).These results suggest that the reduced skeletal muscle protein turnover and altered local muscle metabolism of thyroid hormones may lead to an intracellular diminished T3 availability (21).Therefore, the downregulation of TTR and ALB genes in this study implies that a reduction in the expression of TTR and ALB genes may be a factor in the negative impact of a lysine restricted diet on muscle growth.This could be due to an impact on those proteins responsible for transporting hormones such as thyroid hormones and other signaling molecules.
The expression of several genes coding for the cytochrome P450 family members, such as cytochrome P450 family 1 subfamily A members 1 and 22 (CYP1A1 and CYP1A22), and cytochrome P450 family 2 subfamily E member 1 (CYP2E1), was decreased in pigs fed the lysine-restriction diet (Table 4).It is known that CYP1A1 and CYP2E1 enzymes metabolize polyunsaturated fatty acids to synthesis biologically active intracellular cell-signaling molecules or lipid mediators, such as ω-3 fatty acid derived resolvins and maresins (22).Some of the lipid mediators are well known for their function in maintaining skeletal muscle mass.During nutritional stress (e.g., deficiency), these anti-inflammatory molecules play important roles in reducing skeletal muscle wasting (23).Consequently, the decreased expression of the genes that are associated with polyunsaturated fatty acid metabolism and synthesis of anti-inflammatory molecules, may negatively contribute to skeletal muscle maintenance and homeostasis.
Of the downregulated genes, phenylalanine hydroxylase (PAH) and tyrosine amino transferase (TAT) are associated with AA metabolism.The PAH enzyme catalyzes the conversion of L-phenylalanine to L-tyrosine, which indicates that dietary lysine restriction may suppress the PAH enzyme activity in pigs.It has long been found that a depressed PAH activity is associated with a decreased ability to dispose the phenylalanine load, leading to an impaired tyrosine formation (24).Tyrosine is used not only in protein synthesis but also in the synthesis of some neurotransmitters (e.g., dopamine, epinephrin, and nor-epinephrin) and hormones (e.g., thyroid hormone).In addition, tyrosine acts as a substrate in reaction with oxoglutarate catalyzed by TAT enzyme to produce 4-hydroxyphenylpyruvate and glutamate.Glutamate can act as a precursor for several neurotransmitters (e.g., GABA) (25).Thus, a decrease in the expression of PAH and TAT in response to dietary lysine restriction may indirectly affect the production of some neurotransmitters and hormones associated with animal growth.
As is known, most keratin isoforms are annotated as hair and skin specific in humans, a contamination of our muscle samples by hair or skin tissue is most unlikely due to our rigorous sample preparation protocol.Furthermore, some previous studies had reported the presence of several type I and type II intermediate filament proteins (known as keratins) in developing and striated muscle (28, 29).For example, keratin 19 (K19) and keratin 8 (K8) are prominent keratins expressed in mature striated muscle, where they play a pivotal role in upholding the structural integrity of muscle fibers (30).
There is no previous research, however, that demonstrated a regulatory effect of lysine on keratin gene transcription and translation in vivo.It is of interesting in finding that dietary lysine may play a role in keratin structural activity, which warrants further investigation.
The ECM of skeletal muscle -a complex meshwork consisting of collagens, glycoproteins, proteoglycans, and elastin -maintains the tissue structure and integrity (31).The mRNA abundance of several genes that encode collagens, such as collagen type II alpha 1 chain (COL2A1) and collagen type IX alpha 1 chain (COL9A1), proteoglycans, such as hyaluronan and proteoglycan link protein 1 (HAPLN1) and aggrecan (ACAN), and glycoproteins, such as tenascin N (TNN), was increased in this study.The collagen is the major protein in ECM.It has lower fractional synthetic rate than the myofibrillar proteins, and its concentration in the skeletal muscle is only 15 to 20% of the myofibrillar proteins.The suboptimal nutritional condition, such as AA deficiency, may increase in the collagen fractional synthesis rate while decreasing the syntheses of myofibrillar and sarcoplasmic proteins (32).Therefore, the dietary lysine restriction may be a suboptimal nutritional condition that upregulates the collagen synthesis in pigs.
Cartilage proteoglycan link proteins, on the other hand, are the major non-collagenous components of the ECM cartilage (33).In this study, the expression of aggrecan (ACAN) and hyaluronan and proteoglycan link protein 1 (HAPLN1) genes was increased in response to dietary lysine restriction.ACAN is best known as a major, essential, and defining cartilage component quantitatively.According to Mundlos et al. (34), the gene expression patterns of aggrecan and link protein are identical in different tissues, which agrees with our finding.The regulatory mechanism that promotes such drastic alteration in mRNA expression in response to lysine restriction remain uncertain and will be subject to further studies.
The de novo biosynthesis of serine and glycine can generate carbon units that satisfy many metabolic demands including nucleotide precursors, redox maintenance, and substrates for methylation reaction (35,36).Of the several enzymatic reaction steps in serine biosynthesis pathway, phosphoglycerate dehydrogenase (PHGDH) catalyzes the conversion of 3-phosphoglycerate into 3-phosphohydroxypyruvate.Phosphoserine aminotransferase 1 (PSAT1) is a pivotal enzyme that regulates the production of two metabolites, serine, and α-ketoglutarate (α-KG), which are involved in one carbon metabolism and the TCA cycle, respectively (1,37).This study shows that both PHGDH and PSAT1 genes were upregulated in the longissimus dorsi muscle of the pigs fed the lysine restricted diet.This is in agreement with one of our previous studies, which revealed that dietary lysine deficiency decreased the expression of PHGDH and PSAT1 in the longissimus dorsi muscle of finishing pigs (2).Overexpression of these genes indicate that dietary lysine level could affect the balance of one carbon metabolism and its associated metabolic pathways.
The widely expressed activating transcription factor (ATF5) regulates cell cycle, differentiation, homeostasis, and survival (38,39).The results of this study indicate that the expression of ATF5 was increased in pigs fed the LDD diet compared to the LAD diet.Similar findings in previous studies reported that ATF5 expression was increased in response to cellular stress conditions such as AA limitation, heat stress, or oxidative stress (39,40).Increased activation of ATF5 may upregulate the protein kinases phosphorylation of EIF2α and lead to a global reduction in protein translation.
Six BP terms were significantly enriched under the lysine restriction status (Figure 3), of which protein polymerization (66-fold), fibrinolysis (79-fold), and plasminogen activation (99-fold) were highly enriched.Based on these terms, it can be seen that some DEGs are mostly associated with the homeostasis, wound healing, inflammation, angiogenesis, and several other biological functions (41,42).For example, fibrinogen, a proteolytic glycoprotein, composed of three closely linked polypeptides coded from FGA, FGB, and FGG genes.Dietary lysine restriction downregulated these genes that may affect the rate of fibrinogen synthesis and result in abnormalities in blood circulation (42), which in turn may indirectly affect the transportation of oxygenated blood or nutrients to the muscle.
Secondly, there were 4 MF terms that were significantly enriched (Figure 3), of which, the highly enriched one was endopeptidase inhibitor activity (62-fold).Endopeptidases, including trypsin, chymotrypsin, elastase, and pepsin, are proteolytic enzymes that break the peptide bonds of non-terminal AAs (43,44), and are associated with a wide range of biological roles, such as developmental processes, digestion, fertilization, blood coagulation, apoptosis, fibrinolysis, and immune defense (45).The negative effect of dietary lysine restriction on pig growth performance may be also associated with endopeptidase inhibition.
Thirdly, there were 8 CC terms that were significantly enriched (Figure 3), of which the fibrinogen complex was enriched by 166-fold.The CC terms in GO analysis demonstrate the subcellular location or the structure of macromolecular complex of gene products (46).As shown in Figure 3, a high percentage of gene products among the DEGs are part of fibrinogen complex, intermediate filament, and ECM.Fibrinogen, an important component of ECM, is a large complex protein composed of three polypeptide chains, such as A, B, and G that are encoded by FGA, FGB, and FGG genes, respectively (47,48).The DEG result of this study indicates that the dietary lysine restriction decreased the fibrinogen gene expression (Table 4) that may lead to a reduction in the amount of fibrinogen available to be incorporated into ECM, leading to a weaker ECM and a reduction in the angiogenesis process, which could affect the integrity of ECM and muscle growth (49,50).
Finally, the results of the GO terms of KEGG pathway enrichment indicate that the negative impact of feeding lysine restricted diet to pigs on the growth may be mediated via the signaling and metabolic pathways that include AA biosynthesis, cholesterol metabolism, receptor interaction, and complement and coagulation cascades (Figure 3).For example, Carbamoyl l-phosphate synthase 1 (CPS1) gene encode an enzyme involved in nitrogen metabolism that removes nitrogen waste products from the body (51).The decrease in CPS1 expression was observed in this study suggesting that it may cause disruption in the balance of muscle protein synthesis and degradation, leading to muscle wasting.It is important to note that AA metabolism is a very complex metabolic process whose regulatory mechanism warrants further investigation.
At the end, authors would like to point out that the mRNA abundance does not always correlate with protein abundance, even though many scientists use mRNA abundance to indicate protein abundance because in many cases more mRNA would translate out more protein.Also, it should be kept in mind that the global gene expression data in general are not easy to be interpreted simply because of the complexities not only in the data themselves but also in the aspect of animal biological processes.

Conclusion
A transcriptome analysis using RNA-Seq technology revealed that dietary lysine restriction altered the gene expression profile in longissimus dorsi muscle of growing pigs.We found that 80 genes were differentially expressed, of which 46 genes were down-regulated and 34 genes were up-regulated.Some of the down-regulated genes (e.g., TTR, ALB) are associated with transportation of signaling molecules (e.g., hormones and neurotransmitters) and AA metabolism (e.g., PAH, TAT).On the other hand, the down-regulated genes, such as keratin, ACAN, and COL isoforms are associated with the cell structural molecule activities and AA metabolism (e.g., PHGDH and PSAT1).The GO analysis results suggest that the dietary lysine restriction may influences some biological processes, such as blood coagulation and inflammation, some molecular functions, such as endopeptidase inhibitory activity; and some KEGG pathways, such as the complement and coagulation pathway.Collectively, the results generated from this study have provided some critical novel insight regarding the molecular mechanisms of muscle growth that are associated with dietary lysine supply.

FIGURE 1 Volcano
FIGURE 1Volcano (scatter) plot for significantly up-and down-regulated genes.X-axis and Y-axis denote the Log 2 Fold Change (FC) and −log 10 of adjusted pvalues (padj), respectively; where −1 ≥ log 2 FC ≥ 1.0 and padj < 0.05 were considered significant changes, which are indicated in violet (upregulated) and blue (downregulated) color.

FIGURE 2
FIGURE 2Summary of differentially expressed genes (DEGs) in the skeletal muscle of growing pigs affected by dietary lysine restriction (padj ≤ 0.05).

FIGURE 3
FIGURE 3The significant Gene Ontology (GO) terms enriched by differentially expressed genes (DEGs) in the longissimus dorsi muscle of growing pigs fed a lysine restricted vs. a lysine adequate diet.BP, Biological Processes; MF, Molecular Functions; CC, Cellular Components; KEGG, Kyoto Encyclopedia of Genes and Genomes.

TABLE 1
Composition of the two experimental diets fed to the growing pigs (as-fed basis).

TABLE 2
Analyzed nutrient composition (%, or as indicated) of the two experimental diets fed to the growing pigs (as-fed basis).

TABLE 3
The mapping results of the RNA sequencing (RNA-Seq) reads.

TABLE 4
The top 25 downregulated genes in the longissimus dorsi muscle of the growing pigs fed a lysine-deficient vs. a lysine-adequate diet.
1 log2 Fold Change (FC) indicates the level of expression: the negative (−) values indicate the downregulation.2Adjustedp-values (padj) were obtained by using the Benjamin-Hochberg multiple testing method provided in the CLC Genomic Workbench (Version 11; Qiagen, Germantown, MD).

TABLE 5
The top 25 upregulated genes in the longissimus dorsi muscle of the growing pigs fed a lysine-restricted vs. a lysine-adequate diet.
1 log2 Fold Change (FC) indicates the level of expression: the positive values indicate the upregulation.2