Bi-allelic mutations in MYL1 cause a severe congenital myopathy

Gene regulatory network analysis identifies MYL1, MDH2, GLS, and TRIM28 as the principal proteins in the response of mesenchymal stem cells to Mg2+ ions

Magnesium (Mg)-based implants have emerged as a promising alternative for orthopedic applications, owing to their bioactive properties and biodegradability. As the implants degrade, Mg2+ ions are released, influencing all surrounding cell types, especially mesenchymal stem cells (MSCs). MSCs are vital for bone tissue regeneration, therefore, it is essential to understand their molecular response to Mg2+ ions in order to maximize the potential of Mg-based biomaterials. In this study, we conducted a gene regulatory network (GRN) analysis to examine the molecular responses of MSCs to Mg2+ ions. We used time-series proteomics data collected at 11 time points across a 21-day period for the GRN construction. We studied the impact of Mg2+ ions on the resulting networks and identified the key proteins and protein interactions affected by the application of Mg2+ ions. Our analysis highlights MYL1, MDH2, GLS, and TRIM28 as the primary targets of Mg2+ ions in the response of MSCs during 1-21 days phase. Our results also identify MDH2-MYL1, MDH2-RPS26, TRIM28-AK1, TRIM28-SOD2, and GLS-AK1 as the critical protein relationships affected by Mg2+ ions. By offering a comprehensive understanding of the regulatory role of Mg2+ ions on MSCs, our study contributes valuable insights into the molecular response of MSCs to Mg-based materials, thereby facilitating the development of innovative therapeutic strategies for orthopedic applications.

A knowledge graph approach to predict and interpret disease-causing gene interactions

BACKGROUND

Understanding the impact of gene interactions on disease phenotypes is increasingly recognised as a crucial aspect of genetic disease research. This trend is reflected by the growing amount of clinical research on oligogenic diseases, where disease manifestations are influenced by combinations of variants on a few specific genes. Although statistical machine-learning methods have been developed to identify relevant genetic variant or gene combinations associated with oligogenic diseases, they rely on abstract features and black-box models, posing challenges to interpretability for medical experts and impeding their ability to comprehend and validate predictions. In this work, we present a novel, interpretable predictive approach based on a knowledge graph that not only provides accurate predictions of disease-causing gene interactions but also offers explanations for these results.

RESULTS

We introduce BOCK, a knowledge graph constructed to explore disease-causing genetic interactions, integrating curated information on oligogenic diseases from clinical cases with relevant biomedical networks and ontologies. Using this graph, we developed a novel predictive framework based on heterogenous paths connecting gene pairs. This method trains an interpretable decision set model that not only accurately predicts pathogenic gene interactions, but also unveils the patterns associated with these diseases. A unique aspect of our approach is its ability to offer, along with each positive prediction, explanations in the form of subgraphs, revealing the specific entities and relationships that led to each pathogenic prediction.

CONCLUSION

Our method, built with interpretability in mind, leverages heterogenous path information in knowledge graphs to predict pathogenic gene interactions and generate meaningful explanations. This not only broadens our understanding of the molecular mechanisms underlying oligogenic diseases, but also presents a novel application of knowledge graphs in creating more transparent and insightful predictors for genetic research.

Circ2388 regulates myogenesis and muscle regeneration

The formation of skeletal muscle is a complex process that is coordinated by many regulatory factors, such as myogenic factors and noncoding RNAs. Numerous studies have proved that circRNA is an indispensable part of muscle development. However, little is known about circRNAs in bovine myogenesis. In this study, we discovered a novel circRNA, circ2388, formed by reverse splicing of the fourth and fifth exons of the MYL1 gene. The expression of circ2388 was different between fetal and adult cattle muscle. This circRNA is 99% homologous between cattle and buffalo and is localized in the cytoplasm. Thoroughly, we proved that circ2388 had no effect on cattle and buffalo myoblast proliferation but promotes myoblast differentiation and myotube fusion. Furthermore, circ2388 in vivo stimulated skeletal muscle regeneration in mouse muscle injury model. Taken together, our findings suggest that circ2388 promotes myoblast differentiation and promotes the recovery and regeneration of damaged muscles.

Decabromodiphenyl Ethane Mainly Affected the Muscle Contraction and Reproductive Endocrine System in Female Adult Zebrafish

The novel brominated flame retardant decabromodiphenyl ethane (DBDPE) has become a widespread environmental pollutant. However, the target tissue and toxicity of DBDPE are still not clear. In the current study, female zebrafish were exposed to 1 and 100 nM DBDPE for 28 days. Chemical analysis revealed that DBDPE tended to accumulate in the brain other than the liver and gonad. Subsequently, tandem mass tag-based quantitative proteomics and parallel reaction monitoring verification were performed to screen the differentially expressed proteins in the brain. Bioinformatics analysis revealed that DBDPE mainly affected the biological process related to muscle contraction and estrogenic response. Therefore, the neurotoxicity and reproductive disruptions were validated via multilevel toxicological endpoints. Specifically, locomotor behavioral changes proved the potency of neurotoxicity, which may be caused by disturbance of muscular proteins and calcium homeostasis; decreases of sex hormone levels and transcriptional changes of genes related to the hypothalamic-pituitary-gonad-liver axis confirmed reproductive disruptions upon DBDPE exposure. In summary, our results suggested that DBDPE primarily accumulated in the brain and evoked neurotoxicity and reproductive disruptions in female zebrafish. These findings can provide important clues for a further mechanism study and risk assessment of DBDPE.

Single-cell transcriptome analysis of the zebrafish embryonic trunk

During embryonic development, cells differentiate into a variety of distinct cell types and subtypes with diverse transcriptional profiles. To date, transcriptomic signatures of different cell lineages that arise during development have been only partially characterized. Here we used single-cell RNA-seq to perform transcriptomic analysis of over 20,000 cells disaggregated from the trunk region of zebrafish embryos at the 30 hpf stage. Transcriptional signatures of 27 different cell types and subtypes were identified and annotated during this analysis. This dataset will be a useful resource for many researchers in the fields of developmental and cellular biology and facilitate the understanding of molecular mechanisms that regulate cell lineage choices during development.

Reorganization of chromatin architecture during prenatal development of porcine skeletal muscle

Myofibres (primary and secondary myofibre) are the basic structure of muscle and the determinant of muscle mass. To explore the skeletal muscle developmental processes from primary myofibres to secondary myofibres in pigs, we conducted an integrative three-dimensional structure of genome and transcriptomic characterization of longissimus dorsi muscle of pig from primary myofibre formation stage [embryonic Day 35 (E35)] to secondary myofibre formation stage (E80). In the hierarchical genomic structure, we found that 11.43% of genome switched compartment A/B status, 14.53% of topologically associating domains are changed intradomain interactions (D-scores) and 2,730 genes with differential promoter–enhancer interactions and (or) enhancer activity from E35 to E80. The alterations of genome architecture were found to correlate with expression of genes that play significant roles in neuromuscular junction, embryonic morphogenesis, skeletal muscle development or metabolism, typically, NEFL, MuSK, SLN, Mef2D and GCK. Significantly, Sox6 and MATN2 play important roles in the process of primary to secondary myofibres formation and increase the regulatory potential score and genes expression in it. In brief, we reveal the genomic reorganization from E35 to E80 and construct genome-wide high-resolution interaction maps that provide a resource for studying long-range control of gene expression from E35 to E80.

A ceRNA-associated risk model predicts the poor prognosis for head and neck squamous cell carcinoma patients

Head and neck squamous cell carcinoma (HNSCC) is one of the most malignant cancers with poor prognosis worldwide. Emerging evidence indicates that competing endogenous RNAs (ceRNAs) are involved in various diseases, however, the regulatory mechanisms of ceRNAs underlying HNSCC remain unclear. In this study, we retrieved differentially expressed long non-coding RNAs (DElncRNAs), messenger RNAs (DEmRNAs) and microRANs (DEmiRNAs) from The Cancer Genome Atlas database and constructed a ceRNA-based risk model in HNSCC by integrated bioinformatics approaches. Functional enrichment analyses showed that DEmRNAs might be involved in extracellular matrix related biological processes, and protein–protein interaction network further selected out prognostic genes, including MYL1 and ACTN2. Importantly, co-expressed RNAs identified by weighted co-expression gene network analysis constructed the ceRNA networks. Moreover, AC114730.3, AC136375.3, LAT and RYR3 were highly correlated to overall survival of HNSCC by Kaplan–Meier method and univariate Cox regression analysis, which were subsequently implemented multivariate Cox regression analysis to build the risk model. Our study provides a deeper understanding of ceRNAs on the regulatory mechanisms, which will facilitate the expansion of the roles on the ceRNAs in the tumorigenesis, development and treatment of HNSCC.

MYBPC2 and MYL1 as Significant Gene Markers for Rhabdomyosarcoma

BACKGROUND

Rhabdomyosarcoma is the most common soft tissue tumor in children. Rhabdomyosarcoma commonly results in pain and bleeding caused by tumor compression and is prone to early metastasis and recurrence, which can seriously affect the therapeutic outcomes and long-term prognosis. Up to 37.7% of rhabdomyosarcomas may metastasize. Therefore, the molecular mechanisms underlying rhabdomyosarcoma must be explored to identify an effective target for its early diagnosis and specific treatment.

METHODS

A dataset of 18 rhabdomyosarcoma tissue samples and 6 healthy skeletal muscle samples was downloaded. Differentially expressed genes between rhabdomyosarcoma and healthy tissue samples were identified by GEO2R. Kyoto Encyclopedia of Genes and Genomes and gene ontology pathway enrichment analyses were performed. A protein-protein interaction network was constructed, and hub genes were identified. Expression and survival analyses of hub genes were performed. Additionally, 30 patients with rhabdomyosarcoma were recruited, and overall survival information and samples were collected. Reverse transcription quantitative real-time polymerase chain reaction assays were performed to verify the expression of MYBPC2 and MYL1 in rhabdomyosarcoma tumor tissues. The Kaplan-Meier method was used to explore overall survival based on our clinical data.

RESULTS

In total, 164 genes were up-regulated and 394 were down-regulated in rhabdomyosarcoma tumor tissues. Gene ontology analysis revealed that variations were predominantly enriched in the cell cycle, muscle contraction, muscle system processes, cytoskeleton, nucleotide binding, and cytoskeletal protein binding. The protein-protein interaction network revealed 3274 edges, and 441 nodes were constructed. Ten hub genes were identified; of these, MYBPC2 and MYL1 were significantly up-regulated in rhabdomyosarcoma. Compared with the healthy group, patients with rhabdomyosarcoma exhibiting high expression of MYBPC2 and MYL1 exhibited significantly worse overall survival.

CONCLUSIONS

We found differentially expressed genes between rhabdomyosarcoma and healthy tissue samples. MYBPC2 and MYL1 may be involved in the pathogenesis of rhabdomyosarcoma and therefore deserve further exploration.

Circular RNA circMYL1 Inhibit Proliferation and Promote Differentiation of Myoblasts by Sponging miR-2400

Circular RNAs (circRNAs) are a class of endogenous non-coding RNAs (ncRNAs) involved in regulating skeletal muscle development by sponging miRNAs. In this study, we found that the circMYL1 expression was down-regulated during myoblast proliferation, while gradually up-regulated in myoblast differentiation. The potential role of circMYL1 was identified in the proliferation of bovine myoblast through mRNA and protein expression of proliferation marker genes (PCNA, CyclinD1, and CDK2), cell counting kit-8 assay, flow cytometry analysis, and 5-ethynyl 2′-deoxyuridine (EdU) assay. Analysis of the expression of differentiation marker genes (MyoD, MyoG, and MYH2) and immunofluorescence of Myosin heavy chain (MyHC) was used to assess cell differentiation. The proliferation analysis revealed that circMYL1 inhibited the proliferation of bovine primary myoblast. Furthermore, the differentiation analysis demonstrated that circMYL1 promoted the differentiation of bovine primary myoblast. The luciferase screening and RNA immunoprecipitation (RIP) assays found that circMYL1 could have interaction with miR-2400. Additionally, we demonstrated that miR-2400 promoted proliferation and inhibited differentiation of bovine primary myoblast, while circMYL1 may eliminate the effects of miR-2400, as showed by rescue experiments. Together, our results revealed that a novel circular RNA of circMYL1 could inhibit proliferation and promote differentiation of myoblast by sponging miR-2400.

Identification and Characterization of circRNA in Longissimus Dorsi of Different Breeds of Cattle

Shandong black cattle is a new breed of cattle that is developed by applying modern biotechnology, such as somatic cloning, and conventional breeding methods to Luxi cattle. It is very important to study the function and regulatory mechanism of circRNAs in muscle differentiation among different breeds to improve meat quality and meat production performance and to provide new ideas for beef cattle meat quality improvements and new breed development. Therefore, the goal of this study was to sequence and identify circRNAs in muscle tissues of different breeds of cattle. We used RNA-seq to identify circRNAs in the muscles of two breeds of cattle (Shandong black and Luxi). We identified 14,640 circRNAs and found 655 differentially expressed circRNAs. We also analyzed the classification and characteristics of circRNAs in muscle tissue. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were used on the parental genes of circRNAs. They were mainly involved in a variety of biological processes, such as muscle fiber development, smooth muscle cell proliferation, bone system morphogenesis, tight junctions and the MAPK, AMPK, and mTOR signaling pathways. In addition, we used miRanda to predict the interactions between 14 circRNAs and 11 miRNAs. Based on the above assays, we identified circRNAs (circ0001048, circ0001103, circ0001159, circ0003719, circ0003424, circ0003721, circ0003720, circ0001519, circ0001530, circ0005011, circ0014518, circ0000181, circ0000190, circ0010558) that may play important roles in the regulation of muscle growth and development. Using real-time quantitative PCR, 14 circRNAs were randomly selected to verify the real circRNAs. Luciferase reporter gene system was used to verify the binding site of miR-1 in circ0014518. Our results provide more information about circRNAs regulating muscle development in different breeds of cattle and lay a solid foundation for future experiments.

Neurogenetic fetal akinesia and arthrogryposis: genetics, expanding genotype-phenotypes and functional genomics

BACKGROUND

Fetal akinesia and arthrogryposis are clinically and genetically heterogeneous and have traditionally been refractive to genetic diagnosis. The widespread availability of affordable genome-wide sequencing has facilitated accurate genetic diagnosis and gene discovery in these conditions.

METHODS

We performed next generation sequencing (NGS) in 190 probands with a diagnosis of arthrogryposis multiplex congenita, distal arthrogryposis, fetal akinesia deformation sequence or multiple pterygium syndrome. This sequencing was a combination of bespoke neurogenetic disease gene panels and whole exome sequencing. Only class 4 and 5 variants were reported, except for two cases where the identified variants of unknown significance (VUS) are most likely to be causative for the observed phenotype. Co-segregation studies and confirmation of variants identified by NGS were performed where possible. Functional genomics was performed as required.

RESULTS

Of the 190 probands, 81 received an accurate genetic diagnosis. All except two of these cases harboured class 4 and/or 5 variants based on the American College of Medical Genetics and Genomics guidelines. We identified phenotypic expansions associated with CACNA1S, CHRNB1, GMPPB and STAC3. We describe a total of 50 novel variants, including a novel missense variant in the recently identified gene for arthrogryposis with brain malformations-SMPD4.

CONCLUSIONS

Comprehensive gene panels give a diagnosis for a substantial proportion (42%) of fetal akinesia and arthrogryposis cases, even in an unselected cohort. Recently identified genes account for a relatively large proportion, 32%, of the diagnoses. Diagnostic-research collaboration was critical to the diagnosis and variant interpretation in many cases, facilitated genotype-phenotype expansions and reclassified VUS through functional genomics.

Targeted gene panel use in 2249 neuromuscular patients: the Australasian referral center experience

Objective

To develop, test, and iterate a comprehensive neuromuscular targeted gene panel in a national referral center.

Methods

We designed two iterations of a comprehensive targeted gene panel for neuromuscular disorders. Version 1 included 336 genes, which was increased to 464 genes in Version 2. Both panels used TargetSeq^TM probe‐based hybridization for target enrichment followed by Ion Torrent sequencing. Targeted high‐coverage sequencing and analysis was performed on 2249 neurology patients from Australia and New Zealand (1054 Version 1, 1195 Version 2) from 2012 to 2015. No selection criteria were used other than referral from a suitable medical specialist (e.g., neurologist or clinical geneticist). Patients were classified into 15 clinical categories based on the clinical diagnosis from the referring clinician.

Results

Six hundred and sixty‐five patients received a genetic diagnosis (30%). Diagnosed patients were significantly younger that undiagnosed patients (26.4 and 32.5 years, respectively; P = 4.6326E‐9). The diagnostic success varied markedly between disease categories. Pathogenic variants in 10 genes explained 38% of the disease burden. Unexpected phenotypic expansions were discovered in multiple cases. Triage of unsolved cases for research exome testing led to the discovery of six new disease genes.

Interpretation

A comprehensive targeted diagnostic panel was an effective method for neuromuscular disease diagnosis within the context of an Australasian referral center. Use of smaller disease‐specific panels would have precluded diagnosis in many patients and increased cost. Analysis through a centralized laboratory facilitated detection of recurrent, but under‐recognized pathogenic variants.

Altered miRNA and mRNA Expression in Sika Deer Skeletal Muscle with Age

Studies of the gene and miRNA expression profiles associated with the postnatal late growth, development, and aging of skeletal muscle are lacking in sika deer. To understand the molecular mechanisms of the growth and development of sika deer skeletal muscle, we used de novo RNA sequencing (RNA-seq) and microRNA sequencing (miRNA-seq) analyses to determine the differentially expressed (DE) unigenes and miRNAs from skeletal muscle tissues at 1, 3, 5, and 10 years in sika deer. A total of 51,716 unigenes, 171 known miRNAs, and 60 novel miRNAs were identified based on four mRNA and small RNA libraries. A total of 2,044 unigenes and 11 miRNAs were differentially expressed between adolescence and juvenile sika deer, 1,946 unigenes and 4 miRNAs were differentially expressed between adult and adolescent sika deer, and 2,209 unigenes and 1 miRNAs were differentially expressed between aged and adult sika deer. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that DE unigenes and miRNA were mainly related to energy and substance metabolism, processes that are closely associate with the growth, development, and aging of skeletal muscle. We also constructed mRNA–mRNA and miRNA–mRNA interaction networks related to the growth, development, and aging of skeletal muscle. The results show that mRNA (Myh1, Myh2, Myh7, ACTN3, etc.) and miRNAs (miR-133a, miR-133c, miR-192, miR-151-3p, etc.) may play important roles in muscle growth and development, and mRNA (WWP1, DEK, UCP3, FUS, etc.) and miRNAs (miR-17-5p, miR-378b, miR-199a-5p, miR-7, etc.) may have key roles in muscle aging. In this study, we determined the dynamic miRNA and unigenes transcriptome in muscle tissue for the first time in sika deer. The age-dependent miRNAs and unigenes identified will offer insights into the molecular mechanism underlying muscle development, growth, and maintenance and will also provide valuable information for sika deer genetic breeding.

Identification of eQTLs and sQTLs associated with meat quality in beef

BACKGROUND

Transcription has a substantial genetic control and genetic dissection of gene expression could help us understand the genetic architecture of complex phenotypes such as meat quality in cattle. The objectives of the present research were: 1) to perform eQTL and sQTL mapping analyses for meat quality traits in longissimus dorsi muscle; 2) to uncover genes whose expression is influenced by local or distant genetic variation; 3) to identify expression and splicing hot spots; and 4) to uncover genomic regions affecting the expression of multiple genes.

RESULTS

Eighty steers were selected for phenotyping, genotyping and RNA-seq evaluation. A panel of traits related to meat quality was recorded in longissimus dorsi muscle. Information on 112,042 SNPs and expression data on 8588 autosomal genes and 87,770 exons from 8467 genes were included in an expression and splicing quantitative trait loci (QTL) mapping (eQTL and sQTL, respectively). A gene, exon and isoform differential expression analysis previously carried out in this population identified 1352 genes, referred to as DEG, as explaining part of the variability associated with meat quality traits. The eQTL and sQTL mapping was performed using a linear regression model in the R package Matrix eQTL. Genotype and year of birth were included as fixed effects, and population structure was accounted for by including as a covariate the first PC from a PCA analysis on genotypic data. The identified QTLs were classified as cis or trans using 1 Mb as the maximum distance between the associated SNP and the gene being analyzed. A total of 8377 eQTLs were identified, including 75.6% trans, 10.4% cis, 12.5% DEG trans and 1.5% DEG cis; while 11,929 sQTLs were uncovered: 66.1% trans, 16.9% DEG trans, 14% cis and 3% DEG cis. Twenty-seven expression master regulators and 13 splicing master regulators were identified and were classified as membrane-associated or cytoskeletal proteins, transcription factors or DNA methylases. These genes could control the expression of other genes through cell signaling or by a direct transcriptional activation/repression mechanism.

CONCLUSION

In the present analysis, we show that eQTL and sQTL mapping makes possible positional identification of gene and isoform expression regulators.

High-resolution genome-wide expression analysis of single myofibers using SMART-Seq

Skeletal muscle is a heterogeneous tissue. Individual myofibers that make up muscle tissue exhibit variation in their metabolic and contractile properties. Although biochemical and histological assays are available to study myofiber heterogeneity, efficient methods to analyze the whole transcriptome of individual myofibers are lacking. Here, we report on a single-myofiber RNA-sequencing (smfRNA-Seq) approach to analyze the whole transcriptome of individual myofibers by combining single-fiber isolation with Switching Mechanism at 5' end of RNA Template (SMART) technology. Using smfRNA-Seq, we first determined the genes that are expressed in the whole muscle, including in nonmyogenic cells. We also analyzed the differences in the transcriptome of myofibers from young and old mice to validate the effectiveness of this new method. Our results suggest that aging leads to significant changes in the expression of metabolic genes, such as Nos1, and structural genes, such as Myl1, in myofibers. We conclude that smfRNA-Seq is a powerful tool to study developmental, disease-related, and age-related changes in the gene expression profile of skeletal muscle.

Insights into myosin regulatory and essential light chains: a focus on their roles in cardiac and skeletal muscle function, development and disease

The activity of cardiac and skeletal muscles depends upon the ATP-coupled actin-myosin interactions to execute the power stroke and muscle contraction. The goal of this review article is to provide insight into the function of myosin II, the molecular motor of the heart and skeletal muscles, with a special focus on the role of myosin II light chain (MLC) components. Specifically, we focus on the involvement of myosin regulatory (RLC) and essential (ELC) light chains in striated muscle development, isoform appearance and their function in normal and diseased muscle. We review the consequences of isoform switching and knockout of specific MLC isoforms on cardiac and skeletal muscle function in various animal models. Finally, we discuss how dysregulation of specific RLC/ELC isoforms can lead to cardiac and skeletal muscle diseases and summarize the effects of most studied mutations leading to cardiac or skeletal myopathies.

Recent advances in understanding congenital myopathies

By definition, congenital myopathy typically presents with skeletal muscle weakness and hypotonia at birth. Traditionally, congenital myopathy subtypes have been predominantly distinguished on the basis of the pathological hallmarks present on skeletal muscle biopsies. Many genes cause congenital myopathies when mutated, and a burst of new causative genes have been identified because of advances in gene sequencing technology. Recent discoveries include extending the disease phenotypes associated with previously identified genes and determining that genes formerly known to cause only dominant disease can also cause recessive disease. The more recently identified congenital myopathy genes account for only a small proportion of patients. Thus, the congenital myopathy genes remaining to be discovered are predicted to be extremely rare causes of disease, which greatly hampers their identification. Significant progress in the provision of molecular diagnoses brings important information and value to patients and their families, such as possible disease prognosis, better disease management, and informed reproductive choice, including carrier screening of parents. Additionally, from accurate genetic knowledge, rational treatment options can be hypothesised and subsequently evaluated in vitro and in animal models. A wide range of potential congenital myopathy therapies have been investigated on the basis of improved understanding of disease pathomechanisms, and some therapies are in clinical trials. Although large hurdles remain, promise exists for translating treatment benefits from preclinical models to patients with congenital myopathy, including harnessing proven successes for other genetic diseases.