Sexually antagonistic cytonuclear fitness interactions in Drosophila melanogaster.

Theoretical and empirical studies have shown that selection cannot maintain a joint nuclear-cytoplasmic polymorphism within a population except under restrictive conditions of frequency-dependent or sex-specific selection. These conclusions are based on fitness interactions between a diploid autosomal locus and a haploid cytoplasmic locus. We develop a model of joint transmission of X chromosomes and cytoplasms and through simulation show that nuclear-cytoplasmic polymorphisms can be maintained by selection on X-cytoplasm interactions. We test aspects of the model with a "diallel" experiment analyzing fitness interactions between pairwise combinations of X chromosomes and cytoplasms from wild strains of Drosophila melanogaster. Contrary to earlier autosomal studies, significant fitness interactions between X chromosomes and cytoplasms are detected among strains from within populations. The experiment further demonstrates significant sex-by-genotype interactions for mtDNA haplotype, cytoplasms, and X chromosomes. These interactions are sexually antagonistic--i.e., the "good" cytoplasms in females are "bad" in males--analogous to crossing reaction norms. The presence or absence of Wolbachia did not alter the significance of the fitness effects involving X chromosomes and cytoplasms but tended to reduce the significance of mtDNA fitness effects. The negative fitness correlations between the sexes demonstrated in our empirical study are consistent with the conditions that maintain cytoplasmic polymorphism in simulations. Our results suggest that fitness interactions with the sex chromosomes may account for some proportion of cytoplasmic variation in natural populations. Sexually antagonistic selection or reciprocally matched fitness effects of nuclear-cytoplasmic genotypes may be important components of cytonuclear fitness variation and have implications for mitochondrial disease phenotypes that differ between the sexes.

Gene flow versus local adaptation in the northern acorn barnacle, Semibalanus balanoides: insights from mitochondrial DNA variation.

In reciprocal transplant experiments, Bertness and Gaines (1993) found that Semibalanus balanoides juveniles that had settled in an upper Narragansett Bay estuary survived better in that estuary that did juveniles from coastal localities. The observed pattern of survivorship led to the claim that local adaptation may result from a combination of limited gene flow between and strong selection within these habitats. Here we test the hypothesis that limited gene flow has led to habitat-specific population differentiation using sequence and restriction fragment length polymorphism analyses of the mitochondrial DNA D-loop region of S. balanoides. Samples were analyzed from replicated coastal and estuary localities in both Narragansett Bay, Rhode Island, and Damariscotta River, Maine. The patterns of F(ST) indicate that gene flow between coast and estuary is extensive (Nm > 100) and is not lower in the estuary with lower flushing rates (Narragansett Bay). Given the high estimate of genetic exchange, adaptations for unpredictable environments seem more likely than local adaptation in this species because loci that respond to selection in one generation are essentially homogenized by the next seasons' settlement. Nevertheless, these estimates of neutral gene flow can help identify the strength of selection necessary for local adaptation to accumulate in Semibalanus.

Mutation and selection at silent and replacement sites in the evolution of animal mitochondrial DNA.

Two patterns are presented that illustrate the interaction of mutation and selection in the evolution of animal mtDNA: 1) variation among taxa in the ratio of polymorphism to divergence (rpd) at silent and replacement sites in protein-coding genes, and 2) strand-differences in polymorphism and divergence at 'silent' sites that suggest a mutation-selection balance in the evolution of codon usage. Cytochrome b data from GenBank show that about half of the species pairs tested have a significant excess of amino acid polymorphism, relative to divergence. The remaining half of species pairs do not depart from neutrality, but generally do show an excess of amino acid polymorphism. Sequences from Drosophila pseudoobscura displaying a signature of an expanding population show a slight, but non-significant, deficiency of amino acid polymorphism suggestive of recently intensified selection on mildly deleterious mutations. Genes whose reading frames lie on the major coding strand of Drosophila mtDNA show a preponderance of T- > C substitutions, while genes encoded on the minor strand experience more A- > G than T- > C substitutions between species at both silent and replacement sites. However, silent mutations at third codon positions are introduced into the population in proportions opposite to those observed as fixed differences between species (e.g., an excess of T- > C polymorphisms are found at the ND5 gene on the minor coding strand). The high A + T content of insect mtDNAs imposes strong codon usage bias favoring A-ending and T-ending codons resulting in a distinct mutation-selection balance for genes encoded on opposites strands. Thus, at both replacement and silent sites, mutations that appear to be constrained in terms of divergence between species are in excess within species. The data suggest that mildly deleterious mutations are common in mitochondrial genes. A test of this, and a competing, hypothesis is proposed that requires additional sequence surveys of polymorphism and divergence. An important challenge is to tease apart the impact of mutation and selection on levels of polymorphism versus divergence in a genome that does not generally recombine.

Aging, mating, and the evolution of mtDNA heteroplasmy in Drosophila melanogaster.

Heteroplasmy, the presence of more than one type of mtDNA within cells, is common in animals and has been associated with aging and disease in humans. Changes in the frequencies of mtDNA variants between cell and animal generations thus bears on the evolution of mtDNA and the progression of diverse pathologies. We have used densitometry of Southern blots of individual heteroplasmic Drosophila melanogaster to study the effects of age, increased egg production after mating, and temperature on evolution of heteroplasmy within and between generations. The frequency of the longer mtDNA variant consistently increased between early and late cohorts of F1 offspring derived from 18 independent heteroplasmic mothers as they aged. Neither temperature (flies maintained at 25 degrees C and 18 degrees C) nor the holding of flies as virgins for 10 days before mating had significant effects on transmission patterns. However, at the ends of their lives, flies that had laid eggs at 25 degrees C had a greater frequency of the long mtDNA than did their siblings who had laid eggs at 18 degrees C. The evolution of heteroplasmy within a generation was studied in samples of siblings that either were mated or held as virgins, and then scored for mtDNA haplotype frequencies at two different ages (day 2 and day 14). Mated flies showed a significantly greater increase in the frequency of the long mtDNA variant with age than did the unmated flies. This system provides a model for the joint analysis of generational and chronological age in the transmission dynamics of a molecular polymorphism.

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans.

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.

Neutral and non-neutral evolution of Drosophila mitochondrial DNA.

To test hypotheses of neutral evolution of mitochondrial DNA (mtDNA), nucleotide sequences were determined for 1515 base pairs of the NADH dehydrogenase subunit 5 (ND5) gene in the mitochondrial DNA of 29 lines of Drosophila melanogaster and 9 lines of its sibling species Drosophila simulans. In contrast to the patterns for nuclear genes, where D. melanogaster generally exhibits much less nucleotide polymorphism, the number of segregating sites was slightly higher in a global sample of nine ND5 sequences in D. melanogaster (s = 8) than in the nine lines of D. simulans (s = 6). When compared to variation at nuclear loci, the mtDNA variation in D. melanogaster does not depart from neutral expectations. The ND5 sequences in D. simulans, however, show fewer than half the number of variable sites expected under neutrality when compared to sequences from the period locus. While this reduction in variation is not significant at the 5% level, HKA tests with published restriction data for mtDNA in D. simulans do show a significant reduction of variation suggesting a selective sweep of variation in the mtDNA in this species. Tests of neutral evolution based on the ratios of synonymous and replacement polymorphism and divergence are generally consistent with neutral expectations, although a significant excess of amino acid polymorphism within both species is localized in one region of the protein. The rate of mtDNA evolution has been faster in D. melanogaster than in D. simulans and the population structure of mtDNA is distinct in these species. The data reveal how different rates of mtDNA evolution between species and different histories of neutral and adaptive evolution within species can compromise historical inferences in population and evolutionary biology.