MuscleAtlasExplorer: a web service for studying gene expression in human skeletal muscle

Multiomics profiling of DNA methylation, microRNA, and mRNA in skeletal muscle from monozygotic twin pairs discordant for type 2 diabetes identifies dysregulated genes controlling metabolism

BACKGROUND

A large proportion of skeletal muscle insulin resistance in type 2 diabetes (T2D) is caused by environmental factors.

METHODS

By applying multiomics mRNA, microRNA (miRNA), and DNA methylation platforms in biopsies from 20 monozygotic twin pairs discordant for T2D, we aimed to delineate the epigenetic and transcriptional machinery underlying non-genetic muscle insulin resistance in T2D.

RESULTS

Using gene set enrichment analysis (GSEA), we found decreased mRNA expression of genes involved in extracellular matrix organization, branched-chain amino acid catabolism, metabolism of vitamins and cofactors, lipid metabolism, muscle contraction, signaling by receptor tyrosine kinases pathways, and translocation of glucose transporter 4 (GLUT4) to the plasma membrane in muscle from twins with T2D. Differential expression levels of one or more predicted target relevant miRNA(s) were identified for approximately 35% of the dysregulated GSEA pathways. These include miRNAs with a significant overrepresentation of targets involved in GLUT4 translocation (miR-4643 and miR-548z), signaling by receptor tyrosine kinases pathways (miR-607), and muscle contraction (miR-4658). Acquired DNA methylation changes in skeletal muscle were quantitatively small in twins with T2D compared with the co-twins without T2D. Key methylation and expression results were validated in muscle, myotubes, and/or myoblasts from unrelated subjects with T2D and controls. Finally, mimicking T2D-associated changes by overexpressing miR-548 and miR-607 in cultured myotubes decreased expression of target genes, GLUT4 and FGFR4, respectively, and impaired insulin-stimulated phosphorylation of Akt (Ser473) and TBC1D4.

CONCLUSIONS

Together, we show that T2D is associated with non- and epigenetically determined differential transcriptional regulation of pathways regulating skeletal muscle metabolism and contraction.

MyoData: An expression knowledgebase at single cell/nucleus level for the discovery of coding-noncoding RNA functional interactions in skeletal muscle

•

Regulation of gene expression through non-coding RNAs at single myofiber and nucleus resolution.

•

Reinterpretation of KEGG pathways with microRNA and long non-coding RNA activities.

•

miR-149, -214, and let-7e alter mitochondrial shape.

•

The long non-coding RNA Pvt1 is a sponge for miR-27a.

•

miR-208b regulates Sox6; miR-214 regulates both Sox6 and Slc16a3.

Non-coding RNAs represent the largest part of transcribed mammalian genomes and prevalently exert regulatory functions. Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) can modulate the activity of each other. Skeletal muscle is the most abundant tissue in mammals. It is composed of different cell types with myofibers that represent the smallest complete contractile system. Considering that lncRNAs and miRNAs are more cell type-specific than coding RNAs, to understand their function it is imperative to evaluate their expression and action within single myofibers. In this database, we collected gene expression data for coding and non-coding genes in single myofibers and used them to produce interaction networks based on expression correlations. Since biological pathways are more informative than networks based on gene expression correlation, to understand how altered genes participate in the studied phenotype, we integrated KEGG pathways with miRNAs and lncRNAs. The database also integrates single nucleus gene expression data on skeletal muscle in different patho-physiological conditions. We demonstrated that these networks can serve as a framework from which to dissect new miRNA and lncRNA functions to experimentally validate. Some interactions included in the database have been previously experimentally validated using high throughput methods. These can be the basis for further functional studies. Using database information, we demonstrate the involvement of miR-149, -214 and let-7e in mitochondria shaping; the ability of the lncRNA Pvt1 to mitigate the action of miR-27a via sponging; and the regulatory activity of miR-214 on Sox6 and Slc16a3. The MyoData is available at https://myodata.bio.unipd.it.