MBE Advance Access originally published online on November 9, 2005
Molecular Biology and Evolution 2006 23(3):559-564; doi:10.1093/molbev/msj058
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Article |
Tracing Hybrid Incompatibilities to Single Amino Acid Substitutions
Marine Biology Research Division, Scripps Institution of Oceanography, University of California, San Diego
E-mail: jsharrison{at}ucsd.edu.
| Abstract |
|---|
|
|
|---|
Deleterious interactions among genes cause reductions in fitness of interpopulation hybrids (hybrid breakdown). Identifying genes involved in hybrid breakdown has proven difficult, and few studies have addressed the molecular basis of this widespread phenomenon. Because proper function of the mitochondrial electron transport system (ETS) requires a coadapted set of nuclear and mitochondrial gene products, ETS genes present an attractive system for studying the evolution of coadapted gene complexes within isolated populations and the loss of fitness in interpopulation hybrids. Here we show the effects of single amino acid substitutions in cytochrome c (CYC) on its functional interaction with another ETS protein, cytochrome c oxidase (COX) in the intertidal copepod Tigriopus californicus. The individual and pairwise consequences of three naturally occurring amino acid substitutions in CYC are examined by site-directed mutagenesis and found to differentially effect the rates of CYC oxidation by COX variants from different source populations. In one case, we show that interpopulation hybrid breakdown in COX activity can be attributed to a single naturally occurring amino acid substitution in CYC.
Key Words: nuclear/mitochondrial coadaptation hybridization epistasis cytochrome c cytochrome c oxidase Tigriopus californicus
| Introduction |
|---|
|
|
|---|
When two genetically differentiated populations (or agricultural strains) are crossed, F1 hybrids often display some degree of enhanced performance (hybrid vigor) relative to the parental populations. However, when F1's are intercrossed, the F2's frequently show dramatic declines in performance. This phenomenon of F2 hybrid breakdown has broad implications for diverse areas of biology. For example, in evolutionary biology, hybrid breakdown is viewed as a critical step in the speciation process; in agriculture, the phenomenon underlies the common practice of planting F1 (but not F2) hybrid seed. A model first proposed by Dobzhansky (1937)
Although not specifically addressed by the Dobzhansky-Muller model, incompatibilities between nuclear and mitochondrial genes may play a significant role in hybrid breakdown due to some unique characteristics of mtDNA. First, mitochondrial genes typically have higher mutation rates than nuclear genes and commonly carry a higher number of slightly deleterious mutations (Sackton, Haney, and Rand 2003
). Second, clonal inheritance results in the accumulation of a greater number of fixed differences between isolated populations over time. Finally, mitochondrial genes are all critical subunits of electron transport system (ETS) enzyme complexes where they must interact properly with a suite of nuclear genes (see below).
The intertidal copepod Tigriopus californicus is an excellent laboratory system for investigating the consequences of interpopulation hybridization. Genetic divergence among populations from isolated rocky outcrops is high at both nuclear and mitochondrial gene loci (Burton and Lee 1994
; Burton 1998
; Edmands 1999
, 2001
; Edmands and Harrison 2003
). Mitochondrial sequence divergence exceeding 18% between populations is common, yet crosses between these populations produce fertile hybrid offspring in the laboratory (Burton 1990
). Second generation hybrids exhibit partial hybrid breakdown in a number of fitness-related traits including response to osmotic stress (Burton 1986
), developmental time (Burton 1990
), hatching number, metamorphosis number, survivorship (Edmands 1999
), and cytochrome c oxidase (COX) activity (Edmands and Burton 1999
; Rawson and Burton 2002
).
The fact that hybrid breakdown impacts a broad range of traits suggests that its molecular basis may lie in central biochemical pathways with broad physiological consequences. Energy production via the mitochondrial ETS is an attractive candidate system, not only because it meets this criterion but also because extensive genic interactions are required for ETS function (Blier, Dufresne, and Burton 2001
; Sackton, Haney, and Rand 2003
). For example, COX, which catalyzes the oxidation of cytochrome c (CYC) in the terminal step of the ETS, consists of 3 mitochondrial (COXI, COXII, and COXIII) and 10 nuclear gene products. Rawson and Burton (2002)
found that the activity of COX isolated from a Santa Cruz (SC) T. californicus population was significantly higher with CYC from SC than with CYC from a second population (San Diego, SD). Reciprocally, SD COX activity was higher with native SD CYC than with that from SC. These data provide direct evidence of functional coadaptation between CYC and COX within populations. Furthermore, Willett and Burton (2001
, 2003
) have found that the relative viability of CYC genotypes in interpopulation hybrids can differ dramatically in reciprocal crosses, implicating a direct role for cyto-nuclear interactions on fitness. Interactions between CYC and COX could underlie such fitness differences and contribute to hybrid breakdown.
In considering the interactions between CYC and COX within populations and in interpopulation hybrids, we note that COX shows much higher levels of population-specific amino acid substitution than CYC. For example, the COXII subunit alone, encoded in the mitochondrial genome, differs by 14 amino acid substitutions between the SC and SD populations (table 1). COXII is part of the catalytic core of COX, and it has been suggested that interactions between CYC and COXII are important for proper function (Roberts and Pique 1999
). In contrast, CYC genotypes from SC and SD differ by only three substitutions (table 1). Therefore, the observed functional coadaptation must be attributed to the three amino acid substitutions in CYC and their interaction with the large number of substitutions in the COX enzyme.
|
Here we report on the consequence of each CYC amino acid substitution on the observed coadaptive phenotype. Using site-directed mutagenesis, we constructed CYC variants accounting for the three amino acid substitutions and determined their effect on SC and SD COX activity. Each substitution differentially affected the rates of CYC oxidation by COX from different source populations. Furthermore, we surveyed additional natural populations for CYC variants and found a population (Punto Morro [PM], Baja California, Mexico) where one of our constructed CYC mutants occurs naturally. This PM genotype differs from the SD genotype by a single lysine to glutamine (charge changing) amino acid difference. Using population-specific COX assays, we show in this case that CYC/COX coadaptation can be attributed to a single naturally occurring amino acid substitution in CYC and its interaction with COX.
| Materials and Methods |
|---|
|
|
|---|
CYC Cloning and Expression
CYC variants from two natural T. californicus populations (SC and SD) were previously cloned into a CYC expression vector pBTR. This vector facilitates the coexpression of CYC heme lyase which is necessary for the heterologous expression of functional CYC in Escherichia coli (Pollock et al. 1998
CYC Purification
Cells from 1.6 liters of culture for each CYC variant were pelleted and resuspended in 150 ml of 100 mM NaH2PO4 buffer (pH 7.1) with 0.5 g lysozyme and crystals of DNase and RNase. Suspensions were incubated at 4°C, with mixing, for 24 h. Cellular debris was pelleted, and the supernatant was removed and stored at 4°C. This process was repeated two to three times until the pellet was no longer pink in color. The supernatant from each round of digestion was combined and passed through an Ultrafree-15 Centrifugal filter (Millipore) three times, each time resuspending in 0.05 M potassium phosphate buffer (PB) at pH 7.3. After the final pass through the filter, each sample was resuspended in 3 ml PB. CYC from each sample was then reduced with 0.5 mg ascorbic acid. Ascorbate was removed using a PD-10 desalting column (Amersham Biosciences, Piscataway, N.J.) equilibrated with PB.
COX Activity Assay
Homogenates prepared from 30 adult females were used as a source of population-specific COX. Animals were homogenized in 600 µl of PB with 0.05% Tween 20 and then centrifuged at 5000xg for 5 min; pellets were discarded. Activity of COX for each population was measured with each CYC variant at both 18 and 25°C and replicated four times. Reaction rates were determined by monitoring the oxidation of CYC at 550 nm. Reactions were performed by adding 25 µl of the homogenate to 400 µl of CYC (20 µM) and monitored for 2 min. Reaction rates were calculated using the best linear fit to the 2-min trajectories. Rates were normalized to protein concentration of homogenates determined by standard bicinchoninic acid method (Pierce) using bovine serum albumin as a standard.
| Results and Discussion |
|---|
|
|
|---|
Using site-directed mutagenesis, we constructed CYC variants and tested each of the three amino acid substitutions for their effect on activity of SC and SD COXs (table 1). COX activity assays demonstrate significant functional coadaptation between CYC and COX within populations (Rawson and Burton 2002
0.01).
|
Models of CYC/COX interactions predict that the dipole moment or surface charge potential is critical for proper docking of CYC (Maneg et al. 2004
|
The large temperature effect on COX activity observed here has been noted previously (Rawson and Burton 2002
Because we observed a large effect of individual amino acids on COX activity, we extended a previous survey of T. californicus populations to obtain additional natural CYC variants (Rawson, Brazeau, and Burton 2000
). Copepods obtained from PM, Baja California, Mexico (31°52'N, 116°40'W) possessed a CYC variant identical to the experimental CYC mutant K43Q that is differing from SD by the single change found to be most important in the SC/SD comparison. Our finding of multiple T. californicus populations possessing charge differences in CYC is itself remarkable. In contrast, there are almost no charge differences in CYC among divergent lineages of primates (Schmidt et al. 2005
).
Interestingly, PM COXII also possesses one less acidic residue than SD COXII (table 1). Results of intra- and interpopulation COX/CYC activity assays were again consistent with functional breakdown in the interpopulation assay (fig. 3). Three-way ANOVA again shows the main effects of COX, CYC, and temperature, and all two-way interactions are significant (P < 0.001). The activity of PM COX was consistently higher with PM CYC than with SD CYC at both 18 and 25°C (P = 0.002 [18°C], P < 0.001 [25°C], t-test). SD COX activity was significantly higher with SD CYC than with PM CYC regardless of temperature (P = 0.013 [18°C], P = 0.023 [25°C], t-test). The significant two-way interaction (ANOVA, P < 0.001) and the consistently higher activities of within-population CYC/COX pairs strongly suggest that the observed intrapopulation CYC/COX coadaptation can be attributed to this single amino acid difference between PM and SD CYCs and its functional interaction with COX.
|
There are a number of studies that have considered the coevolution of ETS proteins within species. These studies typically include sequence analysis to identify correlated rate changes among lineages (Jobson et al. 2004
Our study was designed to determine which of the three amino acid differences between SD and SC CYCs are required for coadaptation with SD or SC COXs. Results suggest that the K43Q substitution accounts for a large portion of the CYC/COX coadaptation in SC, a result that is consistent with hypotheses concerning charge potential and electrostatic interactions between CYC and COX. In contrast, although SD COX activity is not independent of the CYC source, none of the individual or pairwise amino acid changes explained the decreased activity with SC CYC. Here we suggest that all three amino acid changes in CYC and their interaction with the large number of substitutions in the COX enzyme jointly result in reduced activity of SD COX with SC CYC. Our discovery of the PM population allowed us to explicitly test the consequences of a single amino acid change on interpopulation CYC/COX activity. Again, intrapopulation COX/CYC combinations resulted in higher COX activity than interpopulation combinations. To our knowledge, this is the first case in natural populations, where the functional basis for functional coadaptation at the molecular level can be isolated to a single amino acid change.
| Acknowledgements |
|---|
|
|
|---|
We thank J. Flowers and C. Ellison for helpful comments on the manuscript. We gratefully acknowledge support of a National Science Foundation grant to R.S.B.
| Footnotes |
|---|
Peter Lockhart, Associate Editor
| References |
|---|
|
|
|---|
Blier, P. U., F. Dufresne, and R. S. Burton. 2001. Natural selection and the evolution of mtDNA-encoded peptides: evidence for intergenomic co-adaptation. Trends Genet. 17:400406.[CrossRef][ISI][Medline]
Burton, R. S. 1986. Evolutionary consequences of restricted gene flow among natural populations of the copepod Tigriopus californicus. Bull. Mar. Sci. 39:526535.
. 1990. Hybrid breakdown in physiological response: a mechanistic approach. Evolution 44:18061813.[CrossRef]
. 1998. Intraspecific phylogeography across the point conception biogeographic boundary. Evolution 52:734745.[CrossRef][ISI]
Burton, R. S., and B.-N. Lee. 1994. Nuclear and mitochondrial gene genealogies and allozyme polymorphism across a major phylogeographic break in the copepod Tigriopus californicus. Proc. Natl. Acad. Sci. USA 91:51975201.
Coyne, J. A., and H. A. Orr. 1997. "Patterns of speciation in Drosophila" revisited. Evolution 51:295303.
Dobzhansky, T. 1937. Genetics and the origin of species. Columbia University Press, New York.
Drosou, V., F. Malatesta, and B. Ludwig. 2002. Mutations in the docking site for cytochrome c on the Paracoccus heme aa3 oxidase. Electron entry and kinetic phases of the reaction. Eur. J. Biochem. 269:29802988.[Medline]
Edmands, S. 1999. Heterosis and outbreeding depression in interpopulation crosses spanning a wide range of divergence. Evolution 53:17571765.[CrossRef]
. 2001. Phylogeography of the intertidal copepod Tigriopus californicus reveals substantially reduced population differentiation at northern latitudes. Mol. Ecol. 10:17431750.[CrossRef][Medline]
Edmands, S., and R. S. Burton. 1999. Cytochrome c oxidase activity in interpopulation hybrids of a marine copepod: a test for nuclear-nuclear or nuclear-cytoplasmic coadaptation. Evolution 53:19721978.[CrossRef]
Edmands, S., and J. S. Harrison. 2003. Molecular and quantitative trait variation within and among populations of the intertidal copepod Tigriopus californicus. Evolution 57:22772285.[CrossRef][ISI][Medline]
Jobson, R. W., R. Nielson, L. Laakkonen, M. Wikstrom, and V. A. Albert. 2004. Adaptive evolution of cytochrome c oxidase: infrastructure for a carnivorous plant radiation. Proc. Natl. Acad. Sci. USA 101:1806418068.
Maneg, M., F. Malatesta, B. Ludwig, and V. Drosou. 2004. Interaction of cytochrome c with cytochrome c oxidase: two different docking scenarios. Biochem. Biophys. Acta 1655:274281.[Medline]
Muller, H. J. 1940. Bearing of the "Drosophila" work in systematics. Pp. 185268 in J. Huxley, ed. The new systematics. Oxford University Press, Oxford.
Orr, H. A., J. P. Masly, and D. C. Presgraves. 2004. Speciation genes. Curr. Opin. Genet. Dev. 14:675679.[CrossRef][ISI][Medline]
Osheroff, N., S. H. Speck, E. Margoliash, E. C. Veerman, J. Wilms, B. W. Konig, and A. O. Muijsers. 1983. The reaction of primate cytochromes c with cytochrome c oxidase: analysis of the polarographic assay. J. Biol. Chem. 258:57315738.
Pollock, W. B. R., F. I. Rosell, M. B. Twitchett, M. E. Dumont, and A. G. Mauk. 1998. Bacterial expression of a mitochondrial cytochrome c. Trimethylation of Lys72 in yeast iso-1-cytochrome c and the alkaline conformational transition. Biochemistry 37:61246131.[CrossRef][Medline]
Presgraves, D. C. 2003. A fine-scale genetic analysis of hybrid incompatibilities in Drosophila. Genetics 163:955972.
Rawson, P. D., and R. S. Burton. 2002. Functional coadaptation between cytochrome c and cytochrome c oxidase within allopatric populations of a marine copepod. Proc. Natl. Acad. Sci. USA 99:1295512958.
Rawson, P. D., D. A. Brazeau, and R. S. Burton. 2000. Isolation and characterization of cytochrome c from the marine copepod Tigriopus californicus. Gene 248:1522.[CrossRef][ISI][Medline]
Roberts, V. A., and M. E. Pique. 1999. Definition of the interaction domain for cytochrome c on cytochrome c oxidase. III. Prediction of the docked complex by a complete, systematic search. J. Biol. Chem. 274:3805138060.
Sackton, T. B., R. A. Haney, and D. M. Rand. 2003. Cytonuclear coadaptation in Drosophila: disruption of cytochrome c oxidase activity in backcross genotypes. Evolution 57:23152325.[CrossRef][ISI][Medline]
Schmidt, T. R., D. E. Wildman, M. Uddin, J. C. Opazo, M. Goodman, and L. I. Grossman. 2005. Rapid electrostatic evolution at the binding site for cytochrome c on cytochrome c oxidase in anthropoid primates. Proc. Natl. Acad. Sci. USA 102:63796384.
Turelli, M., N. H. Barton, and J. A. Coyne. 2001. Theory and speciation. Trends Ecol. Evol. 16:330342.[CrossRef][Medline]
Willett, C. S., and R. S. Burton. 2001. Viability of cytochrome c genotypes depends on cytoplasmic backgrounds in Tigriopus californicus. Evolution 55:15921599.[CrossRef][ISI][Medline]
. 2003. Environmental influences on epistatic interactions: viabilities of cytochrome c genotypes in interpopulation crosses. Evolution 57:22862292.[CrossRef][ISI][Medline]
Wu, C.-I., and C.-T. Ting. 2004. Genes and speciation. Nat. Rev. Genet. 5:114127.[CrossRef][ISI][Medline]
Zhen, Y., C. W. Hoganson, G. T. Babcock, and S. Ferguson-Miller. 1999. Definition of the interaction domain for cytochrome c on cytochrome c oxidase. I. Biochemical, spectral, and kinetic characterization of surface mutants in subunit II of Rhodobacter sphaeroides cytochrome aa3. J. Biol. Chem. 274:3803238041.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
D. I. Bolnick, M. Turelli, H. Lopez-Fernandez, P. C. Wainwright, and T. J. Near Accelerated Mitochondrial Evolution and "Darwin's Corollary": Asymmetric Viability of Reciprocal F1 Hybrids in Centrarchid Fishes Genetics, February 1, 2008; 178(2): 1037 - 1048. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. S. Willett Significant Variation for Fitness Impacts of ETS Loci in Hybrids between Populations of Tigriopus californicus J. Hered., January 1, 2008; 99(1): 56 - 65. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Niehuis, A. K. Judson, and J. Gadau Cytonuclear Genic Incompatibilities Cause Increased Mortality in Male F2 Hybrids of Nasonia giraulti and N. vitripennis Genetics, January 1, 2008; 178(1): 413 - 426. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. F. McDaniel, J. H. Willis, and A. J. Shaw A Linkage Map Reveals a Complex Basis for Segregation Distortion in an Interpopulation Cross in the Moss Ceratodon purpureus Genetics, August 1, 2007; 176(4): 2489 - 2500. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Turelli and L. C. Moyle Asymmetric Postmating Isolation: Darwin's Corollary to Haldane's Rule Genetics, June 1, 2007; 176(2): 1059 - 1088. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Burzynski, M. Zbawicka, D. O. F. Skibinski, and R. Wenne Doubly Uniparental Inheritance Is Associated With High Polymorphism for Rearranged and Recombinant Control Region Haplotypes in Baltic Mytilus trossulus Genetics, November 1, 2006; 174(3): 1081 - 1094. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

indicates P < 0.05 when tested against SC. indicates P < 0.05 when tested against SD.


