Next Article in Journal
Antimicrobial Action of Water-Soluble β-Chitosan against Clinical Multi-Drug Resistant Bacteria
Previous Article in Journal
Synthesis of Spongy-Like Mesoporous Hydroxyapatite from Raw Waste Eggshells for Enhanced Dissolution of Ibuprofen Loaded via Supercritical CO2
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 16
Issue 4
10.3390/ijms16047976
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

RNA Sequencing Identifies Upregulated Kyphoscoliosis Peptidase and Phosphatidic Acid Signaling Pathways in Muscle Hypertrophy Generated by Transgenic Expression of Myostatin Propeptide

by
Yuanxin Miao
1,2,
Jinzeng Yang
3,
Zhong Xu
1,2,
Lu Jing
1,2,
Shuhong Zhao
1,2 and
Xinyun Li
1,2,*
1
Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
2
The Cooperative Innovation Center for Sustainable Pig Production, Wuhan 430070, China
3
Department of Human Nutrition, Food and Animal Sciences, University of Hawaii at Manoa, Honolulu, HI 96822, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2015, 16(4), 7976-7994; https://doi.org/10.3390/ijms16047976
Submission received: 23 December 2014 / Revised: 18 March 2015 / Accepted: 30 March 2015 / Published: 9 April 2015
(This article belongs to the Section Biochemistry)
Browse Figures
Versions Notes

Abstract

:
Myostatin (MSTN), a member of the transforming growth factor-β superfamily, plays a crucial negative role in muscle growth. MSTN mutations or inhibitions can dramatically increase muscle mass in most mammal species. Previously, we generated a transgenic mouse model of muscle hypertrophy via the transgenic expression of the MSTN N-terminal propeptide cDNA under the control of the skeletal muscle-specific MLC1 promoter. Here, we compare the mRNA profiles between transgenic mice and wild-type littermate controls with a high-throughput RNA sequencing method. The results show that 132 genes were significantly differentially expressed between transgenic mice and wild-type control mice; 97 of these genes were up-regulated, and 35 genes were down-regulated in the skeletal muscle. Several genes that had not been reported to be involved in muscle hypertrophy were identified, including up-regulated myosin binding protein H (mybph), and zinc metallopeptidase STE24 (Zmpste24). In addition, kyphoscoliosis peptidase (Ky), which plays a vital role in muscle growth, was also up-regulated in the transgenic mice. Interestingly, a pathway analysis based on grouping the differentially expressed genes uncovered that cardiomyopathy-related pathways and phosphatidic acid (PA) pathways (Dgki, Dgkz, Plcd4) were up-regulated. Increased PA signaling may increase mTOR signaling, resulting in skeletal muscle growth. The findings of the RNA sequencing analysis help to understand the molecular mechanisms of muscle hypertrophy caused by MSTN inhibition.

Graphical Abstract

1. Introduction

Myostatin, a member of the transforming growth factor-β superfamily and also known as growth factor and differentiation factor-8 (GDF-8), is an essential muscle endogenous negative regulatory factor secreted by skeletal muscles. Myostatin is mainly expressed in skeletal muscle tissue and negatively regulates muscle growth and development [1,2]. In 1997, MSTN was first confirmed as a muscle growth inhibitor in a knockout mouse model. As reported, the muscles of the myostatin-null mice weighed almost two times as much as those of the controls due to increased muscle fiber number and diameter [3,4,5]. Myostatin was highly conserved, and it can inhibit muscle growth in other mammalian species. The mutation and inhibition of MSTN can dramatically increase muscle mass in most mammal species. A natural mutation of myostatin in double-muscled cattle may account for the double muscling phenotype [6,7]. A genetic study confirmed that one G/A mutation identified in the 3' UTR of the myostatin gene formed a novel binding site for miR-1 and miR-206 in Texel sheep. These miRNAs inhibited the protein expression of MSTN, causing the muscle hypertrophy phenomenon in this sheep breed [8]. Furthermore, a mutant MSTN gene causes muscle hypertrophy in humans [9].
MSTN is synthesized as an inactive precursor protein, which is the predominant form in the muscle extracellular matrix and in circulation [10]. The precursor protein undergoes dimerization and two proteolytic cleavage processes to generate a C-terminal dimer—an active biological molecule. The dimer remains in an inactive, latent complex with the N-terminal pro-peptide and circulates in the blood [5]. The C-terminal dimer (mature myostatin) can be activated when separated with the N-terminal, and mature myostatin can bind to the transmembrane ActRIIA/B receptor, which then recruits the type I receptor, ALK-4/5, and phosphorylates the type I receptor to activate its kinase domains [11,12,13]. After binding the receptors, Smad2 and Smad3 can be activated to form a heterodimer complex with Smad4. The Smad complex then activates the Erk1/2 MAPK pathway, which can prevent the proliferation of myoblasts via the p21/Rb signaling cascade pathway and can promote anti-apoptotic pathways by activating p53 in differentiated cells. Moreover, activated Smad3 can repress the level of the muscle regulatory factors (MRFs) [14]. During the interference of Smad3, MyoD sequesters in the cytoplasm and does not enter the nucleus to activate myogenesis [15]. The activated Smad complex can activate Smad4 and cause it to translocate to the nucleus to block the transactivation of MyoD and activate the expression of Smad7 [11,16,17]. Furthermore, phosphorylated SMAD3 can inhibit the PI3K/Akt/mTOR signaling pathway via the induction of E3-ligase atrogin-1 [18,19]. A previous study also found that MSTN can induce muscle protein degradation by activating the ubiquitin-proteasome to inhibit the PI3K/Akt/mTOR signaling pathway [20].
Blocking the myostatin signaling transduction pathway can induce a highly muscled phenotype [21]. The ability of several myostatin antagonists, such as follistatin, follistatin-related gene, soluble activin receptor IIB and suramin to block MSTN signaling transduction has been tested [22,23,24,25]. Furthermore, some microRNAs can affect the MSTN signaling pathway by targeting the 3' UTR of myostatin [26,27]. We previously generated transgenic mice via the muscle-specific expression of a cDNA sequence (5'-region 886 nucleotides) encoding the MSTN propeptide. The transgene effectively depressed myostatin function, and transgenic mice showed dramatic growth and increases in muscle mass due to muscle fiber hypertrophy [28,29,30]. To further elucidate the mechanisms of muscle hypertrophy caused by MSTN inhibition, we compared the gene expression profiles of the muscle tissue between MSTN propeptide transgenic mice and their littermate wild-type controls. Our results indicate that the Ky/mybph/actins pathway is a novel noncanonical signaling pathway, and phosphatidic acid inhibits mTOR signaling via Dgki, Dgkz and Plcd4 in the hypertrophic skeletal muscle of the MSTN propeptide transgenic mice.

2. Results

2.1. Skeletal Muscle Weights of the Transgenic Mice and Wild-Type Littermate Controls

The mice were sacrificed at four months of age to evaluate the muscle weight and obtain muscle samples. The results were consistent with our previous report; the individual major muscles of the myostatin propeptide transgenic mice were significantly heavier (38% to 95% more muscle) than those of the wild-type littermate controls (Figure 1). The muscle groups that increased in mass predominantly consisted of fast-twitch muscle fibers, including the gastrocnemius, tibialis anterior, biceps femoris, triceps brachii, longissimus dorsi, and semitendinosus.
Ijms 16 07976 g001 1024
Figure 1. The skeletal muscle weights of the transgenic mice and wild-type littermate controls.
Figure 1. The skeletal muscle weights of the transgenic mice and wild-type littermate controls.
Ijms 16 07976 g001

2.2. Solexa Sequencing Data and Reads Assembly Using Tophat Software

In this study, Solexa sequencing was used to detect the mRNA expression profiles of MSTN propeptide transgenic mice and their littermates. Five RNA libraries were constructed from the total RNA isolated from the gastrocnemius muscle of the transgenic and control mice, which generated 9,710,772 (CN148), 10,371,693 (CN164), 9,991,395 (TN126), 10,258,517 (TN135), and 10,076,991 (TN329) raw reads. After filtering the low quality reads and trimming off the adapted fragments, 8,717,967, 9,418,459, 9,054,439, 9,301,124, and 9,032,266 clean reads were obtained, and the ratios of cleans reads in raw reads were 89.78%, 90.81%, 90.62%, 90.67%, and 89.63% (Table 1). The detailed description of the sequence data is shown in Table 1, and the abundance of mRNA reads mapped to each chromosome is shown in Figure S1. The Tophat aligner software was used to align the sequencing reads to the reference genome download from the ENSEMBEL website ftp://ftp.ensembl.org/pub/release-75/fasta/mus_musculus/dna/. The mapping rate of each sample exceeded 97% (Table 2). The statistical alignment data indicated that the sequencing reads were of high quality and the sequencing depth was sufficient for a differential expression analysis between the two groups. Subsequently, HTSeq was selected to control the quality of aligned reads. The quality control result is shown in Figure S2.
Table
Table 1. The overview of raw RNA sequencing data in muscle samples of propeptide transgenic and control mice.
Table 1. The overview of raw RNA sequencing data in muscle samples of propeptide transgenic and control mice.
CN148CN164TN126TN135TN329
Raw reads9,710,77210,371,6939,991,39510,258,51710,076,991
Clean reads8,717,9679,418,4599,054,4399,301,1249,032,266
Ratio89.78%90.81%90.62%90.67%89.63%
Table
Table 2. The clean reads counts and alignment results of each library.
Table 2. The clean reads counts and alignment results of each library.
CN148CN164TN126TN135TN329
Input8,717,9679,418,4599,054,4399,301,1249,032,266
Output8,717,5239,418,2479,054,1989,300,8769,031,776
Qualified99.99%99.99%99.99%99.99%99.99%
Mapped8,583,2239,251,3388,865,8329,118,8898,890,164
Mapped%98.5%98.2%97.9%98%98.4%
Multiple alignments2,249,6971,869,4721,801,5331,789,7152,407,543
Multiple alignments%26.2%20.2%20.3%19.6%27.1%
Min read length3030303030
Max read length5050505050

2.3. Differentially Expressed Genes of Skeletal Muscle between Transgenic Mice and Their Littermate Controls

The cufflinks software was used to assemble the individual transcripts that had been aligned to the genome and quantify the expression level of each transfrag in the sample. After cufflinks treated the reads, cuffmerge was used to parsimoniously merge the assembled transfrags. The annotation reference file of Mus musculus integrated into the merged assembly was downloaded from ENSEMBEL (ftp://ftp.ensembl.org/pub/release-75/gtf/mus_musculus/) [31]. Cuffdiff was used to calculate the differential expression and analyze the significances of observed changes between two groups. A cluster analysis was performed to elaborate the expression patterns of genes in MSTN propeptide transgenic mice and their controls (Figure 2A). The results showed that approximately 70% genes were up-regulated in the transgenic mice. CummeRbund volcano plots (Figure 2B) were available to intuitively observe the genes that were differentially expressed between the two groups. Finally, the expression levels of 132 genes were identified to be significantly different between the transgenic mice and control mice; among these genes, 97 were up-regulated and 35 were down-regulated (FDR < 0.05). Of the 132 genes, 117 were annotated in the database. The differentially expressed genes with more than two-fold changes and an FPKM value more than 5 are listed in Table 3. Several genes have not been reported to be involved in skeletal muscle hypertrophy, including myosin binding protein H (mybph), beta, gamma and alpha actin (Actb, actc1, actg1), protein tyrosine phosphatase receptor type c (Ptprc), zinc metallopeptidase STE24 (zmpste24), and genes related to phosphatidic acid (PA) pathways (Dgki, Dgkz and Plcd4). Moreover, the gene kyphoscoliosis peptidase (Ky), which plays an important role in muscle growth, was also up-regulated in the skeletal muscle of the transgenic mice. The bar plots of the differentially expressed genes are shown in Figure S3, and all differentially expressed genes identified by RNA-seq are summarized in Table S1.
Ijms 16 07976 g002 1024
Figure 2. Differentially expressed genes between MSTN propeptide transgenic mice and their littermate control mice. (A) A heatmap of FPKM expression values in the five samples. Each row represents a gene ID and the column represents a sample; (B) The volcano plots reveal genes that differ significantly between two groups. The black dots represent the expression levels of genes that do not differ and the red dots represent the expression level differences.
Figure 2. Differentially expressed genes between MSTN propeptide transgenic mice and their littermate control mice. (A) A heatmap of FPKM expression values in the five samples. Each row represents a gene ID and the column represents a sample; (B) The volcano plots reveal genes that differ significantly between two groups. The black dots represent the expression levels of genes that do not differ and the red dots represent the expression level differences.
Ijms 16 07976 g002
Table
Table 3. Differentially expressed genes between myostatin propeptide transgenic and littermate control mice. (FPKM ≥ 5; fold change ≥ 2). inf represent the fold change cannot be calculated because the expression level was not detectable in wild-type mice.
Table 3. Differentially expressed genes between myostatin propeptide transgenic and littermate control mice. (FPKM ≥ 5; fold change ≥ 2). inf represent the fold change cannot be calculated because the expression level was not detectable in wild-type mice.
Gene IDGeneOfficial Full NameCN ValueTN ValueFold Changep-Value
Genes Related to Cytoskeleton and Muscle Development
XLOC_010673Myh6Myosin heavy polypeptide 65.490.70−7.830.00010
XLOC_000666MybphMyosin binding protein H 14.7757.603.900.00005
XLOC_003363RmstRhabdomyosarcoma 2 associated transcript0.5315.1128.720.00005
XLOC_006460Actg1Actin gamma cytoplasmic 111.26355.4231.560.00005
XLOC_019664Actc1Actin alpha cardiac muscle 150.30257.645.120.00005
XLOC_034484KyKyphoscoliosis peptidase (CD1)3.268.812.700.00015
XLOC_026810ActbActin beta26.2467.112.560.00020
XLOC_000205Mstn-proMyostatin propeptide32.03409.3512.780.00005
Genes Related to Lipoprotein Metabolic Process and Intracellular Signaling Cascade
XLOC_034105Apoa1Apolipoprotein A-I11.041.76 −6.270.00025
XLOC_019500DgkzDiacylglycerol kinase zeta22.1490.104.070.00025
XLOC_000353Plcd4Phospholipase C delta 43.7617.384.630.00010
XLOC_025320AlbAlbumin83.449.80−8.510.00005
XLOC_002012G0s2G0/G1 switch gene 24.3523.625.420.00005
Genes for Immune Response
XLOC_007714IghmImmunoglobulin heavy constant mu27.2473.442.700.00020
XLOC_007114Tnfaip2Tumor necrosis factor alpha-induced protein 23.9017.934.600.00010
XLOC_031749Trim12cTripartite motif-containing 12C0.3434.88103.840.00025
Genes for Oxidation Reduction
XLOC_006482Cbr2Carbonyl reductase 24.4441.629.370.00010
XLOC_010336Mss51Mitochondrial translational activator77.3816.04−4.830.00005
XLOC_020605Kcnab1Potassium voltage-gated channel shaker-related Subfamily beta member 12.8611.063.870.00020
Genes for Cell Proliferation
XLOC_001651PtprcProtein tyrosine phosphatase receptor type C (CD45)0.817.909.770.00005
XLOC_026557Gm13841Predicted gene 13841 (ribosomal protein L29)0.00 8.32 inf0.00005
Genes of Proteolysis
XLOC_024383Zmpste24Zinc metallopeptidase STE245.6727.674.880.00015
Genes of Carbohydrate Phosphatase Activity
XLOC_018711Pfkfb36-phosphofructo-2-kinase/fructose-2,6-biphosphatase 38.8330.523.460.00005
Genes of Integral to Membrane, Intrinsic to Membrane
XLOC_015439Aqp4Aquaporin 42.406.812.830.00020
Genes of Extracellular Region
XLOC_025292Csn1s2acasein alpha s2-like A0.008.54inf0.00005
Others
XLOC_029580Josd2Josephin domain containing 2109.3522.64−4.830.00005
XLOC_025924Gm15459Predicted gene 15459 (heat shock protein 8 pseudogene)7.4621.772.920.00005
XLOC_003234Odf3l2Outer dense fiber of sperm tails 3-like 214.3957.113.970.00005

2.4. Function Annotation of the Differentially Expressed Genes

To gain insight into the function of the differentially expressed genes in the hypertrophic skeletal muscle of the MSTN propeptide transgenic mice, the DAVID (Database for Annotation, Visualization and Integrated Discovery) website (http://david.abcc.ncifcrf.gov/) was used to identify the gene function with the following parameters: Count = 2 and EASE = 0.01. The gene ontology functional classification is shown in Figure 3. The functions of genes in biological processes, cellular components and molecular processes were annotated based on the GO categories. Most genes were enriched in two types of cellular components: the extracellular region and extracellular space. The genes enriched in biological processes were mainly involved in biological adhesion and cell adhesion. The genes enriched in molecular function were mainly associated with calcium ion binding.
Ijms 16 07976 g003 1024
Figure 3. Overview of enriched gene ontology functional classifications of muscle differential expression genes between MSTN propeptide transgenic and wild-type mice.
Figure 3. Overview of enriched gene ontology functional classifications of muscle differential expression genes between MSTN propeptide transgenic and wild-type mice.
Ijms 16 07976 g003

2.5. Validation of the Differentially Expressed Genes by RT-qPCR

To verify the reliability of the sequencing data, five genes were randomly selected from Table 3 for RT-qPCR analyses (Figure 4). As expected, the gene expression patterns validated by RT-qPCR strongly correlated with the sequencing results (Pearson’s r correlation coefficient is 0.97591, Figure 5). Both RT-qPCR and solexa sequencing indicated that MSTN propeptide, Cbr2, Tnfaip2 and Mybph were up-regulated, while Mss51 was down-regulated in MSTN propeptide transgenic mice. Among the differentially expressed genes, the MSTN propeptide mRNA level increased by 12.78-fold according to solexa sequencing (Figure S4) and 43-fold according to the RT-qPCR analyses (Figure 4A). We also detected the expression level of endogenous MSTN; the endogenous myostatin expression did not significantly differ (p > 0.05) between the transgenic mice and their littermate wild-type controls (Figure 4B).
Ijms 16 07976 g004 1024
Figure 4. Q-PCR validation of genes from RNA-seq results between MSTN propeptide transgenic and control mice. All samples were normalized to 18S RNA. (A) Tnfaip2, Mybph, Cbr2 and MSTN propeptide were highly expressed in transgenic mice (TN) compared with control mice (CN). Mss51 had lower expression in transgenic mice (TN) compared with control mice (CN); (B) Endogenous myostatin (C-terminal mature peptide) mRNA levels in transgenic mice was slightly higher than the controls, but not significantly different from the control (p > 0.05). The error bars show the SD. Two-tailed t-test was used to calculate the significance of differentially expressed genes. Mark * in Figure 4A represent the expression level significantly different between two group (p < 0.05).
Figure 4. Q-PCR validation of genes from RNA-seq results between MSTN propeptide transgenic and control mice. All samples were normalized to 18S RNA. (A) Tnfaip2, Mybph, Cbr2 and MSTN propeptide were highly expressed in transgenic mice (TN) compared with control mice (CN). Mss51 had lower expression in transgenic mice (TN) compared with control mice (CN); (B) Endogenous myostatin (C-terminal mature peptide) mRNA levels in transgenic mice was slightly higher than the controls, but not significantly different from the control (p > 0.05). The error bars show the SD. Two-tailed t-test was used to calculate the significance of differentially expressed genes. Mark * in Figure 4A represent the expression level significantly different between two group (p < 0.05).
Ijms 16 07976 g004
Ijms 16 07976 g005 1024
Figure 5. Line fit plot of Q-PCR results and RNA-Seq data for selected genes difference expression between MSTN propeptide transgenic and control mice. Linear regression model and R-Squared shown in the figure.
Figure 5. Line fit plot of Q-PCR results and RNA-Seq data for selected genes difference expression between MSTN propeptide transgenic and control mice. Linear regression model and R-Squared shown in the figure.
Ijms 16 07976 g005

2.6. Signaling Pathway Analysis

To identify pathways that are involved in the regulation of muscle hypertrophy in the MSTN propeptide transgenic mice, we employed the DAVID website to analyze the differentially expressed genes with the following parameters: Count = 2 and EASE = 0.1. The pathways and related genes are listed in Table 4. This analysis identified the hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy signaling pathways, both of which regulate muscle growth and development. In addition, pathways associated with the maintenance of muscle structure and myocytes were identified, including the ECM-receptor interaction and focal adhesion pathways. Notably, the gene expression of three key enzymes, Dgki, Dgkz and Plcd4, of the phosphatidic acid (PA) signaling system was significantly enhanced in the skeletal muscle of the transgenic mice compared with their littermate controls (Figure 6). The hematopoietic cell lineage genes Itga4, cd3d and cd5 were also identified.
Table
Table 4. The pathways regulated by differentially expressed genes.
Table 4. The pathways regulated by differentially expressed genes.
CategoryTermGene CountsGenes
KEGG_PATHWAYHypertrophic cardiomyopathy (HCM)5Actb, Itga4, Actg1, Myh6, Acta
KEGG_PATHWAYDilated cardiomyopathy5Actb, Itga4, Actg1, Myh6, Acta
KEGG_PATHWAYFocal adhesion5Lama1, Lamc3, Actb, Itga4, Actg1
KEGG_PATHWAYArrhythmogenic right ventricular cardiomyopathy (ARVC)3Actb, Itga4, Actg1
KEGG_PATHWAYPhosphatidic acid signaling system3Dgki, Dgkz, Plcd4
KEGG_PATHWAYComplement and coagulation cascades3C4bp, Proc, F12
KEGG_PATHWAYECM-receptor interaction3Lama1, Lamc3, Itga4
KEGG_PATHWAYHematopoietic cell lineage3Itga4, Cd3d, Cd5
Ijms 16 07976 g006 1024
Figure 6. A theoretical model of the negative regulation of the PA signal pathway by myostatin in skeletal muscle. Plcd4 can catalyze hydrolyzes phosphatidylinositol 4,5-bisphosphate (PIP2) to inositol 1,4,5-trisphosphate (IP3) and DAG [32,33,34]. Diacylglycerol (DAG) can then be phosphorylated by Dgkz to produce phosphatidic acid (PA). Subsequently, PA binds to the FKBP12-rapamycine binding (FRB) domain of mTOR to induce muscle fiber hypertrophy [35,36].
Figure 6. A theoretical model of the negative regulation of the PA signal pathway by myostatin in skeletal muscle. Plcd4 can catalyze hydrolyzes phosphatidylinositol 4,5-bisphosphate (PIP2) to inositol 1,4,5-trisphosphate (IP3) and DAG [32,33,34]. Diacylglycerol (DAG) can then be phosphorylated by Dgkz to produce phosphatidic acid (PA). Subsequently, PA binds to the FKBP12-rapamycine binding (FRB) domain of mTOR to induce muscle fiber hypertrophy [35,36].
Ijms 16 07976 g006

3. Discussion

High-throughput sequencing technologies are currently emerging as a powerful tool to initially elaborate the genetic mechanisms responsible for a phenotype. Previous studies have detected genes related to the hypertrophic muscle observed in myostatin propeptide transgenic mice by employing microarray analysis [30]. The results identified 52 unique genes that are differentially expressed, and myostatin propeptide promotes adult muscle growth and maintenance by up-regulating the expression levels of myogenic regulatory factors and extracellular matrix components and down-regulating protein degradation-related genes as well as mitochondrial ATP synthesis-related genes [30].
To further illustrate the molecular mechanism and identify potential regulatory genes associated with the muscle hypertrophy caused by myostatin inhibition, we employed RNA sequencing technologies to study differences in the expression patterns between MSTN propeptide transgenic mice and wild-type littermate control mice at four months of age; these mice were younger than those used in the microarray study. In the current study, 132 differentially expressed genes were identified in MSTN propeptide transgenic mice compared with control mice based on the RNA sequencing technology and Cufflinks analysis method. Among these 132 genes, most were up-regulated in transgenic mice, and only 35 genes were down-regulated. Functional annotation showed that these differentially expressed genes mainly correlated with muscle development, lipid metabolism, energy metabolism and the immune system. Intriguingly, both of the abovementioned studies revealed that the function of some differentially expressed genes involves the extracellular region. Additionally, pathway analysis found four pathways in both studies: focal adhesion, ECM-receptor interaction, the regulation of the actin cytoskeleton, and cardiac muscle contraction.
In the current study, we found that the expression levels of cytoskeletal and muscle growth and development-related genes were changed in the transgenic mice. Strikingly, Myh6 was down-regulated in the transgenic mice compared with the wild-type controls. A previous study showed that protein synthesis and muscle hypertrophy were enhanced when the expression of myh6 was repressed [37,38,39], suggesting that myh6 is a negative regulator of muscle mass. This finding is consistent with our results, which showed that Myh6 expression was decreased in muscle hypertrophy mice. Moreover, a previous study reported that the mitochondria number, the ratio of mitochondrial DNA to nuclear DNA, and the level of oxidative enzymes used to synthetize ATP all decreased in myostatin-deficient muscle [40]. Other studies showed that Myh6 had higher ATPase activity compared with other myosins [38,41]. Thus, the energy metabolism of myostatin propeptide transgenic mice may decrease due to a decrease in the expression of Myh6, which reduces ATP consumption. Unlike Myh6, the muscle growth-related genes Ky, Actb, Mybph, Actc1, RMST and Actg1 were up-regulated in our study. Among these genes, Ky (kyphoscoliosis peptidase) plays a vital role in muscle growth; the absence of Ky protein leads to muscular dystrophy [42,43]. A Ky mutation results in degenerative myopathy and loss of muscle hypertrophy responses [42]. Ky protein can interact with myosin-binding protein and filamin C [44]. Moreover, Myosin-binding protein H (Mybph) reportedly interacts with Rho kinase 2 (ROCK2) to affect the interaction of actin-myosin [45,46,47]. Additionally, Actb, Actc1 and Actg1 are three types of actins that participate in the maintenance of the cytoskeleton [48]. A functional analysis of these cytoskeletal and muscle development-related genes suggested that myostatin propeptide can enhance muscle mass via the Ky/Mybph/Actins pathway.
Lipids play a critical role in signal transduction; they often act as a secondary messenger to activate downstream pathways [49]. We found that the overexpression of myostatin propeptide activated the phosphatidic acid (PA) signaling pathway, which generates signaling lipids and participates in the cell cycle and cell growth [50,51]. In the current study, we found that Plcd4 and Dgkz, two key enzymes of the PA signaling pathway, were up-regulated in myostatin propeptide transgenic mice. This result suggested a connection between the myostatin signaling and PA signaling pathways. Plcd4 is a member of the delta class of phospholipase C enzymes. Phospholipase C enzymes can catalyze hydrolyzing phosphatidylinositol 4,5-bisphosphate (PIP2) to two intracellular secondary messengers, inositol 1,4,5-trisphosphate (IP3) and diacylglycerol (DAG) [32]. Previous studies reported that DAG could be phosphorylated by Dgkz, and the phosphorylated DAG could be digested to phosphatidic acid (PA). PA can reportedly bind to mTOR to induce muscle fiber hypertrophy [52,53,54,55,56]. Other studies reported that myostatin mediates myoblast differentiation and myotube hypertrophy by inhibiting Akt/mTOR/p70S6 protein signaling [18,20,57]. Overall, we may infer that myostatin can negatively regulate the PA signing pathway to inhibit the mTOR signaling pathway (Figure 4). In addition to activating the mTOR signaling pathway, PA also activates p21-activated kinase 1 (PAK1) to initiate the release of RhoGDI from Rac1, and RhoGD1 subsequently changes the actin dynamics [58].
In biology, redox reactions can frequently store and release biological energy by electron transfer processes. In our study, we found that heme-binding protein Mss51, which is related to redox reactions, was down-regulated in myostatin propeptide transgenic mice. In yeast, Mss51 can regulate the biogenesis of cytochrome c oxidase (COX), the terminal electron transport chain (ETC) oxidase, thereby controlling energy production [59,60]. Previous studies showed that cytochrome c oxidase subunit VIc (Cox6c) expression is decreased in myostatin propeptide transgenic mice [30]. Although we did not detect differential Cox6c expression in our study, the expression of its upstream gene, Mss51, was decreased. Thus, we can infer that myostatin may influence ATP synthesis via Mss51. In addition, the gene Cbr2, which participates in neutralizing excessive reactive oxygen species (ROS), is also regulated by myostatin propeptide [61,62,63]. Myostatin induces the production of ROS, which decreases muscle mass by regulating antioxidant gene expression and mitochondrial biogenesis [64,65,66]. In combination with our results, these findings suggest that myostatin leads to muscle dystrophy by negatively regulating Cbr2 to increase ROS.

4. Experimental Section

4.1. Muscle Tissue Sample Collections

Myostatin propeptide cDNA under the control of rat myosin light chain 1 (MLC1) with its enhancer, SV40 poly A signal sequence, was constructed and served as the transgene construct, and transgenic mice were generated by standard microinjection techniques as previously described [28,29,30]. Male mice (hemizygous genotype for the transgene) from the high-expressing line were mated with B6SJL F1 wild-type females. The mice were housed in cages with a 12-h light/dark cycle, and the room temperature was maintained at 22 °C. The mice were weaned at four weeks of age and given free access to a standard diet (10% kcal fat, ME 3.85 kcal/g). All procedures and animal care were in accordance with the institution guidelines and approved by the Institutional Animal Care and Use Committee of the University of Hawaii, Honolulu, HI, USA. Male mice were sacrificed at four months of age after 8 h of fasting to sample and dissect the muscle tissue. Gastrocnemius muscle samples from both legs were immediately dissected from the carcass, quickly frozen in liquid nitrogen, and stored at −80 °C. Three transgenic male mice and two wild-type littermate male control samples were used for this study.

4.2. RNA Isolation and RNA Sequencing

Total RNA was extracted from freshly frozen muscle tissues using the TRIzol Reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. A magnetic bead homogenizer was used to homogenize the tissue and TRIzol Reagent. The quality and concentration of RNA was detected by NanoDrop ND2000 (Thermo Fisher Scientific, Waltham, MA, USA) spectrophotometry and gel electrophoresis. Solexa sequencing was used to detect the differentially expressed genes in the gastrocnemius of the transgenic and control mice. Five cDNA libraries were constructed using the TruSeq Stranded Total RNA LT Sample Prep Kit (Illumina, Santiago, CA, USA). The library construction and solexa sequencing were performed by Genergy biological technology Limited company (Shanghai, China).

4.3. Quality Control for Raw Sequencing Data

The raw sequencing data were in a FASTQ file, and the data with clean reads were obtained by trimming the adapter contaminants and filtering the low-quality reads. A quality control tool, HTSeq (https://pypi.python.org/pypi/HTSeq), which can provide a quick impression of high throughput data, was utilized in our study. The quality control results are presented in the supplementary information.

4.4. Bioinformatics Analysis

The software Bowtie2 (http://bowtie-bio.sourceforge.net/index.shtml/) was used to merge the chromosome genomes of Mus musculus downloaded from the ENSEMBL database (ftp://ftp.ensembl.org/pub/release-75/fasta/mus_musculus/dna/) and convert them to the FM index, which can store massive genomic data and rapidly be searched. Tophat (http://ccb.jhu.edu/software/tophat/index.shtml), which uses Bowtie as an engine to align the sequencing reads to reference genome, was used to align the clean reads of each individual to the large genomes. Subsequently, the transcriptome was assembled with the Cufflinks software (http://cufflinks.cbcb.umd.edu/), and these assemblies were merged using Cuffmerge in the presence of the reference genome annotation (ftp://ftp.ensembl.org/pub/release-75/gtf/mus_musculus/). Cuffmerge provided a uniform basis to calculate the expression level of transcripts and genes at different conditions.
Cuffdiff, a separate program contained in Cufflink, was used to calculate the expression level of each individual and test the statistical significance between transgenic and control mice. Cuffdiff used FPKM (fragments per kilobase of transcript per million mapped fragments) method to calculate the expression level. In addition, CummeRbund was used to visualize and integrate all results provided by the Cuffdiff analysis. DAVID Bioinformatics Resources were used for gene annotation and pathway analysis.

4.5. Differential Expression Gene Validation by Real Time RT-PCR

Quantitative real-time PCR was performed to validate the differentially expressed genes identified by the Solexa sequencing data. Five differentially expressed genes were randomly selected for quantitative real-time PCR analysis. The Primer Premier 5.0 (PREMIER Biosoft International, Palo Alto, CA, USA) software was used to design the quantitative primers based on the cDNA sequence. RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) was used for RNA reverse transcription following the manufacturer’s instructions, and DNase I and RNase-free (Thermo Fisher Scientific, Waltham, MA, USA) were used to treat the RNA sample to avoid genomic DNA contamination. The qPCR reaction was performed on a LightCycler 480 Real-Time PCR instrument (Roche, Penzberg, Germany) as follows: single cycle of denaturation at 95 °C for 5 min, 40 cycles of denaturation at 95 °C for 15 s, annealing at 60 °C for 15 s and extension at 72 °C for 15 s. SYBR Green (BIO-RAD, Hercules, CA USA) was used for the qPCR reaction, and all reactions were performed in triplicate. The 2−ΔΔCt method was used to calculate the relative expression levels between transgenic and control mice, and all expression levels were normalized to the expression of 18S RNA. Significant differences between groups were assessed with a two-tailed t-test. The difference between groups was considered significant when p < 0.05.

5. Conclusions

In this study, gene expression profiles associated with hypertrophic skeletal muscle in myostatin propeptide transgenic mice were investigated by RNA sequencing. A total of 132 genes were identified and used for function annotation and gene enrichment analysis. The novel up-regulated genes were related to hypertrophic cardiomyopathy, focal adhesion, and PA signaling pathway, and ROS removal processes, suggesting that myostatin directly affects these processes during the regulation of skeletal muscle growth and mass. Two new signaling pathways, the Ky/Mybph/Actins and PA signaling pathways, were uncovered, which may account for myostatin propeptide-induced muscle hypertrophy. Our results provide information to help further understand how myostatin inhibition causes skeletal muscle hypertrophy.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/16/04/7976/s1.

Acknowledgments

This research was supported by an NSFC Outstanding Youth grant (31025026), the priority area project of Ministry of Education (20130146130001), Guangdong Science and Technology Project (2011A020102003), and the Fundamental Research Funds for the Central Universities (2013PY050), Alan M Krassner and Ingeborg v. F. McKee Fund of the Hawaii Community Foundation (ID #: 13P-60310), and the USDA-NIFA funding for HATCH project, administrated by the College of Tropical Agriculture and Human Resources of the University of Hawaii at Manoa.

Author Contributions

Xinyun Li, Shuhong Zhao and Jinzeng Yang conceived and designed the experiments; Yuanxin Miao performed the experiments; Yuanxin Miao, Zhong Xu and Lu Jing analyzed the data; Yuanxin Miao and Jinzeng Yang contributed reagents/materials/analysis tools; Yuanxin Miao, Xinyun Li, Shuhong Zhao and Jinzeng Yang wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, S.J. Regulation of muscle mass by myostatin. Annu. Rev. Cell Dev. Biol. 2004, 20, 61–86. [Google Scholar] [CrossRef] [PubMed]
  2. Parsons, S.A.; Millay, D.P.; Sargent, M.A.; McNally, E.M.; Molkentin, J.D. Age-dependent effect of myostatin blockade on disease severity in a murine model of limb-girdle muscular dystrophy. Am. J. Pathol. 2006, 168, 1975–1985. [Google Scholar] [CrossRef] [PubMed]
  3. McPherron, A.C.; Lawler, A.M.; Lee, S.J. Regulation of skeletal muscle mass in mice by a new tgf-beta superfamily member. Nature 1997, 387, 83–90. [Google Scholar] [CrossRef] [PubMed]
  4. Zimmers, T.A.; Davies, M.V.; Koniaris, L.G.; Haynes, P.; Esquela, A.F.; Tomkinson, K.N.; McPherron, A.C.; Wolfman, N.M.; Lee, S.J. Induction of cachexia in mice by systemically administered myostatin. Science 2002, 296, 1486–1488. [Google Scholar] [CrossRef] [PubMed]
  5. Tsuchida, K. Myostatin inhibition by a follistatin-derived peptide ameliorates the pathophysiology of muscular dystrophy model mice. Acta Myol. 2008, 27, 14–18. [Google Scholar] [PubMed]
  6. McPherron, A.C.; Lee, S.J. Double muscling in cattle due to mutations in the myostatin gene. Proc. Natl. Acad. Sci. USA 1997, 94, 12457–12461. [Google Scholar] [CrossRef] [PubMed]
  7. Kambadur, R.; Sharma, M.; Smith, T.P.; Bass, J.J. Mutations in myostatin (gdf8) in double-muscled belgian blue and piedmontese cattle. Genome Res. 1997, 7, 910–916. [Google Scholar] [PubMed]
  8. Clop, A.; Marcq, F.; Takeda, H.; Pirottin, D.; Tordoir, X.; Bibe, B.; Bouix, J.; Caiment, F.; Elsen, J.M.; Eychenne, F.; et al. A mutation creating a potential illegitimate microrna target site in the myostatin gene affects muscularity in sheep. Nat. Genet. 2006, 38, 813–818. [Google Scholar] [CrossRef] [PubMed]
  9. Williams, M.S. Myostatin mutation associated with gross muscle hypertrophy in a child. N. Engl. J. Med. 2004, 351, 1030–1030. [Google Scholar] [CrossRef] [PubMed]
  10. Guo, T.Q.; Jou, W.; Chanturiya, T.; Portas, J.; Gavrilova, O.; McPherron, A.C. Myostatin inhibition in muscle, but not adipose tissue, decreases fat mass and improves insulin sensitivity. PLoS ONE 2009, 4, e4937. [Google Scholar] [CrossRef] [PubMed]
  11. Bradley, L.; Yaworsky, P.J.; Walsh, F.S. Myostatin as a therapeutic target for musculoskeletal disease. Cell Mol. Life Sci. 2008, 65, 2119–2124. [Google Scholar] [CrossRef] [PubMed]
  12. Attisano, L.; Wrana, J.L. Signal transduction by the TGF-β superfamily. Science 2002, 296, 1646–1647. [Google Scholar] [CrossRef] [PubMed]
  13. Lee, S.J.; McPherron, A.C. Regulation of myostatin activity and muscle growth. Proc Natl. Acad. Sci. USA 2001, 98, 9306–9311. [Google Scholar] [CrossRef] [PubMed]
  14. Joulia-Ekaza, D.; Cabello, G. Myostatin regulation of muscle development: Molecular basis, natural mutations, physiopathological aspects. Exp. Cell Res. 2006, 312, 2401–2414. [Google Scholar] [CrossRef] [PubMed]
  15. Langley, B.; Thomas, M.; Bishop, A.; Sharma, M.; Gilmour, S.; Kambadur, R. Myostatin inhibits myoblast differentiation by down-regulating myod expression. J. Biol. Chem. 2002, 277, 49831–49840. [Google Scholar] [CrossRef] [PubMed]
  16. Miyazawa, K.; Shinozaki, M.; Hara, T.; Furuya, T.; Miyazono, K. Two major smad pathways in TGF-β superfamily signalling. Genes Cells 2002, 7, 1191–1204. [Google Scholar] [CrossRef] [PubMed]
  17. Yan, X.H.; Chen, Y.G. Smad7: Not only a regulator, but also a cross-talk mediator of TGF-β signalling. Biochem. J. 2011, 434, 1–10. [Google Scholar] [CrossRef] [PubMed]
  18. Goodman, C.A.; McNally, R.M.; Hoffmann, F.M.; Hornberger, T.A. Smad3 induces atrogin-1, inhibits mtor and protein synthesis, and promotes muscle atrophy in vivo. Mol. Endocrinol. 2013, 27, 1946–1957. [Google Scholar] [CrossRef] [PubMed]
  19. Lokireddy, S.; McFarlane, C.; Ge, X.J.; Zhang, H.M.; Sze, S.K.; Sharma, M.; Kambadur, R. Myostatin induces degradation of sarcomeric proteins through a smad3 signaling mechanism during skeletal muscle wasting. Mol. Endocrinol. 2011, 25, 1936–1949. [Google Scholar] [CrossRef] [PubMed]
  20. Trendelenburg, A.U.; Meyer, A.; Rohner, D.; Boyle, J.; Hatakeyama, S.; Glass, D.J. Myostatin reduces Akt/TORC1/p70S6K signaling, inhibiting myoblast differentiation and myotube size. Am. J. Physiol. Cell Physiol. 2009, 296, C1258–C1270. [Google Scholar] [CrossRef] [PubMed]
  21. Bogdanovich, S.; Krag, T.O.; Barton, E.R.; Morris, L.D.; Whittemore, L.A.; Ahima, R.S.; Khurana, T.S. Functional improvement of dystrophic muscle by myostatin blockade. Nature 2002, 420, 418–421. [Google Scholar] [CrossRef] [PubMed]
  22. Diel, P.; Schiffer, T.; Geisler, S.; Hertrampf, T.; Mosler, S.; Schulz, S.; Wintgens, K.F.; Adler, M. Analysis of the effects of androgens and training on myostatin propeptide and follistatin concentrations in blood and skeletal muscle using highly sensitive immuno PCR. Mol. Cell. Endocrinol. 2010, 330, 1–9. [Google Scholar] [CrossRef] [PubMed]
  23. Hill, J.J.; Davies, M.V.; Pearson, A.A.; Wang, J.H.; Hewick, R.M.; Wolfman, N.M.; Qiu, Y. The myostatin propeptide and the follistatin-related gene are inhibitory binding proteins of myostatin in normal serum. J. Biol. Chem. 2002, 277, 40735–40741. [Google Scholar] [CrossRef] [PubMed]
  24. Kishioka, Y.; Thomas, M.; Wakamatsu, J.I.; Hattori, A.; Sharma, M.; Kambadur, R.; Nishimura, T. Decorin enhances the proliferation and differentiation of myogenic cells through suppressing myostatin activity. J. Cell. Physiol. 2008, 215, 856–867. [Google Scholar] [CrossRef] [PubMed]
  25. Nozaki, M.; Li, Y.; Zhu, J.; Ambrosio, F.; Uehara, K.; Fu, F.H.; Huard, J. Improved muscle healing after contusion injury by the inhibitory effect of suramin on myostatin, a negative regulator of muscle growth. Am. J. Sports Med. 2008, 36, 2354–2362. [Google Scholar] [CrossRef] [PubMed]
  26. Huang, Z.Q.; Chen, X.L.; Yu, B.; He, J.; Chen, D.W. Microrna-27a promotes myoblast proliferation by targeting myostatin. Biochem. Biophys. Res. Commun. 2012, 423, 265–269. [Google Scholar] [CrossRef] [PubMed]
  27. Miretti, S.; Martignani, E.; Accornero, P.; Baratta, M. Functional effect of mir-27b on myostatin expression: A relationship in piedmontese cattle with double-muscled phenotype. BMC Genomics 2013, 14, 194. [Google Scholar] [CrossRef] [PubMed]
  28. Yang, J.; Zhao, B. Postnatal expression of myostatin propeptide cdna maintained high muscle growth and normal adipose tissue mass in transgenic mice fed a high-fat diet. Mol. Reprod. Dev. 2006, 73, 462–469. [Google Scholar] [CrossRef] [PubMed]
  29. Yang, J.; Ratovitski, T.; Brady, J.P.; Solomon, M.B.; Wells, K.D.; Wall, R.J. Expression of myostatin pro domain results in muscular transgenic mice. Mol. Reprod. Dev. 2001, 60, 351–361. [Google Scholar] [CrossRef] [PubMed]
  30. Zhao, B.; Li, E.J.; Wall, R.J.; Yang, J. Coordinated patterns of gene expressions for adult muscle build-up in transgenic mice expressing myostatin propeptide. BMC Genomics 2009, 10, 305. [Google Scholar] [CrossRef] [PubMed]
  31. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed]
  32. Laxalt, A.M.; Munnik, T. Phospholipid signalling in plant defence. Curr. Opin. Plant Biol. 2002, 5, 332–338. [Google Scholar] [CrossRef] [PubMed]
  33. Fukami, K.; Takenaka, K.; Nagano, K.; Takenawa, T. Growth factor-induced promoter activation of murine phospholipase C δ4 gene. Eur. J. Biochem. 2000, 267, 28–36. [Google Scholar] [CrossRef] [PubMed]
  34. Liu, N.; Fukami, K.; Yu, H.; Takenawa, T. A new phospholipase C δ4 is induced at s-phase of the cell cycle and appears in the nucleus. J. Biol. Chem. 1996, 271, 355–360. [Google Scholar] [CrossRef] [PubMed]
  35. Veverka, V.; Crabbe, T.; Bird, I.; Lennie, G.; Muskett, F.W.; Taylor, R.J.; Carr, M.D. Structural characterization of the interaction of mTOR with phosphatidic acid and a novel class of inhibitor: Compelling evidence for a central role of the FRB domain in small molecule-mediated regulation of mTOR. Oncogene 2008, 27, 585–595. [Google Scholar] [CrossRef] [PubMed]
  36. Fang, Y.; Vilella-Bach, M.; Bachmann, R.; Flanigan, A.; Chen, J. Phosphatidic acid-mediated mitogenic activation of mTOR signaling. Science 2001, 294, 1942–1945. [Google Scholar] [CrossRef] [PubMed]
  37. Gao, X.M.; Kiriazis, H.; Moore, X.L.; Feng, X.H.; Sheppard, K.; Dart, A.; Du, X.J. Regression of pressure overload-induced left ventricular hypertrophy in mice. Am. J. Physiol. Heart Circ. Physiol. 2005, 288, H2702–H2707. [Google Scholar] [CrossRef] [PubMed]
  38. Swynghedauw, B. Developmental and functional adaptation of contractile proteins in cardiac and skeletal muscles. Physiol. Rev. 1986, 66, 710–771. [Google Scholar] [PubMed]
  39. Ojamaa, K.; Petrie, J.F.; Balkman, C.; Hong, C.; Klein, I. Posttranscriptional modification of myosin heavy-chain gene expression in the hypertrophied rat myocardium. Proc. Natl. Acad. Sci. USA 1994, 91, 3468–3472. [Google Scholar] [CrossRef] [PubMed]
  40. Amthor, H.; Macharia, R.; Navarrete, R.; Schuelke, M.; Brown, S.C.; Otto, A.; Voit, T.; Muntoni, F.; Vrbova, G.; Partridge, T.; et al. Lack of myostatin results in excessive muscle growth but impaired force generation. Proc. Natl. Acad. Sci. USA 2007, 104, 1835–1840. [Google Scholar] [CrossRef] [PubMed]
  41. Tyska, M.J.; Warshaw, D.M. The myosin power stroke. Cell Motil Cytoskeleton 2002, 51, 1–15. [Google Scholar] [CrossRef] [PubMed]
  42. Blanco, G.; Coulton, G.R.; Biggin, A.; Grainge, C.; Moss, J.; Barrett, M.; Berquin, A.; Marechal, G.; Skynner, M.; van Mier, P.; et al. The kyphoscoliosis (ky) mouse is deficient in hypertrophic responses and is caused by a mutation in a novel muscle-specific protein. Hum. Mol. Genet. 2001, 10, 9–16. [Google Scholar] [CrossRef] [PubMed]
  43. Vargas, J.D.; Culetto, E.; Ponting, C.P.; Miguel-Aliaga, I.; Davies, K.E.; Sattelle, D.B. Cloning and developmental expression analysis of ltd-1, the caenorhabditis elegans homologue of the mouse kyphoscoliosis (ky) gene. Mech. Dev. 2002, 117, 289–292. [Google Scholar] [CrossRef] [PubMed]
  44. Beatham, J.; Romero, R.; Townsend, S.K.M.; Hacker, T.; van der Ven, P.F.M.; Blanco, G. Filamin C interacts with the muscular dystrophy KY protein and is abnormally distributed in mouse KY deficient muscle fibres. Hum. Mol. Genet. 2004, 13, 2863–2874. [Google Scholar] [CrossRef] [PubMed]
  45. Conti, A.; Riva, N.; Pesca, M.; Iannaccone, S.; Cannistraci, C.V.; Corbo, M.; Previtali, S.C.; Quattrini, A.; Alessio, M. Increased expression of myosin binding protein H in the skeletal muscle of amyotrophic lateral sclerosis patients. Biochim. Biophys. Acta 2014, 1842, 99–106. [Google Scholar] [CrossRef] [PubMed]
  46. Lin, Y.; Ishikawa, R.; Okagaki, T.; Ye, L.H.; Kohama, K. Stimulation of the ATP-dependent interaction between actin and myosin by a myosin-binding fragment of smooth muscle caldesmon. Cell Motil. Cytoskeleton 1994, 29, 250–258. [Google Scholar] [CrossRef] [PubMed]
  47. dos Remedios, C.G.; Chhabra, D.; Kekic, M.; Dedova, I.V.; Tsubakihara, M.; Berry, D.A.; Nosworthy, N.J. Actin binding proteins: Regulation of cytoskeletal microfilaments. Physiol. Rev. 2003, 83, 433–473. [Google Scholar] [CrossRef] [PubMed]
  48. Doherty, G.J.; McMahon, H.T. Mediation, modulation, and consequences of membrane-cytoskeleton interactions. Annu. Rev. Biophys. 2008, 37, 65–95. [Google Scholar] [CrossRef] [PubMed]
  49. Wang, X.; Devaiah, S.P.; Zhang, W.; Welti, R. Signaling functions of phosphatidic acid. Prog. Lipid Res. 2006, 45, 250–278. [Google Scholar] [CrossRef]
  50. Janda, M.; Planchais, S.; Djafi, N.; Martinec, J.; Burketova, L.; Valentova, O.; Zachowski, A.; Ruelland, E. Phosphoglycerolipids are master players in plant hormone signal transduction. Plant Cell Rep. 2013, 32, 839–851. [Google Scholar] [CrossRef] [PubMed]
  51. Munnik, T. Phosphatidic acid: An emerging plant lipid second messenger. Trends Plant Sci. 2001, 6, 227–233. [Google Scholar] [CrossRef] [PubMed]
  52. You, J.S.; Lincoln, H.C.; Kim, C.R.; Frey, J.W.; Goodman, C.A.; Zhong, X.P.; Hornberger, T.A. The role of diacylglycerol kinase zeta and phosphatidic acid in the mechanical activation of mammalian target of rapamycin (mTOR) signaling and skeletal muscle hypertrophy. J. Biol. Chem. 2014, 289, 1551–1563. [Google Scholar] [CrossRef] [PubMed]
  53. Goodman, C.A.; Frey, J.W.; Mabrey, D.M.; Jacobs, B.L.; Lincoln, H.C.; You, J.S.; Hornberger, T.A. The role of skeletal muscle mTOR in the regulation of mechanical load-induced growth. J. Physiol. 2011, 589, 5485–5501. [Google Scholar] [CrossRef] [PubMed]
  54. Hornberger, T.A. Mechanotransduction and the regulation of mTORC1 signaling in skeletal muscle. Int. J. Biochem. Cell Biol. 2011, 43, 1267–1276. [Google Scholar] [CrossRef] [PubMed]
  55. Gharbi, S.I.; Rincon, E.; Avila-Flores, A.; Torres-Ayuso, P.; Almena, M.; Cobos, M.A.; Albar, J.P.; Merida, I. Diacylglycerol kinase zeta controls diacylglycerol metabolism at the immunological synapse. Mol. Biol. Cell 2011, 22, 4406–4414. [Google Scholar] [CrossRef] [PubMed]
  56. Toschi, A.; Lee, E.; Xu, L.; Garcia, A.; Gadir, N.; Foster, D.A. Regulation of mTORC1 and mTORC2 complex assembly by phosphatidic acid: Competition with rapamycin. Mol. Cell. Biol. 2009, 29, 1411–1420. [Google Scholar] [CrossRef] [PubMed]
  57. Hulmi, J.J.; Silvennoinen, M.; Lehti, M.; Kivela, R.; Kainulainen, H. Altered redd1, myostatin, and akt/mtor/foxo/mapk signaling in streptozotocin-induced diabetic muscle atrophy. Am. J. Physiol. Endocrinol. Metab. 2012, 302, E307–E315. [Google Scholar] [CrossRef] [PubMed]
  58. Abramovici, H.; Mojtabaie, P.; Parks, R.J.; Zhong, X.P.; Koretzky, G.A.; Topham, M.K.; Gee, S.H. Diacylglycerol kinase zeta regulates actin cytoskeleton reorganization through dissociation of Rac1 from RhoGDI. Mol. Biol. Cell 2009, 20, 2049–2059. [Google Scholar] [CrossRef] [PubMed]
  59. Soto, I.C.; Fontanesi, F.; Myers, R.S.; Hamel, P.; Barrientos, A. A heme-sensing mechanism in the translational regulation of mitochondrial cytochrome c oxidase biogenesis. Cell Metab. 2012, 16, 801–813. [Google Scholar] [CrossRef] [PubMed]
  60. Fontanesi, F.; Soto, I.C.; Horn, D.; Barrientos, A. Mss51 and Ssc1 facilitate translational regulation of cytochrome c oxidase biogenesis. Mol. Cell. Biol. 2010, 30, 245–259. [Google Scholar] [CrossRef] [PubMed]
  61. Nakagawa, J.; Ishikura, S.; Asami, J.; Isaji, T.; Usami, N.; Hara, A.; Sakurai, T.; Tsuritani, K.; Oda, K.; Takahashi, M.; et al. Molecular characterization of mammalian dicarbonyl/l-xylulose reductase and its localization in kidney. J. Biol. Chem. 2002, 277, 17883–17891. [Google Scholar] [CrossRef] [PubMed]
  62. Liu, S.; Ma, L.; Huang, W.; Shai, Y.; Ji, X.; Ding, L.; Liu, Y.; Yu, L.; Zhao, S. Decreased expression of the human carbonyl reductase 2 gene hcr2 in hepatocellular carcinoma. Cell Mol. Biol. Lett. 2006, 11, 230–241. [Google Scholar] [CrossRef] [PubMed]
  63. Hoffmann, F.; Maser, E. Carbonyl reductases and pluripotent hydroxysteroid dehydrogenases of the short-chain dehydrogenase/reductase superfamily. Drug Metab. Rev. 2007, 39, 87–144. [Google Scholar] [CrossRef] [PubMed]
  64. Barbieri, E.; Sestili, P. Reactive oxygen species in skeletal muscle signaling. J. Signal Transduct. 2012, 2012. [Google Scholar] [CrossRef]
  65. Lee, K.-P.; Shin, Y.J.; Cho, S.C.; Lee, S.-M.; Bahn, Y.J.; Kim, J.Y.; Kwon, E.-S.; Park, S.C.; Rhee, S.G.; Woo, H.A.; et al. Peroxiredoxin 3 has a crucial role in the contractile function of skeletal muscle by regulating mitochondrial homeostasis. Free Radic. Biol. Med. 2014, 77, 298–306. [Google Scholar] [CrossRef] [PubMed]
  66. Sriram, S.; Subramanian, S.; Juvvuna, P.K.; Ge, X.; Lokireddy, S.; McFarlane, C.D.; Wahli, W.; Kambadur, R.; Sharma, M. Myostatin augments muscle-specific ring finger protein-1 expression through an Nf-kB independent mechanism in SMAD3 null muscle. Mol. Endocrinol. 2014, 28, 317–330. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Miao, Y.; Yang, J.; Xu, Z.; Jing, L.; Zhao, S.; Li, X. RNA Sequencing Identifies Upregulated Kyphoscoliosis Peptidase and Phosphatidic Acid Signaling Pathways in Muscle Hypertrophy Generated by Transgenic Expression of Myostatin Propeptide. Int. J. Mol. Sci. 2015, 16, 7976-7994. https://doi.org/10.3390/ijms16047976

AMA Style

Miao Y, Yang J, Xu Z, Jing L, Zhao S, Li X. RNA Sequencing Identifies Upregulated Kyphoscoliosis Peptidase and Phosphatidic Acid Signaling Pathways in Muscle Hypertrophy Generated by Transgenic Expression of Myostatin Propeptide. International Journal of Molecular Sciences. 2015; 16(4):7976-7994. https://doi.org/10.3390/ijms16047976

Chicago/Turabian Style

Miao, Yuanxin, Jinzeng Yang, Zhong Xu, Lu Jing, Shuhong Zhao, and Xinyun Li. 2015. "RNA Sequencing Identifies Upregulated Kyphoscoliosis Peptidase and Phosphatidic Acid Signaling Pathways in Muscle Hypertrophy Generated by Transgenic Expression of Myostatin Propeptide" International Journal of Molecular Sciences 16, no. 4: 7976-7994. https://doi.org/10.3390/ijms16047976

APA Style

Miao, Y., Yang, J., Xu, Z., Jing, L., Zhao, S., & Li, X. (2015). RNA Sequencing Identifies Upregulated Kyphoscoliosis Peptidase and Phosphatidic Acid Signaling Pathways in Muscle Hypertrophy Generated by Transgenic Expression of Myostatin Propeptide. International Journal of Molecular Sciences, 16(4), 7976-7994. https://doi.org/10.3390/ijms16047976

Article Metrics

Back to TopTop