Establishment of a human skeletal muscle-derived cell line: biochemical, cellular and electrophysiological characterization

Conditionally immortalised equine skeletal muscle cell lines for in vitro analysis

Background

Thoroughbred racehorse performance is largely influenced by a major quantitative trait locus at the myostatin (MSTN) gene which determines aptitude for certain race distances due to a promoter region insertion mutation influencing functional phenotypes in skeletal muscle. To develop an in vitro system for functional experiments we established three novel equine skeletal muscle cell lines reflecting the variation in phenotype associated with MSTN genotype (CC/II, CT/IN and TT/NN for SNP g.66493737C > T/SINE insertion 227 bp polymorphism). Primary equine skeletal muscle myoblasts, isolated from Thoroughbred horse gluteus medius, were conditionally immortalised and evaluated to determine whether cell phenotype and metabolic function were comparable to functional characteristics previously reported for ex vivo skeletal muscle isolated from Thoroughbred horses with each genotype.

Results

Primary myoblasts conditionally immortalised with the temperature sensitive SV40TtsA58 lentivirus vector successfully proliferated and could revert to their primary cell phenotype and differentiate into multinucleated myotubes. Skeletal muscle fibre type, MSTN gene expression, mitochondrial abundance, and mitochondrial function of the three MSTN genotype cell lines, were consistent with equivalent characterisation of ex vivo skeletal muscle samples with these genotypes. Furthermore, addition of coenzyme Q₁₀ (CoQ₁₀) to the cell lines improved mitochondrial function, an observation consistent with ex vivo skeletal muscle samples with these genotypes following supplementation with CoQ₁₀ in the diet.

Conclusions

The observation that the phenotypic characteristics and metabolic function of the cells lines are equivalent to ex vivo skeletal muscle indicates that this in vitro system will enable efficient and cost-effective analyses of equine skeletal muscle for a range of different applications including understanding metabolic function, testing of nutritional supplements, drug test development and gene doping test development. In the multi-billion-euro international Thoroughbred horse industry research advances in the biological function of skeletal muscle are likely to have considerable impact. Furthermore, this novel genotype-specific system may be adapted and applied to human biomedicine to improve understanding of the effects of myostatin in human physiology and medicine.

miR-19b-3p is associated with a diametric response to resistance exercise in older adults and regulates skeletal muscle anabolism via PTEN inhibition

Understanding paradoxical responses to anabolic stimulation and identifying the mechanisms for this inconsistency in mobility-limited older adults may provide new targets for the treatment of sarcopenia. Our laboratory has discovered that dysregulation in microRNA (miRNA) that target anabolic pathways is a potential mechanism resulting in age-associated decreases in skeletal muscle mass and function (sarcopenia). The objective of the current study was to assess circulating miRNA expression profiles in diametric response of leg lean mass in mobility-limited older individuals after a 6-mo progressive resistance exercise training intervention (PRET) and determine the influence of differentially expressing miRNA on regulation of skeletal muscle mass. Participants were dichotomized by gain (Gainers; mean +561.4 g, n = 33) or loss (Losers; mean −589.8 g, n = 40) of leg lean mass after PRET. Gainers significantly increased fat-free mass 2.4% vs. −0.4% for Losers. Six miRNA (miR-1-3p, miR-19b-3p, miR-92a, miR-126, miR-133a-3p, and miR-133b) were significantly identified to be differentially expressed between Gainers and Losers, with miR-19b-3p being the miRNA most highly associated with increases in fat-free mass. Using an aging mouse model, we then assessed if miR-19b-3p expression was different in young mice with larger muscle mass compared with older mice. Circulating and skeletal muscle miR-19b-3p expression was higher in young compared with old mice and was positively associated with muscle mass and grip strength. We then used a novel integrative approach to determine if differences in circulating miR-19b-3p potentially translate to augmented anabolic response in human skeletal muscle cells in vitro. Results from this analysis identified that overexpression of miR-19b-3p targeted and downregulated PTEN by 64% to facilitate significant ∼50% increase in muscle protein synthetic rate as measured with SUnSET. The combine results of these three models identify miR-19b-3p as a potent regulator of muscle anabolism that may contribute to an inter-individual response to PRET in mobility-limited older adults.

An integrated approach to identify environmental modulators of genetic risk factors for complex traits

Complex traits and diseases can be influenced by both genetics and environment. However, given the large number of environmental stimuli and power challenges for gene-by-environment testing, it remains a critical challenge to identify and prioritize specific disease-relevant environmental exposures. We propose a framework for leveraging signals from transcriptional responses to environmental perturbations to identify disease-relevant perturbations that can modulate genetic risk for complex traits and inform the functions of genetic variants associated with complex traits. We perturbed human skeletal-muscle-, fat-, and liver-relevant cell lines with 21 perturbations affecting insulin resistance, glucose homeostasis, and metabolic regulation in humans and identified thousands of environmentally responsive genes. By combining these data with GWASs from 31 distinct polygenic traits, we show that the heritability of multiple traits is enriched in regions surrounding genes responsive to specific perturbations and, further, that environmentally responsive genes are enriched for associations with specific diseases and phenotypes from the GWAS Catalog. Overall, we demonstrate the advantages of large-scale characterization of transcriptional changes in diversely stimulated and pathologically relevant cells to identify disease-relevant perturbations.

Nanopattern surface improves cultured human myotube maturation

BACKGROUND

In vitro maturation of human primary myoblasts using 2D culture remains a challenging process and leads to immature fibers with poor internal organization and function. This would however represent a valuable system to study muscle physiology or pathophysiology from patient myoblasts, at a single-cell level.

METHODS

Human primary myoblasts were cultured on 800-nm wide striated surface between two layers of Matrigel, and in a media supplemented with an inhibitor of TGFβ receptor. Gene expression, immunofluorescence, and Ca²⁺ measurements upon electrical stimulations were performed at various time points during maturation to assess the organization and function of the myotubes.

RESULTS

We show that after 10 days in culture, myotubes display numerous functional acetylcholine receptor clusters and express the adult isoforms of myosin heavy chain and dihydropyridine receptor. In addition, the myotubes are internally well organized with striations of α-actinin and STIM1, and occasionally ryanodine receptor 1. We also demonstrate that the myotubes present robust Ca²⁺ responses to repetitive electrical stimulations.

CONCLUSION

The present method describes a fast and efficient system to obtain well matured and functional myotubes in 2D culture allowing thorough analysis of single-cell Ca²⁺ signals.

Predictive network modeling in human induced pluripotent stem cells identifies key driver genes for insulin responsiveness

Insulin resistance (IR) precedes the development of type 2 diabetes (T2D) and increases cardiovascular disease risk. Although genome wide association studies (GWAS) have uncovered new loci associated with T2D, their contribution to explain the mechanisms leading to decreased insulin sensitivity has been very limited. Thus, new approaches are necessary to explore the genetic architecture of insulin resistance. To that end, we generated an iPSC library across the spectrum of insulin sensitivity in humans. RNA-seq based analysis of 310 induced pluripotent stem cell (iPSC) clones derived from 100 individuals allowed us to identify differentially expressed genes between insulin resistant and sensitive iPSC lines. Analysis of the co-expression architecture uncovered several insulin sensitivity-relevant gene sub-networks, and predictive network modeling identified a set of key driver genes that regulate these co-expression modules. Functional validation in human adipocytes and skeletal muscle cells (SKMCs) confirmed the relevance of the key driver candidate genes for insulin responsiveness.

Insulin resistance is characterized by a defective response (“resistance”) to normal insulin concentrations to uptake the glucose present in the blood, and is the underlying condition that leads to type 2 diabetes (T2D) and increases the risk of cardiovascular disease. It is estimated that 25–33% of the US population are insulin resistant enough to be at risk of serious clinical consequences. For more than a decade, large population studies have tried to discover the genes that participate in the development of insulin resistance, but without much success. It is now increasingly clear that the complex genetic nature of insulin resistance requires novel approaches centered in patient specific cellular models. To fill this gap, we have generated an induced pluripotent stem cell (iPSC) library from individuals with accurate measurements of insulin sensitivity, and performed gene expression and key driver analyses. Our work demonstrates that iPSCs can be used as a revolutionary technology to model insulin resistance and to discover key genetic drivers. Moreover, they can develop our basic knowledge of the disease, and are ultimately expected to increase the therapeutic targets to treat insulin resistance and type 2 diabetes.

Derivation and Characterization of Immortalized Human Muscle Satellite Cell Clones from Muscular Dystrophy Patients and Healthy Individuals

In Duchenne muscular dystrophy (DMD) patients, absence of dystrophin causes muscle wasting by impacting both the myofiber integrity and the properties of muscle stem cells (MuSCs). Investigation of DMD encompasses the use of MuSCs issued from human skeletal muscle. However, DMD-derived MuSC usage is restricted by the limited number of divisions that human MuSCs can undertake in vitro before losing their myogenic characteristics and by the scarcity of human material available from DMD muscle. To overcome these limitations, immortalization of MuSCs appears as a strategy. Here, we used CDK4/hTERT expression in primary MuSCs and we derived MuSC clones from a series of clinically and genetically characterized patients, including eight DMD patients with various mutations, four congenital muscular dystrophies and three age-matched control muscles. Immortalized cultures were sorted into single cells and expanded as clones into homogeneous populations. Myogenic characteristics and differentiation potential were tested for each clone. Finally, we screened various promoters to identify the preferred gene regulatory unit that should be used to ensure stable expression in the human MuSC clones. The 38 clonal immortalized myogenic cell clones provide a large collection of controls and DMD clones with various genetic defects and are available to the academic community.

Measurement of Skeletal Muscle Fiber Contractility with High-Speed Traction Microscopy

We describe a technique for simultaneous quantification of the contractile forces and cytosolic calcium dynamics of muscle fibers embedded in three-dimensional biopolymer gels under auxotonic loading conditions. We derive a scaling law for linear elastic matrices such as basement membrane extract hydrogels (Matrigel) that allows us to measure contractile force from the shape of the relaxed and contracted muscle cell and the Young's modulus of the matrix without further knowledge of the matrix deformations surrounding the cell and without performing computationally intensive inverse force reconstruction algorithms. We apply our method to isolated mouse flexor digitorum brevis (FDB) fibers that are embedded in 10 mg/mL Matrigel. Upon electrical stimulation, individual FDB fibers show twitch forces of 0.37 ± 0.15 μN and tetanic forces (100-Hz stimulation frequency) of 2.38 ± 0.71 μN, corresponding to a tension of 0.44 ± 0.25 kPa and 2.53 ± 1.17 kPa, respectively. Contractile forces of FDB fibers increase in response to caffeine and the troponin-calcium stabilizer tirasemtiv, similar to responses measured in whole muscle. From simultaneous high-speed measurements of cell length changes and cytosolic calcium concentration using confocal line scanning at a frequency of 2048 Hz, we show that twitch and tetanic force responses to electric pulses follow the low-pass filtered calcium signal. In summary, we present a technically simple high-speed method for measuring contractile forces and cytosolic calcium dynamics of single muscle fibers. We expect that our method will help to reduce preparation time, costs, and the number of sacrificed animals needed for experiments such as drug testing.

Cellular, biochemical and molecular changes in muscles from patients with X-linked myotubular myopathy due to MTM1 mutations

Centronuclear myopathies are early-onset muscle diseases caused by mutations in several genes including MTM1, DNM2, BIN1, RYR1 and TTN. The most severe and often fatal X-linked form of myotubular myopathy (XLMTM) is caused by mutations in the gene encoding the ubiquitous lipid phosphatase myotubularin, an enzyme specifically dephosphorylating phosphatidylinositol-3-phosphate and phosphatidylinositol-3,5-bisphosphate. Because XLMTM patients have a predominantly muscle-specific phenotype a number of pathogenic mechanisms have been proposed, including a direct effect of the accumulated lipid on the skeletal muscle calcium channel ryanodine receptor 1, a negative effect on the structure of intracellular organelles and defective autophagy. Animal models knocked out for MTM1 show severe reduction of ryanodine receptor 1 mediated calcium release but, since knocking out genes in animal models does not necessarily replicate the human phenotype, we considered it important to study directly the effect of MTM1 mutations on patient muscle cells. The results of the present study show that at the level of myotubes MTM1 mutations do not dramatically affect calcium homeostasis and calcium release mediated through the ryanodine receptor 1, though they do affect myotube size and nuclear content. On the other hand, mature muscles such as those obtained from patient muscle biopsies exhibit a significant decrease in expression of the ryanodine receptor 1, a decrease in muscle-specific microRNAs and a considerable up-regulation of histone deacetylase-4. We hypothesize that the latter events consequent to the primary genetic mutation, are the cause of the severe decrease in muscle strength that characterizes these patients.

Skeletal muscle characteristics are preserved in hTERT/cdk4 human myogenic cell lines

Background

hTERT/cdk4 immortalized myogenic human cell lines represent an important tool for skeletal muscle research, being used as therapeutically pertinent models of various neuromuscular disorders and in numerous fundamental studies of muscle cell function. However, the cell cycle is linked to other cellular processes such as integrin regulation, the PI3K/Akt pathway, and microtubule stability, raising the question as to whether genetic modification related to the cell cycle results in secondary effects that could undermine the validity of these cell models.

Results

Here we subjected five healthy and disease muscle cell isolates to transcriptomic analysis, comparing immortalized lines with their parent primary populations in both differentiated and undifferentiated states, and testing their myogenic character by comparison with non-myogenic (CD56-negative) cells. Principal component analysis of global gene expression showed tight clustering of immortalized myoblasts to their parent primary populations, with clean separation from the non-myogenic reference. Comparison was made to publicly available transcriptomic data from studies of muscle human pathology, cell, and animal models, including to derive a consensus set of genes previously shown to have altered regulation during myoblast differentiation. Hierarchical clustering of samples based on gene expression of this consensus set showed that immortalized lines retained the myogenic expression patterns of their parent primary populations. Of 2784 canonical pathways and gene ontology terms tested by gene set enrichment analysis, none were significantly enriched in immortalized compared to primary cell populations. We observed, at the whole transcriptome level, a strong signature of cell cycle shutdown associated with senescence in one primary myoblast population, whereas its immortalized clone was protected.

Conclusions

Immortalization had no observed effect on the myogenic cascade or on any other cellular processes, and it was protective against the systems level effects of senescence that are observed at higher division counts of primary cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13395-016-0115-5) contains supplementary material, which is available to authorized users.

Functional characterization of orbicularis oculi and extraocular muscles

Facial muscles are skeletal muscles that control facial expression. Sekulic-Jablanovic et al. characterize orbicularis oculi and extraocular muscles and find divergence in the expression of key molecules for muscle function between facial, extraocular, and quadriceps muscles.

The orbicularis oculi are the sphincter muscles of the eyelids and are involved in modulating facial expression. They differ from both limb and extraocular muscles (EOMs) in their histology and biochemistry. Weakness of the orbicularis oculi muscles is a feature of neuromuscular disorders affecting the neuromuscular junction, and weakness of facial muscles and ptosis have also been described in patients with mutations in the ryanodine receptor gene. Here, we investigate human orbicularis oculi muscles and find that they are functionally more similar to quadriceps than to EOMs in terms of excitation–contraction coupling components. In particular, they do not express the cardiac isoform of the dihydropyridine receptor, which we find to be highly expressed in EOMs where it is likely responsible for the large depolarization-induced calcium influx. We further show that human orbicularis oculi and EOMs express high levels of utrophin and low levels of dystrophin, whereas quadriceps express dystrophin and low levels of utrophin. The results of this study highlight the notion that myotubes obtained by explanting satellite cells from different muscles are not functionally identical and retain the physiological characteristics of their muscle of origin. Furthermore, our results indicate that sparing of facial and EOMs in patients with Duchenne muscular dystrophy is the result of the higher levels of utrophin expression.

β-Arrestin scaffolds and signaling elements essential for the obestatin/GPR39 system that determine the myogenic program in human myoblast cells

Obestatin/GPR39 signaling stimulates skeletal muscle repair by inducing the expansion of satellite stem cells as well as myofiber hypertrophy. Here, we describe that the obestatin/GPR39 system acts as autocrine/paracrine factor on human myogenesis. Obestatin regulated multiple steps of myogenesis: myoblast proliferation, cell cycle exit, differentiation and recruitment to fuse and form multinucleated hypertrophic myotubes. Obestatin-induced mitogenic action was mediated by ERK1/2 and JunD activity, being orchestrated by a G-dependent mechanism. At a later stage of myogenesis, scaffolding proteins β-arrestin 1 and 2 were essential for the activation of cell cycle exit and differentiation through the transactivation of the epidermal growth factor receptor (EGFR). Upon obestatin stimulus, β-arrestins are recruited to the membrane, where they functionally interact with GPR39 leading to Src activation and signalplex formation to EGFR transactivation by matrix metalloproteinases. This signalplex regulated the mitotic arrest by p21 and p57 expression and the mid- to late stages of differentiation through JNK/c-Jun, CAMKII, Akt and p38 pathways. This finding not only provides the first functional activity for β-arrestins in myogenesis but also identify potential targets for therapeutic approaches by triggering specific signaling arms of the GPR39 signaling involved in myogenesis.

The metabolomic profile of gamma-irradiated human hepatoma and muscle cells reveals metabolic changes consistent with the Warburg effect

The two human cell lines HepG2 from hepatoma and HMCL-7304 from striated muscle were γ-irradiated with doses between 0 and 4 Gy. Abundant γH2AX foci were observed at 4 Gy after 4 h of culture post-irradiation. Sham-irradiated cells showed no γH2AX foci and therefore no signs of radiation-induced double-strand DNA breaks. Flow cytometry indicated that 41.5% of HepG2 cells were in G2/M and this rose statistically significantly with increasing radiation dose reaching a plateau at ∼47%. Cell lysates from both cell lines were subjected to metabolomic analysis using Gas Chromatography-Mass Spectrometry (GCMS). A total of 46 metabolites could be identified by GCMS in HepG2 cell lysates and 29 in HMCL-7304 lysates, most of which occurred in HepG2 cells. Principal Components Analysis (PCA) showed a clear separation of sham, 1, 2 and 4 Gy doses. Orthogonal Projection to Latent Structures-Discriminant Analysis (OPLS-DA) revealed elevations in intracellular lactate, alanine, glucose, glucose 6-phosphate, fructose and 5-oxoproline, which were found by univariate statistics to be highly statistically significantly elevated at both 2 and 4 Gy compared with sham irradiated cells. These findings suggested upregulation of cytosolic aerobic glycolysis (the Warburg effect), with potential shunting of glucose through aldose reductase in the polyol pathway, and consumption of reduced Glutathione (GSH) due to γ-irradiation. In HMCL-7304 myotubes, a putative Warburg effect was also observed only at 2 Gy, albeit a lesser magnitude than in HepG2 cells. It is anticipated that these novel metabolic perturbations following γ-irradiation of cultured cells will lead to a fuller understanding of the mechanisms of tissue damage following ionizing radiation exposure.

Intracellular Ca(2+)-handling differs markedly between intact human muscle fibers and myotubes

Background

In skeletal muscle, intracellular Ca²⁺ is an important regulator of contraction as well as gene expression and metabolic processes. Because of the difficulties to obtain intact human muscle fibers, human myotubes have been extensively employed for studies of Ca²⁺-dependent processes in human adult muscle. Despite this, it is unknown whether the Ca²⁺-handling properties of myotubes adequately represent those of adult muscle fibers.

Methods

To enable a comparison of the Ca²⁺-handling properties of human muscle fibers and myotubes, we developed a model of dissected intact single muscle fibers obtained from human intercostal muscle biopsies. The intracellular Ca²⁺-handling of human muscle fibers was compared with that of myotubes generated by the differentiation of primary human myoblasts obtained from vastus lateralis muscle biopsies.

Results

The intact single muscle fibers all demonstrated strictly regulated cytosolic free [Ca²⁺] ([Ca²⁺]_i) transients and force production upon electrical stimulation. In contrast, despite a more mature Ca²⁺-handling in myotubes than in myoblasts, myotubes lacked fundamental aspects of adult Ca²⁺-handling and did not contract. These functional differences were explained by discrepancies in the quantity and localization of Ca²⁺-handling proteins, as well as ultrastructural differences between muscle fibers and myotubes.

Conclusions

Intact single muscle fibers that display strictly regulated [Ca²⁺]_i transients and force production upon electrical stimulation can be obtained from human intercostal muscle biopsies. In contrast, human myotubes lack important aspects of adult Ca²⁺-handling and are thus an inappropriate model for human adult muscle when studying Ca²⁺-dependent processes, such as gene expression and metabolic processes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13395-015-0050-x) contains supplementary material, which is available to authorized users.

Characterization of excitation–contraction coupling components in human extraocular muscles

We show that the expression level of RyR1 in human extraocular muscles (EOMs) is low and that these muscles express different levels of proteins involved in excitation–contraction coupling (ECC) compared with leg muscles (LMs).