Abstract
Variation in protein 3D structures reflects genetic variation and contributes to phenotypic diversity, yet its underlying genetic mechanisms remain unclear. To investigate the relationship between protein 3D structure and phenotype, we predict the 3D structures of 795,649 proteins from 26 maize (Zea mays L.) inbred lines using AlphaFold2. Population genetics analysis of these protein 3D structures reveal that buried residues held greater genomic evolutionary rate profiling (GERP) scores than exposed residues, indicating that buried residues are under stronger purifying selection. The design of the maize nested association mapping population makes it possible to utilize haplotype information and protein 3D structural variation to reveal the molecular mechanisms linking genetic diversity and phenotypic variation for a population with about 5000 individuals. Associating protein 3D structure variation with phenotypes (structure-based proteome-wide association study [PWAS]) identifies 14.2% more (96 vs. 84) significant proteins compared with associating protein sequence with phenotypes (sequence-based PWAS) using 32 agronomic traits. Moreover, structure-based PWAS identifies 24 additional significant proteins unique to predicted structures, whereas sequence-based PWAS identifies 12 additional significant proteins. Structure-based proteome-wide predictions (PWPs) improve genomic prediction accuracy by an average of 3.8% compared with sequence-based PWPs. In general, predicted protein 3D structures represent a powerful approach for understanding the natural diversity of protein haplotypes.