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.

Genome-Wide DNA Methylation and Hydroxymethylation Changes Revealed Epigenetic Regulation of Neuromodulation and Myelination in Yak Hypothalamus

Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production

DNA methylation dynamics during yak adipocyte differentiation

Yak Response to High-Altitude Hypoxic Stress by Altering Mrna Expression and DNA Methylation of Hypoxia-Inducible Factors

A comparative analysis of differentially expressed mRNAs, miRNAs and circRNAs provides insights into the key genes involved in the high-altitude adaptation of yaks

Comparative DNA methylation reveals epigenetic adaptation to high altitude in snub-nosed monkeys

Different Gene Expressions Between Cattle and Yak Provide Insights into High-Altitude Adaptation.

Genome-wide DNA methylation profiles in Tibetan and Yorkshire pigs under high-altitude hypoxia

The transcriptomic landscape of yaks reveals molecular pathways for high altitude adaptation

Transcriptome analysis identified long non-coding RNAs involved in the adaption of yak to high-altitude environments

Whole-Transcriptome Analysis of Yak and Cattle Heart Tissues Reveals Regulatory Pathways Associated With High-Altitude Adaptation

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

Single-cell RNA-seq and single-cell bisulfite-sequencing reveal insights into yak preimplantation embryogenesis

A Conjoint Analysis of Renal Structure and Omics Characteristics Reveal New Insight to Yak High-Altitude Hypoxia Adaptation

Genome Methylation and Regulatory Functions for Hypoxic Adaptation in Tibetan Chicken Embryos.

Genome-wide DNA methylation landscape of four Chinese populations and epigenetic variation linked to Tibetan high-altitude adaptation

Favoring Expression of Yak Alleles in Interspecies F1 Hybrids of Cattle and Yak Under High-Altitude Environments

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

Metabolomics Reveals Potential Plateau Adaptability by Regulating Inflammatory Response and Oxidative Stress-Related Metabolism and Energy Metabolism Pathways in Yak

Comparative Analysis of the Microrna Transcriptome Between Yak and Cattle Provides Insight into High-Altitude Adaptation.

Genome-wide DNA methylation status of Mongolians exhibits signs of cellular stress response related to their nomadic lifestyle