Abstract
Death-associated protein (DAP) undergoes substantial changes in expression during turkey skeletal muscle development, decreasing from the 18 day embryonic stage to 1 day posthatch, and again from 1 day posthatch to 16 weeks of age. These changes suggest that DAP plays an important role at critical stages of the developmental process. The objective of this study was to elucidate the role of DAP in muscle development by examining the effect of reduced DAP expression on global gene expression in proliferating and differentiating turkey pectoralis major muscle satellite cells. Small interfering RNA was used to knock down expression of DAP and the transcriptome was subsequently profiled using a turkey skeletal muscle long oligonucleotide microarray. Microarray data were corroborated using quantitative real-time PCR. In proliferating cells, 458 loci, resulting in 378 uniquely annotated genes, showed differential expression (false discovery rate, FDR < 0.05). Pathway analysis highlighted altered eukaryotic translational initiation factors (eIFs) signaling, protein ubiquitination, sirtuin signaling, and mechanistic target of rapamycin (mTOR) signaling as the primary pathways affected in the knockdown proliferating cells. The findings underpinned the potential DAP involvement in cell proliferation of turkey satellite cells through the coordination between protein synthesis and cell cycle. In differentiating cells, 270 loci, accounting for 189 unique genes, showed differential expression (FDR < 0.05). Decreased expression of genes encoding various myofibrillar proteins and proteins involved in sarcoplasmic reticulum calcium flux suggests that DAP may affect regulation of calcium homeostasis and cytoskeleton signaling. This study provides the first evidence that reduced expression of DAP significantly alters the transcriptome profile of pectoralis major muscle satellite cells, thereby reducing proliferation and differentiation.
Introduction
Death-associated protein (DAP) is a highly conserved, 15-kDa, proline-rich protein with two phosphorylation sites (; ). It was first identified by in studies aimed at identifying novel genes with functional relevance in apoptosis. Subsequent cell culture studies found that in nutrient-rich conditions, DAP is a substrate for phosphorylation by the mechanistic target of rapamycin (mTOR) (). When cells are subjected to amino acid starvation, mTOR activity is down-regulated and DAP undergoes dephosphorylation. Analysis of cells grown under nutrient-rich and nutrient-deprived conditions suggests that the phosphorylated form is functionally silent, whereas dephosphorylated DAP is the active form that acts as an inhibitor of cell autophagy (). Subsequently, others have shown that DAP is involved in regulation of autophagy and apoptosis induced by subtilase cytotoxin (SubAB) of some strains of Shiga-toxigenic E. coli (). Reduced expression of DAP is also associated with adverse clinical outcomes for human breast () and colorectal cancer () owing to altered regulation of apoptosis and autophagy.
In previous studies, our group investigated changes in gene expression in turkey pectoralis major muscle as a function of temporal development (). Based on the following criteria, DAP was one of several differentially expressed (DE) genes selected for further studies. First, DAP expression changed substantially, being highest in embryonic muscle undergoing hyperplasia, followed by a five-fold decrease in neonatal muscle undergoing hypertrophy, and decreased further in muscle from 16-wk market-age turkeys (). Second, little was known about the function of DAP, and most importantly, there was no previous evidence for its role in muscle development.
In subsequent studies, knockdown of DAP expression in turkey pectoralis major muscle satellite cells by small interfering RNA (siRNA) resulted in reduction of proliferation by up to 50% and nearly complete inhibition of the differentiation of satellite cells into myotubes (). Chicken satellite cells transfected with DAP cDNA showed that overexpression of DAP resulted in increased proliferation, differentiation, and myotube diameter, but had no effect on satellite cell apoptosis (). These studies were the first to suggest that DAP may play an important role in pectoralis major muscle development.
The objective of the current study was to initiate studies to define the role of DAP in the early stages of muscle growth and development by testing the hypothesis that altered DAP expression modulates expression of other genes that are critical to pectoralis major muscle development. siRNA was used to reduce DAP expression in proliferating and differentiating turkey myogenic satellite cells. Transcriptome differences between control and knockdown cell cultures were measured using a turkey skeletal muscle long oligonucleotide microarray (TSKLMO) (). Our results provide the first evidence of global gene expression changes upon knockdown of DAP expression. Pathway analysis suggests potential involvement of DAP in protein synthesis, calcium homeostasis, and muscle mass accumulation.
Materials and Methods
Isolation of Turkey Myogenic Satellite Cells
Randombred Control Line 2 (RBC2) Pectoralis major muscle satellite cells, previously harvested according to and stored in liquid nitrogen, were used for this study. The RBC2 line represents a 1967 commercial line turkey that has been maintained at The Ohio State University, Agricultural Research and Development Center Poultry Research Unit, without conscious selection for any traits (; ). An ethics approval statement was not required because no animals were sacrificed for the current study. The cells were analyzed for proliferation and differentiation characteristics as previously described (, ).
Cell Culture
Growth of control and DAP-knockdown satellite cells were performed at The Ohio State University as described in detail by . Briefly, satellite cells were transfected with siRNA targeted against DAP or were mock-transfected (without siRNA) 24 h after plating (Figure 1). Proliferation was induced as previously described (, ) and the extent of proliferation was measured by DNA concentration (). After 72 h of proliferation, cell differentiation was induced by changes in growth media as described by (Figure 1), The differentiation assay was adapted from and Yun et al. (1997) as described by . The extent of differentiation was determined as a function of muscle-specific creatine kinase levels at 48 h after induction of differentiation.
FIGURE 1
RNA Isolation
Satellite cell cultures were frozen at −80°C after 72 h of proliferation or 48 h of differentiation and shipped on dry ice to Michigan State University for RNA isolation and subsequent microarray hybridization. Frozen cells were thawed on ice for 5 min prior to isolation of total RNA with TRIReagent (Molecular Research Center, Inc., Cincinnati, OH, United States) following manufacturer instructions. Contaminating genomic DNA was removed by treatment with RNA-free DNase (Ambion Inc., Austin, TX, United States). Total RNA was quantified using a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States) and sample integrity determined with an Agilent 2100 Bioanalyzer (Santa Clara, CA, United States). Samples with an RNA Integrity Number (RIN) of >8.0 were considered acceptable for use with the microarray.
Microarray Experimental Design
Design, construction, and validation of the TSKMLO array are described in . The microarray was constructed based on cDNA libraries gathered at critical stages in turkey skeletal muscle development (). Information regarding the TSKMLO platform can be found at the National Center for Biotechnology (NCBI) Gene Expression Omnibus (GEO) database (platform accession: GPL9788). In the current study, RNAs isolated from DAP knockdown samples were directly compared to those of the control treatment. A total of four replicates was used for each experiment (proliferation or differentiation) and dye swaps were included in the experimental design to minimize dye bias, resulting in eight arrays per experiment and 16 arrays for the entire knockdown study.
RNA Amplification and Microarray Hybridization
Total RNA was amplified, dye-coupled, and purified using the Amino Allyl MessageAmpTM II aRNA Amplification Kit (Ambion, Inc.) according to the manufacturer’s protocol. Due to low levels of RNA isolated from cell culture samples, two rounds of amplification were performed for all samples. Microarray preparation and hybridization were carried out as described in . All spot intensities were exported as GenePix Results (GPR) files for statistical analysis.
Microarray Statistical Analysis and Gene Annotation
The GPR files were subjected to statistical analysis as described in . Briefly, dye intensity bias was normalized using the “normexp” background correction method based on . Normalized data were described as log2 fluorescent intensities ratio (Cy5/Cy3) or M-value, and statistically analyzed with a linear model using LIMMA () based on overall intercept, the dye assignment for samples from the DAP knockdown treatment, and the random effect of cells derived from the treatment. Mismatch and negative control spots on the microarray were used to confirm that hybridization to the array was specific (). Raw Cy5 and Cy3 intensities, M-value, and LOESS-normalized log2 average intensities (A) were submitted with original GPR files to the NCBI GEO (platform accession: GPL9788, series accession: GSE35660).
Differences in expression levels between DAP knockdown and control were considered significant with an estimated false discovery rate (FDR) <0.05. In this experiment, the comparison of interest was overall differences between DAP knockdown treatment and control within each developmental stage (proliferation or differentiation). Fold change was defined as the ratio of expression of particular gene in DAP-knockdown cells relative to control cell culture. Gene annotation of the array probes was conducted by iterative NCBI BLAST searches (blast.ncbi.nlm.nih.gov) of the original cDNA sequences used to design the array probes. Sequences were first queried against the turkey refseq gene list (genome annotation 103), then the NCBI refseq database, followed by the turkey genome (v5.1) and the NCBI NR database. Sequences that failed to identify significant matches were designated as novel transcripts. Several sequences aligned to introns within described genes in the turkey genome BLAST searches. These were assigned to the described genes. Sequences with significant alignment to the turkey genome outside of described genes and without matches to other databases were designated as undescribed transcripts.
Confirmation of Gene Expression by Quantitative Real-Time PCR
Six genes from the proliferation experiment and 9 genes from the differentiation experiment were chosen for gene expression analysis using quantitative real-time PCR (qPCR). Genes were chosen based on function and comparably high levels of differential expression as seen in the microarray analysis. Primers (Supplementary Table S1) were designed using Primer Express 3.0 software (Applied Biosystems, Foster City, CA, United States) and synthesized by Operon, Inc. (Huntsville, AL, United States). Sample Amino Allyl-aRNAs (2 μg), prepared for microarray hybridization using the Amino Allyl MessageAmpTM II aRNA Amplification Kit (Ambion, Inc.), were reverse-transcribed into cDNA using Superscript III (Invitrogen). Synthesized cDNA was quantified using a Nanodrop. Reactions were run using 6–10 ng cDNA template, 1.2 μM primer mix, and Power SYBR Green PCR Master Mix (Applied Biosystems) in an ABI Prism 7900HT Sequence Detection System (Applied Biosystems). Data were analyzed using the 2–ΔΔCt method (). Based on data indicating no change in Ct for the described treatment conditions (data not shown), UNC-5 homolog B (UNC5B) and Integral Membrane Protein 2C (IMP2C) were used as endogenous control genes for proliferation and differentiation, respectively. Statistical analysis was performed using a t-test (SigmaPlot, San Jose, CA, United States) and results are expressed as a fold change of the DAP knockdown relative to the control.
Functional Analysis
DE transcripts (FDR < 0.05) were subjected to pathway analysis using Qiagen Ingenuity Pathway Analysis (IPA) software (Qiagen, Redwood City, CA, United States). Canonical pathways, indicative of cell signaling and metabolic pathways relevant to the DE genes, were obtained from IPA analysis based on algorithms developed by .
Results
Differential Expression of Satellite Cell Genes
A combined total of 647 loci (representing up to 559 genes) was DE between DAP knockdown and control treatments in the turkey satellite cells (FDR < 0.05). In the proliferating cells, 458 loci (378 unique genes) were DE (Supplementary Table S2), whereas 270 loci (189 unique genes) were DE in differentiation (Figure 2). During proliferation, more genes were down-regulated (244) than up-regulated (214) upon knockdown of DAP (Supplementary Table S2).
FIGURE 2
As with knockdown of proliferating satellite cells, more genes in the differentiating satellite cells were down-regulated (179) than up-regulated (91) (Supplementary Table S3). As expected, expression of most of the myofibrillar proteins were decreased. Those include two loci of alpha actin (ACTA1, FC = −19 and −12), various isoforms of tropomyosin (TPM, FC ranging from −4 to −15), troponin C (TNNC, FC ranging between −2 and −18), troponin I (TNNI, FC ranging from −2.6 to −13), troponin T (TNNT, FC ranging from −1.9 to −35), nebulin (LOC109368623, FC = −4), and titin (TTN, FC = −3.5). In addition, the gene for obscurin (OBSCN), a protein linking the myofibril to the sarcoplasmic reticulum (SR), was down-reglated 2.2-fold.
Eighty loci showed differential expression in both proliferation and differentiation knockdown experiments (Table 1). Expression of 72 loci was directionally consistent between both the proliferating and differentiating cells. Differences in fold change between proliferation and differentiation ranged from 0.009 to 23.9 with an average of 2.6 (Standard Deviation = 3.61). The largest FC deviation was observed for eukaryotic translation initiation factor 3 subunit E (eIF3E) with FC = 29.3 and 5.3 in proliferation and differentiation, respectively. This initiation factor is a component of the initiation factor 3 (eIF-3) complex, required in the initiation of protein synthesis. In contrast, FC values were essentially identical for CDV3 (2.37 vs 2.36). Protein CDV3 homolog, also known as carnitine deficiency-associated gene expressed in ventricle 3, encodes the protein Histone H4.
TABLE 1
 |
Differentially expressed transcripts identified in both proliferation and differentiation knockdown experiments.
FC, fold change, positive and negative FC indicates increased and decreased, respectively, abundance in DAP-knockdown satellite cells. FDR, false discovery rate. Green highlighting indicates down-regulated genes; red highlighting indicates up-regulated genes.
Four gene loci that were up-regulated during proliferation were down-regulated in differentiation (Table 1). These loci include heat shock protein 25 (HSP25, FC = 9.1 vs −19.8), tropomyosin 4 (TPM4, FC = 5.2 vs −4.1), 3-hydroxyacyl-CoA dehydratase 1 (HACD1, FC = 4.3 vs −3.3), and A-kinase anchoring protein 13 (AKAP13, FC = 2.8 vs −2.2). Conversely, three gene loci that were down-regulated during proliferation were up-regulated in differentiation. These include deoxyribodipyrimidine photo-lyase-like (LOC100551374, FC = −2.19 vs 2.35), nuclear factor of activated T-cells (NFATC2, FC = −2.20 vs 2.53), and procollagen C-endopeptidase enhancer 2 (PCOLCE2, FC = −6.4 vs 5.1).
Confirmation of Gene Expression
Differential gene expression identified by microarray analysis was corroborated by performing qPCR on six genes from the proliferation experiment and nine genes from the differentiation experiment (Supplementary Table S1). All genes selected show similar overall trends in differential expression in both microarray and in qPCR results for proliferation (Figure 3) and differentiation (Figure 4).
FIGURE 3
FIGURE 4
Functional Pathway Analysis
DE transcripts (FDR < 0.05) were subjected to IPA analysis using their official gene names and canonical pathways considered in interpreting results. The most affected canonical pathways within DAP-knockdown proliferating and differentiation satellite cells are shown in Table 2. Knockdown of DAP in turkey satellite cells during proliferation altered expression of several eukaryotic initiation factor (eIF) genes as well as mTOR signaling genes (Figure 5). During differentiation, the calcium signaling pathway (Table 2, Figure 6) was the primary canonical pathway altered by DAP knockdown. In addition, actin cytoskeleton signaling, integrin-linked kinase (ILK) signaling, regulation of actin-based motility by Rho and protein kinase A are also revealed by IPA as affected canonical pathways (Table 2). IPA calculates activation (z) scores that infer the activation states of genes with predicted functional interactions within molecular networks. All of these enriched pathways exhibited a negative z-score, indicating a predicted inactivation. The impact of differential expression on molecular functions include aspects of cell morphology, cellular organization, and cellular maintenance (Table 3). Many of the genes showing significant changes in expression (FDR < 0.05) are translated into proteins comprising the myofibrillar architecture including myosin, actin, tropomyosin, the troponin subunits, nebulin, and titin corresponding with the observed failure to form myotubes in our previous study on DAP-knockdown in satellite cells ().
TABLE 2
| Stage | Canonical pathways | p-value |
| Proliferation | EIF2 signaling | 9.68 × 10–11 |
| Protein ubiquitination pathway | 9.98 × 10–8 | |
| Sirtuin signaling pathway | 2.69 × 10–7 | |
| Regulation of eIF4 and p70S6K signaling | 1.56 × 10–6 | |
| mTOR signaling | 1.62 × 10–6 | |
| Differentiation | Calcium signaling | 7.40 × 10–15 |
| Actin cytoskeleton signaling | 3.14 × 10–8 | |
| ILK signaling | 4.99 × 10–7 | |
| Regulation of actin-based motility by Rho | 1.05 × 10–6 | |
| Protein kinase A signaling | 1.07 × 10–6 |
Top five canonical pathways affected by DAP knockdown in proliferating and differentiating satellite cells.
FIGURE 5
FIGURE 6
TABLE 3
| Stage | Altered functions | p-value range |
| Proliferation | Protein synthesis | 7.64 × 10–6–6.36 × 10–15 |
| Cell death and survival | 8.69 × 10–6–2.59 × 10–9 | |
| Cellular assembly and organization | 4.35 × 10–8–4.35 × 10–8 | |
| DNA replication, recombination, and repair | 4.35 × 10–8–4.35 × 10–8 | |
| Cellular movement | 9.06 × 10–7–1.11 × 10–7 | |
| Differentiation | Cell morphology | 3.04 × 10–4–5.88 × 10–19 |
| Cellular assembly and organization | 3.07 × 10–4–2.02 × 10–16 | |
| Cellular function and maintenance | 2.77 × 10–4–8.89 × 10–16 | |
| Cellular development | 3.37 × 10–4–3.69 × 10–15 | |
| Cellular growth and proliferation | 3.37 × 10–4–3.69 × 10–15 |
Top five altered molecular and cellular functions affected by DAP knockdown in proliferating and differentiating satellite cells.
Pathway analysis of the 80 DE loci common to both proliferation and differentiation knockdown experiments showed the most affected canonical pathways include nitric oxide signaling in the cardiovascular system, neuregulin signaling, mTOR signaling, glucocorticoid receptor signaling and unfolded protein response (Table 4). In an IPA comparison analysis, the largest activation score was observed for nitric oxide signaling in the cardiovascular system (activated, z = 1.34, p = 8.79 × 10–6). This result is based primarily on expression of heatshock proteins (HSP90B1, HSP90AA1, and HSP25), AKT1, and AMPK (PRKAA1).
TABLE 4
| Canonical pathways | p-value |
| Nitric oxide signaling in the cardiovascular system | 8.79 × 10–6 |
| Neuregulin signaling | 1.56 × 10–4 |
| mTOR signaling | 3.14 × 10–4 |
| Glucocorticoid receptor signaling | 3.63 × 10–4 |
| Unfolded protein response | 5.30 × 10–4 |
Top five canonical pathways with shared genes affected by DAP knockdown in proliferating and differentiating satellite cells based on p-value.
Discussion
Myocyte proliferation and differentiation are the result of tightly regulated, temporal changes in gene expression (). The processes are not only required for skeletal muscle growth and development, but also for muscle regeneration to maintain integrity of mature adult skeletal muscle. Biological pathways associated with muscle development have been intensively studied; however, many genes involved in this process are as of yet unidentified or their specific roles are not fully elucidated. In a previous transcriptome analysis to characterize temporal changes in gene expression, our group found that DAP transcription was highest in embryonic pectoralis major muscle undergoing hyperplasia, followed by a substantial progressive decline in expression as muscle development shifted to hypertrophy in the 1-day-old hatchling (). Transcript levels of DAP continued to decline between 1 day and 16 weeks of age. Subsequently, we performed siRNA knockdown experiments to elucidate the role of DAP in proliferating and differentiating of pectoralis major muscle satellite cells. observed that DAP knockdown resulted in substantial inhibition of both satellite cell proliferation and differentiation in vitro. Similar results were observed in chicken satellite cell cultures (). Conversely, overexpression of DAP in chicken satellite cells resulted in increased myotube diameter and higher rates of proliferation and differentiation, further suggesting a role for DAP in muscle development based on results from this cell culture model.
The current study was designed to elucidate the molecular mechanisms by which DAP mediates skeletal muscle growth and development. Transcriptomes of DAP-knockdown turkey pectoralis major muscle satellite cells were profiled during proliferating and differentiating states against those of non-treated control cells. To our knowledge, these results provide the first evidence that reduced DAP expression dramatically alters global gene expression in both proliferating and differentiating turkey satellite cells.
Knockdown of DAP in proliferating satellite cells, as demonstrated by microarray analysis, showed expression of this gene was reduced approximately three-fold compared to the control cells. In complementary experiments, an approximate five-fold reduction of DAP was observed by qPCR (Figure 3). However, in the differentiating cells, no significant difference in DAP abundance was observed between treated cells and controls by either microarray or qPCR studies. The absence of DAP differential expression during differentiation was likely the result of the transient transfection system used for these experiments. As described previously (), satellite cells were first transfected with siRNA, grown to 65% confluence (about 72 h), and then induced to differentiate over a period of 48 h. The half-life of DAP transcripts remains unknown and the treatment conditions may have provided enough time for the cells to transcribe sufficient levels of DAP such that levels were not significantly different between the knockdown and control cells by the 48 h differentiation timepoint. It is also important to note that changes in transcript abundance due to siRNA transfection do not necessarily result in changes in protein levels.
Proliferation
In proliferating turkey pectoralis major muscle satellite cells, a majority of the affected genes were down-regulated compared to the control. The microarray data were corroborated with qPCR, with all of the examined genes showing similar overall trend in fold change. Pathway analysis (Table 2) suggested that knockdown of DAP in proliferating turkey satellite cells altered the cellular signals regulating post-transcriptional processes, in large part, through the modification of eIFs and protein turnover. The deviation of such canonical signals would affect cellular activities, particularly protein synthesis, and eventually in extreme cases, lead to cell death (Table 3).
The mTOR signaling pathway was one of the most significantly affected canonical pathways in response to DAP knockdown. The mTOR enzyme is a serine-threonine protein kinase that regulates cell growth, proliferation, and development in response to nutritional and environmental cues (; ). This pathway is regulated in part by DAP which is a substrate for phosphorylation (). Studies suggest that mTOR is important in regulating the ultimate growth and size of skeletal muscle (; ; ). This is due, in large part, to control of the underlying mechanisms of protein synthesis (; Zanchi and Lancha, 2008; ).
In this study, several eIFs, essential for initiation of eukaryotic translation (i.e., protein synthesis), showed altered gene expression as a result of DAP knockdown in proliferating satellite cells (Supplementary Table S2 and Figure 4). Fold change of nine initiation factors affected by DAP knockdown (eIF2S1, eIF3A, eIF3D, eIF3E, eIF3L, eIF3J, eIF4E, eIF4G2, and eIF5) ranged from 0.22 to 29.22. Altered expression of these initiation factors can affect satellite cell proliferation and differentiation (; ). For example, eIF4E was down-regulated more than fourfold (FDR < 0.05) in proliferating turkey satellite cells upon DAP knockdown. The translational protein eIF4E recognizes and binds to the 5′cap of mRNA, aiding in the recruitment of mRNA to the ribosomal complex (). A second eIF (eIF4G2), also commonly known as p97 or DAP5, was down-regulated nearly fourfold (FDR < 0.05). Reduced expression of eIF4G2 by siRNA in HEK293 cells led to decreased overall translation rate and inhibition of cell proliferation (). Other eIFs affected by DAP knockdown include the eIF3s that are necessary for formation of the pre-initiation complex with the 40S ribosomal subunit during translational initiation ().
In the present study, differential expression of several eiF3 subunits was observed with DAP knockdown (FDR < 0.05). One subunit (eIF3D) was down-regulated more than twofold while, the other four, including eIF3A, eIF3E, eIF3J, and eIF3L, were up-regulated. Among these four, eIF3E was the most DE in the proliferating cells with FC > 29-fold. Differential expression changes of eIF3A and eIF3L were comparable (∼fivefold) whereas eIF3J was increased more than ninefold. Because all other eIFs showed a decrease in expression upon knockdown of DAP, it is proposed that expression of the up-regulated eIF3s may be a compensatory mechanism for loss of translational control at early stages during the translation initiation process. Consistent with our findings, yeast cells with reduced expression of eIFs display altered cell cycle profile with an accumulation of the cells in G1 phase (Yu et al., 2006).
Another critical gene significantly altered by DAP knockdown during proliferation was AKT1, down-regulated >4-fold relative to control cells (FDR < 0.05). This gene encodes the AKT1 protein (also termed protein kinase B) that plays important roles in cell metabolism, proliferation, and growth. AKT1 exerts an anti-apoptotic role in cell survival, and mediates cell growth through the mTOR pathway (; ). In mice, targeted knockdown of AKT1 in C2 myoblasts blocked differentiation, although with no apparent effect on cell proliferation (). Conversely, found AKT1 to play a critical role in proliferation of C2.7 myoblasts and 3T3 fibroblasts. showed that Akt-knockout mice had reduced muscle mass, grip strength, and contractile force. Thus, our results indicating reduced AKT1 expression associated with the DAP knockdown are consistent with impaired proliferation and development. In addition, as noted by AKT1 inhibits autophagosome formation and lysosome-directed protein degradation. Thus, it is reasonable to hypothesize that the degree of suppression of AKT1 expression in the DAP knockdown cells may result in increased degradation of muscle proteins due to the absence of the inhibitory effect of AKT1.
The current findings highlight the impaired protein biosynthesis, particularly via modified translation initiation and mTOR signaling, in the DAP-knockdown satellite cells. Such molecular alteration might in turn restrict progression of the cells to initiate a new round of cell division, hence diminishing cell proliferation in the knockdown cells ().
Differentiation
Based on pathway analysis, the calcium signaling pathway was the top altered canonical pathway reported by IPA (Table 2 and Figure 6). A total of 13 genes in this pathway displayed differential expression upon DAP knockdown with expression ranging from 34-fold decrease to 2.5-fold increase. Proteins in the calcium signaling pathway include SR calcium regulators that are responsible for calcium uptake and release, and myofibrillar proteins that respond to changes in [Ca2+] to effect muscle contraction. Of the SR calcium regulators, the voltage dependent calcium channel subunit gamma 1 (CACNG1) and both isoforms of avian ryanodine receptors, i.e., RYR1 and RYR3 were down-regulated approximately 28-fold, 2-fold, and 6-fold, respectively. The CACNG1 gene encodes the dihydropyridine receptor (DHPR) gamma subunit, one of the five subunits of the L-type calcium channel (). The DHPR is located in the transverse tubule and associates with RYR1 to elicit calcium release from the SR into the sarcoplasm (; ). The gamma subunit of the DHPR, in particular, may act as a calcium antagonist during times of cell stress (). Both RYR1 and RYR3 are calcium channel proteins that open upon activation to allow diffusion of Ca2+ from SR lumen to the cytosol. The sarco/endoplasmic reticulum-Ca2+ ATPase is responsible for reducing cytosolic Ca2+ by active transport, coupled with the hydrolysis of ATP, back into the SR (). The SERCA gene ATP2A2 showed nearly threefold down-regulation (FDR < 0.05) upon DAP knockdown. Although SERCA2 is more commonly associated with slow-twitch and cardiac muscle, SERCA2a, the more prominent of the SERCA2 isoforms, has been detected in fast-twitch skeletal muscle and neonatal muscle in humans (Zarain-Herzberg and Alvarez-Fernandez, 2002). Proteins that comprise most of the mass of the muscle myofibril were also downregulated, including the structural proteins titin and nebulin, the contractile proteins myosin and actin, and the Ca2+-sensitive regulatory proteins tropomyosin and the three troponin subunits. In addition, obscurin, which may play a role in initiating myofibril assembly (), was also downregulated, potentially reducing myotube development observed in the knockdown experiments ().
It is widely recognized that Ca2+ acts as a second messenger responsible for regulating cascades of signal transduction. Differential expression of the crucial Ca2+-regulatory proteins responsible for release and uptake of cytosolic Ca2 in differentiating satellite cells suggests aberrant cytosolic [Ca2+] homeostasis with profound downstream consequences for cellular development. In myogenesis, precise control of [Ca2+] is required for activating muscle-specific transcription factors that regulate the fusion of myoblasts to form multinucleated myotubes (). Dysregulated [Ca2+] in the DAP-knockdown cells may interfere with numerous cellular processes including the calcium-mediated pathways that are critical for myotube formation.
The other top altered canonical pathways (Table 2) function in a coordinated manner to regulate cytoskeletal reorganization largely required for cell differentiation. The cytoskeleton, recognized as a highly dynamic scaffold, comprises different molecules, e.g., actin and tubulin, arranged into filaments distributed throughout the cytoplasm. Reorganization of such structures induced by extracellular stimuli, activates intracellular signaling, including Rho-family small GTPases through the links with extracellular matrix adhesion molecules such as integrins and syndecans (; ; ). In this study, two main cytoskeletal proteins, including actin and tubulin, were substantially down-regulated. Within the DAP-knockdown cells, decreased transcripts include two loci corresponding to alpha actin (ACTA1, FC = −19 and −12), cardiac actin (ACTC1, FC = −16) and three loci of alpha tubulin (LOC100545462, FC ranging from −5.5 to −9.2). In addition, expression of integrin beta-1 (ITGB1), an integrin complex subunit, and filamin (FLN2), that serves as a linking protein between two actin filaments, were decreased 4.9-fold and 3.5-fold, respectively. In contrast, syndecan-2 (SDC2) increased 4.3-fold. Two loci of the non-regulatory myosin light chain (MYL1) were the most highly down-regulated at >143-fold and >59-fold. Structurally and functionally different from myosin heavy chain, myosin light chain resides in the neck region of myosin and interacts with Rho A effectors, hence mediating assembly of actomyosin complex. With altered cytoskeleton signaling, the formation of multinucleated myotubes could be disrupted, thus preventing the differentiation of the knockdown cells.
There were 80 DE loci common to both proliferation and differentiation. Of the top five canonical pathways affected (Table 4), the most dramatically affected pathway, Nitric Oxide Signaling in the Cardiovascular System, is of particular interest. The DE genes within this pathway include upregulation of the heat shock proteins (HSP90AA1 and HSP90B), the catalytic subunit of AMPK (PRKAA1), and downregulation of AKT1. It is well established that AKT1 plays a critical role in regulating proliferation, differentiation, size, and viability of muscle cells in part through changes in phosphorylation status (Yun and Matts, 2005). The protein HSP90 plays multiple roles in the developing muscle cell. Like most heat shock proteins, it is involved in protein folding and in refolding of denatured proteins (). In addition, HSP90 binds to AKT1 and maintains a balanced phosphorylation state of AKT1 by reducing dephosphorylation through protein phosphatase 2A (PP2A) (; Yun and Matts, 2005). The substantial downregulation of AKT1 coupled with upregulation of HSP90AA1 and HSP90B could modulate this balance. Moreover, HSP90 is a target of S-nitrosylation that results in loss of its ATP-dependent protein folding activity (), and could result in accumulation of misfolded proteins accelerating autophagy or apoptosis. Finally, the upregulation of PRKAA1 is associated with decreased protein synthesis through phosphorylation of mTOR and tuberous sclerosis complex 2 (TSC2), thereby increasing autophagy (). AMPK also phosphorylates key regulators of autophagy including ULK1, ATG9, and beclin-1.
Early studies of DAP in HeLa cells suggest it is an important positive regulator of cytokine-mediated apoptosis (). Subsequently, showed that DAP is phosphorylated under nutrient-rich conditions but functionally silent under starvation conditions where DAP is dephosphorylated and acts as a suppressor of autophagy. The process of autophagy is implicated in skeletal muscle development and renewal as a method for efficiently recycling and removing proteins from cells (). Ca2+, acting as a second messenger, regulates the processes of autophagy and apoptosis through complex spatiotemporal mechanisms that may either promote or inhibit these pathways (; ). Increased or decreased intracellular [Ca2+] may disrupt regulation of these pathways, leading to altered cellular development. The results of the current study suggest that DAP plays a central role in turkey muscle growth and development, and perhaps in other cell types, by regulating autophagy as well as calcium signaling and protein translation.
In conclusion, this study demonstrated global changes in the transcriptome following knockdown of DAP in developing turkey satellite cells. The findings underline the essential roles of DAP in proliferation, potentially through the transduction signals associated with translation initiation and mTOR signaling. During differentiation, DAP knockdown altered the biological pathways responsible for controlling [Ca2+] and cytoskeleton organization. Future studies will focus on mechanisms by which changes in DAP expression, and consequent changes in the level of cellular DAP protein, alter cellular homeostasis.
Statements
Data availability statement
Raw Cy5 and Cy3 intensities, M-value, and LOESS-normalized log2 average intensities (A) were submitted with original GPR files to the NCBI GEO (platform accession: GPL9788, series accession: GSE35660).
Author contributions
GS, SV, and KR designed the experiments and acquired the funding for this project. SV performed the cell culture work including DAP knockdown experiments. KH and KS isolated RNA and performed the microarray and qPCR experiments. RT analyzed the raw microarray data. KR annotated the microarray oligonucleotides and conducted the IPA analysis. KH made the initial draft of the manuscript with major revisions provided by YM, KR, SV, and GS. All authors read, revised, and approved the final manuscript.
Funding
This project was supported by a grant from USDA National Institute of Food and Agriculture Award number 2005-35604-15628 to GS, SV, KR, and RT.
Acknowledgments
We would like to acknowledge Cindy Coy for technical assistance on the in vitro satellite cell culture work and Juan Abrahante for bioinformatic assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.01036/full#supplementary-material
References
1
AitkenC. E.LorschJ. R. (2012). A mechanistic overview of translation initiation in eukaryotes.Nat. Struct. Mol. Biol.19568–576. 10.1038/nsmb.2303
2
AmbrizX.de LanerolleP.AmbrosioJ. R. (2018). The mechanobiology of the actin cytoskeleton in stem cells during differentiation and interaction with biomaterials.Stem Cells Int.2018:2891957. 10.1155/2018/2891957
3
AndronacheZ.UrsuD.LehnertS.FreichelM.FlockerziV.MelzerW. (2007). The auxiliary subunit gamma 1 of the skeletal muscle L-type Ca2+ channel is an endogenous Ca2+ antagonist.Proc. Natl. Acad Sci. U.S.A.10417885–17890. 10.1073/pnas.0704340104
4
BodineS. C.StittT. N.GonzalezM.KlineW. O.StoverG. L.BauerleinR.et al (2001). Akt/mTOR pathway is a crucial regulator of skeletal muscle hypertrophy and can prevent muscle atrophy in vivo.Nat. Cell. Biol.31014–1019. 10.1038/ncb1101-1014
5
BuckinghamM.BajardL.ChangT.DaubasP.HadchouelJ.MeilhacS.et al (2003). The formation of skeletal muscle: from somite to limb.J. Anat.20259–68. 10.1046/j.1469-7580.2003.00139.x
6
CarlssonL.YuJ. G.ThornellL. E. (2008). New aspects of obscurin in human striated muscles.Histochem. Cell. Biol.13091–103. 10.1007/s00418-008-0413-z
7
DeissL. P.FeinsteinE.BerissiH.CohenO.KimchiA. (1995). Identification of a novel serine/threonine kinase and a novel 15-kD protein as potential mediators of the γ interferon-induced cell death.Genes Dev.915–30. 10.1101/gad.9.1.15
8
FloriniJ. R. (1989). Assay of creatine kinase in microtiter plates using thio-NAD to allow monitoring at 405 nm.Anal. Biochem.182399–404. 10.1016/0003-2697(89)90614-3
9
FuX.WangH.HuP. (2015). Stem cell activation in skeletal muscle regeneration.Cell. Mol. Life Sci.721663–1677. 10.1007/s00018-014-1819-5
10
GoncalvesM. D.PistilliE. E.BalduzziA.BirnbaumM. J.LacheyJ.KhuranaT. S.et al (2010). Akt deficiency attenuates muscle size and function but not the response to ActRIIB inhibition.PLoS One.5:e12707. 10.1371/journal.pone.0012707
11
HarrM. W.DistelhorstC. W. (2010). Apoptosis and autophagy: decoding calcium signals that mediate life or death.Cold Spring Harb. Perspect. Biol.2:a005579. 10.1101/cshperspect.a005579
12
HayN.SonenbergN. (2004). Upstream and downstream of mTOR.Genes Dev.181926–1945. 10.1101/gad.1212704
13
Héron-MilhavetL.FranckhauserC.RanaV.BerthenetC.FisherD.HemmingsB. A.et al (2006). Only Akt1 is required for proliferation while Akt2 promotes cell cycle exit through p21 binding.Mol. Cell Biol.268267–8280. 10.1128/MCB.00201-06
14
HuangC. L.PedersenT. H.FraserJ. A. (2011). Reciprocal dihydropyridine and ryanodine receptor interactions in skeletal muscle activation.J. Muscle Res. Cell Motil.32171–202. 10.1007/s10974-011-9262-9
15
JiaY.YeL.JiK.TomsA. M.DaviesM. L.RugeF.et al (2014). Death associated protein 1 is correlated with the clinical outcome of patients with colorectal cancer and has a role in the regulation of cell death.Oncol. Rep.31175–182. 10.3892/or.2013.2866
16
KorenI.ReemE.KimchiA. (2010). DAP1, a novel substrate of mTOR, negatively regulates autophagy.Curr. Biol.201093–1098. 10.1016/j.cub.2010.04.041
17
KrämerA.GreennJ.PollardJ.Jr.TugendreichS. (2014). Causal analysis approaches in ingenuity pathway analysis.Bioinformatics30523–530. 10.1093/bioinformatics/btt703
18
LeeS. H.McCormickF. (2006). p97/DAP5 is a ribosome-associated factor that facilitates protein synthesis and cell proliferation by modulating the synthesis of cell cycle proteins.EMBO J.254008–4019. 10.1038/sj.emboj.7601268
19
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method.Methods25402–408. 10.1006/meth.2001.1262
20
ManningB. D.CantleyL. C. (2007). AKT/PKB signaling: navigating downstream.Cell1291261–1274. 10.1016/j.cell.2007.06.009
21
Martínez-RuizA.VillanuevaL.González de OrduñaC.López-FerrerD.HiguerasM. A.TarínC.et al (2005). S-nitrosylation of Hsp90 promotes the inhibition of its ATPase and endothelial nitric oxide synthase regulatory activities.Proc. Nat. Acad. Sci. U.S.A.1028525–8530. 10.1073/pnas.0407294102
22
McFarlandD. C.DoumitM. E.MinshallR. D. (1988). The turkey myogenic satellite cell: optimization of in vitro proliferation and differentiation.Tissue Cell20899–908. 10.1016/0040-8166(88)90031-6
23
McFarlandD. C.GilkersonK. K.PesallJ. E.WalkerJ. S.YunY. (1995). Heterogeneity in growth characteristics of satellite cell populations.Cytobios8221–27.
24
NestorK. E. (1977). The stability of two randombred control populations of turkeys.Poult. Sci.5654–57. 10.3382/ps.0560054
25
NestorK. E.McCartneyM. G.BachevN. (1969). Relative contributions of genetics and environment to turkey improvement.Poult. Sci.481944–1949. 10.3382/ps.0481944
26
OhannaM.SoberingA. K.LapointeT.LorenzoL.PraudC.PetroulakisE.et al (2005). Atrophy of S6k1(-/-) skeletal muscle cells reveals distinct mTOR effectors for cell cycle and size control.Nat. Cell Biol.7286–294. 10.1038/ncb1231
27
PallafacchinaG.CalabriaE.SerranoA. L.KalhovdeJ. M.SchiaffinoS. (2002). A protein kinase B-dependent and rapamycin-sensitive pathway controls skeletal muscle growth but not fiber type specification.Proc. Natl. Acad. Sci. U.S.A.999213–9218. 10.1073/pnas.142166599
28
PolymenisM.AramayoR. (2015). Translate to divide: control of the cell cycle by protein synthesis.Microb. Cell294–104. 10.15698/mic2015.04.198
29
ReedK. M.MendozaK. M.JunejaB.FahrenkrugS. C.VellemanS.ChaingW.et al (2008). Characterization of expressed sequence tags from turkey skeletal muscle.Anim. Genet.39635–644. 10.1111/j.1365-2052.2008.01787.x
30
RitchieM. E.SilverJ.OshlackA.HolmesM.DiyagamaD.HollowayA.et al (2007). A comparison of background correction methods for two-colour microarrays.Bioinformatics232700–2707. 10.1093/bioinformatics/btm412
31
RossiA. E.DirksenR. T. (2006). Sarcoplasmic reticulum: the dynamic calcium governor of muscle.Muscle Nerve33715–731. 10.1002/mus.20512
32
RotweinP.WilsonE. M. (2009). Distinct actions of Akt1 and Akt2 in skeletal muscle differentiation.J. Cell Physiol.219503–511. 10.1002/jcp.21692
33
SandriM. (2010). Autophagy in skeletal muscle.FEBS Lett.5841411–1416. 10.1016/j.febslet.2010.01.056
34
SatoS.FujitaN.TsuruoT. (2000). Modulation of Akt kinase activity by binding to Hsp90.Proc. Nat. Acad. Sci. U.S.A.9710832–10837. 10.1073/pnas.170276797
35
ShinJ.McFarlandD. C.StrasburgG. M.VellemanS. G. (2013a). Function of death-associated protein 1 in proliferation, differentiation, and apoptosis of chicken satellite cells.Muscle Nerve48777–790. 10.1002/mus.23832
36
ShinJ.McFarlandD. C.VellemanS. G. (2013b). Migration of turkey muscle satellite cells is enhanced by the syndecan-4 cytoplasmic domain through the activation of RhoA.Mol. Cell. Biochem.375115–130. 10.1007/s11010-012-1534-1
37
SmythG. K. (2005). “LIMMA: linear models for microarray data,” in Bioinformatics and Computational Biology Solutions using R and Bioconductor, edsGentlemanR.CareyV.DudoitS.IrizarryR.HuberW. (New York, NY: Springer), 397–420. 10.1007/0-387-29362-0_23
38
SporerK. R. B.ChiangW.TempelmanR. J.ErnstC. W.ReedK. M.VellemanS. G.et al (2011a). Characterization of a 6K oligonucleotide turkey skeletal muscle microarray.Anim. Genet.4275–82. 10.1111/j.1365-2052.2010.02085.x
39
SporerK. R. B.TempelmanR. J.ErnstC. W.ReedK. M.VellemanS. G.StrasburgG. M. (2011b). Transcriptional profiling identifies differentially expressed genes in developing turkey skeletal muscle.BMC Genomics12:143. 10.1186/1471-2164-12-143
40
SunF.XuX.WangX.ZhangB. (2016). Regulation of autophagy by Ca2+.Tumor Biol.3715467–15476. 10.1007/s13277-016-5353-y
41
ThompsonD. M. (2018). The role of AMPK in the regulation of skeletal muscle size, hypertrophy, and regeneration.Int. J. Mol. Sci.19:3125. 10.3390/ijms19103125
42
VellemanS. G.LiuC.CoyC. S.McFarlandD. C. (2006). Effects of glypican-1 on turkey skeletal muscle cell proliferation, differentiation and fibroblast growth factor 2 responsiveness.Dev. Growth Differ.48271–276. 10.1111/j.1440-169X.2006.00860.x
43
VellemanS. G.LiuX.NestorK. E.McFarlandD. C. (2000). Heterogeneity in growth and differentiation characteristics in male and female satellite cells isolated from turkey lines with different growth rates.Comp. Biochem. Physiol. Pt. A Mol. Integr. Physiol.125503–509. 10.1016/s1095-6433(00)00178-1
44
VellemanS. G.SporerK. R. B.ErnstC. W.ReedK. M.StrasburgG. M. (2012). Versican, matrix Gla protein, and death-associated protein expression affect muscle satellite cell proliferation and differentiation.Poult. Sci.911964–1973. 10.3382/ps.2012-02147
45
WazirU.JiangW. G.SharmaA. K.MokbelK. (2012). The mRNA expression of DAP1 in human breast cancer: correlation with clinicopathological parameters.Cancer Genom. Proteom.9199–201.
46
WullschlegerS.LoewithR.HallM. N. (2006). TOR signaling in growth and metabolism.Cell124471–484. 10.1016/j.cell.2006.01.016
47
YahiroK.TsutsukiH.OguraK.NagasawaS.MossJ.NodaM. (2014). DAP1, a negative regulator of autophagy, controls SubAB-mediated apoptosis and autophagy.Infect. Immun.824899–4908. 10.1128/IAI.02213-14
48
YenC. F.WangH. S.LeeC. L.LiaoS. K. (2014). Roles of integrin-linked kinase in cell signaling and its perspectives as a therapeutic target.Gynecol. Minimally Invasive Ther.367–72. 10.1016/j.gmit.2014.06.002
49
YuL.CastilloL. P.MnaimnehS.HughesT. R.BrownG. W. (2006). A survey of essential gene function in the yeast cell division cycle.Mol. Biol. Cell174736–4747. 10.1091/mbc.e06-04-0368
50
YunB.-G.MattsR. L. (2005). Hsp90 functions to balance the phosphorylation state of Akt during C2C12 myoblast differentiation.Cell. Signal.171477–1485. 10.1016/j.cellsig.2005.03.006
51
YunY.McFarlandD. C.PesallJ. E.GilkersonK. K.Vander WalL. S.FerrinN. H. (1997). Variation in response to growth factor stimuli in satellite cell populations.Comp. Biochem. Physiol. Part A. Physiol.117463–470. 10.1016/s0300-9629(96)00404-5
52
ZanchiN. E.LanchaA. H.Jr. (2008). Mechanical stimuli of skeletal muscle: implications on mTOR/p70s6k and protein synthesis.Eur. J. Appl. Physiol.102253–263. 10.1007/s00421-007-0588-3
53
Zarain-HerzbergA.Alvarez-FernandezG. (2002). Sarcoendo plasmic reticulum Ca2+-ATPase-2 gene: structure and transcriptional regulation of the human gene.Sci. World J.21469–1483. 10.1100/tsw.2002.228
Summary
Keywords
satellite cells, muscle development, death-associated protein, global gene expression, microarray, muscle, turkey
Citation
Horton KA, Sporer KRB, Tempelman RJ, Malila Y, Reed KM, Velleman SG and Strasburg GM (2020) Knockdown of Death-Associated Protein Expression Induces Global Transcriptome Changes in Proliferating and Differentiating Muscle Satellite Cells. Front. Physiol. 11:1036. doi: 10.3389/fphys.2020.01036
Received
31 March 2020
Accepted
28 July 2020
Published
14 August 2020
Volume
11 - 2020
Edited by
Walter Gay Bottje, University of Arkansas, United States
Reviewed by
Paweł Tomasz Maækowiak, Poznań University of Life Sciences, Poland; Carl Joseph Schmidt, University of Delaware, United States
Updates
Copyright
© 2020 Horton, Sporer, Tempelman, Malila, Reed, Velleman and Strasburg.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gale M. Strasburg, stragale@msu.edu
†These authors have contributed equally to this work
This article was submitted to Avian Physiology, a section of the journal Frontiers in Physiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.