Simple Y-Autosomal Incompatibilities Cause Hybrid Male Sterility in Reciprocal Crosses Between Drosophila virilis and D. americana

Postzygotic reproductive isolation evolves when hybrid incompatibilities accumulate between diverging populations. Here, I examine the genetic basis of hybrid male sterility between two species of Drosophila, Drosophila virilis and D. americana. From these analyses, I reach several conclusions. First, neither species carries any autosomal dominant hybrid male sterility alleles: reciprocal F₁ hybrid males are perfectly fertile. Second, later generation (backcross and F₂) hybrid male sterility between D. virilis and D. americana is not polygenic. In fact, I identified only three genetically independent incompatibilities that cause hybrid male sterility. Remarkably, each of these incompatibilities involves the Y chromosome. In one direction of the cross, the D. americana Y is incompatible with recessive D. virilis alleles at loci on chromosomes 2 and 5. In the other direction, the D. virilis Y chromosome causes hybrid male sterility in combination with recessive D. americana alleles at a single QTL on chromosome 5. Finally, in contrast with findings from other Drosophila species pairs, the X chromosome has only a modest effect on hybrid male sterility between D. virilis and D. americana.SPECIATION occurs when populations evolve one or more barriers to interbreeding (Dobzhansky 1937; Mayr 1963). One such barrier is intrinsic postzygotic isolation, which typically evolves when diverging populations accumulate different alleles at two or more loci that are incompatible when brought together in hybrid genomes; negative epistasis between these alleles renders hybrids inviable or sterile (Bateson 1909; Dobzhansky 1937; Muller 1942). Classical and recent studies in diverse animal taxa have provided support for two evolutionary patterns that often characterize the genetics of postzygotic isolation (Coyne and Orr 1989a). The first, Haldane''s rule, observes that when there is F₁ hybrid inviability or sterility that affects only one sex, it is almost always the heterogametic sex (Haldane 1922). Over the years, many researchers have tried to account for this pattern, but only two ideas are now thought to provide a general explanation: the “dominance theory,” which posits that incompatibility alleles are generally recessive in hybrids, and the “faster-male theory,” which posits that genes causing hybrid male sterility diverge more rapidly than those causing hybrid female sterility (Muller 1942; Wu and Davis 1993; Turelli and Orr 1995; reviewed in Coyne and Orr 2004). In some cases, however, additional factors might contribute to Haldane''s rule, including meiotic drive, a faster-evolving X chromosome, dosage compensation, and Y chromosome incompatibilities (reviewed in Laurie 1997; Turelli and Orr 2000; Coyne and Orr 2004).The second broad pattern affecting the evolution of postzygotic isolation is the disproportionately large effect of the X chromosome on heterogametic F₁ hybrid sterility (Coyne 1992). This “large X effect” has been documented in genetic analyses of backcross hybrid sterility (e.g., Dobzhansky 1936; Grula and Taylor 1980; Orr 1987; Masly and Presgraves 2007) and inferred from patterns of introgression across natural hybrid zones (e.g., Machado et al. 2002; Saetre et al. 2003; Payseur et al. 2004). However, in only one case has the cause of the large X effect been unambiguously determined: incompatibilities causing hybrid male sterility between Drosophila mauritiana and D. sechellia occur at a higher density on the X than on the autosomes (Masly and Presgraves 2007). Testing the generality of this pattern will require additional high-resolution genetic analyses in diverse taxa (Presgraves 2008). But whatever its causes, there is now general consensus that the X chromosome often plays a special role in the evolution of postzygotic isolation (Coyne and Orr 2004).The contribution of the Y chromosome to animal speciation is less clear. Y chromosomes have far fewer genes than the X or autosomes, and most of these genes are male specific (Lahn and Page 1997; Carvalho et al. 2009). In Drosophila species, the Y chromosome is typically required for male fertility, but not for viability (Voelker and Kojima 1971). How often, then, does the Y chromosome play a role in reproductive isolation? In crosses between Drosophila species, hybrid male sterility is frequently caused by incompatibilities between the X and Y chromosomes (Schafer 1978; Heikkinen and Lumme 1998; Mishra and Singh 2007) or between the Y and heterospecific autosomal alleles (Patterson and Stone 1952; Vigneault and Zouros 1986; Lamnissou et al. 1996). In crosses between D. yakuba and D. santomea, the Y chromosome causes F₁ hybrid male sterility, and accordingly, shows no evidence for recent introgression across a species hybrid zone (Coyne et al. 2004; Llopart et al. 2005). In mammals, reduced introgression of Y-linked loci (relative to autosomal loci) has been shown across natural hybrid zones of mice (Tucker et al. 1992) and rabbits (Geraldes et al. 2008), suggesting that the Y chromosome contributes to reproductive barriers.Here I examine the genetic basis of hybrid male sterility between two species of Drosophila, D. virilis and D. americana. These species show considerable genetic divergence (K_s ∼0.11, Morales-Hojas et al. 2008) and are currently allopatric: D. virilis is a human commensal worldwide with natural populations in Asia, and D. americana is found in riparian habitats throughout much of North America (Throckmorton 1982; McAllister 2002). Nearly 70 years ago, Patterson et al. (1942) showed that incompatibilities between the D. americana Y chromosome and the second and fifth chromosomes from D. virilis cause hybrid male sterility, a result that was confirmed in a more recent study (Lamnissou et al. 1996). Another study suggested that the X chromosome might play the predominant role in causing hybrid male sterility between D. virilis and D. americana (Orr and Coyne 1989). But because previous genetic analyses had to rely on only a few visible markers to map hybrid male sterility, they lacked the resolution to examine the genomic distribution of incompatibility loci.Using the D. virilis genome sequence, I have developed a dense set of molecular markers to investigate the genetic architecture of hybrid male sterility between D. virilis and D. americana. In this study, I perform a comprehensive set of crosses to address several key questions: What is the effect of the X chromosome on hybrid male sterility between D. virilis and D. americana? What is the effect of the Y chromosome? Approximately how many loci contribute to hybrid male sterility between these Drosophila species? Perhaps surprisingly, the answers to these questions differ dramatically from what has been found for other Drosophila species, including the well-studied D. melanogaster group.     

Alternative Splicing Modulates Ubx Protein Function in Drosophila melanogaster

Polyploidy is an important aspect of the evolution of flowering plants. The potential of gene copies to diverge and evolve new functions is influenced by meiotic behavior of chromosomes leading to segregation as a single locus or duplicated loci. Switchgrass (Panicum virgatum) linkage maps were constructed using a full-sib population of 238 plants and SSR and STS markers to access the degree of preferential pairing and the structure of the tetraploid genome and as a step toward identification of loci underlying biomass feedstock quality and yield. The male and female framework map lengths were 1645 and 1376 cM with 97% of the genome estimated to be within 10 cM of a mapped marker in both maps. Each map coalesced into 18 linkage groups arranged into nine homeologous pairs. Comparative analysis of each homology group to the diploid sorghum genome identified clear syntenic relationships and collinear tracts. The number of markers with PCR amplicons that mapped across subgenomes was significantly fewer than expected, suggesting substantial subgenome divergence, while both the ratio of coupling to repulsion phase linkages and pattern of marker segregation indicated complete or near complete disomic inheritance. The proportion of transmission ratio distorted markers was relatively low, but the male map was more extensively affected by distorted transmission ratios and multilocus interactions, associated with spurious linkages.POLYPLOIDY is common among plants (Masterson 1994; Levin 2002) and is an important aspect of plant evolution. Widespread paleopolyploidy in flowering plant lineages suggests that ancient polyploidization events have contributed to the radiation of angiosperms (Soltis et al. 2009; Van de Peer et al. 2009a). Whole genome duplications are thought to be the sources of evolutionary novelty (Osborn et al. 2003; Freeling and Thomas 2006; Chen 2007; Hegarty and Hiscock 2008; Flagel and Wendel 2009; Leitch and Leitch 2008). Other attributes of polyploids considered to promote evolutionary success include increased vigor, masking of recessive alleles, and reproductive barriers arising from loss of one of the duplicate genes (Soltis and Soltis 2000; Comai 2005; Otto 2007; Van de Peer et al. 2009b). Among crop species, polyploidy likely contributed to trait improvement under artificial selection (Paterson 2005; Udall and Wendell 2006; Dubcovsky and Dvorak 2007; Hovav et al. 2008).Disomic inheritance in polyploids, in contrast to polysomic inheritance, presents opportunities for duplicated genes to diverge and evolve new functions. The relative age of whole genome duplications and the extent of homology between subgenomes greatly influence chromosomal pairing at meiosis (Soltis and Soltis 1995; Wolfe 2001; Ramsey and Schemske 2002). Polysomic inheritance resulting from random chromosome pairing is associated with doubling of a single set of chromosomes. Disomic inheritance resulting from preferential pairing is often associated with polyploidy arising from combinations of divergent genomes. The evolutionary process of diploidization leads to a shift from random to preferential pairing that is not well understood but is genetically defined in systems such as Ph1 of wheat (Triticum aestivum) and PrBn of Brassica napus (Riley and Chapman 1958; Vega and Feldman 1998; Jenczewski et al. 2003). The degree of preferential pairing also affects allelic diversity and the ability to detect linkage. Accurate information about chromosome pairing and whole or partial genome duplications is thus important for both evolutionary studies and in linkage analysis.Such information is extremely limited in the C4 panicoid species Panicum virgatum (switchgrass), which is now viewed as a promising energy crop in the United States and Europe (Lewandowski et al. 2003; McLaughlin and Kszos 2005) and is planted extensively for forage and soil conservation (Vogel and Jung 2001). Little is known about either its genome structure or inheritance. Much current bioenergy feedstock development is focused on tetraploid cytotypes (2n = 4x = 36) due to their higher yield potentials, and an initial segregation study indicated a high degree of preferential pairing in a single F1 mapping population (Missaoui et al. 2005). A once-dominant component of the tallgrass prairie in North America, switchgrass is largely self-incompatible (Martinez-Reyna and Vogel 2002) with predominantly tetraploid or octoploid cytotypes (Hultquist et al. 1997; Lu et al. 1998). Limited gene flow appears possible between different cytotypes suggested by DNA content variation within collection sites and seed lots (Nielsen 1944; Hultquist et al. 1997; Narasimhamoorthy et al. 2008). True diploids appear to be rare (Nielsen 1944; Young et al. 2010). Multivalents in meiosis have not been observed in tetraploids or F1 hybrids between upland and lowland tetraploids, although rare univalents occurred (Barnett and Carver 1967; Martinez-Reyna et al. 2001). However, polysomic inheritance may occur with random bivalent pairing (Howard and Swaminathan 1953).Sustainable production of switchgrass for bioenergy to meet the goal of reducing greenhouse gas emissions will require advances in feedstock production that include improvements in yield (Carroll and Somerville 2009). Switchgrass has extensive genetic diversity and potential for genetic improvements, but each cycle of phenotypic selection can take several years (McLaughlin and Kszos 2005; Parrish and Fike 2005; Bouton 2007). Detailed understanding of genome structure to enable efficient marker-assisted selection (MAS) can speed this process considerably. Complete linkage maps are therefore required to both understand chromosome pairing and allow MAS.We report the construction of the first complete linkage maps of two switchgrass genotypes. The linkage maps provide genetic evidence for disomic inheritance in lowland, tetraploid switchgrass. Gene-derived markers enabled a comparative analysis to sorghum, revealing syntenic relationships between the diploid sorghum genome and the tetraploid switchgrass subgenomes. Transmission ratio distortion and multilocus interactions were analyzed in detail to document their potential influence on map accuracy and map-based studies in switchgrass.     

Unraveling the Complex Trait of Crop Yield With Quantitative Trait Loci Mapping in Brassica napus

Yield is the most important and complex trait for the genetic improvement of crops. Although much research into the genetic basis of yield and yield-associated traits has been reported, in each such experiment the genetic architecture and determinants of yield have remained ambiguous. One of the most intractable problems is the interaction between genes and the environment. We identified 85 quantitative trait loci (QTL) for seed yield along with 785 QTL for eight yield-associated traits, from 10 natural environments and two related populations of rapeseed. A trait-by-trait meta-analysis revealed 401 consensus QTL, of which 82.5% were clustered and integrated into 111 pleiotropic unique QTL by meta-analysis, 47 of which were relevant for seed yield. The complexity of the genetic architecture of yield was demonstrated, illustrating the pleiotropy, synthesis, variability, and plasticity of yield QTL. The idea of estimating indicator QTL for yield QTL and identifying potential candidate genes for yield provides an advance in methodology for complex traits.YIELD is the most important and complex trait in crops. It reflects the interaction of the environment with all growth and development processes that occur throughout the life cycle (Quarrie et al. 2006). Crop yield is directly and multiply determined by yield-component traits (such as seed weight and seed number). Yield-related traits (such as biomass, harvest index, plant architecture, adaptation, resistance to biotic and abiotic constraints) may also indirectly affect yield by affecting the yield-component traits or by other, unknown mechanisms. Increasing evidence suggests that “fine-mapped” quantitative trait loci (QTL) or genes identified as affecting crop yield involve diverse pathways, such as seed number (Ashikari et al. 2005; Tian et al. 2006b; Burstin et al. 2007; Xie et al. 2008; Xing et al. 2008; Xue et al. 2008), seed weight (Ishimaru 2003; Song et al. 2005; Shomura et al. 2008; Wang et al. 2008; Xie et al. 2006, 2008; Xing et al. 2008; Xue et al. 2008), flowering time (Cockram et al. 2007; Song et al. 2007; Xie et al. 2008; Xue et al. 2008), plant height (Salamini 2003; Ashikari et al. 2005; Xie et al. 2008; Xue et al. 2008), branching (Clark et al. 2006; Burstin et al. 2007; Xing et al. 2008), biomass yield (Quarrie et al. 2006; Burstin et al. 2007), resistance and tolerance to biotic and abiotic stresses (Khush 2001; Brown 2002; Yuan et al. 2002; Waller et al. 2005; Zhang 2007; Warrington et al. 2008), and root architecture (Hochholdinger et al. 2008).Many experiments have explored the genetic basis of yield and yield-associated traits (yield components and yield-related traits) in crops. Summaries of identified QTL have been published for wheat (MacCaferri et al. 2008), barley (Von Korff et al. 2008), rice, and maize (http://www.gramene.org/). The results show several common patterns. First, QTL for yield and yield-associated traits tend to be clustered in the genome, which suggests that the QTL of the yield-associated traits have pleiotropic effects on yield. Second, this kind of pleiotropy has not been well analyzed genetically. The QTL for yield (complicated factor), therefore, have not been associated with any yield-associated traits (relatively simple factors, such as plant height). Therefore, they are unlikely to predict accurately potential candidate genes for yield. Third, only a few loci (rarely >10) have been found for each of these traits. Thus, the genetic architecture of yield has remained ambiguous. Fourth, trials were carried out in a few environments and how the mode of expression of QTL for these complex traits might respond in different environments is unclear.In this study, the genetic architecture of crop yield was analyzed through the QTL mapping of seed yield and eight yield-associated traits in two related populations of rapeseed (Brassica napus) that were grown in 10 natural environments. The complexity of the genetic architecture of seed yield was demonstrated by QTL meta-analysis. The idea of estimating indicator QTL (QTL of yield-associated traits, which are defined as the potential genetic determinants of the colocalized QTL for yield) for yield QTL in conjunction with the identification of candidate genes is described.     

Genome Changes After Gene Duplication: Haploidy vs. Diploidy

Since genome size and the number of duplicate genes observed in genomes increase from haploid to diploid organisms, diploidy might provide more evolutionary probabilities through gene duplication. It is still unclear how diploidy promotes genomic evolution in detail. In this study, we explored the evolution of segmental gene duplication in haploid and diploid populations by analytical and simulation approaches. Results show that (1) under the double null recessive (DNR) selective model, given the same recombination rate, the evolutionary trajectories and consequences are very similar between the same-size gene-pool haploid vs. diploid populations; (2) recombination enlarges the probability of preservation of duplicate genes in either haploid or diploid large populations, and haplo-insufficiency reinforces this effect; and (3) the loss of duplicate genes at the ancestor locus is limited under recombination while under complete linkage the loss of duplicate genes is always random at the ancestor and newly duplicated loci. Therefore, we propose a model to explain the advantage of diploidy: diploidy might facilitate the increase of recombination rate, especially under sexual reproduction; more duplicate genes are preserved under more recombination by originalization (by which duplicate genes are preserved intact at a special quasi-mutation-selection balance under the DNR or haplo-insufficient selective model), so genome sizes and the number of duplicate genes in diploid organisms become larger. Additionally, it is suggested that small genomic rearrangements due to the random loss of duplicate genes might be limited under recombination.USUALLY genome size becomes larger from haploid to diploid organisms (Lynch and Conery 2003), and so does the number of duplicate genes observed in genomes (Zhang 2003). It is extensively hypothesized that diploidy might facilitate the preservation and accumulation of duplicate genes, but it is still unclear how diploidy supports the evolution of duplicate genes in detail. The superiority of diploidy is classically attributed to preventing expression of deleterious mutations (Crow and Kimura 1965), but it is also argued that the sheltering of deleterious mutations cannot adequately explain the advantages of diploidy (Perrot et al. 1991).Recombination is a common phenomenon in all three kingdoms of life, Bacteria, Eukarya, and Archaea. It has been reported that recombination influences the loss of duplicate genes (Zhang and Kishino 2004; Xue et al. 2010). In diploid organisms, if recombination between the ancestor locus and the newly duplicated locus is free, the rate of recombination is maximally 0.5, which is commonly observed especially when the two loci are located on different chromosomes. Although recombination should not be regarded as an exception in haploid organisms (Fraser et al. 2007), recombination events usually occur more frequently in diploid populations than they do in haploid populations. In other words, diploidy might facilitate the occurrence of recombination. The difference of recombination behaviors between haploid and diploid organisms is an obvious and important feature during genomic evolution.In our recent studies of genomic duplication, we proposed a new possible way of preserving and accumulating duplicate genes in genomes—originalization (Xue and Fu 2009a). As is well known, for a locus in an infinite diploid population, the frequencies of wild-type and degenerative alleles will move to an equilibrium under purifying selection and mutation, which is known as the mutation–selection balance. After genomic duplication, under two simple selective models, double null recessive (DNR, under which valid individuals require at least one active wild-type allele on the ancestor and newly duplicated loci) and haplo-insufficient (HI or partial dominant, under which valid individuals require at least two active wild-type alleles on both loci) models, a special equilibrium of allele frequencies at the ancestor and newly duplicated loci will be reached under recombination, in which the frequency of wild-type allele is kept high at both loci. Under the HI selective model this balance becomes so stable and flexible that the fixation of a degenerative allele at one of these two loci (or the balance being broken) becomes very difficult even in a modest population (Xue and Fu 2009a,b). However, if the two loci are tightly linked (recombination rate r = 0), this balance of allele frequencies does not appear. As r increases, the balance becomes more stable and the frequency of the wild-type allele at two loci becomes higher. High frequency of the wild-type allele at both loci means that duplicate genes are preserved intact in genomes, so this phenomenon was named originalization.Although many duplicate genes originated from genomic duplications in some species, such as yeast, maize, and fish (Li et al. 2005), those from segmental duplications are also very popular (Zhang et al. 2000; Leister 2004). In haploid populations, most duplication events are small segmental duplications. Therefore, to understand genomic evolution comprehensively, it is necessary to explore the evolution of segmental genomic duplication.Lynch et al. (2001) and Tanaka et al. (2009) have studied the evolution of segmental gene duplication in diploid populations theoretically. However, in this study, we further compared the evolution of segmental gene duplication in haploid vs. diploid populations by numerical and simulation approaches under the DNR and HI selective models. We observed that haploid and diploid populations with the same-size gene pool are very similar under the DNR model and the same recombination rate. Recombination enlarges the probability of preservation of duplicate genes in either haploid or diploid populations via originalization, and haplo-insufficiency reinforces this effect. The loss of duplicate genes at the ancestor locus might be limited under recombination, while under complete linkage, the loss of duplicate genes is random at the ancestor and newly duplicated loci. According to these results, we propose a model with which to explain the revolutionary genomic transition from haploidy to diploidy.     

Estimating the Parameters of Selection on Nonsynonymous Mutations in Drosophila pseudoobscura and D. miranda

We present the results of surveys of diversity in sets of >40 X-linked and autosomal loci in samples from natural populations of Drosophila miranda and D. pseudoobscura, together with their sequence divergence from D. affinis. Mean silent site diversity in D. miranda is approximately one-quarter of that in D. pseudoobscura; mean X-linked silent diversity is about three-quarters of that for the autosomes in both species. Estimates of the distribution of selection coefficients against heterozygous, deleterious nonsynonymous mutations from two different methods suggest a wide distribution, with coefficients of variation greater than one, and with the average segregating amino acid mutation being subject to only very weak selection. Only a small fraction of new amino acid mutations behave as effectively neutral, however. A large fraction of amino acid differences between D. pseudoobscura and D. affinis appear to have been fixed by positive natural selection, using three different methods of estimation; estimates between D. miranda and D. affinis are more equivocal. Sources of bias in the estimates, especially those arising from selection on synonymous mutations and from the choice of genes, are discussed and corrections for these applied. Overall, the results show that both purifying selection and positive selection on nonsynonymous mutations are pervasive.SURVEYS of DNA sequence diversity and divergence are shedding light on a number of questions in evolutionary genetics (for recent reviews, see Akey 2009; Sella et al. 2009). Two of the most important questions of this kind concern the distribution of selection coefficients against deleterious mutations affecting protein sequences and the proportion of amino acid sequence differences between related species that have been fixed by positive selection. Several different methods have been proposed for studying each of these questions, using different features of data on polymorphism and divergence at nonsynonymous and silent sites.For example, the parameters of the distribution of selection coefficients against deleterious amino acid mutations have been estimated by contrasting the numbers of nonsynonymous and silent within-species polymorphisms and fixed differences between species (Sawyer and Hartl 1992; Bustamante et al. 2002; Piganeau and Eyre-Walker 2003; Sawyer et al. 2007); by fitting the frequency spectra of nonsynonymous and silent variants to models of selection, mutation, and drift (Akashi 1999; Eyre-Walker et al. 2006; Keightley and Eyre-Walker 2007; Kryukov et al. 2007; Boyko et al. 2008; Eyre-Walker and Keightley 2009); or by comparing levels of nonsynonymous and silent diversities between species with different population sizes (Loewe and Charlesworth 2006; Loewe et al. 2006). The results of these different approaches generally agree in suggesting that there is a wide distribution of selection coefficients against nonsynonymous mutations and that the mean selection coefficient against heterozygous carriers of such mutations is very small. The results imply that a typical individual from a human population carries several hundred weakly deleterious mutations (Eyre-Walker et al. 2006; Kryukov et al. 2007; Boyko et al. 2008); for a typical Drosophila population, with its much higher level of variability, the number is probably an order of magnitude greater (Loewe et al. 2006; Keightley and Eyre-Walker 2007).The presence of this large load of slightly deleterious mutations in human and natural populations, most of which are held at low frequencies by natural selection, has many implications. From the point of view of understanding human genetic disease, it means that we have to face the likelihood that susceptibility to a disease can be influenced by variants at many loci, each with small effects (Kryukov et al. 2007). The pervasive presence of deleterious mutations throughout the genome contributes to inbreeding depression (Charlesworth and Willis 2009) and may mean that the effective population size is reduced by background selection effects, even in regions of the genome with normal levels of genetic recombination (Loewe and Charlesworth 2007). Their presence may contribute so strongly to Hill–Robertson effects (Hill and Robertson 1966; Felsenstein 1974) that they cause severely reduced levels of diversity and adaptation in low-recombination regions of the genome (Charlesworth et al. 2010) and create a selective advantage to maintaining nonzero levels of recombination (Keightley and Otto 2006; Charlesworth et al. 2010). In addition, having an estimate of the distribution of selection coefficients against deleterious nonsynonymous mutations allows their contribution to between-species divergence to be predicted, providing a way of estimating the fraction of fixed nonsynonymous differences caused by positive selection (Loewe et al. 2006; Boyko et al. 2008; Eyre-Walker and Keightley 2009).It is thus important to collect data that shed light on the properties of selection against nonsynonymous mutations in a wide range of systems and also to compare the results from different methods of estimation, since they are subject to different sources of difficulty and biases. In a previous study, we proposed the use of a comparison between two related species with different effective population sizes for this purpose (Loewe and Charlesworth 2006; Loewe et al. 2006), using Drosophila miranda and D. pseudoobscura as material. These are well suited for this type of study, as they are closely related, live together in similar habitats, and yet have very different levels of silent nucleotide diversity, indicating different effective population sizes (N_e). This study was hampered by our inability to compare the same set of loci across the two species and by the small number of loci that could be used. We here present the results of a much larger study of DNA variation at X-linked and autosomal loci for these two species, using D. affinis as a basis for estimating divergence. We compare the results, applying the method of Loewe et al. (2006) with that of Eyre-Walker and Keightley (2009) for estimating the distribution of deleterious selection coefficients and with McDonald–Kreitman test-based methods for estimating the proportion of nonsynonymous differences fixed by positive selection. While broadly confirming the conclusions from earlier studies, we note some possible sources of bias and describe methods for minimizing their effects.     

Identification of the Major Sex-Determining Region of Turbot (Scophthalmus maximus)

Sex determination in fish is a labile character in evolutionary terms. The sex-determining (SD) master gene can differ even between closely related fish species. This group is an interesting model for studying the evolution of the SD region and the gonadal differentiation pathway. The turbot (Scophthalmus maximus) is a flatfish of great commercial value, where a strong sexual dimorphism exists for growth rate. Following a QTL and marker association approach in five families and a natural population, we identified the main SD region of turbot at the proximal end of linkage group (LG) 5, close to the SmaUSC-E30 marker. The refined map of this region suggested that this marker would be 2.6 cM and 1.4 Mb from the putative SD gene. This region appeared mostly undifferentiated between males and females, and no relevant recombination frequency differences were detected between sexes. Comparative genomics of LG5 marker sequences against five model species showed no similarity of this chromosome to the sex chromosomes of medaka, stickleback, and fugu, but suggested a similarity to a sex-associated QTL from Oreochromis spp. The segregation analysis of the closest markers to the SD region demonstrated a ZW/ZZ model of sex determination in turbot. A small proportion of families did not fit perfectly with this model, which suggests that other minor genetic and/or environmental factors are involved in sex determination in this species.SEX ratio is a central demographic parameter directly related to the reproductive potential of individuals and populations (Penman and Piferrer 2008). The phenotypic sex depends on the processes of both sex determination and sex differentiation. Exogenous factors, such as temperature, hormones, or social behavior, can modify the gonad development pathway in fish (Baroiller and D''Cotta 2001; Piferrer and Guiguen 2008). Both genetic (GSD) and environmental sex determination has been reported in this group (Devlin and Nagahama 2002; Penman and Piferrer 2008), although primary sex determination is genetic in most species (Valenzuela et al. 2003). Among GSD, single, multiple, or polygenic sex-determining (SD) gene systems have been documented (Kallman 1984; Matsuda et al. 2002; Lee et al. 2004; Vandeputte et al. 2007).Sex determination in fish can evolve very rapidly (Woram et al. 2003; Peichel et al. 2004; Ross et al. 2009). Different sex determination mechanisms have been reported between congeneric species and even between populations of the same species (Almeida-Toledo and Foresti 2001; Lee et al. 2004; Mank et al. 2006). The evolution of sex chromosomes involves the suppression of recombination between homologous chromosomes probably to maintain sex-related coadapted gene blocks (Charlesworth et al. 2005; Tripathi et al. 2009). The sex determination pathway appears to be less conserved than other developmental processes (Penman and Piferrer 2008). However, differences are more related to the top of the hierarchy in the developmental pathway, while downstream genes are more conserved (Wilkins 1995; Marín and Baker 1998). As a consequence, the SD master gene in fish can vary among related species (Kondo et al. 2003; Tanaka et al. 2007; Alfaqih et al. 2009). In this sense, fish represent an attractive model for studying the evolution of SD mechanisms and sex chromosomes (Peichel et al. 2004; Kikuchi et al. 2007).A low proportion of fish species have demonstrated sex-associated chromosome heteromorphisms (Almeida-Toledo and Foresti 2001; Devlin and Nagahama 2002; Penman and Piferrer 2008). This is congruent with the rapid evolution of the SD region in fish, and thus in most species the male and female version of this chromosome region appears largely undifferentiated. In spite of this, indirect clues related to progenies of sex/chromosome-manipulated individuals or to segregation of morphologic/molecular sex-associated markers indicate that mechanisms of sex determination in fish are similar to other vertebrates (Penman and Piferrer 2008). With the arrival of genomics, large amounts of different genetic markers and genomic information are available for scanning genomes to look for their association with sex determination. Quantitative trait loci (QTL) (Cnaani et al. 2004; Peichel et al. 2004) or marker association (Felip et al. 2005; Chen et al. 2007) approaches have been used to identify the SD regions in some fish species. Also, microarrays constructed from gonadal ESTs have been applied to detect differentially expressed genes in the process of gonadal differentiation (Baron et al. 2005). Further, the increased genomic resources in model and aquaculture species have allowed the development of both comparative genomics (Woram et al. 2003; Kikuchi et al. 2007; Tripathi et al. 2009) and candidate gene (Shirak et al. 2006; Alfaqih et al. 2009) strategies to identify and characterize the SD region in fish. This has permitted the identification of the SD region in eight fish, including both model and aquaculture species (reviewed in Penman and Piferrer 2008).The turbot is a highly appreciated European aquaculture species, whose harvest is expected to increase from the current 9000 tons to >15,000 tons in 2012 (S. Cabaleiro, personal communication). Females of this species reach commercial size 4–6 months before males do, explaining the interest of the industry in obtaining all-female populations. Although some differences between families can be observed in the production process at farms, sex ratio is usually balanced at ∼1:1. Neither mitotic nor meiotic chromosomes have shown sex-associated heteromorphisms in turbot (Bouza et al. 1994; Cuñado et al. 2001). The proportion of sexes observed in triploid and especially gynogenetic progenies moved Cal et al. (2006a,b) to suggest an XX/XY mechanism in turbot with some additional, either environmental or genetic, factor involved. However, Haffray et al. (2009) have recently claimed a ZZ/ZW mechanism on the basis of the analysis of a large number of progenies from steroid-treated parents. These authors also suggested some (albeit low) influence of temperature in distorting sex proportions after the larval period. Finally, hybridizations between brill (Scophthalmus rhombus) and turbot render monosex progenies, depending on the direction of the cross performed, which suggests different SD mechanisms in these congeneric species (Purdom and ThaCker 1980).In this study, we used the turbot genetic map (Bouza et al. 2007, 2008; Martínez et al. 2008) to look for sex-associated QTL in this species. The identification of a major QTL in a specific linkage group (LG) in the five families analyzed prompted us to refine the genetic map at this LG and to perform a comparative genomics approach against model fish species for a precise location and characterization of the putative SD region. Also, sex-associated QTL markers were screened in a large natural population to provide additional support to our findings and to obtain population parameters at sex-related markers that could aid in interpreting the evolution of this genomic region.     

Mek1 Suppression of Meiotic Double-Strand Break Repair Is Specific to Sister Chromatids,Chromosome Autonomous and Independent of Rec8 Cohesin Complexes

During meiosis, recombination is directed to occur between homologous chromosomes to create connections necessary for proper segregation at meiosis I. Partner choice is determined at the time of strand invasion and is mediated by two recombinases: Rad51 and the meiosis-specific Dmc1. In budding yeast, interhomolog bias is created in part by the activity of a meiosis-specific kinase, Mek1, which is localized to the protein cores of condensed sister chromatids. Analysis of meiotic double-strand break (DSB) repair in haploid and disomic haploid strains reveals that Mek1 suppresses meiotic intersister DSB repair by working directly on sister chromatids. Rec8 cohesin complexes are not required, however, either for suppression of intersister DSB repair or for the repair itself. Regulation of DSB repair in meiosis is chromosome autonomous such that unrepaired breaks on haploid chromosomes do not prevent interhomolog repair between disomic homologs. The pattern of DSB repair in haploids containing Dmc1 and/or Rad51 indicates that Mek1 acts on Rad51-specific recombination processes.IN eukaryotes, meiosis is a specialized type of cell division that produces the gametes required for sexual reproduction. In meiosis, one round of DNA replication is followed by two rounds of chromosome segregation, termed meiosis I and II. As a result of the two divisions, four haploid cells are produced, each containing half the number of chromosomes as the diploid parent. Proper segregation at meiosis I requires connections between homologous chromosomes that are created by a combination of sister chromatid cohesion and recombination (Petronczki et al. 2003). In vegetative cells, cohesion is mediated by multisubunit ring-shaped complexes that are removed by proteolysis of the kleisin subunit, Mcd1/Scc1 (Onn et al. 2008). In meiotic cells, introduction of a meiosis-specific kleisin subunit, Rec8, allows for a two-step removal of cohesion with loss of arm cohesion at anaphase I and centromere cohesion at anaphase II (Klein et al. 1999). Missegregation of chromosomes during meiosis causes abnormal chromosome numbers in gametes that may lead to infertility and genetic disorders such as trisomy 21 or Down''s syndrome.In mitotically dividing budding yeast cells, recombination is mediated by an evolutionarily conserved RecA-like recombinase, Rad51, and occurs preferentially between sister chromatids (Kadyk and Hartwell 1992). In contrast, recombination during meiosis is initiated by the deliberate formation of double-strand breaks (DSBs) by an evolutionarily conserved, topoisomerase-like protein, Spo11, and occurs preferentially between homologous chromosomes (Jackson and Fink 1985; Schwacha and Kleckner 1997; Keeney 2001). After DSB formation, the 5′ ends on either side of the breaks are resected, resulting in 3′ single stranded (ss) tails. Rad51, and the meiosis-specific recombinase Dmc1, bind to the 3′ ssDNA tails to form protein/DNA filaments that promote strand invasion of homologous chromosomes. DNA synthesis and ligation result in the formation of double Holliday junctions, which are then preferentially resolved into crossovers (Allers and Lichten 2001; Hunter 2007).The precise roles that the Rad51 and Dmc1 recombinase activities play in meiotic recombination have been unclear because experiments have indicated both overlapping and distinct functions for the two proteins (Sheridan and Bishop 2006; Hunter 2007). While both rad51Δ and dmc1Δ mutants reduce interhomolog recombination, other studies suggest that Rad51, in complex with the accessory protein Rad54, is involved primarily in intersister DSB repair. In contrast, Dmc1, in conjunction with the accessory protein Rdh54/Tid1 (a paralog of Rad54), effects DSB repair in meiotic cells by invasion of nonsister chromatids (Dresser et al. 1997; Schwacha and Kleckner 1997; Shinohara et al. 1997a,b; Arbel et al. 1999; Bishop et al. 1999; Hayase et al. 2004; Sheridan and Bishop 2006).The preference for recombination to occur between homologous chromosomes during meiosis is created in part by Dmc1. DSBs accumulate in dmc1Δ diploids due to a failure in strand invasion (Bishop et al. 1992; Hunter and Kleckner 2001). In the efficiently sporulating SK1 strain background, these unrepaired breaks trigger the meiotic recombination checkpoint, resulting in prophase arrest (Lydall et al. 1996; Roeder and Bailis 2000). In dmc1Δ mutants, Rad51 is present at DSBs, yet there is no strand invasion of sister chromatids (Bishop 1994; Shinohara et al. 1997a). These results suggest that in addition to Dmc1 promoting interhomolog strand invasion, Rad51 activity must also be suppressed.Recent studies have shown that during meiosis Rad51 recombinase activity is inhibited by two different mechanisms that decrease the formation of Rad51/Rad54 complexes: (1) binding of the meiosis-specific Hed1 protein to Rad51, thereby excluding interaction with Rad54, and (2) reduction in the affinity of Rad54 for Rad51 due to phosphorylation of Rad54 by Mek1 (Tsubouchi and Roeder 2006; Busygina et al. 2008; Niu et al. 2009). Mek1 is a meiosis-specific kinase that is activated in response to DSBs (Niu et al. 2005, 2007; Carballo et al. 2008). In addition to phosphorylating Rad54, Mek1 phosphorylation of an as yet undetermined substrate is required to suppress Rad51/Rad54-mediated strand invasion of sister chromatids (Niu et al. 2009).To dissect the mechanism by which Mek1 suppresses meiotic intersister DSB repair, we took advantage of the ability of yeast cells to undergo haploid meiosis. The lack of homologous chromosomes in haploid cells makes it possible to examine sister-chromatid-specific events in the absence of interhomolog recombination. De Massy et al. (1994) previously observed a delay in DSB repair in haploid cells and proposed that this delay was due to a constraint in using sister chromatids. We have shown that this delay is dependent on MEK1 and utilized the haploid system to determine various biological parameters required to suppress meiotic intersister DSB repair. Our results indicate that Rad51 and Dmc1 recombinase activities have distinct roles during meiosis and that interhomolog bias is established specifically on sister chromatids through regulation of Rad51, not Dmc1. rec8Δ diploids exhibit defects in meiotic DSB repair (Klein et al. 1999; Brar et al. 2009). Given that cohesin complexes are specific for sister chromatids, we investigated the role of REC8 in intersister DSB repair and found it is required neither for suppressing intersister DSB repair during meiosis nor for the repair itself.     

Extended Bayesian LASSO for Multiple Quantitative Trait Loci Mapping and Unobserved Phenotype Prediction

The Bayesian LASSO (BL) has been pointed out to be an effective approach to sparse model representation and successfully applied to quantitative trait loci (QTL) mapping and genomic breeding value (GBV) estimation using genome-wide dense sets of markers. However, the BL relies on a single parameter known as the regularization parameter to simultaneously control the overall model sparsity and the shrinkage of individual covariate effects. This may be idealistic when dealing with a large number of predictors whose effect sizes may differ by orders of magnitude. Here we propose the extended Bayesian LASSO (EBL) for QTL mapping and unobserved phenotype prediction, which introduces an additional level to the hierarchical specification of the BL to explicitly separate out these two model features. Compared to the adaptiveness of the BL, the EBL is “doubly adaptive” and thus, more robust to tuning. In simulations, the EBL outperformed the BL in regard to the accuracy of both effect size estimates and phenotypic value predictions, with comparable computational time. Moreover, the EBL proved to be less sensitive to tuning than the related Bayesian adaptive LASSO (BAL), which introduces locus-specific regularization parameters as well, but involves no mechanism for distinguishing between model sparsity and parameter shrinkage. Consequently, the EBL seems to point to a new direction for QTL mapping, phenotype prediction, and GBV estimation.REGULARIZATION or shrinkage methods are gaining increasing recognition as a valuable alternative to variable selection techniques in dealing with oversaturated or otherwise ill-defined regression problems in both the classical and Bayesian frameworks (e.g., O''hara and Sillanpää 2009). Many studies (e.g., Xu 2003; Wang et al. 2005; Zhang and Xu 2005; De los Campos et al. 2009; Usai et al. 2009; Wu et al. 2009; Xu et al. 2009) have documented the potential of shrinkage methods for quantitative trait locus (QTL) mapping and genomic breeding value (GBV) estimation using genome-wide dense sets of markers. Lee et al. (2008) make a clear connection between phenotype prediction and GBV estimation, suggesting that methods developed for one are also applicable to the other. We thus use the two concepts interchangeably throughout this article.Regularized regression methods, such as ridge regression (Hoerl and Kennard 1970) or the least absolute shrinkage and selection operator (LASSO) (Tibshirani 1996), are essentially penalized likelihood procedures, where suitable penalty functions are added to the negative log-likelihood to automatically shrink spurious effects (effects of redundant covariates) toward zero, while allowing relevant effects to take values farther from zero.It has been pointed out that these non-Bayesian shrinkage methods are not suitable for oversaturated models. Zou and Hastie (2005) and Park and Casella (2008) noted that the LASSO cannot select a number of nonzero effects exceeding the sample size. Xu (2003) found that for ridge regression to work, the number of model effects should be in the same order as the number of observations. This is impractical for genomic selection, which capitalizes on the variation due to small-marker effects, the number of which can exceed the sample size, by contrast to QTL mapping where interest lies mostly in a small subset of loci with large effects on the focal phenotype. In connection with the LASSO, the Bayesian LASSO (BL) (Park and Casella 2008; Yi and Xu 2008) has been proposed to overcome this limitation by imposing a selective shrinkage across regression parameters. Xu (2003) also proposed a Bayesian shrinkage method for QTL mapping, which extends ridge regression in a similar fashion.Although the BL has been successfully applied to QTL mapping (e.g., Yi and Xu 2008) and to GBV estimation (e.g., De los Campos et al. 2009), it relies on a single parameter known as the regularization parameter to simultaneously regulate the overall model sparsity and the extent to which individual regression coefficients are shrunken. However, this is unrealistic when dealing with a large number of predictors whose effect sizes may differ by orders of magnitude. It is therefore natural to ask whether this practice can be relaxed and how such an attempt may impinge on the model performance (e.g., Sun et al. 2010).Here we propose an extension to the Bayesian LASSO for QTL mapping and unobserved phenotype prediction. Our method, the extended Bayesian LASSO (EBL), introduces locus-specific regularization parameters and utilizes a parameterization that clearly separates the overall model sparsity from the degree of shrinkage of individual regression parameters. We use simulated data to investigate the performance of the EBL relative to the Bayesian LASSO in mapping QTL and in predicting unobserved phenotypes. We also compare the performance of the EBL to the Bayesian adaptive LASSO (BAL) recently proposed by Sun et al. (2010), which also assumes locus-specific regularization parameters.     

Deletion in a Quantitative Trait Gene qPE9-1 Associated With Panicle Erectness Improves Plant Architecture During Rice Domestication

Rice plant architecture is an important agronomic trait and a major determinant in high productivity. Panicle erectness is the preferred plant architecture in japonica rice, but the molecular mechanism underlying domestication of the erect panicle remains elusive. Here we report the map-based cloning of a major quantitative trait locus, qPE9-1, which plays an integral role in regulation of rice plant architecture including panicle erectness. The R6547 qPE9-1 gene encodes a 426-amino-acid protein, homologous to the keratin-associated protein 5-4 family. The gene is composed of three Von Willebrand factor type C domains, one transmembrane domain, and one 4-disulfide-core domain. Phenotypic comparisons of a set of near-isogenic lines and transgenic lines reveal that the functional allele (qPE9-1) results in drooping panicles, and the loss-of-function mutation (qpe9-1) leads to more erect panicles. In addition, the qPE9-1 locus regulates panicle and grain length, grain weight, and consequently grain yield. We propose that the panicle erectness trait resulted from a natural random loss-of-function mutation for the qPE9-1 gene and has subsequently been the target of artificial selection during japonica rice breeding.THE worldwide explosion of the human population necessitates an increase in grain yield, which poses a substantial challenge (Rosegrant and Cline 2003). Improvement of plant architecture is considered as a viable approach to increase grain yield, because crop plants with desirable architecture are able to produce much higher yields (Wang and Li 2008). The most striking example arose in the late 1950s, when selection for the semi-dwarf stature in rice and wheat greatly improved plant architecture and yield potential (Peng et al. 1999; Monna et al. 2002; Sasaki et al. 2002; Spielmeyer et al. 2002). Tiller, panicle, and leaf morphology also play important roles in shaping high-yield crop architecture. Most plant architecture traits are controlled by quantitative trait loci (QTL) derived from naturally occurring allelic variation. Rice (Oryza sativa L.) is the most important food crop in the world (White 1994). It is the staple of diet for heavily populated Asian countries as well as many African countries. Numerous QTL or major genes controlling plant architecture traits have been identified and several have recently been cloned (Li et al. 2004; Ashikari et al. 2005; Fan et al. 2006; Xie et al. 2006; Song et al. 2007; Yu et al. 2007; Jin et al. 2008; Shomura et al. 2008; Tan et al. 2008; Xing et al. 2008; Xue et al. 2008). Cloning and functional characterization of these genes not only addresses fundamental questions in plant development, but also facilitates bridging the gap between gene identification and breeding application by improving the precision and efficiency of selection.Rice panicle architecture not only contributes to grain yield, but also to the ecological conditions of cultivated populations and the physicochemical properties of different varieties (Xu et al. 1996; Yuan 1997; Chen et al. 2001). Presently, most japonica rice varieties cultivated in China exhibit the panicle erectness (PE) type of inflorescence (Zhang et al. 2002b). PE varieties typically bear short, erect panicles and leaves, which benefit ventilation and light penetration. As a result, populations of PE varieties show higher photosynthetic rates and material production capacity (Liu et al. 2001; Zhang et al. 2002a; Chen et al. 2007). Additionally, PE rice varieties show increased lodging and fertilizer resistance due to decreased plant height (Xu et al. 1995). Therefore, PE is the preferred plant architecture for high-yield japonica rice.Since development of the first rice PE variety, Guihuahuang, in the early 1960s, a large number of PE japonica varieties have been released in China, including the most well known, Liaojing 5. PE varieties have increased yield potential compared to panicle drooping varieties and therefore are the preferred. PE serves as the most suitable morphological index, and has subsequently been brought into super high-yield breeding. The development and cultivation of PE varieties is considered the third landmark trait after dwarf and hybrid rice in the history of Chinese rice breeding (Zhang et al. 2002b).The genetic mechanisms controlling PE have received some attention. Initially, PE was reported to be governed by a recessive gene (Zhu and Gu 1979), while other studies suggested a major gene with dominant or additive effects, and polygenic modifications serving to regulate PE (Xu et al. 1995; Wang et al. 1997). Chen et al. (2006) proposed that panicle angle was controlled by two major genes with additive-dominance-epistatic effects and also polygenes with additive-dominance-epistatic influences, using major gene–polygene mixed inheritance models and a joint analysis method. Pedigree analysis of PE varieties indicated that two-thirds of the varieties possessed genes from the Italian Balilla variety and shared a close relationship with Liaojing 5 (Zhang et al. 2002b). The dominant EP gene was first reported from chromosome 9, between the two SSR markers RM5833-11 and RM5686-23, at a genetic distance of 1.5 and 0.9 cM, respectively (Kong et al. 2007). In a previous study, we identified and characterized qPE9-1, a major QTL on chromosome 9 responsible for the erect panicle trait using a double-haploid (DH) population derived from a cross between Wuyunjing 8 and Nongken 57 varieties (Yan et al. 2007).However, despite some progresses of the molecular mechanisms governing rice PE, the complexity of the trait results in substantial gaps in our understanding of its regulation. Here we report on a major QTL, qPE9-1, which encodes a keratin-associated protein 5-4, regulates rice PE, and plays pleiotropic roles in an array of plant architecture and yield traits.     

The Conserved miR-51 microRNA Family Is Redundantly Required for Embryonic Development and Pharynx Attachment in Caenorhabditis elegans

A major question about cytokinesis concerns the role of the septin proteins, which localize to the division site in all animal and fungal cells but are essential for cytokinesis only in some cell types. For example, in Schizosaccharomyces pombe, four septins localize to the division site, but deletion of the four genes produces only a modest delay in cell separation. To ask if the S. pombe septins function redundantly in cytokinesis, we conducted a synthetic-lethal screen in a septin-deficient strain and identified seven mutations. One mutation affects Cdc4, a myosin light chain that is an essential component of the cytokinetic actomyosin ring. Five others cause frequent cell lysis during cell separation and map to two loci. These mutations and their dosage suppressors define a signaling pathway (including Rho1 and a novel arrestin) for repairing cell-wall damage. The seventh mutation affects the poorly understood RNA-binding protein Scw1 and severely delays cell separation when combined either with a septin mutation or with a mutation affecting the septin-interacting, anillin-like protein Mid2, suggesting that Scw1 functions in a pathway parallel to that of the septins. Taken together, our results suggest that the S. pombe septins participate redundantly in one or more pathways that cooperate with the actomyosin ring during cytokinesis and that a septin defect causes septum defects that can be repaired effectively only when the cell-integrity pathway is intact.THE fission yeast Schizosaccharomyces pombe provides an outstanding model system for studies of cytokinesis (McCollum and Gould 2001; Balasubramanian et al. 2004; Pollard and Wu 2010). As in most animal cells, successful cytokinesis in S. pombe requires an actomyosin ring (AMR). The AMR begins to assemble at the G₂/M transition and involves the type II myosin heavy chains Myo2 and Myp2 and the light chains Cdc4 and Rlc1 (Wu et al. 2003). Myo2 and Cdc4 are essential for cytokinesis under all known conditions, Rlc1 is important at all temperatures but essential only at low temperatures, and Myp2 is essential only under stress conditions. As the AMR constricts, a septum of cell wall is formed between the daughter cells. The primary septum is sandwiched by secondary septa and subsequently digested to allow cell separation (Humbel et al. 2001; Sipiczki 2007). Because of the internal turgor pressure of the cells, the proper assembly and structural integrity of the septal layers are essential for cell survival.Septum formation involves the β-glucan synthases Bgs1/Cps1/Drc1, Bgs3, and Bgs4 (Ishiguro et al. 1997; Le Goff et al. 1999; Liu et al. 1999, 2002; Martín et al. 2003; Cortés et al. 2005) and the α-glucan synthase Ags1/Mok1 (Hochstenbach et al. 1998; Katayama et al. 1999). These synthases are regulated by the Rho GTPases Rho1 and Rho2 and the protein kinase C isoforms Pck1 and Pck2 (Arellano et al. 1996, 1997, 1999; Nakano et al. 1997; Hirata et al. 1998; Calonge et al. 2000; Sayers et al. 2000; Ma et al. 2006; Barba et al. 2008; García et al. 2009b). The Rho GTPases themselves appear to be regulated by both GTPase-activating proteins (GAPs) and guanine-nucleotide-exchange factors (GEFs) (Nakano et al. 2001; Calonge et al. 2003; Iwaki et al. 2003; Tajadura et al. 2004; Morrell-Falvey et al. 2005; Mutoh et al. 2005; García et al. 2006, 2009a,b). In addition, septum formation and AMR function appear to be interdependent. In the absence of a normal AMR, cells form aberrant septa and/or deposit septal materials at random locations, whereas a mutant defective in septum formation (bgs1) is also defective in AMR constriction (Gould and Simanis 1997; Le Goff et al. 1999; Liu et al. 1999, 2000). Both AMR constriction and septum formation also depend on the septation initiation network involving the small GTPase Spg1 (McCollum and Gould 2001; Krapp and Simanis 2008). Despite this considerable progress, many questions remain about the mechanisms and regulation of septum formation and its relationships to the function of the AMR.One major question concerns the role(s) of the septins. Proteins of this family are ubiquitous in fungal and animal cells and typically localize to the cell cortex, where they appear to serve as scaffolds and diffusion barriers for other proteins that participate in a wide variety of cellular processes (Longtine et al. 1996; Gladfelter et al. 2001; Hall et al. 2008; Caudron and Barral 2009). Despite the recent progress in elucidating the mechanisms of septin assembly (John et al. 2007; Sirajuddin et al. 2007; Bertin et al. 2008; McMurray and Thorner 2008), the details of septin function remain obscure. However, one prominent role of the septins and associated proteins is in cytokinesis. Septins concentrate at the division site in every cell type that has been examined, and in Saccharomyces cerevisiae (Hartwell 1971; Longtine et al. 1996; Lippincott et al. 2001; Dobbelaere and Barral 2004) and at least some Drosophila (Neufeld and Rubin 1994; Adam et al. 2000) and mammalian (Kinoshita et al. 1997; Surka et al. 2002) cell types, the septins are essential for cytokinesis. In S. cerevisiae, the septins are required for formation of the AMR (Bi et al. 1998; Lippincott and Li 1998). However, this cannot be their only role, because the AMR itself is not essential for cytokinesis in this organism (Bi et al. 1998; Korinek et al. 2000; Schmidt et al. 2002). Moreover, there is no evidence that the septins are necessary for AMR formation or function in any other organism. A further complication is that in some cell types, including most Caenorhabditis elegans cells (Nguyen et al. 2000; Maddox et al. 2007) and some Drosophila cells (Adam et al. 2000; Field et al. 2008), the septins do not appear to be essential for cytokinesis even though they localize to the division site.S. pombe has seven septins, four of which (Spn1, Spn2, Spn3, and Spn4) are expressed in vegetative cells and localize to the division site shortly before AMR constriction and septum formation (Longtine et al. 1996; Berlin et al. 2003; Tasto et al. 2003; Wu et al. 2003; An et al. 2004; Petit et al. 2005; Pan et al. 2007; Onishi et al. 2010). Spn1 and Spn4 appear to be the core members of the septin complex (An et al. 2004; McMurray and Thorner 2008), and mutants lacking either of these proteins do not assemble the others at the division site. Assembly of a normal septin ring also depends on the anillin-like protein Mid2, which colocalizes with the septins (Berlin et al. 2003; Tasto et al. 2003). Surprisingly, mutants lacking the septins are viable and form seemingly complete septa with approximately normal timing. These mutants do, however, display a variable delay in separation of the daughter cells, suggesting that the septins play some role(s) in the proper completion of the septum or in subsequent processes necessary for cell separation (Longtine et al. 1996; An et al. 2004; Martín-Cuadrado et al. 2005).It is possible that the septins localize to the division site and yet are nonessential for division in some cell types because their role is redundant with that of some other protein(s) or pathway(s). To explore this possibility in S. pombe, we screened for mutations that were lethal in combination with a lack of septins. The results suggest that the septins cooperate with the AMR during cytokinesis and that, in the absence of septin function, the septum is not formed properly, so that an intact system for recognizing and repairing cell-wall damage becomes critical for cell survival.     

Recurrent Selection on the Winters sex-ratio Genes in Drosophila simulans

Selfish genes, such as meiotic drive elements, propagate themselves through a population without increasing the fitness of host organisms. X-linked (or Y-linked) meiotic drive elements reduce the transmission of the Y (X) chromosome and skew progeny and population sex ratios, leading to intense conflict among genomic compartments. Drosophila simulans is unusual in having a least three distinct systems of X chromosome meiotic drive. Here, we characterize naturally occurring genetic variation at the Winters sex-ratio driver (Distorter on the X or Dox), its progenitor gene (Mother of Dox or MDox), and its suppressor gene (Not Much Yang or Nmy), which have been previously mapped and characterized. We survey three North American populations as well as 13 globally distributed strains and present molecular polymorphism data at the three loci. We find that all three genes show signatures of selection in North America, judging from levels of polymorphism and skews in the site-frequency spectrum. These signatures likely result from the biased transmission of the driver and selection on the suppressor for the maintenance of equal sex ratios. Coalescent modeling indicates that the timing of selection is more recent than the age of the alleles, suggesting that the driver and suppressor are coevolving under an evolutionary “arms race.” None of the Winters sex-ratio genes are fixed in D. simulans, and at all loci we find ancestral alleles, which lack the gene insertions and exhibit high levels of nucleotide polymorphism compared to the derived alleles. In addition, we find several “null” alleles that have mutations on the derived Dox background, which result in loss of drive function. We discuss the possible causes of the maintenance of presence–absence polymorphism in the Winters sex-ratio genes.MEIOTIC drive can leave signatures in the genome similar to positive natural selection without increasing the fitness of an organism (Lyttle 1993). Drive elements are preferentially transmitted during meiosis by disrupting the development or function of sperm carrying the homologous chromosome (Zimmering et al. 1970, meiotic drive sensu lato), or by true chromosome segregation defects during meiosis (Sandler and Novitski 1957, meiotic drive sensu stricto; Tao et al. 2007a). While drive elements may arise on any chromosome, sex-linked drivers have higher population invasion probabilities than autosomal drivers and are more easily detected due to their impact on progeny sex ratios (Hurst and Pomiankowski 1991). To survive, a driver must maintain tight linkage with an insensitive target locus lest it drive against itself, a condition ensured by the lack of recombination between sex chromosomes (Charlesworth and Hartl 1978). Because of the impact drive elements have on sex ratios, sex-linked drivers are often referred to as “sex-ratio distorters” and the phenotype of skewed progeny sex ratios is termed “sex-ratio.” The mere transmission advantage of a driver, unless balanced by some detrimental fitness effect or masked by a suppressor, can cause it to sweep through a population in a manner similar to a positively selected mutation (Edwards 1961; Vaz and Carvalho 2004).Obviously, a complete sweep of a sex-linked driver dooms a male-less (or female-less) population to extinction (Hamilton 1967), and natural selection strongly favors genetic factors that suppress drive and restore Mendelian segregation. Fisher (1930) presented a qualitative argument for the maintenance of an equal sex ratio, which predicts selection on any heritable variant that increases the production of the rarer sex. Fisher''s principle has been formalized mathematically and demonstrated empirically (e.g., Bodmer and Edwards 1960; Carvalho et al. 1998). Suppressors have been identified in a wide variety of meiotic drive systems and are predicted to be strongly favored by natural selection for the maintenance of equal sex ratios (reviewed by Jaenike 2001). Furthermore, the evolution of linked enhancer genes may enable drivers to evade suppression, setting off another bout of Fisherian selection for equal sex ratios (Hartl 1975).Meiotic drive is widespread, with systems identified in mammals, insects, and plants (Jaenike 2001). Drosophila is the most extensively studied insect taxon, and sex-chromosome meiotic drive systems have been identified in more than a dozen species (Jaenike 2001). Cryptic (i.e., suppressed) distorters may be identified when the association between driver and suppressor is lost, such as in hybrids between species or populations that do not share meiotic drive systems (Mercot et al. 1995). The coevolutionary arms race between drivers and suppressors likely contributes to Haldane''s rule (the preferential sterility or inviability of heterogametic hybrids) and is a leading explanation for the importance of X-linked loci in causing hybrid male sterility (Frank 1991; Hurst and Pomiankowski 1991; Tao et al. 2007b; Presgraves 2008). Indeed, two recently characterized hybrid male sterility factors are also sex-ratio distorters—direct evidence of a link between meiotic drive and speciation (Tao et al. 2001; Orr and Irving 2005; Phadnis and Orr 2009).The three X-linked drive systems of Drosophila simulans are genetically distinct and have been termed Paris, Durham, and Winters (Tao et al. 2007a). Here, we focus on the Winters sex-ratio (SR), whose driver and suppressor have been mapped to the gene level and whose molecular and cellular features have been elucidated (Tao et al. 2007a,b). Distortion requires two genes, Distorter on the X (Dox) and Mother of Dox (MDox); Dox is a duplicate copy of MDox (Tao et al. 2007a; Y. Tao, personal communication). The dominant suppressor, Not Much Yang (Nmy), is a retrotransposed copy of Dox on chromosome 3R (Tao et al. 2007b). Nmy likely suppresses Dox through an RNA interference mechanism by forming a double stranded RNA with homology to the distorter RNAs (Tao et al. 2007b). The genes of the Winters sex-ratio are not found in D. melanogaster, which diverged from D. simulans ∼2.3 million years ago (Li et al. 1999). Initial surveys of the genes in the simulans clade indicate that a functional Nmy gene is present in D. mauritiana (Tao et al. 2007b). Thus, the Winters genes are >250,000 years old, the speciation time of D. simulans, D. mauritiana, and D. sechellia (McDermott and Kliman 2008).Signatures of positive selection have been previously detected at genomic regions linked to Drosophila sex-ratio distorters. However, this study represents the first evidence of selection acting directly on a sex-ratio distorter gene and its suppressor gene. In D. recens, driving X chromosomes show reduced nucleotide and haplotype variability relative to standard (nondriving) X chromosomes, and linkage disequilibrium extends over 130 cM of the driving chromosome (Dyer et al. 2007). The Paris driver has been localized to a pair of duplicated loci 150 kb apart; recent work shows reduced haplotype diversity and linkage disequilibrium between variants associated with drive (Derome et al. 2008). In this study, we characterize patterns of genetic variation in natural populations of North American D. simulans and find signatures of recent and strong positive selection at all three genes of the Winters sex-ratio.