Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Conditions for the Maintenance of Prezygotic and Zygotic Isolation in a Continent-Island Model

A continent-island model is studied in which differences are initially present at a gene locus with assortative mating (which causes prezygotic reproductive isolation) and at one or two unlinked selected loci, (which cause zygotic reproductive isolation). This model simulates secondary contact between two populations which have diverged allopatrically at loci which control isolating mechanisms. Two different models of prezygotic isolation and three different models of zygotic isolation are studied. The analysis is focused on those conditions under which differentiation is not lost through the swamping effect of gene now. In each case the critical migration rate (m_c) below which differences can be maintained is found. Results depend on the degree of assortative mating, on the level of selection on the genes involved, and on the particular models of reproductive isolation considered. However, in all the models, for high values of such parameters the two isolating mechanisms interact strongly, and consequently the stability of the difference at each locus (as measured by m_c) is noticeably influenced by the presence of the difference at the other locus (or loci). The nature and strength of the effects of this interaction differ according to whether or not there is selection on the assortative mating locus and whether selection against hybrids or directional selection is hypothesized.     

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation () and introgression (genomic cline center [α] and rate [β]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of , α, and β were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated , α, and β were not strongly co-localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.     

How chromosomal rearrangements shape adaptation and speciation: Case studies in Drosophila pseudoobscura and its sibling species Drosophila persimilis

The gene arrangements of Drosophila have played a prominent role in the history of evolutionary biology from the original quantification of genetic diversity to current studies of the mechanisms for the origin and establishment of new inversion mutations within populations and their subsequent fixation between species supporting reproductive barriers. This review examines the genetic causes and consequences of inversions as recombination suppressors and the role that recombination suppression plays in establishing inversions in populations as they are involved in adaptation within heterogeneous environments. This often results in the formation of clines of gene arrangement frequencies among populations. Recombination suppression leads to the differentiation of the gene arrangements which may accelerate the accumulation of fixed genetic differences among populations. If these fixed mutations cause incompatibilities, then inversions pose important reproductive barriers between species. This review uses the evolution of inversions in Drosophila pseudoobscura and D. persimilis as a case study for how inversions originate, establish and contribute to the evolution of reproductive isolation.     

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural ‘replicate’ hybrid zones between two divergent Populus species via locus‐specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between‐species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited ‘hybrid habitats’, our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.     

Morphological stasis in the Eurytemora affinis species complex (Copepoda: Temoridae)

Morphological stasis has long been regarded as one of the most challenging problems in evolutionary biology. This study focused on the copepod species complex, Eurytemora affinis, as a model system to determine pattern and degree of morphological stasis. This study revealed discordant rates of morphological differentiation, molecular evolution, and reproductive isolation, where speciation was accompanied by lack of morphological differentiation in secondary sex characters. Comparisons were made among phylogenies based on morphometrics, nuclear (allozyme) loci, and mitochondrial DNA (mtDNA) sequences from cytochrome oxidase I, for a total of 43 populations within the complex. These systematic relationships were also compared to patterns of reproductive isolation. In addition, genetic subdivision of nuclear molecular (allozyme) markers (G _ST) and quantitative (morphological) characters (Q _ST) were determined to infer evolutionary forces driving morphological differentiation. The morphometric phylogeny revealed that all clades, excluding the European clade, were morphologically undifferentiated and formed a polytomy (multifurcation). Morphometric distances were not correlated with mtDNA distances, or with patterns of reproductive isolation. In contrast, nuclear and mtDNA phylogenies were mostly congruent. Reproductive isolation proved to be the most sensitive indicator of speciation, given that two genetically and morphologically proximate populations showed evidence of hybrid breakdown. Quantitative genetic (morphological) subdivision (Q _ST = 0.162) was lower than nuclear genetic subdivision (G _ST = 0.617) for four laboratory-reared North American populations, indicating retarded evolution of morphological characters. This result contrasts with most other species, where Q _ST typically exceeds G _ST as a result of directional selection. Thus, in all but the European populations, evolution of the secondary sex characters was marked by morphological stasis, even between reproductively-isolated populations.     

Interspecific Tests of Allelism Reveal the Evolutionary Timing and Pattern of Accumulation of Reproductive Isolation Mutations

Despite extensive theory, little is known about the empirical accumulation and evolutionary timing of mutations that contribute to speciation. Here we combined QTL (Quantitative Trait Loci) analyses of reproductive isolation, with information on species evolutionary relationships, to reconstruct the order and timing of mutations contributing to reproductive isolation between three plant (Solanum) species. To evaluate whether reproductive isolation QTL that appear to coincide in more than one species pair are homologous, we used cross-specific tests of allelism and found evidence for both homologous and lineage-specific (non-homologous) alleles at these co-localized loci. These data, along with isolation QTL unique to single species pairs, indicate that >85% of isolation-causing mutations arose later in the history of divergence between species. Phylogenetically explicit analyses of these data support non-linear models of accumulation of hybrid incompatibility, although the specific best-fit model differs between seed (pairwise interactions) and pollen (multi-locus interactions) sterility traits. Our findings corroborate theory that predicts an acceleration (‘snowballing’) in the accumulation of isolation loci as lineages progressively diverge, and suggest different underlying genetic bases for pollen versus seed sterility. Pollen sterility in particular appears to be due to complex genetic interactions, and we show this is consistent with a snowball model where later arising mutations are more likely to be involved in pairwise or multi-locus interactions that specifically involve ancestral alleles, compared to earlier arising mutations.     

RAPD analysis of genetic diversity and population genetic structure of Stipa krylovii reshov. in Inner Mongolia steppe

Random amplified polymorphic DNA (RAPD) analysis was used to characterize the genetic diversity and population genetic structure of Stipa krylovii populations in Inner Mongolia steppe of North China. Thirteen 10-bp oligonucleotide primers, which generated 237 RAPD bands, were used to analyze 90 plants of five populations from three regions, meadow steppe, typical steppe and desert steppe, from the east to the west. The genetic diversity of Stipa krylovii that was revealed by observed number of alleles (na), expected number of alleles (ne), Nei’s diversity index (h), Shannon’s diversity index (H), amplificated loci, polymorphic loci and the percentage of polymorphic loci (PPB) increased from the east to the west. The Pearson’s correlation analysis between genetic diversity parameters and ecological parameters indicated that the genetic diversity of Stipa krylovii was associated with precipitation and cumulative temperature variations along the longitude (humidity were calculated by precipitation and cumulative temperature). Dendrogram based on Jaccard’s genetic distance showed that the individuals from the same population formed a single subgroup. Although most variation (56.85%) was within populations, there was high genetic differentiation among populations of Stipa krylovii, high differentiation within and between regions by AMOVA analysis. Either Nei’s unbiased genetic distance (G _ST) or gene flow (Nm) among pairwise populations was not correlated with geographical distance by Mantel’s test (P > 0.05), suggesting that there was no consistency with the isolation by distance model in these populations. Natural selection may have played a role in affecting the genetic diversity and population structure, but habitat destruction and degradation in northem grassland in China may be the main factor responsible for high genetic differentiation among populations, within and among regions. The text was submitted by the authors in English.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

GENETIC HITCHHIKING AND THE DYNAMIC BUILDUP OF GENOMIC DIVERGENCE DURING SPECIATION WITH GENE FLOW

A major issue in evolutionary biology is explaining patterns of differentiation observed in population genomic data, as divergence can be due to both direct selection on a locus and genetic hitchhiking. “Divergence hitchhiking” (DH) theory postulates that divergent selection on a locus reduces gene flow at physically linked sites, facilitating the formation of localized clusters of tightly linked, diverged loci. “Genome hitchhiking” (GH) theory emphasizes genome‐wide effects of divergent selection. Past theoretical investigations of DH and GH focused on static snapshots of divergence. Here, we used simulations assessing a variety of strengths of selection, migration rates, population sizes, and mutation rates to investigate the relative importance of direct selection, GH, and DH in facilitating the dynamic buildup of genomic divergence as speciation proceeds through time. When divergently selected mutations were limiting, GH promoted divergence, but DH had little measurable effect. When populations were small and divergently selected mutations were common, DH enhanced the accumulation of weakly selected mutations, but this contributed little to reproductive isolation. In general, GH promoted reproductive isolation by reducing effective migration rates below that due to direct selection alone, and was important for genome‐wide “congealing” or “coupling” of differentiation (F_ST) across loci as speciation progressed.     

Chromosomal rearrangements directly cause underdominant F1 pollen sterility in Mimulus lewisii–Mimulus cardinalis hybrids

Chromosomal rearrangements can contribute to the evolution of postzygotic reproductive isolation directly, by disrupting meiosis in F₁ hybrids, or indirectly, by suppressing recombination among genic incompatibilities. Because direct effects of rearrangements on fertility imply fitness costs during their spread, understanding the mechanism of F₁ hybrid sterility is integral to reconstructing the role(s) of rearrangements in speciation. In hybrids between monkeyflowers Mimulus cardinalis and Mimulus lewisii, rearrangements contain all quantitative trait loci (QTLs) for both premating barriers and pollen sterility, suggesting that they may have facilitated speciation in this model system. We used artificial chromosome doubling and comparative mapping to test whether heterozygous rearrangements directly cause underdominant male sterility in M. lewisii–M. cardinalis hybrids. Consistent with a direct chromosomal basis for hybrid sterility, synthetic tetraploid F₁s showed highly restored fertility (83.4% pollen fertility) relative to diploids F₁s (36.0%). Additional mapping with Mimulus parishii–M. cardinalis and M. parishii–M. lewisii hybrids demonstrated that underdominant male sterility is caused by one M. lewisii specific and one M. cardinalis specific reciprocal translocation, but that inversions had no direct effects on fertility. We discuss the importance of translocations as causes of reproductive isolation, and consider models for how underdominant rearrangements spread and fix despite intrinsic fitness costs.     

Isolation and characterization of microsatellite loci from the apple maggot fly Rhagoletis pomonella (Diptera: Tephritidae)

We report the isolation and development of 81 novel primers for amplifying microsatellite loci in the Rhagoletis pomonella sibling species complex, and the sequencing, characterization and analysis of basic population genetic parameters for nine of these genes. We also report the successful cross‐species amplification of several of these loci. The R. pomonella sibling species complex is a textbook example of genetic differentiation in sympatry via host‐plant shifting. Microsatellite markers can be useful for mapping host‐plant‐associated adaptations in Rhagoletis that generate reproductive isolation and facilitate speciation, as well as for resolving the genetic structure and evolutionary history of fly populations.     

On the evolution of heavy-metal tolerant populations in Armeria maritima: evidence from allozyme variation and reproductive barriers

The process of ecological differentiation leading to the evolution of heavy-metal tolerant populations in Armeria maritima was studied by comparing population genetic structure and pattern of gene flow between populations growing on heavy-metal contaminated against non-contaminated sites using allozyme markers. In addition the evolution of reproductive isolation among populations was studied by measuring pollen fertility in interpopulational hybrids. The allozyme data suggested that in A. maritima multiple independent evolutionary origins of heavy-metal tolerant populations have occurred in the absence of strong genetic bottlenecks. The pattern of gene flow among populations was consistent with the model of isolation by distance with considerable gene flow between neighbor populations, and no reduction of gene flow between tolerant and non-tolerant populations. Hence it appears that substantial gene flow has not hampered genetic differentiation, probably because of the high selection pressure for heavy-metal tolerance. The pattern of reproductive isolation among populations suggests that evolution of heavy-metal tolerant populations has not triggered the development of reproductive barriers against non-tolerant populations. However, partial reproductive isolation has occurred under geographic separation.     

Reproductive barriers between two sympatric beetle species specialized on different host plants

Knowledge on interspecific pre‐ and post‐zygotic isolation mechanisms provides insights into speciation patterns. Using crosses (F₁ and backcrosses) of two closely related flea beetles species, Altica fragariae and A. viridicyanea, specialized on different hosts in sympatry, we measured: (a) the type of reproductive isolation and (b) the inheritance mode of preference and host‐specific performance, using a joint‐scaling test. Each species preferred almost exclusively its host plant, creating strong prezygotic isolation between them, and suggesting that speciation may occur at least partly in sympatry. Reproductive isolation was intrinsic between females of A. fragariae and either A. viridicyanea or F₁ males, whereas the other crosses showed ecologically dependent reproductive isolation, suggesting ecological speciation. The genetic basis of preference and performance was at least partially independent, and several loci coded for preference, which limits the possibility of sympatric speciation. Hence, both ecological and intrinsic factors may contribute to speciation between these species.     

The evolution of intrinsic reproductive isolation in the genus Cakile (Brassicaceae)

In theory, adaptive divergence can increase intrinsic post‐zygotic reproductive isolation (RI), either directly via selection on loci associated with RI, or indirectly via linkage of incompatibility loci with loci under selection. To test this hypothesis, we measured RI at five intrinsic post‐zygotic reproductive barriers between 18 taxa from the genera Cakile and Erucaria (Brassicaceae). Using a comparative framework, we tested whether the magnitude of RI was associated with genetic distance, geographic distance, ecological divergence and parental mating system. Early stages of post‐zygotic RI related to F₁ viability (i.e. initial seed set) tended to be stronger than later stages related to F₁ fecundity (i.e. flower number, fruit number). Mating system significantly influenced early stages of RI, such that RI was lowest when the mother was selfing and father was outcrossing, consistent with an imbalance between sink strength and resistance to provisioning. We found little evidence that adaptive divergence accelerates the evolution of intrinsic post‐zygotic RI, consistent with a nonecological model of evolution that predicts the nonlinear accumulation of RI and RI asymmetry with time (i.e. genetic distance), irrespective of adaptive divergence. Thus, although certain aspects of ecological divergence do not appear to have contributed strongly to the evolution of RI in this system, divergence in mating system actually reduced RI, suggesting that mating system evolution may play a significant role in speciation dynamics.     

Is the between-population variance negligible in the total variance of heterozygosity? Case of a finite number of loci subject to the infinite-allele model in finite monoecious populations

The decomposition of the variance of the average heterozygosity into variances between and within populations is studied in the general case of a finite number of loci. These loci are assumed randomly distributed over chromosome pairs having a non-interference recombination scheme, and independently subject to mutation according to the infinite-allele model. The equilibrium behavior of that decomposition is discussed in the monoecious mating case with regard to each parameter of the model: mutation rate per gene per generation (u), population size (N), number of loci (n), map length of chomosome pairs (L). It is shown that the proportion Q of the between-population variability in the total variance of the average heterozygosity is decreasing as either the mean heterozygosity (θ = 4Nu/(1 + 4Nu)) or the mean number of mutations per gamete per generation (v = nu) is increasing. Moreover, even if Q is always smaller than for this model, it is not negligible unless θ is close to one or v is much larger than one for L long enough.     

Weak premating isolation between Clitarchus stick insect species despite divergent male and female genital morphology

Documenting natural hybrid systems builds our understanding of mate choice, reproductive isolation and speciation. The stick insect species Clitarchus hookeri and C. tepaki differ in their genital morphology and hybridize along a narrow peninsula in northern New Zealand. We utilize three lines of evidence to understand the role of premating isolation and species boundaries: (a) genetic differentiation using microsatellites and mitochondrial DNA; (b) variation in 3D surface topology of male claspers and 2D morphometrics of female opercular organs; and (c) behavioural reproductive isolation among parental and hybrid populations through mating crosses. The genetic data show introgression between the parental species and formation of a genetically variable hybrid swarm. Similarly, the male and female morphometric data show genital divergence between the parental species as well as increased variation within the hybrid populations. This genital divergence has not resulted in reproductive isolation between species, instead weak perimating isolation has enabled the formation of a hybrid swarm. Behavioural analysis demonstrates that the entire mating process influences the degree of reproductive isolation between species undergoing secondary contact. Mechanical isolation may appear strong, whereas perimating isolation is weak.     

Genetic Divergence of Chars of the Genus <Emphasis Type="Italic">Salvelinus</Emphasis> from Kronotsky Lake (Kamchatka Peninsula)

Chars of the genus Salvelinus, inhabiting lakes and lake-river systems, comprise morphological and ecological forms whose taxonomic status is under dispute. In the present work, we have examined genetic variation and divergence in various chars from the Kronotsky lake basin: the lacustrine chars (white, nose, and longhead) and riverine Dolly Varden Salvelinus malma. The study was conducted using analysis of allozyme and microsatellite loci, myogens, RAPD, and restriction analysis of two mtDNA segments. The estimates of heterozygosity at allozyme and microsatellite loci were similar to the corresponding parameters in populations of northern Dolly Varden and Arctic char. Heterozygote deficit was recorded in both samples of separate forms, and in the total sample of all chars from Kronotsky Lake. For allozyme and microsatellite loci, appreciable genetic differentiation among the samples of different char forms was found, which was comparable to that among the spatially isolated populations of northern Dolly Varden. This result indicates reproductive isolation among the char forms examined. However, this isolation is not complete, because no fixed differences between the forms by any of the genetic systems analyzed were found. The genetic differentiation among different forms of lacustrine chars, which corresponds to the interpopulation rather than interspecific level, is thought to be explained by their comparatively recent divergence.__________Translated from Genetika, Vol. 41, No. 8, 2005, pp. 1096–1107.Original Russian Text Copyright © 2005 by Salmenkova, Omel’chenko, Radchenko, Gordeeva, Rubtsova, Romanov.     

Genomic evidence for convergent evolution of a key trait underlying divergence in island birds

Reproductive isolation can be initiated by changes in one or a few key traits that prevent random mating among individuals in a population. During the early stages of speciation, when isolation is often incomplete, there will be a heterogeneous pattern of differentiation across regions of the genome between diverging populations, with loci controlling these key traits appearing the most distinct as a result of strong diversifying selection. In this study, we used Illumina‐sequenced ddRAD tags to identify genomewide patterns of differentiation in three recently diverged island populations of the Monarcha castaneiventris flycatcher of the Solomon Islands. Populations of this species have diverged in plumage colour, and these differences in plumage colour, in turn, are used in conspecific recognition and likely important in reproductive isolation. Previous candidate gene sequencing identified point mutations in MC1R and ASIP, both known pigmentation genes, to be associated with the difference in plumage colour between islands. Here, we show that background levels of genomic differentiation based on over 70,000 SNPs are extremely low between populations of distinct plumage colour, with no loci reaching the level of differentiation found in either candidate gene. Further, we found that a phylogenetic analysis based on these SNPs produced a taxonomy wherein the two melanic populations appear to have evolved convergently, rather than from a single common ancestor, in contrast to their original classification as a single subspecies. Finally, we found evidence that the pattern of low genomic differentiation is the result of both incomplete lineage sorting and gene flow between populations.     

The population biology of Gehyra (Gekkonidae). III. Patterns of microgeographic variation

The distribution of allozymic and chromosomal polymorphisms was examined among central Australian populations of the chromosomally variable genus Gehyra to assess whether they typically have the small deme size and low gene flow levels required by some models of chromosomal speciation. Particular attention was given to comparisons between rock-specialists (Gehyra nana) and habitat-generalists (Gehyra variegata) to investigate whether the former have more restricted gene flow. Both allozyme and chromosome data sets showed greater among population differentiation in the rock-specialists than the habitat-generalists, consistent with predictions from a previous ecological study (Moritz, 1987), although this pattern could also be due to historical effects. This was evident from the relationships between genetic and geographic distance, the conditional frequency of alleles, and F-statistics. However, both taxa appear to have substantial levels of gene flow. This indicates that Gehyra populations typically do not meet the stringent conditions for the fixation of strongly underdominant chromosome rearrangements through strong genetic drift. A consistent deficiency of heterozygotes does, however, suggest the possibility of inbreeding which would increase the likelihood of the establishment of underdominant rearrangements.