Next Article in Journal
Population Structure of a Worldwide Collection of Tropical Japonica Rice Indicates Limited Geographic Differentiation and Shows Promising Genetic Variability Associated with New Plant Type
Next Article in Special Issue
MicroRNA-181a Regulates the Proliferation and Differentiation of Hu Sheep Skeletal Muscle Satellite Cells and Targets the YAP1 Gene
Previous Article in Journal
The Influence of 5-HTTLPR, BDNF Rs6265 and COMT Rs4680 Polymorphisms on Impulsivity in Bipolar Disorder: The Role of Gender
Previous Article in Special Issue
Integrated Hair Follicle Profiles of microRNAs and mRNAs to Reveal the Pattern Formation of Hu Sheep Lambskin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Single Nucleotide Polymorphisms and Haplotype of MYF5 Gene Are Associated with Body Measurements and Ultrasound Traits in Grassland Short-Tailed Sheep

1
Genetic Resources Center, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Animal Disease Prevention and Control Center, Ewenki Autonomous Banner, Hulunbuir 021000, China
*
Author to whom correspondence should be addressed.
These authors contribute equally to this work.
Genes 2022, 13(3), 483; https://doi.org/10.3390/genes13030483
Submission received: 11 February 2022 / Revised: 7 March 2022 / Accepted: 7 March 2022 / Published: 9 March 2022

Abstract

:
Myogenic factor 5 plays active roles in the regulation of myogenesis. The aim of this study is to expose the genetic variants of the MYF5 and its association with growth performance and ultrasound traits in grassland short-tailed sheep (GSTS) in China. The combination technique of sequencing and SNaPshot revealed seven SNPs in ovine MYF5 from 533 adult individuals (male 103 and female 430), four of which are novel ones located at g.6838G > A, g.6989 G > T, g.7117 C > A in the promoter region and g.9471 T > G in the second intron, respectively. Genetic diversity indexes showed the seven SNPs in low or intermediate level, but each of them conformed HWE (p > 0.05) in genotypic frequencies. Association analysis indicated that g.6838G > A, g.7117 C > A, g.8371 T > C, g.9471 T > G, and g.10044 C > T had significant effects on growth performance and ultrasound traits. The diplotypes of H1H3 and H2H3 had higher body weight and greater body size, and haplotype H3 had better performance on meat production than the others. In addition, the dual-luciferase reporter assay showed that there are two active regions in the MYF5 promoter located at −1799~−1197 bp and −514~−241 bp, respectively, but g.6838G > A and g.7117 C > A were out of the region, suggesting these two SNPs influence the phenotype by other pathway. The results suggest that the MYF5 gene might be applied as a promising candidate of functional genetic marker in GSTS breeding.

1. Introduction

The myogenic regulatory factors (MRFs) are essential for regulating the development of muscle fiber and then meat production capacity in mammals [1,2], which include four conserved basic helic-loop-helic (bHLH) molecules functioning as transcription factors as myogenic regulatory factor D (MyoD), myogenic regulatory factor G (MyoG), myogenic factor 5 (MYF5), and myogenic factor 6 (MYF6). All of these four factors affect the determination and maturation of muscle fibers [3]. Variants of these genes can therefore be associated with phenotype of growth, and carcass and meat traits in livestock. Among the four MRFs genes, MYF5 is the earliest one expressed during embryonic muscle development [4]. Experiments evidently lacking MYF5 delayed differentiation and modestly impaired proliferation of satellite cells in MYF5-null knock-out mice [5], and the fetuses exhibited skeletal defects and died immediately after birth in another experiment [6].
Research about the effect of variants in MYF5 on animal production traits has been explored in pigs and cattle. In pigs, variation of MYF5 has been found to be associated with meat quality including loin-eye area, tenderness, total fiber number, water moisture content, and pH value of longissimus dorsi [7,8,9,10]. While in cattle, variation in MYF5 affects growth and meat quality traits including birth, fat, carcass weight, and body length [11,12,13]. In addition, variants of MYF5 have been reported to be associated with growth and reproductive traits in chicken [14], and showed differences in carcass and meat quality traits in pigeon [15] and rabbit [16].
Ovine MYF5 maps to chromosome 3 (OAR3) [17]. Currently, there is a report [1] on variants of ovine MYF5, which are associated with carcass traits in Romney sheep. No other study exists on ovine MYF5 describing its variation in performance, especially in body measurements and ultrasound traits.
Given that the association of MYF5 polymorphism is similar to other livestock in growth and meat quality traits, we hypothesized that there are polymorphisms in MYF5 in ovine that are associated with the traits of growth and meat quality. The objective of this study was therefore to screen for variations in the promoter, 5′UTR, exons, introns, and 3′UTR regions in the MYF5 gene of grassland short-tailed sheep (GSTS) in Inner Mongolia, China, and to explore the presence of associations between variants of MYF5 and growth (body weight and body measurements) and ultrasound traits for revealing its molecular mechanism preliminarily.

2. Materials and Methods

2.1. Experimental Animals

A total of 533 adult GSTS (430 ewes and 103 rams) in Inner Mongolia, China were randomly selected from their breeding farm in Ewenki Autonomous Banner, Inner Mongolia Autonomous Region. These sheep were grazed in the steppe all year round, moved freely, and were in good health.

2.2. Data Collection and Sampling

The blood samples were collected from the jugular vein of sheep (5 mL/sheep) by a qualified veterinarian and stored in blood collection vessels containing EDTA. All operations followed the guideline of Ai Yi Ti Breeding farm of grassland short-tailed sheep, Inner Mongolia Autonomous Region (approval number: AYT-2020112003). Genomic DNA was extracted using the TIANamp Genomic DNA Kit (TIANGEN, Beijing, China) following the manufacturer’s instruction. The DNA purity was tested by Spectrophotometer and dissolved in TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8.0), and stored at −20 °C.
Growth traits including body weight (BW), body length (BL), withers height (WH), chest depth (CD), chest circumference (CC), chest width (CW), cannon bone circumference (CBC), and hip width (HW) were determined as is routine method. Ultrasound traits including eye muscle area (EMA) and backfat thickness (BFT) were measured by the animal specific B ultrasonic instrument whilst the animals were alive (Dawei, Xuzhou, China). Animals to be tested were restrained in a standing position after shaving of the wool and cleaning of the skin on the right side of the region between the 12th and 13th thoracic vertebrae; the ultrasound images of the longissimus muscle were captured using the procedure described by Rozanski et al. [18]. At the end of the measuring, all images of each individual were analyzed by the same person who operated the ultrasound apparatus, using the software fixed in the instrument in order to eliminate systematical error.

2.3. Identifying SNPs and Genotyping

According to the ovine MYF5 sequence provided in Genbank in NCBI (Ovis aries. NC_040254.1), seven pairs of sequential overlapping primers were designed using Primer Premier 5, each pair of primers can amplify about 800 bp in length, and are named MYF5-1, MYF5-2, and MYF5-3, etc., in numerical order, covering the open reading frame (ORF), promoter, 5’UTR, and 3’UTR. The information of the primers is shown in Supplementary Materials Table S1.
For SNPs detection, twenty genomic DNA samples were selected randomly as a pool to identify if mutation of the gene was harbored. PCR was performed in a 25 μL reaction volume, including 2 × Taq PCR MasterMix (12.5 μL), forward and reverse primer (10 μM, 1 μL for each), DNA (1 μL), and ddH2O (9.5 μL). The reaction of thermal cycling included: 95 °C for 5 min, followed by 35 cycles: 95 °C for 30 s, 60 °C for 40 s, and 72 °C for 1 min; and 72 °C for 7 min for the final extension, then the samples were stored at 4 °C. Bidirectional sequencing was operated by Sangon (Shanghai, China). The Chromas software was employed to analyze and determine the location of the SNPs.
After sequencing detection, the technique of SNaPshot was used to type each specific location of SNPs found and then the other seven primers were designed according to the description of Turner et al. (2002), labeled as SNP1, SNP2, and SNP3 etc., in numerical order. The information about the primer pairs is shown in Supplementary Materials Table S2. The PCR reaction mixture was 25.0 μL, including 1.1 × T3 Super PCR Mix 22 μL (TSINGKE, Beijing, China), 10 μM primer F 1 μL, 10 μM Primer R 1 μL, and gDNA 1 μL. The thermal cycling procedure was 98 °C for 2 min, followed by 35 cycles at 98 °C for 10 s, 60 °C for 10 s, and 72 °C for 20 s; then 72 °C for 5 min, and the products were stored in 4 °C.

2.4. The Determination of Transcriptional Regulation of SNPs in the Promoter Region

2.4.1. Cell Culture

The 293T cell line (HEK293T) was purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA), and the cells were cultured in Dulbecco’s modified Eagle medium (DMEM) with 10% fetal bovine serum (FBS) (Ioway, Beijing, China) and 1% penicillin and streptomycin (Invitrogen, Carlsbad, CA, USA). Sheep embryonic fibroblasts cell line (SEF) was courtesy of Professor Ma RZ from the University of Chinese Academy of Science, and the selection of the sheep embryonic fibroblasts was conducted with DMEM/F12 medium. Both of the cells were inoculated into 24-well plates at the density of 1 × 105 cells/well and cultured in 5% CO2 incubator at 37 °C. Cells were treated with 0.1 μM Plicamycin as SP1 inhibitor (Cat#HY-A0122, MedChemExpress, Monmouth Junction, NJ, USA) followed by transfection for 6 h.

2.4.2. Construction of Dual-Luciferase Reporter Expression Vectors

According to CE Design of Vazyme primer software (http://www.vazyme.com, accessed on 22 March 2021), five fragment primers starting in +59 bp with successive deletions of −1799, −1197, −866, −514, and −241 in the promoter region of MYF5 were designed (the primers’ information is shown in Supplementary Materials Table S3), and double restriction sites of XhoIand HandIII were added at each ends of the fragment as PGL4.2-Basic vector request, which were named as PGL4.2-1799, PGL4.2-1197, PGL4.2-866, PGL4.2-514, and PGL4.2-241, respectively. Each insertion of fragments with the same 25 base pairs at 5′ and 3′ ends corresponding to the PGL4.2-Basic was homologously recombined using Clon Express® Ultra One Step Cloning Kit (Vazyme, Nanjing, China), and the ligated products were transformed into DH5α competent cells.

2.4.3. Dual-Luciferase Reporter Assay

The 293. T and SEF cells were inoculated into a 24-well plate at a density of 1 × 105 cells per well initially, and cultured in a 5% CO2 incubator at 37 °C. The cells were transfected when the fusion reached approximately 80% within 24 h. A total of six recombinant plasmids (including one empty plasmid) with 1000 ng and TK-Renilla plasmid with 15 ng were co-transfected into 293T and SEF cells, and each of experiments were repeated three times. After 48 h transfection the cells were collected, and the firefly luciferase (F) and sea kidney luciferase activity (R) were determined by the Dual-Luciferase Reporter® Assay System (Promega, Madison, WI, USA) and detected by GloMax 20/20 Luminometer (Promega, Madison, WI, USA). The firefly luciferase activity levels in each well were normalized by those for Renilla luciferase. All of the experiments were replicated three times and each assay was performed thrice.

2.4.4. Prediction of Transcription Factors in the Promoter Region

Alibaba 2.1 (http://gene-regulation.com/pub/programs/alibaba2/index. Html, accessed on 13 July 2021) online tool was paired with the input complete sequence of MYF5 containing SNP1, 2, and 3 from NCBI firstly, the specific parameters set as: the predicted length is 50 bp, the longest binding site in length is 10 bp, and the shortest is 4 bp, with the confidence above 80%. Then, the predicted transcription factors were further analyzed on the JASPER website (http://jaspar.genereg.net/, accessed on 13 July 2021), both with highest scores were selected.

2.4.5. Construction Luciferase Reporter Plasmid with Mutation

To investigate the difference of luciferase activity among different genotypes of SNP2, a 48 bp (−695 to −742 bp) wild-type sequence containing SNP2 was inserted into PGL4.2-Basic as PGL4.2-SNP2-W, then the site G in wild-type was mutated into T by PCR method to construct a mutant plasmid as PGL4.2-SNP2-M. The constructed plasmid was verified by double digestion and Sanger sequencing. A blank control and the experimental groups were tested for three times in SEF cells. The information of the base mutation primers is shown in Supplementary Materials Table S4.

2.5. Statistical Analysis

The genotype and allele frequency were calculated directly. The polymorphism information content (PIC) was used to estimate the MYF5 in population.
PIC = 1 i = 1 m p i 2 i = 1 m 1 j = i + 1 m 2 p i 2 p j 2 .
p i and p j are the frequencies of i and j alleles, respectively, and n is the number of multiple alleles. PIC > 0.5 indicates that the locus is highly polymorphic, when PIC < 0.25 means the locus is in low polymorphism condition, and 0.25 < PIC < 0.5 suggests that the locus is moderately polymorphic [19].
The Hardy–Weinberg equilibrium (HWE) was evaluated by the Chi-squared test in Haploview 4.2 (Whitehead Institute for Biomedical Research, Cambridge, MA, USA). The r2 and D’ linkage disequilibrium (LD) and haplotype analysis were also calculated with Haploview 4.2.
In this experiment, the general linear model (GLM) as Yij = μ + Gi + Sj + eij was used for analyzing the association between SNPs and phenotypes by SAS 9.4 (SAS Institute Inc., Cary, NC, USA), where Yij was the phenotypic observations; μ was the averaged values; Gi was the fixed effect of genotype; Sj was the fixed effect of sex; and eij was the residual effect. All values were expressed as the mean ± standard deviation. The results with p < 0.05 were considered statistically significant.

3. Results

3.1. SNP Detection and Genetic Parameters

A total of 7 SNPs have been detected in this study by direct sequencing the pool PCR products from seven sequential primers for MYF5 (Figure 1), SNP1, SNP2, and SNP3 locate in the promoter, SNP4 and SNP5 in the first intron, SNP6 in the second intron and SNP7 in the 3′UTR region of the gene; the detailed information is shown in Table 1 and Figure 1. Based on the information above, 533 samples were genotyped by SNaPshot technology. Genetic parameters of genotype and allele frequency, observed heterozygosity (ObsHe), predicted heterozygosity (PredHe), χ2 test for Hardy–Weinberg equilibrium (HWE), and the polymorphism information content (PIC) were calculated for all 7 SNPs and shown in Table 1. All SNPs were in Hardy–Weinberg equilibrium and the genetic diversity of SNP1, SNP5, SNP6, and SNP7 were in low polymorphism condition (PIC < 0.25), and SNP2, SNP3, and SNP4 were in moderate status (0.25 < PIC < 0.5).

3.2. Linkage Disequilibrium and Haplotype of the SNPs

The linkage disequilibrium coefficients of D’ and r2 values (D’ = 1 indicates full linkage, r2 > 0.33 indicates strong linkage) showed that between SNP4 and SNP5 (D’ = 0.97, r2 = 0.73), SNP6 and SNP7 (D’ = 0.82, r2 = 0.64) and SNP1 and SNP2 (D’ = 0.98, r2 = 0.22) were in strong linkages states, and the degree of linkage went from high to low. As Figure 1c displayed, these 7 SNPs formed three haplotype modules in STGS flock.
The haplotype analysis indicated there were eight haplotypes with the frequencies more than 1% in the sampled sheep flock, and H1(-GTATGTC-) was the most abundance one with the frequency of 49.7%, followed by H2 (-GGCTGTC-) and H3 (-AGACAGT-) with 22.2% and 12.6%, respectively (Table 2). Based on the pairwise combination of the eight haplotypes, four major combinations of diplotypes with the frequencies above 5% were identified, which account for 66.3% of the sampled sheep covering 373 individuals, and the other combinations were too small to be considered. Among four combinations, H1H1 was the most common one (GTATGTC/GTATGTC, n = 144, 25.6%) in the tested flock, followed by H1H2 (GTATGTC/GGCTGTC, n = 128, 22.4%), H1H3 (GTATGTC/AGACAGT, n = 59, 10.5%), and H2H3 (GGCTGTC/AGACAGT, n = 42, 7.5%) (Table 2 and Table 3).

3.3. Association of the SNPs and Haplotype Combinations with Phenotype Traits

Through genotyping, three mutations were determined as SNP1, SNP2, and SNP3, which existed in the promoter region. For the SNP1, the sheep with AA genotype has significantly higher BW, WH, CC, and SC than GA and GG ones (p < 0.05). Moreover, individuals of AA genotype had significantly higher CW than that of GG genotype (p < 0.05). For SNP3, the EMA of AC genotype individuals was significantly higher than that of CC (p < 0.05), while there was no difference between CC and AA (p > 0.05). However, there was no significant association detected between the mutants and the phenotypes for SNP2 (Table 4).
The other SNPs, namely SNP4, SNP5, and SNP6, are located in the intron regions of MYF5 and had complicated association with phenotype traits. SNP4 is located in the first intron, and the heterozygous genotype CT had more advantages in BW, BL, WH, CC, CBC, and EMA than CC genotypes (p < 0.05). TT genotype also had higher mean values of BW, CBC, and EMA than CC genotype. Meanwhile, CC, CT genotype had significantly higher WH than CC and TT genotypes (p < 0.05), but there was no significant difference between CC and TT genotypes (p > 0.05).
SNP5 was also located in the first intron, but there is no significant difference among genotypes (p > 0.05). For SNP6, it was located in the second intron of MYF5 gene, and the individuals with GG genotype had significantly higher values in BW, WH, and CBC than individuals with TT genotype (p < 0.05). Genotype GT was in moderate level in all phenotype traits, except that CBC was significantly lower than GG genotype.
For the 3′UTR mutation of SNP7, the sheep with TT genotype had higher BW and greater WH than CT and CC genotype (p < 0.05). Besides, TT genotype had significant higher mean values of chest circumference (CC) than CC genotype, and there was no difference (p > 0.05) between CT and CC genotype.
The result of association analysis amongst the four diplotypes and phenotypic traits showed that the mean value of BW, BL, WH, CC, CW, CBC, and HW of H1H3 and H2H3 were significantly (p < 0.05) higher than those of H1H1 and H1H2. However, there were no significant differences (p > 0.05) among the four diplotypes in two ultrasound traits (EMA and BFT) (Table 5).

3.4. An Optimal Haplotype for Meat Production Breeding

According to the results above, it seems that individuals with haplotype H3 had a better performance than the others, so all individuals were divided into two groups, which were animals with H3 haplotype (H3+) and animals without H3 haplotype (H3−). The result of contrastive analysis showed that H3+ group had significant (p < 0.05) advantages over H3- group in BW, BL, WH, and CW (Table 6).

3.5. The Position and Binding Factor of Three SNPs Functioned in Ovine MYF5 Promoter

After successfully co-transfecting the plasmids of successive truncating fragments of ovine MYF5 promoter with Renilla plasmids into 293T or SEF cells in 24 h, the dual-luciferase activity detecting results tended to be consistent in the two cell lines (only the result of SEF is shown), which indicated that all five truncation fragments are in active, and the fragment from −514 to −241 bp showed the highest activity (p < 10−5) (Figure 2a); this indicated, therefore, that the functional core region of ovine MYF5 gene promoter is located between −514 to −241 bp. However, it is a little surprising that SNP1, 2, and 3 are out of the core region. Certainly, the fragment containing SNP 1, 2, and 3 (−1799 to −866 bp) also had a significant effect (p < 0.001) with the promoter activity (Figure 2a). It is logical to assume that these three SNPs play an important role in regulating the promoter activity.
The results of prediction for the binding sites to transcription factor on the three SNPs showed that SP1 is located in the binding site of SNP2, and there are no binding sites on SNP1 and 3. Hence, the SP1 inhibitor assay conducted, and the result displayed that the luciferase activity decreased significantly (p < 0.001) when SP1 inhibitor added to the culture of SEF cells transfected with PGL4.2-SNP2-W plasmid (Figure 2b). Further, the mutant of PGL4.2-SNP2-M plasmid was transfected, and the activity was significantly (p < 0.001) lower compared to the wild-type (Figure 2c). This suggests that the binding of SP1 to SNP2 plays an important role in exerting the promoter activity of MYF5.

4. Discussion

There were few studies on variations in ovine MYF5, especially the variations associated with the production traits in sheep. Sousa et al. [20] had genotyped MyoD family genes including MYF5 in 192 Santa Inês sheep and also analyzed association with meat quality traits, and had not found an association with MYF5. Phua and Wood [21] firstly reported two variations in ovine MYF5 in 24 sheep samples using Southern hybridization RFLP method. Wang et al. [1] detected two SNPs in exon 3 and an indel in intron 2 from five breeds with total of 645 sheep samples using PCR-SSCP analysis, and found genotype AA from Romney sheep associated with carcass lean meat yield. More MYF5 sequence variation has been described in cattle [12,13,22] and pigs [10,23,24].
In the study, the ovine MYF5 gene was genotyped with seven primers spanning 4606 bp within an indigenous Chinese breed named GSTS and seven SNPs were determined from 533 individuals located in the promoter (SNP1, SNP2, and SNP3), intron 1 (SNP4 and SNP5), intron 2 (SNP6), and 3′UTR (SNP7) of MYF5, respectively. Among them, except for three SNPs in the intron 1 and 3′UTR, which were reported previously in Santa Inês sheep [20], the SNP1, SNP2, SNP3, and SNP6 were detected firstly in Ovis aries. Except for SNP2, this study also revealed significant association (p < 0.05) between SNPs, body measurements, and ultrasound traits in a relatively large amount of sheep samples. This might be the reason that previous studies failed to find the variations in ovine MYF5 gene and the association with measurements of phenotype. The results above further confirmed that MYF5 gene had an effect on the body measurements traits of domestic animals.
Usually, a single SNP has little effect on a phenotype, while haplotype could provide more abundant information, because haplotype composite heterozygous SNPs are located on single chromosomes, play an important role in the mining of disease genes, and search for new strategies for disease treatment [25]. In the present study, we detected eight haplotypes consisting of seven mutant sites with a frequency above 1%, and the most common haplotype is H1 (0.497) in the GSTS population. Through analyzing the association between different haplotypes with the body weight, body size, and ultrasound traits in GSTS flock, we also found that haplotype H3 had a better performance in body weight, body size, and ultrasound traits, hence the combinations of H1H3 and H2H3 showed a higher mean value than other major diplotypes, suggesting that H3 was a beneficial haplotype for body measurements and carcass traits, which laid the foundation for formulating the breeding regime. It is anticipated that strengthening the selection of H3 haplotype would increase the frequency of H3 in the population of the sheep flock, resulting in desirable body measurements traits for improving the meat production performance at an accelerated rate for GSTS.
In the study, seven mutants of SNP were detected in the flock, but no SNPs were found in exon, suggesting that the exon region of MYF5 in GSTS tends to be relatively conserved. It was believed that mutations occurring in intron could not affect biological traits. With the progressing of research, it is found that introns are involved in lots of gene regulation such as splicing in eukaryotes. The variants in introns may lead to change splicing accuracy or efficiency, resulting in changes in amino acid coding and ultimately change in the phenotype of eukaryotes [26]. For example, the variation in the third intron of the IGF2 may affect its binding to the transcription factor ZBED6, which would then improve the meat yield of Chinese Bama pigs [27]. In the study, we detected three SNPs in the intron regions of ovine MYF5, and CC homozygous genotype of the SNP4 had a negative effect on body weight, body measurement traits, and eye muscle area, but CT heterozygous had the highest mean value in all measuring traits among the three genotypes. Comparatively, the mutation of SNP6 had a positive effect on body weight, withers height, and cannon bone circumference. These results supplied clues for further investigation in the mechanisms of transcriptional regulation in gene introns.
Generally, the region of 3′UTR involved the binding site of mRNA and therefore regulated nuclear transport, polyadenylation, and subcellular targeting, as well as rates of localization, translation, and mRNA-specific degradation [28,29]. The function of miRNAs appears to be negative regulators for transcription by interacting with 3′UTR of target gene [30]. Thus, sequence alterations in the region of 3′UTR may involve an inappropriate expression of genes and has been conclusively proven in a lot of studies [29,31,32]. In the current study, we found the sheep with TT genotype at SNP7 in 3′UTR had a higher mean value of body weight, withers height, and chest circumference than CT and CC genotypes. From the routine perceiving, the mutant T at 10,044 bp of MYF5 might be a beneficial allele for muscle development in terms of appearance. Whether this mutation would be favorable to muscle growth without other side effects needs more investigation.
Among the seven SNPs found in the present study, the SNP1, SNP2, and SNP3 are located in the promoter region of ovine MYF5. The luciferase assay determined that SNP2 interacted with transcription factor SP1, SNP1 was significantly associated with the body weight and measurement traits, and SNP3 was significantly associated with eye muscle area. A promoter contains important information about the network of the gene expression regulation [33], and single nucleotide polymorphisms in the promoter regions can affect the gene’s regulation activity by altering the binding efficiency of the transcript factors [34,35]. For instance, Kostek et al. reported that a polymorphism in the IGF1 promoter is significantly associated with total human fat mass [36]. From our study, the results of related SNPs in MYF5 promoter suggest that three SNPs may affect the transcriptional efficiency by changing its promoter regulating activity, and SNP2 specifically may affect the SP1 binding efficiency to regulate the promoter’s activity. Of course, more precise assays such as immunoprecipitation need to be conducted to ensure the interaction between SP1 and the binding site containing SNP2.

5. Conclusions

This is the first report describing the variation of ovine MYF5 in its whole length, and the association between variations and body weight, and body measurements and ultrasound traits in sheep. Through the investigation, seven SNPs were identified in a Chinese indigenous sheep breed and four of them were novel in identification. In addition, g.6838G > A, g.7117 C > A, g.8371 T > C, g.9471 T > G, and g.10044 C > T and diplotypes of H1H3 and H2H3 in GSTS might influence body weight, body measurements, and longissimus dosi ultrasound traits, and haplotype of H3 would be beneficial for improving meat production. These findings laid the foundation for formulating the breeding regime for grassland small-tailed sheep and provided a deep insight into offering a functional genetic marker for improving economically valuable traits in sheep breeding.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13030483/s1, Table S1: Primers and amplification conditions for identifying SNPs; Table S2: SNaPshot single base extension primers; Table S3: Homologous recombination primer designed for construction of expression vector; Table S4: The primer information of the base mutation of SNP2.

Author Contributions

Design of the study, Z.D. and Z.Z.; experimental work, C.L. and W.H.; data collection, Z.Z., C.L., W.Y., S.A. and Y.Z.; data analysis, Z.Z. C.L. and W.H.; data interpretation, Z.D. Z.Z. and Y.Z.; writing—original draft preparation, Z.Z. and C.L.; writing—review and editing, Z.D.; supervision, Z.D.; project administration, Z.D.; funding acquisition, Z.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Strategic Priority Research Program of the Chinese Academy of Sciences, Grant No. (XDA26040305) and the key project of “Prospering Cities with Science and Technology” of Hulunbuir, Grant No. (2021hzzx08).

Institutional Review Board Statement

The study was conducted according to the guidelines of Ai Yi Ti Breeding farm of grassland short-tailed sheep Inner Mongolia Autonomous Region (Approval number: AYT-2020112003).

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, J.Q.; Zhou, H.T.; Forrest, R.H.J.; Hu, J.; Liu, X.; Li, S.B.; Luo, Y.Z.; Hickford, J.G.H. Variation in the ovine MYF5 gene and its effect on carcass lean meat yield in New Zealand Romney sheep. Meat Sci. 2017, 131, 146–151. [Google Scholar] [CrossRef] [PubMed]
  2. Olson, E.N. MyoD family—A paradigm for development—Comment. Genes Dev. 1990, 4, 1454–1461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Hughes, S.M.; Chi, M.M.Y.; Lowry, O.H.; Gundersen, K. Myogenin induces a shift of enzyme activity from glycolytic to oxidative metabolism in muscles of transgenic mice. J. Cell Biol. 1999, 145, 633–642. [Google Scholar] [CrossRef] [PubMed]
  4. Ott, M.O.; Bober, E.; Lyons, G.; Arnold, H.; Buckingham, M. Early expression of the myogenic regulatory gene, MYF-5, in precursor cells of skeletal-muscle in the mouse embro. Development 1991, 111, 1097–1107. [Google Scholar] [CrossRef]
  5. Gayraud-Morel, B.; Chretien, F.; Flarnant, P.; Gomes, D.; Zammit, P.S.; Tajbakhsh, S. A role for the myogenic determination gene Myf5 in adult regenerative myogenesis. Dev. Biol. 2007, 312, 13–28. [Google Scholar] [CrossRef] [Green Version]
  6. Braun, T.; Rudnicki, M.A.; Arnold, H.H.; Jaenisch, R. Targeted inactivation of the muscle regulatory gene MYF-5 results in abnormal rib development and perinatal death. Cell 1992, 71, 369–382. [Google Scholar] [CrossRef]
  7. Liu, M.; Peng, J.; Xu, D.Q.; Zheng, R.; Li, F.E.; Li, J.L.; Zuo, B.; Lei, M.G.; Xiong, Y.Z.; Deng, C.Y.; et al. Association analyses of polymorphisms in porcine MYF5 and MYOD1 genes with carcass traits. Aust. J. Agric. Res. 2007, 58, 1040–1045. [Google Scholar] [CrossRef]
  8. Liu, M.; Peng, J.; Xu, D.Q.; Zheng, R.; Li, F.E.; Li, J.L.; Zuo, B.; Lei, M.G.; Xiong, Y.Z.; Deng, C.Y.; et al. Association of MYF5 and MYOD1 Gene Polymorphisms and Meat Quality Traits in Large White x Meishan F2 Pig Populations. Biochem. Genet. 2008, 46, 720–732. [Google Scholar] [CrossRef]
  9. Wang, W.; Xue, W.; Zhou, X.; Zhang, L.; Wu, J.; Qu, L.; Jin, B.; Zhang, X.; Ma, F.; Xu, X. Effects of candidate genes’ polymorphisms on meat quality traits in pigs. Acta Agric. Scand. Sect. A-Anim. Sci. 2012, 62, 120–126. [Google Scholar] [CrossRef]
  10. Ryu, Y.C.; Lee, E.A.; Chai, H.H.; Park, J.E.; Kim, J.M. Effects of a Novel p.A41P Mutation in the Swine Myogenic factor 5 (MYF5) Gene on Protein Stabilizing, Muscle Fiber Characteristics and Meat Quality. Korean J. Food Sci. Anim. Resour. 2018, 38, 711–717. [Google Scholar] [CrossRef]
  11. Seong, J.; Oh, J.D.; Cheong, I.C.; Lee, K.W.; Lee, H.K.; Suh, D.S.; Jeon, G.J.; Park, K.D.; Kong, H.S. Association between polymorphisms of Myf5 and POU1F1 genes with growth and carcass traits in Hanwoo (Korean cattle). Genes Genom. 2011, 33, 425–430. [Google Scholar] [CrossRef]
  12. Sarti, F.M.; Lasagna, E.; Ceccobelli, S.; Di Lorenzo, P.; Filippini, F.; Sbarra, F.; Giontella, A.; Pieramati, C.; Panella, F. Influence of single nucleotide polymorphisms in the myostatin and myogenic factor 5 muscle growth-related genes on the performance traits of Marchigiana beef cattle. J. Anim. Sci. 2014, 92, 3804–3810. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, C.P.; Raza, S.H.A.; Khan, R.; Sabek, A.; Khan, S.; Ullah, I.; Memon, S.; Abd El-Aziz, A.H.; Shah, M.A.; Li, S.J.; et al. Genetic variants in MYF5 affected growth traits and beef quality traits in Chinese Qinchuan cattle. Genomics 2020, 112, 2804–2812. [Google Scholar] [CrossRef] [PubMed]
  14. Tang, Y.; Zhang, T.; Zhang, G.X.; Wang, J.Y.; Fan, Q.C.; Chen, X.S.; Wei, Y.; Han, K.P.; Wang, Y.J. Eight SNPs of the Myf5 gene and diplotypes associated with growth and reproductive traits in Jinghai yellow chicken. Mol. Biol. Rep. 2014, 41, 6837–6844. [Google Scholar] [CrossRef] [PubMed]
  15. Yin, Z.Z.; Dong, X.Y.; Dong, D.J.; Ma, Y.Z. Association of MYF5 and KLF15 gene polymorphisms with carcass traits in domestic pigeons (Columba livia). Br. Poult. Sci. 2016, 57, 612–618. [Google Scholar] [CrossRef]
  16. Wang, J.; Hu, Y.S.; Elzo, M.A.; Shi, Y.; Jia, X.B.; Chen, S.Y.; Lai, S.J. Genetic effect of Myf5 gene in rabbit meat quality traits. J. Genet. 2017, 96, 673–679. [Google Scholar] [CrossRef]
  17. Broad, T.E.; Burkin, D.J.; Jones, C.; Lewis, P.E.; Ansari, H.A.; Pearce, P.D. Mapping of MYF5, C1R, MYHL, TPI1, IAPP, A2MR and RNR onto sheep chromosome-3Q. Anim. Genet. 1993, 24, 415–419. [Google Scholar] [CrossRef]
  18. Rozanski, S.; Vivian, D.R.; Kowalski, L.H.; Prado, O.R.; Fernandes, S.R.; de Souza, J.C.; de Freitas, J.A. Carcass and meat traits, and non-carcass components of lambs fed ration containing increasing levels of urea. Semin.-Cienc. Agrar. 2017, 38, 1587–1603. [Google Scholar] [CrossRef]
  19. Nei, M.; Roychoudhury, A.K. SAMPLING VARIANCES OF HETEROZYGOSITY AND GENETIC DISTANCE. Genetics 1974, 76, 379–390. [Google Scholar] [CrossRef]
  20. Sousa, L.P.B.; Meira, A.N.; Azevedo, H.C.; Muniz, E.N.; Coutinho, L.L.; Mourao, G.B.; Leao, A.G.; Pedrosa, V.B.; Pinto, L.F.B. Variants in myostatin and MyoD family genes are associated with meat quality traits in Santa Ines sheep. Anim. Biotechnol. 2020, 1–13. [Google Scholar] [CrossRef]
  21. Phua, S.H.; Wood, N.J. The myogenic factor 5 (MYF5) locus in sheep carries a rare Eco RI RFLP. Animal Genetics 1993, 24, 220. [Google Scholar] [CrossRef]
  22. Robakowska-Hyzorek, D.; Oprzadek, J.; Zelazowska, B.; Olbromski, R.; Zwierzchowski, L. Effect of the g.-723G -> T Polymorphism in the Bovine Myogenic Factor 5 (Myf5) Gene Promoter Region on Gene Transcript Level in the Longissimus Dorsi Muscle and on Meat Traits of Polish Holstein-Friesian Cattle. Biochem. Genet. 2010, 48, 450–464. [Google Scholar] [CrossRef]
  23. Kunhareang, S.; Zhou, H.; Hickford, J.G.H. Allelic Variation in the Porcine MYF5 Gene Detected by PCR-SSCP. Mol. Biotechnol. 2009, 41, 208–212. [Google Scholar] [CrossRef]
  24. Verner, J.; Humpolicek, P.; Knoll, A. Impact of MYOD family genes on pork traits in Large White and Landrace pigs. J. Anim. Breed. Genet. 2007, 124, 81–85. [Google Scholar] [CrossRef] [PubMed]
  25. Riazulhaq; Chytil, F. Expression of nuclear retinoic acid receptors in bat adipose tissue. Biochem. Biophys. Res. Commun. 1991, 176, 1539–1544. [Google Scholar] [CrossRef]
  26. Greenwood, T.A.; Kelsoe, J.R. Promoter and intronic variants affect the transcriptional regulation of the human dopamine transporter gene. Genomics 2003, 82, 511–520. [Google Scholar] [CrossRef]
  27. Xiang, G.H.; Ren, J.L.; Hai, T.; Fu, R.; Yu, D.W.; Wang, J.; Li, W.; Wang, H.Y.; Zhou, Q. Editing porcine IGF2 regulatory element improved meat production in Chinese Bama pigs. Cell. Mol. Life Sci. 2018, 75, 4619–4628. [Google Scholar] [CrossRef] [PubMed]
  28. Hesketh, J. 3 ′-Untranslated regions are important in mRNA localization and translation: Lessons from selenium and metallothionein. Biochem. Soc. Trans. 2004, 32, 990–993. [Google Scholar] [CrossRef] [Green Version]
  29. Reamon-Buettner, S.M.; Cho, S.-H.; Borlak, J. Mutations in the 3′-untranslated region of GATA4 as molecular hotspots for congenital heart disease (CHD). BMC Med. Genet. 2007, 8, 38. [Google Scholar] [CrossRef]
  30. Dumortier, O.; Van Obberghen, E. MicroRNAs in pancreas development. Diabetes Obes. Metab. 2012, 14, 22–28. [Google Scholar] [CrossRef]
  31. Shibayama, A.; Cook, E.H.; Feng, J.N.; Glanzmann, C.; Yan, J.; Craddock, N.; Jones, I.R.; Goldman, D.; Heston, L.L.; Sommer, S.S. MECP2 structural and 3′-UTR variants in schizophrenia, autism and other psychiatric diseases: A possible association with autism. Am. J. Med. Genet. Part B-Neuropsychiatr. Genet. 2004, 128B, 50–53. [Google Scholar] [CrossRef]
  32. Yuan, M.; Zhan, Q.; Duan, X.; Song, B.; Zeng, S.; Chen, X.; Yang, Q.; Xia, J. A functional polymorphism at miR-491-5p binding site in the 3 ′-UTR of MMP-9 gene confers increased risk for atherosclerotic cerebral infarction in a Chinese population. Atherosclerosis 2013, 226, 447–452. [Google Scholar] [CrossRef] [PubMed]
  33. Hu, H.; Jia, Q.; Hou, S.; Liu, J.; Zhao, S.; Zhang, W.; Zhang, J.; Zhang, J.; Bian, H. Cloning and Activity Analysis of the Promoter Region of Swine DKK1 Gene. Acta Vet. Et Zootech. Sin. 2017, 48, 1150–1157. [Google Scholar]
  34. Wray, G.A.; Hahn, M.W.; Abouheif, E.; Balhoff, J.P.; Pizer, M.; Rockman, M.V.; Romano, L.A. The evolution of transcriptional regulation in eukaryotes. Mol. Biol. Evol. 2003, 20, 1377–1419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Savinkova, L.; Drachkova, I.; Arshinova, T.; Ponomarenko, P.; Ponomarenko, M.; Kolchanov, N. An Experimental Verification of the Predicted Effects of Promoter TATA-Box Polymorphisms Associated with Human Diseases on Interactions between the TATA Boxes and TATA-Binding Protein. PLoS ONE 2013, 8, e54626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Kostek, M.C.; Devaney, J.M.; Gordish-Dressman, H.; Harris, T.B.; Thompson, P.D.; Clarkson, P.M.; Angelopoulos, T.J.; Gordon, P.M.; Moyna, N.M.; Pescatello, L.S.; et al. A polymorphism near IGF1 is associated with body composition and muscle function in women from the Health, Aging, and Body Composition Study. Eur. J. Appl. Physiol. 2010, 110, 315–324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Structure of the MYF5 gene (a); Sequence variants of the MYF5 gene in GSTS (b); Linkage disequilibrium coefficients (D’ left, R2 right) among MYF5 SNPs (c).
Figure 1. Structure of the MYF5 gene (a); Sequence variants of the MYF5 gene in GSTS (b); Linkage disequilibrium coefficients (D’ left, R2 right) among MYF5 SNPs (c).
Genes 13 00483 g001
Figure 2. Promoter activity analysis of ovine MYF5 gene in sheep embryonic fibroblast (a); The effect of SP1 inhibitor to the PGL 4.2-SNP2-W plasmid (b); Comparison of the relative dual-luciferase activity between wild-type and mutant PGL4.2-SNP2 plasmid (c). Results are shown as mean ± standard deviation and the data are representative of at least three independent assays. The statistically significant difference between groups was tested by independent sample t-test. *** p < 10−3, ***** p < 10−5.
Figure 2. Promoter activity analysis of ovine MYF5 gene in sheep embryonic fibroblast (a); The effect of SP1 inhibitor to the PGL 4.2-SNP2-W plasmid (b); Comparison of the relative dual-luciferase activity between wild-type and mutant PGL4.2-SNP2 plasmid (c). Results are shown as mean ± standard deviation and the data are representative of at least three independent assays. The statistically significant difference between groups was tested by independent sample t-test. *** p < 10−3, ***** p < 10−5.
Genes 13 00483 g002
Table 1. Genetic parameters, genotypic, and allelic frequencies of MYF5.
Table 1. Genetic parameters, genotypic, and allelic frequencies of MYF5.
IDLocusPositionGenotypic Frequency $Allele Frequencyχ2 (HWE)ObsHePredHePIC
SNP1g.6838G > APromoterGG
400/0.750
GA
123/0.231
AA
10/0.019
G
0.866
A
0.134
p > 0.050.2300.2320.205
SNP2g.6989G > TPromoterGG
96/0.180
GT
237/0.445
TT
200/0.375
G
0.402
T
0.598
p > 0.050.4440.4810.365
SNP3g.7106C > APromoterCC
59/0.111
CA
230/0.432
AA
244/0.457
C
0.326
A
0.674
p > 0.050.4310.4400.343
SNP4g.8371T > CIntron1TT
351/0.659
CT
166/0.311
CC
16/0.030
T
0.814
C
0.186
p > 0.050.3110.3020.257
SNP5g.8561G > AIntron1GG
383/0.719
GA
140/0.263
AA
10/0.019
G
0.850
A
0.150
p > 0.050.2620.2550.223
SNP6g.9471T > GIntron2TT
376/0.705
GT
143/0.268
GG
14/0.026
T
0.840
G
0.160
p > 0.050.2660.2680.233
SNP7g.10044C > T3′UTRCC
366/0.687
CT
156/0.293
TT
11/0.066
C
0.833
T
0.167
p > 0.050.2920.2780.240
$ The number before the slash represents the quantity of individuals and after it the frequency; HWE, Hardy–Weinberg equilibrium; ObsHe, observation heterozygosity; PredHe, prediction heterozygosity; PIC, polymorphism information content; 3′UTR, 3′ untranslated region.
Table 2. Haplotypes of MYF5 gene and the frequencies in GSTS flock.
Table 2. Haplotypes of MYF5 gene and the frequencies in GSTS flock.
HaplotypeSNP1SNP2SNP3SNP4SNP5SNP6SNP7Frequency
g.6838G > Ag.6989G > Tg.7106C > Ag.8371T > Cg.8561G > Ag.9471T > Gg.10044C > T
H1GTATGTC0.497
H2GGCTGTC0.222
H3AGACAGT0.126
H4GTACGTC0.038
H5GTCTGTC0.036
H6GGCTGTT0.025
H7GGCTGGC0.019
H8GTCCATC0.017
Table 3. Diplotypes of MYF5 gene and the frequencies in GSTS flock.
Table 3. Diplotypes of MYF5 gene and the frequencies in GSTS flock.
Diplotypes SNP1SNP2SNP3SNP4SNP5SNP6SNP7Sample SizeFrequency
g.6838G > Ag.6989G > Tg.7106C > Ag.8371T > Cg.8561G > Ag.9471T > Gg.10044C > T
H1H1G GT TA AT TG GT TC C1440.256
H1H2G GT GA CT TG GT TC C1280.224
H1H3G AT GA AT CG AT GC T590.105
H2H3G AG GC AT CG AT GC T420.075
Table 4. Association between mutations in MYF5 with body measurements and ultrasound traits.
Table 4. Association between mutations in MYF5 with body measurements and ultrasound traits.
LocusGenotypeBW (kg)BL (cm)WH (cm)CD (cm)CC (cm)CW (cm)CBC (cm)HW (cm)EMA (mm2)BFT (mm)
SNP1AA(10)55.88 ± 15.75 a72.10 ± 5.2376.80 ± 5.75 a40.60 ± 2.6398.6 ± 10.05 a20.40 ± 1.71 a10.15 ± 0.94 a24.00 ± 2.491040.04 ± 376.835.29 ± 2.36
GA(122)51.85 ± 9.59 b73.36 ± 5.4275.31 ± 4.16 b40.85 ± 2.0896.09 ± 6.52 b19.85 ± 1.50 a,b9.84 ± 0.80 b23.31 ± 1.951030.53 ± 294.144.61 ± 1.62
GG(400)50.34 ± 9.88 b72.12 ± 5.4374.52 ± 4.01 b40.42 ± 2.2395.33 ± 6.85 b19.51 ± 1.67 b9.77 ± 0.75 b23.12 ± 1.96995.38 ± 281.724.50 ± 1.57
SNP2GG(96)51.13 ± 10.2872.69 ± 5.0874.82 ± 3.7240.67 ± 2.0995.42 ± 6.9119.59 ± 1.479.79 ± 0.8023.06 ± 1.87979.20 ± 293.434.56 ± 1.65
GT(236)50.66 ± 10.0972.49 ± 5.6074.81 ± 4.1640.51 ± 2.3695.65 ± 6.8719.61 ± 1.699.76 ± 0.8023.17 ± 2.011020.04 ± 308.514.55 ± 1.61
TT(200)50.78 ± 9.7372.16 ± 5.4374.62 ± 4.2040.47 ± 2.0895.54 ± 6.8319.60 ± 1.669.83 ± 0.7223.26 ± 1.97997.74 ± 254.154.51 ± 1.57
SNP3CC(59)50.23 ± 10.4972.64 ± 5.9274.49 ± 3.6240.45 ± 2.3195.13 ± 6.7519.37 ± 1.519.73 ± 0.7222.98 ± 1.72953.17 ± 279.64 b4.46 ± 1.69
AC(229)51.02 ± 9.9272.53 ± 5.3474.84 ± 4.1440.59 ± 2.2495.48 ± 6.7719.58 ± 1.629.75 ± 0.8223.16 ± 1.991018.16 ± 310.42 a4.54 ± 1.59
AA(244)50.71 ± 9.9372.22 ± 5.4374.71 ± 4.1740.48 ± 2.1695.75 ± 6.9619.68 ± 1.689.85 ± 0.7223.25 ± 2.011003.65 ± 263.16 a,b4.55 ± 1.60
SNP4CC(16)47.17 ± 3.92 b70.56 ± 2.7 b74.00 ± 1.67 b40.56 ± 2.4293.43 ± 4.41 b19.68 ± 1.309.46 ± 0.38 b22.62 ± 1.40874.80 ± 158.68 b4.13 ± 0.87
CT(166)52.23 ± 10.0 4 a73.19 ± 5.3 2 a75.54 ± 4.41 a40.92 ± 2.1496.16 ± 6.92 a19.81 ± 1.579.89 ± 0.82 a23.35 ± 1.981036.16 ± 291.39 a4.65 ± 1.66
TT(351)50.29 ± 10.0 4 a72.11 ± 5.55 a,b74.40 ± 3.96 b40.33 ± 2.2195.38 ± 6.88 ab19.50 ± 1.679.76 ± 0.75 a23.13 ± 1.98994.65 ± 286.43 a4.50 ± 1.60
SNP5AA(10)53.33 ± 14.9971.20 ± 5.0775.50 ± 4.8340.40 ± 2.5097.00 ± 9.0920.50 ± 1.64 a10.00 ± 0.9123.70 ± 2.40976.52 ± 349.504.55 ± 1.60
GA(140)52.01 ± 9.5573.21 ± 5.3075.52 ± 4.3640.99 ± 2.1596.03 ± 6.6319.79 ± 1.52 a,b9.84 ± 0.8123.30 ± 1.951034.40 ± 287.754.68 ± 1.66
GG(383)50.29 ± 9.9472.13 ± 5.4874.44 ± 3.9440.36 ± 2.2095.36 ± 6.8619.51 ± 1.67 b9.77 ± 0.7523.13 ± 1.96993.57 ± 283.874.48 ± 1.58
SNP6GG(14)53.60 ± 14.46 a71.78 ± 4.9175.92 ± 5.34 a40.64 ± 2.2797.35 ± 9.4219.85 ± 1.7410.17 ± 0.82 a23.28 ± 2.461046.58 ± 330.024.75 ± 2.14
GT(142)52.04 ± 9.60 a,b73.26 ± 5.4275.37 ± 4.05 a,b40.89 ± 2.1996.09 ± 6.6719.82 ± 1.559.85 ± 0.79 b23.34 ± 1.951025.65 ± 291.464.74 ± 1.67
TT(376)50.26 ± 9.86 b72.11 ± 5.4474.47 ± 4.03 b40.39 ± 2.2095.31 ± 6.8119.51 ± 1.669.76 ± 0.75 b23.12 ± 1.96994.31 ± 283.064.45 ± 1.55
SNP7TT(11)54.97 ± 15.25 a72.27 ± 5.0076.63 ± 5.48 a40.63 ± 2.5098.00 ± 9.74 a20.18 ± 1.7710.04 ± 0.9623.63 ± 2.651037.14 ± 357.625.04 ± 2.39
CT(156)51.20 ± 9.58 b72.96 ± 5.4675.07 ± 4.03 b40.68 ± 2.1995.62 ± 6.59 a,b19.72 ± 1.539.82 ± 0.7923.19 ± 1.871009.6 ± 290.994.60 ± 1.62
CC(366)50.50 ± 9.93 b72.16 ± 5.4374.54 ± 4.06 b40.45 ± 2.2195.47 ± 6.86 b19.53 ± 1.679.77 ± 0.7523.16 ± 1.991000.6 ± 282.574.49 ± 1.57
Values are shown as the least squares means ± standard deviation; a,b Means with different superscripts are significantly different (p < 0.05); body weight (BW), body length (BL), withers height (WH), chest depth (CD), chest circumference (CC), chest width (CW), cannon bone circumference (CBC), hip width (HW). Ultrasound traits including eye muscle area (EMA) and backfat thickness (BFT).
Table 5. Associations of combined genotypes with body measurements and ultrasound traits in GSTS flock.
Table 5. Associations of combined genotypes with body measurements and ultrasound traits in GSTS flock.
Combined GenotypesBW (kg)BL (cm)WH (cm)CD (cm)CC (cm)CW (cm)CBC (cm)HW (cm)EMA (mm2)BFT (mm)
H1H1(144)49.90 ± 9.48 b71.61 ± 5.21 b74.20 ± 3.92 b40.29 ± 2.0195.18 ± 6.91 b19.43 ± 1.72 b9.79 ± 0.68 a,b23.10 ± 2.02 b986.39 ± 253.014.43 ± 1.48
H1H2(128)49.60 ± 9.40 b71.77 ± 5.19 b74.12 ± 3.79 b40.25 ± 2.3295.00 ± 6.58 b19.45 ± 1.72 b9.69 ± 0.77 b23.00 ± 1.97 b1006.77 ± 317.834.50 ± 1.59
H1H3(59)52.49 ± 10.23 a73.67 ± 6.06 a75.72 ± 4.65 a40.91 ± 2.2997.01 ± 6.86 a20.16 ± 1.60 a9.93 ± 0.80 a23.72 ± 2.05 a1047.53 ± 283.504.84 ± 1.77
H2H3(42)51.90 ± 8.25 a73.54 ± 4.55 a75.02 ± 3.47 a40.88 ± 1.8395.50 ± 6.20 a,b19.69 ± 1.35 a,b9.76 ± 0.82 a,b23.11 ± 1.72 b1033.34 ± 289.564.57 ± 1.49
a,b Means with different super-scripts are significantly different (p < 0.05).
Table 6. The effect of haplotype H3 to the body measurements traits and ultrasound traits in GSTS flock.
Table 6. The effect of haplotype H3 to the body measurements traits and ultrasound traits in GSTS flock.
Combined GenotypesBW (kg)BL (cm)WH (cm)CD (cm)CC (cm)CW (cm)CBC (cm)HW (cm)EMA (mm2)BFT (mm)
H3+(131)52.29 ± 10.2273.26 ± 5.4675.47 ± 4.3240.85 ± 2.1296.34 ± 6.8919.89 ± 1.549.87 ± 0.8123.37 ± 2.011037.31 ± 300.134.69 ± 1.68
H3-(402)50.27 ± 9.872.14 ± 5.4174.52 ± 3.9940.43 ± 2.2295.33 ± 6.8319.5 ± 1.669.77 ± 0.7523.12 ± 1.96993.77 ± 281.724.49 ± 1.57
p0.0440.0420.0220.0580.1420.0190.1690.2060.1330.208
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Z.; Liu, C.; Hao, W.; Yin, W.; Ai, S.; Zhao, Y.; Duan, Z. Novel Single Nucleotide Polymorphisms and Haplotype of MYF5 Gene Are Associated with Body Measurements and Ultrasound Traits in Grassland Short-Tailed Sheep. Genes 2022, 13, 483. https://doi.org/10.3390/genes13030483

AMA Style

Zhang Z, Liu C, Hao W, Yin W, Ai S, Zhao Y, Duan Z. Novel Single Nucleotide Polymorphisms and Haplotype of MYF5 Gene Are Associated with Body Measurements and Ultrasound Traits in Grassland Short-Tailed Sheep. Genes. 2022; 13(3):483. https://doi.org/10.3390/genes13030483

Chicago/Turabian Style

Zhang, Zhichao, Cheng Liu, Wenjing Hao, Weiwen Yin, Sitong Ai, Yanfang Zhao, and Ziyuan Duan. 2022. "Novel Single Nucleotide Polymorphisms and Haplotype of MYF5 Gene Are Associated with Body Measurements and Ultrasound Traits in Grassland Short-Tailed Sheep" Genes 13, no. 3: 483. https://doi.org/10.3390/genes13030483

APA Style

Zhang, Z., Liu, C., Hao, W., Yin, W., Ai, S., Zhao, Y., & Duan, Z. (2022). Novel Single Nucleotide Polymorphisms and Haplotype of MYF5 Gene Are Associated with Body Measurements and Ultrasound Traits in Grassland Short-Tailed Sheep. Genes, 13(3), 483. https://doi.org/10.3390/genes13030483

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop