Skip to main content

Cloning and promoter analysis of palladin 90-kDa, 140-kDa, and 200-kDa isoforms involved in skeletal muscle cell maturation

Abstract

Objective

Palladin is a ubiquitous phosphoprotein expressed in vertebrate cells that works as a scaffolding protein. Several isoforms deriving from alternative splicing are originated from the palladin gene and involved in mesenchymal and muscle cells formation, maturation, migration, and contraction. Recent studies have linked palladin to the invasive spread of cancer and myogenesis. However, since its discovery, the promoter region of the palladin gene has never been studied. The objective of this study was to predict, identify, and measure the activity of the promoter regions of palladin gene.

Results

By using promoter prediction programs, we successfully identified the transcription start sites for the Palld isoforms and revealed the presence of a variety of transcriptional regulatory elements including TATA box, GATA, MyoD, myogenin, MEF, Nkx2-5, and Tcf3 upstream promoter regions. The transcriptome profiling approach confirmed the active role of predicted transcription factors in the mouse genome. This study complements the missing piece in the characterization of palladin gene and certainly contributes to understanding the complexity and enrollment of palladin regulatory factors in gene transcription.

Introduction

Palladin is a cytoskeletal associated protein that plays a fundamental role in human and animal cell morphology and adhesion. It was simultaneously and independently identified and characterized by two groups of researchers in early 2000 [1, 2]. At that time, the protein is known to organize and stabilize the actin cytoskeleton and focal adhesions at the cell–cell and cell–matrix junctions in embryonic and fibroblast cells and tissues [2]. Subsequently, the mouse palladin gene (Palld) was described as a complex structure that spans approximately 400 kb and is composed of 25 exons coding for several isoforms [3, 4] that highly interact with various cytoskeletal associated proteins [1, 5,6,7,8]. The contribution of palladin to various biological, physiological processes and signaling pathways within cells is not undeniable [9,10,11] and has been partially elucidated in our laboratory.

From online databases, we found that the human palladin gene generates about 14 isoforms against four for the mouse palladin gene. Investigations are still carried out to classify proteins that regulate or directly interact with palladin with a substantial impact on cell structure and function. However, no study has ever identified the activity of promoter regions for the three well-characterized palladin isoforms (90-kDa, 140-kDa, and 200-kDa). Here, we used bioinformatics tools to predict the promoter region of Palld and describe potential palladin regulatory factors.

Main text

Methods

Cell culture

C2C12 cells, a mouse myogenic cell line, were purchased from the American Type Culture Collection (ATCC® CRL-1772™) and maintained in complete growth medium (GM) containing Dulbecco’s modified Eagle’s medium (DMEM, Gibco) with 10% fetal bovine serum (Gibco), and 1% penicillin/streptomycin (Gibco) in a humidified atmosphere of 5% CO2 at 37 °C. Differentiation of myoblasts into myotubes was induced by exchanging GM to differentiation medium (DM) containing DMEM and 2% horse serum (Gibco). Fresh DM was replaced every day after differentiation induction from Day 0 (D0) to Day 5 (D5).

Promoter analysis of mouse gene Palld

DNA sequences of Palld (Gene ID: 72333) were obtained from the National Center for Biotechnology Information (NCBI) and Mouse Genome Informatics (MGI) databases. Identification of the Transcription Start Site (TSS) of the three palladin isoforms 90-kDa, 140-kDa, and 200-kDa was performed using the Ensembl genome browser. DNA sequences of approximately − 3 kb from the start codon of each isoform were analyzed by using three online promoter prediction programs including Promoter 2.0 Prediction Server, Neural Network Promoter Prediction, and PROMOTER SCAN. To compare these gene sequences with multiple vertebrate genomes, an Evolutionary Conserved Regions (ECR) Browser database was used. Prediction of Transcription Factor Binding Sites in DNA sequences was performed by using MATCH and PROMO tools.

DNA construction of luciferase reporter and DNA transfection

Multiple DNA fragments of predicted Palld promoter regions were amplified from mouse genomic DNA (C57BL/6 J) and subcloned into the pGL4.17 [luc2/neo] vector (Promega, USA) with a downstream tagged luciferase gene. DNA fragments of palladin isoform 90-kDa, 140-kDa, and 200-kDa were respectively inserted into SacI and BglII, KpnI and EcoRV, and KpnI and HindIII sites of the multiple cloning region of the vector.

C2C12 cells were transfected with various palladin promoter constructs tagged with a luciferase gene or pGL4.74 [hRluc, TK] empty vector (Promega, USA) as a control plasmid by using Lipofectamine 2000 reagent (Invitrogen) and further incubated for 4 days following differentiation induction.

Luciferase reporter assay

The Dual-Glo Luciferase Reporter Assay System Kit (Promega) was used according to the manufacturer’s protocol. Luminescence of each sample was quantified on a Tecan Infinite F200 Pro ELISA reader (USA). Values were normalized to the control vector, and the results were expressed as the ratio of firefly luciferase activity to Renilla luciferase activity for each construct.

RNA isolation and sequencing

Total RNA from C2C12 cells was extracted by using High Pure RNA Isolation Kit (Roche) according to the manufacturer’s instructions. The RNA concentration was determined using NanoDrop 2000 spectrophotometer (Thermo fisher) and RNA purity and integrity was evaluated using Experion™ Automated Electrophoresis System. The samples were submitted to Yourgene Bioscience (Taiwan) for sequencing. The strand specific Poly-A mRNA libraries were generated and sequenced by using Illumina HiSeq™ 4000 (NovogeneAIT, Singapore) platform and procedure. Quality trimming was performed and paired-end reads were mapped to the mouse reference genome GRCm38.p6 using TopHat software. Assembled transcripts were acquired by using Cufflinks and transcriptome was assembled via Cuffmerge. Gene expression levels were normalized by computing the Fragments Per Kilobase of transcript per Million mapped reads (FPKM). Because of the downstream usage of sequencing results not discussed in the present work, the data have not yet been made available in a public repository. For sequencing data partially used to support our analysis, raw data can be accessed on request.

Quantitative PCR

The same RNA pools used for RNA seq were employed for qPCR. Complementary DNA (cDNA) was synthesized using iScript cDNA Synthesis Kit (Bio-Rad). The primers used were listed in Additional file 1. Each reaction tubes contained Power SYBR Green PCR Master Mix (Applied Biosystems), Forward and Reverse primer at the final concentration of 500 nM, 1 µl of synthetized cDNA, and PCR grade water to a total volume of 10 µl. Each sample was tested in duplicate using the StepOnePlus Real-Time PCR System (Applied Biosystems) at the following conditions: 95 °C for 10 min, 40 cycles of 95 °C for 15 s, and 60 °C for 1 min. Amplification specificity was assessed by melting curve. The Ct values obtained were normalized against the Adaptor Related Protein Complex 3 Subunit Delta 1 (Ap3d1) endogenous control and expression of mRNA calculated using the 2−∆Ct, a derivative of Livak method [12].

Statistical analysis

Two-way repeated-measures analysis of variance (ANOVA) with Bonferroni corrections for multiple group comparisons at various time point were calculated. Values of p ≤ 0.05 were considered statistically significant.

Results

Identification of predictive Palld promoter regions

To predict and identify the nature of regulatory elements controlling the activation or repression of palladin isoforms, gene fragments of 3000 bp upstream the first exon were generated and the palladin gene promoter region analyzed. The promoter regions were predicted to be within the − 2917 bp/+ 83 bp, − 2721 bp/+ 279 bp, and − 2871 bp/+ 129 bp, respectively for palladin isoform 90-kDa, 140-kDa, and 200-kDa. Promoter 2.0 Prediction Server data indicate the presence of at least one putative eukaryotic polymerase II binding site within each isoform (Additional file 2: Appendix A). Through the Neural Network Promoter Prediction program, various polymerase II binding sites were observed at different positions along with TATA box sequences (Additional file 2: Appendix B). While the existence of TATA box-binding sites implies that of a polymerase II binding site, it did not necessarily refer to the exact transcription initiation location. Web Promoter Scan program revealed potential transcription factor binding sites for the 90-kDa isoform but failed for the prediction of 140-kDa and 200-kDa (Additional file 2: Appendix C).

To further identify conserved elements in Palld between mouse and other animal species, DNA sequences were aligned and submitted to ECR browser. A schematic representation of the mouse palladin gene indicated the location of isoform transcripts and the promoter region of each transcript (Fig. 1a). Several conserved regions between mouse and human were identified. There were no similarities between mouse (Mus musculus) and opossum (Monodelphis domestica) for palladin 90-kDa isoform however the mouse shared three conserved regions with the human (Homo sapiens) and six intronic regions with the chimpanzee (Pan troglodytes) (Fig. 1b). Six intronic regions of 140-kDa isoform were conserved between mouse, human, and chimpanzee. The 200-kDa isoform presented five common conserved regions with a unique region shared between mouse and opossum.

Fig. 1
figure1

Promotor prediction region of palladin isoforms. a Schematic presentation of the murine palladin gene structure. Blue closed boxes show translated regions of different splice variants and red elliptical shapes show predictive localization of each splice variant promoter region. b visualization of Conserved Regions (ECRs) within the predicted promoter region of palladin isoforms of different species. Mouse palladin reference gene was compared to palladin genes from Opossum, Human, and Chimpanzee

By using MATCH and PROMO software, the promoter regions were restricted up to 1.5 kb for each isoform. For the targeted constructs, substantial high levels of luciferase activity were observed 4 days after C2C12 cell transfection. Promoter regions were delineated based on promoter activity increase in comparison to the activity of the promoter-free construct. Relative firefly/Renilla luciferase activity was measured for the palladin 90-kDa isoform (Fig. 2a), palladin 140-kDa isoform (Fig. 2b), and palladin 200-kDa isoform (Fig. 2c) with their respective gene representations and landmarks. These outcomes indicate multiple potential TSS for palladin isoforms with sequences located at the region between − 1418 and + 83 bp of exon 12 for the 90-kDa isoform, -300 and − 1 bp upstream exon 3 for the 140-kDa isoform and − 300 to + 129 bp of exon 1 for the 200-kDa isoform.

Fig. 2
figure2

Palld gene predictive promoter regions location. Luciferase reporter activities at distinct differentiation days (D0, D1, D2, D4) of various Palld gene constructs comprising potential promoter region and representative gene regulatory sequences of palladin. a 90-kDa isoform, b 140-kDa isoform and c 200-kDa isoform. Luciferase activity measured in wells without cells was used as control (Ctrl). Gene fragments were designed with positive and negative numbers respectively representing number of bp upstream and downstream the first coding region of each isoform. The variation of Firefly/Renilla luciferase ratio between constructs were not statistically different considering the mean of three repeats

RNA-Seq transcriptome analysis

To investigate mRNA dynamics inherent to myoblast cell ongoing differentiation for the specific palladin isoforms and some myogenic genes, C2C12 RNA transcripts were analyzed. The expression of palladin isoforms was determined according to the base pair count which differentiate multiple predicted and known palladin isoforms. The expression of genes involved in cell maturation differentially and significantly changed during myogenesis. The highest expressing transcript was palladin 90-kDa mRNA compared to the other isoforms (Fig. 3a). The expression of 90-kDa and 200-kDa isoforms increased during cell differentiation, which was opposite to that of the 140-kDa isoform. Genes encoding muscle hypertrophic and sarcomeric proteins, for example myogenin (Fig. 3b) and myosin heavy chain (Fig. 3c) were upregulated while myf5 (Fig. 3d) was downregulated.

Fig. 3
figure3

Palladin and myogenic genes mRNA modulation during C2C12 differentiation. Six samples generated from C2C12 differentiated from D0 to D5 were submitted for RNA-Seq analysis and results presented as Fragments Per Kilobase of transcript per Million mapped reads (FPKM) ± SEM, n = 3 independent biological replicates. a Palladin isoforms, b myogenin, c Myosin heavy chain were significantly upregulated and d Myf5 was downregulated during cell differentiation. (*) p value < 0.05

Additionally, RNA transcripts analysis confirmed the presence and the activity of predicted regulatory elements such as Myocyte enhancer factor–2 (MEF2)-Mef2a, Mef2b, Mef2c, homeobox protein Nkx2-5, and Tcf3 in an inclusive way (Additional file 3). The qPCR validated and correlated with transcriptome findings in terms of mRNA expression trend (Additional file 4).

Discussion

Here, we address the gap of knowledge about the promoter sites driving the expression of three major forms of palladin isoforms. Analysis using Promoter 2.0 Prediction Server, Neural Network Promoter Prediction, and PROMOTER SCAN software disclosed different outcomes that might be caused by different computational algorithms used by each promoter prediction software which primarily relies on specific motifs included in the core promoter [13,14,15].

MATCH and PROMO software elicited multiple binding sites of seven transcription factors including a; TATA box (core sequence 5′-TATAAA-3′), GATA, MyoD, myogenin, MEF, Nkx2-5 and E47 sites. Most of the predicted transcription factors are involved in the maturation process of myoblast cells into myofibers. As reported previously, C2C12 cells myoblast remodeling during differentiation requires activation of myogenic regulatory factors (MRF) for example myogenic termination gene (MYOD), myogenin, MRF4, and myogenic factor 5 (MYF5) [16,17,18]. In studies of the molecular mechanism underlying heart development, Nkx2.5, MEF family, GATA 4, and GATA 6 were found to be the key regulatory elements in cardiac muscle [19,20,21]. MEF2 site is essential for the amplification of Mef2c involved in skeletal muscle development [22]. Additionally, analysis of the promoters of the actin-related group of serum response factor (SRF) with scaffolding or contractile functions has indicated that approximately 43% of the involved genes in the mouse contains TATA box sequences within the core promoter region [23]. In a cascade regulation process induced by the cytokine transforming growth factor beta (TGF-β), E2A immunoglobulin enhancer‐binding factors E12 and E47 regulate palladin cytoskeletal protein [24]. These findings confirmed our results regarding the identification of active transcription regulation elements on targeted Palld gene sequences. Luciferase assay showed various activities of the putative promoter core region for the three isoforms which were highly triggered during differentiation. Yet, discrepancies between palladin 140-kDa mRNA measured and the transcript level in luciferase activity were assumed to be normal since luciferase activity could not evaluate protein absolute concentration. Palld has indeed different promoter regions that drive the expression of each isoform and that carry TSS sites and transcriptional regulator factors. The downstream activity of Palld ultimately modulates cells’ differentiation or myogenesis within skeletal muscle cells supporting mammalian embryogenesis and development. The ability to control the expression of the palladin gene could then regulate biological events in which the protein is involved.

Limitations

Transcriptome analysis of C2C12 cells’ ongoing differentiation revealed the functional role of elements predicted to regulate palladin isoform gene. However, at this point in the investigation, we are unable to experimentally attest the presence of the predictive regulatory elements within the promoter regions of the Palld gene. Thus, these results might be considered as a starting point for a more explorative study.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Abbreviations

DMEM:

Dulbecco’s modified Eagle’s medium

GM:

Growth medium

DM:

Differentiation medium

mRNA:

Messenger Ribonucleic acid

TSS:

Transcription start site

FPKM:

Fragments Per Kilobase of transcript per Million mapped reads

Ap3d1:

Adaptor Related Protein Complex 3 Subunit Delta 1

MEF2:

Myocyte enhancer factor-2

MRF:

Myogenic regulatory factors

MYOD:

Myogenic termination gene

Myf5:

Myogenic factor 5

References

  1. 1.

    Mykkänen OM, Grönholm M, Rönty M, Lalowski M, Salmikangas P, Suila H, et al. Characterization of human palladin, a microfilament-associated protein. Mol Biol Cell. 2001;12(10):3060–73.

    PubMed  PubMed Central  Article  Google Scholar 

  2. 2.

    Parast MM, Otey CA. Characterization of palladin, a novel protein localized to stress fibers and cell adhesions. J Cell Biol. 2000;150(3):643–55.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Otey CA, Rachlin A, Moza M, Arneman D, Carpen O. The palladin/myotilin/myopalladin family of actin-associated scaffolds. Int Rev Cytol. 2005;246:31–58.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Wang HV, Moser M. Comparative expression analysis of the murine palladin isoforms. Dev Dyn. 2008;237(11):3342–51.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Boukhelifa M, Moza M, Johansson T, Rachlin A, Parast M, Huttelmaier S, et al. The proline-rich protein palladin is a binding partner for profilin. FEBS J. 2005;273(1):26–33.

    Article  Google Scholar 

  6. 6.

    Boukhelifa M, Parast MM, Bear JE, Gertler FB, Otey CA. Palladin is a novel binding partner for Ena/VASP family members. Cell Motil Cytoskeleton. 2004;58(1):17–29.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Rachlin AS, Otey CA. Identification of palladin isoforms and characterization of an isoform-specific interaction between Lasp-1 and palladin. J Cell Sci. 2006;119(Pt 6):995–1004.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Otey CA, Carpen O. Alpha-actinin revisited: a fresh look at an old player. Cell Motil Cytoskeleton. 2004;58(2):104–11.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Nguyen NU, Liang VR, Wang HV. Actin-associated protein palladin is required for migration behavior and differentiation potential of C2C12 myoblast cells. Biochem Biophys Res Commun. 2014;452(3):728–33.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Nguyen NU, Wang HV. Dual roles of palladin protein in in vitro myogenesis: inhibition of early induction but promotion of myotube maturation. PLoS ONE. 2015;10(4):e0124762.

    PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Umegaki Y, Brotons AM, Nakanishi Y, Luo Z, Zhang H, Bonni A, et al. Palladin is a neuron-specific translational target of mTOR signaling that regulates axon morphogenesis. J Neurosci. 2018;38(21):4985–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Solovyev V, Umarov R. Prediction of prokaryotic and eukaryotic promoters using convolutional deep learning neural networks. arXiv:161000121 [Q-BioGN]. 2016.

  14. 14.

    Abeel T, Saeys Y, Bonnet E, Rouzé P, Van de Peer Y. Generic eukaryotic core promoter prediction using structural features of DNA. Genome Res. 2008;18(2):310–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Bedoya O, Bustamante S. CNN-promoter, new consensus promoter prediction program based on neural networks. Revista EIA. 2011;8(15):153–64.

    Google Scholar 

  16. 16.

    Jiayong T, Aihua H, Hui Y, Gang J, Guangmang L, Xiaoling C, et al. Damage to the myogenic differentiation of C2C12 cells by heat stress is associated with up-regulation of several selenoproteins. Sci Rep. 2018;8(1):10601.

    Article  Google Scholar 

  17. 17.

    Hernández-Hernández JM, García-González EG, Brun CE, Rudnicki MA. The myogenic regulatory factors, determinants of muscle development, cell identity and regeneration. Semin Cell Dev Biol. 2017;72:10–8.

    PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Ferri P, Barbieri E, Burattini S, Guescini M, D’Emilio A, Biagiotti L, et al. Expression and subcellular localization of myogenic regulatory factors during the differentiation of skeletal muscle C2C12 myoblasts. J Cell Biochem. 2009;108(6):1302–17.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Branda T. Heart development: molecular insights into cardiac specification and early morphogenesis. Dev Biol. 2003;258(1):1–19.

    Article  Google Scholar 

  20. 20.

    Stefanovic S, Barnett P, van Duijvenboden K, Weber D, Gessler M, Christoffels VM. GATA-dependent regulatory switches establish atrioventricular canal specificity during heart development. Nat Commun. 2014;5:3680.

    PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Desjardins CA, Naya FJ. The function of the MEF2 family of transcription factors in cardiac development, cardiogenomics, and direct reprogramming. J Cardiovasc Dev Dis. 2016;3(3):26.

    PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Wang DZ, Valdez MR, McAnally J, Richardson J, Olson EN. The Mef2c gene is a direct transcriptional target of myogenic bHLH and MEF2 proteins during skeletal muscle development. Development. 2001;128(22):4623–33.

    CAS  PubMed  Google Scholar 

  23. 23.

    Xu M, Gonzalez-Hurtado E, Martinez E. Core promoter-specific gene regulation: TATA box selectivity and Initiator-dependent bi-directionality of serum response factor-activated transcription. Biochim Biophys Acta. 2016;1859(4):553–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Baulida J. Epithelial-to-mesenchymal transition transcription factors in cancer-associated fibroblasts. Mol Oncol. 2017;11(7):847–59.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

The authors deeply appreciate Dr. Michael W. Hughes for the scientific and language editing as well as Dr. Ngoc-Uyen-Nhi Nguyen and Dr. Xiao-Lei Jin revisions.

Funding

No external assistance has been received.

Author information

Affiliations

Authors

Contributions

All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by BEFO, TYL, CYL, PYC, and CLH. The draft of the manuscript was written by BEFO and was substantively revised by CLH, CCL, YCC, and HVW. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yu-Chung Chiang or Hao-Ven Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Oligonucleotides used for qPCR.

Additional file 2: Table S2.

Transcription start sites (TSS) prediction by Promoter 2.0 Prediction Server.

Additional file 3: Figure S1.

Transcription factors general expression.

Additional file 4: Figure S2.

Expression of palladin isoforms and myogenic genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ouali, B.E.F., Liu, T., Lu, C. et al. Cloning and promoter analysis of palladin 90-kDa, 140-kDa, and 200-kDa isoforms involved in skeletal muscle cell maturation. BMC Res Notes 13, 321 (2020). https://doi.org/10.1186/s13104-020-05152-9

Download citation

Keywords

  • Palladin
  • Cytoskeleton
  • Promoter analysis
  • RNA isoform