The purpose of this study was to analyze the transcriptome of MyoD1 gene knockout MDBK cells (bovine kidney cells) using high-throughput sequencing. For the first time, CRISPR/CAS9 technology was used to construct a MyoD1 knockout in MDBK cells and transcriptome sequence analysis was used to examine MyoD1-related target gene expression. Transcriptome sequencing indicated the presence of 723 differentially expressed genes (DEGs) by comparing wild type and MyoD1 knockout MDBK cells and included 178 upregulated and 72 downregulated genes. The DEGs are mainly enriched in Pl-3-kinase and AKT, p53 signaling pathways. Quantitative RT-PCR confirmed that PDE1B, ADAMTS1, DPT, and CCND2 were highly expressed in the leg muscle, longissimus dorsi, and shoulder of Guanling cattle, and CCND2 was inhibited after MyoD1 knockout, suggesting it may be a key downstream gene of MyoD1 and associated with muscle formation and differentiation in Guanling cattle. This provides experimental data for subsequent studies on the regulatory mechanisms of muscle differentiation in Guanling cattle.

Identifying genetic markers of economically valuable traits has practical benefits for the meat goat industry. To better understand the genomic variations influencing body conformation traits, a genome-wide association study was performed on Tashi goats, an indigenous Chinese goat breed. A total of 155 Tashi goats were phenotyped for eight body conformation traits: body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Then, 100 Tashi goats were randomly selected for whole-genome sequencing and genotyped. We obtained 1676.4 Gb of raw data with an average sequencing depth of 6.2X. Clean reads were aligned to the ARS1.2 reference genome, and 11,257,923 single nucleotide polymorphisms (SNPs) were identified. The structure analysis showed that these Tashi goats were almost not genetically related. The 109, 20, 52, 14, 62, 51, 70, and 7 SNPs were significantly associated with body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Within the ±500 kb region of significant SNPs, 183 genes were annotated. The most significantly enriched KEGG pathway was “olfactory transduction”, and the most significantly enriched gene ontology (GO) terms were “cellular process”, “cellular anatomical entity”, and “molecular transducer activity”. Interestingly, we found several SNPs on chromosomes 10 and 11 that have been identified multiple times for all eight body conformation traits located in two fragments (114 kb and 1.03 Mb). In chr.10:25988403-26102739, the six SNPs were tightly linked, the TACTAG genotype was the highest at 91.8%, and the FNTB (Farnesyltransferase, CAAX Box Beta) and CHURC1 (Churchill Domain Containing 1) genes were located. In chr.11:88216493-89250659, ten SNPs were identified with several dependent linkage disequilibrium (LD) blocks, and seven related genes were annotated, but no significant SNP was located in them. Our results provide valuable biological information for improving growth performance with practical applications for genomic selection in goats.