- Research article
- Open Access
The SNPs in myoD gene from normal muscle developing individuals have no effect on muscle mass
BMC Geneticsvolume 20, Article number: 72 (2019)
Myogenic Differentiation 1 (MyoD) is a crucial master switch in regulating muscle-specific gene transcription. Forced expression of myoD is equipped to induce several cell lineages into myoblast, which then differentiate and fuse into myotube. Pig is one of the most significant livestock supplying meat, and has been classified into lean, fat and miniature pig breeds. However, the mechanisms underlying muscle mass variation among different pig breeds have remained unclear. Considering the important effect of MyoD on muscle development, it remains to be investigated whether the difference in muscle mass is caused by its single nucleotide polymorphisms (SNPs) which are the major differences among pig breeds at DNA level.
In this study, we identified the locations of porcine myoD regulatory regions including proximal regulatory region (PRR), distal regulatory region (DRR), and core enhancer (CE) region. There are 8 SNPs in the regulatory regions and 6 SNPs in gene body region, which were identified from lean, fat and miniature pig populations. However, these SNPs have no effects on its temporal expression and transcriptional activity which might lead to the distinction in postnatal muscle mass. In addition, overexpression of myoD clones across from amphibious to mammals including xenopus tropicalis, chicken, mouse and pig whose gene identities vary from 68 to 84%, could promote myogenesis in NIH3T3 fibroblasts cells.
These results proved that myoD nucleotide variations from different pig populations have no effect on muscle mass, suggesting that the function of myoD is highly conserved not only among different pig breeds, but also across different species. Thus, it would be futile to discover SNPs affecting muscle mass in pig populations with normal muscle development.
Different pig breeds vary in muscle mass because of genetic differences and intensive selection. Development of skeletal muscle determines muscle mass. Numerous researches on mechanisms of skeletal muscle development have been reported, and many genes regulating myogenesis have been found [1,2,3]. Among them, MyoD, containing basic helix-loop-helix (bHLH) domain, can bind to muscle-specific genes and promote their expression robustly [4, 5]. Previous study revealed that earlier myoD expression exhibits precocious myogenic differentiation and causes severe muscle hypotrophy . In addition, our earlier research found that myoD showed differential expression at 35 days-post-coitus in Landrace and Lantang pigs . These suggest the expression time and level of MyoD have an important effect on muscle mass. However, there are many factors affecting myoD expression, the primary one is its regulatory elements. So far, three regulatory regions have been identified to regulate myoD expression: proximal regulatory region (PRR), distal regulatory region (DRR), and core enhancer (CE) region. The CE region and DRR are essential to regulate myoD expression [8, 9]. Recent findings have found that myoD transcripts corresponding to CE and DRR enhancers can promote chromatin accessibility and RNA polymerase II recruitment at myoD and Myogenin loci, respectively [10,11,12]. In this case, is the difference in muscle mass among different pig breeds related to the SNPs in myoD gene? Therefore, we obtained myoD SNPs in its regulatory regions and gene body by resequencing, and analyzed them to explore whether these SNPs caused alteration in muscle mass in pigs. This study not only helps us to better understand the mechanisms of different muscle mass among pigs, but also provides new clues to study these mechanisms in the future.
Determination of myoD regulatory regions in pigs
The location of myoD regulatory regions in pigs have not been identified before. Depending on the CE, DRR and PRR of human and mouse myoD reported in documents [13,14,15], we finally identified the locations of porcine myoD CE regions (Fig. 1a), DRR (Fig. 1b) and PRR (Fig. 1c) by comparing nucleotide sequence similarities using BLAST. The CE region of porcine myoD was localized to a 258 bp fragment approximately 21,924 to 22,182 bp upstream of myoD transcriptional start site (TSS), which shows 94% sequence homology with human (Additional file 3: Table S3, Additional file 4: Table S4.). The DRR was located at 5 kb upstream relative to the TSS and the PRR was located between -197 bp and -485 bp (Additional file 4: Table S4).
In these regions, it was found that the CE region contains three E-boxes (Fig. 1a, d). In addition, the DRR contains consensus sequences for three E-boxes, two MEF2 binding sites and so on (Fig. 1b, d), which are nearly identical to human and mice [13, 15, 16]. The PRR contains consensus sequences for an E-box, two SP1 sites, an AP2 binding sites, etc. (Fig. 1c, d), which are necessary elements for mouse PRR to regulate muscle-specific transcription [15, 16].
The SNPs in regulatory regions of myoD gene have no effect on muscle mass
Considering the significant role of myoD regulatory region in itself expression, we hypothesized that the polymorphisms of myoD are associated with its expression time and expression level resulting in the divergence of pig muscle mass. To find out whether the SNPs in myoD gene affect muscle mass, whole myoD genome sequence including regulatory regions and gene body were obtained from resequencing results. One SNP existed in all breeds in the CE region (C41446175T) listed in Table 1. However, this SNP also existed within the pig breed. Therefore, this SNP was unlikely to be the major cause that leads to the different muscle mass. Another SNP at locus 41,446,192 in the CE region was found only in individuals, but no breed specificity. The same results were also observed in DRR and PRR (Table 1). Moreover, these SNPs in the three regulatory regions did not appear in their binding motifs which were marked in black boxes (Fig. 1a, b and c). Therefore, these results indicate that SNPs in myoD regulatory region obtained from normally developed pig populations do not contribute to muscle mass.
The mutation (Arg76Pro) has no effect on muscle mass
Six SNPs in myoD coding region (CDS) were discovered through whole-genome resequencing (Table 2). Among them, four SNPs do not result in amino acids alterations, whereas the other two SNPs at locus 41,424,009 (A → G) and 41,424,010 (C → G) synchronously encoded an amino acid, which resulted in the substitution of amino acid from arginine (Arg) to proline (Pro) at 76th amino acid of MyoD protein sequence (Fig. 2a). Of note, the substitution is mainly led by the alteration at locus 41,424,010 from C to G, which was appeared only in the lean pigs.
To determine whether this mutation can affect MyoD expression and then result in distinction in muscle mass, we constructed a myoD expressing plasmid with a C-terminal FLAG tag (PCDNA3.1-flag-MyoDWT) and a mutagenized plasmid to substitute Arg76 with Pro (PCDNA3.1-flag -MyoDR76P). These plasmids were separately transfected into NIH3T3 cells which cannot differentiate into myotubes because of lacking MyoD expression. When MyoD is forced to be expressed, NIH3T3 cell can differentiate into myotubes. So, it is considered as an ideal cell model for testing exogenous MyoD function. As a result, MyoD expression was up-regulated significantly both at mRNA level and protein level (Fig. 2b, c), which remarkably increased the mRNA levels of Ckm, Cdh15 and Myh3 (Fig. 2d). Actually, Ckm, Cdh15 and Myog are the targets of MyoD and promote myogenic differentiation. However, there was no significant difference between MyoDWT and MyoDR76P plasmid in the expression levels of Ckm, Cdh15 and Myh3 (Fig. 2d). MyoDWT and MyoDR76P plasmids were respectively co-transfected into 293 T cells with a 4Rtk-luc reporter plasmid which was used as MyoD-responsive reporter  (Fig. 2e). The transfection of MyoDWT and MyoDR76P both could activate 4Rtk-luc transcription whereas there was no difference between them (Fig. 2e). Immunofluorescence assay showed that overexpression of MyoDWT or MyoDR76P in NIH3T3 cells equally promoted MyHC expression and generated myotubes (Fig. 2f, g). Collectively, these data made it clear that the mutant protein has neither effect on MyoD transcriptional activity in vitro, nor effect on myogenesis.
MyoD clones across from amphibious to mammals can promote myogenesis
Based on the above-mentioned studies, we demonstrated that these SNPs in CDS and regulatory regions of myoD do not result in alteration of its expression and function. It should be noted that none of these mutations occurred in bHLH domain which presented high conservation among these pig populations. This phenomenon makes us have to believe that the bHLH domain of MyoD is conserved in function among normally growing animals, even different species.
To validate the hypothesis, FLAG-tagged myoD plasmids were constructed using pig, mouse, chicken and xenopus tropicalis, respectively. The sequence homologies of MyoD protein among these four species vary from 66 to 88% (Table 3, Additional file 5: Figure S1.), whereas their bHLH regions are highly-conserved (98.08% identity) (Fig. 3a). When these myoD plasmids were transfected into NIH3T3 cells, the mRNA levels of Ckm, Cdh15 and Myh3 were all significantly increased (Fig. 3d, e). However, from the result, it was indicated that there was no significant disparity among the four species at the expression levels of Ckm, and Myh3 regulated by MyoD (Fig. 3f, g). However, the expression level of Cdh15 in mice was significantly higher than that in other three species (Fig. 3f), Which may be related to the fact that NIH3T3 is a mouse-derived cells. In order to assess the transcriptional activity of MyoD from different species, myoD plasmids was co-transfected into 293 T cells with 4Rtk-luc reporter plasmid. As expected, there is no difference of luciferase activity among them (Fig. 3h). Immunofluorescence assay also showed that overexpression of myoD clones promoted MyHC expression and resulted in the similar phenotype in NIH3T3 cells (Fig. 3i).
Together, these results indicated that all myoD clones across from amphibious to mammals can promote myogenesis of NIH3T3 fibroblasts cells owing to the highly-conserved bHLH domain, even though the rest amino acid residues are different.
The ongoing studies about swine have received extensive attention owning to its association to agricultural, clinical and dietary needs. Also, as one of the most significant livestock worldwide, pigs have been classified into lean, fat and miniature categories according to the differences in muscle and fat proportion, and it is well known that there are genetic differences among these three types of pig breeds. Lean pigs, such as Landrace, Duroc and Pietrain, are characterized by high lean meat percentage and fast-growing muscle [18,19,20]. In contrast, fat pigs including Lantang, Luchuan and Laiwu Blacks, as China indigenous pig breeds, are recognized by high subcutaneous fat content, slow-growing muscle as well as low lean percentage [18, 21]. As miniature pigs, Bamaxiang and Wuzhishan are not only relative slow in growth rate and feed conversion rate, but also less in muscle mass . However, the genetic mechanism underlying the difference of muscle mass has remained unclear.
The distinct muscle mass is not relevant to the SNPs in myoD gene
All SNPs in regulatory regions of myoD gene from 10 pig populations were listed in Tables 1 and 2. After analysis, it was inferred that the initiate up-regulated time and level of myoD expression should be not affected by the SNPs of myoD gene. The initiate expression of myoD occurred in myotome in mice and had crucial role in the specification of progenitor cells [3, 23]. Myotube formation was delayed in MyoD-deficient skeletal muscle [1, 24]. In particular, loss of MyoD facilitated adipogenic trans-differentiation of myoblasts . So importantly myoD acts, it is not workable to expect myoD SNPs cause the alteration of myoD expression characteristics in normal muscle developing individuals. The mutations that affect muscle mass should generally exist in diseased individuals or individuals with phenotypic abnormalities. Although a novel variation of amino acid (Arg76Pro) was found in MyoD CDS, and only appeared in lean pigs (Table 2), this mutation did not alter the characteristics of myoD expression (Fig. 2a). Besides, this mutation also occurred in different species including pig, mouse, chicken and xenopus tropicalis (Additional file 5: Figure S1) whose MyoD transcriptional activity were similar in our experiment (Fig. 3). Therefore, this mutation (Arg76Pro) does not affect the transcriptional activity of MyoD, which was verified in NIH3T3 cells (Fig. 2). These data explained detailly that the distinct muscle mass among pig breeds is not relevant to the myoD SNPs which present in different pig populations with normal muscle development.
No mutation occurs in the bHLH domain of MyoD in animals with normal muscle development
MyoD protein contains a bHLH motif, which can form dimers with the E-protein to promote MyoD transcriptional activity . Previous studies have revealed that the bHLH domain of MyoD is sufficient to convert C3H10T1/2 cells to myoblasts . What’s more, when the bHLH domain of MyoD is replaced by the bHLH domain of NeuroD2, MyoD will become a vital monitor of neurogenesis, other than a master regulator of myogenesis . These results suggest that the bHLH domain of MyoD confers its own lineage determination potential. In addition, Myf5, a bHLH protein highly related with MyoD, is also expressed in skeletal muscle and has a crucial role in muscle cell specification. The mutation (p.Arg95Cys) located at the bHLH domain of Myf5 can impair Myf5 nuclear localization and transcriptional activity, which leads to external ophthalmoplegia, rib, and vertebral anomalies in humans . Notably, the bHLH domain of MyoD is highly-conserved among pig breeds, even species in our study, indicating its importance in transcriptional activation. Since this domain is very important in myogenesis, it is unrealistic to expect functional mutations in this domain for normally developing individuals. In addition, it also indicates that difference in muscle mass among breeds is not triggered by the SNPs in bHLH domain, which was confirmed by our sequencing results.
The prospect for exploring the mechanisms of different muscle mass
The difference in muscle mass among pig breeds is not due to the SNPs in myoD gene, which makes us more strongly believe that it is futile to explore the causes of the differences in muscle mass only based on SNPs in populations with normal muscle development. In general, if there is such a mutation that affects muscle mass, there will be abnormal muscle development, such as double muscle rump. After exclusion of SNPs at DNA level, the probable reason resulting in the difference of muscle mass could be attributed to epigenetic modifications that can also be inherited [9, 30, 31]. These modifications may affect the expression of MyoD by affecting its upstream regulatory factors and then influence the differentiation time of myoblasts, eventually leading to differences in muscle mass.
Overall, our study excludes a SNP-based strategy for exploring the mechanisms of different muscle mass among pigs in myoD or other important myogenic genes. At the same time, it also opens another door for researchers.
In summary, our study has showed that the distinction of muscle mass among pig breeds is not caused by the SNPs in myoD gene, but might be attributed by epigenetic modifications, which conduce to understand the mechanisms underlying muscle mass among different pig breeds.
Sample collection and preparation
The pig populations used for SNPs in myoD gene were Bamaxiang (6, miniature), Wuzhishan (6, miniature), Duroc (6, lean), Landrace (5, lean), Large white (6, lean), Pietrain (6, lean), Guangdong Small-ear Spotted (3, fat), Lantang (3, fat), Luchuan (6, fat) and Laiwu Black (6, fat) pigs. The porcine ear tissues were snap-frozen in liquid nitrogen and stored until further use.
DNA library construction and whole-genome resequencing
Porcine genome DNA was obtained using standard phenol chloroform method. One DNA library was constructed for each sample (Guangdong Small-ear Spotted and Lantang pigs) and then sequenced on an Illumina HiSeq X Ten system by Beijing Novogene (commercial service). We also downloaded other whole genome resequencing data of Bamaxiang, Wuzhishan, Duroc, Landrace, Large white, Pietrain, Luchuan and Laiwu Black pigs from GenBank. All the information of these samples was listed in Additional file 1: Table S1. Sequencing reads were mapped with porcine myoD in Sscrofa11.1 version by BWA-MEM  with default parameters and the polymorphisms of myoD were analyzed by GATK3.7 . The sequencing depth of whole genome resequencing for each pig breed was listed in Table 4.
Cell culture and differentiation
Mouse NIH3T3 cell line was purchased from ATCC, and were cultured in high-glucose DMEM with 10% fetal bovine serum (growth medium, GM) until confluence. NIH3T3 cells were switched into DMEM with 2% horse serum (differentiation medium, DM) when cells reached confluence. All cells were cultured at 37 °C in 5% CO2 incubator.
Plasmids and transfection
The cDNA of myoD from mice, pigs, chicken and xenopus tropicalis were separately cloned into pcDNA3.1-flag vector (Invitrogen, Shanghai, China). NIH3T3 cells were transfected with myoD plasmid or control vector using Lipofectamine 3000 (Invitrogen) according to the manufacturer’s instruction. NIH3T3 cells were passaged to 6-well or 12-well plates 12 h before transfection. All transfections were performed in triplicate for each experiment.
RNA isolation and quantitative real-time PCR
Total RNA was isolated from NIH3T3 cells using TRIzol® Reagent (Invitrogen, Shanghai, China) and cDNA was synthesized from total RNA by Reverse Transcription Kit (Promega, Shanghai, China). Quantitative real-time PCR (qPCR) assays was performed on LightCycler 480 II system (Roche, Basel, Switzerland) by using a SYBR Green qPCR Kit (Genestar, Beijing, China). GAPDH is used as an internal control for normalization. The sequences of qPCR primers are shown in Additional file 2: Table S2. The experimental data were analyzed using 2-∆∆CT method.
NIH3T3 cells were treated with cell Lysis Buffer to completely release total protein. Cell extracts were separated by SDS-PAGE and the proteins were transferred to 0.45 μm PVDF membrane (Bio-Rad, Shanghai, China). Then the PVDF membranes were blocked with 3% BSA for 1 h at room temperature and then incubated at 4 °C overnight with primary antibodies. After washing, the membranes were labeled with proper secondary antibodies. Blots were visualized using a commercial enhanced chemiluminescene detection kit (Thermo Scientific, Beijing, China). GAPDH was used as the internal control. Primary antibodies used in this study included Flag (#8146S, CST) and GAPDH (#AP0063, Bioworld). Secondary antibodies included either anti-rabbit HRP-linked (#7074 S, CST) or anti-mouse HRP-linked (#7076 S, CST) antibodies.
Luciferase reporter assays
293 T Cells were seeded into a 24-well plate well and each plate of cells was transfected with 100 ng of 4R-TK-Luc and 100 ng of MyoD expression vectors or control vectors. Reporter activity was measured at 48 h later with a dual luciferase assay kit (Promega) and BioTek Synergy2. The luciferase activity was normalized by Renilla activity and the total transfection amount was normalized by an empty expression vector.
For myosin heavy-chain staining, cells in culture medium were rinsed with PBS, fixed with 4% paraformaldehyde, permeabilized with 0.5% Triton X-100 and blocked for 1 h with 3% BSA. Then incubated at 4 °C overnight with the anti-MyHC antibody (#ab51263, Abcam), washed three times with PBS and incubated for 1 h with secondary antibodies (#8940S, CST). After three washes by PBS, cells were stained with DAPI for 2 min. Images were captured by fluorescent reverse microscopy (ZEISS, Heidenheim, Germany).
Data are presented as mean ± SE from three independent experiments. Statistical significance was determined by the Student’s t-test, and P < 0.05 was considered as significance (*P < 0.05; **P < 0.01; ***P < 0.001).
Availability of data and materials
The datasets supporting the results of this article are included within the article (and its additional files). The myoD sequence data of Small-ear Spotted and Lantang pigs have been deposited in GenBank database (Accession number: PRJNA530874).
- bHLH domain:
Basic helix-loop-helix domain
Distal regulatory region
Myogenic Differentiation 1
Proximal regulatory region
Quantitative real-time PCR
Single nucleotide polymorphisms
Transcriptional start site
White JD, Scaffidi A, Davies M, McGeachie J, Rudnicki MA, Grounds MD. Myotube formation is delayed but not prevented in MyoD-deficient skeletal muscle: studies in regenerating whole muscle grafts of adult mice. J Histochem Cytochem. 2000;48(11):1531–44.
Andersson L, Georges M. Domestic-animal genomics: deciphering the genetics of complex traits. Nat Rev Genet. 2004;5(3):202–12.
Hettmer S, Wagers AJ. Muscling in: uncovering the origins of rhabdomyosarcoma. Nat Med. 2010;16(2):171–3.
Weintraub H, Tapscott SJ, Davis RL, Thayer MJ, Adam MA, Lassar AB, Miller AD. Activation of muscle-specific genes in pigment, nerve, fat, liver, and fibroblast cell lines by forced expression of MyoD. Proc Natl Acad Sci U S A. 1989;86(14):5434–8.
Singh K, Dilworth FJ. Differential modulation of cell cycle progression distinguishes members of the myogenic regulatory factor family of transcription factors. FEBS J. 2013;280(17):3991–4003.
Schuster-Gossler K, Cordes R, Gossler A. Premature myogenic differentiation and depletion of progenitor cells cause severe muscle hypotrophy in Delta1 mutants. Proc Natl Acad Sci U S A. 2007;104(2):537–42.
Zhao X, Mo D, Li A, Gong W, Xiao S, Zhang Y, Qin L, Niu Y, Guo Y, Liu X, et al. Comparative analyses by sequencing of transcriptomes during skeletal muscle development between pig breeds differing in muscle growth rate and fatness. PLoS One. 2011;6(5):e19774.
Chen JC, Goldhamer DJ. The core enhancer is essential for proper timing of MyoD activation in limb buds and branchial arches. Dev Biol. 2004;265(2):502–12.
Scionti I, Hayashi S, Mouradian S, Girard E, Esteves de Lima J, Morel V, Simonet T, Wurmser M, Maire P, Ancelin K, et al. LSD1 controls timely MyoD expression via MyoD Core enhancer transcription. Cell Rep. 2017;18(8):1996–2006.
Mueller AC, Cichewicz MA, Dey BK, Layer R, Reon BJ, Gagan JR, Dutta A. MUNC, a long noncoding RNA that facilitates the function of MyoD in skeletal myogenesis. Mol Cell Biol. 2015;35(3):498–513.
Garstang MG, Madapura PM. An enhancer-derived RNA muscles in to regulate Myogenin in trans. Mol Cell. 2018;71(1):3–5.
Tsai PF, Dell'Orso S, Rodriguez J, Vivanco KO, Ko KD, Jiang K, Juan AH, Sarshad AA, Vian L, Tran M, et al. A Muscle-Specific Enhancer RNA Mediates Cohesin Recruitment and Regulates Transcription In trans. Mol Cell. 2018;71(1):129–141.e128.
Chen JC, Love CM, Goldhamer DJ. Two upstream enhancers collaborate to regulate the spatial patterning and timing of MyoD transcription during mouse development. Dev Dyn. 2001;221(3):274–88.
Goldhamer DJ, Brunk BP, Faerman A, King A, Shani M, Emerson CP Jr. Embryonic activation of the myoD gene is regulated by a highly conserved distal control element. Development. 1995;121(3):637–49.
Tapscott SJ, Lassar AB, Weintraub H. A novel myoblast enhancer element mediates MyoD transcription. Mol Cell Biol. 1992;12(11):4994–5003.
Asakura A, Lyons GE, Tapscott SJ. The regulation of MyoD gene expression: conserved elements mediate expression in embryonic axial muscle. Dev Biol. 1995;171(2):386–98.
Huang Y, Chen B, Ye M, Liang P, Zhangfang Y, Huang J, Liu M, Songyang Z, Ma W. Ccndbp1 is a new positive regulator of skeletal myogenesis. J Cell Sci. 2016;129(14):2767–77.
Li JQ, Chen ZM, Liu DW, Liu XH, Sun BL, Ling F, Zhang H, Chen YS. Genetic effects of IGF-1 gene on the performance in landrace x Lantang pig resource population. Yi chuan xue bao =. Acta Genet Sin. 2003;30(9):835–9.
Newcom DW, Stalder KJ, Baas TJ, Goodwin RN, Parrish FC, Wiegand BR. Breed differences and genetic parameters of myoglobin concentration in porcine longissimus muscle. J Anim Sci. 2004;82(8):2264–8.
Tang Z, Li Y, Wan P, Li X, Zhao S, Liu B, Fan B, Zhu M, Yu M, Li K. LongSAGE analysis of skeletal muscle at three prenatal stages in Tongcheng and landrace pigs. Genome Biol. 2007;8(6):R115.
Suzuki A, Kojima N, Ikeuchi Y, Ikarashi S, Moriyama N, Ishizuka T, Tokushige H. Carcass composition and meat quality of Chinese purebred and European x Chinese crossbred pigs. Meat Sci. 1991;29(1):31–41.
Freilich M, Wen B, Shafer D, Schleier P, Dard M, Pendrys D, Ortiz D, Kuhn L. Implant-guided vertical bone growth in the mini-pig. Clin Oral Implants Res. 2012;23(6):751–7.
Rudnicki MA, Schnegelsberg PN, Stead RH, Braun T, Arnold HH, Jaenisch R. MyoD or Myf-5 is required for the formation of skeletal muscle. Cell. 1993;75(7):1351–9.
Law C, Cheung P. Expression of non-acetylatable H2A.Z in myoblast cells blocks myoblast differentiation through disruption of MyoD expression. J Biol Chem. 2015;290(21):13234–49.
Wang C, Liu W, Nie Y, Qaher M, Horton HE, Yue F, Asakura A, Kuang S. Loss of MyoD promotes fate Transdifferentiation of myoblasts into Brown adipocytes. EBioMedicine. 2017;16:212–23.
Tapscott SJ. The circuitry of a master switch: Myod and the regulation of skeletal muscle gene transcription. Development. 2005;132(12):2685–95.
Tapscott SJ, Davis RL, Thayer MJ, Cheng PF, Weintraub H, Lassar AB. MyoD1: a nuclear phosphoprotein requiring a Myc homology region to convert fibroblasts to myoblasts. Science (New York, NY). 1988;242(4877):405–11.
Fong AP, Yao Z, Zhong JW, Johnson NM, Farr GH 3rd, Maves L, Tapscott SJ. Conversion of MyoD to a neurogenic factor: binding site specificity determines lineage. Cell Rep. 2015;10(12):1937–46.
Di Gioia SA, Shaaban S, Tuysuz B, Elcioglu NH, Chan WM, Robson CD, Ecklund K, Gilette NM, Hamzaoglu A, Tayfun GA, et al. Recessive MYF5 mutations cause external Ophthalmoplegia, rib, and vertebral anomalies. Am J Hum Genet. 2018;103(1):115–24.
Harada A, Maehara K, Sato Y, Konno D, Tachibana T, Kimura H, Ohkawa Y. Incorporation of histone H3.1 suppresses the lineage potential of skeletal muscle. Nucleic Acids Res. 2015;43(2):775–86.
Luo D, de Morree A. Deltex2 represses MyoD expression and inhibits myogenic differentiation by acting as a negative regulator of Jmjd1c. Proc Natl Acad Sci. 2017;114(15):E3071–e3080.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics (Oxford, England). 2009;25(14):1754–60.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
We would like to thank Zuyong He for his some useful suggestions on this study.
This research was supported by the National Key R&D Program of China (2018YFD0501200), National Natural Science Foundation of China (31772565), Sailing Plan of Guangdong Province (2014YT02H042) and National Swine Industry Technology System (CARS–36).
All the animal procedures were carried out in accordance with China Council on Animal Care and the protocols we used were approved by the Animal Care and Use Committee of Guangdong Province, China. The approval ID or permit numbers are SCXK (Guangdong) 2011–0029 and SYXK (Guangdong) 2011–0112.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The information of all samples. (XLSX 295 kb)
Table S2. The primer sequences for q-PCR. (DOCX 13 kb)
Table S3. Summary of sequence similarity of myoD regulatory regions. (DOCX 15 kb)
Table S4. Sequence size and location of myoD regulatory regions. (DOCX 13 kb)
Figure S1. The similarity comparison of MyoD protein among different species. (DOCX 144 kb)