Evolution of the Drosophila nuclear pore complex results in multiple hybrid incompatibilities

Chromosomes, conflict, and epigenetics: chromosomal speciation revisited

Since Darwin first noted that the process of speciation was indeed the "mystery of mysteries," scientists have tried to develop testable models for the development of reproductive incompatibilities-the first step in the formation of a new species. Early theorists proposed that chromosome rearrangements were implicated in the process of reproductive isolation; however, the chromosomal speciation model has recently been questioned. In addition, recent data from hybrid model systems indicates that simple epistatic interactions, the Dobzhansky-Muller incompatibilities, are more complex. In fact, incompatibilities are quite broad, including interactions among heterochromatin, small RNAs, and distinct, epigenetically defined genomic regions such as the centromere. In this review, we will examine both classical and current models of chromosomal speciation and describe the "evolving" theory of genetic conflict, epigenetics, and chromosomal speciation.

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.

Conditions for mutation-order speciation

Two models for speciation via selection have been proposed. In the well-known model of 'ecological speciation', divergent natural selection between environments drives the evolution of reproductive isolation. In a second 'mutation-order' model, different, incompatible mutations (alleles) fix in different populations adapting to the same selective pressure. How to demonstrate mutation-order speciation has been unclear, although it has been argued that it can be ruled out when gene flow occurs because the same, most advantageous allele will fix in all populations. However, quantitative examination of the interaction of factors influencing the likelihood of mutation-order speciation is lacking. We used simulation models to study how gene flow, hybrid incompatibility, selective advantage, timing of origination of new mutations and an initial period of allopatric differentiation affect population divergence via the mutation-order process. We find that at least some population divergence can occur under a reasonably wide range of conditions, even with moderate gene flow. However, strong divergence (e.g. fixation of different alleles in different populations) requires very low gene flow, and is promoted when (i) incompatible mutations have similar fitness advantages, (ii) less fit mutations arise slightly earlier in evolutionary time than more fit alternatives, and (iii) allopatric divergence occurs prior to secondary contact.

A genetic model for the female sterility barrier between Asian and African cultivated rice species

S(1) is the most important locus acting as a reproductive barrier between Oryza sativa and O. glaberrima. It is a complex locus, with factors that may affect male and female fertility separately. Recently, the component causing the allelic elimination of pollen was fine mapped. However, the position and nature of the component causing female sterility remains unknown. To fine map the factor of the S(1) locus affecting female fertility, we developed a mapping approach based on the evaluation of the degree of female transmission ratio distortion (fTRD) of markers. Through implementing this methodology in four O. sativa x O. glaberrima crosses, the female component of the S(1) locus was mapped into a 27.8-kb (O. sativa) and 50.3-kb (O. glaberrima) region included within the interval bearing the male component of the locus. Moreover, evidence of additional factors interacting with S(1) was also found. In light of the available data, a model where incompatibilities in epistatic interactions between S(1) and the additional factors are the cause of the female sterility barrier between O. sativa and O. glaberrima was developed to explain the female sterility and the TRD mediated by S(1). According to our model, the recombination ratio and allelic combinations between these factors would determine the final allelic frequencies observed for a given cross.

Step-by-step evolution of neo-sex chromosomes in geographical populations of wild silkmoths, Samia cynthia ssp

Geographical subspecies of wild silkmoths, Samia cynthia ssp. (Lepidoptera: Saturniidae), differ considerably in sex chromosome constitution owing to sex chromosome fusions with autosomes, which leads to variation in chromosome numbers. We cloned S. cynthia orthologues of 16 Bombyx mori genes and mapped them to chromosome spreads of S. cynthia subspecies by fluorescence in situ hybridization (FISH) to determine the origin of S. cynthia neo-sex chromosomes. FISH mapping revealed that the Z chromosome and chromosome 12 of B. mori correspond to the Z chromosome and an autosome (A₁) of S. c. ricini (Vietnam population, 2n=27, Z0 in female moths), respectively. B. mori chromosome 11 corresponds partly to another autosome (A₂) and partly to a chromosome carrying nucleolar organizer region (NOR) of this subspecies. The NOR chromosome of S. c. ricini is also partly homologous to B. mori chromosome 24. Furthermore, our results revealed that two A₁ homologues each fused with the W and Z chromosomes in a common ancestor of both Japanese subspecies S. c. walkeri (Sapporo population, 2n=26, neo-Wneo-Z) and S. cynthia subsp. indet. (Nagano population, 2n=25, neo-WZ₁Z₂). One homologue, corresponding to the A₂ autosome in S. c. ricini and S. c. walkeri, fused with the W chromosome in S. cynthia subsp. indet. Consequently, the other homologue became a Z₂ chromosome. These results clearly showed a step-by-step evolution of the neo-sex chromosomes by repeated autosome-sex chromosome fusions. We suggest that the rearrangements of sex chromosomes may facilitate divergence of S. cynthia subspecies towards speciation.

Multiple molecular mechanisms cause reproductive isolation between three yeast species

Incompatibility between nuclear and mitochondrial genomes in yeast species may represent a general mechanism of reproductive isolation during yeast evolution.

Nuclear-mitochondrial conflict (cytonuclear incompatibility) is a specific form of Dobzhansky-Muller incompatibility previously shown to cause reproductive isolation in two yeast species. Here, we identified two new incompatible genes, MRS1 and AIM22, through a systematic study of F2 hybrid sterility caused by cytonuclear incompatibility in three closely related Saccharomyces species (S. cerevisiae, S. paradoxus, and S. bayanus). Mrs1 is a nuclear gene product required for splicing specific introns in the mitochondrial COX1, and Aim22 is a ligase encoded in the nucleus that is required for mitochondrial protein lipoylation. By comparing different species, our result suggests that the functional changes in MRS1 are a result of coevolution with changes in the COX1 introns. Further molecular analyses demonstrate that three nonsynonymous mutations are responsible for the functional differences of Mrs1 between these species. Functional complementation assays to determine when these incompatible genes altered their functions show a strong correlation between the sequence-based phylogeny and the evolution of cytonuclear incompatibility. Our results suggest that nuclear-mitochondrial incompatibility may represent a general mechanism of reproductive isolation during yeast evolution.

Hybrids between species are usually inviable or sterile, possibly due to functional incompatibility between genes from the different species. Incompatible genes are hypothesized to encode interacting components that cannot function properly when paired with alleles from another species. To understand how incompatible gene pairs result in hybrid sterility or inviability, it is important to identify these genes and reconstruct their evolutionary history. A previous study has shown that incompatibility between nuclear and mitochondrial genomes (cytonuclear incompatibility) causes hybrid sterility between two yeast species. To expand on these findings, we screened three yeast species for genes involved in cytonuclear incompatibility, discovering two nuclear genes, MRS1 and AIM22, which encode proteins that are unable to support full mitochondrial function in the hybrids. Of these two genes, Mrs1 is required for removing a specific intron in the mitochondrial COX1 gene. By comparing different yeast species, we find a clear coevolutionary relationship between Mrs1 function and the COX1 intron pattern. We also show that changes in three amino acids in the Mrs1 RNA-binding domain are sufficient to make Mrs1 incompatible in hybrids. Our results suggest that cytonuclear incompatibility may represent a general mechanism of reproductive isolation during yeast evolution.

Introgression of Drosophila simulans nuclear pore protein 160 in Drosophila melanogaster alone does not cause inviability but does cause female sterility

We have been analyzing genes for reproductive isolation by replacing Drosophila melanogaster genes with homologs from Drosophila simulans by interspecific backcrossing. Among the introgressions established, we found that a segment of the left arm of chromosome 2, Int(2L)S, carried recessive genes for hybrid sterility and inviability. That nuclear pore protein 160 (Nup160) in the introgression region is involved in hybrid inviability, as suggested by others, was confirmed by the present analysis. Male hybrids carrying an X chromosome of D. melanogaster were not rescued by the Lethal hybrid rescue (Lhr) mutation when the D. simulans Nup160 allele was made homozygous or hemizygous. Furthermore, we uniquely found that Nup160 is also responsible for hybrid sterility. Females were sterile when D. simulans Nup160 was made homozygous or hemizygous in the D. melanogaster genetic background. Genetic analyses indicated that the D. simulans Nup160 introgression into D. melanogaster was sufficient to cause female sterility but that other autosomal genes of D. simulans were also necessary to cause lethality. The involvement of Nup160 in hybrid inviability and female sterility was confirmed by transgene experiment.

Determinants of divergent adaptation and Dobzhansky-Muller interaction in experimental yeast populations

Divergent adaptation can be associated with reproductive isolation in speciation [1]. We recently demonstrated the link between divergent adaptation and the onset of reproductive isolation in experimental populations of the yeast Saccharomyces cerevisiae evolved from a single progenitor in either a high-salt or a low-glucose environment [2]. Here, whole-genome resequencing and comparative genome hybridization of representatives of three populations revealed 17 mutations, six of which explained the adaptive increases in mitotic fitness. In two populations evolved in high salt, two different mutations occurred in the proton efflux pump gene PMA1 and the global transcriptional repressor gene CYC8; the ENA genes encoding sodium efflux pumps were overexpressed once through expansion of this gene cluster and once because of mutation in the regulator CYC8. In the population from low glucose, one mutation occurred in MDS3, which modulates growth at high pH, and one in MKT1, a global regulator of mRNAs encoding mitochondrial proteins, the latter recapitulating a naturally occurring variant. A Dobzhansky-Muller (DM) incompatibility between the evolved alleles of PMA1 and MKT1 strongly depressed fitness in the low-glucose environment. This DM interaction is the first reported between experimentally evolved alleles of known genes and shows how reproductive isolation can arise rapidly when divergent selection is strong.

Speciation genetics: current status and evolving approaches

The view of species as entities subjected to natural selection and amenable to change put forth by Charles Darwin and Alfred Wallace laid the conceptual foundation for understanding speciation. Initially marred by a rudimental understanding of hereditary principles, evolutionists gained appreciation of the mechanistic underpinnings of speciation following the merger of Mendelian genetic principles with Darwinian evolution. Only recently have we entered an era where deciphering the molecular basis of speciation is within reach. Much focus has been devoted to the genetic basis of intrinsic postzygotic isolation in model organisms and several hybrid incompatibility genes have been successfully identified. However, concomitant with the recent technological advancements in genome analysis and a newfound interest in the role of ecology in the differentiation process, speciation genetic research is becoming increasingly open to non-model organisms. This development will expand speciation research beyond the traditional boundaries and unveil the genetic basis of speciation from manifold perspectives and at various stages of the splitting process. This review aims at providing an extensive overview of speciation genetics. Starting from key historical developments and core concepts of speciation genetics, we focus much of our attention on evolving approaches and introduce promising methodological approaches for future research venues.

Hybrid incompatibility genes: remnants of a genomic battlefield?

Hybrid incompatibility (including sterility, lethality, and less extreme negative effects) interests evolutionary biologists because of its role in speciation as a reproductive isolating barrier. It also has unusual genetic properties, being mainly due to interactions between at least two genes. Recent studies have identified some of the interacting genes that underlie hybrid incompatibility. These genes represent a wide array of functions, including those involved in oxidative respiration, nuclear trafficking, DNA-binding, and plant defense. Accumulating evidence suggests genomic conflict frequently drives the divergence causing incompatibilities in hybrids. The evidence bearing on this genomic conflict hypothesis is assessed and ways to test it conclusively are suggested.

Molecular mechanisms of polyploidy and hybrid vigor

Hybrids such as maize (Zea mays) or domestic dog (Canis lupus familiaris) grow bigger and stronger than their parents. This is also true for allopolyploids such as wheat (Triticum spp.) or frog (i.e. Xenopus and Silurana) that contain two or more sets of chromosomes from different species. The phenomenon, known as hybrid vigor or heterosis, was systematically characterized by Charles Darwin (1876). The rediscovery of heterosis in maize a century ago has revolutionized plant and animal breeding and production. Although genetic models for heterosis have been rigorously tested, the molecular bases remain elusive. Recent studies have determined the roles of nonadditive gene expression, small RNAs, and epigenetic regulation, including circadian-mediated metabolic pathways, in hybrid vigor, which could lead to better use and exploitation of the increased biomass and yield in hybrids and allopolyploids for food, feed, and biofuels.

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.

Mitochondrial gene in the nuclear genome induces reproductive barrier in rice

Hybrid incompatibility in F(1) hybrids or later generations is often observed as sterility or inviability. This incompatibility acts as postzygotic reproductive isolation, which results in the irreversible divergence of species. Here, we show that the reciprocal loss of duplicated genes encoding mitochondrial ribosomal protein L27 causes hybrid pollen sterility in F(1) hybrids of the cultivated rice Oryza sativa and its wild relative O. glumaepatula. Functional analysis revealed that this gene is essential for the later stage of pollen development, and distribution analysis suggests that the gene duplication occurred before the divergence of the AA genome species. On the basis of these results, we discuss the possible contribution of the "founder effect" in establishing this reproductive barrier.

A novel method to detect proteins evolving at correlated rates: identifying new functional relationships between coevolving proteins

Interacting proteins evolve at correlated rates, possibly as the result of evolutionary pressures shared by functional groups and/or coevolution between interacting proteins. This evolutionary signature can be exploited to learn more about protein networks and to infer functional relationships between proteins on a genome-wide scale. Multiple methods have been introduced that detect correlated evolution using amino acid distances. One assumption made by these methods is that the neutral rate of nucleotide substitution is uniform over time; however, this is unlikely and such rate heterogeneity would adversely affect amino acid distance methods. We explored alternative methods that detect correlated rates using protein-coding nucleotide sequences in order to better estimate the rate of nonsynonymous substitution at each branch (d(N)) normalized by the underlying synonymous substitution rate (d(S)). Our novel likelihood method, which was robust to realistic simulation parameters, was tested on Drosophila nuclear pore proteins, which form a complex with well-documented physical interactions. The method revealed significantly correlated evolution between nuclear pore proteins, where members of a stable subcomplex showed stronger correlations compared with those proteins that interact transiently. Furthermore, our likelihood approach was better able to detect correlated evolution among closely related species than previous methods. Hence, these sequence-based methods are a complementary approach for detecting correlated evolution and could be applied genome-wide to provide candidate protein-protein interactions and functional group assignments using just coding sequences.

The nuclear pore complex has entered the atomic age

Nuclear pore complexes (NPCs) perforate the nuclear envelope and represent the exclusive passageway into and out of the nucleus of the eukaryotic cell. Apart from their essential transport function, components of the NPC have important, direct roles in nuclear organization and in gene regulation. Because of its central role in cell biology, it is of considerable interest to determine the NPC structure at atomic resolution. The complexity of these large, 40-60 MDa protein assemblies has for decades limited such structural studies. More recently, exploiting the intrinsic modularity of the NPC, structural biologists are making progress toward understanding this nanomachine in molecular detail. Structures of building blocks of the stable, architectural scaffold of the NPC have been solved, and distinct models for their assembly proposed. Here we review the status of the field and lay out the challenges and the next steps toward a full understanding of the NPC at atomic resolution.

Along the speciation continuum in sticklebacks

Speciation can be viewed as a continuum, potentially divisible into several states: (1) continuous variation within panmictic populations, (2) partially discontinuous variation with minor reproductive isolation, (3) strongly discontinuous variation with strong but reversible reproductive isolation and (4) complete and irreversible reproductive isolation. Research on sticklebacks (Gasterosteidae) reveals factors that influence progress back and forth along this continuum, as well as transitions between the states. Most populations exist in state 1, even though some of these show evidence of disruptive selection and positive assortative mating. Transitions to state 2 seem to usually involve strong divergent selection coupled with at least a bit of geographic separation, such as parapatry (e.g. lake and stream pairs and mud and lava pairs) or allopatry (e.g. different lakes). Transitions to state 3 can occur when allopatric or parapatric populations that evolved under strong divergent selection come into secondary contact (most obviously the sympatric benthic and limnetic pairs), but might also occur between populations that remained in parapatry or allopatry. Transitions to state 4 might be decoupled from these selective processes, because the known situations of complete, or nearly complete, reproductive isolation (Japan Sea and Pacific Ocean pair and the recognized gasterosteid species) are always associated with chromosomal rearrangements and environment-independent genetic incompatibilities. Research on sticklebacks has thus revealed complex and shifting interactions between selection, adaptation, mutation and geography during the course of speciation.

Genome-wide investigation of reproductive isolation in experimental lineages and natural species of Neurospora: identifying candidate regions by microarray-based genotyping and mapping

Inherent incompatibilities between genetic components from genomes of different species may cause intrinsic reproductive isolation. In evolution experiments designed to instigate speciation in laboratory populations of the filamentous fungus Neurospora, we previously discovered a pair of incompatibility loci (dfe and dma) that interact negatively to cause severe defects in sexual reproduction. Here we show that the dfe-dma incompatibility also is a significant cause of genetic isolation between two naturally occurring species of Neurospora (N. crassa and N. intermedia). The strong incompatibility interaction has a simple genetic basis (two biallelic loci) and antagonistic epistasis occurs between heterospecific alleles only, consistent with the Dobzhansky-Muller model of genic incompatibility. We developed microarray-based, restriction-site associated DNA (RAD) markers that identified approximately 1500 polymorphisms between the genomes of the two species, and constructed the first interspecific physical map of Neurospora. With this new mapping resource, the approximate genomic locations of the incompatibility loci were determined using three different approaches: genome scanning, bulk-segregant analyses, and introgression. These population, quantitative, and classical genetics methods concordantly identified two candidate regions, narrowing the search for each incompatibility locus to only approximately 2% of the nuclear genome. This study demonstrates how advances in high-throughput, genome-wide genotyping can be applied to mapping reproductive isolation genes and speciation research.

Chromosomal inversions and species differences: when are genes affecting adaptive divergence and reproductive isolation expected to reside within inversions?

Many factors can promote speciation, and one which has received much attention is chromosomal inversions. A number of models propose that the recombination suppressing effects of inversions facilitate the maintenance of differences between interbreeding populations in genes affecting adaptive divergence and reproductive isolation. These models predict that such genes will disproportionately reside within inversions, rather than in collinear regions. This hypothesis has received some support, but exceptions exist. Additionally, the effects of known low levels of recombination within inversions on these models are uninvestigated. Here, simulations are used to compare the maintenance of genetic differences between populations following secondary contact and hybridization in different inversion models. We compare regions with no recombination within them to regions with low recombination and to collinear regions with free recombination. Our most general finding is that the low levels of recombination within an inversion often result in the loss of accentuated divergence in inverted regions compared to collinear ones. We conclude that inversions can facilitate the maintenance of species differences under some conditions, but that large or qualitative differences between inverted and collinear regions need not occur. We also find that strong selection facilitates maintenance of divergence in a manner analogous to inversions.

Genetics and ecological speciation

Species originate frequently by natural selection. A general mechanism by which this occurs is ecological speciation, defined as the evolution of reproductive isolation between populations as a result of ecologically-based divergent natural selection. The alternative mechanism is mutation-order speciation in which populations fix different mutations as they adapt to similar selection pressures. Although numerous cases now indicate the importance of ecological speciation in nature, very little is known about the genetics of the process. Here, we summarize the genetics of premating and postzygotic isolation and the role of standing genetic variation in ecological speciation. We discuss the role of selection from standing genetic variation in threespine stickleback (Gasterosteus aculeatus), a complex of species whose ancestral marine form repeatedly colonized and adapted to freshwater environments. We propose that ecological speciation has occurred multiple times in parallel in this group via a "transporter" process in which selection in freshwater environments repeatedly acts on standing genetic variation that is maintained in marine populations by export of freshwater-adapted alleles from elsewhere in the range. Selection from standing genetic variation is likely to play a large role in ecological speciation, which may partly account for its rapidity.

Parental squabbles and genome expression: lessons from the polyploids

The merger of evolutionarily diverged genomes to form a new polyploid genetic system can involve extensive remodeling of gene regulation. A recent paper in BMC Biology provides important insights into regulatory events that have affected the evolution of allopolyploid cotton.

Genetics and lineage-specific evolution of a lethal hybrid incompatibility between Drosophila mauritiana and its sibling species

The Dobzhansky-Muller model posits that intrinsic postzygotic reproductive isolation--the sterility or lethality of species hybrids--results from the evolution of incompatible epistatic interactions between species: favorable or neutral alleles that become fixed in the genetic background of one species can cause sterility or lethality in the genetic background of another species. The kind of hybrid incompatibility that evolves between two species, however, depends on the particular evolutionary history of the causative substitutions. An allele that is functionally derived in one species can be incompatible with an allele that is functionally derived in the other species (a derived-derived hybrid incompatibility). But an allele that is functionally derived in one species can also be incompatible with an allele that has retained the ancestral state in the other species (a derived-ancestral hybrid incompatibility). The relative abundance of such derived-derived vs. derived-ancestral hybrid incompatibilities is unknown. Here, we characterize the genetics and evolutionary history of a lethal hybrid incompatibility between Drosophila mauritiana and its two sibling species, D. sechellia and D. simulans. We show that a hybrid lethality factor(s) in the pericentric heterochromatin of the D. mauritiana X chromosome, hybrid lethal on the X (hlx), is incompatible with a factor(s) in the same small autosomal region from both D. sechellia and D. simulans, Suppressor of hlx [Su(hlx)]. By combining genetic and phylogenetic information, we infer that hlx-Su(hlx) hybrid lethality is likely caused by a derived-ancestral incompatibility, a hypothesis that can be tested directly when the genes are identified.