Abstract:DNA methylation modifications are implicated in many biological processes. As the most common epigenetic mechanism DNA methylation also affects muscle growth and development. The majority of previous studies have focused on different varieties of yak, but little is known about the epigenetic regulation mechanisms in different age groups of animals. The development of muscles in the different stages of yak growth remains unclear. In this study, we selected the longissimus dorsi muscle tissue at three different growth stages of the yak, namely, 90-day-old fetuses (group E), six months old (group M), and three years old (group A). Using RNA-Seq transcriptome sequencing and methyl-RAD whole-genome methylation sequencing technology, changes in gene expression levels and DNA methylation status throughout the genome were investigated during the stages of yak development. Each group was represented by three biological replicates. The intersections of expression patterns of 7694 differentially expressed genes (DEGs) were identified (padj < 0.01, |log2FC| > 1.2) at each of the three developmental periods. Time-series expression profile clustering analysis indicated that the DEGs were significantly arranged into eight clusters which could be divided into two classes (padj < 0.05), class I profiles that were downregulated and class II profiles that were upregulated. Based on this cluster analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that DEGs from class I profiles were significantly (padj < 0.05) enriched in 21 pathways, the most enriched pathway being the Axon guidance signaling pathway. DEGs from the class II profile were significantly enriched in 58 pathways, the pathway most strongly enriched being Metabolic pathway. After establishing the methylation profiles of the whole genomes, and using two groups of comparisons, the three combinations of groups (M-vs.-E, M-vs.-A, A-vs.-E) were found to have 1344, 822, and 420 genes, respectively, that were differentially methylated at CCGG sites and 2282, 3056, and 537 genes, respectively, at CCWGG sites. The two sets of data were integrated and the negative correlations between DEGs and differentially methylated promoters (DMPs) analyzed, which confirmed that TMEM8C, IGF2, CACNA1S and MUSTN1 were methylated in the promoter region and that expression of the modified genes was negatively correlated. Interestingly, these four genes, from what was mentioned above, perform vital roles in yak muscle growth and represent a reference for future genomic and epigenomic studies in muscle development, in addition to enabling marker-assisted selection of growth traits.

Molecular Characterization and Expression Analysis of Myf6 Gene in Yak (bos Grunniens)

Polymorphism and Sequence Analysis of Partial Fragment of BMP6 Gene in Five Sheep Breeds

Whole-transcriptome analysis of longissimus dorsi muscle in cattle-yaks reveals the regulatory functions of ADAMTS6 gene in myoblasts

Transcriptome and DNA Methylation Analyses of the Molecular Mechanisms Underlying with Longissimus dorsi Muscles at Different Stages of Development in the Polled Yak

microRNA Temporal-Specific Expression Profiles Reveal longissimus dorsi Muscle Development in Tianzhu White Yak

Integrative ATAC-seq and RNA-seq Analysis of the Longissimus Dorsi Muscle of Gannan Yak and Jeryak

Integration of ATAC-Seq and RNA-Seq Analysis to Identify Key Genes in the Longissimus Dorsi Muscle Development of the Tianzhu White Yak

Cloning and sequencing of GH, H-FABP and HSP72 genes in yak

Gelatinase MMP Expression Is Correlated with Muscle Fiber Growth in Maiwa Yak Gluteus Maximus

Comparative Transcriptome Analysis of mRNA and miRNA during the Development of Longissimus Dorsi Muscle of Gannan Yak and Tianzhu White Yak

Complete characterization of the yak testicular development using accurate full-length transcriptome sequencing

Copy Number Variations of KLF6 Modulate Gene Transcription and Growth Traits in Chinese Datong Yak (Bos Grunniens)

Cloning and Molecular Characteristics of SRY Gene in Yak

Organization and Complexity of the Yak (Bos Grunniens) Immunoglobulin Loci

Mutations in the FOXO3 Gene and Their Effects on Meat Traits in Gannan Yaks

Genetic structure analysis of yak breeds and their response to adaptive evolution

Transcriptomic analysis reveals differentially expressed genes and key immune pathways in the spleen of the yak ( Bos grunniens ) at different growth stage

N6-Methyladenosine Modification-Related Genes Express Differentially in Sterile Male Cattle-Yaks

Delineating Genetic Variability: Comparative Analysis of Yak MHC-DQB Gene through cDNA Cloning

Whole-genome resequencing reveals genetic diversity, differentiation, and selection signatures of yak breeds/populations in southwestern China

Metallothionein in Yak Characterization of Metallothionein-III in Yak (bos Grunniens)