Integrated GWAS and methylation analysis identify DNMT3A as an important regulator of growth in rabbits

The parameters of individual growth curve can serve as pseudo-phenotype for genetic evaluation in livestock. In this study, we compared five nonlinear growth models using post-weaning body weights of 706 New Zealand White rabbits. Under the best-fitting model, two parameters of mature weight and maturity rate were subjected to GWAS through single-step genomic BLUP framework that integrated phenotypic records from non-genotyped animals with 41,359 SNPs genotyped in 198 individuals. Association analysis identified 147 relevant genomic regions, and also highlighted DNMT3A as a promising candidate gene for further functional investigation. siRNA-mediated knockdown of DNMT3A significantly impaired myoblast proliferation. Whole-genome bisulfite sequencing of DNMT3A-knockdown myoblasts identified 69,480 differentially methylated regions (DMRs). Integrative analyses revealed substantial overlap between DMR-associated genes and GWAS candidate genes, with significant enrichment in vitamin B6 and tyrosine metabolism pathways. These findings suggest that DNMT3A may regulate rabbit growth via mediating DNA methylation of downstream genes.