Emergence of Complex Haplotypes from Microevolutionary Variation in Sequence and Structure of <Emphasis Type="Italic">Colias</Emphasis> Phosphoglucose Isomerase

A molecular evolutionary explanation of natural genetic variation requires analysis of specific variants’ evolutionary dynamics. To pursue this for phosphoglucose isomerase (PGI) of Colias butterflies, whose polymorphism is maintained by strong natural selection, we assembled a large data set of wild haplotypes, highly variable at the amino acid and DNA levels. The most common electrophoretic, i.e., charge macrostate, allele class, 3, is conserved in its pattern of charged amino acid residues. The next most common macrostate, 4, has multiple patterns of charge, i.e., microstates, while less common (1, 2, 5, 6) macrostates are very diverse. Macrostate 4 shows significant linkage disequilibrium (LD) among its variants, especially for two groups of five haplotypes each. We find extensive intragenic recombination among all haplotypes except the two high-LD groups of macrostate 4, which display none. Phyletic relations among haplotypes are largely reticulate, again except for the high-LD groups of macrostate 4, which form clades with strong bootstrap support. Charge-changing and linked charge-neutral amino acid variants occur in diverse parts of PGI’s sequence. Homology-based modeling of PGI’s structure shows that these regions are related spatially in ways suggesting functional interaction. The high-LD groups of macrostate 4 display parallel amino acid variation in these regions. This pattern of haplotype clades with high LD among multiple varying sites, emerging from chaotically recombining variation, may be a “signature” of refinement of complex adaptive sequences by recombination and selection. It should be tested further in this study system and others as a possibly general feature of the evolution of living complexity.