Hierarchical genetic structure and gene flow in three sympatric species of Amazonian rodents

The population genetic structure of three species of Amazonian rodents ( Oligoryzomys microtis, Oryzomys capito , and Mesomys hispidus ) is examined for mtDNA sequence haplotypes of the cytochrome b gene by hierarchical analysis of variance and gene flow estimates based on fixation indices ( N _ST) and coalescence methods. Species samples are from the same localities along 1000 km of the Rio Juruá in western Amazonian Brazil, but each species differs in important life history traits such as population size and reproductive rate. Average haplotype differentiation, hierarchical haplotype apportionment, and gene flow estimates are contrasted in discussing the current and past population structure. Two species exhibit isolation by distance patterns wherein gene flow is largely limited to geographically adjacent localities. Mesomys exhibits this pattern throughout its range along the river. More than 75% of haplotype variation is apportioned among localities and regions, and estimates of Nm for pair-wise comparisons are nearly always less than 1. Oligoryzomys shows weak isolation by distance, but only over the largest geographical distances. Nm values for this species are nearly always above 1 and most (about 80%) of haplotype variation is contained within local populations. In contrast, Oryzomys exhibits no genetic structure throughout its entire distribution; Nm values average 17 and nearly 90% of the total haplotype variance is contained within local populations. Although gene flow estimates are high, the pattern of Nm as a function of geographical distance suggests that this species experienced a more recent invasion of the region and is still in genetic disequilibrium under its current demographic conditions.     

Estimates of Gene Flow among Populations, Geographic Races, and Species in the Ipomopsis Aggregata Complex

Interpopulational gene flow within a species can reduce population differentiation due to genetic drift, whereas genetic exchange among taxa can impede speciation. We used allozyme data to estimate gene flow within and among geographic races and species of perennial herbs in the Ipomopsis aggregata complex (Polemoniaceae). Estimates of interpopulational gene flow within taxa from two methods (F statistics and private alleles) were correlated with one another. Gene flow among populations within each geographic race (subspecies) of I. aggregata was relatively high (Nm greater than approximately 1.0). Gene flow was also high among populations of I. arizonica and among four northern populations of I. tenuituba. However, gene flow was low (Nm less than 1.0) for I. tenuituba when a population representing subsp. macrosiphon was included. This is consistent with previous findings that subsp. macrosiphon has had an independent origin and is reproductively, as well as geographically, isolated. A recently developed model, based on hierarchical F statistics, was employed to estimate genetic exchange among taxa. Gene flow estimates were generally high among races of I. aggregata (dNmrace greater than 1.0) but were low among subspecies of I. tenuituba (dNmrace less than 1.0). Consistent with morphological evidence, estimates of interspecific gene flow were moderate between I. aggregata and I. tenuituba, which hybridize in several areas. However, contrary to morphological evidence, we estimated relatively high levels of interspecific gene flow involving I. arizonica. Our results suggest that I. arizonica has hybridized with other species without the transfer of morphological traits.(ABSTRACT TRUNCATED AT 250 WORDS)     

Microsatellite variation within and among populations of Oryza officinalis (Poaceae), an endangered wild rice from China

Oryza officinalis Wall. ex Watt. is an agriculturally important but seriously endangered species of wild rice. To obtain more accurate estimates of population structure for improved conservation planning of the species, genetic variability at 14 microsatellite DNA loci was examined in population samples covering most of the species' range in China. Considerable genetic variability (overall Na = 1.886, P = 62%, HO = 0.056, HE = 0.216, and HS = 0.277) was detected at the 14 loci in 442 individuals of the 18 natural populations. The evaluation of partitioning of genetic variability (FST = 0.442) suggested high genetic differentiation among the Chinese O. officinalis populations. An overall value of Nm = 0.316 suggested limited gene flow occurred among the sampled populations. Most of the populations showed heterozygote deficits in tests of Hardy-Weinberg equilibrium and significantly positive FIS values. This could be due to some inbreeding occurring in this predominantly outcrossing species. For effective in situ conservation and restoration genetics, maintenance of significant historical processes is particularly important, including high outbreeding, considerable gene flow, and large population effective sizes. The high FST values detected among populations in this study are instructive for adopting a conservation plan that includes representative populations with the greatest genetic variation for either in situ conservation management or germplasm collection expeditions.     

Genetic diversity within and among sinai populations of three Ballota species (Lamiaceae)

Ballota undulata, Ballota kaiseri, and Ballota saxatilis are very rare (and endemic--B. kaiseri), threatened species growing in St. Catherine Protectorate, southern Sinai, Egypt. They are subjected to a number of threats that have caused populations to decline in both number and size. For the long-term survival of these species, an appropriate conservation strategy for the maintenance of their genetic variation should be developed. This study measures genetic diversity within and among populations of these Ballota species and determines the conservation implications of the results. The genetic analyses demonstrated that the three Ballota species maintain relatively high levels of genetic diversity (He = 0.195-0.317) and that most of the their genetic diversity was found within populations (GST = 0.045-0.099). Indirect estimates of historical gene flow for B. undulata and B. saxatilis were relatively high (Nm(W) = 5.25 and 3.37, respectively) but suggest that there is somewhat less gene movement among B. kaiseri populations (Nm(W) = 2.29). The levels of genetic diversity maintained within populations of the three Ballota species indicate that an appropriate sampling design for ex situ safeguarding should capture the majority of the genetic diversity found within these taxa.     

GENETIC VARIATION AND ESTIMATES OF GENE FLOW IN CLARKIA SPECIOSA SUBSP. POLYANTHA (ONAGRACEAE)

Estimates of interpopulational gene flow and the levels and distribution of genetic variation in Clarkia speciosa subsp. polyantha were obtained using enzyme electrophoresis. Eight enzymes encoded by 17 loci were analyzed. Nei's mean genetic identity was 0.96, indicating little genetic divergence among populations. Gene diversity statistics also suggest little heterogeneity among populations. Interpopulational gene flow, estimated according to Slatkin (1985), was fairly high, Nm = 3.9, probably accounting for the lack of differentiation among populations.     

Population genetic analysis of the genes APOE, APOB(3'VNTR) and ACE in some black and Amerindian communities from Colombia

A population genetic study was carried out with the APOE, APOB and ACE loci in 17 Colombian human populations. Ten of them were Amerindian communities coming from the northeastern part of Colombia, Pacific region, Eastern Plains and Amazonia. Six were black populations from Providence Island, Caribbean and Pacific coasts. Finally, the Mestizo population of Bogota was studied as well. The APOE and ACE loci were in Hardy-Weinberg equilibrium, whereas the APOB locus was not studied in all populations. The genetic heterogeneity was substantially greater among the Amerindian populations (G(ST) = 0.059) than in the Afrocolombian populations (G(ST) = 0.009). Also the gene flow population pair estimates were so much higher among the Afrocolombian populations (Nm = 49.08 +/- 43.07) than among Amerindian populations (Nm = 9.66 +/- 18.04). Different phylogenetic and multivariant analyses showed that the Amerindian populations analyzed were clustered in three different arrays: one constituted by the Colombian northeastern and Pacific populations, the second one by the two Amazon populations (Coreguaje and Nukak) and the last one by the Yuco (the unique Caribbe-speaking population among those studied). The latter population was highly divergent from a genetic point of view from the remainder Amerindian populations studied. By using the Mantel test, the existence of a positive and significant correlation between the genetic and geographical distances found among Amerindian populations was demonstrated. This fact was not observed among the Afrocolombian populations. Nevertheless, an isolation-by-distance Slatkin analysis test did not show a significant clear structure of this special pattern among the Indian tribes studied.     

Patterns of variation within and between Greek populations of Ceratitis capitata suggest extensive gene flow and latitudinal clines

Four natural Greek populations of Mediterranean fruit fly, Ceratitis capitata (Wiedemann), was studied for genetic variability at 25 enzyme loci. The comparison of polymorphism within and between populations shows two loci with high between-population heterozygosity (HT) and very high fixation index (F(ST)) values, suggesting the presence of balancing selection. The gradual decline of common allele frequency of the polymorphic loci tested indicated that latitudinal clines are present in Greece. Indirect estimates of gene flow based both on Wright's method (Nm*) and on the Slatkin's method (Nm*), which depends on the frequencies of rare alleles found in only one population, revealed a substantial amount of gene flow (Nm = 3.493 and Nm* = 3.197). These estimates of gene flow may well explain why the "introduced" Greek populations of C. capitata, in spite of their low genetic variability, display the same polymorphic loci. Gene flow in combination with natural selection and genetic drift may have played an important role to genetic differentiation in this species in Greece.     

A study of genetic distances and variability in several species pf the gemis Ctenomys (Rodentia: Octodontidae) with special references to a probable causal role of chromosomes in speciation

Genetic variability measured by allozymic electrophoresis has been studied in several species of the subterranean rodent genus Ctenomys (Octodontidae). The study was carried out with the main purpose of analysing a special group known as the 'Corrientes group' which inhabits the province of Corrientes in Argentina. The members of the group are, with high probability, isolated reproductively due to their karyotypic differences diploid numbers are between 42 and 70). To evaluate whether the chromosomal differences were the causes of speciation or if they arose a posteriori, we compared the level of the genetic distance among taxa within the orrientes group with all those measured between intraspecific populations in Ctenomys . The results indicate that the distances among the populations of the Corrientes group are at or below the level of those measured among traspecific populations in Ctenomys . Thus, it is unlikely that genetic differentiation triggered speciation. It is shown that this low level of genetic differentiation is not in contradiction with the high levels expected for species originating from one of the possible mechanisms of chromosomal speciation named, chromosomal transilience. Although for geographical reasons it seems obviousthat gene flow is precluded among the members of the Corrientes group, the Slatkin method for estimating Nm values was also used. Because high values of Nm exist but no isolation by distance could be detected, it is suggested that reproductive and geographical isolation are very ecent. Finally, analysis of population variability suggests that the high levels of heterozygosity observed (a) can be explained by the population structure, and (b) are within the range of expected values if bottlenecks have occurred in the recent history of the Corrientes group. As a general conclusion, the results indicate that in the Corrientes group the genetic data support a putative causal role for chromosomes in speciation.     

Oligoryzomys flavescens (Rodentia, Muridae): gene flow among populations from central-eastern Argentina

In species acting as hosts of infectious agents, the extent of gene flow between populations is of particular interest because the expansion of different infectious diseases is usually related to the dispersal of the host. We have estimated levels of gene flow among populations of the sigmodontine rodent Oligoryzomys flavescens, in which high titers of antibodies have been detected for a Hantavirus in Argentina that produces a severe pulmonary syndrome. Enzyme polymorphism was studied by means of starch gel electrophoresis in 10 populations from the area where human cases of Hantavirus have occurred. Genetic differentiation between populations was calculated from F_ST values with the equation Nm = [(1/F_ST−1]/4. To assess the relative importance of current gene flow and historical associations between populations, the relationship of population pairwise log Nm and log geographic distance was examined. Low F_ST (mean = 0.038) and high Nm (15.27) values suggest high levels of gene flow among populations. The lack of an isolation by distance pattern would indicate that this species has recently colonized the area. The northernmost population, located on the margin of a great river, shows very high levels of gene flow with the downstream populations despite the large geographic distances. Passive transport of animals down the river by floating plants would promote unidirectional gene flow. This fact and the highest mean heterozygosity of that northernmost population suggest it is a center of dispersal within the species' range. This revised version was published online in July 2006 with corrections to the Cover Date.     

RAPD variation and population genetic structure in Prunus mahaleb (Rosaceae), an animal-dispersed tree

We examined the patterns of random amplified polymorphic DNA (RAPD) variation among seven Prunus mahaleb (Rosaceae) populations extending over approximately 100 km2 to examine local differentiation in relation to spatial isolation due to both geographical distance and differences in elevation. No less than 51. 4% of the RAPD loci were polymorphic, but very few were fixed and among-population variation accounted for 16.46% of variation in RAPD patterns. Mean gene diversity was 0.1441, with mean Nei's genetic diversity for individual populations ranging between 0.089 and 0.149. Mean GST value across loci was 0.1935 (range, 0.0162-0.4685), giving an average estimate for Nm of 1.191. These results suggest extensive gene flow among populations, but higher GST and lower Nm values relative to other outcrossing, woody species with endozoochorous dispersal, also suggest a process of isolation by distance. The combined effect of both geographical and elevation distances and nonoverlapping flowering and fruiting phenophases on the GST matrix was partially significant, revealing only marginal isolation of the P. mahaleb populations. The matrix correlation between estimated Nm values among populations and the geographical + elevation distance matrices (r = -0.4623, P = 0.07), suggests a marginal trend for more isolated populations to exchange less immigrants. Long-distance seed dispersal by efficient medium-sized frugivorous birds and mammals is most likely associated to the high levels of within-population genetic diversity. However, vicariance factors and demographic bottlenecks (high postdispersal seed and seedling mortality) explain comparatively high levels of local differentiation.     

Genetic diversity and structure of Cordyceps sinensis populations from extensive geographical regions in China as revealed by inter-simple sequence repeat markers

Cordyceps sinensis is one of the most valuable medicinal caterpillar fungi native to China. However, its productivity is extremely limited and the species is becoming endangered. The genetic diversity of eighteen C. sinensis populations across its major distributing regions in China was evaluated by inter-simple sequence repeat (ISSR) markers. A total of 141 markers were produced in 180 individuals from the 18 populations, of which 99.3% were polymorphic. The low average of Shannon (0.104) and Nei index (0.07) of the 18 populations indicates that there are little genetic variations within populations. For all 18 populations, estimates of total gene diversity (H_T), gene diversity within populations (H_S), coefficient of genetic differentiation (G_ST), and gene flow (Nm) were 0.170, 0.071, 0.583, and 0.357, respectively. This pattern suggests that the genetic diversity of C. sinensis is low and most of the ISSR variations are found among populations with little gene exchange. The 18 populations are divided into five groups based on the genetic distance and the grouping pattern matches with the geographic distribution along the latitudinal gradient. The five groups show obvious difference in the G_ST and Nm values. Therefore, the genetic diversification of C. sinensis populations may be determined by geographic isolation and the combined effects of life history characters and the interaction with host insect species. The information illustrated by this study is useful for selecting in situ conservation sites of C. sinensis.     

EXTINCTION AND RECOLONIZATION: THEIR EFFECTS ON THE GENETIC DIFFERENTIATION OF LOCAL POPULATIONS

In this paper, we use a model by Slatkin (1977) to investigate the genetic effects of extinction and recolonization for a species whose population structure consists of an array of local demes with some migration among them. In particular, we consider the conditions under which extinction and recolonization might enhance or diminish gene flow and increase or decrease the rate of genetic differentiation relative to the static case with no extinctions. We explicitly take into account the age-structure that is established within the array of populations by the extinction and colonization process. We also consider two different models of the colonization process, the so-called “migrant pool” and “propagule pool” models. Our theoretical studies indicate that the genetic effects of extinction and colonization depend upon the relative magnitudes of K, the number of individuals founding new colonies, and 2Nm, twice the number of migrants moving into extant populations. We find that these genetic effects are surprisingly insensitive to the extinction rate. We conclude that, in order to assess the genetic effects of the population dynamics, we must first answer an important empirical question that is essentially ecological: is colonization a behavior distinct from migration?     

Estimates of gene flow, genetic substructure and population heterogeneity in bracken (Pteridium aquilinum)

Bracken [ Pteridium aquilinum (L.) Kuhn] is a cosmopolitan species and is a noxious weed in many areas. Because of its abundance, particularly in Britain, bracken affords an ideal system for investigating various aspects of population genetics and evolution. High mobility of dispersal units (spores) suggests that rates of gene flow among distant populations should be high. Gene flow is a major evolutionary force that influences the genetic structure of populations. To examine the effects of gene flow on population heterogeneity and population substructuring in bracken, starch gel electrophoresis of enzymes was used to provide the necessary genetic database. Allele frequency data at 21 loci were obtained for seven British populations, one Majorcan and one from the eastern United States. A model was employed to estimate the amount of gene flow ( Nm ) at several levels. Gene flow among British populations was extremely high ( Nm = 36.51), one of the highest estimates reported for plants. Among eight European populations gene flow was lower (but still considered high) at Nm = 2.47. Trans-Atlantic gene flow was low ( Nm = 0.0926).
F -statistics were used to assess population heterogeneity and substructuring. The data indicate that, compared with other species, there is very little genetic differentiation among British populations of bracken. Indeed, it appears that the whole island is behaving as a single randommating population. This result is consistent with high levels of gene flow. Only one population (on the Isle of Arran) showed statistically significant genetic substructuring. Habitat heterogeneity on the island and age structure are hypothesized as possible causes of this result.
The data reported here support previous studies demonstrating that bracken is genetically polymorphic and is an outcrossing species.     

Population genetic study of allozyme variation in natural populations of Drosophila antonietae (Insecta, Diptera)

Drosophila antonietae is an endemic South American cactophilic species found in relictual xerophytic vegetation, mostly associated with Cereus hildmaniannus cactus. Low differentiation among populations of this species has been detected using several markers. In this work, we performed an allozyme genetic variability analysis of 11 natural populations of D. antonietae and included a discussion about the possible influences of several evolutionary processes that might be acting to maintain the pattern observed. The genetic variability of 14 isoenzyme loci was analysed and showed a high genetic diversity (average observed heterozygosity = 0.319) and a moderate genetic differentiation among populations ( F statistics = 0.0723). A correlation between genetic and geographical and ecological distances was detected among pairs of populations and the regional equilibrium analysis was thus applied. This analysis resulted in Nm (number of migrants) of approximately 3.21, indicating that moderate levels of both gene flow and genetic drift occur in this species, with gene flow overlapping genetic drift. However, considering ecological features of drosophilids, we propose a hypothesis to explain the moderate differentiation encountered as a result of three different processes, or a combination of them: (1) gene flow; (2) a short period of differentiation, i.e. maintenance of ancestral polymorphism; and (3) action of natural selection. Moreover, if gene flow is present, the high genetic diversity compared with other cactophilic and non-cactophilic species could be due to differential selection in different populations followed by gene exchange among them. These factors are discussed in the light of D. antonietae 's historical and evolutionary association with the host cactus.     

Allopolyploid origin and population genetics of the rare orchid Spiranthes diluvialis

The process of becoming and the attributes of being polyploid play a major role in the development and maintenance of genetic variation in allopolyploid species. A genetic survey employing protein electrophoresis on 12 populations of S. diluvialis, as well as on populations of eight congeneric species, was conducted to assess the putative allopolyploid origin of S. diluvialis and to determine the genetic variability within and among populations. Genetic identity values indicated S. diluvialis was more similar to S. magnicamporum (0.619) and S. romanzoffiana (0.727) than to any of the other congeneric species assayed. Similar to most allopolyploids, S. diluvialis showed high levels of fixed, or nearly fixed, heterozygosity and a high percentage of polymorphic loci (57.1-71.4%). The mean number of alleles per polymorphic locus in populations of S. diluvialis (2.6-3.3), however, was similar to mean values for both animal-pollinated, outcrossing, diploid species, and geographically restricted, diploid species (2.6 and 2.5, respectively). Genetic divergence among populations (mean Fst = 0.083) was low, leading to relatively high estimates of interpopulational gene flow (mean Nm = 5.41). Thus, each population harbors most of the genetic variability found within the species. The genetic variation observed within S. diluvialis supports the occurrence of at least two separate hybridization events giving rise to S. diluvialis.     

Nuclear and chloroplast microsatellites reveal extreme population differentiation and limited gene flow in the Aegean endemic Brassica cretica (Brassicaceae)

Nuclear and chloroplast microsatellite markers were used to study population structure and gene flow among seven Cretan populations of the Aegean endemic plant species Brassica cretica (Brassicaceae). Both nuclear and chloroplast markers revealed exceptionally high levels of population differentiation (overall F(ST)=0.628 and 1.000, respectively) and relatively little within-population diversity (overall H(S)=0.211 and 0.000, respectively). Maximum-likelihood estimates of directional migration rates were low among all pairs of populations (average Nm=0.286). There was no evidence that differences in flower colour between populations had any influence on historical levels of gene flow. In addition, a haplotype network showed that all five chloroplast haplotypes found in the sample were closely related. Together, these results suggest that current patterns of diversification in B. cretica are mainly a result of genetic drift during the last half million years. The main conclusions from the present study are consistent with the prevailing hypothesis that plant diversification in the Aegean region is driven by random rather than adaptive differentiation among isolated populations.     

Gene flow and effective population size in Lepanthes (Orchidaceae): a case for genetic drift

Genetic drift can play an important role in population differentiation, particularly when effective population sizes are small and gene flow is limited. Such conditions are suspected to be common in the species-rich Orchidaceae. We investigated the likelihood of genetic drift in natural populations of three endemic species of Lepanthes (Orchidaceae) from Puerto Rico. We estimated effective population size, Ne, using three ecologically based methods. Two of the three estimates were based on variance in reproductive potential and the third was based on coalescence time. All estimates of Ne were usually <40% of the standing population size, resulting in values of <20 individuals per population. Based on starch gel electrophoresis of isozymes, Nm estimates suggest restricted gene flow among populations in the range of one or less successful migrant per generation. Genetic differentiation among populations is expected under such conditions from random genetic drift. Indeed we observed high genetic differentiation among populations (L. rubripetala, F_ST> G_ST, 6; θ.248, 0.266, 0.293; L. rupestris, 0.148, 0.169, 0.138; L. eltoroensis, 0.251, 0.219, 0.218, respectively). Genetic drift is likely to be important for population differentiation in Lepanthes as a result of small effective population sizes and restricted gene flow.     

A standardized method for quantifying unidirectional genetic introgression

Genetic introgression of domesticated to wild conspecifics is of great concern to the genetic integrity and viability of the wild populations. Therefore, we need tools that can be used for monitoring unidirectional gene flow from domesticated to wild populations. A challenge to quantitation of unidirectional gene flow is that both the donor and the recipient population may be genetically substructured and that the subpopulations are subjected to genetic drift and may exchange migrants between one another. We develop a standardized method for quantifying and monitoring domesticated to wild gene flow and demonstrate its usefulness to farm and wild Atlantic salmon as a model species. The challenge of having several wild and farm populations was circumvented by in silico generating one analytical center point for farm and wild salmon, respectively. Distributions for the probability that an individual is wild were generated from individual‐based analyses of observed wild and farm genotypes using STRUCTURE. We show that estimates of proportions of the genome being of domesticated origin in a particular wild population can be obtained without having a historical reference sample for the same population. The main advantages of the method presented are the standardized way in which genetic processes within and between populations are taken into account, and the individual‐based analyses giving estimates for each individual independent of other individuals. The method makes use of established software, and as long as genetic markers showing generic genetic differences between domesticated and wild populations are available, it can be applied to all species with unidirectional gene flow. Results from our method are easy to interpret and understand, and will serve as a powerful tool for management, especially because there is no need for a specific historical wild reference sample.     

七筋菇自然居群的遗传结构分析

采用ISSR分子标记,对七筋菇(Clintonia udensis)17个居群的遗传多样性与遗传结构进行了研究。结果表明:七筋菇不同居群的多态位点百分率PPB为11.90%~59.52%,总的多态位点百分率PPB为98.8%,具有高的遗传多样性。Shannon多样性指数(0.6903)和基因分化系数(GST=0.6944)均揭示出七筋菇居群间存在明显的遗传差异,AMOVA分析结果也显示遗传变异主要发生在居群之间(81.47%),而居群内部的遗传变异仅为18.53%。七筋菇居群间的遗传距离从0.1871~0.6632,平均为0.3838,大于同一物种居群间的平均遗传距离值(0.05),同样表明七筋菇居群间的遗传多样性存在较大差异。七筋菇居群间的基因流Nm=0.2200,远远低于一般广布种植物的基因流(Nm=1.881)。Mantel检测显示居群间的遗传距离与地理距离之间没有显著相关性(r=0.029,P=0.3196)。七筋菇分布范围广以及其进化历史是其具有高遗传多样性的原因;居群间存在较高遗传变异可能是由于七筋菇本身的生物学特性、有限的基因流以及遗传漂变等原因造成的。     

Speciation despite gene flow when developmental pathways evolve

Abstract.— Evolutionary biologists assume that species formation requires a drastic reduction in gene exchange between populations, but the rate sufficient to prevent speciation is unknown. To study speciation, we use a new class of population genetic models that incorporate simple developmental genetic rules, likely present in all organisms, to construct the phenotype. When we allow replicate populations to evolve in parallel to a new, shared optimal phenotype, often their hybrids acquire poorly regulated phenotypes: Dobzhansky-Muller incompatibilities arise and postzygotic reproductive isolation evolves. Here we show that, although gene exchange does inhibit this process, it is the proportion of migrants exchanged ( m ) rather than the number of migrants ( Nm ) that is critical, and rates as high as 16 individuals exchanged per generation still permit the evolution of postzygotic isolation. Stronger directional selection counters the inhibitory effect of gene flow, increasing the speciation probability. We see similar results when populations in a standard two-locus, two-allele Dobzhansky-Muller model are subject to simultaneous directional selection and gene flow. However, in developmental pathway models with more than two loci, gene flow is more able to impede speciation. Genetic incompatibilities arise as frequent by-products of adaptive evolution of traits determined by regulatory pathways, something that does not occur when phenotypes are modeled using the standard, additive genetic framework. Development therefore not only constrains the microevolutionary process, it also facilitates the interactions among genes and gene products that make speciation more likely–even in the face of strong gene flow.