*Corresponding Author:
Youfeng Zhou
Department of Pediatrics, Fujian Maternity and Child Health Hospital, College of Clinical Medicine for Obstetrics & Gynecology and Pediatrics, Fujian Medical University, Fuzhou, Fujian 350001, China
E-mail:
qwerrtyy014725@163.com
This article was originally published in a special issue, “New Research Outcomes in Drug and Health Sciences”
Indian J Pharm Sci 2023:85(6) Spl Issue “66-73”

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License, which allows others to remix, tweak, and build upon the work non-commercially, as long as the author is credited and the new creations are licensed under the identical terms

Abstract

The purpose of this study was to simulate muscular dystrophy disease by constructing a laminin-alpha 2 zebrafish model, and to study the changes of muscle transcriptome using transcriptomic methods. Using clustered regularly interspaced short palindromic repeats associated protein 9 gene editing technology; the zebrafish laminin-alpha 2 gene was knocked out or mutated. First, the appropriate guide RNA sequence was designed, and then the guide RNA and clustered regularly interspaced short palindromic repeats associated protein 9 proteins were synthesized and injected into zebrafish embryos. The lamininalpha 2 model fish is screened and identified by detecting the mutation of DNA or messenger RNA in the progeny fish genome. The muscle phenotype identification results of zebrafish model 4 showed that the muscle phenotypes of normal and unable swimming individuals in laminin-alpha 2 knockout group were abnormal, while the muscle phenotypes of clustered regularly interspaced short palindromic repeats associated protein 9 control group were not. The muscular dystrophy characteristics of laminin-alpha 2 model fish were confirmed. Functional enrichment and pathway analysis of differentially expressed genes were performed to further reveal the biological processes and regulatory pathways associated with muscular dystrophy in laminin-alpha 2 model fish. Transcriptome sequencing analysis revealed the difference in transcriptome expression between the muscle tissue of laminin-alpha 2 zebrafish model and wild-type zebrafish. These differentially expressed genes are involved in biological processes such as muscle development, metabolism and signal transduction, providing important clues for further understanding of the pathogenesis of muscular dystrophy. Through functional enrichment and pathway analysis, the key signaling pathways and molecular mechanisms involved in the pathogenesis of muscular dystrophy were identified. These findings provide theoretical basis for the discovery of new therapeutic targets and formulation of corresponding therapeutic strategies, and provide new ideas and methods for the treatment of muscle diseases.

Keywords

Laminin-alpha 2 zebrafish model, muscular dystrophy, transcriptomics, clustered regularly interspaced short palindromic repeats associated protein 9, metabolism

Muscle disease is a kind of disease that seriously affects human health and quality of life, its etiology and pathogenesis have not been fully clarified. Muscular dystrophy is a common muscle disease characterized by muscle weakness, loss of muscle mass, and abnormal muscle metabolism. Understanding the pathogenesis of muscle diseases is of great significance for diagnosis and treatment [1-3].

Muscle is one of the most important tissues in the human body and plays a key role in maintaining posture, movement ability and metabolic balance[4,5]. The occurrence of muscle diseases is related to a variety of factors, including genetic factors, environmental factors, metabolic abnormalities, neuromuscular connection disorders, etc.,[6-9]. Muscular dystrophy is one of the common types. Patients with muscular dystrophy often show symptoms such as decreased muscle strength, muscle atrophy, fatigue and uneven muscle hypertrophy, which seriously affect patients' quality of life[10-13].

However, the etiology and pathogenesis of muscle diseases are still not fully understood, which limits progress in their diagnosis and treatment. Traditional research methods are often restricted by problems such as difficulty in obtaining samples, limited technical means and limited research objects. Therefore, the development of new research methods and the establishment of suitable animal models become the key to the study of muscle diseases.

In recent years, transcriptomics has been widely applied in life science research, providing a powerful tool for studying the molecular mechanism of diseases[14]. Transcriptomics can comprehensively reveal all the transcript information in cells or tissues, so as to help us understand the changes in gene expression, the regulation of signaling pathways and the occurrence and development mechanism of related diseases. For muscle diseases, transcriptomic studies can reveal changes in gene expression related to muscle development, metabolism and signal transduction, providing important clues for disease prevention and treatment[15-17].

Transcriptomic studies of muscle diseases have made remarkable progress in the last few years. However, due to the difficulty in acquiring muscle samples and the heterogeneity of human muscle tissue, it is limited to study muscle transcriptomics directly in human body. Therefore, finding suitable animal models to replace human studies has become an effective strategy[18-20]. As a commonly used experimental model organism, zebrafish Danio rerio has advantages such as strong fecundity, short reproductive cycle, and clear genetic background, making it an ideal model for studying muscle diseases[21].

In this study, we will use zebrafish as a model organism to construct a zebrafish model of Laminin- Alpha 2 (LAMA2) muscular dystrophy, and use transcriptomic methods to conduct an in-depth study of its muscle transcriptome. Zebrafish have a highly conserved genome, and their muscle development and physiological processes are very similar to those of humans, so they are highly comparable and reliable.

In order to construct the LAMA2 zebrafish model, we will use Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing technology to knockout or mutate the LAMA2 gene of zebrafish. CRISPR/ Cas9 technology is an important breakthrough in genome editing in recent years because it is efficient, accurate and economical, and can realize the targeted editing of specific genes. Through the establishment of this model, we will simulate muscular dystrophy disease and deeply study the changes of transcriptome in its muscle tissue, so as to reveal the pathogenesis of muscular dystrophy[22-25].

In transcriptomics studies, we will apply highthroughput transcriptome sequencing techniques, such as RNA Sequencing (RNA-Seq), to sequence the collected muscle tissue samples of zebrafish. The RNA-Seq technique enables the quantification and identification of transcripts on a genome-wide scale, providing a comprehensive understanding of transcriptome changes. By analyzing the sequencing data, we will identify and compare the transcriptome expression differences in muscle tissue between LAMA2 model fish and wild-type zebrafish, and identify differentially expressed genes related to muscle development, metabolism and signal transduction. Through the construction of LAMA2 zebrafish model and transcriptomic studies, we can deeply understand the development process of muscle diseases, discover new therapeutic targets and formulate corresponding therapeutic strategies, and make contributions to improving the quality of life of patients. We will further promote the application of bioinformatics technology in the biomedical field, and open up new ways for the prevention and treatment of muscle diseases.

Materials and Methods

RNA-Seq data source:

Muscle tissue sample collection: Muscle tissue samples were respectively collected from LAMA2 model fish and wild-type zebrafish. During the collection process, we tried to select muscle tissues with the same anatomical position and similar size, and through rapid freezing and preservation, to maintain the integrity of RNA.

RNA extraction and library preparation: Total RNA was extracted from the collected muscle tissue, and the RNA was purified and quality controlled using appropriate reagents and methods. We then used the RNA-Seq library preparation kit to construct the RNA-Seq library according to the manufacturer's recommended procedures, including RNA reverse transcription, fragmentation, and sequencing splicing.

RNA-Seq: Constructed RNA-Seq libraries will be sent to high-throughput sequencing platforms, such as Illumina HiSeq or NextSeq, for sequencing. Through high throughput sequencing technology, we will obtain sequencing readings of all transcripts in muscle tissue.

Experimental methods:

Construction of the LAMA2 zebrafish model: Using CRISPR/Cas9 gene editing technology, we designed appropriate primers and probes to trigger LAMA2 gene knockout or mutation in zebrafish fertilized eggs by introducing Cas9 protein and the single guide RNA (sgRNA) of the LAMA2 gene. Through screening and identification, the zebrafish individuals with LAMA2 gene mutation were obtained, and the LAMA2 muscular dystrophy zebrafish model was established.

Collecting muscle tissue samples of zebrafish: Muscle tissue samples were collected from LAMA2 model fish and wild-type zebrafish. Using appropriate techniques and tools, we remove the muscle tissue from the fish body and maintain the integrity and purity of the sample as much as possible. The collected samples will be used for subsequent transcriptomic analysis.

Differential expression analysis: Differential expression analysis of two conditions/groups (two biological replicates per condition) was performed using the DESeq2 R package (1.16.1). DESeq2 provide statistical routines for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. The resulting p values were adjusted using the Benjamini and Hochberg’s approach for controlling the false discovery rate. Genes with an adjusted p<0.05 found by DESeq2 were assigned a differentially expressed gene[26].

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis: GO enrichment analysis of differentially expressed genes was implemented by the top GO R package, in which gene length bias was corrected[27]. GO terms with corrected p<0.05 were considered significantly enriched by differential expressed genes.

Protein-Protein Interaction (PPI) analysis of differentially expressed genes: PPI analysis of differentially expressed genes was based on the Search Tool for the Retrieval of Interacting Genes (STRING) database, which known and predicted PPIs[28].

Results and Discussion

The establishment of the LAMA2 zebrafish model was successfully verified through molecular biology and genetic analysis. The experiment was divided into two groups; LAMA2 knockdown group and Cas9 control group. In LAMA2 knockdown group, muscle phenotypes of juvenile fish in normal swimming state and unable to swim state were selected for identification. Juvenile fish in Cas9 control group had no abnormal state and could swim normally, so individuals in Cas9 group were randomly selected as muscle phenotype control. Detailed phenotypic analysis of the model was performed, including assessment of muscle structure and function. The differences between the model and wild-type zebrafish were compared. The muscle phenotype identification results of zebrafish model 4 showed that the muscle phenotypes of normal and unable swimming individuals in LAMA2 knockout group were abnormal, while the muscle phenotypes of Cas9 control group were not. Muscular dystrophy was identified in LAMA2 model fish (fig. 1A and fig. 1B).

pharmaceutical-sciences-group

Fig. 1: (A): In Cas9 control group, juvenile fish swam normally without any abnormality and (B): In LAMA2 group, juvenile fish could not swim normally

Using the transcriptome sequencing technology to LAMA2 model fish and wild type zebrafish muscle tissue samples were analyzed, and won the 3054 DEGs, 1201 DEGs of up and 1853 DEGs of down (fig. 2). In hierarchal clustering analysis, the regions with different colors represent different clustering information and have similar gene expression patterns in the same group. They may have similar functions or participate in the same biological processes (fig. 3).

pharmaceutical-sciences-logs

Fig. 2: Volcanic map analysis of muscle tissue samples of LAMA2 model fish and wild-type zebrafish Note: (Image): No difference; (Image): Down-regulated genes and (Image): Up-regulated genes

pharmaceutical-sciences-tissue

Fig. 3: Heat map analysis of muscle tissue samples of LAMA2 model fish and wild-type zebrafish Note: (Image): Cas9 and (Image): LAMA2, blue represents down-regulated genes and red represent up-regulated genes

Functional enrichment and pathway analysis of differentially expressed genes were performed to further reveal the biological processes and regulatory pathways associated with muscular dystrophy in LAMA2 model fish. Identify key signaling pathways and molecular mechanisms that may be involved in the pathogenesis of muscular dystrophy (fig. 4). Through differential expression gene analysis, MYH1, MYH2, MYH3, MYH4, ACTA1, TNNT1, TNNT3, TNNI1, MYOD1, MYOG, MYF5, MYF6, PAX3, PAX7, MRF4, MYH1, MYH2, MYH3, MYH4, Acta1, TNnt1, TNnt3, TNNI1, myod1, myog, myF5, myF6,IGF1, IGF2, FGF2, FGF6, related to metabolism are AMPK, PGC1α, PPARγ, SIRT1, GLUT1, GLUT4, G6PD, HK2, FASN, ACACA, PPARG, LPL, MAFbx, MuRF1, AKT1, mTOR and signal transduction related genes are WNT1, WNT3, CTNNB1, AXIN2, TGFB1, TGFBR1, SMAD2, SMAD3, MAPK1, MAPK3, MAPK8, MAPK14,PIK3CA, AKT1, FOXO1 and GSK3B (Table 1).

pharmaceutical-sciences-muscle

Fig. 4: GO/KEGG analysis of muscle tissue samples of LAMA2 model fish and wild-type zebrafish

Category Gene
Genes related to muscle development MYH1, MYH2, MYH3, MYH4, ACTA1, TNNT1, TNNT3, TNNI1, MYOD1, MYOG, MYF5, MYF6, PAX3, PAX7, MRF4, IGF1, IGF2, FGF2 and FGF6
Metabolism-related genes AMPK, PGC1α, PPARγ, SIRT1, GLUT1, GLUT4, G6PD, HK2, FASN, ACACA, PPARG, LPL, MAFbx, MuRF1, AKT1 and mTOR
Genes involved in signal transduction WNT1, WNT3, CTNNB1, AXIN2, TGFB1, TGFBR1, SMAD2, SMAD3, MAPK1, MAPK3, MAPK8, MAPK14, PIK3CA, AKT1, FOXO1 and GSK3B

Table 1: Functional Enrichment and Pathway Analysis Through Differential Expression Genes

The interaction relationship in the STRING protein interaction database (http://string-db.org/) was applied. For the species contained in the database, the interaction relationship of target gene set (such as differential gene list) was directly extracted from the database to construct the network (fig. 5). The filter was conducted according to the expression levels of all genes, and the Standard Deviation (SD) values of all samples were sorted from largest to smallest. The top 25% genes were selected for co-expression analysis (fig. 6).

pharmaceutical-sciences-network

Fig. 5: Schematic diagram of protein interaction network

pharmaceutical-sciences-core

Fig. 6: Schematic diagram of core presentation network Note: Each rectangular node in the figure represents a gene, and the associations between nodes represent coexpression relationships in which the correlation calculated by coexpression is higher than the threshold. Di fferent colors represent different network modules analyzed by WGCNA

The aim of this study was to construct a zebrafish model of LAMA2 muscular dystrophy and to reveal its underlying molecular mechanism through muscular transcriptomics studies[29-32]. Through the analysis and discussion of the experimental results, we can deeply understand the characteristics of muscle development and function abnormalities in the LAMA2 model, further elucidate the pathophysiological process of the disease, and provide new clues for the development of therapeutic strategies[33-36].

In this study, we successfully constructed a LAMA2 zebrafish model to simulate the clinical phenotype of muscular dystrophy. Using single-cell sequencing and RNA-Seq techniques, we obtained and comprehensively analyzed the transcriptome data of muscle tissue in the LAMA2 model. In single-cell sequencing data, we identified and classified multiple muscle cell subtypes and identified changes in gene expression associated with muscle development and function. RNASeq data analysis revealed a large number of differentially expressed genes and pathways in the LAMA2 model, which are involved in biological processes such as muscle structure, metabolic regulation, and cellular stress. Through differential expression gene analysis, MYH1, MYH2, MYH3, MYH4, ACTA1, TNNT1, TNNT3, TNNI1, MYOD1, MYOG, MYF5, MYF6, PAX3, PAX7, MRF4, MyH1, MyH2, MyH3, MyH4, Acta1, TNnt1, TNnt3, TNNI1, myod1, myog, myF5, myF6, IGF1, IGF2, FGF2, FGF6, related to metabolism are AMPK, PGC1α, PPARγ, SIRT1, GLUT1, GLUT4, G6PD, HK2, FASN, ACACA, PPARG, LPL, MAFbx, MuRF1, AKT1, mTOR and signal transduction related genes are WNT1, WNT3, CTNNB1, AXIN2, TGFB1, TGFBR1, SMAD2, SMAD3, MAPK1, MAPK3, MAPK8, MAPK14,PIK3CA, AKT1, FOXO1 and GSK3B[37-42].

Further discussing the results of the zebrafish model of LAMA2 muscular dystrophy, we can compare the study results with existing clinical and molecular genetic data. Although zebrafish models cannot fully simulate human diseases, we can find some common pathological features and molecular mechanisms through comparison. For example, we observed muscle degeneration and atrophy in the LAMA2 model, which is consistent with clinical manifestations of muscular dystrophy in humans[29-32]. In addition, by analyzing differentially expressed genes and pathways, we can identify some potential therapeutic targets and drug candidates. For example, targeting the differentially expressed gene family of muscle structural proteins[33-35]. We can explore the efficacy of related drugs to promote muscle regeneration and functional recovery.

Although we achieved a series of meaningful results in this study, there are some research limitations and future directions for improvement. First of all, this study focuses on transcriptomic analysis, which can be combined with proteomics and metabolomics to obtain more comprehensive molecular information in the future. Second, although we have verified the changes of some differentially expressed genes, the function and regulatory mechanism of these genes still need to be further studied. In addition, further functional experiments and gene knockout can be used to confirm the importance of differentially expressed genes in muscle development and function.In conclusion, this study provides an in-depth understanding of LAMA2 muscular dystrophy zebrafish through the construction of the model and muscle transcriptomics studies. We reveal differentially expressed genes and pathways associated with abnormal muscle development and function in the LAMA2 model, providing important clues for further study of the pathophysiological mechanisms and potential therapeutic targets of this disease. This study provides new ideas and directions for the treatment and clinical management of muscular dystrophy and is expected to make important contributions to the health and well-being of patients.

Successful construction of the LAMA2 zebrafish model; successfully simulated the clinical phenotype of LAMA2 muscular dystrophy, including muscle degeneration and atrophy. This provides a reliable model system for studying the pathophysiological process of the disease and seeking treatment strategies.

Muscle transcriptomic studies reveal differentially expressed genes and pathways. Using single-cell sequencing and RNA-Seq techniques, we identified changes in gene expression in muscle cells in the LAMA2 model. We identified multiple muscle cell subtypes and identified differentially expressed genes and pathways related to muscle development and function.

The molecular mechanisms underlying abnormal muscle development and function are revealed. Our results indicate abnormalities in the LAMA2 model involved in biological processes such as cytoskeletal recombination, mitochondrial function impairment, and inflammatory responses. These changes may be closely related to pathological features such as muscle degeneration and atrophy.

Providing clues to therapeutic targets and new drug development through analysis of differentially expressed genes and pathways, we identified a number of potential therapeutic targets and drug candidates. In particular, the study of differentially expressed gene families of muscle structural proteins is expected to provide new directions for therapeutic strategies for muscle regeneration and functional recovery.

Funding:

This article was supported by the scientific innovation start-up foundation of Fujian Maternity and Child Health Hospital (YCXM20-23).

Conflicts of interests:

The authors declare no conflict of interests.

References