Supporting data for "Identifying candidate genetic variants for egg number by analyzing over 1000 fully sequenced layers"

Egg production over a long laying cycle until 700 days of age is fancied for modern layer chickens breeding. It is influenced by the onset of laying, stability during the peak period, and persistence at late laying stages. Conventional single-single nucleotide polymorphisms (SNP) association analyses have identified additive loci, but few studies have explored dominance effects or integrated multi-omics data to investigate the genetic basis of egg production traits from the onset to 700 days of age. A full diallel cross of 1,004 chickens was·subjected to whole-genome sequencing. Transcriptome data from the ovary was available for a subset of 120 chickens. A genome-wide association study (GWAS) was conducted using an additive-dominance model for cumulative egg number and egg number at different stages. Expression quantitative trait loci (eQTL) mapping was applied to investigate associations between SNPs and gene expression. A transcriptome-wide association study (TWAS) was conducted to explore the associations between gene expression and egg production traits to identify candidate genes. <br>The additive-dominance model identified 5,892 significant SNPs, comprising 805 additive SNPs and 360 dominance SNPs shared between two or more traits. By integrating loci identified through GWAS with eQTL-mapping, the expression level of 27 genes were found to be associated with significant SNPs. Further integration with TWAS results revealed four novel candidate genes. For the loci with significant SNP effects, we found a positive but insignificant correlation between the ratios of dominance to additive effects and observed heterosis. Observed heterosis was positively correlated with heterosis predicted based on dominance effects and allele frequencies of all SNPs. <br>We identified candidate genetic variants for egg production traits by analyzing 1,004 fully sequenced layers. Detection benefited from incorporating dominance into the GWAS model. Traits with higher heterosis tended to be more affected by genes with dominant mode of action. Moreover, multi-omics data allowed to contribute to deciphering genetic mechanisms underlying egg production by establishing connections between genetic variants, gene expression, and egg number.