Transcriptomic profiles of muscular dystrophy with myositis (mdm) in extensor digitorum longus, psoas, and soleus muscles from mice.

BACKGROUND: Titinopathies are inherited muscular diseases triggered by genetic mutations in the titin gene. Muscular dystrophy with myositis (mdm) is one such disease caused by a LINE repeat insertion, leading to exon skipping and an 83-amino acid residue deletion in the N2A-PEVK region of mouse titin. This region has been implicated in a number of titin-titin ligand interactions, hence are important for myocyte signaling and health. Mice with this mdm mutation develop a severe and progressive muscle degeneration. The range of phenotypic differences observed in mdm mice shows that the deletion of this region induces a cascade of transcriptional changes extending to numerous signaling pathways affected by the titin filament. Previous research has focused on correlating phenotypic differences with muscle function in mdm mice. These studies have provided understanding of the downstream physiological effects resulting from the mdm mutation but only provide insights on processes that can be physiologically observed and measured. We used differential gene expression (DGE) to compare the transcriptomes of extensor digitorum longus (EDL), psoas and soleus muscles from wild-type and mdm mice to develop a deeper understand of these tissue-specific responses. RESULTS: The overall expression pattern observed shows a well-differentiated transcriptional signature in mdm muscles compared to wild type. Muscle-specific clusters observed within the mdm transcriptome highlight the level of variability of each muscle to the deletion. Differential gene expression and weighted gene co-expression network analysis showed a strong directional response in oxidative respiration-associated mitochondrial genes, which aligns with the poor shivering and non-shivering thermogenesis previously observed. Sln, which is a marker associated with shivering and non-shivering thermogenesis, showed the strongest expression change in fast-fibered muscles. No drastic changes in MYH expression levels were reported, which indicated an absence of major fiber-type switching events. Overall expression shifts in MYH isoforms, MARPs, and extracellular matrix associated genes demonstrated the transcriptional complexity associated with mdm mutation. The expression alterations in mitochondrial respiration and metabolism related genes in the mdm muscle dominated over other transcriptomic changes, and likely account for the late stage cellular responses in the mdm muscles. CONCLUSIONS: We were able to demonstrate that the complex nature of mdm mutation extends beyond a simple rearrangement in titin gene. EDL, psoas and soleus exemplify unique response modes observed in skeletal muscles with mdm mutation. Our data also raises the possibility that failure to maintain proper energy homeostasis in mdm muscles may contribute to the pathogenesis of the degenerative phenotype in mdm mice. Understanding the full disease-causing molecular cascade is difficult using bulk RNA sequencing techniques due to intricate nature of the disease. The development of the mdm phenotype is temporally and spatially regulated, hence future studies should focus on single fiber level investigations.

Bridging the muscle genome to phenome across multiple biological scales.

Muscle is highly hierarchically organized, with functions shaped by genetically controlled expression of protein ensembles with different isoform profiles at the sarcomere scale. However, it remains unclear how isoform profiles shape whole-muscle performance. We compared two mouse hindlimb muscles, the slow, relatively parallel-fibered soleus and the faster, more pennate-fibered tibialis anterior (TA), across scales: from gene regulation, isoform expression and translation speed, to force-length-velocity-power for intact muscles. Expression of myosin heavy-chain (MHC) isoforms directly corresponded with contraction velocity. The fast-twitch TA with fast MHC isoforms had faster unloaded velocities (actin sliding velocity, Vactin; peak fiber velocity, Vmax) than the slow-twitch soleus. For the soleus, Vactin was biased towards Vactin for purely slow MHC I, despite this muscle's even fast and slow MHC isoform composition. Our multi-scale results clearly identified a consistent and significant dampening in fiber shortening velocities for both muscles, underscoring an indirect correlation between Vactin and fiber Vmax that may be influenced by differences in fiber architecture, along with internal loading due to both passive and active effects. These influences correlate with the increased peak force and power in the slightly more pennate TA, leading to a broader length range of near-optimal force production. Conversely, a greater force-velocity curvature in the near-parallel fibered soleus highlights the fine-tuning by molecular-scale influences including myosin heavy and light chain expression along with whole-muscle characteristics. Our results demonstrate that the individual gene, protein and whole-fiber characteristics do not directly reflect overall muscle performance but that intricate fine-tuning across scales shapes specialized muscle function.

Contributions of Titin and Collagen to Passive Stress in Muscles from mdm Mice with a Small Deletion in Titin's Molecular Spring.

Muscular dystrophy with myositis (mdm) is a naturally occurring mutation in the mouse Ttn gene that results in higher passive stress in muscle fibers and intact muscles compared to wild-type (WT). The goal of this study was to test whether alternative splicing of titin exons occurs in mdm muscles, which contain a small deletion in the N2A-PEVK regions of titin, and to test whether splicing changes are associated with an increase in titin-based passive tension. Although higher levels of collagen have been reported previously in mdm muscles, here we demonstrate alternative splicing of titin in mdm skeletal muscle fibers. We identified Z-band, PEVK, and C-terminus Mex5 exons as splicing hotspots in mdm titin using RNA sequencing data and further reported upregulation in ECM-associated genes. We also treated skinned mdm soleus fiber bundles with trypsin, trypsin + KCl, and trypsin + KCL + KI to degrade titin. The results showed that passive stress dropped significantly more after trypsin treatment in mdm fibers (11 ± 1.6 mN/mm2) than in WT fibers (4.8 ± 1 mN/mm2; p = 0.0004). The finding that treatment with trypsin reduces titin-based passive tension more in mdm than in WT fibers supports the hypothesis that exon splicing leads to the expression of a stiffer and shorter titin isoform in mdm fibers. After titin extraction by trypsin + KCl + KI, mdm fibers (6.7 ± 1.27 mN/mm2) had significantly higher collagen-based passive stress remaining than WT fibers (2.6 ± 1.3 mN/mm2; p = 0.0014). We conclude that both titin and collagen contribute to higher passive tension of mdm muscles.

Comparative analysis of the transcriptomes of EDL, psoas, and soleus muscles from mice.

BACKGROUND: Individual skeletal muscles have evolved to perform specific tasks based on their molecular composition. In general, muscle fibers are characterized as either fast-twitch or slow-twitch based on their myosin heavy chain isoform profiles. This approach made sense in the early days of muscle studies when SDS-PAGE was the primary tool for mapping fiber type. However, Next Generation Sequencing tools permit analysis of the entire muscle transcriptome in a single sample, which allows for more precise characterization of differences among fiber types, including distinguishing between different isoforms of specific proteins. We demonstrate the power of this approach by comparing the differential gene expression patterns of extensor digitorum longus (EDL), psoas, and soleus from mice using high throughput RNA sequencing. RESULTS: EDL and psoas are typically classified as fast-twitch muscles based on their myosin expression pattern, while soleus is considered a slow-twitch muscle. The majority of the transcriptomic variability aligns with the fast-twitch and slow-twitch characterization. However, psoas and EDL exhibit unique expression patterns associated with the genes coding for extracellular matrix, myofibril, transcription, translation, striated muscle adaptation, mitochondrion distribution, and metabolism. Furthermore, significant expression differences between psoas and EDL were observed in genes coding for myosin light chain, troponin, tropomyosin isoforms, and several genes encoding the constituents of the Z-disk. CONCLUSIONS: The observations highlight the intricate molecular nature of skeletal muscles and demonstrate the importance of utilizing transcriptomic information as a tool for skeletal muscle characterization.

Microfluidics-free single-cell genomics with templated emulsification.

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.