Pervasive Adaptive Evolution among Interactors

Pervasive Adaptive Evolution among Interactors of the Drosophila Hybrid
Inviability Gene, Nup96
Daven C. Presgraves* and Wolfgang Stephan*
*Section of Evolutionary Biology, Biocenter, University of Munich, Planegg-Martinsried, Germany; and
Department of Biology, University of Rochester
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.
Introduction
New species arise through the evolution of reproductive isolation (Dobzhansky 1937; Coyne and Orr 2004).
Among diverging allopatric Drosophila species, prezygotic
isolation (e.g., sexual isolation) and intrinsic postzygotic
isolation (e.g., hybrid sterility or inviability) evolve gradually and do so at roughly equal rates (Coyne and Orr 1989,
1997). The evolution of intrinsic postzygotic isolation corresponds to the accumulation of alleles that have deleterious
epistatic interactions in species hybrids—so-called hybrid
incompatibilities (Dobzhansky 1937; Muller 1940, 1942;
Orr 1995). Many classical genetic analyses have confirmed
that alleles fixed in 1 species are often functionally incompatible with the genetic backgrounds of closely related species (Coyne and Orr 2004). However, the normal functions
and evolutionary histories of hybrid incompatibility genes
have been determined in only a few cases, in part because
the particular loci involved have proven difficult to identify
(Orr et al. 2004; Wu and Ting 2004; Orr 2005).
In Drosophila, just 4 hybrid incompatibility genes have
been identified. Ods is a duplicated homeobox transcription factor that is misexpressed in the testes of Drosophila
simulans–Drosophila mauritiana hybrids, causing male
sterility (Ting et al. 1998; Barbash and Ashburner 2003;
Barbash et al. 2003; Ting et al. 2004). Hmr is a single-copy
gene encoding a Myb-related protein with putative DNAbinding properties, which causes lethality and female sterility in F1 Drosophila melanogaster–D. simulans hybrids
(Barbash and Ashburner 2003; Barbash et al. 2003). The
D. simulans allele of Nup96, a single-copy gene, encodes
a nuclear pore protein that causes recessive lethality when
combined with a hemizygous D. melanogaster X chromosome (Presgraves et al. 2003). And, JYAlpha, which encodes
a male fertility–essential Na1/K1 ATPase, was recently
found to cause sterility due to its change in genetic map position: JYAlpha is on the fourth chromosome in D. melanogaster but on the third chromosome in D. simulans; hybrids
Key words: Drosophila, hybrid inviability, hybrid incompatibility,
speciation.
E-mail: [email protected].
Mol. Biol. Evol. 24(1):306–314. 2007
doi:10.1093/molbev/msl157
Advance Access publication October 20, 2006
Ó The Author 2006. Published by Oxford University Press on behalf of
the Society for Molecular Biology and Evolution. All rights reserved.
For permissions, please e-mail: [email protected]
homozygous D. melanogaster at all major chromosomes
except the fourth are therefore sterile as they lack JYalpha
(Masly et al. 2006).
It is unclear what role, if any, these particular genes
played in reproductive isolation in nature, but they nevertheless serve as important models whose molecular properties and evolutionary histories are informative about the
evolution and genetics of hybrid incompatibilities in general. For example, a striking, but still tentative, pattern
has already emerged from 3 of these genes: each has a history of recurrent positive selection at its protein-coding sequence (Ting et al. 1998; Presgraves et al. 2003; Barbash
et al. 2004). (The remaining gene, JYAlpha, causes hybrid
sterility not by its functional divergence but by a change in
map position [Lynch and Force 2000]). Adaptive divergence at Ods is confined to the D. mauritiana lineage (Ting
et al. 1998), whereas Hmr and Nup96 evolved adaptively
in both the D. melanogaster and the D. simulans lineages
(Presgraves et al. 2003; Barbash et al. 2004).
Hybrid incompatibility genes are not unconditionally
deleterious but instead depend on alleles at interacting loci.
By identifying specific loci, we can make inferences about
the interactions that lead to hybrid sterility and inviability.
For example, the D. simulans allele of Nup96 causes lethality in hybrids that are hemizygous for the D. melanogaster X
chromosome but not in hybrids hemizygous for the D. simulans X (Presgraves 2003; Presgraves et al. 2003). This is
not surprising as the D. simulans allele of Nup96 is not expected to be incompatible with loci on the D. simulans X
chromosome. On the contrary, the functional divergence
of Nup96 in D. simulans must have coevolved with functional divergence at 1 (or more) of its interacting genes. Because Nup96 does not encode a transcription factor and
because the relevant divergence in Nup96 seems to be concentrated in its protein-coding region (Presgraves et al.
2003), its hybrid inviability likely involves a protein–
protein incompatibility rather than a protein–DNA incompatibility (Landry et al. 2005).
Fortunately, a good deal is known about Nup96’s function and potential protein interactors. Nup96 is 1 of ;30 different nucleoporins that together form a single nuclear pore
complex (NPC) (Allen et al. 2000, 2002; Rout et al. 2000;
Tran and Wente 2006). Thousands of NPCs perforate the
nuclear envelopes of eukaryotic cells, mediating molecular
Adaptive Evolution of Nucleoporins 307
traffic into and out of the nucleus. The NPC is largely conserved throughout eukaryotic evolution in its modular
architecture, complement of nucleoporins, and interactions
among particular subsets of nucleoporins (Allen et al. 2000,
2002; Rout et al. 2000; Vasu and Forbes 2001; Bapteste et al.
2005; Devos et al. 2006; Tran and Wente 2006). Nup96 is 1
of 8 known proteins that form the Nup107 subcomplex
(known as the Nup84 subcomplex in yeast), a stable structural component of the NPC, along with Nup75, Nup107,
Nup133, Nup160, Seh1, Nup37, and Nup43 (Siniossoglou
et al. 1996; Rappsilber et al. 2000; Siniossoglou et al. 2000;
Vasu et al. 2001; Allen et al. 2002; Lutzmann et al. 2002;
Walther et al. 2003; Devos et al. 2004). From what is now
known, Nup96 occupies a central position in the Nup107
subcomplex and interacts directly with Nup75 and
Nup107 (Lutzmann et al. 2002). In addition to these stable
Nup107 subcomplex proteins, 3 dynamic nucleoporins
(Nup98, Nup153, and Sec13) are known to shuttle between
the nuclear interior and the NPC, where they physically dock
at the Nup107 subcomplex (Griffis et al. 2002). All 3 of these
dynamic nucleoporins interact directly with Nup96 (Vasu
et al. 2001; Hodel et al. 2002; Enninga et al. 2003);
Nup153 and Nup98 are also known to bind with Nup107,
Nup133, and Nup160 (Vasu et al. 2001).
Because Nup96 is a hub for multiple interactions in the
Nup107 subcomplex—and because Nup96 has experienced
repeated bouts of adaptive evolution in the D. melanogaster
and the D. simulans lineages—we investigated the possibility of adaptive coevolution between Nup96 and some of its
interacting proteins. We surveyed DNA sequence divergence and polymorphism in D. melanogaster and D. simulans at 5 additional nucleoporins: Nup75, Nup107, Nup133,
Nup98, and Nup153. Here we report that divergence at every
1 of these genes was driven by adaptive evolution.
Materials and Methods
Genes and Fly Stocks
We identified Nup75 (CG5733), Nup107 (CG6743),
Nup133 (CG6958), Nup98 (CG10198), and Nup153
(CG4453) using Blast searches of the D. melanogaster
and D. simulans whole-genome sequences using mouse,
human, Xenopus, and yeast homologs as queries. All of
these nucleoporins appear to be encoded by single-copy
genes in the 2 Drosophila genomes. For the 2 large shuttling
nucleoporins, we sequenced 2.4 kb from the 5# end Nup153
(positions 27–2387 of 9246 bp) and 1.3 kb from 3# end of
Nup98 (positions 1906–3223 of 3223 bp) as these are
known to encode the protein domains that mediate interactions with the Nup107 subcomplex (Enarson et al. 1998;
Vasu et al. 2001; Hodel et al. 2002). For the 3 Nup107
subcomplex nucleoporins, we sequenced their entire protein-coding regions. For completeness, we also included
the previously sequenced entire protein-coding region of
Nup96 in our analyses below (Presgraves et al. 2003).
For each gene, we sampled 11–15 chromosomes from
Zimbabwean populations of both D. melanogaster and D.
simulans. The lines were kindly provided by Dr Charles
Aquadro (Cornell University). We studied sequences from
Zimbabwe to minimize the confounding effects of the recent demographic expansion associated with the out-of-
Africa events of both species (Glinka et al. 2003; Baudry
et al. 2004; Ometto et al. 2005; Thornton and Andolfatto
2006). Toextract and isogenize autosomesinD.melanogaster,
we used standard crosses involving balancer chromosomes.
For those lines for which autosomes could not be made homozygous (e.g., due to the presence of balanced recessive
lethals), we extracted genomic DNA from deficiency heterozygotes in which the deficiency deletes a small chromosomal region including the relevant nucleoporin (i.e., from
flies hemizygous for a Zimbabwe chromosome at the locus
of interest). Details about the particular stocks and crosses
used for each nucleoporin are available upon request.
We obtained Drosophila sechellia and Drosophila yakuba sequences for each gene using Blast searches of the
complete genome sequences available from the Broad Institute and the Washington University Genome Sequencing
Center, respectively.
Sequencing and Analysis
We extracted genomic DNA from ;20 females from
each line using Puregene DNA isolation kits (Gentra Systems, Minneapolis, MN). Polymerase chain reaction and internal sequencing primers were designed using the
available D. melanogaster (Adams et al. 2000) and D. simulans (Washington University Genome Sequencing Center)
genome sequences as guides. We purified PCR products
with EXOSAP-IT (USB, Cleveland, OH) and sequenced
both strands. Some sequencing was done at the University
of Munich, using the DYEnamic ET terminator cycle sequencing kit (Amersham Biosciences, Buckinghamshire,
UK), and some at the University of Rochester, using
ABI BigDye terminator chemistry (PerkinElmer, Wellesley, MA). Sequences were assembled and edited using
Sequencher v. 4.5 (Gene Codes, Ann Arbor, MI) and
manually aligned with Se-Al v. 2.0 (Rambaut 1996). All
sequences have been deposited in GenBank, under accession numbers EF057923–EF058045.
We used DnaSP v. 4.0 for most population genetic
analyses (Rozas et al. 2003). To estimate probabilities for
Tajima’s D (Tajima 1989) and Fay and Wu’s H (Fay and
Wu 2000) statistics (with and without recombination), we
simulated 5,000 neutral genealogies, conditioning on the
observed hw (Watterson 1975), using J. Fay’s web application (http://www.genetics.wustl.edu/jflab/htest.html). The
population recombination rate, q 5 4Nr (Hudson 1987),
used in these simulations was estimated from the data.
Results
Adaptive Evolution of the Nup107 Subcomplex
We used MK tests (McDonald and Kreitman 1991) to
determine if patterns of protein polymorphism and divergence at 6 nucleoporins, Nup75, Nup107, Nup133,
Nup98, Nup153, and, for completeness, Nup96 (Presgraves
et al. 2003), are consistent with histories of neutral mutation
and drift. Under a neutral model of molecular evolution, in
which only neutral mutations contribute to variation within
and divergence between species, the ratio of replacement to
synonymous polymorphisms will equal the ratio of replacement to synonymous fixed differences (Kimura 1983;
308 Presgraves and Stephan
Table 1
MK Tests for Nucleoporins in Drosophila melanogaster and Drosophila simulans
Polymorphisms
Gene
Nup75
Nup96
Nup107
Nup133
Nup98
Nup153
Pooled
Fixed Differences
R
S
R/S Ratio
R
S
R/S Ratio
G
P Value
24
27
29
29
14
10
133
95
108
105
127
44
32
511
0.253
0.250
0.276
0.228
0.318
0.313
0.260
27
27
35
34
14
49
186
37
34
41
33
17
49
211
0.730
0.794
0.854
1.030
0.824
1.000
0.882
9.764
11.888
13.365
22.813
4.054
8.647
78.192
0.0018
0.0006
0.0003
,0.0001
0.0441
0.0033
,10 15
McDonald and Kreitman 1991). Table 1 shows that this null
neutral model can be rejected for all 6 genes. Each—the
4 stable Nup107 subcomplex proteins (Nup75, Nup96,
Nup107, and Nup133) and the 2 dynamic shuttling nucleoporins (Nup98 and Nup153)—shows a significant excess
of amino acid replacement fixations between species, consistent with recurrent bouts of adaptive protein evolution.
Sliding window analyses of replacement/synonymous
(R/S) polymorphism and divergence are provided as supplementary figures, Supplementary Material online.
Lineage-Specific Adaptive Evolution
To determine whether the adaptive evolution of any
particular nucleoporin gene was restricted to the D. melanogaster lineage or to the D. simulans lineage, or occurred
in both lineages, we performed lineage-specific MK tests.
Table 2 shows the results of 3 of these tests for each gene.
For the first test, D. simulans and D. yakuba sequences were
used as near and distant outgroups, respectively, to polarize
mutations in the D. melanogaster lineage. For the second
test, D. melanogaster and D. yakuba were used as near and
distant outgroups, respectively, to polarize mutations in the
D. simulans lineage. These 2 tests include fixations that occurred over the 3–5 Myr of independent evolution in D.
melanogaster and in D. simulans, respectively (Hey and
Kliman 1993; Li 1997; Tamura et al. 2004). For the third
test, D. sechellia and D. melanogaster sequences were used
as near and distant outgroups, respectively, to polarize substitutions in the most recent ;400,000 years of evolution in
the D. simulans lineage (Kliman et al. 2000). The latter test,
therefore, asks about patterns of substitution in the more
recent history of D. simulans. (Speciation between D. simulans and D. sechellia was sufficiently recent that some of
the history of the alleles in our sample could predate the
species split [Kliman et al. 2000], possibly confounding polymorphism and divergence. The D. sechellia–D. simulans
MK tests should therefore be interpreted with caution.) It
is important to note that lineage-specific MK tests suffer reduced statistical power, for 2 reasons. First, mapping mutations onto each lineage reduces the number of observations
in each cell of the MK table. Second, mutations that cannot
be accurately polarized, due to, for example, multiple substitutions or alignment gaps that make inference of ancestral and derived states ambiguous, are excluded from the
analysis.
Table 2
Lineage-Specific MK Tests for Nucleoporins in Drosophila melanogaster and Drosophila simulans
Polymorphisms
Gene
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
Nup75
Nup96
Nup107
Nup133
Nup98
Nup153
Pooled
Contrast
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
D.
melanogaster lineage
simulans lineage
simulans lineage (recent)
melanogaster lineage
simulans lineage
simulans lineage (recent)
melanogaster lineage
simulans lineage
simulans lineage (recent)
melanogaster lineage
simulans lineage
simulans lineage (recent)
melanogaster lineage
simulans lineage
simulans lineage (recent)
melanogaster lineage
simulans lineage
simulans lineage (recent)
melanogaster lineage
simulans lineage
simulans lineage (recent)
* Fisher’s exact probability.
Fixed Differences
R
S
R/S Ratio
R
S
R/S Ratio
G
P Value
2
17
15
4
17
16
16
10
8
12
9
8
6
7
7
6
2
2
46
62
56
21
59
51
35
58
56
37
60
59
49
53
51
18
18
18
12
16
18
172
264
253
0.095
0.288
0.294
0.114
0.293
0.286
0.432
0.167
0.136
0.245
0.170
0.157
0.333
0.389
0.389
0.500
0.125
0.111
0.267
0.235
0.221
9
10
1
12
10
0
6
19
4
10
13
2
4
4
3
12
28
7
53
84
17
19
14
1
19
7
1
27
12
0
15
16
4
9
9
0
26
19
1
115
77
7
0.474
0.714
1.000
0.632
1.429
0.000
0.222
1.583
—
0.667
0.813
0.500
0.444
0.444
—
0.462
1.474
7.000
0.461
1.091
2.429
4.427
3.258
0.682
8.082
8.046
—
1.159
22.314
—
3.664
9.409
1.342
0.141
0.032
—
0.017
13.744
16.133
5.397
54.732
28.816
0.0354
0.0711
0.4090
0.0045
0.0046
1.0000*
0.2071
,0.0001
0.0005*
0.0556
0.0022
0.2467
0.7076
0.8587
0.0366*
0.8957
0.0002
,0.0001
0.0202
,0.0001
,0.0001
Adaptive Evolution of Nucleoporins 309
Table 3
DNA Polymorphism and Divergence at Nucleoporins in Drosophila melanogaster and Drosophila simulans Populations
D. melanogaster
Nup75
Nup96
Nup107
Nup133
Nup98
Nup153
Mean
Map Positiona
bp
n
pb
55B (2R)
95B (3R)
32A (2L)
94C (3R)
95B (3R)
14F (X)
2286
3398
2560
4185
1308
2680
2736
12
15
12
11
15
11
12.7
0.0151
0.0145
0.0163
0.0258
0.0266
0.0104
0.0181
D. simulans
Dc
0.202
0.638
1.189
0.291
0.471
0.994
0.022
n
pb
12
15
11
13
15
11
12.8
0.0426
0.0236
0.0374
0.0221
0.0207
0.0132
0.0266
Dc
0.547
0.746
0.079
0.350
0.615
0.460
0.466
Kad
KSd
Ka/Ks
0.022
0.017
0.024
0.016
0.018
0.030
0.021
0.156
0.101
0.141
0.084
0.108
0.107
0.116
0.138
0.168
0.169
0.191
0.170
0.281
0.186
a
Cytogical position.
p for silent sites, with X-linked loci corrected by multiplying by 4/3.
c
Tajima’s D for silent sites.
d
Ka and Ks are Jukes–Cantor corrected.
b
In the D. melanogaster lineage, 4 of 6 nucleoporins
have experienced an excess of replacement fixations (i.e.,
R/S for fixations is greater than R/S for polymorphism).
Two of 4, Nup75 and Nup96, show individually significant
excesses of replacement fixations consistent with adaptive
evolution (table 2). A third gene, Nup133, shows a marginally nonsignificant excess of replacement fixations.
In the D. simulans lineage, 4 of 6 nucleoporins show individually significant evidence for adaptive evolution
(Nup96, Nup107, Nup133, and Nup153; table 2). Similar histories of adaptive evolution cannot be excluded for the remaining 2 genes because Nup75 shows excess amino acid
divergence that is marginally nonsignificant (table 2, line
15); and Nup98 shows a significant excess of replacement
fixations since the split from D. sechellia (table 2, line 15).
Thus, at least 5 (nearly 6) of 6 nucleoporins show some evidence of recurrent adaptive protein evolution in D. simulans.
Nup107 and Nup153 have experienced recurrent bouts
of positive selection strictly limited to the D. simulans lineage (table 2, lines 8 and 17), whereas none of the genes
show D. melanogaster–limited bouts of positive selection.
Only 1 gene, Nup96, shows strong evidence for parallel
bouts of positive selection in both species (table 2, lines
4 and 5). The data for Nup75 and Nup133 are suggestive
of similar parallel bouts of positive selection, but in 1
of the species (D. simulans and D. melanogaster, respectively), the lineage-specific tests are marginally nonsignificant. Overall, the Nup107 subcomplex appears to have
diverged more in D. simulans than in the D. melanogaster
(table 2, lines 19–21).
inferred larger effective population size (Aquadro et al.
1988; Akashi 1996), D. simulans shows ;50% higher
levels of silent diversity, on average, than D. melanogaster
(table 3). None of the nucleoporin genes in either species
exhibits the dramatically reduced variability expected after
a recent, strong selective sweep.
The frequency spectra of silent polymorphisms, summarized by Tajima’s D statistic, are more skewed toward
rare alleles in D. simulans than in D. melanogaster for 5
of 6 genes. A comparison of the distribution of all silent
polymorphisms pooled across loci differs significantly between species (Kolmogorov–Smirnov test, P 5 0.0002).
This species difference is consistent with either a stronger
signal of population expansion or stronger selective constraints on silent sites in D. simulans than in D. melanogaster (Akashi 1995). The sole exception to this pattern is
Nup107: in D. melanogaster, silent variability at Nup107
shows a negative Tajima’s D (table 3; without and with recombination, P 5 0.092 and P 5 0.013, respectively) and
a negative Fay and Wu’s H (without and with recombination, P 5 0.114 and P 5 0.057, respectively). These patterns of variability are consistent with a fairly old selective
sweep at Nup107 in D. melanogaster (relative to the approximate detection limit of ;0.1Ne generations in Drosophila), but it is difficult to rule out alternative (e.g.,
demographic) explanations. Thus, neither the levels nor
the frequency spectra of silent variability at Nup107 subcomplex genes shows particularly compelling evidence
for recent, strong selective sweeps in either species.
Discussion
Little Evidence for Recent Selective Sweeps
Recent selective sweeps caused by the fixation of new,
favorable mutations can leave characteristic signatures in
the level and distribution of linked neutral nucleotide variability (Maynard Smith and Haigh 1974; Kaplan et al.
1989; Nielsen 2005), for example, by depressing silent variability and causing an excess of rare variants. Table 3
shows summary statistics of the levels of silent (synonymous and noncoding) variability and of the frequency spectra of silent mutations for each of the 6 nucleoporins. Levels
of silent polymorphism at each locus are not unusual for
African populations of these species (Andolfatto 2001;
Mousset and Derome 2004). Consistent with its previously
The structure and function of NPCs are largely conserved among eukaryotes. It is therefore surprising that
all 4 Nup107 subcomplex nucleoporins and both interacting
shuttling nucleoporins studied have histories of recurrent
adaptive evolution during the last 3–5 Myr of divergence
between D. melanogaster and D. simulans. Several recent
analyses of polymorphism and divergence data across loci
in Drosophila have revealed an overall excess of amino
acid divergence consistent widespread adaptive evolution
(Fay et al. 2002; Smith and Eyre-Walker 2002; Sawyer
et al. 2003). However, even in such multilocus surveys,
only 5–10% of genes show individually significant evidence
for adaptive evolution (Smith and Eyre-Walker 2002;
310 Presgraves and Stephan
Presgraves 2005). The fact that each of 6 genes in our survey
shows individually significant evidence for adaptive evolution is therefore striking. This finding—along with the fact
that these proteins are known interactors—suggests that the
bouts of adaptation at these 6 genes are not independent but
rather reflect histories of molecular coevolution. We note,
however, that our evidence is circumstantial as only functional tests can confirm coevolution. These tests require that
the particular functional domains of these nucleoporins become better characterized.
What might cause such a large fraction of the proteins
of the otherwise conserved NPC to diverge so dramatically?
The pattern of divergence observed here suggests that a
simple scenario of ecological adaptation to the extrinsic
environment—say, adaptation to a new host plant—is unlikely. D. melanogaster and D. simulans are thought to have
evolved in allopatry, in southeastern sub-Saharan Africa and
in Madagascar, respectively (Lachaise et al. 1988; Dean and
Ballard 2004; Lachaise and Silvain 2004; Lachaise et al.
2004). For an ecological scenario to account for the divergence seen here, these 2 historically geographically isolated
species must have independently experienced an ecological
challenge—at about the same time in their evolutionary histories—and responded to that challenge with multiple beneficial substitutions at the same (or related) sets of genes.
Although we cannot rule out such a scenario, an alternative
class of explanation seems more plausible.
In particular, the parallel bouts of adaptive evolution at
nucleoporins in both D. melanogaster and D. simulans suggest that these 2 species may have inherited an unresolved
conflict involving NPCs from their common ancestor. Under this scenario, targeted nucleoporins in these 2 species
would have evolved in response to a conflict using different
sets of substitutions, with interacting nucleoporins coevolving. This adaptive divergence would then incidentally lead
to a hybrid incompatibility. There are at least 3 types of
genetic conflicts in which NPCs—the ‘‘gatekeepers of
the nucleus’’ (Wente 2000)—might plausibly be involved.
actions between flies and viruses (but see Dostert et al. 2005;
Wang et al. 2006), it is easy to imagine that fly pathogens
have evolved similar strategies, leading to conflicts that
drive the divergence of Drosophila nucleoporins.
Genetic Conflict over Centromeric Drive
Meiotic drive of selfish centromeres during oogenesis
represents a second genetic conflict potentially fueling the
divergence of the Nup107 subcomplex. Of the 4 meiotic
products of oogenesis, only 1 becomes the primary oocyte
and the other 3 are discarded as polar bodies. This inherent
asymmetry in oogenesis thus creates an opportunity for centromeres, the sites of kinetochore assembly and spindle attachment, to compete for access to the primary oocyte
(Zwick et al. 1999; Henikoff et al. 2001). There is direct
evidence for such centromeric drive in plants and animals
(Pardo-Manuel de Villena and Sapienza 2001; Fishman and
Willis 2005), as well as some indirect evidence suggested
by the rapid evolution of both centromeric heterochromatin
and centromere-binding proteins in plants and animals, including Drosophila (Malik and Henikoff 2001; Malik et al.
2002; Talbert et al. 2002). Centromeric drive is predicted to
have deleterious side effects for both female and male fertility, including nondisjunction (Zwick et al. 1999) and
other meiotic defects (Henikoff and Malik 2002; Malik
2005). Centromere-binding proteins may therefore evolve
rapidly in order to suppress selfish centromeres. How might
the Nup107 subcomplex be involved in centromeric drive?
One of the surprises to emerge from recent work on the
Nup107 subcomplex is its localization during cell division:
during mitosis the nuclear envelope breaks down and NPCs
disassemble; however, the Nup107 subcomplex remains intact, specifically localizes to the kinetochores of chromosomes, and promotes spindle assembly (Belgareh et al.
2001; Loiodice et al. 2004). The members of the Nup107
subcomplex therefore represent another class of rapidly
evolving, centromere-associated proteins potentially involved in suppressing centromeric drive.
Genetic Conflict between Host and Pathogens
Genetic Conflict over Segregation Distortion
First, there is ample evidence that nucleoporins are involved in antagonistic interactions with viruses and viruslike particles of retrotransposons (Cronshaw and Matunis
2004; Smith and Helenius 2004; Fontoura et al. 2005; Irwin
et al. 2005). Some viruses (e.g., adenovirus and herpes simplex virus) are too large to pass through the nuclear pore and
instead dock on the cytoplasmic face of NPCs to inject their
genomes into the nucleus (Cronshaw and Matunis 2004;
Smith and Helenius 2004; Fontoura et al. 2005). Other
viruses directly target particular nucleoporins for disruption.
In human cells, the matrix protein of the vesicular stomatitus
virus inhibits mRNA export by specifically disrupting
Nup98 (Enninga et al. 2002). In response, host cells upregulate expression of Nup98 and Nup96 via an interferonmediated immune response. Similarly, poliovirus (Gustin
and Sarnow 2001) and rhinovirus (Gustin and Sarnow
2002) impair specific nuclear import pathways by proteolyzing Nup153. Furthermore, HIV-1 and Tf1 retrotransposons
also interact with Nup153 (Varadarajan et al. 2005). Although relatively little is known about the molecular inter-
Segregation distortion represents a third genetic conflict that might involve the NPC. In D. melanogaster,
a male-specific meiotic drive complex called Segregation
Distorter (SD) resides on 1–5% of second chromosomes
in natural populations (Hartl and Hiraizumi 1976; Kusano
et al. 2003). During spermatogenesis in SD/SD1 males,
SD1-bearing sperm are incapacitated so that virtually all
progeny inherit SD instead of the 50% expected under Mendelian inheritance. The 3 major loci of the SD complex are
1) the Segregation distorter locus (Sd), the main driver that
encodes a truncated duplicate of RanGAP (Merrill et al.
1999); 2) Enhancer of SD [E(SD)], an upward modifier
of the strength of distortion; and 3) Responder (Rsp), the
target of distortion. Rsp is an array of AT-rich 120-bp
repeats for which repeat copy number is correlated with sensitivity to distortion—sensitive Rsp alleles comprise 700–
2500 repeats, whereas insensitive alleles comprise a few
dozen to ;200 (Wu et al. 1988). Wildtype RanGAP localizes
to the cytoplasm where it functions as a major regulator of
nuclear transport; in contrast, truncated Sd-RanGAP was
Adaptive Evolution of Nucleoporins 311
recently shown to mislocalize to the nuclear interior, where it
causes distortion by disrupting nuclear transport during spermatogenesis (Kusano et al. 2001). Genes able to influence
nuclear transport—including nucleoporins—could in principle suppress or compensate for the deleterious effects of distortion and might therefore evolve in response to this genetic
conflict. We have argued elsewhere that RanGAP-mediated
distortion may be an ancient genetic conflict, predating the
D. melanogaster–D. simulans split, which has driven rapid
divergence of several nuclear transport genes, including RanGAP (Presgraves DC, unpublished data).
Thus, the rapid molecular evolution of the Nup107
subcomplex and its interactors could plausibly reflect histories of 1) host–pathogen conflict, 2) conflict over centromeric drive, or 3) conflict over nuclear transport–related
segregation distortion. Although admittedly speculative,
such genetic conflict–driven scenarios seem likely to explain the simultaneous bouts of recurrent adaptive substitutions concentrated in the Nup107 subcomplex genes of
D. melanogaster and D. simulans.
Coevolutionary Consequences for Hybrid
Incompatibilities
Population genetic models of the evolution of hybrid
incompatibilities usually assume that substitutions occur independently (Orr 1995). For instance, in the simple 2-locus,
2-allele Dobzhansky–Muller model, an ancestral aabb genotype first becomes AAbb in 1 of 2 lineages; then, with equal
probability, a B substitution can occur either in the same lineage yielding AABB (scenario 1) or in the sister lineage yielding aaBB (scenario 2). The pattern of molecular evolution
among the Nup107 subcomplex genes implies, however,
that the histories of interacting genes are not independent—the Nup107 subcomplex genes appear to coevolve.
Such nonindependence among coevolving interactors has
2 consequences for the evolution of hybrid incompatibilities.
First, coevolution accelerates the rate of evolution of hybrid
incompatibilities relative to a model with independent substitutions because substitutions at 1-locus increase the probability of substitutions at interacting loci—within a lineage
an a/A substitution can drive a b/B substitution.
Second, coevolution decreases the probability of
evolving incompatibilities between 2 derived alleles
(derived–derived) versus incompatibilities between 1 derived and once ancestral allele (derived–ancestral). Under
scenario 1 in the model above, the AABB 3 aabb species
cross will produce AaBb hybrids exposed to a possible
incompatibility between the derived B allele and the ancestral a allele (as ‘‘B’’ has never been tested by natural selection on an ‘‘a’’ genetic background). Similarly, under
scenario 2, the AAbb 3 aaBB species cross will produce
AaBb hybrids exposed to incompatibilities between the
derived A and derived B alleles. Coevolution among interacting loci will cause substitutions to be concentrated
in 1 lineage (scenario 1), increasing the probability of
derived–ancestral incompatibilities relative to a model with
independent substitutions (Orr 1995). Thus, the nonindependence of substitutions in coevolving genes should increase the rate of evolution of derived–ancestral hybrid
incompatibilities.
Coevolution between Nup96 and a Partner on the X
Chromosome?
Genetic analysis shows that the D. simulans allele of
Nup96 causes hybrid inviability via an epistatic interaction
with an incompatible ‘‘partner’’ gene (or genes) on the X
chromosome of D. melanogaster (Presgraves 2003; Presgraves et al. 2003). Of the 8 known members of the
Nup107 subcomplex and the three shuttling nucleoporins
that dock with it, only 1 is encoded by a gene on the X chromosome: the shuttling nucleoporin, Nup153. Interestingly,
Nup153 shows evidence of adaptive evolution only in the
D. simulans lineage, suggesting that the D. simulans allele
of Nup96 is functional in the context of the newly evolved
allele of Nup153 in D. simulans. If true, then this putative hybrid incompatibility could represent a derived–
ancestral incompatibility—the derived allele of Nup96
(from D. simulans) is incompatible with the ancestral allele
of Nup153 (from D. melanogaster). It is important to note,
however, that Nup153 remains only a candidate partner
gene until functional tests demonstrate its role in the
Nup96 incompatibility.
Supplementary Material
Supplementary figures and Accession Numbers
EF057923–EF058045 are available at Molecular Biology
and Evolution online (http://www.mbe.oxfordjournals.org/).
Acknowledgments
We thank Andrea Betancourt, Victoria Cattani, J.P.
Masly, Allen Orr, and Shanwu Tang for comments on
the manuscript. This work was supported by an Alexander
von Humboldt Research Fellowship and funding from the
University of Rochester to D.C.P. and by a grant from the
VW-Foundation (I/78815) to W.S.
Literature Cited
Adams MD, Celniker SE, Holt RA, et al. (195 co-authors). 2000.
The genome sequence of Drosophila melanogaster. Science.
287:2185–2195.
Akashi H. 1995. Inferring weak selection from patterns of polymorphism and divergence at ‘‘silent’’ sites in Drosophila DNA.
Genetics. 139:1067–1076.
Akashi H. 1996. Molecular evolution between Drosophila melanogaster and D. simulans: reduced codon bias, faster rates of
amino acid substitution, and larger proteins in D. melanogaster.
Genetics. 144:1297–1307.
Allen TD, Cronshaw JM, Bagley S, Kiseleva E, Goldberg MW.
2000. The nuclear pore complex: mediator of translocation between nucleus and cytoplasm. J Cell Sci. 113:1651–1659.
Allen NPC, Patel SS, Huang L, Chalkley RJ, Burlingame A,
Lutzmann M, Hurt EC, Rexach M. 2002. Deciphering networks
of protein interactions at the nuclear pore complex. Mol Cell
Proteomics. 1:930–946.
Andolfatto P. 2001. Contrasting patterns of X-linked and autosomal nucleotide variation in Drosophila melanogaster and Drosophila simulans. Mol Biol Evol. 18:279–290.
Aquadro CF, Lado KM, Noon WA. 1988. The rosy region of Drosophila melanogaster and Drosophila simulans. I. Contrasting
levels of naturally occurring DNA restriction map variation
and divergence. Genetics. 119:875–888.
312 Presgraves and Stephan
Bapteste E, Charlebois RL, MacLeod D, Brochier C. 2005. The
two tempos of nuclear pore complex evolution: highly adapting proteins in an ancient frozen structure. Genome Biol.
6:R85.
Barbash DA, Ashburner M. 2003. A novel system of fertility rescue in Drosophila hybrids reveals a link between hybrid lethality and female sterility. Genetics. 163:217–226.
Barbash DA, Awadalla P, Tarone AM. 2004. Functional divergence caused by ancient positive selection of a Drosophila hybrid incompatibility locus. PLoS Biol. 2:839–848.
Barbash DA, Siino DF, Tarone AM, Roote J. Forthcoming 2003.
A rapidly evolving Myb-related protein causes species isolation in Drosophila. Proc Natl Acad Sci USA. 100:5302–
5307.
Baudry E, Viginier B, Veuille M. 2004. Non-African populations
of Drosophila melanogaster have a unique origin. Mol Biol
Evol. 21:1482–1491.
Belgareh N, Rabut G, Bai SW, et al. (12 co-authors). 2001. An
evolutionarily conserved NPC subcomplex, which redistributes in part to kinetochores in mammalian cells. J Cell Biol.
154:1147–1160.
Coyne JA, Orr HA. 1989. Patterns of speciation in Drosophila.
Evolution. 43:362–381.
Coyne JA, Orr HA. 1997. ‘‘Patterns of speciation in Drosophila’’
revisited. Evolution. 51:295–303.
Coyne JA, Orr HA. 2004. Speciation. Sunderland (MA): Sinauer.
Cronshaw JM, Matunis MJ. 2004. The nuclear pore complex: disease associations and functional correlations. Trends Endocrinol Metab. 15:34–39.
Dean MD, Ballard JWO. 2004. Linking phylogenetics with population genetics to reconstruct the geographic origin of a species. Mol Phylogenet Evol. 32:998–1009.
Devos D, Dokudovskaya S, Alber F, Williams R, Chait BT, Sali
A, Rout MP. 2004. Components of coated vesicles and nuclear
pore complexes share a common molecular architecture. PLoS
Biol. 2:2085–2093.
Devos D, Dokudovskaya S, Williams R, Alber F, Eswar N, Chait
BT, Rout MP, Sali A. 2006. Simple fold composition and modular architecture of the nuclear pore complex. Proc Natl Acad
Sci USA. 103:2172–2177.
Dobzhansky T. 1937. Genetics and the origin of species. New
York: Columbia University Press.
Dostert C, Jouanguy E, Irving P, Troxler L, Galiana-Arnoux D,
Hetru C, Hoffman JA, Imler JL. 2005. The Jak-STAT signaling
pathway is required but not sufficient for the antiviral response
of Drosophila. Nat Immunol. 6:946–953.
Enarson P, Enarson M, Bastos R, Burke B. 1998. Amino-terminal
sequences that direct nucleoporin Nup153 to the inner surface
of the nuclear envelope. Chromosoma. 107:228–236.
Enninga J, Levay A, Fontoura BMA. 2003. Sec13 Shuttles between the nucleus and the cytoplasm and stably interacts with
Nup96 at the nuclear pore complex. Mol Biol Cell. 23:
7271–7284.
Enninga J, Levy DE, Blobel G, Fontoura BMA. 2002. Role of
nucleoporin induction in releasing an mRNA nuclear export
block. Science. 295:1523–1525.
Fay JC, Wu C-I. 2000. Hitchhiking under positive Darwinian selection. Genetics. 155:1405–1413.
Fay JC, Wyckoff GJ, Wu C-I. 2002. Testing the neutral theory of
molecular evolution with genomic data from Drosophila. Nature. 415:1024–1026.
Fishman L, Willis JN. 2005. A novel meiotic drive locus almost
completely distorts segregation in mimulus (monkeyflower)
hybrids. Genetics. 169:347–353.
Fontoura BMA, Faria PA, Nussenzveig DR. 2005. Viral interactions with the nuclear transport machinery: discovering and
disrupting pathways. Life. 57:65–72.
Glinka S, Ometto L, Mousset S, Stephan W, DeLorenzo D. 2003.
Demography and natural selection have shaped genetic variation in Drosophila melanogaster: a multi-locus approach. Genetics. 165:1269–1278.
Griffis ER, Altan N, Lippincott-Schwartz J, Powers MA. 2002.
Nup98 is a mobile nucleoporin with transcription dependent
dynamics. Mol Biol Cell. 13:1282–1297.
Gustin KE, Sarnow P. 2001. Effects of poliovirus infection on
nucleo-cytoplasmic trafficking and nuclear pore complex composition. EMBO J. 20:240–249.
Gustin KE, Sarnow P. 2002. Inhibition of nuclear import and alteration of nuclear pore complex composition by rhinovirus.
J Virol. 76:8787–8796.
Hartl DL, Hiraizumi Y. 1976. Segregation distortion. In: Novitski
E, Ashburner M, editors. The genetics and biology of Drosophila melanogaster. New York: Academic Press.
Henikoff S, Ahmad K, Malik HS. 2001. The centromere paradox:
stable inheritance with rapidly evolving DNA. Science.
293:1098–1102.
Henikoff S, Malik HS. 2002. Centromeres: selfish drivers. Nature.
417:227.
Hey J, Kliman RM. 1993. Population genetics and phylogenetics
of DNA sequence variation at multiple loci within the Drosophila melanogaster species complex. Mol Biol Evol.
10:804–822.
Hodel AE, Hodel MR, Griffis ER, Hennig KA, Ratner GA, Xu S,
Powers MA. 2002. The three-dimensional structure of the
auteoproteolytic, nuclear pore-targeting domain of the human
nucleoporin Nup98. Mol Cell. 10:347–358.
Hudson RR. 1987. Estimating the recombination parameter of a
finite population-model without selection. Genet Res. 50:
245–250.
Irwin B, Aye M, Baldi P, Beliakova-Bethell N, Cheng H, Dou Y,
Liou W, Sandmeyer S. 2005. Retroviruses and yeast retrotransposons use overlapping sets of host genes. Genome Res.
15:641–654.
Kaplan NL, Hudson RR, Langley CH. 1989. The ‘‘hitch-hiking
effect’’ revisited. Genetics. 123:887–899.
Kimura M. 1983. The neutral theory of molecular evolution. Cambridge: Cambridge University Press.
Kliman RM, Andolfatto P, Coyne JA, Depaulis F, Kreitman M,
Berry AJ, McCarter J, Wakeley J, Hey J. 2000. The population
genetics of the origin and divergence of the Drosophila simulans complex species. Genetics. 156:1913–1931.
Kusano A, Staber C, Chan HYE, Ganetzky B. 2003. Closing the
(Ran)GAP on segregation distortion in Drosophila. BioEssays.
25:108–115.
Kusano A, Staber, Ganetzky B. 2001. Nuclear mislocalization of
enzymatically active RanGAP causes segregation distortion in
Drosophila. Dev Cell. 1:351–361.
Lachaise D, Capy P, Cariou M-L, Joly D, Lemeunier F, David JR.
2004. Nine relatives from one African ancestor: population biology and evolution of the Drosophila melanogaster subgroup
species. In: Singh RS, Uyenoyama M, editors. The evolution of
population biology. New York: Cambridge University Press.
Lachaise D, Cariou M-L, David JR, Lemeunier F, Tsacas L,
Ashburner M. 1988. Historical biogeography of the Drosophila
melanogaster species subgroup. Evol Biol. 22:159–225.
Lachaise D, Silvain J-F. 2004. How two Afrotropical endemics
made two cosmopolitan human commensals: the Drosophila
melanogaster-D. simulans palaeogeographic riddle. Genetica.
120:17–39.
Landry CR, Wittkopp PJ, Taubes CH, Ranz JM, Clark AG, Hartl
DL. 2005. Compensatory cis-trans evolution and the dysregulation of gene expression in interspecific hybrids of Drosophila. Genetics. 171:1813–1822.
Li W-H. 1997. Molecular evolution. Sunderland (MA): Sinauer.
Adaptive Evolution of Nucleoporins 313
Loiodice I, Alves A, Rabut G, Van Overbeek M, Ellenberg J,
Sibarita JB, Doye V. 2004. The entire Nup107-160 complex,
including three new members, is targeted as one entity to kinetochores in mitosis. Mol Cell Biol. 15:3333–3344.
Lutzmann M, Kunze R, Buerer A, Aebi U, Hurt E. 2002. Modular
self-assembly of a Y-shaped multiprotein complex from seven
nucleoporins. EMBO J. 21:387–397.
Lynch M, Force AG. 2000. The origin of interspecific genomic
incompatibility via gene duplication. Am Nat. 156:
590–605.
Malik HS. 2005. Mimulus finds centromeres in the driver’s seat.
Trends Ecol Evol. 20:151–154.
Malik HS, Henikoff S. 2001. Adaptive evolution of Cid, a
centromere-specific histone in Drosophila. Genetics. 157:
1293–1298.
Malik HS, Vermaak D, Henikoff S. 2002. Recurrent evolution of
DNA-binding motifs in the Drosophila centromeric histone.
Proc Natl Acad Sci USA. 99:1449–1454.
Masly JP, Jones CD, Noor MAF, Locke J, Orr HA. 2006. Gene
transposition as a novel cause of hybrid male sterility. Science.
313:1448–1450.
Maynard Smith J, Haigh J. 1974. The hitch-hiking effect of
a favourable gene. Genet Res. 23:23–35.
McDonald JH, Kreitman M. 1991. Adaptive protein evolution at
the Adh locus in Drosophila. Nature. 351:652–654.
Merrill C, Bayraktaroglu L, Kusano A, Ganetzky B. 1999. Truncated RanGAP encoded by the Segregation Distorter locus of
Drosophila. Science. 283:1742–1745.
Mousset S, Derome N. 2004. Molecular polymorphism in Drosophila melanogaster and D. simulans: what have we learned
from recent studies? Genetica. 120:79–86.
Muller HJ. 1940. Bearing of the Drosophila work on systematics.
In: Huxley JS, editor. The new systematics. Oxford: Clarendon
Press. p. 185–268.
Muller HJ. 1942. Isolating mechanisms, evolution, and temperature. Biol Symp. 6:71–125.
Nielsen R. 2005. Molecular signatures of natural selection. Annu
Rev Genet. 39:197–218.
Ometto L, Glinka S, DeLorenzo D, Stephan W. 2005. Inferring the
effects of demography and selection on Drosophila melanogaster populations from a chromosome-wide scan of
DNA variation. Mol Biol Evol. 22:2119–2130.
Orr HA. 1995. The population genetics of speciation: the evolution of hybrid incompatibilities. Genetics. 139:1805–1813.
Orr HA. 2005. The genetic basis of reproductive isolation: insights
from Drosophila. Proc Natl Acad Sci USA. 102:6522–6526.
Orr HA, Irving S. 2005. Segregation distortion in hybrids between
the Bogota and USA subspecies of Drosophila pseudoobscura.
Genetics. 169:671–682.
Orr HA, Masly JP, Presgraves DC. 2004. Speciation genes. Curr
Opin Genet Dev. 14:675–679.
Pardo-Manuel de Villena F, Sapienza C. 2001. Female meiosis
drive karyotypic evolution in mammals. Genetics. 159:
1179–1189.
Presgraves DC. 2003. A fine-scale genetic analysis of hybrid incompatibilities in Drosophila. Genetics. 163:955–972.
Presgraves DC. 2005. Recombination enhances protein adaptation
in Drosophila melanogaster. Curr Biol. 15:1651–1656.
Presgraves DC, Balagopalan L, Abmayr SM, Orr HA. 2003.
Adaptive evolution drives divergence of a hybrid inviability
gene between two species of Drosophila. Nature. 423:
715–719.
Rambaut A. 1996. Se-Al: sequence alignment editor. Available at
http://evolve.zoo.ox.ac.uk.
Rappsilber J, Siniossoglou S, Hurt EC, Mann M. 2000. A generic
strategy to analyze the spatial organization of multi-protein
complexes by cross-linking and mass spectrometry. Anal
Chem. 72:267–275.
Rout MP, Aitchison JD, Suprapto A, Hjertaas K, Zhao Y, Chait
BT. 2000. The yeast nuclear pore complex: composition,
architecture and transport mechanism. J Cell Biol. 148:
635–651.
Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R. 2003.
DnaSP, DNA polymorphism analyses by the coalescent and
other methods. Bioinformatics. 19:2496–2497.
Sawyer SA, Kulathinal RJ, Bustamante CD, Hartl DL. 2003.
Bayesian analysis suggests that most amino acid replacements
in Drosophila are driven by positive selection. J Mol Evol.
57:S154–S164.
Siniossoglou S, Lutzmann M, Santos-Rosa H, Leonard K, Mueller
S, Aebi U, Hurt EC. 2000. Structure and assembly of the
Nup84p complex. J Cell Biol. 149:41–53.
Siniossoglou S, Wimmer C, Rieger M, Doye V, Tekotte H, Weise
C, Emig S, Segref A, Hurt EC. 1996. A novel complex of
nucleoporins, which includes Sec13p and a Sec13p homolog,
is essential for normal nuclear pores. Cell. 84:265–275.
Smith N, Eyre-Walker A. 2002. Adaptive protein evolution in
Drosophila. Nature. 415:1022–1024.
Smith AE, Helenius A. 2004. How viruses enter animal cells. Science. 304:237–242.
Tajima F. 1989. Statistical method for testing the neutral mutation
hypothesis by DNA polymorphism. Genetics. 123:585–595.
Talbert PB, Masuelli R, Tyagi AP, Comai L, Henikoff S. 2002.
Centromeric localization and adaptive evolution of an Arabidopsis histone H3 variant. Plant Cell. 14:1053–1066.
Tamura K, Subramanian S, Kumar S. 2004. Temporal patterns of
fruitfly (Drosophila) evolution revealed by mutation clocks.
Mol Biol Evol. 21:36–44.
Thornton K, Andolfatto P. 2006. Approximate Bayesian inference
reveals evidence for a recent, severe, bottleneck in non-African
populations of Drosophila melanogaster. Genetics. 172:1607–
1619.
Ting C-T, Tsauer SC, Sun S, Brown WE, Chen YC, Patel NH, Wu
C-I. 2004. Gene duplication and speciation in Drosophila: evidence from the Odysseus locus. Proc Natl Acad Sci.
101:12232–12235.
Ting C-T, Tsaur S-C, Wu M-L, Wu C-I. 1998. A rapidly evolving
homeobox at the site of a hybrid sterility gene. Science.
282:1501–1504.
Tran EJ, Wente SR. 2006. Dynamic nuclear pore complexes: life
on the edge. Cell. 125:1041–1053.
Varadarajan P, Mahalingam S, Liu P, Ng SBH, Gandotra S,
Dorairajoo DSK, Balasundaram D. 2005. The functionally
conserved nucleoporins Nup124p from fission yeast and the
human Nup153 mediate nuclear import and activity of the
Tf1 retrotransposon and HIV-1 vpr. Mol Biol Cell.
16:1823–1838.
Vasu SK, Forbes DJ. 2001. Nuclear pores and nuclear assembly.
Curr Opin Cell Biol. 13:363–375.
Vasu SK, Shah S, Orjalo A, Park M, Fischer WH, Forbes DJ.
2001. Novel vertebrate nucleoporins Nup133 and Nup160 play
a role in mRNA export. J Cell Biol. 155:339–353.
Walther TC, Askjaer P, Gentzel M, Haberman A, Griffiths G,
Wilm M, Mattaj IW, Hetzer M. 2003. RanGTP mediates
nuclear pore complex assembly. Nature. 424:689–694.
Wang XH, Aliyari R, Li WX, Li WH, Kim K, Carthew R,
Atkinson P, Ding SW. 2006. RNA interference directs innate
immunity against viruses in adult Drosophila. Science. 312:
452–454.
Watterson GA. 1975. On the number of segregating sites in genetical models without recombination. Theor Popul Biol.
7:256–276.
314 Presgraves and Stephan
Wente SR. 2000. Gatekeepers of the nucleus. Science. 288:1374–
1377.
Wu C-I, Lyttle TW, Wu M-L, Lin GF. 1988. Association between
DNA satellite sequences and the responder of segregation distortion in D. melanogaster. Cell. 54:179–189.
Wu C-I, Ting C-T. 2004. Genes and speciation. Nat Rev Genet.
5:114–122.
Zwick ME, Salstrom JL, Langley CH. 1999. Genetic variation in
rates of nondisjunction: association of two naturally occurring
polymorphisms in the chromokinesin nod with increased rates
of nondisjunction in Drosophila melanogaster. Genetics.
152:1605–1614.
Jianzhi Zhang, Associate Editor
Accepted October 17, 2006