- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, China
Skeletal muscle growth and development is a complicated process that is regulated at multiple steps and by numerous myogenesis genes. RNA editing represents one of the events at the post-transcriptional level, which contributes to the diversity of transcriptome and proteome by altering the nucleotides of RNAs. However, RNA editing events in the skeletal muscle of yaks are still not well defined. This study conducted whole-genome RNA-editing identification in skeletal muscle of yaks at embryonic stage (ES) and adult stage (AS). We found a total of 11,168 unique RNA editing sites, most of which were detected in the intergenic region. After annotation, we totally identified 2,718 editing sites within coding regions, among which 858 were missense changes. Moreover, totally 322 editing sites in the 3′ untranslated regions (UTR) were also predicted to alter the set of miRNA target sites, indicating that RNA editing may be involved in translational repression or mRNA degradation. We found 838 RNA editing sites (involving 244 common genes) that are edited differentially in ES as compared to AS. According to the KEGG enrichment analysis, these differentially edited genes were mainly involved in pathways highly related to skeletal muscle development and myogenesis, including MAPK, AMPK, Wnt, and PI3K-Akt signaling pathways. Altogether, our work presents the first characterization of RNA editing sites within yak skeletal muscles on a genome-wide scale and enhances our understanding of the mechanism of skeletal muscle development and myogenesis.
Introduction
Yak (Bos grunniens) is a vital species of livestock prevalent across the Qinghai-Tibet plateau and the adjacent alpine regions (1). It plays a crucial role in promoting the local economy by providing products such as meat, milk, hair, transport, and fuel for the residents. Yak meat is highly nutritious, being rich in protein and amino acids, and having low-fat content (2). Nonetheless, compared to cattle, the growth rate of yak is slow due to the lack of an efficient yak breeding program to improve the growth traits. Therefore, improving the growth rate is one of the significant breeding objectives in the yak industry.
Skeletal muscle accounts for ~40% of the bodyweight and is the main meat-producing tissue. In livestock, the development of skeletal muscle has a significant influence on the growth rate and meat yield (3). Vertebrate skeletal muscle has an essential role in locomotion and metabolism, and it originates from paraxial mesoderm in the embryonic stages (4). In prenatal period, increasing fiber numbers is the main characteristic of skeletal muscle growth. After birth, muscle development mainly depends on muscle fiber hypertrophy. Skeletal muscle development (myogenesis) is a complicated but orderly process involving commitment of multipotent precursor cells into myoblasts, myoblast proliferation, differentiation, fusion to myofibers, and myotubes with multiple nuclei, as well as their eventual accommodation into slow-and fast-twitch muscle fibers (5). This process is precisely orchestrated via several well-known transcription factors (TFs), such as the myocyte enhancer factor 2 (MEF2) family members (6), paired box protein 3/7 (Pax3/7) (7), and members of the myogenic regulatory factors (MRFs) family (8). Although many protein-encoding genes and non-coding RNAs (ncRNA) have been shown to make important contributions to the intricate process of muscle development, the distinct mechanisms have not yet been elucidated completely.
RNA editing is one of the most important post-transcriptional events altering the nucleotide composition of a transcript through the insertion, deletion, or substitution of nucleotides. It leads to differences bcetween a genomic DNA sequence and its corresponding mRNA sequence. RNA editing may occur within coding genes or non-coding regions, thereby resulting in non-synonymous substitutions, regulating alternative splicing (AS), and stability of RNAs (9). Notably, the A-to-I RNA editing is a frequently occurring RNA editing type within many animals, which is mediated by the adenosine deaminases acting on RNA (ADAR family) (10). In double-stranded RNAs (dsRNAs), ADAR family members can deaminize adenosine (A) to inosine. The latter is read by the cell machinery as guanosine (G). Additionally, cytidine-to-uridine (C-to-U) editing is also a common RNA editing type observed in animals, which is triggered by enzymes of the apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) family (11). Recently, many bioinformatic tools and next-generation sequencing (NGS) techniques have been developed for detecting RNA editing events, which has enabled us to discover genome-wide RNA editing events and their molecular function. Presently, many studies for identification of RNA editing events have been performed in human (12), cattle (13), pigs (14), sheep (15), and chicken (16). Many studies have shown that dysregulation of RNA editing is associated with various human diseases, such as the presence of tumors (17), epilepsy (18), Alzheimer's disease (19), and sporadic amyotrophic lateral sclerosis (ALS) (20). The importance of RNA editing in myogenesis was demonstrated in a recent study by Yang et al. (21). In this study, the authors outlined the landscape of the RNA editome in porcine skeletal muscle across 27 developmental stages and identified a series of RNA editing sites related to myogenesis, indicating the crucial role of RNA editing in the process of skeletal muscle development. However, there is limited information available on the role of RNA editing in skeletal muscle development in yak. In this work, we present a genome-wide landscape of the RNA editome in yak skeletal muscle at embryonic stage (ES) and adult stage (AS) and evaluate the potential role of RNA editing events in muscle development. Our study largely extends the list of RNA editing sites in yak and provides valuable insights for better understanding of the regulatory mechanisms of muscle development in yak.
Materials and Methods
Data Resource
In order to investigate the RNA editing sites in muscle tissues of yak, six specific strand RNA-seq data (PRJNA550017) from our previous study (22) were downloaded from NCBI SRA (https://www.ncbi.nlm.nih.gov/sra/). This dataset was developed from longissimus dorsi muscle tissues of 3 embryonic (90 days old) and 3 adult (3 years old) yaks. The Illumina HiSeq 2500 platform was employed to sequence cDNA libraries of these samples (2 × 150-bp paired-end read length).
Whole Genome Sequencing (WGS)
Blood samples were gathered from yaks having the same genetic background as the study samples used for RNA-seq. All experiments and animal care procedures were carried out according to the guidelines of the Animal Administration and Ethics Committee of Lanzhou Institute of Husbandry and Pharmaceutical Sciences of CAAS (Permit No. SYXK-2014–0002). DNA extraction kit (Tiangen Biotech, Beijing, China) was utilized to extract genomic DNA from blood in line with the specific instructions. DNA libraries were constructed with the use of the MGIEasy Universal DNA Library Prep Kit (MGI, Shenzhen, China) following the BGI's standard preparation protocol. Paired-end sequencing (2 × 150-bp) was carried out on an MGISEQ2000 (MGI, Shenzhen, China) at Frasergen Bioinformatics Co., Ltd (Wuhan, China).
Quality Control and Reads Mapping
Reads including ploy-N, adapter or low-quality sequences were filtered out to obtain clean reads. The Burrows-Wheeler Aligner (BWA mem; version v0.7.17) (23) was used to align clean reads of WGS data to domestic yak genome sequence (LU_ Bosgru_v3.0). Samtools (Version 1.9) (24) was employed for sorting mapped reads, and duplicated reads were eliminated with Picard tools v. 2.13.2 (https://broadinstitute.github.io/picard/) using the function MarkDuplicate. We then mapped the clean reads from RNA-seq data to the genome sequence of domestic yak (LU_Bosgru_v3.0) using the HISAT2 (2.1.0) (25). The resulting SAM files were converted into BAM files using Samtools (Version 1.9) (24).
RNA Editing Detection
In order to detect the RNA editing sites, the REDItoolDnaRna.py in REDItools v1.0.4 (26) was utilized. Parameters recommended by the researcher were chosen after mild modifications (minimal quality score (-q) 25, minimal mapping quality score (-m) 25, minimal editing frequency (-n) 0.05, minimal homopolymeric length (-O) 5, minimal read coverage (-c) 10, and minimal read number that supports variation (-v) 3) (27). Finally, only the edited sites were retained for further study, which were present in at least two samples. The SnpEff (v4.3t) (28) was employed to annotate RNA editing sites according to Ensembl-based gene annotation (release 90).
Validation of RNA Editing Sites Through Sanger Sequencing
To validate the reliability of RNA editing sites identified in this study, three editing sites were randomly selected for PCR validation. Six independent samples used for PCR validation were collected at the same developmental stage as those used for RNA-Seq. Total RNA was extracted using the animal tissue RNA isolation kit (ZDGSY, Beijing, China) following the manufacturer's protocol. Genomic DNA was isolated using the animal tissue Genomic DNA Kit (ZDGSY, Beijing, China). Total RNA was used for reverse transcription by PrimeScript™ RT reagent kit with gDNA Eraser (TaKaRa, Dalian, China). Equal amounts of cDNA were mixed. The primer details are presented in Supplementary Table S1. The 20 μl PCR reactions contained 1 μl of mixed cDNA, 1 μl of each sense and anti-sense, 10 μl of GoTaq® Green Master Mix (Promega, Madison, WI, USA) and 7 μl of ddH2O. The PCR program was set as follows: 95°C, 2 min, (95°C, 30s; 58°C, 30 s; 72°C, 1 min) 30 cycles and 72°C, 5 min. The PCR products were then sequenced by Sanger sequencing.
Analyzing the Effects of RNA Editing on MiRNA Regulation
We used RNAhybrid (-b 1 -c -f 2,8 -m 100,000 -u 1 -v 1 -e−10) (29) and miRanda (sc 140 -en−10 -scale 4 -strict) (30) to detect binding targets of the miRNAs on the edited and reference sequences for editing sites in the 3'untranslated region (UTR). According to the predicted miRNA-mRNA interactions supported by both tools, the gain of miRNA binding targets was defined as the interactions existing in the edited sequences but not in the reference sequences. On the contrary, the loss of miRNA binding targets was defined as the interactions missing in the edited sequences but not in the reference sequences.
Differential RNA-Editing Analysis
For an RNA editing site in a given sample, the RNA editing level was calculated as the ratio of the reads supporting the edited base to the total number of reads detected on this site. To identify RNA-edited sites related to muscle development, we performed Tukey's Honest Significant Difference approach to screen significantly differently edited sites between ES and AS. Sites with P-value < 0.05 were labeled as significant.
Annotation and Enrichment Analysis
The genes with RNA editing sites located in exonic regions were extracted for functional enrichment analysis. To explore the functions of these putative RNA editing sites, we carried out gene ontology (GO) analysis based on the G: Profile web approach (31). KOBAS 3.0 software was employed for Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis by performing a hypergeometric test (32). The cut-off value of significant GO terms and KEGG pathway was P-value < 0.05.
Conservation Analysis of RNA Editing Sites
For detecting A-to-I editing events with a high conservation degree, editing sites discovered in the present work were compared with those identified in humans based on the REDIportal database (http://srv00.recas.ba.infn.it/atlas/) (33). In total, 4,627,557 editing sites in humans were adopted to conduct conservation analyses. We extracted flanking regions of 50 bp at the editing sites in domestic yak genome sequence (LU_Bosgru_v3.0) and blasted them against 50 bp flanking regions of the identified A-to I editing events in humans using the Nucleotide Basic Local Alignment Search Tool (Blastn). Hits with Expect (E) values ≤ 1e−5 and identity ≥ 85% indicated editing events with high conservation. The E-value served as a measure of the number of hits one can “expect” to see by chance when searching a database.
Results
RNA Editing Detection and Validation
To determine RNA editing sites at the genomic level in the longissimus dorsi muscle, we collected six strand-specific RNA sequencing samples from both the ES and AS. After trimming the adaptor sequence and low-quality reads, we obtained 282,913,382 total clean reads for ES and 363,479,946 total clean reads for AS. The total mapped ratio between the reads and the reference genome of all samples ranged from 88.19 to 92.58% (Table 1). We identified a total of 31,244 raw editing sites using strict filtering criteria (Supplementary Table S2). After potential SNPs filtering, 11,168 high-confidence RNA editing events were kept (Supplementary Table S3). The ES1 individual showed the fewest, whereas the AS1 individual showed the highest number of editing events. Additionally, these RNA editing events showed a non-uniform location within the yak chromosomes (Figure 1). The highest number of editing sites (670) was observed on chromosome 8, while the least number of editing sites (126) was found on chromosome 29. As our expectation, the number of RNA editing sites was different between ES and AS. Our analysis indicates that 6,829 RNA editing sites were shared between the two groups, 1,075 editing sites were specific to ES, and 3,264 editing sites were specific to AS (Figure 2A). A total of 12 different types of RNA editing were identified, each with a proportion >1%. Over 40% of these sites correspond to A-to-G and C-to-T type, consistent with A-to-I and C-to-U editing (Figure 2B). To validate the predicted RNA editing sites, three sites were randomly selected for PCR and Sanger sequencing. The sites were considered to be successfully verified if the cDNA sequence was heterozygous while the corresponding DNA sequencing was homozygous. After aligning the cDNA and DNA sequences, selected editing sites were consistent with the prediction (Supplementary Figure S1).
Figure 2. Characteristics of the editome in the skeletal muscle. (A) Venn diagram showing the number of shared editing sites between AS and ES; (B) Proportion of 12 types of RNA editing; (C) Neighbor sequence preferences of A-to-G RNA editing; (D) Distribution of the identified RNA editing sites throughout the genome.
Sequence Preference and Annotation of RNA Editing Sites in Skeletal Muscle
As shown in Figure 2C, G is underrepresented at −1 position and overrepresented at +1 position of the A-to-G editing sites discovered in our study. The identified sequence preference for A-to-G editing sites corresponds to the findings of previous reports on mammalian ADAR (34, 35). To characterize the distribution of RNA editing sites, the SnpEff tool was applied to identify each RNA editing site that corresponded to the annotated gene. Totally RNA editing sites overlapped with 2,584 annotated genes. The RNA editing sites identified in this study were annotated to 7 types of genomic locations: intergenic regions, introns, coding sequences (CDS), 3′ UTRs, 5′ UTRs, ncRNA, and splice sites (Supplementary Table S4). Most RNA editing sites were localized in the intergenic region (38.15%), followed by introns (28.38%) (Figure 2D). Furthermore, a total of 858 editing events residing in CDS were identified as missense editing events, which could influence the functions of protein.
Cross-Species Analysis Between Yak and Human
To investigate the conservation of RNA editomes between human and yak, upstream and downstream sequence of editing sites was aligned using nucleotide BLAST tools. With a rigorous filtering threshold (identity >85%, E-value <1e−5), totally eight conserved RNA editing sites were identified (Supplementary Table S5). Of these, editing site in the COPI coat complex subunit alpha (COPA) gene resulted in missense substitution.
Impact of RNA Editing on MiRNA Regulation
RNA editing events residing in miRNA binding sites could change the miRNA-mRNA interactions and modulate translational suppression or mRNA degradation. For the RNA editing sites residing in 3' UTRs, TargetScan and miRanda were adopted for detecting the miRNA binding targets for each reference sequence. By predicting miRNA targets, 232 editing sites were presumed to generate 223 novel miRNA binding sites and abolish 172 primitive ones (Figure 3; Supplementary Table S6). For a better understanding of the functional roles of these 190 target genes, GO and KEGG pathway enrichment analyses were performed. As shown in Supplementary Figure S2, these target genes were categorized into two categories: biological process (BP) and cellular component (CC). The five most significantly enriched terms in the BP category were macromolecule localization, cellular macromolecule localization, cellular localization, localization, and protein localization. According to KEGG analyses, the impacted target genes were mostly enriched in pathways associated with muscle development, such as mTOR, Notch, and insulin signaling pathways (Supplementary Figure S3).
Differential RNA-Editing Analysis in Skeletal Muscle at Different Periods
The average editing level for each sample ranged from 0.23 to 0.32 (Figure 4A). Based on hierarchical clustering, it could be found that the discrepancies in the editing levels within groups were less than those between groups (Figure 4B). These results suggest that RNA editing sites at the genome-wide level can be adopted for characterizing the developmental stage of skeletal muscle. To identify the muscle development-associated RNA editing events, this work utilized Tukey's Honest Significant Difference approach to search for the differential editing sites between ES and AS. We discovered a total of 838 differential editing events (Supplementary Table S7). GO analysis indicated that these genes with differential editing levels were involved in 120 terms (Supplementary Table S8). In the BP category, muscle cell development, muscle cell differentiation, and striated muscle cell differentiation were the most abundant terms (Figure 5A). KEGG pathway analyses demonstrated that these differentially edited genes were significantly enriched in muscle development-related pathways, such as AMPK signaling pathway, focal adhesion, MAPK signaling pathway, insulin signaling pathways, ECM-receptor interaction, Wnt signaling pathway and PI3K-Akt signaling pathway (Figure 5B; Supplementary Table S9).
Figure 4. Characteristics of RNA editing levels within and between groups. (A) Distribution of editing levels across samples; (B) Hierarchical clustering of RNA editing levels across samples.
Figure 5. GO and KEGG enrichment of genes with differential editing levels. (A) GO enrichment terms; (B) KEGG pathway enrichment terms. Rich Factor is the ratio of differentially edited gene numbers annotated in this pathway term to all gene numbers annotated in this pathway term.
Discussion
RNA editing represents one of the important modifications at the post-transcriptional level, and it alters RNA nucleotide sequences, thereby impacting the mRNA structure and proteomic diversity. It was first identified in trypanosomes (36). Thanks to the continuous advancements in NGS technologies, many RNA editing events have been discovered in animals (37) and plants (38). To date, however, RNA editing in yak has been far less studied. To explore how RNA editing affects the muscle development processes in yak, we identified RNA editing sites at ES and AS at the whole-sequence level.
The present study discovered 11,168 editing events in yak longissimus dorsi muscle. It is worth noting that there are lesser editing sites at embryonic stage as compared with the adult stage, which suggests the decreased activity of RNA editing at embryonic stage. Notably, the two canonical RNA editing events, viz., A-to-I and C-to-U, occupy up to 41.12% of those putative editing events, indicating that the accuracy of our results is high. Additionally, eight RNA editing sites were conserved between yak and human. Consistent with the results obtained for previous studies (13, 21), the low level of overlap between editing sites identified in this study and the reported human editing sites suggests that most RNA editing sites are poorly conserved through evolution.
Unlike other studies in which the intronic region includes the most RNA editing sites (39), our study indicated that the intergenic region (38.15%) had the most RNA editing sites, followed by the intronic region (28.38%). Two possible explanations may account for this discrepancy. The first is that the annotation of the yak reference genome was inaccurate and incomplete compared with the human reference genome. The second possibility is that our RNA-seq libraries which were prepared with RiboZero rRNA Removal Kit, contain a large number of intergenic lncRNAs. However, our results indicate a potentially relationship between RNA editing events that reside in non-coding regions and the regulation of gene expression.
The recoding RNA editing sites lead to non-synonymous replacements that increase proteomic diversity (40). It is reported that numerous recoding RNA editing sites are preserved and have functional and evolutionary importance (41, 42). Altogether 858 recoded RNA editing sites were identified within 650 genes, some of which were related to skeletal muscle development. For example, myosin heavy chain 3 (MYH3) belongs to the myosin heavy chain (MYH) family and is predominantly expressed in distinct muscle developmental stages. Missense mutations in the MYH3 gene have been reported to give rise to human muscle development disorders (43). Sad1 And UNC84 Domain Containing 1 (SUN1) belongs to the linker of nucleoskeleton and cytoskeleton (LINC) complex, which plays a critical role in myotube formation (44). A previous study revealed that loss of SUN1 leads to myofibers with a smaller diameter, and this slows down the adult skeletal muscle regeneration (45). Myosin 18B (MYO18B) is the new nontraditional myosin heavy chain that is mostly expressed in human cardiac and skeletal muscle (46). Disruption of MYO18B could suppress the proliferation and differentiation of C2C12 mouse myoblasts (47). Therefore, these recoding RNA editing sites may have a crucial role in the skeletal muscle development in the yak. Except for these recoding editing sites, some RNA editing sites within the 3′ UTR of mRNAs could regulate the expression of target genes by creating or destroying miRNA binding sites (40). In the present study, we discovered 253 RNA editing sites that might alter target gene expression by changing the power of miRNA binding. The KEGG enrichment analysis indicated that the impacted target genes were enriched in some critical pathways associated with muscle development, including hippo (48), mTOR, and insulin signaling pathways (49). These results suggest that RNA editing events within miRNA binding sites might regulate skeletal muscle development via translational repression or mRNA degradation.
In this study, a total of 838 sites were found to be differentially edited between the AS and ES groups. According to KEGG analyses, genes that showed different editing levels were mostly enriched in several important pathways associated with muscle development, e.g., MAPK, AMPK, Wnt and PI3K-Akt signaling pathways and ECM-receptor interaction. The MAPK signaling pathway was considered to be a crucial regulator of skeletal muscle growth and development (50). Myocyte enhancer factor 2C (MEF2C) and voltage-dependent, alpha-2/delta subunit 1(CACNA2D1) are involves in MAPK signaling pathway. MEF2C belongs to the MADS-box transcription enhancer factor 2 (MEF2) family, which is known to be involved in myogenesis. It was found that miR-204–5p can inhabit the myoblast differentiation through repression MEF2C gene (51). CACNA2D1 belongs to the alpha-2/delta subunit family, which is associated with voltage-gated calcium channels. Variations of the CACNA2D1 gene are significantly associated with bovine carcass traits (52). Laminin 2 (LAMA2) encodes an important extracellular matrix (ECM) protein which is mainly expressed in the basement membrane of skeletal muscle (53). Variants in LAMA2 were identified to be the causal mutations of congenital merosin-deficient muscular dystrophy (54). These observations suggest that RNA editing may influence skeletal muscle development and myogenesis.
Conclusion
In conclusion, the present study comprehensively analyzes the RNA editome in the skeletal muscle of yaks at ES and AS. We identified 11,168 high-confidence RNA editing sites. Of these, many RNA editing sites with different editing levels may potentially contribute to myogenesis and muscle development. However, their biological functions and regulatory mechanisms need to be further investigated in depth. Our research expands the list of RNA editing sites in yak and offers profounder insight into understanding the mechanism of yak muscle development.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
The animal study was reviewed and approved by Animal Administration and Ethics Committee of Lanzhou Institute of Husbandry and Pharmaceutical Sciences of CAAS (Permit No. SYXK-2014–0002).
Author Contributions
XW, PY, and CL conceived the research and drafted the manuscript with comments from all authors. XW, MC, and XM conducted the experiments and analyses. JP, LX, and XG participated in experiments and revised the manuscript. All authors approved the final version.
Funding
This study was supported by the National Natural Science Foundation of China (32102500), Agricultural Science and Technology Innovation Program (25-LZIHPS-01), Gansu Province Science Foundation for Youths (20JR5RA576), Foundation for Innovation, Groups of Basic Research in Gansu Province (20JR5RA580), and Central Public-interest Scientific Institution Basal Research Fund (Y2022LM10).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.871814/full#supplementary-material
References
1. Wiener G, Han J, Long R. The Yak, 2nd ed. FAO Regional office for Asia and the Pacific: Bangkok, Thailand (2003).
2. Tian C. Han, L Yu, Q, Shi X, Wang W. Changes in tenderness and cathepsins activity during post mortem ageing of yak meat. Can J Anim Sci. (2013) 93:321–28. doi: 10.4141/cjas2012-102
3. Rehfeldt C, Te P M, Wimmers K, Brameld J M, Nissen P M, Berri C, et al. Advances in research on the prenatal development of skeletal muscle in animals in relation to the quality of muscle-based food. II–Genetic factors related to animal performance and advances in methodology. Animal. (2011) 5:718–30. doi: 10.1017/S1751731110002454
4. Buckingham M. Skeletal muscle formation in vertebrates. Curr Opin Genet Dev. (2001) 11:440–48. doi: 10.1016/S0959-437X(00)00215-X
5. Bentzinger C F, Wang Y X, Rudnicki M A. Building muscle: Molecular regulation of myogenesis. Cold Spring Harb Perspect Biol. (2012) 4:a8342. doi: 10.1101/cshperspect.a008342
6. Edmondson D G, Lyons G E, Martin J F, Olson E N. Mef2 gene expression marks the cardiac and skeletal muscle lineages during mouse embryogenesis. Development. (1994) 120:1251–63. doi: 10.1242/dev.120.5.1251
7. Relaix F, Montarras D, Zaffran S, Gayraud-Morel B, Rocancourt D, Tajbakhsh S, et al. Pax3 and Pax7 have distinct and overlapping functions in adult muscle progenitor cells. J Cell Biol. (2006) 172:91–102. doi: 10.1083/jcb.200508044
8. Hernández-Hernández J M, García-González E G, Brun C E, Rudnicki M A. The myogenic regulatory factors, determinants of muscle development, cell identity and regeneration. Semin Cell Dev Biol. (2017) 72:10–8. doi: 10.1016/j.semcdb.2017.11.010
9. Gott J M, Emeson R B. Functions and mechanisms of RNA editing. Annu Rev Genet. (2000) 34:499–531. doi: 10.1146/annurev.genet.34.1.499
10. Nishikura K. A-to-I editing of coding and non-coding RNAs by ADARs. Nat Rev Mol Cell Biol. (2016) 17:83–96. doi: 10.1038/nrm.2015.4
11. Schmidt J W, Wehde B L, Sakamoto K, Triplett A A, Anderson S M, Tsichlis P N, et al. Stat5 regulates the phosphatidylinositol 3-kinase/Akt1 pathway during mammary gland development and tumorigenesis. Mol Cell Biol. (2014) 34:1363–77. doi: 10.1128/MCB.01220-13
12. Tan M H, Li Q, Shanmugam R, Piskol R, Kohler J, Young A N, et al. Dynamic landscape and regulation of RNA editing in mammals. Nature. (2017) 550:249–54. doi: 10.1038/nature24041
13. Bakhtiarizadeh M R, Salehi A, Rivera R M. Genome-wide identification and analysis of A-to-I RNA editing events in bovine by transcriptome sequencing. PLoS ONE. (2018) 13:e193316. doi: 10.1371/journal.pone.0193316
14. Martínez-Montes A M, Fernández A, Pérez-Montarelo D, Alves E, Benítez R M, Nuñez Y, et al. Using RNA-Seq SNP data to reveal potential causal mutations related to pig production traits and RNA editing. Anim Genet. (2017) 48:151–65. doi: 10.1111/age.12507
15. Zhang Y, Han D, Dong X, Wang J, Chen J, Yao Y, et al. Genome-wide profiling of RNA editing sites in sheep. J Anim Sci Biotechnol. (2019) 10:31. doi: 10.1186/s40104-019-0331-z
16. Frésard L, Leroux S, Roux P F, Klopp C, Fabre S, Esquerré D, et al. Genome-wide characterization of RNA editing in chicken embryos reveals common features among vertebrates. PLoS ONE. (2015) 10:e126776. doi: 10.1371/journal.pone.0126776
17. Baysal B E, Sharma S, Hashemikhabir S, Janga S C, RNA. editing in pathogenesis of cancer. Cancer Res. (2017) 77:3733–39. doi: 10.1158/0008-5472.CAN-17-0520
18. Srivastava P K, Bagnati M, Delahaye-Duriez A, Ko J H, Rotival M, Langley S R, et al. Genome-wide analysis of differential RNA editing in epilepsy. Genome Res. (2017) 27:440–50. doi: 10.1101/gr.210740.116
19. Wu S, Yang M, Kim P, Zhou X. ADeditome provides the genomic landscape of A-to-I RNA editing in Alzheimer's disease. Brief Bioinform. (2021) 22:bbaa384. doi: 10.1093/bib/bbaa384
20. Hosaka T, Tsuji H, Kwak S, RNA. editing: A new therapeutic target in amyotrophic lateral sclerosis and other neurological diseases. Int J Mol Sci. (2021) 22:10958. doi: 10.3390/ijms222010958
21. Yang Y, Zhu M, Fan X, Yao Y, Yan J, Tang Y, et al. Developmental atlas of the RNA editome in Sus scrofa skeletal muscle. DNA Res. (2019) 26:261–72. doi: 10.1093/dnares/dsz006
22. Ma X, Jia C, Chu M, Fu D, Lei Q, Ding X, et al. Transcriptome and DNA methylation analyses of the molecular mechanisms underlying with longissimus dorsi muscles at different stages of development in the polled yak. Genes (Basel). (2019) 10:970. doi: 10.3390/genes10120970
23. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. (2009) 25:1754–60. doi: 10.1093/bioinformatics/btp324
24. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. (2009) 25:2078–79. doi: 10.1093/bioinformatics/btp352
25. Kim D, Langmead B, Salzberg S L, HISAT. A fast spliced aligner with low memory requirements. Nat Methods. (2015) 12:357–60. doi: 10.1038/nmeth.3317
26. Picardi E, Pesole G. REDItools: High-throughput RNA editing detection made easy. Bioinformatics. (2013) 29:1813–14. doi: 10.1093/bioinformatics/btt287
27. Picardi E, D'Erchia A M, Montalvo A, Pesole G. Using REDItools to detect RNA editing events in NGS datasets. Curr Protoc Bioinform. (2015) 49:12. doi: 10.1002/0471250953.bi1212s49
28. Cingolani P, Platts A, Wang L L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; Iso-2; Iso-3. Fly (Austin). (2012) 6:80–92. doi: 10.4161/fly.19695
29. McGeary S E, Lin K S, Shi C Y, Pham T M, Bisaria N, Kelley G M, et al. The biochemical basis of microRNA targeting efficacy. Science. (2019) 366:eaav1741. doi: 10.1126/science.aav1741
30. John B, Enright A J, Aravin A, Tuschl T, Sander C, Marks D S. Human MicroRNA targets. PLoS Biol. (2004) 2:e363. doi: 10.1371/journal.pbio.0020363
31. Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, et al. G:Profiler: A web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. (2019) 47:W191–98. doi: 10.1093/nar/gkz369
32. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S., et al. KOBAS 20: A web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. (2011) 39:W316–22. doi: 10.1093/nar/gkr483
33. Picardi E, D'Erchia AM, Lo GC, Pesole G. REDIportal: A comprehensive database of A-to-I RNA editing events in humans. Nucleic Acids Res. (2017) 45:D750–57. doi: 10.1093/nar/gkw767
34. Chen JY, Peng Z, Zhang R, Yang XZ, Tan BC, Fang H, et al. RNA editome in rhesus macaque shaped by purifying selection. PLoS Genet. (2014) 10:e1004274. doi: 10.1371/journal.pgen.1004274
35. Carmi S, Borukhov I, Levanon EY. Identification of widespread ultra-edited human RNAs. PLoS Genet. (2011) 7:e1002317. doi: 10.1371/journal.pgen.1002317
36. Benne R, Van den Burg J, Brakenhoff JP, Sloof P, Van Boom JH, Tromp MC. Major transcript of the frameshifted coxII gene from trypanosome mitochondria contains four nucleotides that are not encoded in the DNA. Cell. (1986) 46:819–26. doi: 10.1016/0092-8674(86)90063-2
37. Danecek P, Nellåker C, McIntyre RE, Buendia-Buendia JE, Bumpstead S, Ponting CP, et al. High levels of RNA-editing site conservation amongst 15 laboratory mouse strains. Genome Biol. (2012) 13:26. doi: 10.1186/gb-2012-13-4-r26
38. Lo G C, Hernández I, Ceci L R, Pesole G, Picardi E, RNA. editing in plants: A comprehensive survey of bioinformatics tools and databases. Plant Physiol Biochem. (2019) 137:53–61. doi: 10.1016/j.plaphy.2019.02.001
39. Hwang T, Park CK, Leung AK, Gao Y, Hyde TM, Kleinman JE, et al. Dynamic regulation of RNA editing in human brain development and disease. Nat Neurosci. (2016) 19:1093–99. doi: 10.1038/nn.4337
40. Eisenberg E, Levanon EY. A-to-I RNA editing-immune protector and transcriptome diversifier. Nat Rev Genet. (2018) 19:473–90. doi: 10.1038/s41576-018-0006-1
41. Burns CM, Chu H, Rueter SM, Hutchinson LK, Canton H, Sanders-Bush E, et al. Regulation of serotonin-2C receptor G-protein coupling by RNA editing. Nature. (1997) 387:303–08. doi: 10.1038/387303a0
42. Köhler M, Burnashev N, Sakmann B, Seeburg PH. Determinants of Ca2+ permeability in both TM1 and TM2 of high affinity kainate receptor channels: Diversity by RNA editing. Neuron. (1993) 10:491–500. doi: 10.1016/0896-6273(93)90336-P
43. Toydemir RM, Rutherford A, Whitby FG, Jorde LB, Carey JC, Bamshad MJ. Mutations in embryonic myosin heavy chain (MYH3) cause Freeman-Sheldon syndrome and Sheldon-Hall syndrome. Nat Genet. (2006) 38:561–65. doi: 10.1038/ng1775
44. Mattioli E, Columbaro M, Capanni C, Maraldi NM, Cenni V, Scotlandi K, et al. Prelamin A-mediated recruitment of SUN1 to the nuclear envelope directs nuclear positioning in human muscle. Cell Death Differ. (2011) 18:1305–15. doi: 10.1038/cdd.2010.183
45. Loo TH, Ye X, Chai RJ, Ito M, Bonne G, Ferguson-Smith AC, et al. The mammalian LINC complex component SUN1 regulates muscle regeneration by modulating drosha activity. Elife. (2019) 8:e49485. doi: 10.7554/eLife.49485.023
46. Salamon M, Millino C, Raffaello A, Mongillo M, Sandri C, Bean C, et al. Human MYO18B, a novel unconventional myosin heavy chain expressed in striated muscles moves into the myonuclei upon differentiation. J Mol Biol. (2003) 326:137–49. doi: 10.1016/S0022-2836(02)01335-9
47. Tang X, Wang J, Zhou S, Zhou J, Jia G, Wang H, et al. MiR-760 regulates skeletal muscle proliferation in rheumatoid arthritis by targeting Myo18b. Mol Med Rep. (2019) 20:4843–54. doi: 10.3892/mmr.2019.10775
48. Gnimassou O, Francaux M, Deldicque L. Hippo pathway and skeletal muscle mass regulation in mammals: A controversial relationship. Front Physiol. (2017) 8:190. doi: 10.3389/fphys.2017.00190
49. Bassel-Duby R, Olson E N. Signaling pathways in skeletal muscle remodeling. Annu Rev Biochem. (2006) 75:19–37. doi: 10.1146/annurev.biochem.75.103004.142622
50. Keren A, Tamir Y, Bengal E. The p38 MAPK signaling pathway: A major regulator of skeletal muscle development. Mol Cell Endocrinol. (2006) 252:224–30. doi: 10.1016/j.mce.2006.03.017
51. Cheng X, Du J, Shen L, Tan Z, Jiang D, Jiang A, et al. MiR-204-5p regulates C2C12 myoblast differentiation by targeting MEF2C and ERRγ. Biomed Pharmacother. (2018) 101:528–35. doi: 10.1016/j.biopha.2018.02.096
52. Yuan Z R, Xu S Z. Novel SNPs of the bovine CACNA2D1 gene and their association with carcass and meat quality traits. Mol Biol Rep. (2011) 38:365–70. doi: 10.1007/s11033-010-0117-z
53. Gawlik K I, Durbeej M. Skeletal muscle laminin and MDC1A: Pathogenesis and treatment strategies. Skelet Muscle. (2011) 1:9. doi: 10.1186/2044-5040-1-9
Keywords: RNA editome, skeletal muscle, development, yak, transcriptome
Citation: Wu X, Chu M, Ma X, Pei J, Xiong L, Guo X, Liang C and Yan P (2022) Genome-Wide Identification of RNA Editing Sites Affecting Muscle Development in Yak. Front. Vet. Sci. 9:871814. doi: 10.3389/fvets.2022.871814
Received: 08 February 2022; Accepted: 07 June 2022;
Published: 28 June 2022.
Edited by:
Kui Li, Agricultural Genomics Institute at Shenzhen (CAAS), ChinaReviewed by:
Jiabo Wang, Southwest Minzu University, ChinaHui Wang, Southwest Minzu University, China
Copyright © 2022 Wu, Chu, Ma, Pei, Xiong, Guo, Liang and Yan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chunnian Liang, Y2h1bm5pYW4yMDA2JiN4MDAwNDA7MTYzLmNvbQ==; Ping Yan, cGluZ3lhbmx6JiN4MDAwNDA7MTYzLmNvbQ==