Search

Organismic evolution requires that variation at distinct hierarchical levels and attributes be coherently integrated, often in the face of disparate environmental and genetic pressures. A central part of the evolutionary analysis of biological systems remains to decipher the causal connections between organism-wide (or genome-wide) attributes (e.g., mRNA abundance, protein length, codon bias, recombination rate, genomic position, mutation rate, etc) as well as their role-together with mutation, selection, and genetic drift-in shaping patterns of evolutionary variation in any of the attributes themselves. Here we combine genome-wide evolutionary analysis of protein and gene expression data to highlight fundamental relationships among genomic attributes and their associations with the evolution of both protein sequences and gene expression levels. Our results show that protein divergence is positively coupled with both gene expression polymorphism and divergence. We show moreover that although the number of protein-protein interactions in Drosophila is negatively associated with protein divergence as well as gene expression polymorphism and divergence, protein-protein interactions cannot account for the observed coupling between regulatory and structural evolution. Furthermore, we show that proteins with higher rates of amino acid substitutions tend to have larger sizes and tend to be expressed at lower mRNA abundances, whereas genes with higher levels of gene expression divergence and polymorphism tend to have shorter sizes and tend to be expressed at higher mRNA abundances. Finally, we show that protein length is negatively associated with both number of protein-protein interactions and mRNA abundance and that interacting proteins in Drosophila show similar amounts of divergence. We suggest that protein sequences and gene expression are subjected to similar evolutionary dynamics, possibly because of similarity in the fitness effect (i.e., strength of stabilizing selection) of disruptions in a gene's protein sequence or its mRNA expression. We conclude that, as more and better data accumulate, understanding the causal connections among biological traits and how they are integrated over time to constrain or promote structural and regulatory evolution may finally become possible.

We report multiple cases in which disruption of hsp70 regulatory regions by transposable element (TE) insertions underlies natural variation in expression of the stress-inducible molecular chaperone Hsp70 in Drosophila melanogaster. Three D. melanogaster populations from different continents are polymorphic for jockey or P element insertions in the promoter of the hsp70Ba gene. All three TE insertions are within the same 87-bp region of hsp70Ba promoter, and we suggest that the distinctive promoter architecture of hsp genes may make them vulnerable to TE insertions. Each of the TE insertions reduces Hsp70 levels, and RNase protection assays demonstrate that such insertions can reduce transcription of the hsp70Ba gene. In addition, the TEs alter two measures of organismal fitness, inducible thermotolerance and female reproductive success. Thus, transposition can create quantitative genetic variation in gene expression within populations, on which natural selection can act.

To determine how the modern copy number (5) of hsp70 genes in Drosophila melanogaster evolved, we localized the duplication events that created the genes in the phylogeny of the melanogaster group, examined D. melanogaster genomic sequence to investigate the mechanisms of duplication, and analyzed the hsp70 gene sequences of Drosophila orena and Drosophila mauritiana. The initial two-to-four hsp70 duplication occurred 10--15 MYA, according to fixed in situ hybridization to polytene chromosomes, before the origin and divergence of the melanogaster and five other species subgroups of the melanogaster group. Analysis of more than 30 kb of flanking sequence surrounding the hsp70 gene clusters suggested that this duplication was likely a retrotransposition. For the melanogaster subgroup, Southern hybridization and an hsp70 restriction map confirmed the conserved number (4) and arrangement of hsp70 genes in the seven species other than D. melanogaster. Drosophila melanogaster is unique; tandem duplication and gene conversion at the derived cluster yielded a fifth hsp70 gene. The four D. orena hsp70 genes are highly similar and concertedly evolving. In contrast, the D. mauritiana hsp70 genes are divergent, and many alleles are nonfunctional. The proliferation, concerted evolution, and maintenance of functionality in the D. melanogaster hsp70 genes is consistent with the action of natural selection in this species.