Contributions of alternative splicing to muscle type development and function

Rbm24-mediated post-transcriptional regulation of skeletal and cardiac muscle development, function and regeneration

RNA-binding proteins are critically involved in the post-transcriptional control of gene expression during embryonic development and in adult life, contributing to regulating cell differentiation and maintaining tissue homeostasis. Compared to the relatively well documented functions of transcription factors, the regulatory roles of RNA-binding proteins in muscle development and function remain largely elusive. However, deficiency of many RNA-binding proteins has been associated with muscular defects, neuromuscular disorders and heart diseases, such as myotonic dystrophy, amyotrophic lateral sclerosis, and cardiomyopathy. Rbm24 is highly conserved among vertebrates and is one of the best characterized RNA-binding proteins with crucial implication in the myogenic and cardiomyogenic programs. It presents the distinctive particularity of displaying highly restricted expression in both skeletal and cardiac muscles, with changes in subcellular localization during the process of differentiation. Functional analyses using different vertebrate models have clearly demonstrated its requirement for skeletal muscle differentiation and regeneration as well as for myocardium organization and cardiac function, by regulating the expression of both common and distinct target genes in these tissues. The challenge remains to decipher the dynamic feature of post-transcriptional circuits regulated by Rbm24 during skeletal myogenesis, cardiomyogenesis, and muscle repair. This review discusses current understanding of its function in striated muscles and its possible implication in human disease, with the aim of identifying research gaps for future investigation.

The m6A modification regulates the composition of myofiber types in chicken skeletal muscle

As a widespread epigenetic RNA modification, N6-methyladenosine (m6A) plays essential regulatory roles in multiple biological processes. However, its function in maintaining and modulating myofiber-type properties remains largely unknown. To investigate the post-transcriptional modification underlying the myofiber type diversity in chicken skeletal muscle, we evaluated the m6A methylation levels of chicken skeletal muscles with different phenotypic traits, and profiled a transcriptome-wide m6A map in the oxidative and glycolytic skeletal muscles by methylated RNA immunoprecipitation sequencing (MeRIP-seq). Our results showed that the levels of m6A methylation in chicken skeletal muscles were closely related to the composition of myofiber types. The m6A methylation level of anterior latissimus dorsi (ALD, typical oxidative skeletal muscle) was the highest among the three muscles and significantly higher than that of the pectoralis major (PM, typical glycolytic skeletal muscle) (P<0.05). We found that about 24.77 % and 33.50 % of genes were modified by m6A methylation in the PM and ALD, respectively, and identified 6,530 and 9,965 m6A peaks, which were mainly located in the coding sequence (CDS) and stop codon. About 3.14 % of m6A modified genes showed significantly differential methylation levels between these two muscles. Intriguingly, the myofiber type-related genes, such as MYOT, TPM3, TPM1, PDK1, MBNL1, and MYH1G, showed differences in m6A methylation and mRNA expression. Further analysis revealed that the m6A methylation was positively correlated with gene expression homeostasis. It is exciting we found that the expression level of ALKBH5 mRNA and protein, was closely related to the composition of myofiber types. ALKBH5 over-expression could regulate the expression levels of genes related to muscle contraction and metabolism, including MYH1E, MYH1G, MYH7B, PDK1, and TPM1, suggesting the effect of ALKBH5 on the formation of myofiber-type properties in chicken skeletal muscle. Our results contribute to a better understanding of epigenetic factors involved in forming chicken myofiber-type properties and provide new targets for further investigation into chicken's growth development and meat quality.

Differential inclusion of NEB exons 143 and 144 provides insight into NEB-related myopathy variant interpretation and disease manifestation

Biallelic pathogenic variants in the gene encoding nebulin (NEB) are a known cause of congenital myopathy. We present two brothers with congenital myopathy and compound heterozygous variants (NC_000002.12:g.151692086G>T; NM_001271208.2: c.2079C>A; p.(Cys693Ter) and NC_000002.12:g.151533439T>C; NM_001271208.2:c.21522+3A>G) in NEB. Transcriptomic sequencing on affected individual muscles revealed that the extended splice variant c.21522+3A>G causes exon 144 skipping. Nebulin isoforms containing exon 144 are known to be mutually exclusive with isoforms containing exon 143, and these isoforms are differentially expressed during development and in adult skeletal muscles. Affected individuals' MRI patterns of muscle involvement were compared with the known pattern of relative abundance of these two isoforms in muscle. We propose that the pattern of muscle involvement in these affected individuals better fits the distribution of exon 144-containing isoforms in muscle than with previously published MRI findings in NEB-related disease due to other variants. Our report introduces disease pathogenesis and manifestation as a result of alteration of isoform distributions in muscle.

Alternative splicing perspective to prey preference of environmentally friendly biological agent Cryptolaemus montrouzieri

BACKGROUND

Cryptolaemus montrouzieri (Coccinellidae) is widely utilized as biological control agents in modern agriculture. A comprehensive understanding of its food preference can help guide mass rearing and safety management during field application of pest control. Although some studies have paid attentions to the impacts of prey shift on C. montrouzieri, little is known regarding the role of post-transcriptional regulations in its acclimation to unnatural preys.

RESULTS

We performed a genome-wide investigation on alternative splicing dynamics in C. montrouzieri in response to the predation transition from natural prey to unnatural ones. When feeding on undesired diets, 402-764 genes were differentially alternative spliced in C. montrouzieri. It is noteworthy that the majority of these genes (> 87%) were not differentially expressed, and these differentially spliced genes regulated distinct biological processes from differentially expressed genes, such as organ development and morphogenesis, locomotory behavior, and homeostasis processes. These suggested the functionally nonredendant role of alternative splicing in modulating physiological and metabolic responses of C. montrouzieri to the shift to undesired preys. In addition, the individuals feeding on aphids were subject to a lower level of changes in splicing than other alternative diets, which might be because of the similar chemical and microbial compositions. Our study further suggested a putative coupling of alternative splicing and nonsense-mediated decay (AS-NMD), which may play an important role in fine-tuning the protein repertoire of C. montrouzieri, and promoting its acclimation to predation changes.

CONCLUSION

These findings highlight the key role of alternative splicing in modulating the acclimation of ladybirds to prey shift and provide new genetic clues for the future application of ladybirds in biocontrol.

Mining alternative splicing patterns in scRNA-seq data using scASfind

Single-cell RNA-seq (scRNA-seq) is widely used for transcriptome profiling, but most analyses focus on gene-level events, with less attention devoted to alternative splicing. Here, we present scASfind, a novel computational method to allow for quantitative analysis of cell type-specific splicing events using full-length scRNA-seq data. ScASfind utilizes an efficient data structure to store the percent spliced-in value for each splicing event. This makes it possible to exhaustively search for patterns among all differential splicing events, allowing us to identify marker events, mutually exclusive events, and events involving large blocks of exons that are specific to one or more cell types.

Undetected Neuromuscular Disease in Patients after Heart Transplantation

(1) Heart transplantation (HTX) improves the overall survival and functional status of end-stage heart failure patients with cardiomyopathies (CMPs). The majority of CMPs have genetic causes, and the overlap between CMPs and inherited myopathies is well documented. However, the long-term outcome in skeletal muscle function and possibility of an undiagnosed underlying genetic cause of both a cardiac and skeletal pathology remain unknown. (2) Thirty-nine patients were assessed using open and standardized interviews on muscle function, a quality-of-life (EuroQol EQ-5D-3L) questionnaire, and a physical examination (Medical Research Council Muscle scale). Whole-exome sequencing was completed in three stages for those with skeletal muscle weakness. (3) Seven patients (17.9%) reported new-onset muscle weakness and motor limitations. Objective muscle weakness in the upper and lower extremities was seen in four patients. In three of them, exome sequencing revealed pathogenic/likely pathogenic variants in the genes encoding nexilin, myosin heavy chain, titin, and SPG7. (4) Our findings support a positive long-term outcome of skeletal muscle function in HTX patients. However, 10% of patients showed clinical signs of myopathy due to a possible genetic cause. The integration of genetic testing and standardized neurological assessment of motor function during the peri-HTX period should be considered.

Alternative splicing plays a nonredundant role in greater amberjack (Seriola dumerili) in acclimation to ambient salinity fluctuations

Alternative splicing (AS) is an important post-transcriptional mechanism for adaptation of fish to environmental stress. Here, we performed a genome-wide investigation to AS dynamics in greater amberjack (Seriola dumerili), an economical marine teleost, in response to hypo- (10 ppt) and hyper-salinity (40 ppt) stresses. Totally, 2267-2611 differentially spliced events were identified in gills and kidney upon the exposure to undesired salinity regimes. In gills, genes involved in energy metabolism, stimulus response and epithelial cell differentiation were differentially spliced in response to salinity variation, while sodium ion transport and cellular amide metabolism were enhanced in kidney to combat the adverse impacts of salinity changes. Most of these differentially spliced genes were not differentially expressed, and AS was found to regulate different biological processes from differential gene expression, indicative of the functionally nonredundant role of AS in modulating salinity acclimation in greater amberjack. Together, our study highlights the important contribution of post-transcriptional mechanisms to the adaptation of fish to ambient salinity fluctuations and provides theoretical guidance for the conservation of marine fishery resources against increasingly environmental challenges.

Divergent Regulatory Roles of Transcriptional Variants of the Chicken LDB3 Gene in Muscle Shaping

LIM domain binding 3 (LDB3) serves as a striated muscle-specific Z-band alternatively spliced protein that plays an important role in mammalian skeletal muscle development, but its regulatory role and molecular mechanism in avian muscle development are still unclear. In this study, we reanalyzed RNA sequencing data sets of 1415 samples from 21 chicken tissues published in the NCBI GEO database. First, three variants (LDB3-X, LDB3-XN1, and LDB3-XN2) generated by alternative splicing of the LDB3 gene were identified in chicken skeletal muscle, among which LDB3-XN1 and LDB3-XN2 are novel variants. LDB3-X and LDB3-XN1 are derived from exon skipping in chicken skeletal muscle at the E18-D7 stage and share three LIM domains, but LDB3-XN2 lacks a LIM domain. Our results preliminarily suggest that the formation of three variants of LDB3 is regulated by RBM20. The three splice isomers have divergent functions in skeletal muscle according to in vitro and in vivo assays. Finally, we identified the mechanism by which different variants play different roles through interactions with IGF2BP1 and MYHC, which promote the proliferation and differentiation of chicken myoblasts, in turn regulating chicken myogenesis. In conclusion, this study revealed the divergent roles of three LDB3 variants in chicken myogenesis and muscle remodeling and demonstrated their regulatory mechanism through protein-protein interactions.

Transcriptome analysis of Apis mellifera antennae reveals molecular divergence underlying the division of labour in worker bees

The olfactory system plays a fundamental role in mediating insect behaviour. Worker bees exhibit an age-dependent division of labour, performing discrete sets of behaviours throughout their lifespan. The behavioural states of bees rely on their sense of the environment and chemical communication via their olfactory system, the antennae. However, the olfactory adaptation mechanism of worker bees during their behavioural development remains unclear. In this study, we conducted a comprehensive and quantitative analysis of antennal gene expression in the Apis mellifera of newly emerged workers, nurses, foragers and defenders using RNA-seq. We found that the antenna tissues of honey bees continued developing after transformation from newly emerged workers to adults. Additionally, we identified differentially expressed genes associated with bee development and division of labour. We validated that major royal jelly protein genes are highly and specifically expressed in nurse honey bee workers. Furthermore, we identified and validated significant alternative splicing events correlated with the development and division of labour. These findings provide a comprehensive transcriptome profile and a new perspective on the molecular mechanisms that may underlie the worker honey bee division of labour.

Differential inclusion of NEB exons 143 and 144 provides insight into NEB-related myopathy variant interpretation and disease manifestation

Biallelic pathogenic variants in the gene encoding nebulin (NEB) are a known cause of congenital myopathy. We present two individuals with congenital myopathy and compound heterozygous variants (NM_001271208.2: c.2079C>A; p.(Cys693Ter) and c.21522+3A>G ) in NEB. Transcriptomic sequencing on patient muscle revealed that the extended splice variant c.21522+3A>G causes exon 144 skipping. Nebulin isoforms containing exon 144 are known to be mutually exclusive with isoforms containing exon 143, and these isoforms are differentially expressed during development and in adult skeletal muscles. Patients MRIs were compared to the known pattern of relative abundance of these two isoforms in muscle. We propose that the pattern of muscle involvement in these patients better fits the distribution of exon 144-containing isoforms in muscle than with previously published MRI findings in NEB-related disease due to other variants. To our knowledge this is the first report hypothesizing disease pathogenesis through the alteration of isoform distributions in muscle.

PINK1 regulated mitophagy is evident in skeletal muscles

PINK1, mutated in familial forms of Parkinson's disease, initiates mitophagy following mitochondrial depolarization. However, it is difficult to monitor this pathway physiologically in mice as loss of PINK1 does not alter basal mitophagy levels in most tissues. To further characterize this pathway in vivo, we used mito-QC mice in which loss of PINK1 was combined with the mitochondrial-associated POLGD257A mutation. We focused on skeletal muscle as gene expression data indicates that this tissue has the highest PINK1 levels. We found that loss of PINK1 in oxidative hindlimb muscle significantly reduced mitophagy. Of interest, the presence of the POLGD257A mutation, while having a minor effect in most tissues, restored levels of muscle mitophagy caused by the loss of PINK1. Although our observations highlight that multiple mitophagy pathways operate within a single tissue, we identify skeletal muscle as a tissue of choice for the study of PINK1-dependant mitophagy under basal conditions.

Peripheral thickening of the sarcomeres and pointed end elongation of the thin filaments are both promoted by SALS and its formin interaction partners

During striated muscle development the first periodically repeated units appear in the premyofibrils, consisting of immature sarcomeres that must undergo a substantial growth both in length and width, to reach their final size. Here we report that, beyond its well established role in sarcomere elongation, the Sarcomere length short (SALS) protein is involved in Z-disc formation and peripheral growth of the sarcomeres. Our protein localization data and loss-of-function studies in the Drosophila indirect flight muscle strongly suggest that radial growth of the sarcomeres is initiated at the Z-disc. As to thin filament elongation, we used a powerful nanoscopy approach to reveal that SALS is subject to a major conformational change during sarcomere development, which might be critical to stop pointed end elongation in the adult muscles. In addition, we demonstrate that the roles of SALS in sarcomere elongation and radial growth are both dependent on formin type of actin assembly factors. Unexpectedly, when SALS is present in excess amounts, it promotes the formation of actin aggregates highly resembling the ones described in nemaline myopathy patients. Collectively, these findings helped to shed light on the complex mechanisms of SALS during the coordinated elongation and thickening of the sarcomeres, and resulted in the discovery of a potential nemaline myopathy model, suitable for the identification of genetic and small molecule inhibitors.

Analysis of alternative splicing in chicken macrophages transfected with overexpression/knockdown of RIP2 gene

Receptor-interacting protein 2 (RIP2) plays a critical role in the transduction of many signaling pathways and is associated with many diseases. Alternative splicing (AS) is an essential and ubiquitous regulatory mechanism of gene expression that contributes to distinct transcript variants and many different kinds of proteins. In this present study, we characterized genome-wide AS events in wild-type chicken macrophages (WT) and RIP2 overexpression/knockdown chicken macrophages (oeRIP2/shRIP2) by high-throughput RNA sequencing technology. A total of 1901, 2061, and 817 differentially expressed (DE) AS genes were identified in the comparison of oeRIP2 vs. WT, oeRIP2 vs. shRIP2, and shRIP2 vs. WT, respectively. These DE AS genes participated in many important KEGG pathways, including regulation of autophagy, Wnt signaling pathway, Ubiquitin mediated proteolysis, MAPK signaling pathway, and Focal adhesion, etc. In conclusion, this research provided a broad atlas of the genome-wide scale of the AS event landscape in RIP2 overexpression/knockdown and wild-type chicken macrophages. This research also provides the theoretical basis of the gene network related to RIP2.

Alternative Transcripts Diversify Genome Function for Phenome Relevance to Health and Diseases

Manipulation using alternative exon splicing (AES), alternative transcription start (ATS), and alternative polyadenylation (APA) sites are key to transcript diversity underlying health and disease. All three are pervasive in organisms, present in at least 50% of human protein-coding genes. In fact, ATS and APA site use has the highest impact on protein identity, with their ability to alter which first and last exons are utilized as well as impacting stability and translation efficiency. These RNA variants have been shown to be highly specific, both in tissue type and stage, with demonstrated importance to cell proliferation, differentiation and the transition from fetal to adult cells. While alternative exon splicing has a limited effect on protein identity, its ubiquity highlights the importance of these minor alterations, which can alter other features such as localization. The three processes are also highly interwoven, with overlapping, complementary, and competing factors, RNA polymerase II and its CTD (C-terminal domain) chief among them. Their role in development means dysregulation leads to a wide variety of disorders and cancers, with some forms of disease disproportionately affected by specific mechanisms (AES, ATS, or APA). Challenges associated with the genome-wide profiling of RNA variants and their potential solutions are also discussed in this review.

A roadmap for delivering a human musculoskeletal cell atlas

Advances in single-cell technologies have transformed the ability to identify the individual cell types present within tissues and organs. The musculoskeletal bionetwork, part of the wider Human Cell Atlas project, aims to create a detailed map of the healthy musculoskeletal system at a single-cell resolution throughout tissue development and across the human lifespan, with complementary generation of data from diseased tissues. Given the prevalence of musculoskeletal disorders, this detailed reference dataset will be critical to understanding normal musculoskeletal function in growth, homeostasis and ageing. The endeavour will also help to identify the cellular basis for disease and lay the foundations for novel therapeutic approaches to treating diseases of the joints, soft tissues and bone. Here, we present a Roadmap delineating the critical steps required to construct the first draft of a human musculoskeletal cell atlas. We describe the key challenges involved in mapping the extracellular matrix-rich, but cell-poor, tissues of the musculoskeletal system, outline early milestones that have been achieved and describe the vision and directions for a comprehensive musculoskeletal cell atlas. By embracing cutting-edge technologies, integrating diverse datasets and fostering international collaborations, this endeavour has the potential to drive transformative changes in musculoskeletal medicine.

Comparative analyses of dynamic transcriptome profiles highlight key response genes and dominant isoforms for muscle development and growth in chicken

BACKGROUND

Modern breeding strategies have resulted in significant differences in muscle mass between indigenous chicken and specialized broiler. However, the molecular regulatory mechanisms that underlie these differences remain elusive. The aim of this study was to identify key genes and regulatory mechanisms underlying differences in breast muscle development between indigenous chicken and specialized broiler.

RESULTS

Two time-series RNA-sequencing profiles of breast muscles were generated from commercial Arbor Acres (AA) broiler (fast-growing) and Chinese indigenous Lushi blue-shelled-egg (LS) chicken (slow-growing) at embryonic days 10, 14, and 18, and post-hatching day 1 and weeks 1, 3, and 5. Principal component analysis of the transcriptome profiles showed that the top four principal components accounted for more than 80% of the total variance in each breed. The developmental axes between the AA and LS chicken overlapped at the embryonic stages but gradually separated at the adult stages. Integrative investigation of differentially-expressed transcripts contained in the top four principal components identified 44 genes that formed a molecular network associated with differences in breast muscle mass between the two breeds. In addition, alternative splicing analysis revealed that genes with multiple isoforms always had one dominant transcript that exhibited a significantly higher expression level than the others. Among the 44 genes, the TNFRSF6B gene, a mediator of signal transduction pathways and cell proliferation, harbored two alternative splicing isoforms, TNFRSF6B-X1 and TNFRSF6B-X2. TNFRSF6B-X1 was the dominant isoform in both breeds before the age of one week. A switching event of the dominant isoform occurred at one week of age, resulting in TNFRSF6B-X2 being the dominant isoform in AA broiler, whereas TNFRSF6B-X1 remained the dominant isoform in LS chicken. Gain-of-function assays demonstrated that both isoforms promoted the proliferation of chicken primary myoblasts, but only TNFRSF6B-X2 augmented the differentiation and intracellular protein content of chicken primary myoblasts.

CONCLUSIONS

For the first time, we identified several key genes and dominant isoforms that may be responsible for differences in muscle mass between slow-growing indigenous chicken and fast-growing commercial broiler. These findings provide new insights into the regulatory mechanisms underlying breast muscle development in chicken.

RUNX1T1 modulates myogenic differentiation by regulating the calcium signaling pathway and the alternative splicing of ROCK2

RUNX1T1 (Runt-related transcription factor 1, translocated to 1) plays a wide-ranging and diverse role in cellular development, including hematopoiesis and adipogenesis. However, little is known about the function of RUNX1T1 in the skeletal muscle development. Here, we assessed the impact of RUNX1T1 on the proliferation and myogenic differentiation of goat primary myoblasts (GPMs). It was observed that RUNX1T1 is highly expressed during the early stages of myogenic differentiation and the fetal stage. Moreover, the knockdown of RUNX1T1 promotes the proliferation and inhibits myogenic differentiation and mitochondrial biogenesis of GPMs. RNA sequencing analysis revealed that significantly differentially expressed genes in RUNX1T1 knockdown cells were enriched in the calcium signaling pathway. Additionally, we discovered that RUNX1T1 regulates alternative splicing (AS) events involved in myogenesis. We also show that silencing RUNX1T1 blocked the Ca²⁺ -CAMK signaling pathway and reduced the expression levels of muscle-specific isoforms of recombinant rho associated coiled coil containing crotein kinase 2 (ROCK2) during myogenic differentiation, partially explaining why RUNX1T1 deficiency leads to the impairment of myotube formation. These findings suggest that RUNX1T1 is a novel regulator of myogenic differentiation that regulates the calcium signaling pathway and AS of ROCK2. Overall, our results highlight the critical role of RUNX1T1 in myogenesis and broaden our understanding of myogenic differentiation.

The oxoglutarate dehydrogenase complex is involved in myofibril growth and Z-disc assembly in Drosophila

Myofibrils are long intracellular cables specific to muscles, composed mainly of actin and myosin filaments. The actin and myosin filaments are organized into repeated units called sarcomeres, which form the myofibrils. Muscle contraction is achieved by the simultaneous shortening of sarcomeres, which requires all sarcomeres to be the same size. Muscles have a variety of ways to ensure sarcomere homogeneity. We have previously shown that the controlled oligomerization of Zasp proteins sets the diameter of the myofibril. Here, we looked for Zasp-binding proteins at the Z-disc to identify additional proteins coordinating myofibril growth and assembly. We found that the E1 subunit of the oxoglutarate dehydrogenase complex localizes to both the Z-disc and the mitochondria, and is recruited to the Z-disc by Zasp52. The three subunits of the oxoglutarate dehydrogenase complex are required for myofibril formation. Using super-resolution microscopy, we revealed the overall organization of the complex at the Z-disc. Metabolomics identified an amino acid imbalance affecting protein synthesis as a possible cause of myofibril defects, which is supported by OGDH-dependent localization of ribosomes at the Z-disc.

Mechanism of elastic energy storage of honey bee abdominal muscles under stress relaxation

Energy storage of passive muscles plays an important part in frequent activities of honey bee abdomens due to the muscle distribution and open circulatory system. However, the elastic energy and mechanical properties of structure in passive muscles remain unclear. In this article, stress relaxation tests on passive muscles from the terga of the honey bee abdomens were performed under different concentrations of blebbistatin and motion parameters. In stress relaxation, the load drop with the rapid and slow stages depending on stretching velocity and stretching length reflects the features of myosin-titin series structure and cross-bridge-actin cyclic connections in muscles. Then a model with 2 parallel modules based on the 2 feature structures in muscles was thus developed. The model described the stress relaxation and stretching of passive muscles from honey bee abdomen well for a good fitting in stress relaxation and verification in loading process. In addition, the stiffness change of cross-bridge under different concentrations of blebbistatin is obtained from the model. We derived the elastic deformation of cross-bridge and the partial derivatives of energy expressions on motion parameters from this model, which accorded the experimental results. This model reveals the mechanism of passive muscles from honey bee abdomens suggesting that the temporary energy storage of cross-bridge in terga muscles under abdomen bending provides potential energy for springback during the periodic abdomen bending of honey bee or other arthropod insects. The finding also provides an experimental and theoretical basis for the novel microstructure and material design of bionic muscle.

Comprehensive analysis of pre-mRNA alternative splicing regulated by m6A methylation in pig oxidative and glycolytic skeletal muscles

BACKGROUND

Different types of skeletal myofibers exhibit distinct physiological and metabolic properties that are associated with meat quality traits in livestock. Alternative splicing (AS) of pre-mRNA can generate multiple transcripts from an individual gene by differential selection of splice sites. N6-methyladenosine (m6A) is the most abundant modification in mRNAs, but its regulation for AS in different muscles remains unknown.  RESULTS: We characterized AS events and m6A methylation pattern in pig oxidative and glycolytic muscles. A tota1 of 1294 differential AS events were identified, and differentially spliced genes were significantly enriched in processes related to different phenotypes between oxidative and glycolytic muscles. We constructed the regulatory network between splicing factors and corresponding differential AS events and identified NOVA1 and KHDRBS2 as key splicing factors. AS event was enriched in m6A-modified genes, and the methylation level was positively correlated with the number of AS events in genes. The dynamic change in m6A enrichment was associated with 115 differentially skipping exon (SE-DAS) events within 92 genes involving in various processes, including muscle contraction and myofibril assembly. We obtained 23.4% SE-DAS events (27/115) regulated by METTL3-meditaed m6A and experimentally validated the aberrant splicing of ZNF280D, PHE4DIP, and NEB. The inhibition of m6A methyltransferase METTL3 could induce the conversion of oxidative fiber to glycolytic fiber in PSCs.

CONCLUSION

Our study suggested that m6A modification could contribute to significant difference in phenotypes between oxidative and glycolytic muscles by mediating the regulation of AS. These findings would provide novel insights into mechanisms underlying muscle fiber conversion.

The insect perspective on Z-disc structure and biology

Myofibrils are the intracellular structures formed by actin and myosin filaments. They are paracrystalline contractile cables with unusually well-defined dimensions. The sliding of actin past myosin filaments powers contractions, and the entire system is held in place by a structure called the Z-disc, which anchors the actin filaments. Myosin filaments, in turn, are anchored to another structure called the M-line. Most of the complex architecture of myofibrils can be reduced to studying the Z-disc, and recently, important advances regarding the arrangement and function of Z-discs in insects have been published. On a very small scale, we have detailed protein structure information. At the medium scale, we have cryo-electron microscopy maps, super-resolution microscopy and protein-protein interaction networks, while at the functional scale, phenotypic data are available from precise genetic manipulations. All these data aim to answer how the Z-disc works and how it is assembled. Here, we summarize recent data from insects and explore how it fits into our view of the Z-disc, myofibrils and, ultimately, muscles.

Stretching muscle cells induces transcriptional and splicing transitions and changes in SR proteins

Alternative splicing is an RNA processing mechanism involved in skeletal muscle development and pathology. Muscular diseases exhibit splicing alterations and changes in mechanobiology leading us to investigate the interconnection between mechanical forces and RNA processing. We performed deep RNA-sequencing after stretching muscle cells. First, we uncovered transcriptional changes in genes encoding proteins involved in muscle function and transcription. Second, we observed that numerous mechanosensitive genes were part of the MAPK pathway which was activated in response to stretching. Third, we revealed that stretching skeletal muscle cells increased the proportion of alternatively spliced cassette exons and their inclusion. Fourth, we demonstrated that the serine and arginine-rich proteins exhibited stronger transcriptional changes than other RNA-binding proteins and that SRSF4 phosphorylation is mechanosensitive. Identifying SRSF4 as a mechanosensitive RNA-binding protein that might contribute to crosstalk between mechanotransduction, transcription, and splicing could potentially reveal novel insights into muscular diseases, particularly those with unknown etiologies.

Hereditary Basis of Coat Color and Excellent Feed Conversion Rate of Red Angus Cattle by Next-Generation Sequencing Data

Angus cattle have made remarkable contributions to the livestock industry worldwide as a commercial meat-type breed. Some evidence supported that Angus cattle with different coat colors have different feed-to-meat ratios, and the genetic basis of their coat color is inconclusive. Here, genome-wide association study was performed to investigate the genetic divergence of black and red Angus cattle with 63 public genome sequencing data. General linear model analysis was used to identify genomic regions with potential candidate variant/genes that contribute to coat color and feed conversion rate. Results showed that six single nucleotide polymorphisms (SNPs) and two insertion−deletions, which were annotated in five genes (ZCCHC14, ANKRD11, FANCA, MC1R, and LOC532875 [AFG3-like protein 1]), considerably diverged between black and red Angus cattle. The strongest associated loci, namely, missense mutation CHIR18_14705671 (c.296T > C) and frameshift mutation CHIR18_12999497 (c.310G>-), were located in MC1R. Three consecutive strongly associated SNPs were also identified and located in FANCA, which is widely involved in the Fanconi anemia pathway. Several SNPs of highly associated SNPs was notably enriched in ZCCHC14 and ANKRD11, which are related to myofiber growth and muscle development. This study provides a basis for the use of potential genetic markers to be used in future breeding programs to improve cattle selection in terms of coat color and meat phenotype. This study is also helpful to understand the hereditary basis of different coat colors and meat phenotypes. However, the putative candidate genes or markers identified in this study require further investigation to confirm their phenotypic causality and potential effective genetic relationships.

Rbfox1 is required for myofibril development and maintaining fiber type-specific isoform expression in Drosophila muscles

Protein isoform transitions confer muscle fibers with distinct properties and are regulated by differential transcription and alternative splicing. RNA-binding Fox protein 1 (Rbfox1) can affect both transcript levels and splicing, and is known to contribute to normal muscle development and physiology in vertebrates, although the detailed mechanisms remain obscure. In this study, we report that Rbfox1 contributes to the generation of adult muscle diversity in Drosophila Rbfox1 is differentially expressed among muscle fiber types, and RNAi knockdown causes a hypercontraction phenotype that leads to behavioral and eclosion defects. Misregulation of fiber type-specific gene and splice isoform expression, notably loss of an indirect flight muscle-specific isoform of Troponin-I that is critical for regulating myosin activity, leads to structural defects. We further show that Rbfox1 directly binds the 3'-UTR of target transcripts, regulates the expression level of myogenic transcription factors myocyte enhancer factor 2 and Salm, and both modulates expression of and genetically interacts with the CELF family RNA-binding protein Bruno1 (Bru1). Rbfox1 and Bru1 co-regulate fiber type-specific alternative splicing of structural genes, indicating that regulatory interactions between FOX and CELF family RNA-binding proteins are conserved in fly muscle. Rbfox1 thus affects muscle development by regulating fiber type-specific splicing and expression dynamics of identity genes and structural proteins.

Diversification of muscle types in Drosophila embryos

Drosophila embryonic somatic muscles represent a simple and tractable model system to study the gene regulatory networks that control diversification of cell types. Somatic myogenesis in Drosophila is initiated by intrinsic action of the mesodermal master gene twist, which activates a cascade of transcriptional outputs including myogenic differentiation factor Mef2, which triggers all aspects of the myogenic differentiation program. In parallel, the expression of a combinatorial code of identity transcription factors (iTFs) defines discrete particular features of each muscle fiber, such as number of fusion events, and specific attachment to tendon cells or innervation, thus ensuring diversification of muscle types. Here, we take the example of a subset of lateral transverse (LT) muscles and discuss how the iTF code and downstream effector genes progressively define individual LT properties such as fusion program, attachment and innervation. We discuss new challenges in the field including the contribution of posttranscriptional and epitranscriptomic regulation of gene expression in the diversification of cell types.

The Hox transcription factor Ultrabithorax binds RNA and regulates co-transcriptional splicing through an interplay with RNA polymerase II

Transcription factors (TFs) play a pivotal role in cell fate decision by coordinating gene expression programs. Although most TFs act at the DNA layer, few TFs bind RNA and modulate splicing. Yet, the mechanistic cues underlying TFs activity in splicing remain elusive. Focusing on the Drosophila Hox TF Ultrabithorax (Ubx), our work shed light on a novel layer of Ubx function at the RNA level. Transcriptome and genome-wide binding profiles in embryonic mesoderm and Drosophila cells indicate that Ubx regulates mRNA expression and splicing to promote distinct outcomes in defined cellular contexts. Our results demonstrate a new RNA-binding ability of Ubx. We find that the N51 amino acid of the DNA-binding Homeodomain is non-essential for RNA interaction in vitro, but is required for RNA interaction in vivo and Ubx splicing activity. Moreover, mutation of the N51 amino acid weakens the interaction between Ubx and active RNA Polymerase II (Pol II). Our results reveal that Ubx regulates elongation-coupled splicing, which could be coordinated by a dynamic interplay with active Pol II on chromatin. Overall, our work uncovered a novel role of the Hox TFs at the mRNA regulatory layer. This could be an essential function for other classes of TFs to control cell diversity.

Transcriptome-Wide Analyses Identify Dominant as the Predominantly Non-Conservative Alternative Splicing Inheritance Patterns in F1 Chickens

Transcriptome analysis has been used to investigate many economically traits in chickens; however, alternative splicing still lacks a systematic method of study that is able to promote proteome diversity, and fine-tune expression dynamics. Hybridization has been widely utilized in chicken breeding due to the resulting heterosis, but the dynamic changes in alternative splicing during this process are significant yet unclear. In this study, we performed a reciprocal crossing experiment involving the White Leghorn and Cornish Game chicken breeds which exhibit major differences in body size and reproductive traits, and conducted RNA sequencing of the brain, muscle, and liver tissues to identify the inheritance patterns. A total of 40 515 and 42 612 events were respectively detected in the brain and muscle tissues, with 39 843 observed in the liver; 2807, 4242, and 4538 events significantly different between two breeds were identified in the brain, muscle, and liver tissues, respectively. The hierarchical cluster of tissues from different tissues from all crosses, based on the alternative splicing profiles, suggests high tissue and strain specificity. Furthermore, a comparison between parental strains and hybrid crosses indicated that over one third of alternative splicing genes showed conserved patterns in all three tissues, while the second prevalent pattern was non-additive, which included both dominant and transgressive patterns; this meant that the dominant pattern plays a more important role than suppression. Our study provides an overview of the inheritance patterns of alternative splicing in layer and broiler chickens, to better understand post-transcriptional regulation during hybridization.

Mechanobiology of muscle and myofibril morphogenesis

Muscles generate forces for animal locomotion. The contractile apparatus of muscles is the sarcomere, a highly regular array of large actin and myosin filaments linked by gigantic titin springs. During muscle development many sarcomeres assemble in series into long periodic myofibrils that mechanically connect the attached skeleton elements. Thus, ATP-driven myosin forces can power movement of the skeleton. Here we review muscle and myofibril morphogenesis, with a particular focus on their mechanobiology. We describe recent progress on the molecular structure of sarcomeres and their mechanical connections to the skeleton. We discuss current models predicting how tension coordinates the assembly of key sarcomeric components to periodic myofibrils that then further mature during development. This requires transcriptional feedback mechanisms that may help to coordinate myofibril assembly and maturation states with the transcriptional program. To fuel the varying energy demands of muscles we also discuss the close mechanical interactions of myofibrils with mitochondria and nuclei to optimally support powerful or enduring muscle fibers.

RNA-Binding Proteins in the Post-transcriptional Control of Skeletal Muscle Development, Regeneration and Disease

Embryonic myogenesis is a temporally and spatially regulated process that generates skeletal muscle of the trunk and limbs. During this process, mononucleated myoblasts derived from myogenic progenitor cells within the somites undergo proliferation, migration and differentiation to elongate and fuse into multinucleated functional myofibers. Skeletal muscle is the most abundant tissue of the body and has the remarkable ability to self-repair by re-activating the myogenic program in muscle stem cells, known as satellite cells. Post-transcriptional regulation of gene expression mediated by RNA-binding proteins is critically required for muscle development during embryogenesis and for muscle homeostasis in the adult. Differential subcellular localization and activity of RNA-binding proteins orchestrates target gene expression at multiple levels to regulate different steps of myogenesis. Dysfunctions of these post-transcriptional regulators impair muscle development and homeostasis, but also cause defects in motor neurons or the neuromuscular junction, resulting in muscle degeneration and neuromuscular disease. Many RNA-binding proteins, such as members of the muscle blind-like (MBNL) and CUG-BP and ETR-3-like factors (CELF) families, display both overlapping and distinct targets in muscle cells. Thus they function either cooperatively or antagonistically to coordinate myoblast proliferation and differentiation. Evidence is accumulating that the dynamic interplay of their regulatory activity may control the progression of myogenic program as well as stem cell quiescence and activation. Moreover, the role of RNA-binding proteins that regulate post-transcriptional modification in the myogenic program is far less understood as compared with transcription factors involved in myogenic specification and differentiation. Here we review past achievements and recent advances in understanding the functions of RNA-binding proteins during skeletal muscle development, regeneration and disease, with the aim to identify the fundamental questions that are still open for further investigations.

A Candidate RNAi Screen Reveals Diverse RNA-Binding Protein Phenotypes in Drosophila Flight Muscle

The proper regulation of RNA processing is critical for muscle development and the fine-tuning of contractile ability among muscle fiber-types. RNA binding proteins (RBPs) regulate the diverse steps in RNA processing, including alternative splicing, which generates fiber-type specific isoforms of structural proteins that confer contractile sarcomeres with distinct biomechanical properties. Alternative splicing is disrupted in muscle diseases such as myotonic dystrophy and dilated cardiomyopathy and is altered after intense exercise as well as with aging. It is therefore important to understand splicing and RBP function, but currently, only a small fraction of the hundreds of annotated RBPs expressed in muscle have been characterized. Here, we demonstrate the utility of Drosophila as a genetic model system to investigate basic developmental mechanisms of RBP function in myogenesis. We find that RBPs exhibit dynamic temporal and fiber-type specific expression patterns in mRNA-Seq data and display muscle-specific phenotypes. We performed knockdown with 105 RNAi hairpins targeting 35 RBPs and report associated lethality, flight, myofiber and sarcomere defects, including flight muscle phenotypes for Doa, Rm62, mub, mbl, sbr, and clu. Knockdown phenotypes of spliceosome components, as highlighted by phenotypes for A-complex components SF1 and Hrb87F (hnRNPA1), revealed level- and temporal-dependent myofibril defects. We further show that splicing mediated by SF1 and Hrb87F is necessary for Z-disc stability and proper myofibril development, and strong knockdown of either gene results in impaired localization of kettin to the Z-disc. Our results expand the number of RBPs with a described phenotype in muscle and underscore the diversity in myofibril and transcriptomic phenotypes associated with splicing defects. Drosophila is thus a powerful model to gain disease-relevant insight into cellular and molecular phenotypes observed when expression levels of splicing factors, spliceosome components and splicing dynamics are altered.

Benzo[a]pyrene diol epoxide-induced transformed cells identify the significance of hsa_circ_0051488, a ERCC1-derived circular RNA in pulmonary squamous cell carcinoma

ERCC1 is a gene for repairing DNA damage whose function is related to carcinogenic-induced tumorigenesis and the effectiveness of platinum therapies. Circular RNAs (circRNAs) are products of posttranscriptional regulation with pleiotropic effects on the pathogenesis of lung cancer. We aim to identify that specific circRNAs derived from ERCC1 can regulate key biological processes involved in the development of lung cancer. We performed bioinformatics analysis, in vitro experiments, and analyzed clinical samples, to determine the biological features of a certain ERCC1-derived circRNA termed as hsa_circ_0051488 in benzo[a]pyrene diol epoxide-induced malignant transformed cell and lung cancer cell. The well-established model of transformed cells provided an ideal platform for analyzing the molecular characteristics of this circRNA in the malignant transformation of lung epithelial cell, which supports that hsa_circ_0051488 functions in the onset and growth of lung squamous cell carcinoma (LUSC). Further analysis indicates that the absence of hsa_circ_0051488 promoted the proliferation of cells with the malignant phenotype. Extensive experiments confirm that hsa_circ_0051488 is present in the cytoplasm and functioned as a competing endogenous RNA. In particular, hsa_circ_0051488 binds to mir-6717-5p, thereby modulating the expression of SATB2 gene, a lung cancer suppressor. Furthermore, our in silico experiments indicate that SATB2 can inhibit multiple tumor pathways and its expression positively correlated with the tumor suppressor gene CRMP1. These findings suggest a possible regulatory mechanism of hsa_circ_0051488 in LUSC, and that the newly discovered hsa_circ_0051488/miR-6717-5p/SATB2 axis may be a potential route for therapeutic intervention of LUSC.

CRISPR gene editing in pluripotent stem cells reveals the function of MBNL proteins during human in vitro myogenesis

Alternative splicing has emerged as a fundamental mechanism for the spatiotemporal control of development. A better understanding of how this mechanism is regulated has the potential not only to elucidate fundamental biological principles, but also to decipher pathological mechanisms implicated in diseases where normal splicing networks are misregulated. Here, we took advantage of human pluripotent stem cells to decipher during human myogenesis the role of muscleblind-like (MBNL) proteins, a family of tissue-specific splicing regulators whose loss of function is associated with myotonic dystrophy type 1 (DM1), an inherited neuromuscular disease. Thanks to the CRISPR/Cas9 technology, we generated human-induced pluripotent stem cells (hiPSCs) depleted in MBNL proteins and evaluated the consequences of their losses on the generation of skeletal muscle cells. Our results suggested that MBNL proteins are required for the late myogenic maturation. In addition, loss of MBNL1 and MBNL2 recapitulated the main features of DM1 observed in hiPSC-derived skeletal muscle cells. Comparative transcriptomic analyses also revealed the muscle-related processes regulated by these proteins that are commonly misregulated in DM1. Together, our study reveals the temporal requirement of MBNL proteins in human myogenesis and should facilitate the identification of new therapeutic strategies capable to cope with the loss of function of these MBNL proteins.

Myofibril and mitochondria morphogenesis are coordinated by a mechanical feedback mechanism in muscle

Complex animals build specialised muscles to match specific biomechanical and energetic needs. Hence, composition and architecture of sarcomeres and mitochondria are muscle type specific. However, mechanisms coordinating mitochondria with sarcomere morphogenesis are elusive. Here we use Drosophila muscles to demonstrate that myofibril and mitochondria morphogenesis are intimately linked. In flight muscles, the muscle selector spalt instructs mitochondria to intercalate between myofibrils, which in turn mechanically constrain mitochondria into elongated shapes. Conversely in cross-striated leg muscles, mitochondria networks surround myofibril bundles, contacting myofibrils only with thin extensions. To investigate the mechanism causing these differences, we manipulated mitochondrial dynamics and found that increased mitochondrial fusion during myofibril assembly prevents mitochondrial intercalation in flight muscles. Strikingly, this causes the expression of cross-striated muscle specific sarcomeric proteins. Consequently, flight muscle myofibrils convert towards a partially cross-striated architecture. Together, these data suggest a biomechanical feedback mechanism downstream of spalt synchronizing mitochondria with myofibril morphogenesis.

The Emerging Role of the RBM20 and PTBP1 Ribonucleoproteins in Heart Development and Cardiovascular Diseases

Alternative splicing is a regulatory mechanism essential for cell differentiation and tissue organization. More than 90% of human genes are regulated by alternative splicing events, which participate in cell fate determination. The general mechanisms of splicing events are well known, whereas only recently have deep-sequencing, high throughput analyses and animal models provided novel information on the network of functionally coordinated, tissue-specific, alternatively spliced exons. Heart development and cardiac tissue differentiation require thoroughly regulated splicing events. The ribonucleoprotein RBM20 is a key regulator of the alternative splicing events required for functional and structural heart properties, such as the expression of TTN isoforms. Recently, the polypyrimidine tract-binding protein PTBP1 has been demonstrated to participate with RBM20 in regulating splicing events. In this review, we summarize the updated knowledge relative to RBM20 and PTBP1 structure and molecular function; their role in alternative splicing mechanisms involved in the heart development and function; RBM20 mutations associated with idiopathic dilated cardiovascular disease (DCM); and the consequences of RBM20-altered expression or dysfunction. Furthermore, we discuss the possible application of targeting RBM20 in new approaches in heart therapies.