Abstract
Blue-green eggs exhibit unique shell color; however, compared to commercial layers, blue-green eggshell chickens have lower egg production and lack uniform shell colors. Aiming to confirm the molecular mechanisms that affect shell color and egg production, this study collected the uteruses of 12 blue-green eggshell chickens (BG group) and six Hy-Line layers (Brown group), which had significantly different shell color indexes (SCI) and egg numbers at 300 days of age (EN300). Transcriptome sequencing and comparative analyses were subsequently performed. BG hens were divided into two groups for comparative analysis (BGblue vs. BGgreen and BGlow vs. BGhigh, respectively), based on the differences in SCI and EN300, respectively. The result of weighted gene co-expression network (WGCNA) analysis showed that the sequenced and mapped 16,785 genes were clustered into 18 modules, among which six modules with a total of 4270 genes were highly correlated with SCI and EN300 traits. Five hundred and eleven differentially expressed genes (DEGs) belonged to the six key modules. Through KEGG mapping, GO enrichment, and Cytoscape network analysis, nine Hub genes were tightly associated with SCI and EN300. The up-regulated genes were CCR2, CCR8, CD40LG, IL18RAP, INHBA, and P2RY13, while the down-regulated genes were ABCA13, ADCY2, and GRM1. Co-analyses with the results of comparisons between the BG subgroups revealed that the expression of solute carrier (SLC) proteins and ABC transporters were highly related to eggshell color, while cytokine-cytokine receptor interactions and neuroactive ligand-receptor interactions were key pathways affecting egg production. The expression of extracellular cytokines and membrane receptors were significantly up-regulated in low-yield chickens. The candidate DEGs and pathways found in the study will assist in clarifying the molecular mechanisms affecting shell color and egg production, and improve the breeding of blue-green eggshell chickens.