Understanding the population structure and genetic bases of well-adapted cattle breeds to local environments is one of the most essential tasks to develop appropriate genetic improvement programs. We performed a comprehensive study to investigate the population structure, divergence and selection signatures at genome-wide level in diverse Chinese local cattle using Bovine HD SNPs array, including two breeds from North China, one breed from Northwest China, three breeds from Southwest China and two breeds from South China. Population genetic analyses revealed the genetic structures of these populations were mostly related to the geographic locations. Notably, we detected 294 and 1263 candidate regions under selection using the di and iHS approaches, respectively. A series of group-specific and breed-specific candidate genes were identified, which are involved in immune response, sexual maturation, stature related, birth and bone weight, embryonic development, coat colors and adaptation. Furthermore, haplotype diversity and network pattern for candidate genes, including LPGAT1, LCORL, PPP1R8, RXFP2 and FANCA, suggest that these genes have been under differential selection pressures in various environmental conditions. Our results shed insights into diverse selection during breed formation in Chinese local cattle. These findings may promote the application of genome-assisted breeding for well-adapted local breeds with economic and ecological importance.

Average daily gain (ADG) is an important trait that contributes to the production efficiency and economic benefits in the beef cattle industry. The molecular mechanisms of ADG have not yet been fully explored because most recent association studies for ADG are based on SNPs or haplotypes. We reported a systematic CNV discovery and association analysis for ADG in Chinese Simmental beef cattle. Our study identified 4912 nonredundant CNVRs with a total length of ~ 248.7 Mb, corresponding to ~ 8.9% of the cattle genome. Using probe-based CNV association, we identified 24 and 12 significant SNP probes within five deletions and two duplications for ADG, respectively. Among them, we found one common deletion with 89 kb imbedded in LHFPL Tetraspan Subfamily Member 6 (LHFPL6) at 22.9 Mb on BTA12, which has high frequency (12.9%) dispersing across population. CNV selection test using VST statistic suggested this common deletion may be under positive selection in Chinese Simmental cattle. Moreover, this deletion was not overlapped with any candidate SNP for ADG compared with previous SNPs-based association studies, suggesting its important role for ADG. In addition, we identified one rare deletion near gene Growth Factor Receptor-bound Protein 10 (GRB10) at 5.1 Mb on BTA4 for ADG using both probe-based association and region-based approaches. Our results provided some valuable insights to elucidate the genetic basis of ADG in beef cattle, and these findings offer an alternative perspective to understand the genetic mechanism of complex traits in terms of copy number variations in farm animals.

Myogenesis is a complex process controlled by several coding and non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) that are known to function as endogenous microRNAs (miRNAs) sponges. Cerebellar Degeneration-Related protein 1 antisense (CDR1as) is the most spotlighted circRNA that is known as an miR-7 sponge, which has bloomed circRNAs’ research in animal disease and physiology. Here, we screened for miRNAs and mRNA associated with CDR1as and further characterized their regulatory function during muscle differentiation. We found that a total of 43 miRNAs (including miR-107-3p, miR-125b-5p, miR-140-5p, miR-29a-3p, and miR-27a-3p upregulated) and 789 mRNAs (including ANGPT1, E2F2, CCN1, FGFR1, and MEF2C downregulated) were differentially expressed in goat skeletal muscle satellite cells (SMSCs). Further, knockdown of CDR1as and ANGPT1 inhibited SMSCs differentiation. miR-27a-3p was differentially upregulated after the knockdown of CDR1as in SMSCs. Overexpressed miR-27a-3p decreased SMSCs differentiation. Via RNAhybrid and luciferase, miR-27a-3p was identified to regulate ANGPT1. We discovered that miR-27a-3p has an inverse relationship with CDR1as and decreases the expression level of ANGPT1 during SMSCs differentiation. In summary, our study demonstrates that siCDR1as inhibits myoblast differentiation by downregulating ANGPT1 mRNA via miR-27a-3p in SMSCs.

The use of next-generation probiotics (NGP) in pigs for combating diseases has been subject to limited research. Here we explored the potential of a well-known NGP candidate Akkermansia muciniphila targeting pig gut health. In the first screening experiment, we found that the abundance of A. muciniphila peaked at 14 d old but decreased at weaning (21 d old; P < 0.05), suggesting the weaning period may be an effective window for A. muciniphila intervention. Following that, 48 crossbred weaned pigs at 28 d old were randomly assigned to five groups: control (CON), high/low live A. muciniphila (HA/LA), and high/low heat-killed A. muciniphila (HIA/LIA). From 1 to 28 d old, the CON group received gastric infusion of anaerobic sterile saline every other day; the HA and LA groups were gavaged every other day with 1 × 10¹⁰ CFU/5 mL and 5 × 10⁸ CFU/5 mL live A. muciniphila, respectively; and the HIA and LIA groups were gavaged every other day with 1 × 10¹⁰ CFU/5 mL and 5 × 10⁸ CFU/5 mL heat-killed A. muciniphila, respectively. At d 29, pigs in the CON group were randomly and equally divided into two groups, one of which was named the enterotoxigenic Escherichia coli (ETEC) group, and all groups except CON received a 5-d ETEC challenge. The supplementation of A. muciniphila numerically reduced the diarrhea rate of weaned pigs compared to the pigs that only received the ETEC challenge (P = 0.57), but the LIA group had a higher diarrhea rate than the CON group (P < 0.05). Consistent with this, the supplementation of A. muciniphila improved the small intestinal morphology and structure, proportion of CD4⁺ T lymphocytes in the blood, as well as the expression of genes related to intestinal barrier and antioxidant indices of pigs with ETEC challenge, especially for the LA group (P < 0.05). Meanwhile, A. muciniphila supplementation reduced the expression of ETEC virulence factor genes in the ileum and colon of pigs challenged by ETEC (P < 0.05). Therefore, A. muciniphila may protect the intestinal health of weaned piglets from damage caused by ETEC infection, but the effect may vary depending on the concentration and activity of A. muciniphila.

Based on a pangenome graph platform, we simultaneously analyzed the impacts of SNPs and SVs in the population structure and phenotypic formation of global cattle using 2,409 individuals from 82 breeds. We demonstrated that SVs, like SNPs, effectively explain the population structure of global cattle. Genomic regions under strong selection, identified using both SNPs and SVs, consistently revealed footprints associated with human-mediated selection of economic traits in European improved cattle or natural selection of geographical adaptations. Notably, we detected that ∼40.14% of SVs were not tagged (LD, r2 < 0.6) by nearby SNPs. These “orphan” SVs may uncover new genetic signals and represent recent mutations associated with specific selection pressures or local environmental adaptation. Selected SVs tagged by SNPs also play causal or dominant roles in regions under selection. For example, our single-cell RNA sequencing has demonstrated that a notable SNP-tagged SV functions as an enhancer of the IGFBP7 gene, regulating fat deposition through IGFBP7+ cells. In conclusion, these SV-related mechanisms likely have caused some differences in economic traits and local adaptability across global cattle populations. Our integrated approaches highlight the unique and indispensable roles of SVs in shaping genetic diversity, offering novel insights into adaptation, selection, and strategies for improving cattle populations.