DNA methylation is an important epigenetic mechanism involved in expression of genes in many biological processes including muscle growth and development. Its effects on economically important traits are evinced from reported significant differences in meat quality traits between Japanese black (Wagyu) and Chinese Red Steppes cattle, thus presenting a unique model for analyzing the effects of DNA methylation on these traits. In the present study, we performed whole genome DNA methylation analysis in the two breeds by whole genome bisulfite sequencing (WGBS). Overall, 23150 differentially methylated regions (DMRs) were identified which were located in 8596 genes enriched in 9922 GO terms, of which 1046 GO terms were significantly enriched (p<0.05) including lipid translocation (GO: 0034204) and lipid transport (GO: 0015914). KEGG analysis showed that the DMR related genes were distributed among 276 pathways. Correlation analysis found that 331 DMRs were negatively correlated with the expression levels of differentially expressed genes (DEGs) with 21 DMRs located in promoter regions. Our results identified novel candidate DMRs and DEGs correlated with meat quality traits, which will be valuable for future genomic and epigenomic studies of muscle development and for marker assisted selection of meat quality traits.

Milk fat is a key factor affecting milk quality and is also a major trait targeted in dairy cow breeding. To determine how the synthesis and the metabolism of lipids in bovine milk is regulated at the miRNA level, primary mammary epithelial cells (pMEC) derived from two Chinese Holstein dairy cows that produced extreme differences in milk fat percentage were cultured by the method of tissue nubbles culture. Small RNA libraries were constructed from each of the two pMEC groups, and Solexa sequencing and bioinformatics analysis were then used to determine the abundance of miRNAs and their differential expression pattern between pMECs. Target genes and functional prediction of differentially expressed miRNAs by Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes analysis illustrated their roles in milk fat metabolism. Results show that a total of 292 known miRNAs and 116 novel miRNAs were detected in both pMECs. Identification of known and novel miRNA candidates demonstrated the feasibility and sensitivity of sequencing at the cellular level. Additionally, 97 miRNAs were significantly differentially expressed between the pMECs. Finally, three miRNAs including bta-miR-33a, bta-miR-152 and bta-miR-224 whose predicted target genes were annotated to the pathway of lipid metabolism were screened and verified by real-time qPCR and Western-blotting experiments. This study is the first comparative profiling of the miRNA transcriptome in pMECs that produce different milk fat content.

DNA methylation is known to play an important role in regulating gene expression during biological development and tissue differentiation in eukaryotes. In this study, we used the fluorescence-labeled methylation-sensitive amplified polymorphism (F-MSAP) method to assess the extent and pattern of cytosine methylation in muscle, heart, liver, spleen, lung, kidney and stomach from the swine strain Laiwu, and we also examined specific methylation patterns in the seven tissues. In total, 96,371 fragments, each representing a recognition site cleaved by either or both EcoRI + HpaII and EcoRI + MspI, the HpaII and MspI are isoschizomeric enzymes, were amplified using 16 pairs of selective primers. A total of 50,094 sites were found to be methylated at cytosines in seven tissues. The incidence of DNA methylation was approximately 53.99% in muscle, 51.24% in the heart, 50.18% in the liver, 53.31% in the spleen, 51.97% in the lung, 51.15% in the kidney and 53.39% in the stomach, as revealed by the incidence of differential digestion. Additionally, differences in DNA methylation levels imply that such variations may be related to specific gene expression during tissue differentiation, growth and development. Three types of bands were generated in the F-MSAP profile, the total numbers of these three types of bands in the seven tissues were 46,277, 24,801 and 25,293, respectively. In addition, different methylation patterns were observed in seven tissues from pig, and almost all of the methylation patterns detected by F-MSAP could be confirmed by Southern analysis using the isolated amplified fragments as probes. The results clearly demonstrated that the F-MSAP technique can be adapted for use in large-scale DNA methylation detection in the pig genome.