help button home button Genetics J Virology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Originally published as Genetics Published Articles Ahead of Print on April 19, 2006.

Genetics, Vol. 173, 1465-1477, July 2006, Copyright © 2006
doi:10.1534/genetics.105.051011

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
genetics.105.051011v1
173/3/1465    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Willett, C. S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Willett, C. S.

Deleterious Epistatic Interactions Between Electron Transport System Protein-Coding Loci in the Copepod Tigriopus californicus

Christopher S. Willett1

Department of Biology, University of North Carolina, Chapel Hill, North Carolina 27599-3280

1 Address for correspondence: Department of Biology, University of North Carolina, CB 3280 Coker Hall, Chapel Hill, NC 27599-3280.
E-mail: willett4{at}email.unc.edu

Manuscript received September 12, 2005. Accepted for publication April 12, 2006.


   ABSTRACT
TOP
 ABSTRACT
MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
The nature of epistatic interactions between genes encoding interacting proteins in hybrid organisms can have important implications for the evolution of postzygotic reproductive isolation and speciation. At this point very little is known about the fitness differences caused by specific closely interacting but evolutionarily divergent proteins in hybrids between populations or species. The intertidal copepod Tigriopus californicus provides an excellent model in which to study such interactions because the species range includes numerous genetically divergent populations that are still capable of being crossed in the laboratory. Here, the effect on fitness due to the interactions of three complex III proteins of the electron transport system in F2 hybrid copepods resulting from crosses of a pair of divergent populations is examined. Significant deviations from Mendelian inheritance are observed for each of the three genes in F2 hybrid adults but not in nauplii (larvae). The two-way interactions between these genes also have a significant impact upon the viability of these hybrid copepods. Dominance appears to play an important role in mediating the interactions between these loci as deviations are caused by heterozygote/homozygote deleterious interactions. These results suggest that the fitness consequences of the interactions of these three complex III-associated genes could influence reproductive isolation in this system.

EPISTASIS contributes to a wide range of evolutionary phenomena such as reproductive isolation, the potential for the evolution of sex, and the nature of adaptation (WHITLOCK and PHILLIPS 1995). However, specific interactions between individual loci and the impact of these interactions on fitness are largely unknown for naturally occurring variation in most animal species. In the context of speciation, epistasis plays a primary role in postzygotic reproductive isolation, which is thought to be caused primarily by negative fitness interactions between loci. These deleterious interactions are known as Dobzhansky–Muller (D–M) incompatibilities (the model was originally proposed by DOBZHANSKY 1936 and MULLER 1942). D–M incompatibilities are hypothesized to result from the breakup of positive epistasis that develops between interacting genes within distinct evolutionary lineages (coadapted gene complexes). These interactions are generally thought to be between two or more nuclear-encoded genes; however, rapidly evolving mitochondrial-encoded genes in many taxa provide another potential partner in deleterious interactions. Despite the crucial role of D–M incompatibilities in postzygotic reproductive isolation, there are very few examples of the types of specific molecular interactions that are involved. Several cases are known from Drosophila, where one gene involved in an interaction has been identified but the interacting partner(s) and the physiological process that are disrupted remain to be elucidated (TING et al. 1998; BARBASH et al. 2003, 2004; PRESGRAVES et al. 2003). Here I look for evidence that nuclear and mitochondrial-encoded proteins in the electron transport system (ETS) in Tigriopus californicus could be involved in deleterious interactions and contribute to reproductive isolation in hybrids from population crosses of this species.

In taxa with partial reproductive isolation that can be crossed, patterns of segregation of molecular markers or chromosomal regions can be examined in hybrid offspring to determine candidate regions for D–M incompatibilities. Studies of crosses between species often reveal widespread detrimental interactions (where nonparental genotype combinations have a negative effect on fitness) and occasional positive interactions (conversely, when nonparental genotype combinations have increased fitness) (PALOPOLI and WU 1994; RIESEBERG et al. 1996; BURKE et al. 1998; COYNE et al. 1998; FISHMAN et al. 2001; PRESGRAVES 2003). Other studies have shown evidence for deleterious epistatic interactions in crosses between divergent populations of the same species (DOBZHANSKY 1950; ARMBRUSTER et al. 1997; LI et al. 1997; GALLOWAY and FENSTER 1999; FENSTER and GALLOWAY 2000a,b; HALL and WILLIS 2005). One caveat for many of these segregation studies is that the observed patterns of non-Mendelian inheritance could result from meiotic drive rather than from the breakdown of coadapted gene complexes; this has been shown to play a role for at least one well-characterized locus in an interspecific plant cross (FISHMAN and WILLIS 2005). Studies of segregation of genotypes in hybrids will sometimes allow the proportion of heterozygous/homozygous (where one locus is homozygous for the same allele and the second locus is heterozygous for two alleles) to homozygous/homozygous (both loci are homozygous in a two-locus comparison) incompatibilities to be estimated. The ratio of these types of incompatibilities can have significant effects on the results obtained from models of postzygotic isolation (TURELLI and ORR 2000).

The intertidal copepod T. californicus has proven to be a useful system in which to examine the initial stages of allopatric speciation. Geographically distinct populations of this copepod from the Pacific coast of North America are dramatically divergent in mitochondrial DNA (mtDNA) even over relatively short distances (BURTON and LEE 1994; BURTON 1998; BURTON et al. 1999; EDMANDS 2001; WILLETT and BURTON 2004). The changes in mtDNA include a large number of replacement differences leading to substantial divergence in these mtDNA-encoded proteins. Crosses between these copepod populations typically show that F2 hybrids (but not F1 hybrids) have lower fitness by a variety of measures when each class is compared to the parental populations, suggesting that coadapted gene complexes are disrupted in these hybrids (BURTON 1986, 1987, 1990; EDMANDS 1999; EDMANDS and DEIMLER 2004). Some genetic regions that contribute to this reduction in F2 hybrid fitness have been identified by segregation studies of molecular markers in crosses of divergent copepod populations (BURTON 1987; HARRISON and EDMANDS 2006).

The intimate interactions between mtDNA-encoded proteins and nuclear-encoded proteins in the physiologically indispensable ETS, coupled with the high rate of divergence in mtDNA in T. californicus, suggested this system as a candidate system in which to examine the evolution of epistasis and its potential role in F2 hybrid breakdown. Several studies have suggested that mitochondrial/nuclear coadaptation has diverged in these genetically isolated copepod populations: EDMANDS and BURTON (1999) found evidence for mitochondrial/nuclear coadaptation in backcross hybrids between populations. RAWSON and BURTON (2002) uncovered functional coadaptation between nuclear-encoded cytochrome c (CYC) and complex IV (which contains three mtDNA-encoded proteins) in studies of enzymatic rates. WILLETT and BURTON (2001, 2003) demonstrated that the fitness of CYC genotypes in hybrids is influenced by the cytoplasmic type of the cross and one specific cross yielded results consistent with coadaptation between CYC and mtDNA-encoded proteins in both F2 and advanced generation hybrids. In this cross [between populations from Santa Cruz (SC) and Abalone Cove (AB)] the selection against non-mtDNA type CYC genotypes occurred during the juvenile life stages (postzygote), thereby ruling out meiotic drive as an explanation for skewed ratios. Results obtained in the prior studies of coadaptation between mtDNA-encoded and nuclear-encoded proteins appear to be highly influenced by the particular pair of populations crossed (and the direction of the cross), suggesting that independent evolution has occurred in these coadaptation systems in each lineage.

Complex III of the ETS (also known as ubiquinol/cytochrome c oxioreductase) includes both the mtDNA-encoded cytochrome b (CYTB) and a set of nuclear-encoded subunits (10 in mammals), 2 of which, cytochrome c1 (CYC1) and the rieske iron–sulfur protein (RISP), contain redox centers (Figure 1) and are found in the bacterial complex III homolog (IWATA et al. 1998). These three functionally important components of complex III have been characterized from several T. californicus populations and have been shown to have differences in amino acid sequence in comparisons across populations (WILLETT and BURTON 2004). Much like the pattern seen for other mtDNA genes (BURTON and LEE 1994; BURTON 1998; BURTON et al. 1999; EDMANDS 2001), CYTB was remarkably divergent between T. californicus populations with synonymous pairwise divergence exceeding 65% for comparisons between SC, AB, and San Diego (SD) copepod populations. In contrast, synonymous divergences for the CYC1 and RISP genes were ~5% for comparisons of these same sets of populations (WILLETT and BURTON 2004). Complex III interacts closely with CYC (a soluble protein that transfers electrons from complex III to complex IV; SARASTE 1999); CYC has also diverged between copepod populations, divergence that has potentially been driven by positive selection (RAWSON et al. 2000; WILLETT and BURTON 2004).


Figure 1
View larger version (15K):
[in this window]
[in a new window]
[Download PPT slide]
  		FIGURE 1.— Schematic of complex III and its position in the electron transport system. Complex III is known to function as a dimer in the inner membrane of the mitochondrion as depicted in the diagram. Within each half, the cytochrome b protein (CYTB) is mtDNA encoded, while the rieske iron–sulfur protein (RISP), cytochrome c1 (CYC1) proteins, and the soluble protein cytochrome c (CYC) are all encoded in the nucleus. Ubiquinol (QH2) is a lipid soluble compound that moves within the membrane and is oxidized by complex III. The diagram is derived from information in SARASTE (1999) and LANGE and HUNTE (2002).

 
Given the importance of these ETS proteins for energy generation and the presence of changes in the amino acid sequences between copepod populations, an unanswered question is whether these differences in proteins will lead to detrimental interactions within the ETS that will negatively influence fitness in hybrids from genetically divergent populations. An examination of genotypic viabilities in F2 hybrids from T. californicus populations for these three interacting proteins (CYC, CYC1, and RISP) and their combined impact on viability affords the opportunity to uncover the nature of nuclear/nuclear interactions as well as nuclear/mtDNA interactions for this system. The degree to which heterozygote/homozygote vs. homozygote/homozygote types of interactions are occurring is also examined. Additional comparisons of levels of epistasis between these three proteins with two non-ETS proteins (malic enzyme homologs, ME1 and ME2) will help to illuminate the level of background interactions expected for arbitrarily chosen proteins. The results presented in this article suggest that for a cross between two Southern California copepod populations, Abalone Cove and San Diego, interactions between these three complex III nuclear proteins are having an impact upon F2 hybrid fitness. These interactions appear to be nuclear/nuclear and primarily heterozygote/homozygote in type and could be responsible for deleterious interactions that act as D–M incompatibilities in these hybrid copepods.


   MATERIALS AND METHODS
TOP
 ABSTRACT
 MATERIALS AND METHODS
RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Crosses and culturing conditions:
T. californicus were collected from intertidal rock pools at two sites in southern California, SD (32°44.74' N, 117°15.30' W, San Diego County) and AB (33°44.26' N, 118°22.52' W, Los Angeles County). The copepods were maintained in mass culture in artificial seawater (Instant Ocean, Aquarium Systems) in 400-ml beakers at 20° with a 12:12 L:D photoperiod. Cultures were maintained at a concentration of 35 parts per thousand seawater and fed with commercial flake fish food (O.S.I. Spirulina, Ocean Science International), but copepods were also allowed to consume natural algal growth.

To test for the existence of epistatic interactions between ETS proteins, crosses were set up using copepods from the AB and SD populations that had been maintained under the culturing conditions described above for at least 1 year (multiple generations are likely to have occurred under these laboratory conditions). Virgin females were obtained by separating clasped pairs (BURTON 1985) and placing these females in dishes with males from the other population. Crosses were done in petri dishes with 15 males and 15 females, and two plates were set up for each of the reciprocal crosses (ABf x SDm and ABm x SDf). Female copepods were moved to new dishes when copepodids of the next generation were observed (if copepodids molt to the final adult stage, then generations will be indistinguishable). Males were removed and discarded because females mate only once in their lifetime and store sufficient sperm to produce multiple clutches. F2 progeny were obtained by allowing F1 male and female copepods to mate after replicate crosses were mixed (to minimize the chance of inbreeding). Nucleotide variability (which should be a proxy for functionally important genetic variation) is much lower within these populations than between them, reflecting the apparent long-term restrictions on gene flow between these populations. F1 females were then transferred to new petri dishes when F2 copepodids were observed. This crossing design has the potential disadvantage of not allowing cross-to-cross variability to be uncovered but it does permit a large number of F2 progeny to be collected, thereby increasing power to detect epistatic interactions.

F2 copepods were collected as both nauplii and adults: First-stage nauplii were obtained by incubating mature egg sacs (red in color) from F1 females in seawater at which point nauplii would hatch out and were immobilized on filter paper and then collected into 10 µl of proteinase-K cell-lysis buffer (HOELZEL and GREEN 1992). Adults were separated according to sex and then collected and placed into 20 µl of proteinase-K cell-lysis buffer. Both adults and nauplii were placed in 96-well polymerase chain reaction (PCR) plates and incubated at 65° for 1 hr and 95° for 15 min to prepare DNA for subsequent PCR reactions. Two 96-well plates each of nauplii, adult males, and adult females were collected and stored at –20°.

ETS gene marker scoring and data analysis:
PCR-based markers were developed for each of the three ETS proteins to allow the scoring of the genotypes of hybrids by length differences, i.e., allele-specific PCR products of different sizes. Genotyping assays for the gene CYC were described in WILLETT and BURTON (2001). Primers for RISP and CYC1 took advantage of length differences in introns between AB and SD populations uncovered in these sequences in WILLETT and BURTON (2004). The primers for RISP had the following sequences: RISPcon.f 5'-TGGCAGATCCAAGGCAAATACTAACAAT-3' and RISPcon.r 5'-TGGTTAGATTTGAATACGGAACCTTCA-3'. The CYC1 primers were CYC1ex3con.f 5'-GACTCGTTTGACCATGCTTCCATTCG-3' and CYC1ex4con.r 5'-CGGCATGGGTAGATAGTCGGTCAAT-3'. PCR reactions were run under standard conditions (with 55° annealing temperatures for both sets of reactions) and genotypes were determined by visualizing differences in the size of PCR products on 1.5% agarose gels stained with ethidium bromide.

Each of the complex III markers was tested for departures from homogeneity between sexes, using contingency tests. Deviations from Mendelian 1:2:1 ratios were then examined for the sexes separately or combined. This test can be used to test the significance of differences in genotypic viability that can then be compared across reciprocal crosses (with different mtDNA backgrounds) to determine whether mitochondrial/nuclear interactions are influencing F2 hybrid fitness. Relative viabilities of homozygous genotypes and their standard deviations were calculated according to HALDANE's (1956) formulas. For the estimated viability value this is simply wAA = 2NAA/(NAa + 1), where NAA is the number of one homozygous class and NAa is the number of heterozygotes. The standard deviation of this value is

Formula

To test for epistatic interactions between nuclear-encoded gene products (nuclear/nuclear interactions), a multiplicative model was used in which the null hypothesis was that single-locus effects are mutually independent to generate expected numbers without epistasis. The expected number of each two-locus genotype class was calculated by using observed single-locus genotype frequencies to correct for deviations associated with each of the two single loci. This was done for each genotypic combination (nine total) by multiplying the two observed genotypic frequencies of the two loci and then multiplying these by the total number of copepods scored for both markers. This simplifies to

Formula
where Exp NAABB is the expected number of a two-locus genotypic class with corresponding numbers of NAA and NBB observed for each locus individually and Ntotal equal to the total number of copepods scored for both loci. Three-way epistatic interactions between RISP, CYC1, and CYC were detected using a procedure analogous to that described above for two-way interactions. The expected number for each of the 27 three-locus classes was derived from the average of each of the three two-locus classes multiplied by the excluded locus class. This simplifies to the formula

Formula
where Exp NAABBCC is the expected number of a three-locus genotypic class; the observed number in each two-locus class is NAABB, NAACC, or NBBCC; the observed number of each single-locus class is NAA, NBB, or NCC; and NTotal is the total sample size for which all three genotypes were obtained. The significance of deviations from independence in single-, two-, and three-locus comparisons was tested with a {chi}2-test. A sequential Bonferroni correction was used to determine a critical P-value for each table that adjusted for multiple tests. The significance of interactions with the functionally distantly related malic enzyme homologs (ME1 and ME2) was calculated using genotypes collected from the same F2 individuals scored herein (genotypes from J. BERKOWITZ and C. S. WILLETT, unpublished data). These markers (in addition to two-way interaction results for histone H1 from WILLETT and BURTON 2001) will provide some insight into the background level of epistatic interactions in these crosses.


   RESULTS
TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Deviations from Mendelian inheritance for ETS markers:
To test the potential contribution of nuclear-encoded ETS genes to hybrid breakdown in population crosses in T. californicus the genotypic-specific viabilities of the three ETS loci were examined individually. The genotypes for CYC, CYC1, and RISP were determined both from hybrid F2 adults and from first-stage nauplii from two reciprocal crosses between AB and SD populations (Table 1). There are significant deviations from Mendelian inheritance for each of these markers in both of the reciprocal crosses for F2 adults but no significant deviations are observed for F2 nauplii, suggesting that fitness differences between genotypic classes are leading to selection during development from nauplii to adult for each of these three markers in these hybrids. The ABf x SDm cross shows significant differences between the sexes for RISP and CYC1 genotypic ratios (only for CYC1 when significance values are corrected for multiple comparisons). For all three ETS markers there is a marked reduction in viability for both homozygote classes in adults when compared to the heterozygous class (Table 2).WILLETT and BURTON (2001) performed an analogous cross under similar conditions and scored AB x SD F2 adults for CYC; these results are summarized in Table 2 for comparison.


View this table:
[in this window]
[in a new window]

  		TABLE 1 Observed genotypic numbers for three ETS genes in F2 progeny from interpopulation crosses of Tigriopus californicus between Abalone Cove (AB) and San Diego (SD)

 

View this table:
[in this window]
[in a new window]

  		TABLE 2 Relative viabilities of ETS markers for F2 hybrids in AB x SD population crosses

 
The potential for CYC coadaptation with mtDNA-encoded proteins has been suggested by a number of studies (EDMANDS and BURTON 1999; WILLETT and BURTON 2001, 2003); however, these studies have focused largely on a different population cross (Abalone Cove and Santa Cruz in central California). Both RISP and CYC1 are found in the same protein complex as CYTB, which sets up the potential for the coadaptation of these nucDNA-encoded proteins with the mtDNA-encoded CYTB. To test for nuclear/mitochondrial coadaptation it is instructive to look at the representation of the alternate homozygous genotypic classes with the expectation that there will be a 1:1 ratio between them (or equivalent viabilities when compared to the heterozygous genotypic class). An examination of the viabilities in Table 2 does not lend support to a large role of mitochondrial/nuclear coadaptation in determining F2 hybrid fitness in the AB x SD crosses for these three markers. There is not widespread agreement with a simple model of ETS gene/mtDNA coadaptation when comparing bias in alternate homozygous classes with the maternal parent in each cross for each of the three genes. With maternal inheritance of mtDNA, nuclear/mtDNA coadaptation would predict that the viability of the homozygous class not matching the female parent should be reduced. For the current results, only three of nine possible comparisons between homozygous classes in Table 2 (underlined) show biases matching this prediction for the current results.

Epistatic interactions between ETS genes:
Given the close interactions of CYC, CYC1, and RISP proteins in complex III of the ETS, hybridization could result in the breakdown of coadaptation between these proteins in this complex. These nuclear/nuclear gene interactions within complex III could then be a potential source of hybrid inviability. To test for the presence of this coadaptation, I have examined the patterns of epistasis observed between the markers in the three complex III-associated genes in the F2 hybrid populations examined above. There are substantial deviations from independence for two-way comparisons of these three markers in adult F2 hybrids but no evidence of departures in F2 nauplii (Table 3). The departures from multiplicative independence were calculated by using expected values that correct for the fitness effects of each of the single loci. The two comparisons involving CYC1 display significant deviations for the adults summed across the reciprocal crosses and large deviations (approaching significance only when correcting for multiple tests) for at least one of each of the reciprocal crosses. For the RISP/CYC comparison, only the total female sample has a marginally significant deviation from independence (approaching the critical sequential Bonferroni-corrected P-value).


View this table:
[in this window]
[in a new window]

  		TABLE 3 Two-way epistatic interactions between ETS loci in F2 hybrids

 
The nature of deviations from independence could yield insights into the types of epistasis occurring between these genes. The simplest pattern that could be generated by D–M incompatibilities would be an excess of parental-type double homozygotes and a deficit of non-parental-type double homozygotes. However, more complex patterns of homozygote/heterozygote incompatibilities are also possible (modeled by TURELLI and ORR 2000). Figure 2 separates the components of the {chi}2-deviations into the nine different two-locus genotypic classes for CYC and CYC1. This provides a qualitative assessment of which interactions are overrepresented vs. underrepresented for these two genes in these hybrid copepods. For the CYC and CYC1 genotypic combinations, the homozygote/heterozygote classes are all underrepresented, suggesting their potential involvement in deleterious interactions. Doubly parental SD homozygotes (e.g., CYCSD/SD/CYC1SD/SD) show large positive deviations (as does the nonparental type, CYCSD/SD/CYC1AB/AB).


Figure 2
View larger version (26K):
[in this window]
[in a new window]
[Download PPT slide]
  		FIGURE 2.— Relative contributions to {chi}2-deviations for genotype-by-genotype classes in F2 hybrids of a cross between SD and AB for T. californicus for CYC and CYC1. The magnitude of the deviation in {chi}2 is shown by the length of the bars on the y-axis while the direction indicates more (positive) or less (negative) in each category than expected. The genotypic combinations are indicated at the top while the reciprocal cross or combined result is indicated by shading as shown in the inset. Boxed classes are heterozygous/homozygous comparisons while the two underlined classes are nonparental homozygous/homozygous comparisons. The {chi}2-value for each cross and associated P-value are given below each cross name in the inset.

 
A similar pattern is found in the RISP/CYC1 comparison for deviations from multiplicative independence (Figure 3). Once again the heterozygous/homozygous classes are underrepresented. The large positive effect of the RISPSD/SD/CYC1AB/AB class and the similar result seen for a nonparental genotype with CYC and CYC1 are not consistent with completely negative interactions between genes from different populations. The deviations for RISP/CYC (for which the overall {chi}2-values were largely marginally significant or insignificant) are shown in Figure 4. For these two genes the largest contribution to the overall {chi}2-deviation comes from positive deviations for both parental-type double homozygous classes. The heterozygous/homozygous genotypic combinations are once again all underrepresented from RISP/CYC.


Figure 3
View larger version (25K):
[in this window]
[in a new window]
[Download PPT slide]
  		FIGURE 3.— Relative contributions to {chi}2-deviations for genotype-by-genotype classes in F2 hybrids of a cross between SD and AB for T. californicus for CYC1 and RISP. The magnitude, direction, and cross type are all as indicated in Figure 1. Boxed classes are heterozygous/homozygous comparisons while the two underlined classes are nonparental homozygous/homozygous comparisons. The {chi}2-value for each cross and associated P-value are given below each cross name in the inset. Note the difference in scale of the y-axis from Figure 2.

 

Figure 4
View larger version (25K):
[in this window]
[in a new window]
[Download PPT slide]
  		FIGURE 4.— Relative contributions to {chi}2-deviations for genotype-by-genotype classes in F2 hybrids of a cross between SD and AB for T. californicus for CYC and RISP. The magnitude, direction, and cross type are all as indicated in Figure 1. Boxed classes are heterozygous/homozygous comparisons while the two underlined classes are nonparental homozygous/homozygous comparisons. The {chi}2-value for each cross and associated P-value are given below each cross name in the inset. Note the difference in scale of the y-axis from Figures 2 and 3.

 
There is limited evidence for significant three-way interactions between RISP, CYC, and CYC1 when the effects of the two-way interactions are factored out (Table 4). When all adults in both reciprocal crosses are summed (total adults in Table 4), they deviate substantially from independence between the three ETS genes but do not reach the significance threshold with a sequential Bonferroni correction. The total adult genotypic class deviates from independence by having an observed excess of triple heterozygotes, and, in addition, several double-heterozygote/single-homozygote classes are underrepresented (results not shown). For the three-way interactions as well as the two-way interactions previously, the female classes show greater deviations from independence than do the male classes. Although the interactions between these three ETS genes appear to be complex, it appears that their interactions have significant effects upon hybrid fitness in this population cross.


View this table:
[in this window]
[in a new window]

  		TABLE 4 Three-way interactions between RISP, CYC1, and CYC genes in F2 hybrids

 
Interactions of ETS proteins with ME homologs:
Markers have been developed for two different malic enzyme homologs from T. californicus, and these have been used to genotype the same F2 hybrids genotyped for the three nuclear-encoded ETS genes in this study (J. BERKOWITZ and C. S. WILLETT, unpublished results). There were significant departures from 1:2:1 Mendelian expectations for these two loci when examined singly that involved greatly reduced numbers of both homozygous classes (results not shown). We can gain some insight into the levels of background epistasis by looking at two-way interactions between the ME proteins and the three ETS proteins in the F2 hybrids. The same multiplicative model of independent gene action was tested for the comparisons of each of these two ME markers with the three ETS genes (Table 5). None of the deviations were significant when corrected for multiple tests using the sequential Bonferroni procedure. The ABf x SDm cross for ME/CYC had a large departure from independence that approached the sequential Bonferroni-adjusted P-value. This departure was driven by an excess of individuals in the CYCSD/SD/MEAB/AB class that led to a significant deviation for the F2 adult females; F2 males of the reciprocal cross had a nonsignificant trend for a deviation in the same genotypic class. In comparison with the results involving the three ETS proteins, these data suggest fewer fitness interactions in hybrids between these two ME homologs with the ETS proteins (and linked genes).


View this table:
[in this window]
[in a new window]

  		TABLE 5 Two-way epistatic interactions between ETS loci and malic enzyme homologs in F2 hybrids

 

   DISCUSSION
TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Significant epistatic interactions between ETS proteins:
The three genes, CYC, CYC1, and RISP, encode proteins crucial for the function of complex III and are involved in epistatic interactions that have deleterious impacts upon hybrid fitness in crosses between AB and SD populations of T. californicus. The deleterious interactions for these crosses appear to be largely nuclear/nuclear in nature despite the extensive interactions of these proteins with mtDNA-encoded proteins (the limited evidence for mitochondrial/nuclear coadaptation for these markers in this cross is discussed later). For each of these markers the heterozygote/homozygote genotypic classes are underrepresented in comparison to numbers expected with no epistasis (Figures 2–4GoGo). This congruence in deviations from independence across reciprocal crosses is a strong argument for the pattern's significance despite the lack of replication of each reciprocal cross individually in the experimental design (all F2 adults were pooled to maximize sample size). These interactions between the genes encoding CYC, CYC1, and RISP could then be an example of both partners of a D–M incompatibility leading to lowered hybrid fitness. Below, I explore other potential explanations for this pattern and argue that they are unlikely to be leading to the observed pattern.

If any of these three genes are linked to another, this could generate a significant pattern of nonindependence; however, other work examining patterns of segregation in backcrosses with these markers has demonstrated that they are on three separate chromosomes (H. B. LEISY and C. S. WILLETT, unpublished results). The lack of significant deviations from independence in two-way comparisons for first-stage nauplii also demonstrates that these genes are not linked to one another. Another possibility is that genes linked to these three markers are involved in the interactions rather than the three ETS genes. Given the close nature of interactions between each of these proteins in complex III they are clearly strong candidates for causing the observed epistatic interactions.

An examination of background levels of epistasis between genes that are not functionally interacting would provide a measure of the uniqueness of the observed interactions. Epistatic interactions approaching significance were much less common for other two-way comparisons using these markers and the two ME homologs or histone H1 (WILLETT and BURTON 2001) or for the ME homologs to each other (J. BERKOWITZ and C. S. WILLETT, unpublished results). Although malic enzyme is not directly in the ETS, most animals have an isoform of this protein that functions in the mitochondrion (KOCHAN et al. 1995; DOLEZAL et al. 2004) and could lead to functional interactions with complex III. HARRISON and EDMANDS (2006) have shown that a combination of microsatellite and RAPD markers on different chromosomes only rarely shows significant pairwise epistatic interactions in F2 backcrosses between SD and a population closely related to AB (Royal Palms). Only 6 of 132 possible pairwise comparisons were significant in their data set (only 2 after Bonferroni correction); however, measures across all chromosomes for this cross in their study do suggest the presence of more complex interactions. Combined, these results suggest that the three nuclear-encoded ETS proteins show many more significant pairwise epistatic interactions than other non-ETS genes and that high numbers of significant epistatic interactions are not ubiquitous across markers in crosses of these copepod populations.

Another argument for the epistatic interactions being caused by the deleterious interaction of the protein products of the CYC, CYC1, and RISP genes is that the magnitude of interactions is related to the degree of structural and functional interactions within complex III (and to CYC). The soluble CYC protein is bound to complex III at the CYC1 protein during the electron-transfer phase from complex III to CYC (LANGE and HUNTE 2002). In comparison with the 26-amino-acid differences between the AB and SD CYTB proteins, there is only a single (CYC1) or two (RISP, CYC) amino acid differences between these three nuclear-encoded proteins. Although some of the amino acid changes between T. californicus populations in these proteins fall near the interacting regions, none of them are the interacting residues themselves (WILLETT and BURTON 2004). If, as this study suggests, these changes have important fitness consequences, this would imply that the basic interactions such as the residues responsible for binding have remained unchanged while the residues that have changed will cause smaller modifications to the complex III enzyme. The RISP protein interacts both in structure and in function with the CYC1 and CYTB proteins but has no direct interactions with CYC (IWATA et al. 1998). For two-way comparisons between these three genes, the comparisons between RISP and CYC had smaller deviations from independence as would be predicted by the more limited structural/functional interactions. If we use yeast or mammal complex III as a model there are either six or seven (respectively) other nuclear-encoded subunits in complex III. Of these subunits, QCR6p (corresponding to subunit 8 in mammals) and QCR9p (corresponding to subunit 10) are close in the structure to the active sites of CYC1 and RISP (IWATA et al. 1998; ZHANG et al. 1998; LANGE and HUNTE 2002). Although complex III can function without these two subunits, it appears that they are important for complex assembly and may modify functional interactions as well (KIM et al. 1987; SCHMITT and TRUMPOWER 1991; ZARA et al. 2004). To gain a complete understanding of the evolution of interactions in this complex for these hybrid copepods, these additional subunits should be characterized as well.

For the F2 adult hybrids of the AB and SD populations, dominance appears to play a large role in structuring the interactions between these genes. The two-way genotypes with homozygous/heterozygous combinations are all underrepresented for comparisons between the three ETS gene markers (Figures 2–4GoGo). This pattern suggests that deleterious interactions involving components of complex III may involve one heterozygous locus and one homozygous locus. In models of postzygotic isolation the ratio of heterozygous/homozygous to homozygous/homozygous incompatibilities influences the ability of reproductive isolation to evolve and can help explain the phenomenon known as the large X effect (TURELLI and ORR 2000). They conclude that a small ratio of heterozygote/homozygote incompatibilities (h1 in their terminology) to homozygote/homozygote incompatibilities (h2) would make dominance a good explanation for the observation of a large X effect for differences between species. T. californicus does not have differentiated sex chromosomes so this large X effect is not applicable to this species, but the results for epistasis at complex III could indicate that heterozygous/homozygous deleterious incompatibilities can be uncoupled from homozygous/homozygous interaction effects. In results seen to date from Drosophila and wasp studies h1 is generally much lower than h2 (BREEUWER and WERREN 1995; HOLLOCHER and WU 1996; TRUE et al. 1996; TURELLI and ORR 2000) although the results of TAO et al. (2003), which looked at hybrid male sterility in Drosophila simulans and D. mauritiana using introgressions, suggest that dominance/epistatic relationships are very complex for loci influencing this trait.

It is not clear what the proximate cause of high frequency of heterozygote/homozygote deleterious interactions would be for the genes in this study. Complex III is a stable dimer with each monomer containing one each of the CYC1, RISP, and CYTB proteins (SARASTE 1999) (see Figure 1). There is the potential for both parental alleles to be present in the same complex III dimer, perhaps contributing to the observed interactions. It is also possible that the observed patterns of genotypic deviations derive from the influence of genes linked to the ETS markers combined with the interactions of the ETS proteins themselves. For this study F2 hybrids were examined, which allows for only one generation of recombination; therefore, there will be fairly large numbers of linked genes with the potential to influence the viability of the ETS markers scored in this study. Future studies are planned to isolate these three ETS markers from linked genes, which may help to clear up the nature of the interactions between these genes from the effects of background loci.

The homozygous–homozygous classes do not always deviate in a manner predicted by the D–M incompatibility model; in several cases the nonparental-type double-homozygote class (e.g., RISPSD/SD/CYC1AB/AB in Figure 3) is overrepresented. Clearly the nuclear/nuclear interactions between the three ETS proteins (or linked genes) have a complex relationship to fitness in the hybrid copepods of this population cross. If the observed epistasis is caused solely by the interactions of the ETS proteins, these results would suggest that both positive and negative fitness interactions can result even from single pairs of genes. Whether these loci with both negative and positive epistatic interactions will function as D–M incompatibilities overall will depend upon the relative strength of the effects and whether the negative outweighs the positive. The results of the epistatic interactions observed in this study suggest that D–M incompatibilities are very complex even when broken down to the level of individual loci.

Mitochondrial/nuclear coadaptation:
For the crosses between copepods from the AB and SD populations, coadaptation between the nuclear-encoded CYC1 and RISP genes and the mtDNA-encoded proteins does not appear to play a significant role in determining the fitness of F2 hybrids. RISP and CYC1 are part of the same complex with CYTB (Figure 1). CYTB is evolving more rapidly than the nuclear-encoded RISP and CYC1 proteins with 26-amino-acid differences between the AB and SD proteins (WILLETT and BURTON 2004). For complex III the expectation from structural and functional considerations would be that RISP would have the closest interactions with CYTB (IWATA et al. 1998); however, the deviations seen in the F2 hybrids are not generally consistent with the deviations expected from mtDNA/nuclear coadaptation for this cross (Table 2). It is possible that such coadaptation could be obscured if the neighboring regions of the genome contain opposing interactions that are stronger in effect. These interactions may be detectable in advanced generation hybrids if they exist. Functional studies on complex III such as those performed by RAWSON and BURTON (2002) for CYC and complex IV could also help illuminate any potential coadaptation.

For the CYC gene, only the ABm x SDf cross for CYC has a lowered viability of the non-mtDNA type homozygous class for adults (SD mtDNA-type in this cross) and a nonsignificant difference between the sexes (Tables 1 and 2). The reciprocal cross does not lend additional support to the suggestion that CYC is involved in mitochondrial/nuclear coadaptation in hybrids of these two populations, however, as the same direction of bias toward SD homozygotes is observed in the two reciprocal crosses. The CYC protein interacts with complex IV (cytochrome c reductase), a complex primarily composed of three mtDNA-encoded proteins (SARASTE 1999), in addition to the close interaction of CYC with complex III. In a cross between AB and a copepod population from SC, CYC did show fitness deviations consistent with mitochondrial/nuclear coadaptation in previous studies (WILLETT and BURTON 2001, 2003). A physiological study has also suggested the possibility of coadaptation for CYC and complex IV enzymatic rates for the AB and SC populations (EDMANDS and BURTON 1999). Combined, these studies suggest that if mitochondrial/nuclear coadaptation is contributing to hybrid breakdown in survival in these copepod hybrids it may involve interactions between CYC and complex IV. In other systems, few studies to date have looked at the contribution of individual ETS proteins to organismal fitness due to genomic coadaptation. A number of studies have used crosses between populations or species to examine the magnitude of mitochondrial/nuclear coadaptation on fitness (CLARK and LYCKEGAARD 1988; HUTTER and RAND 1995; RAND et al. 2001, 2006; JAMES and BALLARD 2003; PERROT-MINNOT et al. 2004). Significant impacts of mitochondrial/nuclear coadaptation upon fitness have been found in these studies, demonstrating that such interactions might not be uncommon in diverging populations or species (along with significant effects of nuclear/nuclear interactions in some studies as well). The specific genetic interactions leading to these mitochondrial/nuclear interactions remain largely unknown at this point.

The results obtained in this study can be compared to results for CYC obtained from a previous study using the AB and SD populations (WILLETT and BURTON 2001), and this comparison shows similar overall trends but a few larger differences as well (Table 2). The direction of bias between the AB/AB and SD/SD alternate homozygote classes in this 2001 experiment is qualitatively similar to the results obtained in the present study; the largest discrepancy for CYC is for adults in the ABf x SDm cross that show little difference between the viabilities of the alternate homozygotes in the current results but large differences in opposite directions between males and females in the 2001 results. This variation in results suggests that there can be differences from cross to cross in the magnitude and expression of hybrid incompatibilities even when conducted under similar conditions (the largest controlled difference was the use of artificial seawater for the current experiment vs. natural seawater for the previous experiment). It is conceivable that other experimental factors could influence the repeatability of these fitness effects; in particular, relative levels of competition (different numbers of nauplii will be in each of the F2 petri dishes) and algal growth in the cultures are factors that vary to some extent but not in a systematic way between experiments. It is also possible that the differences between studies result from genetic variation segregating within each of the populations. The copepods used in these two sets of crosses were from outbred lab populations sampled at different times from the field. It should be noted, however, that the magnitude of genetic variation within each population is much lower than the divergence seen between each of these three populations (typically around an order of magnitude lower for nuclear genes if variation at silent sites is compared within and between these two populations; WILLETT and BURTON 2004). In addition, 30 females and 30 males were used to set up each cross, which should have the effect of averaging over any genetic variation present within each of the populations. Clearly, repeatability is an issue for interpreting the single-locus effects but could potentially reflect an underlying biological reality that many factors can influence the relative fitness of copepods, and at these relatively low levels of reproductive isolation, specific deleterious interactions might be somewhat transient in nature.

In support of the changeable nature of genetic interactions, the environment appeared to play a significant role in the expression of potential incompatibilities in a previous study of interpopulation hybrids of T. californicus. WILLETT and BURTON (2003) found that temperature environment affected whether there was selection at CYC in a cross between AB and Santa Cruz. Temperature interactions would not be surprising given that these copepods are ectotherms and the ETS should have evolved to function best under the range of conditions that particular populations of copepods have experienced in their past. The stability of protein–protein interactions is often highly dependent upon temperature and may then lead to the widespread involvement of temperature in the expression of D–M incompatibilities. Replication of the current study under a different temperature environment could help reveal if nuclear/nuclear interactions detected in this study for ETS genes are influenced as dramatically by the environmental conditions.


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
LITERATURE CITED
 
I thank E. Hoddeson, Q. Qin, and J. Berkowitz for their help with the copepod crosses and sample preparation. Helpful suggestions on the manuscript were provided by R. Burton, M. Servedio, T. Vision, A. Bouck, and S. Hartman; by the University of North Carolina (UNC) evolgen reading group; and by two anonymous reviewers. I thank M. Servedio for help with analyses. C.W. was supported by money from UNC and National Science Foundation grant DEB-0516139.


   LITERATURE CITED
TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 

ARMBRUSTER, P., W. E. BRADSHAW and C. M. HOLZAPPEL, 1997 Evolution of the genetic architecture underlying fitness in the pitcher-plant mosquito, Wyeomyia smithii. Evolution 51: 451–458.[CrossRef]

BARBASH, D. A., D. F. SIINO, A. M. TARONE and J. ROOTE, 2003 A rapidly evolving MYB-related protein causes species isolation in Drosophila. Proc. Natl. Acad. Sci. USA 100: 5302–5307.[Abstract/Free Full Text]

BARBASH, D. A., P. AWADALLA and A. M. TARONE, 2004 Functional divergence caused by ancient positive selection of a Drosophila hybrid incompatibility locus. PloS Biol. 2: 839–848.

BREEUWER, J. A., and J. H. WERREN, 1995 Hybrid breakdown between two haplodiploid species: the role of nuclear and cytoplasmic genes. Evolution 49: 705–717.[CrossRef]

BURKE, J. M., T. J. VOSS and M. L. ARNOLD, 1998 Genetic interactions and natural selection in Louisiana iris hybrids. Evolution 52: 1304–1310.[CrossRef]

BURTON, R. S., 1985 Mating system of the intertidal copepod Tigriopus californicus. Mar. Biol. 86: 247–252.[CrossRef]

BURTON, R. S., 1986 Evolutionary consequences of restricted gene flow in the intertidal copepod Tigriopus californicus. Bull. Mar. Sci. 39: 526–535.

BURTON, R. S., 1987 Differentiation and integration of the genome in populations of Tigriopus californicus. Evolution 41: 504–513.[CrossRef]

BURTON, R. S., 1990 Hybrid breakdown in developmental time in the copepod Tigriopus californicus. Evolution 44: 1814–1822.[CrossRef]

BURTON, R. S., 1998 Intraspecific phylogeography across the Point Conception biogeographic boundary. Evolution 52: 734–745.[CrossRef]

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: 5197–5201.[Abstract/Free Full Text]

BURTON, R. S., P. D. RAWSON and S. EDMANDS, 1999 Genetic architecture of physiological phenotypes: empirical evidence for coadapted gene complexes. Am. Zool. 39: 451–462.

CLARK, A. G., and E. M. S. LYCKEGAARD, 1988 Natural selection with nuclear and cytoplasmic transmission. III. Joint analysis of segregation and mtDNA in Drosophila melanogaster. Genetics 118: 471–481.[Abstract/Free Full Text]

COYNE, J. A., S. SIMEONIDIS and P. ROONEY, 1998 Relative paucity of genes causing inviability in hybrids between Drosophila melanogaster and D. simulans. Genetics 150: 1091–1103.[Abstract/Free Full Text]

DOBZHANSKY, T., 1936 Studies on hybrid sterility. II. Localization of sterility factors in Drosophila pseudoobscura hybrids. Genetics 21: 113–135.[Free Full Text]

DOBZHANSKY, TH., 1950 Genetics of natural populations. XIX. Origin of heterosis through natural selection in populations of Drosophila pseudoobscura. Genetics 35: 288–302.[Free Full Text]

DOLEZAL, P., S. VANACOVA, J. TACHEZY and I. HRDY, 2004 Malic enzymes of Trichomonas vaginalis: two enzyme families, two distinct origins. Gene 329: 81–92.[CrossRef][Medline]

EDMANDS, S., 1999 Heterosis and outbreeding depression in interpopulation crosses spanning a wide range of divergence. Evolution 53: 1757–1768.[CrossRef]

EDMANDS, S., 2001 Phylogeography of the intertidal copepod Tigriopus californicus reveals substantially reduced population differentiation at northern latitudes. Mol. Ecol. 10: 1743–1750.[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: 1972–1978.[CrossRef]

EDMANDS, S., and J. K. DEIMLER, 2004 Local adaptation, intrinsic coadaptation and the effects of environmental stress on interpopulation hybrids in the copepod Tigriopus californicus. J. Exp. Mar. Biol. Ecol. 303: 183–196.[CrossRef]

FENSTER, C. B., and L. F. GALLOWAY, 2000a Inbreeding and outbreeding depression in natural populations of Chamaecrista fasciculata (Fabaceae). Conserv. Biol. 14: 1406–1412.[CrossRef]

FENSTER, C. B., and L. F. GALLOWAY, 2000b Population differentiation in an annual legume: genetic architecture. Evolution 54: 1157–1172.[CrossRef][Medline]

FISHMAN, L., and J. H. WILLIS, 2005 A novel meiotic drive locus completely distorts segregation in Mimulus (monkeyflower) hybrids. Genetics 169: 347–353.[Abstract/Free Full Text]

FISHMAN, L., A. KELLY, E. MORGAN and J. H. WILLIS, 2001 A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions. Genetics 159: 1701–1716.[Abstract/Free Full Text]

GALLOWAY, L. F., and C. B. FENSTER, 1999 The role of nuclear and cytoplasmic factors in the adaptive evolution of populations of Chamaecrista fasciculata (Fabaceae). Evolution 53: 1734–1743.[CrossRef]

HALDANE, J. B. S., 1956 The estimation of viabilities. J. Genet. 59: 29–36.

HALL, M. C., and J. H. WILLIS, 2005 Transmission ratio distortion in intraspecific hybrids of Mimulus guttatus: implications for genomic divergence. Genetics 170: 375–386.[Abstract/Free Full Text]

HARRISON, J. S., and S. EDMANDS, 2006 Chromosomal basis of viability differences in Tigriopus californicus interpopulation hybrids. J. Evol. Biol. (in press).

HOELZEL, A. R., and A. GREEN, 1992 Analysis of population-level variation by sequencing PCR-amplified DNA, pp. 159–187 in Practical Approach Series: Molecular Genetic Analysis of Populations, edited by A. R. HOELZEL. Oxford University Press, New York.

HOLLOCHER, H., and C.-I. WU, 1996 The genetics of reproductive isolation in the Drosophila simulans clade: X vs. autosomal effects and male vs. female effects. Genetics 143: 1243–1255.

HUTTER, C. M., and D. M. RAND, 1995 Competition between mitochondrial haplotypes in distinct nuclear genetic environments: Drosophila pseudoobscura vs. D. persimilis. Genetics 140: 537–548.[Abstract]

IWATA, S., J. L. LEE, K. OKADA, J. K. LEE, M. IWATA et al., 1998 Complete structure of the 11-subunit bovine mitochondrial cytochrome bc1 complex. Science 281: 64–71.[Abstract/Free Full Text]

JAMES, A. C., and J. W. O. BALLARD, 2003 Mitochondrial genotype affects fitness in Drosophila simulans. Genetics 164: 187–194.[Abstract/Free Full Text]

KIM, C. H., C. BALNY and T. E. KING, 1987 Role of the hinge protein in the electron transfer between cardiac cytochrome c1 and c. J. Biol. Chem. 262: 8103–8108.[Abstract/Free Full Text]

KOCHAN, Z., J. KARBOWSKA, F. BUKATO and J. SWIERCZYNSKI, 1995 Comparative studies on NADP+-linked malic enzyme in the central nervous system of ectothermic and endothermic animals. Comp. Biochem. Physiol. 110B: 309–314.[CrossRef]

LANGE, C., and C. HUNTE, 2002 Crystal structure of the yeast cytochrome bc1 complex with its bound substrate cytochrome c. Proc. Natl. Acad. Sci. USA 99: 2800–2805.[Abstract/Free Full Text]

LI, Z., S. R. M. PINSON, A. H. PATERSON, W. D. PARK and J. W. STANSEL, 1997 Genetics of hybrid sterility and hybrid breakdown in an intrasubspecific rice (Oryza sativa L.) population. Genetics 145: 1139–1148.[Abstract]

MULLER, H. J., 1942 Isolating mechanisms, evolution, and temperature. Biol. Symp. 6: 71–125.

PALOPOLI, M. F., and C.-I WU, 1994 Genetics of hybrid male sterility between Drosophila sibling species: a complex web of epistasis is revealed in interspecific studies. Genetics 138: 329–341.[Abstract]

PERROT-MINNOT, M.-J., A. MIGEON and M. NAVAJAS, 2004 Intergenomic interactions affect female reproduction: evidence from introgression and inbreeding depression in a haplodiploid mite. Heredity 93: 551–558.[CrossRef][Medline]

PRESGRAVES, D. C., 2003 A fine-scale genetic analysis of hybrid incompatibilities in Drosophila. Genetics 163: 955–972.[Abstract/Free Full Text]

PRESGRAVES, D. C., L. BALAGOPALAN, S. M. ABMAYR and H. A. ORR, 2003 Adaptive evolution drives divergence of a hybrid inviability gene between two species of Drosophila. Nature 423: 715–719.[CrossRef][Medline]

RAND, D. M., A. G. CLARK and L. M. KANN, 2001 Sexually antagonistic cytonuclear fitness interactions in Drosophila melanogaster. Genetics 159: 173–187.[Abstract/Free Full Text]

RAND, D. M., A. FRY and L. SHELDAHL, 2006 Nuclear-mitochondrial epistasis and Drosophila aging: introgression of Drosophila simulans mtDNA modifies longevity in D. melanogaster nuclear backgrounds. Genetics 172: 329–341.[Abstract/Free Full Text]

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: 12955–12958.[Abstract/Free Full Text]

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: 15–22.[CrossRef][Medline]

RIESEBERG, L. H., B. SINERVO, C. R. LINDER, M. C. UNGERER and D. M. ARIAS, 1996 Role of gene interactions in hybrid speciation: evidence from ancient and experimental hybrids. Science 272: 741–745.[Abstract]

SARASTE, M., 1999 Oxidative phosphorylation at the fin de siecle. Science 283: 1488–1493.[Abstract/Free Full Text]

SCHMITT, M. E., and B. L. TRUMPOWER, 1991 The petite phenotype resulting from a truncated copy of subunit 6 results from loss of assembly of the cytochrome bc1 complex and can be suppressed by overexpression of subunit 9. J. Biol. Chem. 266: 14958–14963.[Abstract/Free Full Text]

TAO, Y., Z-B. ZENG, J. LI, D. L. HARTL and C. C. LAURIE, 2003 Genetic dissection of hybrid incompatibilities between Drosophila simulans and D. mauritiana. II. Mapping hybrid male sterility loci on the third chromosome. Genetics 164: 1399–1418.[Abstract/Free Full Text]

TING, C.-T., S.-C. TSAUR, M.-L. WU and C.-I WU, 1998 A rapidly evolving homeobox at the site of a hybrid sterility gene. Science 282: 1501–1504.[Abstract/Free Full Text]

TRUE, J. R., B. S. WEIR and C. C. LAURIE, 1996 A genome-wide survey of hybrid incompatibility factors by the introgression of marked segments of Drosophila mauritiana chromosomes into D. simulans. Genetics 142: 819–837.[Abstract]

TURELLI, M., and H. A. ORR, 2000 Dominance, epistasis, and the genetics of postzygotic isolation. Genetics 154: 1663–1679.[Abstract/Free Full Text]

WHITLOCK, M. C., and P. C. PHILLIPS, 1995 Multiple fitness peaks and epistasis. Annu. Rev. Ecol. Syst. 26: 601–629.[CrossRef]

WILLETT, C. S., and R. S. BURTON, 2001 Viability of cytochrome c depends on cytoplasmic background in Tigriopus californicus. Evolution 55: 1592–1599.[CrossRef][Medline]

WILLETT, C. S., and R. S. BURTON, 2003 Environmental influences on epistatic interactions: viabilities of cytochrome c genotypes in interpopulation crosses. Evolution 57: 2286–2292.[CrossRef][Medline]

WILLETT, C. S., and R. S. BURTON, 2004 Evolution of interacting proteins in the mitochondrial electron transport system in a marine copepod. Mol. Biol. Evol. 21: 443–453.[Abstract/Free Full Text]

ZARA, V., I. PALMISANO, L. CONTE and B. L. TRUMPOWER, 2004 Further insights into the assembly of the yeast cytochrome bc1 complex based on analysis of single and double mutants lacking supernumerary subunits and cytochrome b. Eur. J. Biochem. 271: 1209–1218.[Medline]

ZHANG, Z., J. HUANG, V. M. SCULMEISTER, Y.-I. CHI, K. K. KIM et al., 1998 Electron transfer by domain movement in cytochrome bc1. Nature 392: 677–684.[CrossRef][Medline]

Communicating editor: D. M. RAND




This article has been cited by other articles:


Home page
 GeneticsHome page
W. J. Dickinson
Synergistic Fitness Interactions and a High Frequency of Beneficial Changes Among Mutations Accumulated Under Relaxed Selection in Saccharomyces cerevisiae
Genetics, March 1, 2008; 178(3): 1571 - 1578.
[Abstract] [Full Text] [PDF]

Home page
 J HeredHome page
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]

Home page
 GeneticsHome page
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]

Home page
 GeneticsHome page
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]