Copy number variation (CNV) is a structural variant with significant impact on genetic diversity. CNV has been widely used in breeding for growth traits, meat production or quality, and coat color. SRY-like box genes (SOXs) are a class of transcription factors that play a regulatory role in cell fate specification and differentiation. SOX5 and SOX8 belong to subgroups D and E of the SOXs, respectively. Previous studies have shown that SOX5 and SOX8 are essential in the development of bones. In this study, we explored the association between the growth traits and CNVs of SOX5 and SOX8 in 326 Ashidan yaks and detected mRNA expression levels in different tissues. Our results illustrated that CNVs of SOX5 and SOX8 were significantly associated with withers height at 18 months of age and chest girth at 30 months of age (p < 0.05). The CNV combination of SOX5 and SOX8 was significantly associated with withers height at 18 months of age (p < 0.01). SOX5 expression in the lung was significantly higher than in the heart, spleen, kidney, and muscle (p < 0.05). SOX8 expression in the lung was significantly higher than in the liver and muscle (p < 0.05). Our results provide evidence that the CNVs of SOX5 and SOX8 genes could be used as new markers for the selection of yak growth traits.

Long-haired individuals in the Tianzhu white yak population are a unique genetic resource, and have important landscape value. Copy number variation (CNV) is an important source of phenotypic variation in mammals. In this study, we used resequencing technology to detect the whole genome of 10 long-haired Tianzhu white yaks (LTWY) and 10 normal-haired Tianzhu white yaks (NTWY), and analyzed the differences of CNV in the genome of LTWYs and NTWYs. A total of 110268 CNVs were identified, 2006 CNVRs were defined, and the distribution map of these CNVRs on chromosomes was constructed. The comparison of LTWYs and NTWYs identified 80 differential CNVR-harbored genes, which were enriched in lipid metabolism, cell migration and other functions. Notably, some differential genes were identified as associated with hair growth and hair-follicle development (e.g., ASTN2, ATM, COL22A1, GK5, SLIT3, PM20D1, and SGCZ). In general, we present the first genome-wide analysis of CNV in LTWYs and NTWYs. Our results can provide new insights into the phenotypic variation of different hair lengths in Tianzhu white yaks.

Understanding animal growth plays an important role in improving animal genetics and breeding. In order to explore the early growth and development law of Ashidan yak, the body weight (BW), wither height (WH), body oblique length (BL) and chest girth (CG) of 260 female Ashidan yaks were measured. These individuals grew under grazing conditions, and growth traits were measured at 6, 12, 18 and 30 months of age. Then the absolute growth and relative growth of Ashidan yak were calculated, and five nonlinear models (Logistic model, Gompertz model, Brody model, von Bertalanffy model and Richards model) were used to fit the growth curve of Ashidan yak. The fitting effect of the model was evaluated according to MSE, AIC and BIC. The results showed that the growth rate of Ashidan yak was the fastest from 12 to 18 months old, and the growth was slow or even stagnant from 6 to 12 months old. The AIC and BIC values of the Richards model were the lowest among the five models, with an AIC value of 4543.98 and a BIC value of 4563.19. The Richards model estimated body weight at 155.642 kg. In summary, the growth rate of female Ashidan yak changes with the seasons, growing faster in warm seasons and slower in cold seasons. Richards model is the best model to describe the growth curve of female Ashidan yak in five nonlinear models.

Although significant advancement has been made in the study of adipogenesis, knowledge about how chromatin accessibility regulates yak adipogenesis is lacking. We here described genome-wide dynamic chromatin accessibility in preadipocytes and adipocytes by using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), and thus revealed the unique characteristics of open chromatin during yak adipocyte differentiation. The chromatin accessibility of preadipocytes and adipocytes exhibited a similar genomic distribution, displaying a preferential location within the intergenic region, intron, and promoter. The pathway enrichment analysis identified that genes with differential chromatin accessibility were involved in adipogenic metabolism regulation pathways, such as the peroxisome proliferator activated receptor-γ (PPAR) signaling pathway, wingless-type MMTV integration site (Wnt) signaling pathway, and extracellular matrix-receptor (ECM–receptor) interaction. Integration of ATAC-seq and mRNA-seq revealed that genes with a high expression were associated with high levels of chromatin accessibility, especially within 1 kb upstream and downstream of the transcription start site. In addition, we identified a series of transcription factors (TFs) related to adipogenesis and created the TF regulatory network, providing the possible interactions between TFs during yak adipogenesis. This study is crucial for advancing the understanding of transcriptional regulatory mechanisms of adipogenesis and provides valuable information for understanding the adaptation of plateau species to high-altitude environments by maintaining whole body homeostasis through fat metabolism.

The growth hormone receptor (GHR) is a member of the cytokine/hematopoietic factor receptor superfamily, which plays an important role in the growth and development, immunity, and metabolism of animals. This study identified a 246 bp deletion variant in the intronic region of the GHR gene, and three genotypes, including type II, type ID, and type DD, were observed. Genotype analysis of structural variation (SV) was performed on 585 individuals from 14 yak breeds, and it was found that 246 bp deletion was present in each breed. The II genotype was dominant in all yak breeds except for SB yak. The association analysis of gene polymorphisms and growth traits in the ASD yak population showed that the 246 bp SV was significantly associated with body length at 6 months (p < 0.05). GHR messenger RNA (mRNA) was expressed in all the tested tissues, with significantly higher levels in the liver, muscle, and fat than in other organs. The results of transcription activity showed that the luciferase activity of the pGL4.10-DD vector was significantly higher than that of the pGL4.10-II vector (p < 0.05). Additionally, the transcription-factor binding prediction results showed that the SV in the runt-related transcription factor 1 (Runx1) transcription-factor binding site may affect the transcriptional activity of the GHR gene, regulating yak growth and development. This study showed that the novel SV of the GHR gene could be used as a candidate molecular marker for the selection of the early growth trait in ASD yak.

In the Suidae family, warthogs show significant survival adaptability and trait specificity. This study offers a comparative genomic analysis between the warthog and other Suidae species, including the Luchuan pig, Duroc pig, and Red River hog. By integrating the four genomes with sequences from the other four species, we identified 8868 single-copy orthologous genes. Based on 8868 orthologous protein sequences, phylogenetic assessments highlighted divergence timelines and unique evolutionary branches within suid species. Warthogs exist on different evolutionary branches compared to DRCs and LCs, with a divergence time preceding that of DRC and LC. Contraction and expansion analyses of warthog gene families have been conducted to elucidate the mechanisms of their evolutionary adaptations. Using GO, KEGG, and MGI databases, warthogs showed a preference for expansion in sensory genes and contraction in metabolic genes, underscoring phenotypic diversity and adaptive evolution direction. Associating genes with the QTLdb-pigSS11 database revealed links between gene families and immunity traits. The overlap of olfactory genes in immune-related QTL regions highlighted their importance in evolutionary adaptations. This work highlights the unique evolutionary strategies and adaptive mechanisms of warthogs, guiding future research into the distinct adaptability and disease resistance in pigs, particularly focusing on traits such as resistance to African Swine Fever Virus.