Integrative ATAC-seq and RNA-seq analysis reveals lactation performance between Sewa sheep and East Friesian sheep

Lactation performance is a pivotal economic trait in sheep production, yet its underlying epigenetic regulatory mechanisms remain poorly understood. In the present study, we integrated ATAC-seq and RNA-seq to compare chromatin accessibility landscapes and transcriptomic in mammary gland tissues from Sewa sheep (SWS) and East Friesian sheep (EFS). Histological characterization revealed that SWS exhibited significantly smaller mammary acini area, smaller lipid droplet area, and reduced lipid droplet diameter compared to EFS. ATAC-seq analysis identified 15,902 differentially accessible regions (DARs) between the two breeds, with motif enrichment analysis uncovering key transcription factors potentially governing lactation traits. RNA-seq analysis revealed 1,163 differentially expressed genes (DEGs), which were involved in lactation regulation. Integrated analysis identified 441 overlapping genes, and enriched in glycolysis/gluconeogenesis (e.g., PGAM1, ENO1) and pyruvate metabolism (e.g., ACACA, ACSS1, ACYP1). Collectively, our study provides new insights into the epigenetic regulatory mechanisms underlying lactation performance differences in sheep.