Whole-exome sequencing identified rare variants associated with body length and girth in cattle

Body measurements can be used in determining body size to monitor the cattle growth and examine the response to selection. Despite efforts putting into the identification of common genetic variants, the mechanism understanding of the rare variation in complex traits about body size and growth remains limited. Here, we firstly performed GWAS study for body measurement traits in Simmental cattle, however there were no SNPs exceeding significant level associated with body measurements. To further investigate the mechanism of growth traits in beef cattle, we conducted whole exome analysis of 20 cattle with phenotypic differences on body girth and length, representing the first systematic exploration of rare variants on body measurements in cattle. By carrying out a three-phase process of the variant calling and filtering, a sum of 1158, 1151, 1267, and 1303 rare variants were identified in four phenotypic groups of two growth traits, higher/ lower body girth (BG_H and BG_L) and higher/lower body length (BL_H and BL_L) respectively. The subsequent functional enrichment analysis revealed that these rare variants distributed in 886 genes associated with collagen formation and organelle organization, indicating the importance of collagen formation and organelle organization for body size growth in cattle. The integrative network construction distinguished 62 and 66 genes with different co-expression patterns associated with higher and lower phenotypic groups of body measurements respectively, and the two sub-networks were distinct. Gene ontology and pathway annotation further showed that all shared genes in phenotypic differences participate in many biological processes related to the growth and development of the organism. Together, these findings provide a deep insight into rare genetic variants of growth traits in cattle and this will have a promising application in animal breeding.