Identification of Goat Supernumerary Teat Phenotype Using Wide-Genomic Copy Number Variants

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 (<i>V</i>_ST and <i>F</i>_ST) 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% <i>V</i>_ST 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% <i>F</i>_ST 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 <i>LOC102185621</i>, <i>LOC102190481</i>, 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.