Review. The strength and genetic basis of reproductive isolating barriers in flowering plants

Speciation is characterized by the evolution of reproductive isolation between two groups of organisms. Understanding the process of speciation requires the quantification of barriers to reproductive isolation, dissection of the genetic mechanisms that contribute to those barriers and determination of the forces driving the evolution of those barriers. Through a comprehensive analysis involving 19 pairs of plant taxa, we assessed the strength and patterns of asymmetry of multiple prezygotic and postzygotic reproductive isolating barriers. We then reviewed contemporary knowledge of the genetic architecture of reproductive isolation and the relative role of chromosomal and genic factors in intrinsic postzygotic isolation. On average, we found that prezygotic isolation is approximately twice as strong as postzygotic isolation, and that postmating barriers are approximately three times more asymmetrical in their action than premating barriers. Barriers involve a variable number of loci, and chromosomal rearrangements may have a limited direct role in reproductive isolation in plants. Future research should aim to understand the relationship between particular genetic loci and the magnitude of their effect on reproductive isolation in nature, the geographical scale at which plant speciation occurs, and the role of different evolutionary forces in the speciation process.

Speciation with gene flow in the large white-headed gulls: does selection counterbalance introgression?

We investigated the role of selection in generating and maintaining species distinctness in spite of ongoing gene flow, using two zones of secondary contact between large gull species in Europe (Larus argentatus and Larus cachinnans) and North America (Larus glaucescens and Larus occidentalis). We used the pattern of neutral genetic differentiation at nine microsatellite loci (F(ST)) as an indicator of expected changes under neutral processes and compared it with phenotypic differentiation (P(ST)) for a large number of traits (size, plumage melanism and coloration of bare parts). Even assuming very low heritability, interspecific divergence between L. glaucescens and L. occidentalis in plumage melanism and orbital ring colour clearly exceeded neutral differentiation. Similarly, melanism of the central primaries was highly divergent between L. argentatus and L. cachinnans. Such divergence is unlikely to have arisen randomly and is therefore attributed to spatially varying selection. Variation in plumage melanism in both transects agrees with Gloger's rule, which suggests that latitude (and associated sun and humidity gradients) could be the selective pressure shaping differentiation in plumage melanism. We suggest that strong species differentiation in orbital ring colour results from sexual selection. We conclude that these large gull species, along with other recently diverged species that hybridize after coming into secondary contact, may differ only in restricted regions of the genome that are undergoing strong disruptive selection because of their phenotypic effects.

Searching for candidate speciation genes using a proteomic approach: seminal proteins in field crickets

In many animals, male seminal proteins influence gamete interactions and fertilization ability and are probably involved in barriers to gene flow between diverging lineages. Here we use a proteomic approach to identify seminal proteins that are transferred to females during copulation and that may be involved in fertilization barriers between two hybridizing field crickets (Gryllus firmus and Gryllus pennsylvanicus). Analyses of patterns of divergence suggest that much of the field cricket genome has remained undifferentiated following the evolution of reproductive isolation. By contrast, seminal protein genes are highly differentiated. Tests of selection reveal that positive selection is likely to be responsible for patterns of differentiation. Together, our observations suggest that some of the loci encoding seminal proteins may indeed play a role in fertilization barriers in field crickets.

Mean time to resolution of gene duplication

The mean time to resolution of gene duplication (T(r)) is studied in this paper under the double null recessive (DNR) and haplo-insufficient (HI) models within the same analytical and simulation framework. We show that when population size is not too small (more precisely Nmu > 0.1), T(r) for unlinked duplication is usually larger than that for linked and T(r) for unlinked duplication under the HI model might be greatly prolonged, which were consistent with previous observations. Furthermore, by analytical approach we here indicate the primary underlying mechanism is that the frequency of the original (or wild-type) chromosomal haplotype of the linked duplication decreases nearly exponential to zero with time while that of the unlinked decreases quickly to an quasi-equilibrium; and this phenomenon is particularly profound under the HI model, because the quasi-equilibrium frequency of the original chromosomal haplotype (x(0)) under the HI model is higher than that under the DNR model. These results suggest that recombination and HI model might jointly contribute to the marked prolongation of T(r) even in a modest population. The prolonged T(r) and higher quasi-equilibrium frequency of the original allele at both duplicated loci might have offered more opportunities for the emergence of novel genes.

Ecological reproductive isolation of coast and inland races of Mimulus guttatus

Adaptive divergence due to habitat differences is thought to play a major role in formation of new species. However it is rarely clear the extent to which individual reproductive isolating barriers related to habitat differentiation contribute to total isolation. Furthermore, it is often difficult to determine the specific environmental variables that drive the evolution of those ecological barriers, and the geographic scale at which habitat-mediated speciation occurs. Here, we address these questions through an analysis of the population structure and reproductive isolation between coastal perennial and inland annual forms of the yellow monkeyflower, Mimulus guttatus. We found substantial morphological and molecular genetic divergence among populations derived from coast and inland habitats. Reciprocal transplant experiments revealed nearly complete reproductive isolation between coast and inland populations mediated by selection against immigrants and flowering time differences, but not postzygotic isolation. Our results suggest that selection against immigrants is a function of adaptations to seasonal drought in inland habitat and to year round soil moisture and salt spray in coastal habitat. We conclude that the coast and inland populations collectively comprise distinct ecological races. Overall, this study suggests that adaptations to widespread habitats can lead to the formation of reproductively isolated species.

Recurrent positive selection of the Drosophila hybrid incompatibility gene Hmr

Lethality in hybrids between Drosophila melanogaster and its sibling species Drosophila simulans is caused in part by the interaction of the genes Hybrid male rescue (Hmr) and Lethal hybrid rescue (Lhr). Hmr and Lhr have diverged under positive selection in the hybridizing species. Here we test whether positive selection of Hmr is confined only to D. melanogaster and D. simulans. We find that Hmr has continued to diverge under recurrent positive selection between the sibling species D. simulans and Drosophila mauritiana and along the lineage leading to the melanogaster subgroup species pair Drosophila yakuba and Drosophila santomea. Hmr encodes a member of the Myb/SANT-like domain in ADF1 (MADF) family of transcriptional regulators. We show that although MADF domains from other Drosophila proteins have predicted ionic properties consistent with DNA binding, the MADF domains encoded by different Hmr orthologs have divergent properties consistent with binding to either the DNA or the protein components of chromatin. Our results suggest that Hmr may be functionally diverged in multiple species.

The molecular forms of Anopheles gambiae: a phenotypic perspective

The African malaria mosquito Anopheles gambiae is undergoing speciation, being split into the M and S molecular forms. Speciation is the main process promoting biological diversity, thus, new vector species might complicate disease transmission. Genetic differentiation between the molecular forms has been extensively studied, but phenotypic differences between them, the evolutionary forces that generated divergence, and the mechanisms that maintain their genetic isolation have only recently been addressed. Here, we review recent studies suggesting that selection mediated by larval predation and competition promoted divergence between temporary and permanent freshwater habitats. These differences explain the sharp discontinuity in distribution of the molecular forms between rice fields and surrounding savanna, but they can also explain the concurrent cline between humid and arid environments due to the dependence on permanent habitats in the latter. Although less pronounced, differences in adult body size, reproductive output, and longevity also suggest that the molecular forms have adapted to distinct niches. Reproductive isolation between the molecular forms is achieved by spatial swarm segregation, although within-swarm mate recognition appears to play a role in certain locations. The implications of these results to disease transmission and control are discussed and many of the gaps in our understanding are highlighted.

Molecular demographic history of the annual sunflowers Helianthus annuus and H. petiolaris--large effective population sizes and rates of long-term gene flow

Hybridization between distinct species may lead to introgression of genes across species boundaries, and this pattern can potentially persist for extended periods as long as selection at some loci or genomic regions prevents thorough mixing of gene pools. However, very few reliable estimates of long-term levels of effective migration are available between hybridizing species throughout their history. Accurate estimates of divergence dates and levels of gene flow require data from multiple unlinked loci as well as an analytical framework that can distinguish between lineage sorting and gene flow and incorporate the effects of demographic changes within each species. Here we use sequence data from 18 anonymous nuclear loci in two broadly sympatric sunflower species, Helianthus annuus and H. petiolaris, analyzed within an "isolation with migration" framework to make genome-wide estimates of the ages of these two species, long-term rates of gene flow between them, and effective population sizes and historical patterns of population growth. Our results indicate that H. annuus and H. petiolaris are approximately one million years old and have exchanged genes at a surprisingly high rate (long-term N(ef)m estimates of approximately 0.5 in each direction), with somewhat higher rates of introgression from H. annuus into H. petiolaris than vice versa. In addition, each species has undergone dramatic population expansion since divergence, and both species have among the highest levels of genetic diversity reported for flowering plants. Our results provide the most comprehensive estimate to date of long-term patterns of gene flow and historical demography in a nonmodel plant system, and they indicate that species integrity can be maintained even in the face of extensive gene flow over a prolonged period.

Evidence for contrasting modes of selection at interacting globin genes in the European rabbit (Oryctolagus cuniculus)

In hybrid zones between genetically differentiated populations, variation in locus-specific rates of introgression may reflect adaptation to different environments or adaptation to different genetic backgrounds. The European rabbit, Oryctolagus cuniculus, is well-suited to studies of such hybrid zone dynamics because it is composed of two genetically divergent subspecies that hybridize in a zone of secondary contact in central Iberia. A species-wide survey of allozyme variation revealed a broad range of locus-specific divergence levels (F(ST) ranged from 0 to 0.54, mean F(ST)=0.16). Interestingly, the two loci that fell at opposite ends of the distribution of F(ST) values, haemoglobin alpha-chain (HBA) and haemoglobin beta-chain (HBB), encode interacting subunits of the haemoglobin protein. The contrasting patterns of spatial variation at these two loci could not be reconciled under a neutral model of population structure. The HBA gene exhibited higher-than-expected levels of population differentiation, consistent with a history of spatially varying selection. The HBB gene exhibited lower-than-expected levels of population differentiation, consistent with some form of spatially uniform selection. Patterns of linkage disequilibrium and allele frequency variation do not appear to fit any simple model of two-locus epistatic selection.

Population genetic structure across dissolved oxygen regimes in an African cichlid fish

Ecological isolation is a process whereby gene flow between selective environments is reduced due to selection against maladapted dispersers, migrant alleles, or hybrids. Although ecological isolation has been documented in several systems, gene flow can often be high among selective regimes, and more studies are thus needed to better understand the conditions under which ecological gradients or divergent selective regimes should influence population structure. We test for ecological isolation in a system in which high plasticity occurs with respect to traits that are adaptive in alternate forms under different environmental conditions. Pseudocrenilabrus multicolor victoriae is a widespread haplochromine cichlid fish in East Africa that exploits both normoxic (normal oxygen) rivers/lakes and hypoxic (low oxygen) swamps. Here, we examine population structure, using mitochondrial DNA and microsatellites, to determine if genetic divergence is significantly increased between dissolved oxygen regimes relative to within them, while controlling for geographical structure. Our results indicate that geographical separation influences population structure, while no effects of divergent selection with respect to oxygen regimes were detected. Specifically, we document (i) genetic clustering according to geographical region, but no clustering according to oxygen regime; (ii) higher genetic variation among than within regions, but no effect of oxygen regime on genetic variation; (iii) isolation by distance within one region; and (iv) decreasing genetic variability with increasing geographical distance from Lake Victoria. We speculate that plasticity may be facilitating gene flow between oxygen regimes in this system.

Reproductive isolation among two interfertile Rhododendron species: low frequency of post-F1 hybrid genotypes in alpine hybrid zones

Hybrids between the acid-loving species Rhododendron ferrugineum and the basic soil species Rhododendron hirsutum occur on soils of intermediate pH in the European Alps. Material from two hybrid zones approximately 500 m apart, and also nearby populations of each parent species, was surveyed for presence/absence of 31 random amplified polymorphic DNA markers that distinguish parents. Based on morphological assessment, the material comprised 51 putative hybrids, 18 putative R. ferrugineum individuals and 26 putative R. hirsutum plants. RAPD data were analysed using a Bayesian approach implemented by the program newhybrids, and also by principal coordinates analysis. The identity of all R. ferrugineum plants examined was confirmed; however, of the putative R. hirsutum individuals examined, two were certainly and 11 possibly hybrid derivatives. Among all hybrid derivatives examined, about half were designated as F1s or a similar class, otherwise backcrosses to R. hirsutum appeared to be common whereas other hybrid classes were rare and backcrosses to R. ferrugineum possibly absent. Despite this, artificially generated seed of F2 class and backcrosses in each direction showed greater viability than one parent (R. hirsutum). Introgression from R. ferrugineum was also detected in a population that from morphology appeared to contain only R. hirsutum. Hence, the direction of backcrossing might be highly asymmetric within hybrid zones, causing unidirectional gene flow from R. ferrugineum into R. hirsutum. Conversely, the rarity of backcrosses to R. ferrugineum, F2s and later hybrid generations, which might be due to phenology effects and habitat-mediated selection, could play a part in restricting gene flow towards R. ferrugineum.

Ecological niche dimensionality and the evolutionary diversification of stick insects

The degree of phenotypic divergence and reproductive isolation between taxon pairs can vary quantitatively, and often increases as evolutionary divergence proceeds through various stages, from polymorphism to population differentiation, ecotype and race formation, speciation, and post-speciational divergence. Although divergent natural selection promotes divergence, it does not always result in strong differentiation. For example, divergent selection can fail to complete speciation, and distinct species pairs sometimes collapse (‘speciation in reverse’). Widely-discussed explanations for this variability concern genetic architecture, and the geographic arrangement of populations. A less-explored possibility is that the degree of phenotypic and reproductive divergence between taxon pairs is positively related to the number of ecological niche dimensions (i.e., traits) subject to divergent selection. Some data supporting this idea stem from laboratory experimental evolution studies using Drosophila, but tests from nature are lacking. Here we report results from manipulative field experiments in natural populations of herbivorous Timema stick insects that are consistent with this ‘niche dimensionality’ hypothesis. In such insects, divergent selection between host plants might occur for cryptic colouration (camouflage to evade visual predation), physiology (to detoxify plant chemicals), or both of these niche dimensions. We show that divergent selection on the single niche dimension of cryptic colouration can result in ecotype formation and intermediate levels of phenotypic and reproductive divergence between populations feeding on different hosts. However, greater divergence between a species pair involved divergent selection on both niche dimensions. Although further replication of the trends reported here is required, the results suggest that dimensionality of selection may complement genetic and geographic explanations for the degree of diversification in nature.

Ancient hybridization and mitochondrial capture between two species of chipmunks

Models that posit speciation in the face of gene flow are replacing classical views that hybridization is rare between animal species. We use a multilocus approach to examine the history of hybridization and gene flow between two species of chipmunks (Tamias ruficaudus and T. amoenus). Previous studies have shown that these species occupy different ecological niches and have distinct genital bone morphologies, yet appear to be incompletely isolated reproductively in multiple areas of sympatry. We compared data from four sequenced nuclear loci and from seven microsatellite loci to published cytochrome b sequences. Interspecific gene flow was primarily restricted to introgression of the T. ruficaudus mitochondrial genome into a sympatric subspecies of T. amoenus, T. a. canicaudus, with the four sequenced nuclear loci showing little to no interspecific allele sharing. Microsatellite data were consistent with high levels of differentiation between the species and also showed no current gene flow between broadly sympatric populations of T. a. canicaudus and T. ruficaudus. Coalescent analyses date the mtDNA introgression event from the mid-Pleistocene to late Pliocene. Overall, these data indicate that introgression has had a minimal impact on the nuclear genomes of T. amoenus and T. ruficaudus despite multiple independent hybridization events. Our findings challenge long-standing assumptions on patterns of reproductive isolation in chipmunks and suggest that there may be other examples of hybridization among the 23 species of Tamias that occur in western North America.

A population genomic analysis of species boundaries: neutral processes, adaptive divergence and introgression between two hybridizing plant species

Little is known about the nature of species boundaries between closely related plant species and about the extent of introgression as a consequence of hybridization upon secondary contact. To address these topics we analyzed genome-wide differentiation between two closely related Silene species, Silene latifolia and S. dioica, and assessed the strength of introgression in sympatry. More than 300 AFLP markers were genotyped in three allopatric and three sympatric populations of each species. Outlier analyses were performed separately for sympatric and allopatric populations. Both positive and negative outlier loci were found, indicating that divergent and balancing selection, respectively, have shaped patterns of divergence between the two species. Sympatric populations of the two species were found to be less differentiated genetically than allopatric populations, indicating that hybridization has led to gene introgression. We conclude that differentiation between S. latifolia and S. dioica has been shaped by a combination of introgression and selection. These results challenge the view that species differentiation is a genome-wide phenomenon, and instead support the idea that genomes can be porous and that species differentiation has a genic basis.

To be or not to be a hamlet pair in sympatry

Recent adaptive radiations, such as the fringillid finches on the Galapagos archipelago or the cichlid fishes in the East African Great Lakes, are invaluable model systems for evolutionary and ecological research. Puebla et al., in this issue of Molecular Ecology, have established a group of colorful coral reef fishes from the Caribbean sea for studying the early phases of species formation in a marine adaptive radiation. It appears that local evolutionary processes are important in this system, which might have even triggered in situ speciation.

Population structure and genetic diversity in two species of Hawaiian picture-winged Drosophila

Over the last several decades many picture-winged Drosophila have become less common in both geographical distribution and local population size (pers. obs., Foote pers. comm., Montgomerey pers. comm.). Here we report on a study of two Hawaiian Drosophila species, D. engyochracea, and D. hawaiiensis, to determine the impact that changes in population sizes over the past thirty years have had on the genetic diversity of these species. D. engyochracea is known from only two locations on the Island of Hawai'i (Kipuka Ki and Kipuka Pua'ulu), while D. hawaiiensis is currently more wide spread across Hawai'i Island. We collected 65 D. hawaiiensis and 66 D. engyochracea from two forest patches (kipuka) isolated by a 400 year old volcanic ash deposit. DNA sequence data for 515 bases of the mitochondrial gene COII was analyzed for both species to estimate relative total genetic diversity as well as inter-kipuka gene flow. The more wide spread species, D. hawaiiensis, has more genetic diversity (23 vs. 11 unique haplotypes) than the rarer species, D. engyochracea. The distribution of haplotypes in the kipuka is consistent with more gene flow in D. engyochracea than in D. hawaiiensis. Phylogenetic analysis indicates a small number of individuals morphologically identified as one species but have DNA sequence diagnostic for the other species. These results are consistent with these individuals being descendant from hybrids between species.

Oxidative damage in gill of Mytilus edulis from Merseyside, UK, and reversibility after depuration

Mussels were collected from the urban/industrialized site of New Brighton, Merseyside and the relatively non-industrial site of Llandudno, North Wales. All mussels were identified as Mytilus edulis by PCR amplification of Mefp1. DNA single strand breaks and 8-oxo-7,8-dihydro-2'-deoxyguanosine were measured in gill within 24h of collection, using the COMET assay, both with and without formamidopyrimidine glycosylase. Gill lipid peroxidation was also measured within 24h. No difference between sites was found for frank SSB and malonaldehyde levels, however 8-oxo-dG and 4-hydroxynonenal were significantly greater in New Brighton mussels compared to Llandudno mussels. After 1-month laboratory maintenance, lipid peroxidation and 8-oxo-dG levels were lower. In contrast, frank SSB were higher. This could reflect enhanced DNA repair excision, though we cannot exclude the possibility of other non-oxidative DNA damage. The results suggest that laboratory maintenance allows recovery from environmentally induced oxidative damage, which was more extensive at Merseyside.

Conservation implications of hybridization in Hawaiian picture-winged Drosophila

In this review, we discuss the importance of hybridization among species for the conservation of Hawaiian picture-winged Drosophila. Hybridization can be a positive evolutionary process that creates new species and increases the adaptation of populations and species through the spread of adaptive alleles and traits. Conversely, hybridization can disrupt the genetic integrity of species or populations and this may be most detrimental among taxa that are recently hybridizing due to recent ecological changes. The loss of biodiversity in Hawaiian Drosophila through hybridization may be facilitated by habitat alteration and introduced species that reduce population sizes and alter geographic distributions of native species. We briefly review the evidence for hybridization in the genus Drosophila and then focus on hybridization in the Hawaiian picture-winged Drosophila. We examine three general approaches for identifying hybrids and for assessing the factors that appear to contribute to hybridization and the potential ecological and evolutionary outcomes of hybrids in natural populations. Overall, the potential for hybridization among species will likely increase the risk of extinction for Hawaiian picture-winged Drosophila species. Thus, it is important to consider the potential for hybridization among species when developing plans for the conservation of Hawaiian Drosophila.

Heterogeneous genomic differentiation between walking-stick ecotypes: "isolation by adaptation" and multiple roles for divergent selection

Genetic differentiation can be highly variable across the genome. For example, loci under divergent selection and those tightly linked to them may exhibit elevated differentiation compared to neutral regions. These represent "outlier loci" whose differentiation exceeds neutral expectations. Adaptive divergence can also increase genome-wide differentiation by promoting general barriers to neutral gene flow, thereby facilitating genomic divergence via genetic drift. This latter process can yield a positive correlation between adaptive phenotypic divergence and neutral genetic differentiation (described here as "isolation-by-adaptation"). Here, we examine both these processes by combining an AFLP genome scan of two host plant ecotypes of Timema cristinae walking-sticks with existing data on adaptive phenotypic divergence and ecological speciation in these insects. We found that about 8% of loci are outliers in multiple population comparisons. Replicated comparisons between population-pairs using the same versus different host species revealed that 1-2% of loci are subject to host-related selection specifically. Locus-specific analyses revealed that up to 10% of putatively neutral (nonoutlier) AFLP loci exhibit significant isolation-by-adaptation. Our results suggest that selection may affect differentiation directly, via linkage, or by facilitating genetic drift. They thus illustrate the varied and sometimes nonintuitive contributions of selection to heterogeneous genomic differentiation.

Quantifying the constraining influence of gene flow on adaptive divergence in the lake-stream threespine stickleback system

The constraining effect of gene flow on adaptive divergence is often inferred but rarely quantified. We illustrate ways of doing so using stream populations of threespine stickleback (Gasterosteus aculeatus) that experience different levels of gene flow from a parapatric lake population. In the Misty Lake watershed (British Columbia, Canada), the inlet stream population is morphologically divergent from the lake population, and presumably experiences little gene flow from the lake. The outlet stream population, however, shows an intermediate phenotype and may experience more gene flow from the lake. We first used microsatellite data to demonstrate that gene flow from the lake is low into the inlet but high into the outlet, and that gene flow from the lake remains relatively constant with distance along the outlet. We next combined gene flow data with morphological and habitat data to quantify the effect of gene flow on morphological divergence. In one approach, we assumed that inlet stickleback manifest well-adapted phenotypic trait values not constrained by gene flow. We then calculated the deviation between the observed and expected phenotypes for a given habitat in the outlet. In a second approach, we parameterized a quantitative genetic model of adaptive divergence. Both approaches suggest a large impact of gene flow, constraining adaptation by 80-86% in the outlet (i.e., only 14-20% of the expected morphological divergence in the absence of gene flow was observed). Such approaches may be useful in other taxa to estimate how important gene flow is in constraining adaptive divergence in nature.

Hybridization and introgression in a mixed population of the intertidal seaweeds Fucus evanescens and F. serratus

The introduced Fucus evanescens (hermaphroditic) and the native F. serratus (dioecious) have been in secondary contact along the Danish coast of the Kattegat Sea for 60-100 years and dioecious hybrids have been observed at Blushøj for several years. Hybridization in Fucus is unusual because it appears to always involve a hermaphroditic and a dioecious parental pair. We determined the degree and spatial patterns of introgression for 286 individuals using 10 microsatellite loci and cpDNA. Hybrids accounted for nearly 13% of the population, yet parental species were well differentiated (F(ST) = 0.633). The presence of F. evanescens chloroplasts in 100% of F(1) hybrids revealed asymmetrical hybridization. Fucus evanescens cpDNA was observed in 50% of introgressed and 5.4% of pure F. serratus, but no F. serratus cpDNA was found in F. evanescens. In contrast, nuclear DNA introgression was symmetrical with an equal amount ( approximately 1.5%) of genes introgressed into each parental species. Survivorship and viability data suggest selection against hybrids in the hybrid zone.

Patterns of reproductive isolation in Mediterranean deceptive orchids

The evolution of reproductive isolation is of central interest in evolutionary biology. In plants, this is typically achieved by a combination of pre- and postpollination mechanisms that prevent, or limit, the amount of interspecific gene flow. Here, we investigated and compared two ecologically defined groups of Mediterranean orchids that differ in pollination biology and pollinator specificity: sexually deceptive orchids versus food-deceptive orchids. We used experimental crosses to assess the strength of postmating prezygotic, and postzygotic reproductive isolation, and a phylogenetic framework to determine their relative rates of evolution. We found quantitative and qualitative differences between the two groups. Food-deceptive orchids have weak premating isolation but strong postmating isolation, whereas the converse situation characterizes sexually deceptive orchids. Only postzygotic reproductive isolation among food-deceptive orchids was found to evolve in a clock-like manner. Comparison of evolutionary rates, within a common interval of genetic distance, showed that the contribution of postmating barriers was more relevant in the food-deceptive species than in the sexually deceptive species. Asymmetry in prezygotic isolation was found among food-deceptive species. Our results indicate that postmating barriers are most important for reproductive isolation in food-deceptive orchids, whereas premating barriers are most important in sexually deceptive orchids. The different rate of evolution of reproductive isolation and the relative strength of pre- and postmating barriers may have implication for speciation processes in the two orchid groups.

Population genetics of speciation in nonmodel organisms: I. Ancestral polymorphism in mangroves

The level of DNA polymorphism in the ancestral species at the time of speciation can be estimated using DNA sequences from many loci sampled from 2 or more extant species. The comparison between ancestral and extant polymorphism can be informative about the population genetics of speciation. In this study, we collected and analyzed DNA sequences of approximately 60 genes from 4 species of Sonneratia, a common genus of mangroves on the Indo-Pacific coasts. We found that the 3 ancestral species were comparable to each other in terms of level of polymorphism. However, the ancestral species at the time of speciation were substantially more polymorphic than the extant geographical populations. This ancestral polymorphism is in fact larger than, or at least equal to, the level of polymorphism of the entire species across extant geographical populations. The observations are not fully compatible with speciation by strict allopatry. We suggest that, at the time of speciation, the ancestral species consisted of interconnected but strongly divided geographical populations. This population structure would give rise to high level of polymorphism across species range. This approach of studying the speciation history by genomic means should be applicable to nonmodel organisms.

Diversity and divergence patterns in regulatory genes suggest differential gene flow in recently derived species of the Hawaiian silversword alliance adaptive radiation (Asteraceae)

The impact of gene flow and population size fluctuations in shaping genetic variation during adaptive radiation, at both the genome-wide and gene-specific levels, is very poorly understood. To examine how historical population size and gene flow patterns within and between loci have influenced lineage divergence in the Hawaiian silversword alliance, we have investigated the nucleotide sequence diversity and divergence patterns of four floral regulatory genes (ASAP1-A, ASAP1-B, ASAP3-A, ASAP3-B) and a structural gene (ASCAB9). Levels and patterns of molecular divergence across these five nuclear loci were estimated between two recently derived species (Dubautia ciliolata and Dubautia arborea) which are presumed to be sibling species. This multilocus analysis of genetic variation, haplotype divergence and historical demography indicates that population expansion and differential gene flow occurred subsequent to the divergence of these two lineages. Moreover, contrasting patterns of allele- sharing for regulatory loci vs. a structural locus between these two sibling species indicate alternative histories of genetic variation and partitioning among loci where alleles of the floral regulatory loci are shared primarily from D. arborea to D. ciliolata and alleles of the structural locus are shared in both directions. Taken together, these results suggest that adaptively radiating species can exhibit contrasting allele migration rates among loci such that allele movement at specific loci may supersede genetic divergence caused by drift and that lineage divergence during adaptive radiation can be associated with population expansion.

Influences of gene flow on adaptive speciation in theDubautia arborea-D. ciliolatacomplex

Mechanisms of reproductive isolation during plant speciation are often unclear because distinct species often experience high levels of gene flow and hybridization. Adaptive radiations such as the Hawaiian silversword alliance (HSA) provide unique opportunities to study the interactions of selection, gene flow and isolating mechanisms during the speciation process. We examined patterns of phenotypic and genetic differentiation in Dubautia arborea and Dubautia ciliolata, two parapatric HSA taxa that show marked morphological divergence but evidence of weak molecular differentiation, in order to estimate genome-wide differentiation and gene flow patterns. We scored 166 amplified fragment length polymorphism markers in a set of 89 plants from two populations each of D. arborea and D. ciliolata and phenotypically D. arborea-like and D. ciliolata-like plants from a natural hybrid zone. Analyses of population subdivision showed low levels of differentiation between the two species (F(ST) = 0.089) and evidence that the phenotypically parental hybrid zone plants were largely of parental species rather than of hybrid origin. A Bayesian analysis of population ancestry identified a number of plants with admixed D. arborea and D. ciliolata ancestry, even in nonhybrid-zone populations. These results suggest that genome-wide low levels of differentiation between D. arborea and D. ciliolata are in part due to gene flow, and favour models of genic speciation and collective evolution in which gene flow has different effects on selected loci vs. nonselected genomic regions. We discuss ecological and climatic factors that may have shaped patterns of differentiation in this species complex.

A new approach to estimate parameters of speciation models with application to apes

How populations diverge and give rise to distinct species remains a fundamental question in evolutionary biology, with important implications for a wide range of fields, from conservation genetics to human evolution. A promising approach is to estimate parameters of simple speciation models using polymorphism data from multiple loci. Existing methods, however, make a number of assumptions that severely limit their applicability, notably, no gene flow after the populations split and no intralocus recombination. To overcome these limitations, we developed a new Markov chain Monte Carlo method to estimate parameters of an isolation-migration model. The approach uses summaries of polymorphism data at multiple loci surveyed in a pair of diverging populations or closely related species and, importantly, allows for intralocus recombination. To illustrate its potential, we applied it to extensive polymorphism data from populations and species of apes, whose demographic histories are largely unknown. The isolation-migration model appears to provide a reasonable fit to the data. It suggests that the two chimpanzee species became reproductively isolated in allopatry approximately 850 Kya, while Western and Central chimpanzee populations split approximately 440 Kya but continued to exchange migrants. Similarly, Eastern and Western gorillas and Sumatran and Bornean orangutans appear to have experienced gene flow since their splits approximately 90 and over 250 Kya, respectively.

Cyto-nuclear genomic dissociation and the African elephant species question

Studies of skull morphology and of nuclear DNA have strongly concluded that African elephants comprise two species. Nonetheless, Debruyne (2005) has suggested a single-species model for Loxodonta based on the polyphyly of a single genetic locus, mitochondrial DNA (mtDNA). Discordant patterns between mitochondrial and nuclear DNA markers were subsequently reported in some African savanna elephant populations, further supporting a two-species model, and prompting us to re-examine here the geographic distribution of different elephant morphotypes and their relationship to nuclear and mtDNA phylogeographic patterns. We used exact tests to compare the distribution of forest elephant-typical and savanna elephant-typical characteristics across eight published datasets containing morphological, mtDNA or nuclear DNA data for African elephants. Among the elephants examined by Debruyne (2005), we found that patterns of forest vs. savanna characteristics were significantly different (p < 10(-5)) between mtDNA and morphology, suggesting the presence of cyto-nuclear genomic dissociation. We show that the eight African elephant continent-wide datasets compared, including that of Debruyne (2005), together support a two-species model with cyto-nuclear genomic dissociation rather than a one-species model, and together indicate that Africa harbors two species of elephant.

Estimating a geographically explicit model of population divergence

Patterns of genetic variation can provide valuable insights for deciphering the relative roles of different evolutionary processes in species differentiation. However, population-genetic models for studying divergence in geographically structured species are generally lacking. Since these are the biogeographic settings where genetic drift is expected to predominate, not only are population-genetic tests of hypotheses in geographically structured species constrained, but generalizations about the evolutionary processes that promote species divergence may also be potentially biased. Here we estimate a population-divergence model in montane grasshoppers from the sky islands of the Rocky Mountains. Because this region was directly impacted by Pleistocene glaciation, both the displacement into glacial refugia and recolonization of montane habitats may contribute to differentiation. Building on the tradition of using information from the genealogical relationships of alleles to infer the geography of divergence, here the additional consideration of the process of gene-lineage sorting is used to obtain a quantitative estimate of population relationships and historical associations (i.e., a population tree) from the gene trees of five anonymous nuclear loci and one mitochondrial locus in the broadly distributed species Melanoplus oregonensis. Three different approaches are used to estimate a model of population divergence; this comparison allows us to evaluate specific methodological assumptions that influence the estimated history of divergence. A model of population divergence was identified that significantly fits the data better compared to the other approaches, based on per-site likelihood scores of the multiple loci, and that provides clues about how divergence proceeded in M. oregonensis during the dynamic Pleistocene. Unlike the approaches that either considered only the most recent coalescence (i.e., information from a single individual per population) or did not consider the pattern of coalescence in the gene genealogies, the population-divergence model that best fits the data was estimated by considering the pattern of gene lineage coalescence across multiple individuals, as well as loci. These results indicate that sampling of multiple individuals per population is critical to obtaining an accurate estimate of the history of divergence so that the signal of common ancestry can be separated from the confounding influence of gene flow-even though estimates suggest that gene flow is not a predominant factor structuring patterns of genetic variation across these sky island populations. They also suggest that the gene genealogies contain information about population relationships, despite the lack of complete sorting of gene lineages. What emerges from the analyses is a model of population divergence that incorporates both contemporary distributions and historical associations, and shows a latitudinal and regional structuring of populations reminiscent of population displacements into multiple glacial refugia. Because the population-divergence model itself is built upon the specific events shaping the history of M. oregonensis, it provides a framework for estimating additional population-genetic parameters relevant to understanding the processes governing differentiation in geographically structured species and avoids the problems of relying on overly simplified and inaccurate divergence models. The utility of these approaches, as well as the caveats and future improvements, for estimating population relationships and historical associations relevant to genetic analyses of geographically structured species are discussed.

Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera

As a model system, Arabidopsis thaliana and its wild relatives have played an important role in the study of genomics and evolution in plants. In this study, we examined the genetic diversity of the chalcone synthase (Chs) gene, which encodes a key enzyme of the flavonoid pathway and is located on chromosome five, as well as two Chs-like genes on the first and fourth chromosomes of Arabidopsis. The objectives of the study are to determine if natural selection operates differentially on the paralogs of the Chs gene family in A. thaliana and Arabidopsis halleri ssp. gemmifera. The mode of selection was inferred from Tajima's D values from noncoding and coding regions, as well as from the ratio of nonsynonymous to synonymous substitutions. Both McDonald-Kreitman and HKA tests revealed the effects of selection on the allelic distribution, except for the chromosome 1 paralog in ssp. gemmifera. The Chs gene on chromosome 5 was under purifying selection in both species. Significant, negative Tajima's D values at synonymous sites and positive Fay and Wu's H values within coding region, plus reduced genetic variability in introns, indicated effects of background selection in shaping the evolution of this gene region in A. thaliana. The Chs paralog on chromosome 1 was under positive selection in A. thaliana, while interspecific introgression and balancing selection determined the fates of the paralog and resulted in high heterogeneity in ssp. gemmifera. Local adaptation differentiated populations of Japan and China at the locus. In contrast, the other Chs-paralog of chromosome 4 was shaped by purifying selection in A. thaliana, while under positive selection in ssp. gemmifera, as indicated by dn/ds>1. Moreover, these contrasting patterns of selection have likely resulted in functional divergence in Arabidopsis, as indicated by radical amino acid substitutions at the chalcone synthase/stilbene synthase motif of the Chs genes. Unlike previous studies of the evolutionary history of A. thaliana, the high levels of genetic diversity in most gene regions of Chs paralogs and nonsignificant Tajima's D in the intron sequences of the Chs gene family in A. thaliana did not reflect the effects of a recent demographic expansion.

DISCORDANT DIVERGENCE TIMES AMONG Z-CHROMOSOME REGIONS BETWEEN TWO ECOLOGICALLY DISTINCT SWALLOWTAIL BUTTERFLY SPECIES

We investigate multilocus patterns of differentiation between parental populations of two swallowtail butterfly species that differ at a number of ecologically important sex-linked traits. Using a new coalescent-based approach, we show that there is significant heterogeneity in estimated divergence times among five Z-linked markers, rejecting a purely allopatric speciation model. We infer that the Z chromosome is a mosaic of regions that differ in the extent of historical gene flow, potentially due to isolating barriers that prevent the introgression of species-specific traits that result in hybrid incompatibilities. Surprisingly, a candidate region for a strong barrier to introgression, Ldh, does not show a significantly deeper divergence time than other markers on the Z chromosome. Our approach can be used to test alternative models of speciation and can potentially assign chronological order to the appearance of factors contributing to reproductive isolation between species.

Widespread Genealogical Nonmonophyly in Species of Pinus Subgenus Strobus

Phylogenetic relationships among Pinus species from subgenus Strobus remain unresolved despite combined efforts based on nrITS and cpDNA. To provide greater resolution among these taxa, a 900-bp intron from a late embryogenesis abundant (LEA)-like gene (IFG8612)was sequenced from 39 pine species, with two or more alleles representing 33 species. Nineteen of 33 species exhibited allelic nonmonphyly in the strict consensus tree, and 10 deviated significantly from allelic monophyly based on topology incongruence tests. Intraspecific nucleotide diversity ranged from 0.0 to 0.0211, and analysis of variance shows that nucleotide diversity was strongly associated (P < 0.0001)with the degree of species monophyly. Although species nonmonophyly complicates phylogenetic interpretations, this nuclear locus offers greater topological support than previously observed for cpDNA or nrITS. Lacking evidence for hybridization, recombination, or imperfect taxonomy, we feel that incomplete lineage sorting remains the best explanation for the polymorphisms shared among species. Depending on the species, coalescent expectations indicate that reciprocal monophyly will be more likely than paraphyly in 1.71 to 24.0 x 10(6) years, and that complete genome-wide coalescence in these species may require up to 76.3 x 10(6) years. The absence of allelic coalescence is a severe constraint in the application of phylogenetic methods in Pinus, and taxa sharing similar life history traits with Pinus are likely to show species nonmonophyly using nuclear markers.

Recurrent replacement of mtDNA and cryptic hybridization between two sibling bat species Myotis myotis and Myotis blythii

The two sibling bat species Myotis myotis and Myotis blythii occur in sympatry over wide areas of Southern and Central Europe. Morphological, ecological and previous genetic evidence supported the view that the two species constitute two well-differentiated groups, but recent phylogenetic analyses have shown that the two species share some mtDNA haplotypes when they occur in sympatry. In order to see whether some genetic exchange has occurred between the two species, we sequenced a highly variable segment of the mitochondrial control region in both species living in sympatry and in allopatry. We also analysed the nuclear diversity of 160 individuals of both species found in two mixed nursery colonies located north and south of the Alps. MtDNA analysis confirmed that European M. blythii share multiple, identical or very similar haplotypes with M. myotis. Since allopatric Asian M. blythii presents mtDNA sequences that are very divergent from those of the two species found in Europe, we postulate that the mitochondrial genome of the European M. blythii has been replaced by that of M. myotis. The analysis of nuclear diversity shows a strikingly different pattern, as both species are well differentiated within mixed nursery colonies (F(ST) = 0.18). However, a Bayesian analysis of admixture reveals that the hybrids can be frequently observed, as about 25% of sampled M. blythii show introgressed genes of M. myotis origin. In contrast, less than 4% of the M. myotis analysed were classified as non-parental genotypes, revealing an asymmetry in the pattern of hybridization between the two species. These results show that the two species can interbreed and that the hybridization is still ongoing in the areas of sympatry. The persistence of well-differentiated nuclear gene pools, in spite of an apparent replacement of mitochondrial genome in European M. blythii by that of M. myotis, is best explained by a series of introgression events having occurred repeatedly during the recent colonization of Europe by M. blythii from Asia. The sharp contrast obtained from the analysis of mitochondrial and nuclear markers further points to the need to cautiously interpret results based on a single class of genetic markers.

Molecular identification of mosquito species in the Anopheles annularis group in southern Asia

The Anopheles annularis group of subgenus Cellia Theobald (Diptera: Culicidae) includes five currently recognized species in southern Asia: An. annularis Van der Wulp, Anopheles nivipes (Theobald) and Anopheles philippinensis Ludlow, which are widespread in the region, Anopheles pallidus Theobald, which is known in Sri Lanka, India and Myanmar, and Anopheles schueffneri Stanton, which occurs in Java and Sumatra. Identification of the four mainland species based on morphology is problematic. In view of the fact that the three widespread species are variously involved in malaria transmission in different parts of the region, we developed a species-specific polymerase chain reaction assay based on rDNA internal transcribed spacer 2 (ITS2) sequences to facilitate entomological and epidemiological studies of the four species. The method proved to be reliable when tested over a wide geographical area.

Rampant gene exchange across a strong reproductive barrier between the annual sunflowers, Helianthus annuus and H. petiolaris

Plant species may remain morphologically distinct despite gene exchange with congeners, yet little is known about the genomewide pattern of introgression among species. Here we analyze the effects of persistent gene flow on genomic differentiation between the sympatric sunflower species Helianthus annuus and H. petiolaris. While the species are strongly isolated in testcrosses, genetic distances at 108 microsatellite loci and 14 sequenced genes are highly variable and much lower (on average) than for more closely related but historically allopatric congeners. Our analyses failed to detect a positive association between levels of genetic differentiation and chromosomal rearrangements (as reported in a prior publication) or proximity to QTL for morphological differences or hybrid sterility. However, a significant increase in differentiation was observed for markers within 5 cM of chromosomal breakpoints. Together, these results suggest that islands of differentiation between these two species are small, except in areas of low recombination. Furthermore, only microsatellites associated with ESTs were identified as outlier loci in tests for selection, which might indicate that the ESTs themselves are the targets of selection rather than linked genes (or that coding regions are not randomly distributed). In general, these results indicate that even strong and genetically complex reproductive barriers cannot prevent widespread introgression.

Taxonomic considerations in listing subspecies under the U.S. Endangered Species Act

The U.S. Endangered Species Act (ESA) allows listing of subspecies and other groupings below the rank of species. This provides the U.S. Fish and Wildlife Service and the National Marine Fisheries Service with a means to target the most critical unit in need of conservation. Although roughly one-quarter of listed taxa are subspecies, these management agencies are hindered by uncertainties about taxonomic standards during listing or delisting activities. In a review of taxonomic publications and societies, we found few subspecies lists and none that stated standardized criteria for determining subspecific taxa. Lack of criteria is attributed to a centuries-old debate over species and subspecies concepts. Nevertheless, the critical need to resolve this debate for ESA listings led us to propose that minimal biological criteria to define disjunct subspecies (legally or taxonomically) should include the discreteness and significance criteria of distinct population segments (as defined under the ESA). Our subspecies criteria are in stark contrast to that proposed by supporters of the phylogenetic species concept and provide a clear distinction between species and subspecies. Efforts to eliminate or reduce ambiguity associated with subspecies-level classifications will assist with ESA listing decisions. Thus, we urge professional taxonomic societies to publish and periodically update peer-reviewed species and subspecies lists. This effort must be paralleled throughout the world for efficient taxonomic conservation to take place.

Reduced recombination rate and genetic differentiation between the M and S forms of Anopheles gambiae s.s

Genetic differentiation between the largely sympatric molecular forms M and S of Anopheles gambiae appears mostly limited to division 6 and part of division 5 of the X chromosome. This region is adjacent to the centromere and includes the rDNA that was used to define these forms. This localized differentiation between populations that experience gene flow strongly suggests that this region contains genes responsible for reproductive isolation. Regions adjacent to centromeres are known to experience less recombination in several species and it has recently been suggested that low recombination rates can facilitate the accumulation and maintenance of isolation genes in partially isolated populations. Therefore, we measured the recombination rate in division 5D/6 directly and estimate that it is at least 16-fold reduced across this region compared to the remainder of the X chromosome. Additionally, sequence data from four loci from field-collected mosquitoes from several West African countries show very strong differentiation between the molecular forms in division 5D/6, whereas none was observed in two loci elsewhere on the X chromosome. Furthermore, genetic variation was substantially lower in division 5D/6 compared to the two reference loci, and the inferred genealogies of the division 5D/6 genes show patterns consistent with selective sweeps. This suggests that the reduced recombination rate has increased the effect of selection on this region and that our data are consistent with the hypothesis that reduced recombination rates can play a role in the accumulation of isolation genes in the face of gene flow.

Drosophila flies in "Evolution Canyon" as a model for incipient sympatric speciation

The genetic basis of population divergence leading to adaptive radiation and speciation is a major unresolved problem of evolutionary biology. Molecular elucidation of "speciation genes" advanced recently, yet it remains without clear identification of the gene complexes participating in reproductive isolation between natural populations, particularly, in sympatry. Genetic divergence was discovered between Drosophila melanogaster populations inhabiting ecologically contrasting, opposite slopes in "Evolution Canyon" (EC), Mt. Carmel, Israel. Interslope migration of flies is easy and verified. Nevertheless, significant interslope D. melanogaster population divergence was established at EC involving habitat choice, mate choice, thermal and drought tolerances, adaptive genes, and mobile elements. Parallel patterns of stress tolerance, habitat choice, and mate choice were demonstrated in Drosophila simulans at EC, although on a smaller scale. However, some tests for interslope genetic differentiation in Drosophila, derived from the opposite EC slopes, gave somewhat controversial results. Here we present new empirical data on interslope genetic divergence of Drosophila at EC, and summarize previous supporting and controversial results. We suggest that Drosophila populations at EC represent a rare example, demonstrating how selection overrides migration, and propose an ad hoc ecological model of incipient sympatric divergence.

Shifting distributions and speciation: species divergence during rapid climate change

Questions about how shifting distributions contribute to species diversification remain virtually without answer, even though rapid climate change during the Pleistocene clearly impacted genetic variation within many species. One factor that has prevented this question from being adequately addressed is the lack of precision associated with estimates of species divergence made from a single genetic locus and without incorporating processes that are biologically important as populations diverge. Analysis of DNA sequences from multiple variable loci in a coalescent framework that (i) corrects for gene divergence pre-dating speciation, and (ii) derives divergence-time estimates without making a priori assumptions about the processes underlying patterns of incomplete lineage sorting between species (i.e. allows for the possibility of gene flow during speciation), is critical to overcoming the inherent logistical and analytical difficulties of inferring the timing and mode of speciation during the dynamic Pleistocene. Estimates of species divergence that ignore these processes, use single locus data, or do both can dramatically overestimate species divergence. For example, using a coalescent approach with data from six loci, the divergence between two species of montane Melanoplus grasshoppers is estimated at between 200,000 and 300,000 years before present, far more recently than divergence estimates made using single-locus data or without the incorporation of population-level processes. Melanoplus grasshoppers radiated in the sky islands of the Rocky Mountains, and the analysis of divergence between these species suggests that the isolation of populations in multiple glacial refugia was an important factor in promoting speciation. Furthermore, the low estimates of gene flow between the species indicate that reproductive isolation must have evolved rapidly for the incipient species boundaries to be maintained through the subsequent glacial periods and shifts in species distributions.

Adaptive radiation and hybridization in Wallace's Dreamponds: evidence from sailfin silversides in the Malili Lakes of Sulawesi

Adaptive radiations are extremely useful to understand factors driving speciation. A challenge in speciation research is to distinguish forces creating novelties and those relevant to divergence and adaptation. Recently, hybridization has regained major interest as a potential force leading to functional novelty and to the genesis of new species. Here, we show that introgressive hybridization is a prominent phenomenon in the radiation of sailfin silversides (Teleostei: Atheriniformes: Telmatherinidae) inhabiting the ancient Malili Lakes of Sulawesi, correlating conspicuously with patterns of increased diversity. We found the most diverse lacustrine species-group of the radiation to be heavily introgressed by genotypes originating from streams of the lake system, an effect that has masked the primary phylogenetic pattern of the flock. We conclude that hybridization could have acted as a key factor in the generation of the flock's spectacular diversity. To our knowledge, this is the first empirical evidence for massive reticulate evolution within a complex animal radiation.