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.

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

Population Structure, and Selection Signatures Underlying High-Altitude Adaptation Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Cattle

Copy number variation of the HPGDS gene in the Ashidan yak and its associations with growth traits

Analysis of genomic copy number variations through whole-genome scan in Yunling cattle

Association analysis of KMT2D copy number variation as a positional candidate for growth traits

Impact of STAT5A- CNVs on growth traits in Chinese beef cattle breeds

Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep

Genome-wide Analysis Reveals Differential Selection Involved with Copy Number Variation in Diverse Chinese Cattle.

Genomic sequencing analysis reveals copy number variations and their associations with economically important traits in beef cattle

Effect of copy number variation of PLA2G2A gene to growth traits in Chinese cattle

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

Identification of Candidate Genes Associated with Yak Body Size Using a Genome-Wide Association Study and Multiple Populations of Information

The Identification of Goat KCNJ15 Gene Copy Number Variation and Its Association with Growth Traits

Whole-Genome Resequencing Highlights the Unique Characteristics of Kecai Yaks

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

Genome–Wide Copy Number Variation Analysis in Chinese Hainan Yellow Cattle Using Whole–Genome Re–Sequencing Data

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

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

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

Copy number variation of bovine SHH gene is associated with body conformation traits in Chinese beef cattle

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