This study aimed to identify the genes related to the body size of pigs by conducting genome-wide selection analysis (GWSA). We performed a GWSA scan on 50 pigs belonging to four small-bodied pig populations (Diannan small-eared pig, Bama Xiang pig, Wuzhishan pig, and Jeju black pig from South Korea) and 124 large-bodied pigs. We used the genetic parameters of the pairwise fixation index (F_ST) and π ratio (case/control) to screen candidate genome regions and genes related to body size. The results revealed 47,339,509 high-quality SNPs obtained from 174 individuals, while 280 interacting candidate regions were obtained from the top 1% signal windows of both parameters, along with 187 genes (e.g., ADCK4, AMDHD2, ASPN, ASS1, and ATP6V0C). The results of the candidate gene (CG) annotation showed that a series of CGs (e.g., MSTN, LTBP4, PDPK1, PKMYT1, ASS1, and STAT6) was enriched into the gene ontology terms. Moreover, molecular pathways, such as the PI3K-Akt, HIF-1, and AMPK signaling pathways, were verified to be related to body development. Overall, we identified a series of key genes that may be closely related to the body size of pigs, further elucidating the heredity basis of body shape determination in pigs and providing a theoretical reference for molecular breeding.

Boer goats, as kemp in meat-type goats, are selected and bred from African indigenous goats under a long period of artificial selection. Their advantages in multiple economic traits, particularly their plump growth, have attracted worldwide attention. The current study displayed the genome-wide selection signature analyses of South African indigenous goat (AF), African Boer (BH), and Australian Boer (AS) to investigate the hereditary basis of artificial selection in different stages. Four methods (principal component analysis, nucleotide diversity, linkage disequilibrium decay, and neighbor-joining tree) implied the genomic diversity changes with different artificial selection intensities in Boer goats. In addition, the θπ, F_ST, and XP-CLR methods were used to search for the candidate signatures of positive selection in Boer goats. Consequently, 339 (BH vs. AF) and 295 (AS vs. BH) candidate genes were obtained from SNP data. Especially, 10 genes (e.g., BMPR1B, DNER, ITGAL, and KIT) under selection in both groups were identified. Functional annotation analysis revealed that these genes are potentially responsible for reproduction, metabolism, growth, and development. This study used genome-wide sequencing data to identify inheritance by artificial selection. The results of the current study are valuable for future molecular-assisted breeding and genetic improvement of goats.

The follicle development (FD) is an important factor determining litter size in animals. Recent studies have found that noncoding RNAs (ncRNAs) play an important role in FD. In particular, the role of the regulatory mechanism of competing endogenous RNAs (ceRNAs) that drive FD has attracted increasing attention. Therefore, this study explored the genetic basis of goat FD by obtaining the complete follicular transcriptome of Dazu black goats at different developmental stages. Results revealed that 128 messenger RNAs (mRNAs), 4 long noncoding RNAs (lncRNAs), 49 microRNAs (miRNAs), and 290 circular RNAs (circRNAs) were significantly differentially expressed (DE) between large and small follicles. Moreover, DEmRNAs were enriched in many signaling pathways related to FD, as well as GO terms related to molecular binding and enzyme activity. Based on the analysis of the ceRNA network (CRN), 34 nodes (1 DElncRNAs, 10 DEcircRNAs, 14 DEmiRNAs, and 9 DEmRNAs) and 35 interactions (17 DEcircRNAs–DEmRNAs, 2 DElncRNAs–DEmiRNAs, and 16 DEmRNA–DEmiRNAs) implied that the CRN could be involved in the FD of goats. In conclusion, we described gene regulation by DERNAs and lncRNA/circRNA–miRNA–mRNA CRNs in the FD of goats. This study provided insights into the genetic basis of FD in precise transcriptional regulation.

Horns are an important breeding trait for sheep. However, no widely recognized viewpoint on the regulatory genes and mechanisms of horns is available, and the genetic basis of the four-horn phenotype (FHP) is unclear. This work conducted a genome-wide association study with 100 sheep genomes from multiple breeds to investigate the genetic basis of the FHP. The results revealed three significant associations (corrected as p < 1.64 × 10⁻⁸) of the InDels (CHR2: g.133,742,709delA, g.133,743,215insC, and g.133,743,940delT) for FHP in the intergenic sequence (IGS) between the MTX2 and the LOC105609047 of CHR2. Moreover, 14 significant associations (corrected as p < 1.42 × 10⁻⁹) of SNPs with the FHP phenotype were identified in CHR2 and CHR16, including five (e.g., CHR16: g.40,351,378G > A and g.40,352,577G > A) located in the intron of the ADAMTS12 gene, eight (e.g., CHR2: g.133,727,513C > T and g.133,732,145T > G) in the IGS between MTX2 and LOC105609047, and only one (CHR2: g.133,930,761A > G) in the IGS between HOXD1 and MTX2. Obvious divergence was also observed in genotype patterns between the FHP and others (two horns and hornless) in the HOXD1 and ADAMTS12 gene regions. An extremely significant linkage also occurred between Loci I and Loci II within 100 individuals (LD = −156.02186, p < 0.00001). In summary, our study indicated that the genomic sequences from CHR2 and CHR16 contributed to the FHP in sheep, specifically the key candidate genes HOXD1 and ADAMTS12. These results improved our understanding of the Mendelian genetic basis of the FHP in sheep.

Although coat color is an important economic phenotype in domesticated yaks (Bos grunniens), its genetic basis is not yet fully understood. In this study, a genome-wide selective sweep and high-frequency runs of homozygosity (ROH) identification were performed on 50 yaks with different coat colors to investigate candidate genes (CDGs) related to coat color. The results suggested that 2263 CDGs were identified from the 5% interaction windows of the F_ST and θπ ratio, along with 2801 and 2834 CDGs from black and brown yaks with iHS, respectively. Furthermore, 648 and 691 CDGs from black and brown yaks, which were widely enriched in pathways related to melanogenesis, melanocyte differentiation, and melanosome organization via GO and KEGG functional enrichment, respectively, were confirmed on the basis of the intersection of three parameters. Additionally, the genome of brown yaks presented more ROH, longer ROH fragments, and higher inbreeding levels than those of black yaks. Specifically, a large number of genes related to melanin synthesis and regulation (e.g., UST, TCF25, and AHRR) from the ROH islands were confirmed to be under strong selection. In summary, the results of this study enhance the understanding of the genetic basis for determining yak coat color.

Supernumerary teats (SNTs) or nipples often emerge around the mammary line. This study performed a genome-wide selective sweep analysis (GWS) at the copy number variant (CNV) level using two selected signal calculation methods (VST and FST) to identify candidate genes associated with SNTs in goats. A total of 12,310 CNVs were identified from 37 animals and 123 CNVs, with the top 1% VST values including 84 candidate genes (CDGs). Of these CDGs, minichromosome maintenance complex component 3, ectodysplasin A receptor associated via death domain, and cullin 5 demonstrated functions closely related to mammary gland development. In addition, 123 CNVs with the top 1% FST values were annotated to 97 CDGs. 5-Hydroxytryptamine receptor 2A, CCAAT/enhancer-binding protein alpha, and the polymeric immunoglobulin receptor affect colostrum secretion through multiple signaling pathways. Two genes, namely, RNA-binding motif protein 46 and β-1,3-galactosyltransferase 5, showed a close relation to mammary gland development. Six CNVs were identified and annotated to five genes by intersecting the top 1% of candidate CNVs with both parameters. These genes include LOC102185621, LOC102190481, and UDP-glucose pyrophosphorylase 2, which potentially affect the occurrence of BC through multiple biological processes, such as cell detoxification, glycogen synthesis, and phospholipid metabolism. In conclusion, we discovered numerous genes related to mammary development and breast cancer (BC) through a GWS, which suggests the mechanism of SNTs in goats and a certain association between mammary cancer and SNTs.