The genetic basis of reproductive isolation: insights from Drosophila

Understanding reproductive isolation based on the rice model

Reproductive isolation is both an indicator of speciation and a mechanism for maintaining species identity. Here we review the progress in studies of hybrid sterility in rice to illustrate the present understanding of the molecular and evolutionary mechanisms underlying reproductive isolation. Findings from molecular characterization of genes controlling hybrid sterility can be summarized with three evolutionary genetic models. The parallel divergence model features duplicated loci generated by genome evolution; in this model, the gametes abort when the two copies of loss-of-function mutants meet in hybrids. In the sequential divergence model, mutations of two linked loci occur sequentially in one lineage, and negative interaction between the ancestral and nascent alleles of different genes causes incompatibility. The parallel-sequential divergence model involves three tightly linked loci, exemplified by a killer-protector system formed of mutations in two steps. We discuss the significance of such findings and their implications for crop improvement.

Sex and speciation: Drosophila reproductive tract proteins- twenty five years later

The protein electrophoresis revolution, nearly fifty years ago, provided the first glimpse into the nature of molecular genetic variation within and between species and showed that the amount of genetic differences between newly arisen species was minimal. Twenty years later, 2D electrophoresis showed that, in contrast to general gene-enzyme variation, reproductive tract proteins were less polymorphic within species but highly diverged between species. The 2D results were interesting and revolutionary, but somewhat uninterpretable because, at the time, rapid evolution and selective sweeps were not yet part of the common vocabulary of evolutionary biologists. Since then, genomic studies of sex and reproduction-related (SRR) genes have grown rapidly into a large area of research in evolutionary biology and are shedding light on a number of phenomena. Here we review some of the major and current fields of research that have greatly contributed to our understanding of the evolutionary dynamics and importance of SRR genes and genetic systems in understanding reproductive biology and speciation.

Evidence for compensatory evolution of ribosomal proteins in response to rapid divergence of mitochondrial rRNA

Rapid evolution of mitochondrial DNA (mtDNA) places intrinsic selective pressures on many nuclear genes involved in mitochondrial functions. Mitochondrial ribosomes, for example, are composed of mtDNA-encoded ribosomal RNAs (rRNAs) and a set of more than 60 nuclear-encoded ribosomal proteins (mRP) distinct from the cytosolic RPs (cRP). We hypothesized that the rapid divergence of mt-rRNA would result in rapid evolution of mRPs relative to cRPs, which respond to slowly evolving nuclear-encoded rRNA. In comparisons of rates of nonsynonymous and synonymous substitutions between a pair of divergent populations of the copepod Tigriopus californicus, we found that mRPs showed elevated levels of amino acid changes relative to cRPs. This pattern was equally strong at the interspecific level, between three pairs of sister species (Nasonia vitripennis vs. N. longicornis, Drosophila melanogaster vs. D. simulans, and Saccharomyces cerevisae vs. S. paradoxus). This high rate of mRP evolution may result in intergenomic incompatibilities between taxonomic lineages, and such incompatibilities could lead to dysfunction of mitochondrial ribosomes and the loss of fitness observed among interpopulation hybrids in T. californicus and interspecific hybrids in other species.

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

Background

Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated (‘allopatric’) populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations.

Results

As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation.

Conclusions

Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of ‘genomic islands of speciation’.

Genomic divergence during speciation: causes and consequences

Speciation is a fundamental process responsible for the diversity of life. Progress has been made in detecting individual 'speciation genes' that cause reproductive isolation. In contrast, until recently, less attention has been given to genome-wide patterns of divergence during speciation. Thus, major questions remain concerning how individual speciation genes are arrayed within the genome, and how this affects speciation. This theme issue is dedicated to exploring this genomic perspective of speciation. Given recent sequencing and computational advances that now allow genomic analyses in most organisms, the goal is to help move the field towards a more integrative approach. This issue draws upon empirical studies in plants and animals, and theoretical work, to review and further document patterns of genomic divergence. In turn, these studies begin to disentangle the role that different processes, such as natural selection, gene flow and recombination rate, play in generating observed patterns. These factors are considered in the context of how genomes diverge as speciation unfolds, from beginning to end. The collective results point to how experimental work is now required, in conjunction with theory and sequencing studies, to move the field from descriptive studies of patterns of divergence towards a predictive framework that tackles the causes and consequences of genome-wide patterns.

Establishment of new mutations under divergence and genome hitchhiking

Theoretical models addressing genome-wide patterns of divergence during speciation are needed to help us understand the evolutionary processes generating empirical patterns. Here, we examine a critical issue concerning speciation-with-gene flow: to what degree does physical linkage (r < 0.5) of new mutations to already diverged genes aid the build-up of genomic islands of differentiation? We used simulation and analytical approaches to partition the probability of establishment for a new divergently selected mutation when the mutation (i) is the first to arise in an undifferentiated genome (the direct effect of selection), (ii) arises unlinked to any selected loci (r = 0.5), but within a genome that has some already diverged genes (the effect of genome-wide reductions in gene flow for facilitating divergence, which we term 'genome hitchhiking'), and (iii) arises in physical linkage to a diverged locus (divergence hitchhiking). We find that the strength of selection acting directly on a new mutation is generally the most important predictor for establishment, with divergence and genomic hitchhiking having smaller effects. We outline the specific conditions under which divergence and genome hitchhiking can aid mutation establishment. The results generate predictions about genome divergence at different points in the speciation process and avenues for further work.

The genomics of speciation-with-gene-flow

The emerging field of speciation genomics is advancing our understanding of the evolution of reproductive isolation from the individual gene to a whole-genome perspective. In this new view it is important to understand the conditions under which 'divergence hitchhiking' associated with the physical linkage of gene regions, versus 'genome hitchhiking' associated with reductions in genome-wide rates of gene flow caused by selection, can enhance speciation-with-gene-flow. We describe here a theory predicting four phases of speciation, defined by changes in the relative effectiveness of divergence and genome hitchhiking, and review empirical data in light of the theory. We outline future directions, emphasizing the need to couple next-generation sequencing with selection, transplant, functional genomics, and mapping studies. This will permit a natural history of speciation genomics that will help to elucidate the factors responsible for population divergence and the roles that genome structure and different forms of hitchhiking play in facilitating the genesis of new biodiversity.

Is speciation accompanied by rapid evolution? Insights from comparing reproductive and nonreproductive transcriptomes in Drosophila

The tempo and mode of evolutionary change during speciation have remained contentious until recently. While much of the evidence claiming speciation is an abrupt and rapid process comes from fossil data, recent molecular phylogenetics show that the background of gradual evolution is often broken by accelerated rates of molecular evolution during speciation. However, what kinds of genes affect or are affected by speciation remains unexplored. Our analysis of 4843 protein-coding genes in five species of the Drosophila melanogaster subgroup shows that while ~70% of genes follow clock-like evolution, between 17-19.67% of loci show signatures of accelerated rates of evolution in recently formed species. These genes show 2-3-fold higher rates of substitution in recently diverged species compared to older species. This fraction of loci affects a diverse range of functions. Only a small proportion of reproductive genes experience speciation-related accelerated changes but many sex-and -reproduction related genes show an interesting pattern of persistent rapid evolution suggesting that sex-and-reproduction related genes are under constant selective pressures. The identification of loci associated with accelerated evolution allows us to address the mechanisms of rapid evolution and speciation, which in our study appears to be a combination of both selection and rapid demographical changes.

Signatures of selection and sex-specific expression variation of a novel duplicate during the evolution of the Drosophila desaturase gene family

The tempo and mode of evolution of loci with a large effect on adaptation and reproductive isolation will influence the rate of evolutionary divergence and speciation. Desaturase loci are involved in key biochemical changes in long-chain fatty acids. In insects, these have been shown to influence adaptation to starvation or desiccation resistance and in some cases act as important pheromones. The desaturase gene family of Drosophila is known to have evolved by gene duplication and diversification, and at least one locus shows rapid evolution of sex-specific expression variation. Here, we examine the evolution of the gene family in species representing the Drosophila phylogeny. We find that the family includes more loci than have been previously described. Most are represented as single-copy loci, but we also find additional examples of duplications in loci which influence pheromone blends. Most loci show patterns of variation associated with purifying selection, but there are strong signatures of diversifying selection in new duplicates. In the case of a new duplicate of desat1 in the obscura group species, we show that strong selection on the coding sequence is associated with the evolution of sex-specific expression variation. It seems likely that both sexual selection and ecological adaptation have influenced the evolution of this gene family in Drosophila.

Genes versus phenotypes in the study of speciation

Despite persistent debate on the nature of species, the widespread adoption of Mayr's biological species concept has led to a heavy emphasis on the importance of reproductive isolation to the speciation process. Equating the origin of species with the evolution of reproductive isolation has become common practice in the study of speciation, coincident with an increasing focus on elucidating the specific genetic changes (i.e.-speciation genes) underlying intrinsic reproductive barriers between species. In contrast, some have recognized that reproductive isolation is usually a signature effect rather than a primary cause of speciation. Here we describe a research paradigm that shifts emphasis from effects to causes in order to resolve this apparent contradiction and galvanize the study of speciation. We identify major elements necessary for a balanced and comprehensive investigation of the origin of species and place the study of so-called "speciation genes" into its appropriate context. We emphasize the importance of characterizing diverging phenotypes, identifying relevant evolutionary forces acting on those phenotypes and their role in the causal origins of reduced gene flow between incipient species, and the nature of the genetic and phenotypic boundaries that results from such processes. This approach has the potential to unify the field of speciation research, by allowing us to make better "historical" predictions about the fate of diverging populations regardless of taxon.

Scanning the European corn borer (Ostrinia spp.) genome for adaptive divergence between host-affiliated sibling species

It has recently been shown that the European corn borer, a major pest of maize crops, is actually composed of two genetically differentiated and reproductively isolated taxa, which are found in sympatry over a wide geographical range in Eurasia. Each taxon is adapted to specific host plants: Ostrinia nubilalis feeds mainly on maize, while O. scapulalis feeds mainly on hop or mugwort. Here, we present a genome scan approach as a first step towards an integrated molecular analysis of the adaptive genomic divergence between O. nubilalis and O. scapulalis. We analysed 609 AFLP marker loci in replicate samples of sympatric populations of Ostrinia spp. collected on maize, hop and mugwort, in France. Using two genome scan methods based on the analysis of population differentiation, we found a set of 28 outlier loci that departed from the neutral expectation in one or the other method (of which a subset of 14 loci were common to both methods), which showed a significantly increased differentiation between O. nubilalis and O. scapulalis, when compared to the rest of the genome. A subset of 12 outlier loci were sequenced, of which 7 were successfully re-amplified as target candidate loci. Three of these showed homology with annotated lepidopteran sequences from public nucleotide databases.

Inference of selection and stochastic effects in the house mouse hybrid zone

We explored the transition of 13 X-linked markers across two separate portions of the house mouse hybrid zone, asking whether such a comparison can distinguish the effects of selection from random factors. A heuristic search in the likelihood landscape revealed more complex likelihood profiles for data sampled in two-dimensional (2D) space relative to data sampled along a linear transect. Randomized resampling of localities analyzed for individual loci showed that deletion of sites away from the zone center can decrease cline width estimates whereas deletion of sites close to the center can significantly increase the width estimates. Deleting localities for all loci resulted in wider clines if the number of samples from the center was limited. The results suggest that, given the great variation in width estimates resulting from inclusion/exclusion of sampling sites, the geographic sampling design is important in hybrid zone studies and that our inferences should take into account measures of uncertainty such as support intervals. The comparison of the two transects indicates cline widths are narrower for loci in the central part of the X chromosome, suggesting selection is stronger in this region and genetic incompatibilities may have at least partly common architecture in the house mouse hybrid zone.

Evolutionary patterns of two major reproduction candidate genes (Zp2 and Zp3) reveal no contribution to reproductive isolation between bovine species

BACKGROUND

It has been established that mammalian egg zona pellucida (ZP) glycoproteins are responsible for species-restricted binding of sperm to unfertilized eggs, inducing the sperm acrosome reaction, and preventing polyspermy. In mammals, ZP apparently represents a barrier to heterospecific fertilization and thus probably contributes to reproductive isolation between species. The evolutionary relationships between some members of the tribe Bovini are complex and highly debatable, particularly, those involving Bos and Bison species for which interspecific hybridization is extensively documented. Because reproductive isolation is known to be a major precursor of species divergence, testing evolutionary patterns of ZP glycoproteins may shed some light into the speciation process of these species. To this end, we have examined intraspecific and interspecific genetic variation of two ZP genes (Zp2 and Zp3) for seven representative species (111 individuals) from the Bovini tribe, including five species from Bos and Bison, and two species each from genera Bubalus and Syncerus.

RESULTS

A pattern of low levels of intraspecific polymorphism and interspecific divergence was detected for the two sequenced fragments each for Zp2 and Zp3. At intraspecific level, none of neutrality tests detected deviations from neutral equilibrium expectations for the two genes. Several haplotypes in both genes were shared by multiple species from Bos and Bison.

CONCLUSIONS

Here we argue that neither ancestral polymorphism nor introgressive hybridization alone can fully account for haplotype sharing among species from Bos and Bison, and that both scenarios have contributed to such a pattern of haplotype sharing observed here. Additionally, codon-based tests revealed strong evidence for purifying selection in the Zp3 coding haplotype sequences and weak evidence for purifying selection in the Zp2 coding haplotype sequences. Contrary to a general genetic pattern that genes or genomic regions contributing to reproductive isolation between species often evolve rapidly and show little or no gene flow between species, these results demonstrate that, particularly, those sequenced exons of the Zp2 and the Zp3 did not show any contribution to reproductive isolation between the bovine species studied here.

Ecological Divergence and the Origins of Intrinsic Postmating Isolation with Gene Flow

The evolution of intrinsic postmating isolation has received much attention, both historically and in recent studies of speciation genes. Intrinsic isolation often stems from between-locus genetic incompatibilities, where alleles that function well within species are incompatible with one another when brought together in the genome of a hybrid. It can be difficult for such incompatibilities to originate when populations diverge with gene flow, because deleterious genotypic combinations will be created and then purged by selection. However, it has been argued that if genes underlying incompatibilities are themselves subject to divergent selection, then they might overcome gene flow to diverge between populations, resulting in the origin of incompatibilities. Nonetheless, there has been little explicit mathematical exploration of such scenarios for the origin of intrinsic incompatibilities during ecological speciation with gene flow. Here we explore theoretical models for the origin of intrinsic isolation where genes subject to divergent natural selection also affect intrinsic isolation, either directly or via linkage disequilibrium with other loci. Such genes indeed overcome gene flow, diverge between populations, and thus result in the evolution of intrinsic isolation. We also examine barriers to neutral gene flow. Surprisingly, we find that intrinsic isolation sometimes weakens this barrier, by impeding differentiation via ecologically based divergent selection.

Insights into species divergence and the evolution of hermaphroditism from fertile interspecies hybrids of Caenorhabditis nematodes

The architecture of both phenotypic variation and reproductive isolation are important problems in evolutionary genetics. The nematode genus Caenorhabditis includes both gonochoristic (male/female) and androdioecious (male/hermaprodite) species. However, the natural genetic variants distinguishing reproductive mode remain unknown, and nothing is known about the genetic basis of postzygotic isolation in the genus. Here we describe the hybrid genetics of the first Caenorhabditis species pair capable of producing fertile hybrid progeny, the gonochoristic Caenorhabditis sp. 9 and the androdioecious C. briggsae. Though many interspecies F(1) arrest during embryogenesis, a viable subset develops into fertile females and sterile males. Reciprocal parental crosses reveal asymmetry in male-specific viability, female fertility, and backcross viability. Selfing and spermatogenesis are extremely rare in XX F(1), and almost all hybrid self-progeny are inviable. Consistent with this, F(1) females do not express male-specific molecular germline markers. We also investigated three approaches to producing hybrid hermaphrodites. A dominant mutagenesis screen for self-fertile F(1) hybrids was unsuccessful. Polyploid F(1) hybrids with increased C. briggsae genomic material did show elevated rates of selfing, but selfed progeny were mostly inviable. Finally, the use of backcrosses to render the hybrid genome partial homozygous for C. briggsae alleles did not increase the incidence of selfing or spermatogenesis relative to the F(1) generation. These hybrid animals were genotyped at 23 loci, and significant segregation distortion (biased against C. briggsae) was detected at 13 loci. This, combined with an absence of productive hybrid selfing, prevents formulation of simple hypotheses about the genetic architecture of hermaphroditism. In the near future, this hybrid system will likely be fruitful for understanding the genetics of reproductive isolation in Caenorhabditis.

Reproductive isolation in hybrid mice due to spermatogenesis defects at three meiotic stages

Early in the process of speciation, reproductive failures occur in hybrid animals between genetically diverged populations. The sterile hybrid animals are often males in mammals and they exhibit spermatogenic disruptions, resulting in decreased number and/or malformation of mature sperms. Despite the generality of this phenomenon, comparative study of phenotypes in hybrid males from various crosses has not been done, and therefore the comprehensive genetic basis of the disruption is still elusive. In this study, we characterized the spermatogenic phenotype especially during meiosis in four different cases of reproductive isolation: B6-ChrX(MSM), PGN-ChrX(MSM), (B6 × Mus musculus musculus-NJL/Ms) F(1), and (B6 × Mus spretus) F(1). The first two are consomic strains, both bearing the X chromosome of M. m. molossinus; in B6-ChrX(MSM), the genetic background is the laboratory strain C57BL/6J (predominantly M. m. domesticus), while in PGN-ChrX(MSM) the background is the PGN2/Ms strain purely derived from wild M. m. domesticus. The last two cases are F(1) hybrids between mouse subspecies or species. Each of the hybrid males exhibited cell-cycle arrest and/or apoptosis at either one or two of three distinct meiotic stages: premeiotic stage, zygotene-to-pachytene stage of prophase I, and metaphase I. This study shows that the sterility in hybrid males is caused by spermatogenic disruptions at multiple stages, suggesting that the responsible genes function in different cellular processes. Furthermore, the stages with disruptions are not correlated with the genetic distance between the respective parental strains.

Speciation genes in plants

BACKGROUND

Analyses of speciation genes--genes that contribute to the cessation of gene flow between populations--can offer clues regarding the ecological settings, evolutionary forces and molecular mechanisms that drive the divergence of populations and species. This review discusses the identities and attributes of genes that contribute to reproductive isolation (RI) in plants, compares them with animal speciation genes and investigates what these genes can tell us about speciation.

SCOPE

Forty-one candidate speciation genes were identified in the plant literature. Of these, seven contributed to pre-pollination RI, one to post-pollination, prezygotic RI, eight to hybrid inviability, and 25 to hybrid sterility. Genes, gene families and genetic pathways that were frequently found to underlie the evolution of RI in different plant groups include the anthocyanin pathway and its regulators (pollinator isolation), S RNase-SI genes (unilateral incompatibility), disease resistance genes (hybrid necrosis), chimeric mitochondrial genes (cytoplasmic male sterility), and pentatricopeptide repeat family genes (cytoplasmic male sterility).

CONCLUSIONS

The most surprising conclusion from this review is that identities of genes underlying both prezygotic and postzygotic RI are often predictable in a broad sense from the phenotype of the reproductive barrier. Regulatory changes (both cis and trans) dominate the evolution of pre-pollination RI in plants, whereas a mix of regulatory mutations and changes in protein-coding genes underlie intrinsic postzygotic barriers. Also, loss-of-function mutations and copy number variation frequently contribute to RI. Although direct evidence of positive selection on speciation genes is surprisingly scarce in plants, analyses of gene family evolution, along with theoretical considerations, imply an important role for diversifying selection and genetic conflict in the evolution of RI. Unlike in animals, however, most candidate speciation genes in plants exhibit intraspecific polymorphism, consistent with an important role for stochastic forces and/or balancing selection in development of RI in plants.

Transmission ratio distortion results in asymmetric introgression in Louisiana Iris

BACKGROUND

Linkage maps are useful tools for examining both the genetic architecture of quantitative traits and the evolution of reproductive incompatibilities. We describe the generation of two genetic maps using reciprocal interspecific backcross 1 (BC1) mapping populations from crosses between Iris brevicaulis and Iris fulva. These maps were constructed using expressed sequence tag (EST)- derived codominant microsatellite markers. Such a codominant marker system allowed for the ability to link the two reciprocal maps, and compare patterns of transmission ratio distortion observed between the two.

RESULTS

Linkage mapping resulted in markers that coalesced into 21 linkage groups for each of the reciprocal backcross maps, presumably corresponding to the 21 haploid chromosomes of I. brevicaulis and I. fulva. The composite map was 1190.0-cM long, spanned 81% of the I. brevicaulis and I. fulva genomes, and had a mean density of 4.5 cM per locus. Transmission ratio distortion (TRD) was observed in 138 (48.5%) loci distributed in 19 of the 21 LGs in BCIB, BCIF, or both BC1 mapping populations. Of the distorted markers identified, I. fulva alleles were detected at consistently higher-than-expected frequencies in both mapping populations.

CONCLUSIONS

The observation that I. fulva alleles are overrepresented in both mapping populations suggests that I. fulva alleles are favored to introgress into I. brevicaulis genetic backgrounds, while I. brevicaulis alleles would tend to be prevented from introgressing into I. fulva. These data are consistent with the previously observed patterns of introgression in natural hybrid zones, where I. fulva alleles have been consistently shown to introgress across species boundaries.

The nuclear component of a cytonuclear hybrid incompatibility in Mimulus maps to a cluster of pentatricopeptide repeat genes

Characterizing the genetic and molecular basis of hybrid incompatibilities is a first step toward understanding their evolutionary origins. We fine mapped the nuclear restorer (Rf) of cytoplasm-dependent anther sterility in Mimulus hybrids by identifying and targeting regions of the Mimulus guttatus genome containing large numbers of candidate pentatricopeptide repeat genes (PPRs). The single Mendelian locus Rf was first isolated to a 1.3-cM region on linkage group 7 that spans the genome's largest cluster of PPRs, then split into two tightly linked loci (Rf1 and Rf2) by <10 recombination events in a large (N = 6153) fine-mapping population. Progeny testing of fertile recombinants demonstrated that a dominant M. guttatus allele at each Rf locus was sufficient to restore fertility. Each Rf locus spans a physical region containing numerous PPRs with high homology to each other, suggesting recent tandem duplication or transposition. Furthermore, these PPRs have higher homology to restorers in distantly related taxa (petunia and rice) than to PPRs elsewhere in the Mimulus genome. These results suggest that the cytoplasmic male sterility (CMS)-PPR interaction is highly conserved across flowering plants. In addition, given our theoretical understanding of cytonuclear coevolution, the finding that hybrid CMS results from interactions between a chimeric mitochondrial transcript that is modified by Rf loci identified as PPRs is consistent with a history of selfish mitochondrial evolution and compensatory nuclear coevolution within M. guttatus.

The impact of the triploid block on the origin and evolution of polyploid plants

Polyploidization, a widespread phenomenon among plants, is considered a major speciation mechanism. Polyploid plants have a high degree of immediate post-zygotic reproductive isolation from their progenitors, as backcrossing to either parent will produce mainly nonviable progeny. This reproductive barrier is called triploid block and it is caused by malfunction of the endosperm. Nevertheless, the main route to polyploid formation is via unreduced gametes and unstable triploid progeny, suggesting that there are ways to overcome the triploid block. Until recently, the mechanistic basis for unreduced gamete formation and the triploid block were completely unknown. Recent developments have revealed genetic pathways leading to unreduced gamete formation as well as the underlying genetic basis for the triploid block in Arabidopsis. These novel findings will provide the basis for a genetic understanding of polyploid formation and subsequent speciation in plants.

Interaction of novel Dobzhansky-Muller type genes for the induction of hybrid lethality between Gossypium hirsutum and G. barbadense cv. Coastland R4-4

Hybrid lethality was identified in interspecific hybrids between two cotton species, Gossypium hirsutum and G. barbadense cv. Coastland R4-4 (R4-4). Genetic analysis indicated that the lethal symptom was controlled by two dominant complementary genes, one from G. hirsutum and another from R4-4. Microsatellite mapping identified the location of the causal gene in G. hirsutum as chromosome D8, while the R4-4 gene was placed on chromosome D11. Our data indicate that these genes conform to the Dobzhansky-Muller model, and are novel for the induction of hybrid lethality in Gossypium. Following the genetic nomenclature, we propose that the two novel Dobzhansky-Muller genes from G. hirsutum and from R4-4 be named Le ( 3 ) and Le ( 4 ), respectively. Given what we know about their inheritance patterns, their genotypes should be Le ( 3 ) Le ( 3 ) le ( 4 ) le ( 4 ) in G. hirsutum, and le ( 3 ) le ( 3 ) Le ( 4 ) Le ( 4 ) in R4-4. Data from this study supported previous information in that expression of the lethal symptom might be affected by the dosage of causal alleles and the environment in which plants are growing.

Identifying footprints of directional and balancing selection in marine and freshwater three-spined stickleback (Gasterosteus aculeatus) populations

Natural selection is expected to leave an imprint on the neutral polymorphisms at the adjacent genomic regions of a selected gene. While directional selection tends to reduce within-population genetic diversity and increase among-population differentiation, the reverse is expected under balancing selection. To identify targets of natural selection in the three-spined stickleback (Gasterosteus aculeatus) genome, 103 microsatellite and two indel markers including expressed sequence tags (EST) and quantitative trait loci (QTL)-associated loci, were genotyped in four freshwater and three marine populations. The results indicated that a high proportion of loci (14.7%) might be affected by balancing selection and a lower proportion (2.8%) by directional selection. The strongest signatures of directional selection were detected in a microsatellite locus and two indel markers located in the intronic regions of the Eda-gene coding for the number of lateral plates. Yet, other microsatellite loci previously found to be informative in QTL-mapping studies revealed no signatures of selection. Two novel microsatellite loci (Stn12 and Stn90) located in chromosomes I and VIII, respectively, showed signals of directional selection and might be linked to genomic regions containing gene(s) important for adaptive divergence. Although the coverage of the total genomic content was relatively low, the predominance of balancing selection signals is in agreement with the contention that balancing, rather than directional selection is the predominant mode of selection in the wild.

Learning about modes of speciation by computational approaches

How often do the early stages of speciation occur in the presence of gene flow? To address this enduring question, a number of recent papers have used computational approaches, estimating parameters of simple divergence models from multilocus polymorphism data collected in closely related species. Applications to a variety of species have yielded extensive evidence for migration, with the results interpreted as supporting the widespread occurrence of parapatric speciation. Here, we conduct a simulation study to assess the reliability of such inferences, using a program that we recently developed MCMC estimation of the isolation-migration model allowing for recombination (MIMAR) as well as the program isolation-migration (IM) of Hey and Nielsen (2004). We find that when one of many assumptions of the isolation-migration model is violated, the methods tend to yield biased estimates of the parameters, potentially lending spurious support for allopatric or parapatric divergence. More generally, our results highlight the difficulty in drawing inferences about modes of speciation from the existing computational approaches alone.

Polymorphism in hybrid male sterility in wild-derived Mus musculus musculus strains on proximal chromosome 17

The hybrid sterility-1 (Hst1) locus at Chr 17 causes male sterility in crosses between the house mouse subspecies Mus musculus domesticus (Mmd) and M. m. musculus (Mmm). This locus has been defined by its polymorphic variants in two laboratory strains (Mmd genome) when mated to PWD/Ph mice (Mmm genome): C57BL/10 (carrying the sterile allele) and C3H (fertile allele). The occurrence of sterile and/or fertile (wild Mmm x C57BL)F1 males is evidence that polymorphism for this trait also exists in natural populations of Mmm; however, the nature of this polymorphism remains unclear. Therefore, we derived two wild-origin Mmm strains, STUS and STUF, that produce sterile and fertile males, respectively, in crosses with C57BL mice. To determine the genetic basis underlying male fertility, the (STUS x STUF)F1 females were mated to C57BL/10 J males. About one-third of resulting hybrid males (33.8%) had a significantly smaller epididymis and testes than parental animals and lacked spermatozoa due to meiotic arrest. A further one-fifth of males (20.3%) also had anomalous reproductive traits but produced some spermatozoa. The remaining fertile males (45.9%) displayed no deviation from values found in parental individuals. QTL analysis of the progeny revealed strong associations of male fitness components with the proximal end of Chr 17, and a significant effect of the central section of Chr X on testes mass. The data suggest that genetic incompatibilities associated with male sterility have evolved independently at the proximal end of Chr 17 and are polymorphic within both Mmd and Mmm genomes.

A mouse speciation gene encodes a meiotic histone H3 methyltransferase

Speciation genes restrict gene flow between the incipient species and related taxa. Three decades ago, we mapped a mammalian speciation gene, hybrid sterility 1 (Hst1), in the intersubspecific hybrids of house mouse. Here, we identify this gene as Prdm9, encoding a histone H3 lysine 4 trimethyltransferase. We rescued infertility in male hybrids with bacterial artificial chromosomes carrying Prdm9 from a strain with the "fertility" Hst1(f) allele. Sterile hybrids display down-regulated microrchidia 2B (Morc2b) and fail to compartmentalize gammaH2AX into the pachynema sex (XY) body. These defects, seen also in Prdm9-null mutants, are rescued by the Prdm9 transgene. Identification of a vertebrate hybrid sterility gene reveals a role for epigenetics in speciation and opens a window to a hybrid sterility gene network.

The phenomics and expression quantitative trait locus mapping of brain transcriptomes regulating adaptive divergence in lake whitefish species pairs (Coregonus sp.)

We used microarrays and a previously established linkage map to localize the genetic determinants of brain gene expression for a backcross family of lake whitefish species pairs (Coregonus sp.). Our goals were to elucidate the genomic distribution and sex specificity of brain expression QTL (eQTL) and to determine the extent to which genes controlling transcriptional variation may underlie adaptive divergence in the recently evolved dwarf (limnetic) and normal (benthic) whitefish. We observed a sex bias in transcriptional genetic architecture, with more eQTL observed in males, as well as divergence in genome location of eQTL between the sexes. Hotspots of nonrandom aggregations of up to 32 eQTL in one location were observed. We identified candidate genes for species pair divergence involved with energetic metabolism, protein synthesis, and neural development on the basis of colocalization of eQTL for these genes with eight previously identified adaptive phenotypic QTL and four previously identified outlier loci from a genome scan in natural populations. Eighty-eight percent of eQTL-phenotypic QTL colocalization involved growth rate and condition factor QTL, two traits central to adaptive divergence between whitefish species pairs. Hotspots colocalized with phenotypic QTL in several cases, revealing possible locations where master regulatory genes, such as a zinc-finger protein in one case, control gene expression directly related to adaptive phenotypic divergence. We observed little evidence of colocalization of brain eQTL with behavioral QTL, which provides insight into the genes identified by behavioral QTL studies. These results extend to the transcriptome level previous work illustrating that selection has shaped recent parallel divergence between dwarf and normal lake whitefish species pairs and that metabolic, more than morphological, differences appear to play a key role in this divergence.

Genetics of intrinsic postzygotic isolation in a circumpolar plant species, Draba nivalis (Brassicaceae)

Sterility barriers, ranging from incomplete to fully developed, were recently demonstrated within taxonomic species of the genus Draba, suggesting the existence of numerous, cryptic biological species. Because these taxa are predominately selfers and of Pleistocene origin, it was concluded that hybrid sterility evolved quickly and possibly by genetic drift. Here we used genetic mapping and QTL analyses to determine the genetic basis of hybrid sterility between geographically distant populations of one of these taxonomic species, Draba nivalis. Fifty microsatellite loci were mapped, and QTL analyses identified five loci underlying seed fertility and two underlying pollen fertility. Four of five seed fertility QTLs reduced fertility in heterozygotes, an observation most consistent with drift-based fixation of underdominant sterility loci. However, several nuclear-nuclear interactions were also found, including two that acted like reciprocal translocations with lowest fitness in double heterozygotes, and two that had a pattern of fitness consistent with Bateson-Dobzhansky-Muller incompatibilities. In contrast, pollen fertility QTLs exhibited additive inheritance, with lowest fertility associated with the paternal allele, a pattern of inheritance suggestive of cytonuclear incompatibilities. The results imply that multiple genetic mechanisms underlie the rapid evolution of reproductive barriers in Draba.

Cytonuclear genic incompatibilities cause increased mortality in male F2 hybrids of Nasonia giraulti and N. vitripennis

The haplodiploid wasp genus Nasonia is a promising model for studying the evolution of genic incompatibilities due to the existence of interfertile species and haploid males. The latter allows for significantly reducing the sample size required to detect and map recessive dysfunctional genic interactions. We exploited these features to study the genetics of intrinsic hybrid inviability in male F2 hybrids of Nasonia giraulti and N. vitripennis. Analyzing marker segregation in 225 hybrid embryos, we inferred a linkage map with 38 framework markers. The markers were tested for marker transmission ratio distortion (MTRD) and interchromosomal linkage disequilibrium in populations of embryonic and adult hybrids. We found evidence for four transmission ratio distorting loci (TRDL). Three TRDL showed a deficit of the N. giraulti allele in hybrids with N. vitripennis cytoplasm. A separate TRDL exhibited a deficiency of the N. vitripennis allele in hybrids with N. giraulti cytoplasm. We ascribe the observed MTRD in adult hybrids to cytonuclear genic incompatibilities causing differential mortality during development since hybrid embryos did not show MTRD. The identified cytonuclear genic incompatibilities in F2 hybrids with N. vitripennis cytoplasm account for most of the intrinsic hybrid inviability in this cross. The high mortality rate in F2 hybrids with N. giraulti cytoplasm cannot be explained by the single identified TRDL alone, however.

Hitchhiking mapping reveals a candidate genomic region for natural selection in three-spined stickleback chromosome VIII

Identification of genes and genomic regions under directional natural selection has become one of the major goals in evolutionary genetics, but relatively little work to this end has been done by applying hitchhiking mapping to wild populations. Hitchhiking mapping starts from a genome scan using a randomly spaced set of molecular markers followed by a fine-scale analysis in the flanking regions of the candidate regions under selection. We used the hitchhiking mapping approach to narrow down a selective sweep in the genomic region flanking a candidate locus (Stn90) in chromosome VIII in the three-spined stickleback (Gasterosteus aculeatus). Twenty-four microsatellite markers were screened in an approximately 800-kb region around the candidate locus in three marine and four freshwater populations. The patterns of genetic diversity and differentiation in the candidate region were compared to those of a putatively neutral set of markers. The Bayesian FST-test indicated an elevated genetic differentiation, deviating significantly from neutral expectations, at a continuous region of approximately 20 kb upstream from the candidate locus. Furthermore, a method developed for an array of microsatellite markers rejected neutrality in a region of approximately 90 kb flanking the candidate locus supporting the selective sweep hypothesis. Likewise, the genomewide pattern of genetic diversity differed from the candidate region in a bottleneck analysis suggesting that selection, rather than demography, explains the reduced genetic diversity at the candidate interval. The neutrality tests suggest that the selective sweep had occurred mainly in the Lake Pulmanki population, but the results from bottleneck analyses indicate that selection might have operated in other populations as well. These results suggest that the narrow interval around locus Stn90 has likely been under directional selection, but the region contains several predicted genes, each of which can be the actual targets of selection. Understanding of the functional significance of this genomic region in an ecological context will require a more detailed sequence analysis.

Mouse consomic strains: exploiting genetic divergence between Mus m. musculus and Mus m. domesticus subspecies

Consomic (chromosome substitution) strains (CSs) represent the most recent addition to the mouse genetic resources aimed to genetically analyze complex trait loci (QTLs). In this study, we report the development of a set of 28 mouse intersubspecific CSs. In each CS, we replaced a single chromosome of the C57BL/6J (B6) inbred strain (mostly Mus m. domesticus) with its homolog from the PWD/Ph inbred strain of the Mus m. musculus subspecies. These two progenitor subspecies diverged less than 1 million years ago and accumulated a large number of genetic differences that constitute a rich resource of genetic variation for QTL analyses. Altogether, the 18 consomic, nine subconsomic, and one conplastic strain covered all 19 autosomes, X and Y sex chromosomes, and mitochondrial DNA. Most CSs had significantly lower reproductive fitness compared with the progenitor strains. CSs homosomic for chromosomes 10 and 11, and the C57BL/6J-Chr X males, failed to reproduce and were substituted by less affected subconsomics carrying either a proximal, central, or distal part of the respective chromosome. A genome-wide scan of 965 DNA markers revealed 99.87% purity of the B6 genetic background. Thirty-three nonsynonymous substitutions were uncovered in the protein-coding regions of the mitochondrial DNA of the B6.PWD-mt conplastic strain. A pilot-phenotyping experiment project revealed a high number of variations among B6.PWD consomics.

First evidence for postzygotic reproductive isolation between two populations of Eurasian perch (Perca fluviatilis L.) within Lake Constance

BACKGROUND

The evolution of reproductive traits, such as hybrid incompatibility (postzygotic isolation) and species recognition (prezygotic isolation), have shown their key role in speciation. Theoretical modeling has recently predicted that close linkage between genes controlling pre- and postzygotic reproductive isolation could accelerate the conditions for speciation. Postzygotic isolation could develop during the sympatric speciation process contributing to the divergence of populations. Using hybrid fitness as a measure of postzygotic reproductive isolation, we empirically studied population divergence in perch (Perca fluviatilis L.) from two genetically divergent populations within a lake.

RESULTS

During spawning time of perch we artificially created parental offspring and F1 hybrids of the two populations and studied fertilization rate and hatching success under laboratory conditions. The combined fitness measure (product of fertilization rate and hatching success) of F1 hybrids was significantly reduced compared to offspring from within population crosses.

CONCLUSION

Our results suggest intrinsic genetic incompatibility between the two populations and indicate that population divergence between two populations of perch inhabiting the same lake may indeed be promoted by postzygotic isolation.

Two alternate mechanisms contribute to the persistence of interdependent lineages in Pogonomyrmex harvester ants

Some populations of Pogonomyrmex harvester ants comprise pairs of highly differentiated lineages with queens mating at random with several males of their own and of the alternate lineage. These queens produce two types of diploid offspring, those fertilized by males of the queens' lineage which develop into new queens and those fertilized by males of the other lineage which mostly develop into functionally sterile workers. This unusual mode of genetic caste determination has been found in 26 populations and a total of four lineage pairs (F(1)-F(2), G(1)-G(2), H(1)-H(2) and J(1)-J(2)) have been described in these populations. Despite the fact that a few interlineage queens are produced, previous studies revealed that there is a complete lack of genetic introgression between lineages. Here we quantify the proportion of interlineage queens produced in each of the four lineage pairs and determine the fate of these queens. In the F(1)-F(2), G(1)-G(2) and H(1)-H(2) lineage pairs, interlineage queens were produced by a minority of colonies. These colonies exclusively produced interlineage queens and workers, suggesting that interlineage eggs can develop into queens in these three pairs of lineages in the absence of competition with pure-lineage brood. An analysis of three key stages of the colony life cycle revealed that colonies headed by interlineage queens failed to grow sufficiently to produce reproductive individuals. In laboratory comparisons, interlineage queens produced fewer viable eggs, with the effect that they raised fewer workers and lost more weight per worker produced than pure-lineage queens. In the J(1)-J(2) lineage pair, we did not find a single interlineage queen, raising the possibility that interlineage eggs have completely lost the ability to develop into queens in this lineage pair. Hence, two distinct mechanisms seem to account for the complete lack of between-lineage gene flow in the F(1)-F(2), G(1)-G(2), H(1)-H(2) and J(1)-J(2) lineage pairs.

The locus of evolution: evo devo and the genetics of adaptation

An important tenet of evolutionary developmental biology ("evo devo") is that adaptive mutations affecting morphology are more likely to occur in the cis-regulatory regions than in the protein-coding regions of genes. This argument rests on two claims: (1) the modular nature of cis-regulatory elements largely frees them from deleterious pleiotropic effects, and (2) a growing body of empirical evidence appears to support the predominant role of gene regulatory change in adaptation, especially morphological adaptation. Here we discuss and critique these assertions. We first show that there is no theoretical or empirical basis for the evo devo contention that adaptations involving morphology evolve by genetic mechanisms different from those involving physiology and other traits. In addition, some forms of protein evolution can avoid the negative consequences of pleiotropy, most notably via gene duplication. In light of evo devo claims, we then examine the substantial data on the genetic basis of adaptation from both genome-wide surveys and single-locus studies. Genomic studies lend little support to the cis-regulatory theory: many of these have detected adaptation in protein-coding regions, including transcription factors, whereas few have examined regulatory regions. Turning to single-locus studies, we note that the most widely cited examples of adaptive cis-regulatory mutations focus on trait loss rather than gain, and none have yet pinpointed an evolved regulatory site. In contrast, there are many studies that have both identified structural mutations and functionally verified their contribution to adaptation and speciation. Neither the theoretical arguments nor the data from nature, then, support the claim for a predominance of cis-regulatory mutations in evolution. Although this claim may be true, it is at best premature. Adaptation and speciation probably proceed through a combination of cis-regulatory and structural mutations, with a substantial contribution of the latter.

Linkage maps of the dwarf and Normal lake whitefish (Coregonus clupeaformis) species complex and their hybrids reveal the genetic architecture of population divergence

Elucidating the genetic architecture of population divergence may reveal the evolution of reproductive barriers and the genomic regions implicated in the process. We assembled genetic linkage maps for the dwarf and Normal lake whitefish species complex and their hybrids. A total of 877 AFLP loci and 30 microsatellites were positioned. The homology of mapped loci between families supported the existence of 34 linkage groups (of 40n expected) exhibiting 83% colinearity among linked loci between these two families. Classes of AFLP markers were not randomly distributed among linkage groups. Both AFLP and microsatellites exhibited deviations from Mendelian expectations, with 30.4% exhibiting significant segregation distortion across 28 linkage groups of the four linkage maps in both families (P < 0.00001). Eight loci distributed over seven homologous linkage groups were significantly distorted in both families and the level of distortion, when comparing homologous loci of the same phase between families, was correlated (Spearman R = 0.378, P = 0.0021). These results suggest that substantial divergence incurred during allopatric glacial separation and subsequent sympatric ecological specialization has resulted in several genomic regions that are no longer complementary between dwarf and Normal populations issued from different evolutionary glacial lineages.

The genetic basis of intrinsic and extrinsic post-zygotic reproductive isolation jointly promoting speciation in the lake whitefish species complex (Coregonus clupeaformis)

Understanding the genetic architecture of reproductive barriers and the evolutionary forces that drove their divergence represents a considerable challenge towards understanding speciation. The objective of this study was to determine the genetic basis of intrinsic and extrinsic post-zygotic isolation in diverging populations of dwarf and normal lake whitefish with allopatric glacial origins. We found that the rate of embryonic mortality was 5.3-6.5 times higher in dwarf-normal hybrid backcrosses during development than in F1 dwarf and normal crosses. When comparing embryos that died during development against larvae that successfully hatched, patterns of Mendelian segregation at 101 loci whose linkage is known identified 13 loci distributed over seven linkage groups that exhibited significant shifts in segregation ratios leading to significant segregation distortion at these loci in the surviving progeny. Controlled crosses and quantitative trait loci analysis revealed a significant genetic basis for developmental time until emergence, a trait critical to fish larval survival in nature. Hatching backcross progeny exhibited asynchronous emergence and transgressive segregation, suggesting that extrinsic post-zygotic isolation may select against hybridization in specific environmental contexts. Evidence of a genetic basis for increased embryonic mortality followed by asynchronous emergence indicated that intrinsic and extrinsic mechanisms are not mutually exclusive in the formation and maintenance of reproductive isolation, but may be jointly promoting population divergence and ultimately speciation.

Macroevolutionary data suggest a role for reinforcement in pollination system shifts

Reproductive isolation can evolve either as a by-product of divergent selection or through reinforcement. We used the Cape flora of South Africa, known for its high level of pollination specialization, as a model system to test the potential role of shifts in pollination system in the speciation process. Comparative analysis of 41 sister-species pairs (representing Geraniaceae, Iridaceae, and Orchidaceae) for which complete pollinator, edaphic, and distribution data are available showed that for sister species with overlapping distribution ranges, pollination system shifts are significantly associated with edaphic shifts. In contrast, there is no significant association between pollination system shifts and edaphic shifts for allopatric sister species. These results are interpreted as evidence for reinforcement.

Pervasive adaptive evolution among interactors of the Drosophila hybrid inviability gene, Nup96

Nup96 is involved in a lethal hybrid incompatibility between 2 fruit fly species, Drosophila melanogaster and Drosophila simulans. Recurrent adaptive evolution drove the rapid functional divergence of Nup96 in both the D. melanogaster and the D. simulans lineages. Functional divergence of Nup96 between these 2 species is unexpected as Nup96 encodes part of the Nup107 subcomplex, an architectural component of nuclear pore complexes, the macromolecular channels in nuclear envelopes that mediate nucleocytoplasmic traffic in all eukaryotes. Here we study the evolutionary histories of 5 of Nup96's protein interactors--3 stable Nup107 subcomplex proteins (Nup75, Nup107, and Nup133) and 2 mobile nucleoporins (Nup98 and Nup153)--and show that all 5 have experienced recurrent adaptive evolution. These results are consistent with selection-driven coevolution among molecular interactors within species causing the incidental evolution of incompatible interactions seen in hybrids between species. We suggest that genetic conflict-driven processes may have contributed to the rapid molecular evolution of Nup107 subcomplex genes.