BACKGROUND: Indigenous Chinese cattle have abundant genetic diversity and a long history of artificial selection, giving local breeds advantages in adaptability, forage tolerance and resistance. The detection of selective sweeps and comparative genome analysis of selected breeds and ancestral populations provide a basis for understanding differences among breeds and for the identification and utilization of candidate genes. We investigated genetic diversity, population structure, and signatures of selection using genome-wide sequencing data for a new breed of Qinchuan cattle (QNC, n = 21), ancestral Qinchuan cattle (QCC, n = 20), and Zaosheng cattle (ZSC, n = 19).
RESULTS: A population structure analysis showed that the ancestry components of QNC and ZSC were similar. In addition, the QNC and ZSC groups showed higher proportions of European taurine ancestry than that of QCC, and this may explain the larger body size of QNC, approaching that of European cattle under long-term domestication and selection. A neighbor-joining tree revealed that QCC individuals were closely related, whereas QNC formed a distinct group. To search for signatures of selection in the QNC genome, we evaluated nucleotide diversity (θπ), the fixation index (FST) and Tajima\'s D. Overlapping selective sweeps were enriched for one KEGG pathway, the apelin signaling pathway, and included five candidate genes (MEF2A, SMAD2, CAMK4, RPS6, and PIK3CG). We performed a comprehensive review of genomic variants in QNC, QCC, and ZSC using whole-genome sequencing data. QCC was rich in novel genetic diversity, while diversity in QNC and ZSC cattle was reduced due to strong artificial selection, with divergence from the original cattle.
CONCLUSIONS: We identified candidate genes associated with production traits. These results support the success of selective breeding and can guide further breeding and resource conservation of Qinchuan cattle.

Tropomyosin 3 (TPM3) plays a significant role as a regulatory protein in muscle contraction, affecting the growth and development of skeletal muscles. Despite its importance, limited research has been conducted to investigate the influence of TPM3 on bovine skeletal muscle development. Therefore, this study revealed the role of TPM3 in bovine myoblast growth and development. This research involved conducting a thorough examination of the Qinchuan cattle TPM3 gene using bioinformatics tools to examine its sequence and structural characteristics. Furthermore, TPM3 expression was evaluated in various bovine tissues and cells using quantitative real-time polymerase chain reaction (qRT-PCR). The results showed that the coding region of TPM3 spans 855 bp, with the 161st base being the T base, encoding a protein with 284 amino acids and 19 phosphorylation sites. This protein demonstrated high conservation across species while displaying a predominant α-helix secondary structure despite being an unstable acidic protein. Notably, a noticeable increase in TPM3 expression was observed in the longissimus dorsi muscle and myocardium of calves and adult cattle. Expression patterns varied during different stages of myoblast differentiation. Functional studies that involved interference with TPM3 in Qinchuan cattle myoblasts revealed a very significantly decrease in S-phase cell numbers and EdU-positive staining (P < 0.01), and disrupted myotube morphology. Moreover, interference with TPM3 resulted in significantly (P < 0.05) or highly significantly (P < 0.01) decreased mRNA and protein levels of key proliferation and differentiation markers, indicating its role in the modulation of myoblast behavior. These findings suggest that TPM3 plays an essential role in bovine skeletal muscle growth by influencing myoblast proliferation and differentiation. This study provides a foundation for further exploration into the mechanisms underlying TPM3-mediated regulation of bovine muscle development and provides valuable insights that could guide future research directions as well as potential applications for livestock breeding and addressing muscle-related disorders.

A genome-wide association study (GWAS) is an effective tool for identifying the dominant genes of complex economic traits in livestock by statistical analysis of genotype data and measured phenotype data. In this study, we rigorously measured 14 body conformation traits in 254 Qinchuan cattle, comprising body weight (BW), body height (BOH), back height (BAH), buttock height (BUH), chest depth (CD), chest width (CW), hip cross height (HCH), body length (BL), hip width (HW), rump length (RL), pin bone width (PBW), chest girth (CG), abdomen circumference (AG), and calf circumference (CC). After quality control, 281,889 SNPs were generated for GWAS with different traits. A total of 250 suggestive SNPs (p < 3.54 × 10-6) were screened and 37 candidate genes were annotated. Furthermore, we performed a linkage disequilibrium analysis of SNP loci and considered published studies, identifying the eight genes (ADAMTS17, ALDH1A3, CHSY1, MAGEL2, MEF2A, SYNM, CNTNAP5, and CTNNA3) most likely to be involved in growth traits. This study provides new insights into the regulatory mechanisms of bovine body size development, which can be very useful in the development of management and breeding strategies.

In the conservation of livestock and poultry resources, population genetic diversity and genetic structure of the conservation population are important factors affecting the effectiveness of conservation. In this study, whole-genome resequencing technology was used to detect genomic variation in 100 Qinchuan cattle (30 bulls and 70 cows). By analyzing population genetic diversity, runs of homozygosity (ROH) distribution features, kinship relationships, and family structure, the conservation effectiveness of Qinchuan cattle was comprehensively evaluated. The results showed that a total of 20,968,017 high-quality SNPs were detected in 100 Qinchuan cattle, the average minimum allele frequency was 0.191±0.124, the average polymorphic information content was 0.279±0.131, and the average observed heterozygosity was 0.275±0.131, the average expected heterozygosity is 0.279±0.131, indicating that the genetic diversity of the Qinchuan cattle conservation population is relatively rich. The average identity by state (IBS) distance of the Qinchuan conservation population was 0.243±0.020, with a value of 0.242±0.021 for the bulls. The results of the kinship G-matrix were consistent with the results of the IBS distance matrix, both showing that some individuals in the conservation population had close kinship. A total of 8258 genomic ROH were detected in 100 Qinchuan cattle, with a total length of 9.64 GB. The average length of ROH fragments was 1.167±1.203 Mb, 69.35% of the ROH were short ROH with a length of 0.5~1 Mb, and the average total length of ROH per individual was 96.40 Mb. The average inbreeding coefficient based on ROH was 0.039±0.039, with a value of 0.044±0.035 for the bulls, indicating that some bulls had a certain degree of inbreeding accumulation. The results of the phylogenetic tree combined with kinship analysis showed that the individuals in the Qinchuan cattle conservation population could be divided into eight families, including seven families with bulls and one family without bulls. This study demonstrated that the genetic diversity of the Qinchuan conservation population is relatively rich, with no significant inbreeding accumulation, but there is a risk of inbreeding among some individuals. Therefore, it is necessary to strengthen selection and mating to ensure the sustainable development of Qinchuan cattle resources.
在畜禽资源保护中，群体遗传多样性和遗传结构是决定保种效果的重要因素。本研究采用全基因组重测序技术检测100头秦川牛(30头公牛、70头母牛)的基因组变异，通过分析群体遗传多样性、连续纯合片段(runs of homozygosity，ROH)分布特征、亲缘关系和家系结构，对秦川牛的保种效果进行了综合评估。结果显示，100头秦川牛共检测到20,968,017个高质量SNPs位点，平均最小等位基因频率为0.191±0.124，平均多态信息含量为0.279±0.131，平均观察杂合度为0.275±0.131，平均期望杂合度为0.279±0.131，表明秦川保种群遗传多样性较为丰富。秦川牛保种群体平均状态同源(identity by state，IBS)遗传距离为0.243±0.020，其中公牛为0.242±0.021，亲缘关系G矩阵结果与IBS距离矩阵结果一致，均显示秦川牛保种群部分个体间亲缘关系较近。100头秦川牛个体共检测到8258个基因组ROH，ROH总长度为9.64 GB，平均ROH长度为1.167±1.203 Mb，69.35%的ROH是长度为0.5~1 Mb的短ROH。个体平均ROH总长度为96.40 Mb。基于ROH的平均近交系数为0.039±0.039，其中30头秦川公牛的平均近交系数为0.044±0.035，表明部分公牛个体存在一定程度的近交积累。进化树结果显示，秦川牛保种群所测个体可分为8个家系，包括7个含公牛家系和1个不含公牛家系。本研究表明，秦川牛保种群的遗传多样性较为丰富，未出现较大程度近交积累，但部分个体间存在近交风险，应强化选配以确保秦川牛资源的可持续发展。.

This study aimed to understand the influence of the non-genetic factors that include breeding year, season, and sex of growth and development traits of Qinchuan cattle and to estimate the heritability of body weight at different growth stages. The Qinchuan cattle measurement records were by the Experiment farm of the National Beef Cattle Improvement Center (Yangling, China) from 2000 to 2017. SPSS and R software were used to analyze the influence of non-genetic factors on body size traits that include body weight (BW), withers height (WH), hip height (HH), body length (BL), chest circumference (CC), abdominal girth (AG), and calf girth (CG), at birth, 6, 12, 18, and 24 months of age. Meanwhile, the single-trait animal model of DMU software was used to estimate the variance component and the heritability of body weight. The results of GLM analysis showed as follows: sex, birth year, and birth season had effects on the body size traits of Qinchuan cattle at different growth stages. Respectively, the heritability of body weight at birth, 6, 12, 18, and 24 months of age were 0.43, 0.32, 0.37, 0.32, and 0.38.

Qinchuan cattle has gradually improved in body shape and growth rate in the long-term breeding process from the draft cattle to beef cattle. As the head of the five local yellow cattle in China, the Qinchuan cattle has been designated as a specialized beef cattle breed. We investigated the selection signatures using whole genome sequencing data in Qinchuan cattle. Based on Fst, we detected hundreds of candidate genes under selection across Qinchuan, Red Angus, and Japanese Black cattle. Through protein-protein interaction analysis and functional annotation of candidate genes, the results revealed that KMT2E, LTBP1 and NIPBL were related to brain size, body characteristics, and limb development, respectively, suggesting that these potential genes may affect the growth and development traits in Qinchuan cattle. ARIH2, DACT1 and DNM2, et al. are related to meat quality. Meanwhile, TBXA2R can be used as a gene associated with reproductive function, and USH2A affect coat color. This provided a glimpse into the formation of breeds and molecular genetic breeding. Our findings will promote genome-assisted breeding to improve animal production and health.

The SMAD1 and SMAD5 genes belong to mothers against decapentaplegic proteins family, which participate in the BMP pathway to control skeletal myogenesis and growth. In the present study, we analyzed the associations between polymorphisms of SMAD1 and SMAD5 genes promoter and important economical traits in Qinchuan cattle. Four SNPs in the SMAD1 gene promoter and three SNPs in the SMAD5 promoter were identified by sequencing of 448 Qinchuan cattles. Allelic and frequency analyses of these SNPs resulted in eight haplotypes both in the promoters of the two genes promoter and identified potential cis-regulatory transcription factor (TF) components. In addition, correlation analysis showed that cattle SMAD1 promoter activity of individuals with Hap4 (P < 0.01) was stronger than that of individuals with Hap2. while the transcriptional activity of individuals with Hap3 within SMAD5 gene promoter was significantly (P < 0.01) higher followed by H2. Uniformly, diplotypes H4-H6 of SMAD1 gene and H1-H3 of SMAD5 gene performed significant (P < 0.01) associations with body measurement and improved carcass quality traits. All these results have indicated that polymorphisms in SMAD1 and SMAD5 genes promoter could impact the transcriptional regulation and then affect muscle content in beef cattle. Moreover, both the SMAD1 and SMAD5 genes were expressed ubiquitously in 10 tissues and had higher expression in the longissimus thoracis tissue from 6-month-old and 12-month-old cattle than in cattle of other ages. We can conclude that SMAD1 and SMAD5 genes may play an important role in muscle growth and development, and the variants mapped within SMAD1 and SMAD5 genes can be utilized in molecular marker-assisted selection for cattle carcass quality and body measurement traits in breed improvement programs of Qinchuan cattle.

The PLIN1 gene produces a phosphorylated protein wrapped in lipid droplets in adipocytes. This phosphorylation assists the mobilization of fat into adipose tissue. The purpose of the experiment was to study the polymorphism of the PLIN1 gene and its relationship with the body and carcass characteristics of Qinchuan cattle to find molecular genetic markers that can be used for breeding. The expression level of the PLIN1 gene was determined in various tissues by qRT-PCR. The results showed that the highest level of PLN1 expression was found in subcutaneous fat, followed by the heart and longissimus muscle, and the lowest level was found in the kidney. Five SNP loci of the PLIN1 gene were identified in 510 Qinchuan cattle, including g.3580T>C (SNP1), g.3898G>A (SNP2), g.8333G>A (SNP3), g.10517T>C (SNP4), and g.10538G>T (SNP5). The results show that SNP1, SNP2, SNP3, and SNP4 were moderately polymorphic (0.25 < PIC < 0.5), while SNP5 was minimally polymorphic (PIC < 0.25). SNP2, SNP3, and SNP5 were within Hardy-Weinberg equilibrium (P > 0.05), but SNP1 and SNP4 were not (P < 0.05). Correlation analysis showed that the five SNPs of the PLIN1 gene were correlated with back-fat depth, intramuscular fat, and chest depth of Qinchuan cattle. The double haplotype H2H4 in Qinchuan beef was associated with body and carcass traits. We conclude that variants mapped within PLIN1 can be used in marker-assisted selection for carcass quality and body traits in breed improvement programs for Qinchuan cattle.

In our study, four single nucleotide polymorphisms (SNPs) were identified in exon 2 of cofilin-1 (CFL1) gene in 488 Chinese Qinchuan (QC) cattle, which included two missense mutations T 2084G and G 2107C, two synonymous mutations T 2052C and T 2169C. Further, we evaluated haplotype frequency and linkage disequilibrium (LD) coefficient of four SNPs. At SNP T 2052C, G 2107C and T 2169C, the QC cattle population belonged to intermediate genetic diversity (0.25 < PIC-value < 0.5), whereas SNP T-2084G belonged to low polymorphism (PIC-value < 0.25). Haplotype analysis showed that 6 different haplotypes (frequency > 0.03). LD analysis showed that SNP G 2107C and T 2169C, SNP G 2107C and T 2084G were high LD, respectively (r2 > 0.33). Association analysis indicated that SNP T 2052C was significantly associated with body length, chest breadth, chest depth and body mass in the QC population (p < 0.01 or p < 0.05). SNP G 2107C was significantly associated with rump length (p < 0.05). SNP T 2169C was significantly associated with chest breadth and chest depth (p < .01 or p < .05). The results of our study suggest that the CFL1 gene may be a strong candidate gene that affects growth traits in the QC cattle breeding program.

Qinchuan cattle is one of the five yellow cattle breeds in China with good performance of meat. The proliferation and differentiation level of muscle and fat are closely related to the growth and development of the organism and are the key factors affecting the quality of meat. In order to study the effect of lncRNA on the fat tissues of Qinchuan cattle, six calf and adult bovine adipose tissues were selected for high-throughput sequencing. We obtained 3,716 lncRNA candidates from calves and adult cattle fat samples, among them 789 lncRNA were annotated and 2,927 lncRNA were novel lncRNA. A number of lncRNAs were highly abundant, and 119 lncRNA were differentially expressed between two developmental stages. We further validated several differentially expressed lncRNAs using qPCR, and the results were consistent with the sequencing data. Therefore, we conclude that lncRNA may play an important role in adipose tissue in different age groups of cattle.