MALADAPTED GENE COMPLEXES WITHIN POPULATIONS OF THE INTERTIDAL COPEPOD TIGRIOPUS CALIFORNICUS?

The prevalence of F2 hybrid breakdown in interpopulation crosses of the marine copepod Tigriopus californicus can be explained by disruption of coadapted gene complexes. This study further dissects the nature of hybrid gene interactions, revealing that parental populations may also harbor maladapted gene complexes. Diagnostic molecular markers (14) were assayed in reciprocal F2 hybrids to test for gene interactions affecting viability. Results showed some evidence of nuclear–nuclear coadaptation. Although there were no significant examples of pairwise linkage disequilibrium between physically unlinked loci, one of the two reciprocal crosses did show an overall excess of parental double homozygotes and an overall dearth of nonparental double homozygotes. In contrast, the nuclear–cytoplasmic data showed a stronger tendency toward maladaptation within the specific inbred lines used in this study. For three out of four loci with significant frequency differences between reciprocal F2, homozygotes were favored on the wrong cytoplasmic background. A separate study of reciprocal backcross hybrids between the same two populations (but different inbred lines) revealed faster development time when the full haploid nuclear genome did not match the cytoplasm. The occurrence of such suboptimal gene complexes may be attributable to effects of genetic drift in small, isolated populations.     

Interpopulation hybrid breakdown maps to the mitochondrial genome

Hybrid breakdown, or outbreeding depression, is the loss of fitness observed in crosses between genetically divergent populations. The role of maternally inherited mitochondrial genomes in hybrid breakdown has not been widely examined. Using laboratory crosses of the marine copepod Tigriopus californicus, we report that the low fitness of F(3) hybrids is completely restored in the offspring of maternal backcrosses, where parental mitochondrial and nuclear genomic combinations are reassembled. Paternal backcrosses, which result in mismatched mitochondrial and nuclear genomes, fail to restore hybrid fitness. These results suggest that fitness loss in T. californicus hybrids is completely attributable to nuclear-mitochondrial genomic interactions. Analyses of ATP synthetic capacity in isolated mitochondria from hybrid and backcross animals found that reduced ATP synthesis in hybrids was also largely restored in backcrosses, again with maternal backcrosses outperforming paternal backcrosses. The strong fitness consequences of nuclear-mitochondrial interactions have important, and often overlooked, implications for evolutionary and conservation biology.     

Deleterious epistatic interactions between electron transport system protein-coding loci in the copepod Tigriopus californicus

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.     

DNA sequences and cross-breeding experiments in the hawthorn spider mite Amphitetranychus viennensis reveal high genetic differentiation between Japanese and French populations

Sequence variation of the complete second internal transcribed spacer (ITS2, 445 bp) of nuclear ribosomal DNA and part of the mitochondrial cytochrome oxidase I gene (COI, 350 bp) was examined in Amphitetranychus viennensis (Zacher) mites (Acari:Tetranychidae) from four French and four Japanese locations. Sequence analysis consistently revealed the separation of the samples in two major groups: French mites differed from Japanese by 3.8–4.1% of the nucleotide divergence in COI sequences. These two groups also displayed distinct ITS2 consensus sequences (2.1% nucleotide divergence). A few variations, not affecting the diagnostic sites around the consensus sequence, were revealed among cloned copies of the same individual. Reciprocal crosses and backcrosses between one French and two Japanese populations disclosed strong reproductive incompatibility. However, fertile hybrid females were obtained, indicating the conspecificity of the tested mites. Despite the presence of Wolbachia in the French strain, but not in the Japanese ones, our crosses did not display the unidirectional incompatibility typically produced by this microorganism, but rather a bidirectional – although asymmetrical – incompatibility pattern. The post-zygotic incompatibilities in A. viennensis cannot be explained by the presence of Wolbachia but to some extent by mite genome divergence resulting from limited gene exchange between allopatric populations. Experiments of Wolbachia elimination by antibiotic treatment and subsequent crosses with cured strains are still needed to fully understand the reproductive incompatibility patterns in this mite species.     

HETEROSIS AND OUTBREEDING DEPRESSION IN INTERPOPULATION CROSSES SPANNING A WIDE RANGE OF DIVERGENCE

The intertidal copepod Tigriopus californicus was used as a model organism to look at effects of crossing distance on fitness and to investigate the genetic mechanisms responsible. Crosses were conducted between 12 pairs of populations spanning a broad range of both geographic distance (5 m to 2007 km) and genetic distance (0.2% to 22.3% sequence divergence for a 606-bp segment of the mitochondrial COI gene). For each pair of populations, three fitness components (hatching number, survivorship number, and metamorphosis number) were measured in up to 16 cohorts including parentals, reciprocal F₁, F₂, F₃, and first-generation backcross hybrids. Comparisons of each set of cohorts allowed estimation of within- and between-locus gene interaction. Relative to parentals, F₁ hybrids showed a trend toward increased fitness, with no correspondence with population divergence, and a decrease in variance, which in some cases correlated with population divergence. In sharp contrast, F₂ hybrids had a decrease in fitness and an increase in variance that both corresponded to population divergence. Genetic interpretation of these patterns suggests that both the beneficial effects of dominance and the detrimental effects of breaking up coadaptation are magnified by increasing evolutionary distance between populations. Because there is no recombination in T. californicus females, effects of recombination can be assessed by comparing F₁ hybrid males and females backcrossed to parentals. Both recombinant and nonrecombinant backcross hybrids showed a decline in fitness correlated with population divergence, indicating that segregation among chromosomes contributes to the breakup of coadaptation. Although there was no difference in mean fitness between the two backcross types, recombinational backcrosses showed greater variance for fitness than nonrecombinational backcrosses, suggesting that the breakup of parental gene ombinations within chromosomes has both beneficial and detrimental effects.     

Viability effects and not meoitic drive cause dramatic departures from Mendelian inheritance for malic enzyme in hybrids of Tigriopus californicus populations

The genetic basis of post-zygotic reproductive isolation is beginning to be untangled in closely related species, but less is known about the genetics of reproductive isolation between divergent populations. Here, two genes encoding malic enzyme (ME) are isolated from the copepod Tigriopus californicus and their influence upon lowered viability in F(2) hybrids of genetically divergent populations is determined. Each ME gene has diverged extensively between T. californicus populations and one gene shows evidence for a recent selective sweep. Segregation patterns of genotypes for both ME genes in adult F(2) hybrids reveal dramatic departures from Mendelian inheritance, deviations that are not seen in F(2) nauplii implying that selection is acting during development based upon the genotype at these ME genes. These results imply that selection against deleterious gene combinations and not aberrant segregation (i.e. meiotic drive) is likely to lead to dramatic departures from Mendelian inheritance observed in these crosses.     

No evidence for faster male hybrid sterility in population crosses of an intertidal copepod (Tigriopus californicus)

Two different forces are thought to contribute to the rapid accumulation of hybrid male sterility that has been observed in many inter-specific crosses, namely the faster male and the dominance theories. For male heterogametic taxa, both faster male and dominance would work in the same direction to cause the rapid evolution of male sterility; however, for taxa lacking differentiated sex chromosomes only the faster male theory would explain the rapid evolution of male hybrid sterility. It is currently unknown what causes the faster evolution of male sterility, but increased sexual selection on males and the sensitivity of genes involved in male reproduction are two hypotheses that could explain the observation. Here, patterns of hybrid sterility in crosses of genetically divergent copepod populations are examined to test potential mechanisms of faster male evolution. The study species, Tigriopus californicus, lacks differentiated, hemizygous sex chromosomes and appears to have low levels of divergence caused by sexual selection acting upon males. Hybrid sterility does not accumulate more rapidly in males than females in these crosses suggesting that in this taxon male reproductive genes are not inherently more prone to disruption in hybrids.     

Genetic consequences of many generations of hybridization between divergent copepod populations

Crosses between populations of the copepod Tigriopus californicus typically result in outbreeding depression. In this study, replicate hybrid populations were initiated with first generation backcross hybrids between two genetically distinct populations from California: Royal Palms (RP) and San Diego (SD). Reciprocal F(1) were backcrossed to SD, resulting in expected starting frequencies of 25% RP/75% SD nuclear genes on either a pure RP cytoplasmic or a pure SD cytoplasmic background. After 1 year of hybridization (up to 15 generations), seven microsatellite loci were scored in two replicates on each cytoplasmic background. Frequencies of the rarer RP alleles increased significantly in all four replicates, regardless of cytoplasmic source, producing a mean hybridity of 0.97 (maximum = 1), instead of the expected 0.50. Explicit tests for heterozygote excess across loci and replicates showed significant deviations. Only the two physically linked markers showed linkage disequilibrium in all replicates. Subsequent fitness assays in parental populations and early generation hybrids revealed lower fitness in RP than SD, and significant F(2) breakdown. Computer simulations showed that selection must be invoked to explain the shift in allele frequencies. Together, these results suggest that hybrid inferiority in early generations gave way to hybrid superiority in later generations.     

Hybrid breakdown weakens under thermal stress in population crosses of the copepod Tigriopus californicus

The outcome of hybridization can be impacted by environmental conditions, which themselves can contribute to reproductive isolation between taxa. In crosses of genetically divergent populations, hybridization can have both negative and positive impacts on fitness, the balance between which might be tipped by changes in the environment. Genetically divergent populations of the intertidal copepod Tigriopus californicus have been shown to differ in thermal tolerance at high temperatures along a latitudinal gradient. In this study, a series of crosses were made between pairs of genetically divergent populations of T. californicus, and the thermal tolerance of these hybrids was tested. In most cases, the first-generation hybrids had relatively high thermal tolerance and the second-generation hybrids were not generally reduced below the less-tolerant parental population for high temperature tolerance. This pattern contrasts with previous studies from crosses of genetically divergent populations of this copepod, which often shows hybrid breakdown in these second-generation hybrids for other measures of fitness. These results suggest that high temperature stress could either increase the positive impacts of hybridization or decrease the negative impacts of hybridization resulting in lowered hybrid breakdown in these population crosses.     

Fitness and morphological outcomes of many generations of hybridization in the copepod Tigriopus californicus

Hybridization between genetically divergent populations is an important evolutionary process, with an outcome that is difficult to predict. We used controlled crosses and freely mating hybrid swarms, followed for up to 30 generations, to examine the morphological and fitness consequences of interpopulation hybridization in the copepod Tigriopus californicus. Patterns of fitness in two generations of controlled crosses were partly predictive of long‐term trajectories in hybrid swarms. For one pair of populations, controlled crosses revealed neutral or beneficial effects of hybridization after the F1 generation, and hybrid swarm fitness almost always equalled or exceeded that of the midparent. For a second pair, controlled crosses showed F2 hybrid breakdown, but increased fitness in backcrosses, and hybrid swarm fitness deviated both above and below that of the parentals. Nevertheless, individual swarm replicates exhibited different fitness trajectories over time that were not related in a simple manner to their hybrid genetic composition, and fixation of fitter hybrid phenotypes was not observed. Hybridization did not increase overall morphological variation, and underlying genetic changes may have been masked by phenotypic plasticity. Nevertheless, one type of hybrid swarm exhibited a repeatable pattern of transgressively large eggsacs, indicating a positive effect of hybridization on individual fecundity. Additionally, both parental and hybrid swarms exhibited common phenotypic trends over time, indicating common selective pressures in the laboratory environment. Our results suggest that, in a system where much work has focused on F2 hybrid breakdown, the long‐term fitness consequences of interpopulation hybridization are surprisingly benign.     

Hybrid Breakdown in Cichlid Fish

Studies from a wide diversity of taxa have shown a negative relationship between genetic compatibility and the divergence time of hybridizing genomes. Theory predicts the main breakdown of fitness to happen after the F1 hybrid generation, when heterosis subsides and recessive allelic (Dobzhansky-Muller) incompatibilities are increasingly unmasked. We measured the fitness of F2 hybrids of African haplochromine cichlid fish bred from species pairs spanning several thousand to several million years divergence time. F2 hybrids consistently showed the lowest viability compared to F1 hybrids and non-hybrid crosses (crosses within the grandparental species), in agreement with hybrid breakdown. Especially the short- and long-term survival (2 weeks to 6 months) of F2 hybrids was significantly reduced. Overall, F2 hybrids showed a fitness reduction of 21% compared to F1 hybrids, and a reduction of 43% compared to the grandparental, non-hybrid crosses. We further observed a decrease of F2 hybrid viability with the genetic distance between grandparental lineages, suggesting an important role for negative epistatic interactions in cichlid fish postzygotic isolation. The estimated time window for successful production of F2 hybrids resulting from our data is consistent with the estimated divergence time between the multiple ancestral lineages that presumably hybridized in three major adaptive radiations of African cichlids.     

Significant variation for fitness impacts of ETS loci in hybrids between populations of Tigriopus californicus

The connections between the genes that cause hybrid incompatibilities and the physiological processes disrupted in hybrids by these incompatibilities are not well understood. The interactions between proteins in the electron transport system (ETS) in the copepod, Tigriopus californicus, have emerged as a potential model system to explore such connections. In this study, the effects on hybrid fitness of 3 different nuclear loci encoding proteins of the ETS are examined in hybrid copepods obtained from crosses of genetically divergent populations of this species. The potential interactions between these genes and mitochondrial-encoded proteins of the ETS are also explored; these interactions have been shown to have diverged functionally between these populations in other studies. Large deviations from Mendelian inheritance are found in genotypic ratios at each of the 3 loci in adults but not in nauplii, demonstrating genotype-based selection during development. The length of developmental time of hybrids appears to influence the pattern of deviations in these loci, likely in conjunction with levels of competition in these crosses. The major finding of this study is that in repeated crosses, the nature of deviations at these ETS loci shows dramatic differences suggesting that slight perturbations in initial conditions can dramatically shift the patterns of selection at these ETS loci in interpopulation hybrids.     

Genetic studies of two sister species in the Drosophila melanogaster subgroup, D. yakuba and D. santomea

We performed genetic analysis of hybrid sterility and of one morphological difference (sex-comb tooth number) on D. yakuba and D. santomea, the former species widespread in Africa and the latter endemic to the oceanic island of S?o Tomé, on which there is a hybrid zone. The sterility of hybrid males is due to at least three genes on the X chromosome and at least one on the Y, with the cytoplasm and large sections of the autosomes having no effect. F1 hybrid females carrying two X chromosomes from either species are perfectly fertile despite their genetic similarity to completely sterile F1 hybrid males. This implies that the appearance of Haldane's rule in this cross is at least partially due to the faster accumulation of genes causing male than female sterility. The larger effects of the X and Y chromosomes than of the autosomes, however, also suggest that the genes causing male sterility are recessive in hybrids. Some female sterility is also seen in interspecific crosses, but this does not occur between all strains. This is seen in pure-species females inseminated by heterospecific males (probably reflecting incompatibility between the sperm of one species and the female reproductive tract of the other) as well as in inseminated F1 and backcross females, probably reflecting genetically based incompatibilities in hybrids that affect the reproductive system. The latter 'innate' sterility appears to involve deleterious interactions between D. santomea chromosomes and D. yakuba cytoplasm. The difference in male sex-comb tooth number appears to involve fairly large effects of the X chromosome. We discuss the striking evolutionary parallels in the genetic basis of sterility, in the nature of sexual isolation, and in morphological differences between the D. santomea/D. yakuba divergence and two other speciation events in the D. melanogaster subgroup involving island colonization.     

Tracing hybrid incompatibilities to single amino acid substitutions

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.     

THE GENETICS OF CENTRAL AND MARGINAL POPULATIONS OF DROSOPHILA SERRATA. II. HYBRID BREAKDOWN IN FITNESS COMPONENTS AS A CORRELATED RESPONSE TO SELECTION FOR DESICCATION RESISTANCE

Multiple-peak epistasis is one of the four premises that underlie Wright's shifting-balance theory of evolution. A selection experiment was conducted in an attempt to push different geographic populations to different fitness peaks as a correlated response to selection for an additively controlled character (desiccation resistance). Four populations of Drosophila serrata, sampled from central and marginal areas of its distribution along a 3000-km stretch of Australia's east coast, underwent selection for desiccation resistance for 14 generations. After selection had ceased, control lines from each of the populations were crossed to determine the amount of hybrid breakdown that existed before selection and selected lines were crossed to determine the amount of hybrid breakdown after selection. Hybrid breakdown was measured in three fitness traits: developmental time, viability, and fecundity. When the individual crosses were examined, virtually no evidence was found for hybrid breakdown between these populations. However, the level of hybrid breakdown in development time in the control lines increased as the distance between the populations in the field increased. This relationship was lost in the selected lines. Therefore, selection for desiccation resistance influenced the level of hybrid breakdown in a fitness trait, although selection may need to be maintained for longer than 14 generations if a new relationship between hybrid breakdown and distance is to be formed.     

Elevated oxidative damage is correlated with reduced fitness in interpopulation hybrids of a marine copepod

Aerobic energy production occurs via the oxidative phosphorylation pathway (OXPHOS), which is critically dependent on interactions between the 13 mitochondrial DNA (mtDNA)-encoded and approximately 70 nuclear-encoded protein subunits. Disruptive mutations in any component of OXPHOS can result in impaired ATP production and exacerbated oxidative stress; in mammalian systems, such mutations are associated with ageing as well as numerous diseases. Recent studies have suggested that oxidative stress plays a role in fitness trade-offs in life-history evolution and functional ecology. Here, we show that outcrossing between populations with divergent mtDNA can exacerbate cellular oxidative stress in hybrid offspring. In the copepod Tigriopus californicus, we found that hybrids that showed evidence of fitness breakdown (low fecundity) also exhibited elevated levels of oxidative damage to DNA, whereas those with no clear breakdown did not show significantly elevated damage. The extent of oxidative stress in hybrids appears to be dependent on the degree of genetic divergence between their respective parental populations, but this pattern requires further testing using multiple crosses at different levels of divergence. Given previous evidence in T. californicus that hybridization disrupts nuclear/mitochondrial interactions and reduces hybrid fitness, our results suggest that such negative intergenomic epistasis may also increase the production of damaging cellular oxidants; consequently, mtDNA evolution may play a significant role in generating postzygotic isolating barriers among diverging populations.     

Viability of cytochrome c genotypes depends on cytoplasmic backgrounds in Tigriopus californicus

Because of their extensive functional interaction, mitochondrial DNA (mtDNA) and nuclear genes may evolve to form coadapted complexes within reproductively isolated populations. As a consequence of coadaptation, the fitness of particular nuclear alleles may depend on mtDNA genotype. Among populations of the copepod Tigriopus californicus, there are high levels of amino acid substitutions in both the mtDNA genes encoding subunits of cytochrome c oxidase (COX) and the nuclear gene encoding cytochrome c (CYC), the substrate for COX. Because of the functional interaction between enzyme and substrate proteins, we hypothesized that the fitness of CYC genotypes would depend on mtDNA genotype. To test this hypothesis, segregation ratios for CYC and a second nuclear marker (histone H1) unrelated to mitochondrial function were scored in F2 progeny of several reciprocal interpopulation crosses. Genotypic ratios at the CYC locus (but not the H1 locus) differed between reciprocal crosses and differed from expected Mendelian ratios, suggesting that CYC genotypic fitnesses were strongly influenced by cytoplasmic (including mtDNA) background. However, in most cases the nature of the deviations from Mendelian ratios and differences between reciprocal crosses are not consistent with simple coevolution between CYC and mtDNA background. In a cross in which both newly hatched larvae and adults were sampled, only the adult sample showed deviations from Mendelian ratios, indicating that genotypic viabilities differed. In two of six crosses, large genotypic ratio differences for CYC were observed between the sexes. These results suggest that significant variation in nuclear-mtDNA coadaptation may exist between T. californicus populations and that the relative viability of specific cytonuclear allelic combinations is somehow affected by sex.     

Population differentiation in the beetle Tribolium castaneum. II. Haldane'S rule and incipient speciation

The heterogametic sex tends to be rare, absent, sterile, or deformed in F1 hybrid crosses between species, a pattern called Haldane's rule (HR). The introgression of single genes or chromosomal regions from one drosophilid species into the genetic background of another have shown that HR is most often associated with fixed genetic differences in inter-specific crosses. However, because such introgression studies have involved species diverged several hundred thousand generations from a common ancestor, it is not clear whether HR attends the speciation process or results from the accumulation of epistatically acting genes postspeciation. We report the first evidence for HR prior to speciation in crosses between two populations of the red flour beetle, Tribolium castaneum, collected 931 km apart in Colombia and Ecuador. In this cross, HR is manifested as an increase in the proportion of deformed males compared to females and the expression of HR is temperature dependent. Neither population, when crossed to a geographically distant population from Japan, exhibits HR at any rearing temperature. Using joint-scaling analysis and additional data from backcrosses and F2's, we find that the hybrid incompatibilities and the emergence of HR are concurrent processes involving interactions between X-linked and autosomal genes. However, we also find many examples of incompatibilities manifest by F2 and backcross hybrids but not by F1 hybrids and most incompatibilities are not sex different in their effects, even when they involve both X-autosomal interactions and genotype-by-environment interactions. We infer that incipient speciation in flour beetles can occur with or without HR and that significant hybrid incompatibilities result from the accumulation of epistatically acting gene differences between populations without differentially affecting the heterogametic sex in F1 hybrids. The temperature dependence of the incompatibilities supports the inference that genotype-by-environment interactions and adaptation to different environments contribute to the genetic divergence important to postzygotic reproductive isolation.     

THE GENOMIC TRAJECTORY OF HYBRID SWARMS: OUTCOMES OF REPEATED CROSSES BETWEEN POPULATIONS OF TIGRIOPUS CALIFORNICUS

Introgressive hybridization between genetically divergent populations is an important evolutionary process. The degree to which repeated hybridization events between the same parental taxa lead to similar genomic outcomes is unknown. This study addressed this question by following genomic trajectories of replicate hybrid swarms of the copepod Tigriopus californicus over many generations of free mating. Swarm composition was determined both by differential reproductive success of founder individuals and subsequent selection on hybrid genotypes. For one cross, between two populations showing differential fitness in the laboratory and no hybrid breakdown, the genetic trajectory was highly repeatable: replicates rapidly became dominated by alleles from the fitter parent. In a second cross, between two populations showing similar fitness and significant F2 hybrid breakdown, alleles from alternative populations dominated different replicates. Swarms exhibited a general temporal trend of decreasing cytonuclear mismatch. Some patterns of differential introgression across the genome were strikingly congruent amongst swarm replicates, both within and between cross types, and reflected patterns of segregation distortion previously observed within controlled crosses between the same parental populations. Differences in heterozygosity between the sexes, and evidence for a previously suspected sex‐distortion locus, suggest that complex interactions between sex and genotype influence hybrid swarm outcome.     

Population differentiation in the beetleTribolium castaneum. I. Genetic architecture

We used joint-scaling analyses in conjunction with rearing temperature variation to investigate the contributions of additive, non-additive, and environmental effects to genetic divergence and incipient speciation among 12 populations of the red flour beetle, Tribolium castaneum, with small levels of pairwise nuclear genetic divergence (0.033 < Nei's D < 0.125). For 15 population pairs we created a full spectrum of line crosses (two parental, two reciprocal F1's, four F2's, and eight backcrosses), reared them at multiple temperatures, and analyzed the numbers and developmental defects of offspring. We assayed a total of 219,388 offspring from 5147 families. Failed crosses occurred predominately in F2's, giving evidence of F2 breakdown within this species. In all cases where a significant model could be fit to the data on offspring number, we observed at least one type of digenic epistasis. We also found maternal and cytoplasmic effects to be common components of divergence among T. castaneum populations. In some cases, the most complex model tested (additive, dominance, epistatic, maternal, and cytoplasmic effects) did not provide a significant fit to the data, suggesting that linkage or higher order epistasis is involved in differentiation between some populations. For the limb deformity data, we observed significant genotype-by-environment interaction in most crosses and pure parent crosses tended to have fewer deformities than hybrid crosses. Complexity of genetic architecture was not correlated with either geographic distance or genetic distance. Our results support the view that genetic incompatibilities responsible for postzygotic isolation, an important component of speciation, may be a natural but serendipitous consequence of nonadditive genetic effects and structured populations.