Whole-genome resequencing-based comparative variant analysis identifies candidate genes associated with cross-beak phenotype in Huiyang Bearded chickens

Cross-beaks are deemed a threat to poultry health, productivity, and animal welfare. Nevertheless, due to sporadic cases, heterogeneity of gene loci and incomplete dominance, the molecular mechanism of cross-beak formation, especially the degree of cross, is not yet clear. Thus, we screen key genes and reveal the possible phenotypic formation mechanism of cross-beak by comparison with different degrees of deformity in Huiyang Bearded chickens by compare whole-genome resequencing-based variant analysis. Comparative analysis between cross-beak and normal-beaked chickens identified differential variants in several candidate genes, including CDH11, CTNNAL1, NRXN3, NRXN1, CDH5, SDC3, and DHFR. Genes harboring these variants were enriched in pathways related to cell adhesion molecules and metabolic processes, with functional annotations involving cell-cell adhesion and neural crest cell migration. Comparative analysis between chickens with severe and slight cross-beak deformities identified additional candidate genes, including MRPL21, NSUN2, DDX55, GNB3, and NFKB2. These genes were associated with enriched terms and pathways related to focal adhesion, amyotrophic lateral sclerosis, steroid 7 -hydroxylase activity, and skin-barrier establishment. These findings provide a preliminary catalogue of genetic variants and candidate genes for future functional studies of cross-beak development and severity in chickens.