Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

An unusual barrier to gene flow: perpetually immature larvae from inter-population crosses in the flour beetle, Tribolium castaneum

We genetically characterize an unusual hybrid incompatibility phenotype manifest in F₁ offspring of crosses between two populations of Tribolium castaneum. Hybrid larvae cease development at the third larval instar, persisting as ‘perpetually immature larvae’ thereafter. Although unable to produce viable adult hybrid offspring with one another, each population produces abundant, fertile hybrids with other populations, indicating a recent origin of the incompatibility and facilitating genetic studies. We mapped the paternal component of the hybrid phenotype to a single region, which exhibits two characteristics common to hybrid incompatibility: marker transmission ratio distortion within crosses and elevated genetic divergence between populations. The incompatible variation and an elevation in between‐population genetic divergence is associated with a region containing the T. castaneum ecdysone receptor homologue, a major regulatory switch, controlling larval moults, pupation and metamorphosis. This contributes to understanding the genetics of speciation in the Coleoptera, one of the most speciose of all arthropod taxa.     

LOCAL ADAPTATION,PHENOTYPIC DIFFERENTIATION,AND HYBRID FITNESS IN DIVERGED NATURAL POPULATIONS OF ARABIDOPSIS LYRATA

Selection for local adaptation results in genetic differentiation in ecologically important traits. In a perennial, outcrossing model plant Arabidopsis lyrata, several differentiated phenotypic traits contribute to local adaptation, as demonstrated by fitness advantage of the local population at each site in reciprocal transplant experiments. Here we compared fitness components, hierarchical total fitness and differentiation in putatively ecologically important traits of plants from two diverged parental populations from different continents in the native climate conditions of the populations in Norway and in North Carolina (NC, U.S.A.). Survival and number of fruits per inflorescence indicated local advantage at both sites and aster life‐history models provided additional evidence for local adaptation also at the level of hierarchical total fitness. Populations were also differentiated in flowering start date and floral display. We also included reciprocal experimental F₁ and F₂ hybrids to examine the genetic basis of adaptation. Surprisingly, the F₂ hybrids showed heterosis at the study site in Norway, likely because of a combination of beneficial dominance effects from different traits. At the NC site, hybrid fitness was mostly intermediate relative to the parental populations. Local cytoplasmic origin was associated with higher fitness, indicating that cytoplasmic genomes also may contribute to the evolution of local adaptation.     

Genetic architecture and phenotypic plasticity of thermally-regulated traits in an eruptive species, <Emphasis Type="Italic">Dendroctonus ponderosae</Emphasis>

Phenotypic plasticity in thermally-regulated traits enables close tracking of changing environmental conditions, and can thereby enhance the potential for rapid population increase, a hallmark of outbreak insect species. In a changing climate, exposure to conditions that exceed the capacity of existing phenotypic plasticity may occur. Combining information on genetic architecture and trait plasticity among populations that are distributed along a latitudinal cline can provide insight into how thermally-regulated traits evolve in divergent environments and the potential for adaptation. Dendroctonus ponderosae feed on Pinus species in diverse climatic regimes throughout western North America, and show eruptive population dynamics. We describe geographical patterns of plasticity in D. ponderosae development time and adult size by examining reaction norms of populations from multiple latitudes. The relative influence of additive and non-additive genetic effects on population differences in the two phenotypic traits at a single temperature is quantified using line-cross experiments and joint-scaling tests. We found significant genetic and phenotypic variation among D. ponderosae populations. Simple additive genetic variance was not the primary source of the observed variation, and dominance and epistasis contributed greatly to the genetic divergence of the two thermally-regulated traits. Hybrid breakdown was also observed in F₂ hybrid crosses between northern and southern populations, further indication of substantial genetic differences among clinal populations and potential reproductive isolation within D. ponderosae. Although it is unclear what maintains variation in the life-history traits, observed plasticity in thermally-regulated traits that are directly linked to rapid numerical change may contribute to the outbreak nature of D. ponderosae, particularly in a changing climate.     

Genetic analysis of larval survival and larval growth of two populations of Leptinotarsa decemlineata on tomato

The genetics of adaptation to tomato in Leptinotarsa decemlineata (Say) were investigated in reciprocal F₁, F₂, and backcross populations generated from crosses between beetles from a tomato adapted population and from a population that was poorly adapted to tomato. Larvae from the parent and test populations were reared on tomato for four days, after which survivorship and larval weights were recorded. Most results indicate that differences in larval growth and survival on tomato between the parent populations are largely determined by autosomal, polygenic mechanisms, the inheritance of which involves a significant dominance component. However, results from F₂ crosses are not consistent with this conclusion. A significant difference in larval weights, but not in survival, between reciprocal F₁ populations in an analysis of combined data from four separate experiments suggests that maternal cytoplasmic effects may contribute to differences in larval performance on tomato between the adapted and unadapted populations. The unusual results obtained from F₂ crosses in this study are not atypical of results from previous studies of the genetics of adaptation to host plants by the Colorado potato beetle. Host plant adaptation by Colorado potato beetles may therefore involve unusual genetic mechanisms that are not easily assessed by classical Mendelian analysis.     

Evolution of developmental homeostasis inDrosophila mojavensis

Summary Variation in life histories among populations of cactophilicDrosophila mojavensis has been hypothesized to be a by-product of a shift to one of two alternate host plants. When cultured on the ancestral and a secondary host cactus, a Baja population expressed shorter development times and smaller thorax sizes than a mainland population, but viability did not differ. Comparisons with all reciprocal F₁ and F₂ crosses between populations revealed that genetic differences in development time and thorax size were largely additive. Homeostasis in these life history traits was population specific, except for viability. Homeostasis in development time was greater in the Baja population than in the other crosses, suggesting dominance for decreased homeostasis in the mainland population. Underdominance in viability homeostasis of the F₁ hybrids suggested some incompatibility between populations. Homeostasis in thorax size was greater in females than in males and differed among parental populations. Maintenance of heritable differences and genetic variation for homeostasis in these traits suggested a role for cactus-specific differences in environmental uncertainty caused by variation in breeding site duration and abundance in nature.     

Landscape heterogeneity and local adaptation define the spatial genetic structure of Pacific salmon in a pristine environment

Identifying the spatial distribution of genetic variation across the landscape is an essential step in informing species conservation. Comparison of closely related and geographically overlapping species can be particularly useful in cases where landscape may similarly influence genetic structure. Congruent patterns among species highlight the importance that landscape heterogeneity plays in determining genetic structure whereas contrasting patterns emphasize differences in species-specific ecology and life-history or the importance of species-specific adaptation to local environments. We examined the interacting roles of demography and adaptation in determining spatial genetic structure in two closely related and geographically overlapping species in a pristine environment. Using single nucleotide polymorphism (SNP) loci exhibiting both neutral and putative adaptive variation, we evaluated the genetic structure of sockeye salmon in the Copper River, Alaska; these data were compared to existing data for Chinook salmon from the same region. Overall, both species exhibited patterns of isolation by distance; the spatial distribution of populations largely determined the distribution of genetic variation across the landscape. Further, both species exhibited largely congruent patterns of within- and among-population genetic diversity, highlighting the role that landscape heterogeneity and historical processes play in determining spatial genetic structure. Potential adaptive differences among geographically proximate sockeye salmon populations were observed when high F_ST outlier SNPs were evaluated in a landscape genetics context. Results were evaluated in the context of conservation efforts with an emphasis on reproductive isolation, historical processes, and local adaptation.     

Effects on Embryo Development Time and Survival of Intercrossing Three Geographically Separate Populations of Southeast Alaska coho salmon, Oncorhynchus kisutch

We investigated adaptive differences among three geographically separate populations of coho salmon (Oncorhynchus kisutch) by forming first generation intercrosses (hybrid lines) and comparing them to parental types (control lines). Broodstock for the experiment came from the Gastineau Hatchery in Juneau, Hidden Falls Hatchery on Baranof Island, and Neets Bay Hatchery near Ketchikan, Alaska. All were isolated hatchery populations separated by 220–400 km and derived 15–20 years previously from single local wild populations. For each population, gametes were taken from 50 mature salmon of each sex and combined to form nine lines (three control and six hybrid); each line had 50 full-sibling families which were assigned to separate cells of an incubator at Gastineau Hatchery. Embryo survival and development times were measured as indicators of locally adapted fitness traits. Two of the control lines had higher survival rates than hybrid lines formed between either of their parental populations and other populations. Differences (p < 0.05) were found between development times for control and hybrid groups, which varied by as many as 20 days between families and as many as 30 days between control lines. The intermediate expression of development time of intercrossed lines is consistent with additive genetic variation of development time between the ancestral populations of coho salmon and indicates that important genetic divergence exists between the populations. A loss of local adaptation through a change in seasonal timing of completion of embryonic development would occur in intercrosses between the populations.     

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene‐targeted PCR‐based assay for next‐generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F _ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single‐gene Mantel tests comparing methylation and pairwise F _ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.     

Recent local adaptation of sockeye salmon to glacial spawning habitats

Salmonids spawn in highly diverse habitats, exhibit strong genetic population structuring, and can quickly colonize newly created habitats with few founders. Spawning traits often differ among populations, but it is largely unknown if these differences are adaptive or due to genetic drift. To test if sockeye salmon (Oncorhynchus nerka) populations are adapted to glacial, beach, and tributary spawning habitats, we examined variation in heritable phenotypic traits associated with spawning in 13 populations of wild sockeye salmon in Lake Clark, Alaska. These populations were commonly founded between 100 and 400 hundred sockeye salmon generations ago and exhibit low genetic divergence at 11 microsatellite loci (F _ST < 0.024) that is uncorrelated with spawning habitat type. We found that mean P _ST (phenotypic divergence among populations) exceeded neutral F _ST for most phenotypic traits measured, indicating that phenotypic differences among populations could not be explained by genetic drift alone. Phenotypic divergence among populations was associated with spawning habitat differences, but not with neutral genetic divergence. For example, female body color was lighter and egg color was darker in glacial than non-glacial habitats. This may be due to reduced sexual selection for red spawning color in glacial habitats and an apparent trade-off in carotenoid allocation to body and egg color in females. Phenotypic plasticity is an unlikely source of phenotypic differences because Lake Clark sockeye salmon spend nearly all their lives in a common environment. Our data suggest that Lake Clark sockeye salmon populations are adapted to spawning in glacial, beach and tributary habitats and provide the first evidence of a glacial spawning ecotype in salmonids. Glacial spawning habitats are often young (i.e., <200 years old) and ephemeral. Thus, local adaptation of sockeye salmon to glacial habitats appears to have occurred recently.     

Fine‐scale temporal adaptation within a salmonid population: mechanism and consequences

We demonstrate a clear example of local adaptation of seasonal timing of spawning and embryo development. The consequence is a population of pink salmon that is segmented into spawning groups that use the same limited habitat. We synthesize published observations with results of new analyses to demonstrate that genetic variation of these traits results in survival differentials related to that variation, and that density‐dependent embryo mortality and seasonally variable juvenile mortality are a mechanism of selection. Most examples of local adaptation in natural systems depend on observed correlations between environments and fitness traits, but do not fully demonstrate local adaptation: that the trait is genetically determined, exhibits different fitness in common environments or across different environments, and its variation is mechanistically connected to fitness differences. The geographic or temporal scales of local adaptation often remain obscure. Here, we show that heritable, fine‐scale differences of timing of reproductive migration in a pink salmon (Oncorhynchus gorbuscha) resulted in temporal structure that persisted several generations; the differences enable a density‐dependent population to pack more spawners into limited spawning habitat, that is, enhance its fitness. A balanced trade‐off of survivals results because embryos from early‐migrating fish have a lower freshwater survival (harsh early physical conditions and disturbance by late spawners), but emigrant fry from late‐migrating fish have lower marine survivals (timing of their vernal emergence into the estuarine environment). Such fine‐scale local adaptations increase the genetic portfolio of the populations and may provide a buffer against the impacts of climate change.     

LOCAL VARIATION IN THE GENETIC BASIS OF PAEDOMORPHOSIS IN THE SALAMANDER AMBYSTOMA TALPOIDEUM

The hypothesis that local isolated populations differed in the genetic basis for life-history traits was tested in the salamander Ambystoma talpoideum. Genetic basis was defined as the specific genetic architecture (additive and nonadditive) that contributes, along with maternal and environmental factors, to the phenotype. All crosses within and between three populations were made to produce nine F₁ populations. Nine within-population crosses produced the F₂ generation. This design does not permit an estimation of the exact nature of the genetic basis (e.g., additive, nonadditive) for any trait within populations. However, hybrid dissimilarity in the F₂ generation was taken as evidence of a different genetic basis for a trait in each population. The genetic basis of life-history pathway (metamorphosis vs. paedomorphosis) and per capita fecundity differed between two populations. The genetic basis of life-history pathway, per capita fecundity, survival, and growth rate was similar between the remaining sets of populations. This study and related ones (Semlitsch and Wilbur, 1989; Semlitsch et al., 1990) suggest that a heterochronic shift that causes rapid morphological evolution between metamorphosis and paedomorphosis (a macroevolutionary pattern) can evolve independently and does not require a macromutation or other nonmicroevolutionary mechanisms.     

Chromosomal fusion and life history‐associated genomic variation contribute to within‐river local adaptation of Atlantic salmon

Chromosomal inversions have been implicated in facilitating adaptation in the face of high levels of gene flow, but whether chromosomal fusions also have similar potential remains poorly understood. Atlantic salmon are usually characterized by population structure at multiple spatial scales; however, this is not the case for tributaries of the Miramichi River in North America. To resolve genetic relationships between populations in this system and the potential for known chromosomal fusions to contribute to adaptation, we genotyped 728 juvenile salmon using a 50 K SNP array. Consistent with previous work, we report extremely weak overall population structuring (Global F_ST = 0.004) and failed to support hierarchical structure between the river's two main branches. We provide the first genomic characterization of a previously described polymorphic fusion between chromosomes 8 and 29. Fusion genomic characteristics included high LD, reduced heterozygosity in the fused homokaryotes, and strong divergence between the fused and the unfused rearrangement. Population structure based on fusion karyotype was five times stronger than neutral variation (F_ST = 0.019), and the frequency of the fusion was associated with summer precipitation supporting a hypothesis that this rearrangement may contribute local adaptation despite weak neutral differentiation. Additionally, both outlier variation among populations and a polygenic framework for characterizing adaptive variation in relation to climate identified a 250‐Kb region of chromosome 9, including the gene six6 that has previously been linked to age‐at‐maturity and run‐timing for this species. Overall, our results indicate that adaptive processes, independent of major river branching, are more important than neutral processes for structuring these populations.     

Effects of genetic distance on heterosis in a <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> model system

Habitat fragmentation and small population sizes can lead to inbreeding and loss of genetic variation, which can potentially cause inbreeding depression and decrease the ability of populations to adapt to altered environmental conditions. One solution to these genetic problems is the implementation of genetic rescue, which re-establishes gene flow between separated populations. Similar techniques are being used in animal and plant breeding to produce superior production animals and plants. To optimize fitness benefits in genetic rescue programs and to secure high yielding domestic varieties in animal and plant breeding, knowledge on the genetic relatedness of populations being crossed is imperative. In this study, we conducted replicated crosses between isogenic Drosophila melanogaster lines from the Drosophila Genetic Reference Panel. We grouped lines in two genetic distance groups to study the effect of genetic divergence between populations on the expression of heterosis in two fitness components; starvation resistance and reproductive output. We further investigated the transgenerational effects of outcrossing by investigating the fitness consequences in both the F₁- and the F₃-generations. High fitness enhancements were observed in hybrid offspring compared to parental lines, especially for reproductive output. However, the level of heterosis declined from the F₁- to the F₃-generation. Generally, genetic distance did not have strong impact on the level of heterosis detected, although there were exceptions to this pattern. The best predictor of heterosis was performance of parental lines with poorly performing parental lines showing higher hybrid vigour when crossed, i.e. the potential for heterosis was proportional to the level of inbreeding depression. Overall, our results show that outcrossing can have very strong positive fitness consequences for genetically depauperate populations.     

Low but significant genetic differentiation underlies biologically meaningful phenotypic divergence in a large Atlantic salmon population

Despite decades of research assessing the genetic structure of natural populations, the biological meaning of low yet significant genetic divergence often remains unclear due to a lack of associated phenotypic and ecological information. At the same time, structured populations with low genetic divergence and overlapping boundaries can potentially provide excellent models to study adaptation and reproductive isolation in cases where high‐resolution genetic markers and relevant phenotypic and life history information are available. Here, we combined single nucleotide polymorphism (SNP)‐based population inference with extensive phenotypic and life history data to identify potential biological mechanisms driving fine‐scale subpopulation differentiation in Atlantic salmon (Salmo salar) from the Teno River, a major salmon river in Europe. Two sympatrically occurring subpopulations had low but significant genetic differentiation (F_ST = 0.018) and displayed marked differences in the distribution of life history strategies, including variation in juvenile growth rate, age at maturity and size within age classes. Large, late‐maturing individuals were virtually absent from one of the two subpopulations, and there were significant differences in juvenile growth rates and size at age after oceanic migration between individuals in the respective subpopulations. Our findings suggest that different evolutionary processes affect each subpopulation and that hybridization and subsequent selection may maintain low genetic differentiation without hindering adaptive divergence.     

Multigenerational outbreeding effects in Chinook salmon (Oncorhynchus tshawytscha)

Outbreeding, mating between genetically divergent individuals, may result in negative fitness consequences for offspring via outbreeding depression. Outbreeding effects are of notable concern in salmonid research as outbreeding can have major implications for salmon aquaculture and conservation management. We therefore quantified outbreeding effects in two generations (F₁ hybrids and F₂ backcrossed hybrids) of Chinook salmon (Oncorhynchus tshawytscha) derived from captively-reared purebred lines that had been selectively bred for differential performance based on disease resistance and growth rate. Parental lines were crossed in 2009 to create purebred and reciprocal hybrid crosses (n = 53 families), and in 2010 parental and hybrid crosses were crossed to create purebred and backcrossed hybrid crosses (n = 66 families). Although we found significant genetic divergence between the parental lines (F_ST = 0.130), reciprocal F₁ hybrids showed no evidence of outbreeding depression (hybrid breakdown) or favorable heterosis for weight, length, condition or survival. The F₂ backcrossed hybrids showed no outbreeding depression for a suite of fitness related traits measured from egg to sexually mature adult life stages. Our study contributes to the current knowledge of outbreeding effects in salmonids and supports the need for more research to better comprehend the mechanisms driving outbreeding depression.     

Spatial and temporal genetic structure in a hybrid cordgrass invasion

Invasive hybrids and their spread dynamics pose unique opportunities to study evolutionary processes. Invasive hybrids of native Spartina foliosa and introduced S. alterniflora have expanded throughout San Francisco Bay intertidal habitats within the past 35 years by deliberate plantation and seeds floating on the tide. Our goals were to assess spatial and temporal scales of genetic structure in Spartina hybrid populations within the context of colonization history. We genotyped adult and seedling Spartina using 17 microsatellite loci and mapped their locations in three populations. All sampled seedlings were hybrids. Bayesian ordination analysis distinguished hybrid populations from parent species, clearly separated the population that originated by plantation from populations that originated naturally by seed and aligned most seedlings within each population. Population genetic structure estimated by analysis of molecular variance was substantial (F_ST=0.21). Temporal genetic structure among age classes varied highly between populations. At one population, the divergence between adults and 2004 seedlings was low (F_ST=0.02) whereas at another population this divergence was high (F_ST=0.26). This latter result was consistent with local recruitment of self-fertilized seed produced by only a few parental plants. We found fine-scale spatial genetic structure at distances less than ∼200 m, further supporting local seed and/or pollen dispersal. We posit a few self-fertile plants dominating local recruitment created substantial spatial genetic structure despite initial long-distance, human dispersal of hybrid Spartina through San Francisco Bay. Fine-scale genetic structure may more strongly develop when local recruits are dominated by the offspring of a few self-fertile plants.     

Signatures of diversifying selection at EST‐SSR loci and association with climate in natural Eucalyptus populations

Understanding the environmental parameters that drive adaptation among populations is important in predicting how species may respond to global climatic changes and how gene pools might be managed to conserve adaptive genetic diversity. Here, we used Bayesian F_ST outlier tests and allele–climate association analyses to reveal two Eucalyptus EST‐SSR loci as strong candidates for diversifying selection in natural populations of a southwestern Australian forest tree, Eucalyptus gomphocephala (Myrtaceae). The Eucalyptus homolog of a CONSTANS‐like gene was an F_ST outlier, and allelic variation showed significant latitudinal clinal associations with annual and winter solar radiation, potential evaporation, summer precipitation and aridity. A second F_ST outlier locus, homologous to quinone oxidoreductase, was significantly associated with measures of temperature range, high summer temperature and summer solar radiation, with important implications for predicting the effect of temperature on natural populations in the context of climate change. We complemented these data with investigations into neutral population genetic structure and diversity throughout the species range. This study provides an investigation into selection signatures at gene‐homologous EST‐SSRs in natural Eucalyptus populations, and contributes to our understanding of the relationship between climate and adaptive genetic variation, informing the conservation of both putatively neutral and adaptive components of genetic diversity.     

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.