The Effect of Recent Admixture on Inference of Ancient Human Population History

Despite the widespread study of genetic variation in admixed human populations, such as African-Americans, there has not been an evaluation of the effects of recent admixture on patterns of polymorphism or inferences about population demography. These issues are particularly relevant because estimates of the timing and magnitude of population growth in Africa have differed among previous studies, some of which examined African-American individuals. Here we use simulations and single-nucleotide polymorphism (SNP) data collected through direct resequencing and genotyping to investigate these issues. We find that when estimating the current population size and magnitude of recent growth in an ancestral population using the site frequency spectrum (SFS), it is possible to obtain reasonably accurate estimates of the parameters when using samples drawn from the admixed population under certain conditions. We also show that methods for demographic inference that use haplotype patterns are more sensitive to recent admixture than are methods based on the SFS. The analysis of human genetic variation data from the Yoruba people of Ibadan, Nigeria and African-Americans supports the predictions from the simulations. Our results have important implications for the evaluation of previous population genetic studies that have considered African-American individuals as a proxy for individuals from West Africa as well as for future population genetic studies of additional admixed populations.STUDIES of archeological and genetic data show that anatomically modern humans originated in Africa and more recently left Africa to populate the rest of the world (Tishkoff and Williams 2002; Barbujani and Goldstein 2004; Garrigan and Hammer 2006; Reed and Tishkoff 2006; Campbell and Tishkoff 2008; Jakobsson et al. 2008; Li et al. 2008). Given the central role Africa has played in the origin of diverse human populations, understanding patterns of genetic variation and the demographic history of populations within Africa is important for understanding the demographic history of global human populations. The availability of large-scale single-nucleotide polymorphism (SNP) data sets coupled with recent advances in statistical methodology for inferring parameters in population genetic models provides a powerful means of accomplishing these goals (Keinan et al. 2007; Boyko et al. 2008; Lohmueller et al. 2009; Nielsen et al. 2009).It is important to realize that studies of African demographic history using genetic data have come to qualitatively different conclusions regarding important parameters. Some recent studies have found evidence for ancient (>100,000 years ago) two- to fourfold growth in African populations (Adams and Hudson 2004; Marth et al. 2004; Keinan et al. 2007; Boyko et al. 2008). Other studies have found evidence of very recent growth (Pluzhnikov et al. 2002; Akey et al. 2004; Voight et al. 2005; Cox et al. 2009; Wall et al. 2009) or could not reject a model with a constant population size (Pluzhnikov et al. 2002; Voight et al. 2005). It is unclear why studies found such different parameter estimates. However, these studies all differ from each other in the amount of data considered, the types of data used (e.g., SNP genotypes vs. full resequencing), the genomic regions studied (e.g., noncoding vs. coding SNPs), and the types of demographic models considered (e.g., including migration vs. not including migration postseparation of African and non-African populations).Another important way in which studies of African demographic history differ from each other is in the populations sampled. Some studies have focused on genetic data from individuals sampled from within Africa (Pluzhnikov et al. 2002; Adams and Hudson 2004; Voight et al. 2005; Keinan et al. 2007; Cox et al. 2009; Wall et al. 2009), while other studies included American individuals with African ancestry (Adams and Hudson 2004; Akey et al. 2004; Marth et al. 2004; Boyko et al. 2008). While there is no clear correspondence between those studies which sampled native African individuals (as opposed to African-Americans) and particular growth scenarios, it is clear from previous studies that African-American populations do differ from African populations in their recent demographic history. In particular, genetic studies suggest that there is wide variation in the degree of European admixture in most African-American individuals in the United States and that they have, on average, ∼80% African ancestry and 20% European ancestry (Parra et al. 1998; Pfaff et al. 2001; Falush et al. 2003; Patterson et al. 2004; Tian et al. 2006; Lind et al. 2007; Reiner et al. 2007; Price et al. 2009; Bryc et al. 2010). Furthermore, both historical records and genetic evidence suggest that the admixture process began quite recently, within the last 20 generations (Pfaff et al. 2001; Patterson et al. 2004; Seldin et al. 2004; Tian et al. 2006). Recent population admixture can alter patterns of genetic variation in a discernible and predictable way. For example, recently admixed populations will exhibit correlation in allele frequencies (i.e., linkage disequilibrium) among markers that differ in frequency between the parental populations. This so-called admixture linkage disequilibrium (LD) (Chakraborty and Weiss 1988) can extend over long physical distances (Lautenberger et al. 2000) and decays exponentially with time the since the admixture process began (i.e., recently admixed populations typically exhibit LD over a longer physical distance than anciently admixed populations).While it is clear that African-American populations have a different recent demographic history than do African populations from within Africa and that admixture tracts can be identified in admixed individuals (Falush et al. 2003; Patterson et al. 2004; Tang et al. 2006; Sankararaman et al. 2008a,b; Price et al. 2009; Bryc et al. 2010), the effect that admixture has on other patterns of genetic variation remains unclear. For example, Xu et al. (2007) found similar LD decay patterns when comparing African-American and African populations. It is also unclear whether the recent admixture affects our ability to reconstruct ancient demographic events (such as expansions that predate the spread of humans out of Africa) from whole-genome SNP data. Most studies of demographic history have summarized the genome-wide SNP data by allele frequency or haplotype summary statistics. If these summary statistics are not sensitive to the recent European admixture, then the African-American samples may yield estimates of demographic parameters that are close to the true demographic parameters for the ancestral, unsampled, African populations. This would suggest that the differences in growth parameter estimates obtained from African populations cannot be explained by certain studies sampling African-American individuals and others sampling African individuals from within Africa. However, if these statistics are sensitive to recent admixture, then they may give biased estimates of growth parameters.Here, we examine the effect of recent admixture on the estimation of population demography. In particular, we estimate growth parameters from simulated data sets using SNP frequencies as well as a recently developed haplotype summary statistic (Lohmueller et al. 2009). We compare the demographic parameter estimates made from the admixed and nonadmixed populations and find that some parameter estimates are qualitatively similar between the two populations when inferred using allele frequencies. Inferences of growth using haplotype-based approaches appear to be more sensitive to recent admixture than inferences based on SNP frequencies. We discuss implications that our results have for interpreting studies of demography in admixed populations.     

Polymorphic Genes of Major Effect: Consequences for Variation,Selection and Evolution in Arabidopsis thaliana

The importance of genes of major effect for evolutionary trajectories within and among natural populations has long been the subject of intense debate. For example, if allelic variation at a major-effect locus fundamentally alters the structure of quantitative trait variation, then fixation of a single locus can have rapid and profound effects on the rate or direction of subsequent evolutionary change. Using an Arabidopsis thaliana RIL mapping population, we compare G-matrix structure between lines possessing different alleles at ERECTA, a locus known to affect ecologically relevant variation in plant architecture. We find that the allele present at ERECTA significantly alters G-matrix structure—in particular the genetic correlations between branch number and flowering time traits—and may also modulate the strength of natural selection on these traits. Despite these differences, however, when we extend our analysis to determine how evolution might differ depending on the ERECTA allele, we find that predicted responses to selection are similar. To compare responses to selection between allele classes, we developed a resampling strategy that incorporates uncertainty in estimates of selection that can also be used for statistical comparisons of G matrices.THE structure of the genetic variation that underlies phenotypic traits has important consequences for understanding the evolution of quantitative traits (Fisher 1930; Lande 1979; Bulmer 1980; Kimura 1983; Orr 1998; Agrawal et al. 2001). Despite the infinitesimal model''s allure and theoretical tractability (see Orr and Coyne 1992; Orr 1998, 2005a,b for reviews of its influence), evidence has accumulated from several sources (artificial selection experiments, experimental evolution, and QTL mapping) to suggest that genes of major effect often contribute to quantitative traits. Thus, the frequency and role of genes of major effect in evolutionary quantitative genetics have been a subject of intense debate and investigation for close to 80 years (Fisher 1930; Kimura 1983; Orr 1998, 2005a,b). Beyond the conceptual implications, the prevalence of major-effect loci also affects our ability to determine the genetic basis of adaptations and species differences (e.g., Bradshaw et al. 1995, 1998).Although the existence of genes of major effect is no longer in doubt, we still lack basic empirical data on how segregating variation at such genes affects key components of evolutionary process (but see Carrière and Roff 1995). In other words, How does polymorphism at genes of major effect alter patterns of genetic variation and covariation, natural selection, and the likely response to selection? The lack of data stems, in part, from the methods used to detect genes of major effect: experimental evolution (e.g., Bull et al. 1997; Zeyl 2005) and QTL analysis (see Erickson et al. 2004 for a review) often detect such genes retrospectively after they have become fixed in experimental populations or the species pairs used to generate the mapping population. The consequences of polymorphism at these genes on patterns of variation, covariation, selection, and the response to selection—which can be transient (Agrawal et al. 2001)—are thus often unobserved.A partial exception to the absence of data on the effects of major genes comes from artificial selection experiments, in which a substantial evolutionary response to selection in the phenotype after a plateau is often interpreted as evidence for the fixation of a major-effect locus (Frankham et al. 1968; Yoo 1980a,b; Frankham 1980; Shrimpton and Robertson 1988a,b; Caballero et al. 1991; Keightley 1998; see Mackay 1990 and Hill and Caballero 1992 for reviews). However, many of these experiments report only data on the selected phenotype (e.g., bristle number) or, alternatively, the selected phenotype and some measure of fitness (e.g., Frankham et al. 1968, Yoo 1980b; Caballero et al. 1991; Mackay et al. 1994; Fry et al. 1995; Nuzhdin et al. 1995; Zur Lage et al. 1997), making it difficult to infer how a mutation will affect variation, covariation, selection, and evolutionary responses for a suite of traits that might affect fitness themselves. One approach is to document how variation at individual genes of major effect affects the genetic variance–covariance matrix (“G matrix”; Lande 1979), which represents the additive genetic variance and covariance between traits.Although direct evidence for variation at major-effect genes altering patterns of genetic variation, covariation, and selection is rare, there is abundant evidence for the genetic mechanisms that could produce these dynamics. A gene of major effect could have these consequences due to any of at least three genetic mechanisms: (1) pleiotropy, where a gene of major effect influences several traits, including potentially fitness, simultaneously, (2) physical linkage or linkage disequilibrium (LD), in which a gene of major effect is either physically linked or in LD with other genes that influence other traits under selection, and (3) epistasis, in which the allele present at a major-effect gene alters the phenotypic effect of other loci and potentially phenotypes under selection. Evidence for these three evolutionary genetic mechanisms leading to changes in suites of traits comes from a variety of sources, including mutation accumulation experiments (Clark et al. 1995; Fernandez and Lopez-Fanjul 1996), mutation induction experiments (Keightley and Ohnishi 1998), artificial selection experiments (Long et al. 1995), and transposable element insertions (Rollmann et al. 2006). For pleiotropy in particular, major-effect genes that have consequences on several phenotypic traits are well known from the domestication and livestock breeding literature [e.g., myostatin mutations in Belgian blue cattle and whippets (Arthur 1995; Grobet et al. 1997; Mosher et al. 2007), halothane genes in pigs (Christian and Rothschild 1991; Fujii et al. 1991), and Booroola and Inverdale genes in sheep (Amer et al. 1999; Visscher et al. 2000)]. While these data suggest that variation at major-effect genes could—and probably does—influence variation, covariation, and selection on quantitative traits, data on the magnitude of these consequences remain lacking.Recombinant inbred line (RIL) populations are a promising tool for investigating the influence of major-effect loci. During advancement of the lines from F₂''s to RILs, alternate alleles at major-effect genes (and most of the rest of the genome) will be made homozygous, simplifying comparisons among genotypic classes. Because of the high homozygosity, individuals within RILs are nearly genetically identical, facilitating phenotyping of many genotypes under a range of environments. In addition, because of recombination, alternative alleles are randomized across genetic backgrounds—facilitating robust comparisons between sets of lines differing at a major-effect locus.Here we investigate how polymorphism at an artificially induced mutation, the erecta locus in Arabidopsis thaliana, affects the magnitude of these important evolutionary genetic parameters under ecologically realistic field conditions. We use the Landsberg erecta (Ler) × Columbia (Col) RIL population of A. thaliana to examine how variation at a gene of major effect influences genetic variation, covariation, and selection on quantitative traits in a field setting. The Ler × Col RIL population is particularly suitable, because it segregates for an artificially induced mutation at the erecta locus, which has been shown to influence a wide variety of plant traits. The Ler × Col population thus allows a powerful test of the effects of segregating variation at a gene—chosen a priori—with numerous pleiotropic effects. The ERECTA gene is a leucine-rich receptor-like kinase (LRR-RLK) (Torii et al. 1996) and has been shown to affect plant growth rates (El-Lithy et al. 2004), stomatal patterning and transpiration efficiency (Masle et al. 2005; Shpak et al. 2005), bacterial pathogen resistance (Godiard et al. 2003), inflorescence and floral organ size and shape (Douglas et al. 2002; Shpak et al. 2003, 2004), and leaf polarity (Xu et al. 2003; Qi et al. 2004).Specifically, we sought to answer the following questions: (1) Is variation at erecta significantly associated with changes to the G matrix? (2) Is variation at erecta associated with changes in natural selection on genetically variable traits? And (3) is variation at erecta associated with significantly different projected evolutionary responses to selection?     

The Drosophila Claudin Kune-kune Is Required for Septate Junction Organization and Tracheal Tube Size Control

The vertebrate tight junction is a critical claudin-based cell–cell junction that functions to prevent free paracellular diffusion between epithelial cells. In Drosophila, this barrier is provided by the septate junction, which, despite being ultrastructurally distinct from the vertebrate tight junction, also contains the claudin-family proteins Megatrachea and Sinuous. Here we identify a third Drosophila claudin, Kune-kune, that localizes to septate junctions and is required for junction organization and paracellular barrier function, but not for apical-basal polarity. In the tracheal system, septate junctions have a barrier-independent function that promotes lumenal secretion of Vermiform and Serpentine, extracellular matrix modifier proteins that are required to restrict tube length. As with Sinuous and Megatrachea, loss of Kune-kune prevents this secretion and results in overly elongated tubes. Embryos lacking all three characterized claudins have tracheal phenotypes similar to any single mutant, indicating that these claudins act in the same pathway controlling tracheal tube length. However, we find that there are distinct requirements for these claudins in epithelial septate junction formation. Megatrachea is predominantly required for correct localization of septate junction components, while Sinuous is predominantly required for maintaining normal levels of septate junction proteins. Kune-kune is required for both localization and levels. Double- and triple-mutant combinations of Sinuous and Megatrachea with Kune-kune resemble the Kune-kune single mutant, suggesting that Kune-kune has a more central role in septate junction formation than either Sinuous or Megatrachea.EPITHELIA are essential for separating physiologically distinct body compartments and regulating trafficking between them. For proper function, it is imperative that epithelia maintain effective barriers against free paracellular diffusion. To this end, epithelial cells contain occluding junctions, which regulate paracellular permeability. In vertebrates, this is accomplished by tight junctions (TJ), structures that are characterized by regions of close membrane apposition between adjacent cells known as “kissing points” (Tsukita and Furuse 2002). While the TJ is made up of at least 40 different components (Schneeberger and Lynch 2004), the core proteins responsible for the paracellular barrier are the claudins (Angelow et al. 2008).Claudins are four-transmembrane domain proteins that form homo- and heterophilic interactions within the same cell (Furuse et al. 1999; Blasig et al. 2006) and with claudins in adjacent cells (Furuse et al. 1999), thereby establishing the paracellular seal. There are 24 members of the claudin family in mammals, many of which display distinct, tissue-specific expression patterns (Kiuchi-Saishin et al. 2002; Angelow et al. 2008). Mutations in several claudins can cause significant paracellular permeability defects in mice. For example, mutations in claudin-14 increase TJ permeability in the organ of Corti and cause deafness (Ben-Yosef et al. 2003), while loss of claudin-1 compromises epidermal barrier function (Furuse et al. 2002).In Drosophila, primary (ectodermally derived) epithelia lack discernable TJs and instead use pleated septate junctions (SJ) for the paracellular barrier (Baumgartner et al. 1996; Lamb et al. 1998; Genova and Fehon 2003; Paul et al. 2003). However, despite sharing a common barrier function, vertebrate TJs and invertebrate SJs differ in several ways. While vertebrate TJs are positioned apical to adherens junctions (AJ) and contain conserved apical polarity proteins, SJs are basal to AJs and contain conserved basolateral polarity proteins (reviewed in Tepass 2003; Wu and Beitel 2004). In addition, SJs do not contain kissing points, but rather ladder-like septa that span the intermembrane space (Lane and Swales 1982; Tepass and Hartenstein 1994).Beyond their general epithelial barrier function, SJs are also required for several tissue-specific processes. Glial cells, for example, ensheath nerve fibers and use SJs to maintain the blood–brain barrier (Auld et al. 1995; Baumgartner et al. 1996; Schwabe et al. 2005). In the embryonic tracheal system, SJs are required for the apical secretion of the lumenal matrix modifying proteins, Vermiform (Verm) and Serpentine (Serp), which act through undefined pathways to restrict tube length (Wang et al. 2006). This secretory pathway appears to be specific for Verm and Serp, since other apical proteins are secreted normally in SJ mutants. SJ proteins have also been shown to play a role in morphogenesis of the heart tube, even though this tissue lacks typical SJ septa (Yi et al. 2008).Although SJs have clear differences from vertebrate TJs, SJs contain at least two claudins, Megatrachea (Mega) and Sinuous (Sinu), both of which are required for the paracellular barrier (Behr et al. 2003; Wu et al. 2004; Stork et al. 2008). In this article, we identify a third claudin, Kune-kune (Kune), that is an integral SJ protein. Like the other claudins, Kune is required for maintaining epithelial paracellular barrier and tracheal tube size control and is not required for apical-basal polarity. We also find that, of all three characterized claudins, Kune has a more severe SJ phenotype, suggesting that it is a more central player in SJ organization and function than previously characterized Drosophila claudins.     

Deletion of a Novel F-Box Protein,MUS-10, in Neurospora crassa Leads to Altered Mitochondrial Morphology,Instability of mtDNA and Senescence

While mitochondria are renowned for their role in energy production, they also perform several other integral functions within the cell. Thus, it is not surprising that mitochondrial dysfunction can negatively impact cell viability. Although mitochondria have received an increasing amount of attention in recent years, there is still relatively little information about how proper maintenance of mitochondria and its genomes is achieved. The Neurospora crassa mus-10 mutant was first identified through its increased sensitivity to methyl methanesulfonate (MMS) and was thus believed to be defective in some aspect of DNA repair. Here, we report that mus-10 harbors fragmented mitochondria and that it accumulates deletions in its mitochondrial DNA (mtDNA), suggesting that the mus-10 gene product is involved in mitochondrial maintenance. Interestingly, mus-10 begins to senesce shortly after deletions are visualized in its mtDNA. To uncover the function of MUS-10, we used a gene rescue approach to clone the mus-10 gene and discovered that it encodes a novel F-box protein. We show that MUS-10 interacts with a core component of the Skp, Cullin, F-box containing (SCF) complex, SCON-3, and that its F-box domain is essential for its function in vivo. Thus, we provide evidence that MUS-10 is part of an E3 ubiquitin ligase complex involved in maintaining the integrity of mitochondria and may function to prevent cellular senescence.THE mus-10 mutant was isolated from a screen aimed at identifying Neurospora crassa strains that were sensitive to MMS and therefore likely to lack proper DNA repair mechanisms (Kafer and Perlmutter 1980). Epistasis analyses involving mus-10 suggested that it belonged to the uvs-6 epistasis group, which functions in recombination repair (Kafer and Perlmutter 1980; Kafer 1983). However, mus-10 did not display several phenotypes common to other members of the uvs-6 epistasis group: chromosomal instability, a high sensitivity to histidine, and the inability to produce viable ascospores in homozygous crosses (Newmeyer et al. 1978; Newmeyer and Galeazzi 1978; Kafer and Perlmutter 1980; Kafer 1981; Schroeder 1986; Watanabe et al. 1997; Handa et al. 2000; Sakuraba et al. 2000). Furthermore, the frequencies of spontaneous and radiation-induced mutation observed in mus-10 were similar to those of a wild-type strain (Kafer 1981). Past efforts to uncover the nature of these discrepancies or the function of the mus-10 gene product have been uninformative.The majority of cellular ATP is produced in mitochondria through aerobic respiration, which couples electron flow through respiratory complexes within the mitochondrial inner membrane with oxidative phosphorylation. Besides their role in ATP synthesis, mitochondria are also involved in many other cellular processes including beta-oxidation (Bartlett and Eaton 2004), calcium homeostasis (Gunter et al. 2004; Rimessi et al. 2008), production of iron-sulfur clusters (Zheng et al. 1998; Gerber and Lill 2002; Lill and Muhlenhoff 2005; Rouault and Tong 2005), and apoptosis (Green 2005; Antignani and Youle 2006; Xu and Shi 2007). Although virtually all mitochondrial proteins are encoded within the nucleus, a small number of proteins are encoded by mitochondrial DNA (mtDNA). The integrity of the mitochondrial genome may affect cell survival as mutations in mtDNA accumulate in patients suffering from severe neurological diseases including Alzheimer''s, Huntington''s and Parkinson''s, as well as several types of cancer (Chatterjee et al. 2006; Higuchi 2007; Krishnan et al. 2007; Reeve et al. 2008). The number of mtDNA mutations also increases with age, suggesting a link between mitochondrial dysfunction and ageing (Cortopassi and Arnheim 1990; Corral-Debrinski et al. 1992; Cortopassi et al. 1992; Simonetti et al. 1992; Reeve et al. 2008). Contrary to the single genome in the nucleus, there are several copies of mtDNA in each mitochondrion. Thus, defects in a few mitochondrial genomes do not necessarily lead to mitochondrial dysfunction. Many patients suffering from mitochondrial diseases exhibit heteroplasmy, a phenomenon in which a mixture of wild-type and mutant mtDNAs exist in a single cell. The ratio of wild-type to mutant mtDNAs is critical in determining the penetrance of the genetic defect, where mutant loads >60% are required to cause respiratory chain dysfunction within an individual cell (Boulet et al. 1992; Chomyn et al. 1992; Sciacco et al. 1994).Even though N. crassa strains are generally deemed immortal if they can be subcultured ∼50 times, a wild-type strain was recently reported to senesce after 12,000 hr of growth, implying that this fungus undergoes natural or programmed ageing (Maheshwari and Navaraj 2008; Kothe et al. 2010). However, replicative life span is also influenced by genetic background as certain mutations can cause progressive deterioration of growth, ultimately leading to death. One such example is the nuclear-encoded natural death (nd), which when mutant causes a senescence phenotype correlating with the accumulation of multiple mtDNA deletions (Sheng 1951; Seidel-Rogol et al. 1989). The deletions of mtDNA in nd occurred between two 70- to 701-bp direct repeats, suggesting that the nd gene product regulates recombination, repair, or replication of mtDNA (Bertrand et al. 1993). Another nuclear mutation, senescence (sen), was isolated from N. intermedia and introgressed into N. crassa (Navaraj et al. 2000). Deletions were also observed in the mtDNA of sen mutants, but unlike those occurring in nd were flanked by 6- to 10-bp repeats typically associated with GC-rich palindromic sequences (D''Souza et al. 2005). The nature of the sequences that flanked the mtDNA deletions in these two mutants supported the existence of two distinct systems of mtDNA recombination in N. crassa: a general system of homologous recombination (system I) and a site-specific mechanism (system II), mediated in part by nd and sen, respectively (Bertrand et al. 1993; D''Souza et al. 2005). The nd and sen mutations have been mapped to linkage groups I and V, respectively, but neither gene has been cloned and the precise function of their gene products remains unclear. Two ultraviolet (UV)-sensitive mutants, uvs-4 and uvs-5, are thought to undergo senescence, but unfortunately, these strains have not been studied in great detail (Schroeder 1970; Perkins et al. 1993; Hausner et al. 2006). Premature senescence has also been observed in cytoplasmic mutants of N. crassa including the E35 and ER-3 stopper mutants that harbor large mtDNA deletions, as well as strains that accumulate mitochondrial plasmids capable of inserting into mtDNA through homologous recombination (de Vries et al. 1986; Akins et al. 1989; Myers et al. 1989; Niagro and Mishra 1989; Court et al. 1991; Alves and Videira 1998).While trying to establish the role of MUS-10 in DNA repair, we discovered that the mus-10 mutant exhibited a shortened life span, an abnormal mitochondrial morphology and mtDNA instability. We cloned the mus-10 gene through its ability to complement the MMS sensitivity of the mus-10 mutant and revealed that it encoded a novel F-box protein. This suggested that MUS-10 is part of an Skp, Cullin, F-box containing (SCF) E3 ubiquitin ligase complex that targets proteins for degradation by the 26S proteasome. The data we present in this article offer proof that an SCF complex can regulate both mitochondrial maintenance and cellular senescence.     

Bayesian Inference of Genetic Parameters Based on Conditional Decompositions of Multivariate Normal Distributions

It is widely recognized that the mixed linear model is an important tool for parameter estimation in the analysis of complex pedigrees, which includes both pedigree and genomic information, and where mutually dependent genetic factors are often assumed to follow multivariate normal distributions of high dimension. We have developed a Bayesian statistical method based on the decomposition of the multivariate normal prior distribution into products of conditional univariate distributions. This procedure permits computationally demanding genetic evaluations of complex pedigrees, within the user-friendly computer package WinBUGS. To demonstrate and evaluate the flexibility of the method, we analyzed two example pedigrees: a large noninbred pedigree of Scots pine (Pinus sylvestris L.) that includes additive and dominance polygenic relationships and a simulated pedigree where genomic relationships have been calculated on the basis of a dense marker map. The analysis showed that our method was fast and provided accurate estimates and that it should therefore be a helpful tool for estimating genetic parameters of complex pedigrees quickly and reliably.MUCH effort in genetics has been devoted to revealing the underlying genetic architecture of quantitative or complex traits. Traditionally, the polygenic model has been used extensively to estimate genetic variances and breeding values of natural and breeding populations, where an infinite number of genes is assumed to code for the trait of interest (Bulmer 1971; Falconer and Mackay 1996). The genetic variance of a quantitative trait can be decomposed into an additive part that corresponds to the effects of individual alleles and a part that is nonadditive because of interactions between alleles. Attention has generally been focused on the estimation of additive genetic variance (and heritability), since additive variation is directly proportional to the response of selection via the breeder''s equation (Falconer and Mackay 1996, Chap. 11). However, to estimate additive genetic variation and heritability accurately, it can be important to identify potential nonadditive sources in genetic evaluations (Misztal 1997; Ovaskainen et al. 2008; Waldmann et al. 2008), especially if the pedigree being analyzed contains a large proportion of full-sibs and clones, as these in particular give rise to nonadditive genetic relationships (Lynch and Walsh 1998, pp. 145). The polygenic model using pedigree and phenotypic information, i.e., the animal model (Henderson 1984), has been the model of choice for estimating genetic parameters in breeding and natural populations (Abney et al. 2000; Sorensen and Gianola 2002; O′Hara et al. 2008).Recent breakthroughs in molecular techniques have made it possible to create genome-wide, single nucleotide polymorphism (SNP) maps. These maps have helped to uncover a vast amount of new loci responsible for trait expression and have provided general insights into the genetic architecture of quantitative traits (e.g., Valdar et al. 2006; Visscher 2008; Flint and Mackay 2009). These insights can help when calculating disease risks in humans, when attempting to increase the yield from breeding programs, and when estimating relatedness in conservation programs. High-density SNPs of many species of importance to science and agriculture can now be scored quickly and relatively cheaply, for example, in mice (Valdar et al. 2006), chickens (Muir et al. 2008), and dairy cattle (VanRaden et al. 2009).In the analysis of populations of breeding stock, the inclusion of dense marker data has improved the predictive ability (i.e., reliability) of genetic evaluations compared to the traditional phenotype model, both in simulations (Meuwissen et al. 2001; Calus et al. 2008; Hayes et al. 2009) and when using real data (Legarra et al. 2008; VanRaden et al. 2009; González-Recio et al. 2009). Meuwissen et al. (2001) suggested that the effect of all markers should first be estimated, and then summed, to obtain genomic estimated breeding values (GEBVs). An alternative procedure, where all markers are used to compute the genomic relationship matrix (in place of the additive polygenic relationship matrix) has also been suggested (e.g., Villanueva et al. 2005; VanRaden 2008; Hayes et al. 2009); this matrix is then incorporated into the statistical analysis to estimate GEBVs. A comparison of both procedures (VanRaden 2008) yielded similar estimates of GEBVs in cases where the effect of an individual allele was small. In addition, if not all pedigree members have marker information, a combined relationship matrix derived from both genotyped and ungenotyped individuals could be computed; this has been shown to increase the accuracy of GEBVs (Legarra et al. 2009; Misztal et al. 2009). Another plausible option to incorporate marker information is to use low-density SNP panels within families and to trace the effect of SNPs from high-density genotyped ancestors, as suggested by Habier et al. (2009) and Weigel et al. (2009). However, fast and powerful computer algorithms, which can use the marker information as efficiently as possible in the analysis of quantitative traits, are needed to obtain accurate GEBVs from genome-wide marker data.This study describes the development of an efficient Bayesian method for incorporating general relationships into the genetic evaluation procedure. The method is based on expressing the multivariate normal prior distribution as a product of one-dimensional normal distributions, each conditioned on the descending variables. When evaluating the genetic parameters of natural and breeding populations, high-dimensional distributions are often used as prior distributions of various genetic effects, such as the additive polygenic effect (Wang et al. 1993), multivariate additive polygenic effects (Van Tassell and Van Vleck 1996), and quantitative trait loci (QTL) effects via the identical-by-decent matrix (Yi and Xu 2000). A Bayesian framework is adopted to obtain posterior distributions of all unknown parameters, estimated by using Markov chain Monte Carlo (MCMC) sampling algorithms in the software package WinBUGS (Lunn et al. 2000, 2009). By performing prior calculations in the form of the factorized product of simple univariate conditional distributions, the computational time of the MCMC estimation procedure is reduced considerably. This feature permits rapid inference for both the polygenic model and the genomic relationship model. Moreover, the decomposition allows for inbreeding of varying degree, since the correct genetic covariance structure can be inferred into the analysis. In this article, we test the method on two previously published pedigree data sets: phenotype data from a large pedigree of Scots pine, incorporation of information on both additive and dominance genetic relationships (Waldmann et al. 2008); and genomic information obtained from a genome-wide scan of a simulated animal population (Lund et al. 2009).     

Genetic Control of Photoperiod Sensitivity in Maize Revealed by Joint Multiple Population Analysis

Variation in maize for response to photoperiod is related to geographical adaptation in the species. Maize possesses homologs of many genes identified as regulators of flowering time in other species, but their relation to the natural variation for photoperiod response in maize is unknown. Candidate gene sequences were mapped in four populations created by crossing two temperate inbred lines to two photoperiod-sensitive tropical inbreds. Whole-genome scans were conducted by high-density genotyping of the populations, which were phenotyped over 3 years in both short- and long-day environments. Joint multiple population analysis identified genomic regions controlling photoperiod responses in flowering time, plant height, and total leaf number. Four key genome regions controlling photoperiod response across populations were identified, referred to as ZmPR1–4. Functional allelic differences within these regions among phenotypically similar founders suggest distinct evolutionary trajectories for photoperiod adaptation in maize. These regions encompass candidate genes CCA/LHY, CONZ1, CRY2, ELF4, GHD7, VGT1, HY1/SE5, TOC1/PRR7/PPD-1, PIF3, ZCN8, and ZCN19.MAIZE (Zea mays L. subsp. mays) was domesticated in southern Mexico and its center of diversity is in tropical Latin America (Goodman 1999; Matsuoka et al. 2002), where precipitation rates and day lengths cycle annually. The presumed ancestor of maize, teosinte (Zea mays L. subsp. parviglumis), likely evolved photoperiod sensitivity to synchronize its reproductive phases to the wetter, short-day growing season (Ribaut et al. 1996; Campos et al. 2006). A critical event in the postdomestication evolution of maize was its spread from tropical to temperate regions of the Americas (Goodman 1988), requiring adaptation to longer day lengths. The result of this adaptation process is manifested today as a major genetic differentiation between temperate and tropical maize (Liu et al. 2003) and substantially reduced photoperiod sensitivity of temperate maize (Gouesnard et al. 2002). Tropical maize exhibits delayed flowering time, increased plant height, and a greater total leaf number when grown in temperate latitudes with daily dark periods <11 hr (Allison and Daynard 1979; Warrington and Kanemasu 1983a,b). Identifying the genes underlying maize photoperiod sensitivity will provide insight into the postdomestication evolution of maize and may reduce barriers to the use of diverse tropical germplasm resources for improving temperate maize production (Holland and Goodman 1995; Liu et al. 2003; Ducrocq et al. 2009).Natural variation at key genes in flowering time pathways is related to adaptation and evolution of diverse plant species (Caicedo et al. 2004; Shindo et al. 2005; Turner et al. 2005; Cockram et al. 2007; Izawa 2007; Slotte et al. 2007). Identification of some of the genes controlling adaptation in numerous plant species relied on regulatory pathways elucidated in Arabidopsis (Simpson and Dean 2002). Many key genes in the Arabidopsis flowering time regulatory pathways are conserved across diverse plant species (Kojima et al. 2002; Hecht et al. 2007; Kwak et al. 2008), but their functions have diverged, resulting in unique regulatory pathways in some phylogenetic groups (Colasanti and Coneva 2009). For example, FRI and FLC control most natural variation for vernalization response in Arabidopsis (Caicedo et al. 2004; Shindo et al. 2005), but wheat and barley appear to lack homologs of these genes and regulate vernalization response with different genes (Yan et al. 2004).Maize exhibits tremendous natural variation for flowering time (Gouesnard et al. 2002; Camus-Kulandaivelu et al. 2006), for which numerous QTL have been identified (Chardon et al. 2004). In contrast, only a few flowering time mutants are known and only a handful of flowering time genes, including DWARF8 (D8), DELAYED FLOWERING1 (DLF1), VEGETATIVE TO GENERATIVE TRANSITION1 (VGT1), and INDETERMINATE GROWTH1 (ID1), have been cloned in maize (Thornsberry et al. 2001; Colasanti et al. 2006; Muszynski et al. 2006; Salvi et al. 2007; Colasanti and Coneva 2009). Variation at or near D8 and VGT1 is related to latitudinal adaptation, but these genes do not appear to regulate photoperiod responses and account for only a limited proportion of the standing flowering time variation in maize (Camus-Kulandaivelu et al. 2006, 2008; Ducrocq et al. 2008; Buckler et al. 2009).Quantitative trait loci (QTL) mapping was a key first step to identifying the genes underlying natural variation for flowering time in Arabidopsis (Koornneef et al. 2004). Photoperiodic QTL have been mapped previously in individual biparental maize mapping populations (Koester et al. 1993; Moutiq et al. 2002; Wang et al. 2008; Ducrocq et al. 2009). Such studies are informative with respect to the parents from which the populations were derived, but often do not reflect the genetic heterogeneity of broader genetic reference populations (Holland 2007).Association mapping (Thornsberry et al. 2001; Ersoz et al. 2007) and combined analysis of multiple biparental crosses (Rebaï et al. 1997; Rebaï and Goffinet 2000; Blanc et al. 2006; Verhoeven et al. 2006; Yu et al. 2008) represent alternative approaches to understanding the variation in genetic control for complex traits among diverse germplasm. Association mapping has limited power to identify genes that affect traits closely associated with population structure, such as flowering time in maize (Camus-Kulandaivelu et al. 2006; Ersoz et al. 2007). In contrast, joint QTL analysis of multiple populations is not hindered by the associations between causal genes and population structure. Combined QTL analysis of multiple mapping populations provides improved power to detect QTL, more precise estimation of their effects and positions, and better understanding of their functional allelic variation and distribution across more diverse germplasm compared to single-population mapping (Rebaï et al. 1997; Wu and Jannink 2004; Jourjon et al. 2005; Blanc et al. 2006; Verhoeven et al. 2006; Yu et al. 2008; Buckler et al. 2009). Joint analysis also provides a direct test of the importance of higher-order epistatic interactions between founder alleles at individual loci with genetic backgrounds (Jannink and Jansen 2001; Blanc et al. 2006). In this study, joint analysis of multiple populations was used to test directly the hypothesis that diverse tropical maize lines carry functionally similar alleles at key photoperiod loci, which would imply genetic homogeneity for a common set of mutations and a shared evolutionary pathway for photoperiod insensitivity.The objective of this study was to integrate candidate gene analyses with photoperiod QTL mapping across multiple maize populations. We tested candidate floral regulators known from other species for associations with natural variation for photoperiod response in maize. We analyzed flowering time in four interrelated recombinant inbred line (RIL) populations, each derived from crosses between temperate and tropical maize parents (Figure 1), in both long- and short-day environments to characterize their responses to distinct photoperiods. Joint population analysis provided high resolution of many QTL positions, permitting robust testing of underlying candidate genes. We directly and indirectly mapped homologs of flowering time candidates genes from Arabidopsis, rice, and barley on a dense consensus genetic map of these four populations, permitting identification of homologs that colocalize with genome regions associated with variation for photoperiod response. These mapping families are being integrated into the maize nested association mapping (NAM) population (Buckler et al. 2009; McMullen et al. 2009) because they were genotyped with the maize NAM map SNP markers, they involve the common parent B73, and their seed and genotypic information (File S1 cont.) are publicly available. Their availability further expands the genetic diversity represented by the maize NAM population and enhances this valuable public community resource.Open in a separate windowFigure 1.—Factorial mating of two temperate (B73 and B97) and two tropical (CML254 and Ki14) inbred maize lines to create four related recombinant inbred line mapping populations.     

In Planta Mutagenesis Determines the Functional Regions of the Wheat Puroindoline Proteins

In planta analysis of protein function in a crop plant could lead to improvements in understanding protein structure/function relationships as well as selective agronomic or end product quality improvements. The requirements for successful in planta analysis are a high mutation rate, an efficient screening method, and a trait with high heritability. Two ideal targets for functional analysis are the Puroindoline a and Puroindoline b (Pina and Pinb, respectively) genes, which together compose the wheat (Triticum aestivum L.) Ha locus that controls grain texture and many wheat end-use properties. Puroindolines (PINs) together impart soft texture, and mutations in either PIN result in hard seed texture. Studies of the PINs'' mode of action are limited by low allelic variation. To create new Pin alleles and identify critical function-determining regions, Pin point mutations were created in planta via EMS treatment of a soft wheat. Grain hardness of 46 unique PIN missense alleles was then measured using segregating F₂:F₃ populations. The impact of individual missense alleles upon PIN function, as measured by grain hardness, ranged from neutral (74%) to intermediate to function abolishing. The percentage of function-abolishing mutations among mutations occurring in both PINA and PINB was higher for PINB, indicating that PINB is more critical to overall Ha function. This is contrary to expectations in that PINB is not as well conserved as PINA. All function-abolishing mutations resulted from structure-disrupting mutations or from missense mutations occurring near the Tryptophan-rich region. This study demonstrates the feasibility of in planta functional analysis of wheat proteins and that the Tryptophan-rich region is the most important region of both PINA and PINB.NATURAL selection has captured a relatively small subset of potentially useful protein sequences. Unraveling the critical features of proteins via understanding the process of their evolution is a powerful approach for proteins present in many diverse species (Bashford et al. 1987; Hampsey et al. 1988). However, this approach is not feasible for the wheat puroindolines (PINs) that are present only in hexaploid wheat and related species (Massa and Morris 2006). The PINs are unique in structure in having a tryptophan-rich domain and are members of the protease inhibitor/seed storage/lipid transfer protein family (PF00234) (Finn et al. 2008).The tryptophan-rich domain has been hypothesized to control PIN function (Giroux and Morris 1997), but there is no unbiased direct evidence for this since previous studies have focused on the tryptophan box alone (Evrard et al. 2008). A nonbiased approach would consist of random mutagenesis followed by functional analysis (Bowie et al. 1990). This approach has been used extensively for proteins that can be expressed in vitro using either random (Tarun et al. 1998; Guo et al. 2004; Smith and Raines 2006; Georgelis et al. 2007) or site-directed mutations (Miyahara et al. 2008; Osmani et al. 2008). However, functional analysis of many plant proteins in vitro may not be comparable to in planta analysis. In the case of puroindolines, there is no in vitro assay that properly mimics the synergistic binding of PINA and PINB to starch granules or is as easy to measure as grain hardness. Therefore, creation and analysis of a large number of new alleles in wheat in planta is an ideal approach to dissect PIN function.The absence of high-throughput transformation and/or functional screening methods in most crop plants is the largest obstacle in the way of in planta protein functional analysis. However, high-throughput in vitro random or targeted mutagenesis followed by functional analysis has been demonstrated in Arabidopsis thaliana (Dunning et al. 2007) and Nicotiana benthamiana (Boter et al. 2007). Traditional in planta mutagenesis followed by analysis of loss-of-function mutations has been used to clone unknown genes (Xiong et al. 2001) or to define function for candidate genes (Haralampidis et al. 2001; Qi et al. 2006). A high-throughput in planta functional approach for PINA and PINB seems attractive for three reasons. First, the EMS mutation rate in wheat is higher than in any other plant (Slade et al. 2005; Feiz et al. 2009a). Second, PINs control the vast majority of variation in grain hardness (Campbell et al. 1999). Finally, a small-scale preliminary study indicated the feasibility of this approach (Feiz et al. 2009a).PINA and PINB are cysteine-rich proteins unique in having a tryptophan-rich domain (Blochet et al. 1993) and together compose the wheat Hardness (Ha) locus (Giroux and Morris 1998; Wanjugi et al. 2007a). Ha is located on chromosome 5DS and is the major determinant of wheat endosperm texture (Mattern et al. 1973; Law et al. 1978; Campbell et al. 1999). Soft texture (Ha) results when both Pin genes are wild type (Pina-D1a, Pinb-D1a) while hard texture (ha) results from mutations in either Pin (Giroux and Morris 1997, 1998). Transgenic studies in rice (Krishnamurthy and Giroux 2001), wheat (Beecher et al. 2002; Martin et al. 2006), and corn (Zhang et al. 2009) have demonstrated that Pin mutations are causative to hard grain texture. PINA and PINB are not functionally interchangeable and control grain hardness via cooperative binding to starch granules (Hogg et al. 2004; Swan et al. 2006; Wanjugi et al. 2007a; Feiz et al. 2009b). PIN binding to starch granules is mediated by polar lipids (Greenblatt et al. 1995) and PIN abundance is correlated with seed polar lipid content (Feiz et al. 2009b). Variation in PIN function affects grain hardness along with nearly all end product quality traits (Hogg et al. 2005; Martin et al. 2007, 2008; Wanjugi et al. 2007b; Feiz et al. 2008). Determining PINs'' function-determining regions could lead to greater knowledge of their mode of action and to wheat quality improvements. Current PIN functional analyses have been limited to in vitro tests of binding to each other (Ziemann et al. 2008) or to yeast membranes (Evrard et al. 2008).Here, we report the creation and functional analysis in planta of new alleles of PINA and PINB. This is the first successful in planta functional analysis of a crop plant protein.     

Distribution of Fitness Effects Caused by Single-Nucleotide Substitutions in Bacteriophage f1

Empirical knowledge of the fitness effects of mutations is important for understanding many evolutionary processes, yet this knowledge is often hampered by several sources of measurement error and bias. Most of these problems can be solved using site-directed mutagenesis to engineer single mutations, an approach particularly suited for viruses due to their small genomes. Here, we used this technique to measure the fitness effect of 100 single-nucleotide substitutions in the bacteriophage f1, a filamentous single-strand DNA virus. We found that approximately one-fifth of all mutations are lethal. Viable ones reduced fitness by 11% on average and were accurately described by a log-normal distribution. More than 90% of synonymous substitutions were selectively neutral, while those affecting intergenic regions reduced fitness by 14% on average. Mutations leading to amino acid substitutions had an overall mean deleterious effect of 37%, which increased to 45% for those changing the amino acid polarity. Interestingly, mutations affecting early steps of the infection cycle tended to be more deleterious than those affecting late steps. Finally, we observed at least two beneficial mutations. Our results confirm that high mutational sensitivity is a general property of viruses with small genomes, including RNA and single-strand DNA viruses infecting animals, plants, and bacteria.MUTATIONAL fitness effects are relevant to many evolutionary processes. For instance, they determine the fraction of mutations that evolves neutrally (Ohta 1992), the amount of genetic variation at the mutation–selection balance (Haldane 1937), processes of fitness decay, such as Muller''s ratchet (Butcher 1995), mutational meltdown (Lynch et al. 1993), or lethal mutagenesis (Bull et al. 2007), the ability of organisms to fix beneficial mutations and evolve novel functions (Wagner 2005), or the origin of sex and recombination (Peck et al. 1997; de Visser et al. 2003). Considerable progress has been made in characterizing mutational fitness effects using model organisms or studying genetic variation in natural populations (Eyre-Walker and Keightley 2007). For instance, mutation–accumulation experiments suggest that the average effect of spontaneous deleterious mutations is 1% or lower (Kibota and Lynch 1996) in Escherichia coli, while roughly 90% of engineered gene knockouts are viable (Baba et al. 2006) and transposon insertions reduce fitness by 3% or less on average (Elena et al. 1998). In yeast, mutation–accumulation and chemical mutagenesis experiments have shown that mutations reduce fitness by 1–4% on average in diploid strains (Zeyl and de Visser 2001; Szafraniec et al. 2003; Joseph and Hall 2004). In nematodes most mutations have fitness effects lower than 1% (Keightley and Caballero 1997; Davies et al. 1999), in Drosophila the average effect of mutations ranges from 0.5 to 3.5% (Mukai et al. 1972; Ohnishi 1977; Fernández and López-Fanjul 1996; Fry et al. 1999), and, in humans, most segregating amino acid substitutions have fitness effects lower than 10% (Eyre-Walker and Keightley 2007).Although mutation–accumulation studies provide valuable information about the average effects of deleterious mutations, their power to infer the entire distribution of mutational effects, including neutral and lethal mutations, is more limited. Also, excluding bias due to selection can be problematic, and the precise location and nature of each mutation is often unknown. On the other hand, studies based on engineering mutations have been generally restricted to large deletions or insertions, which are probably infrequent in nature compared to point mutations. A direct and powerful approach that helps us to solve these difficulties consists of introducing single-nucleotide substitutions by site-directed mutagenesis. Due to their small genome sizes, viruses are excellent systems for achieving this goal. In previous work, this technique has been used for studying mutational fitness effects in several RNA viruses (Sanjuán et al. 2004; Carrasco et al. 2007; Domingo-Calap et al. 2009). However, less is known for DNA viruses—but see Domingo-Calap et al. (2009). Here, we use this approach to characterize the distribution of mutational fitness effects in the bacteriophage f1, an inovirus of the bacteriophage m13 clade, making two important improvements over previous work: first, the number of mutations tested is higher (100) and second, the contribution of experimental error to the observed distribution is explicitly accounted for. We show that one-fifth of single-nucleotide substitutions are lethal, while viable ones reduce fitness by 11% on average and can be described by a heavy-tail two-parameter distribution such as the log-normal. Interestingly, the fraction of beneficial mutations is unexpectedly high. We also compare the average effects of different mutation types and of mutations affecting different genes.     

Autophosphorylation Within the Atg1 Activation Loop Is Required for Both Kinase Activity and the Induction of Autophagy in Saccharomyces cerevisiae

Autophagy is an evolutionarily conserved degradative pathway that has been implicated in a number of physiological events important for human health. This process was originally identified as a response to nutrient deprivation and is thought to serve in a recycling capacity during periods of nutritional stress. Autophagy activity appears to be highly regulated and multiple signaling pathways are known to target a complex of proteins that contains the Atg1 protein kinase. The data here extend these observations and identify a particular phosphorylation event on Atg1 as a potential control point within the autophagy pathway in Saccharomyces cerevisiae. This phosphorylation occurs at a threonine residue, T226, within the Atg1 activation loop that is conserved in all Atg1 orthologs. Replacing this threonine with a nonphosphorylatable residue resulted in a loss of Atg1 protein kinase activity and a failure to induce autophagy. This phosphorylation required the presence of a functional Atg1 kinase domain and two known regulators of Atg1 activity, Atg13 and Atg17. Interestingly, the levels of this modification were found to increase dramatically upon exposure to conditions that induce autophagy. In addition, T226 phosphorylation was associated with an autophosphorylated form of Atg1 that was found specifically in cells undergoing the autophagy process. In all, these data suggest that autophosphorylation within the Atg1 activation loop may represent a point of regulatory control for this degradative process.MACROAUTOPHAGY (hereafter referred to as autophagy) is a highly conserved process of self-degradation that is essential for cell survival during periods of nutrient limitation (Tsukada and Ohsumi 1993). During autophagy, a double membrane grows out from a specific nucleation site, known as the pre-autophagosomal structure, or PAS, in Saccharomyces cerevisiae and the phagophore assembly site in mammals (Suzuki and Ohsumi 2007). This membrane encapsulates bulk protein and other constituents of the cytoplasm and ultimately targets this material to the vacuole/lysosome for degradation (Xie and Klionsky 2007). Recent studies have linked this pathway to a number of processes important for human health, including tumor suppression, innate immunity, and neurological disorders, like Huntington''s disease (Rubinsztein et al. 2007; Levine and Kroemer 2008). Determining how this pathway is regulated is therefore important for our understanding of these processes and our attempts to manipulate autophagy in clinically beneficial ways.Most of the molecular components of the autophagy pathway were initially characterized in the budding yeast, S. cerevisiae, but orthologs of many of these Atg proteins have since been found in other eukaryotes (Tsukada and Ohsumi 1993; Meijer et al. 2007). A complex of proteins that contains the Atg1 protein kinase is of special interest and appears to be a key point of regulatory control within this pathway (Kamada et al. 2000; Budovskaya et al. 2005; He and Klionsky 2009; Stephan et al. 2009). In S. cerevisiae, genetic and biochemical data indicate that this complex is targeted by at least three different signaling pathways. Two of these pathways, involving the Tor and cAMP-dependent protein kinases, inhibit this process, whereas the AMP-activated protein kinase is needed for the full induction of autophagy (Noda and Ohsumi 1998; Wang et al. 2001; Budovskaya et al. 2004; Stephan and Herman 2006; Kamada et al. 2010). The manner in which these signaling pathways regulate Atg1 activity and the precise role of this kinase in the autophagy process are presently matters of intense scrutiny.Although Atg1 kinase activity is required for the induction of autophagy, relatively little is known about how this enzyme is regulated in vivo. Two proteins associated with Atg1, Atg13 and Atg17, have been shown to be required for full Atg1 kinase activity both in vitro and in vivo (Kamada et al. 2000; Stephan et al. 2009). The roles of these proteins appear to be conserved through evolution as functional homologs of both have been identified in fruit flies and/or mammals (Hara et al. 2008; Chan et al. 2009; Chang and Neufeld 2009; Ganley et al. 2009; Hosokawa et al. 2009; Jung et al. 2009; Mercer et al. 2009). However, it is not yet clear precisely how these proteins stimulate Atg1 activity. In this study, we show that Atg1 is autophosphorylated within the activation loop and that this phosphorylation is required for both Atg1 kinase activity and the induction of autophagy. The activation loop is a structurally conserved element within the kinase domain and phosphorylation within this loop is often a necessary prerequisite for efficient substrate binding and/or phosphotransfer in the catalytic site (Johnson et al. 1996; Nolen et al. 2004). This loop generally corresponds to the sequence between two signature elements within the core kinase domain, the DFG and APE motifs (Hanks and Hunter 1995). Phosphorylation within this loop tends to result in a more ordered structure for this region and the proper positioning of key elements within the catalytic core of the kinase domain (Knighton et al. 1991; Johnson and O''reilly 1996; Huse and Kuriyan 2002). We found that Atg1 activation loop phosphorylation was correlated with the onset of autophagy and that replacing the site of phosphorylation with a phosphomimetic residue led to constitutive Atg1 autophosphorylation in vivo. In all, the data here suggest that Atg1 phosphorylation within its activation loop may be an important point of regulation within the autophagy pathway and models that discuss these data are presented.     

Parallel Genetic and Proteomic Screens Identify Msps as a CLASP–Abl Pathway Interactor in Drosophila

Regulation of cytoskeletal structure and dynamics is essential for multiple aspects of cellular behavior, yet there is much to learn about the molecular machinery underlying the coordination between the cytoskeleton and its effector systems. One group of proteins that regulate microtubule behavior and its interaction with other cellular components, such as actin-regulatory proteins and transport machinery, is the plus-end tracking proteins (MT+TIPs). In particular, evidence suggests that the MT+TIP, CLASP, may play a pivotal role in the coordination of microtubules with other cellular structures in multiple contexts, although the molecular mechanism by which it functions is still largely unknown. To gain deeper insight into the functional partners of CLASP, we conducted parallel genetic and proteome-wide screens for CLASP interactors in Drosophila melanogaster. We identified 36 genetic modifiers and 179 candidate physical interactors, including 13 that were identified in both data sets. Grouping interactors according to functional classifications revealed several categories, including cytoskeletal components, signaling proteins, and translation/RNA regulators. We focused our initial investigation on the MT+TIP Minispindles (Msps), identified among the cytoskeletal effectors in both genetic and proteomic screens. Here, we report that Msps is a strong modifier of CLASP and Abl in the retina. Moreover, we show that Msps functions during axon guidance and antagonizes both CLASP and Abl activity. Our data suggest a model in which CLASP and Msps converge in an antagonistic balance in the Abl signaling pathway.COORDINATION of cytoskeletal dynamics is an essential process for the regulation of virtually all aspects of cellular behavior including cell shape changes, cell division, and cell motility (e.g., Rodriguez et al. 2003; Kodama et al. 2004). Not only is the cytoskeletal system coordinated with numerous cellular pathways to control cell behavior, it also functions as a central organizing scaffold for multiple effector protein complexes downstream of signaling pathways and cellular processes such as intracellular transport. Yet, little is known regarding the molecular machinery that governs the integrated coordination of various cytoskeletal components. Evidence suggests that the microtubule (MT) plus-end tracking protein CLASP [cytoplasmic linker protein (CLIP)-associated protein], which has been implicated in mitotic spindle formation (Inoue et al. 2000, 2004) and in linking MT ends to other cell structures such as the cell cortex and kinetochore (Akhmanova et al. 2001; Maiato et al. 2003; Mimori-Kiyosue et al. 2005; Reis et al. 2009), may play a pivotal role in the overall coordination of cytoskeletal networks. Not only does it affect MT dynamics, but CLASP may function as an actin-MT crosslinker as well, as it possesses actin-binding activity (Tsvetkov et al. 2007) and CLASP-bound microtubules appear to track along F-actin bundles in growth cones (Lee et al. 2004).We previously showed that CLASP functions downstream of Abelson (Abl) nonreceptor tyrosine kinase (Lee et al. 2004), which is a key signaling molecule that modulates the cytoskeleton downstream of numerous cell surface receptor inputs and plays essential roles in various contexts including cell motility and human disease (Van Etten 1999; Moresco and Koleske 2003; Bradley and Koleske 2009). While most cytoskeletal-related studies of Abl have focused on its regulation of actin dynamics (Lanier and Gertler 2000; Wills et al. 2002; Hernandez et al. 2004; Bradley and Koleske 2009), few studies have examined its MT effectors such as CLASP and how they may be involved in the coordination of both cytoskeletal networks. Additionally, we previously reported that Drosophila CLASP is necessary for accurate embryonic axon guidance at the central nervous system (CNS) midline where conserved guidance factors (Netrins and Slits) control growth cone navigation (Lee et al. 2004). In embryonic axons, we found that CLASP is required for Abl function (Lee et al. 2004); however, additional partner proteins required for CLASP activity during axon guidance are largely unknown.Recent studies of CLASP function focusing on cell culture and imaging have identified several CLASP-binding proteins, such as additional MT+TIPs (CLIP-170 and EB1) (Akhmanova et al. 2001; Mimori-Kiyosue et al. 2005) and cell cortex-associated proteins (LL5beta and ELKS) (Lansbergen et al. 2006). While investigation of the detailed interactions and functional significance of these types of molecules and their relation to CLASP has been important to understanding the CLASP molecular mechanism, elucidating how CLASP functions in a broader context will require expanding our awareness of the entire CLASP network, or “interactome.” As of yet, there has not been a comprehensive unbiased survey of CLASP functional interactors.A long history of molecular pathway dissection in multiple model systems and biological contexts has shown that genetic and proteomic interactome screens are powerful tools for defining the network of functional partners for any given gene of interest (Xu et al. 1990; Simon et al. 1991; Carthew et al. 1994; Karim et al. 1996; Rebay et al. 2000; St Johnston 2002). Determining the CLASP interaction network can not only define MT+TIP-associated proteins and regulators of MT biology, but it can also reveal new classes of molecules that interact with CLASP. For example, a recent study suggested that CLASPs function as actin-MT crosslinkers because they possess actin-binding activity (Tsvetkov et al. 2007), but few specific actin-binding CLASP interactors have been identified, and the functional relevance of CLASP–actin interaction is still unclear. A systematic approach to define the CLASP interactome has the potential for significantly increasing our ability to understand the CLASP mechanism.Therefore, to expand our knowledge of the CLASP functional mechanism, we have performed a multilevel genetic and proteomic screen for CLASP interactors in Drosophila. The single Drosophila CLASP ortholog has been given multiple names [orbit/multiple asters (MAST)/chromosome bows (chb)] (Fedorova et al. 1997; Inoue et al. 2000; Lemos et al. 2000), but here, we refer to it as CLASP. Our screen has identified novel and specific partners in several functional categories including cytoskeletal components, signaling proteins, and the unanticipated class of translation/RNA regulators. To validate the findings of this screen, we focused our initial investigation on the conserved MT+TIP identified among the cytoskeletal effectors in both the genetic and proteomic screens, Minispindles (Msps, ortholog of the human CKAP5 [cytoskeleton associated protein 5)/TOG (tumor overexpressed gene)/Xenopus Xmap215)]. Msps function and regulation of MT stability has been studied previously in the context of the mitotic spindle (Cullen et al. 1999; Lee et al. 2001; Barros et al. 2005) and in centrosomes (Popov et al. 2002; Cassimeris and Morabito 2004), although it has not been shown to functionally interact with CLASP nor play any role in the nervous system.Here, we report that Msps is an in vivo antagonist of CLASP and interacts strongly with Abl. Furthermore, we show that Msps functions during axon guidance. Our data suggest a model in which CLASP and Msps act antagonistically to provide the growth cone with a rapidly adaptable output for Abl-dependent responses to attractive and repulsive guidance cues.     

Distinct Regulation of Mlh1p Heterodimers in Meiosis and Mitosis in Saccharomyces cerevisiae

Mlh1p forms three heterodimers that are important for mismatch repair (Mlh1p/Pms1p), crossing over during meiosis (Mlh1p/Mlh3p), and channeling crossover events into a specific pathway (Mlh1p/Mlh2p). All four proteins contain highly conserved ATPase domains and Pms1p has endonuclease activity. Studies of the functional requirements for Mlh1p/Pms1p in Saccharomyces cerevisae revealed an asymmetric contribution of the ATPase domains to repairing mismatches. Here we investigate the functional requirements of the Mlh1p and Mlh3p ATPase domains in meiosis by constructing separation of function mutations in Mlh3p. These mutations are analogous to mutations of Mlh1p that have been shown to lead to loss of ATP binding and/or ATP hydrolysis. Our data suggest that ATP binding by Mlh3p is required for meiotic crossing over while ATP hydrolysis is dispensable. This has been seen previously for Mlh1p. However, when mutations that affect ATP hydrolysis by both Mlh3p and Mlh1p are combined within a single cell, meiotic crossover frequencies are reduced. These observations suggest that the function of the Mlh1p/Mlh3p heterodimer requires both subunits to bind ATP but only one to efficiently hydrolyze it. Additionally, two different amino acid substitutions to the same residue (G97) in Mlh3p affect the minor mismatch repair function of Mlh3p while only one of them compromises its ability to promote crossing over. These studies thus reveal different functional requirements among the heterodimers formed by Mlh1p.CROSSING over during meiosis not only generates variation but is also important for providing the necessary interactions between homologous chromosomes that ensure correct segregation at division I of meiosis. Recombination is initiated by the production of programmed double-strand breaks (DSBs), catalyzed by the covalently attached Spo11p (Bergerat et al. 1997; Keeney et al. 1997), aided by a number of proteins (reviewed in Keeney and Neale 2006). DSBs are made at a much higher frequency than crossovers, and designation of only a subset to yield crossovers is thought to occur during early stages of DSB repair (Borner et al. 2004). At least two distinct pathways contribute to the production of crossover events in Saccharomyces cerevisiae. The major pathway is dependent on Msh4p/Msh5p and the mismatch repair proteins Mlh1p and Mlh3p (Ross-MacDonald and Roeder 1994; Hollingsworth et al. 1995; Hunter and Borts 1997; Wang et al. 1999; Abdullah et al. 2004) and the second pathway is dependent on Mus81p/Mms4p endonuclease (de los Santos et al. 2001, 2003).Mitotic mismatch repair (MMR) is the process by which mutations that arise during DNA replication and recombination are recognized and removed (reviewed in Kolodner 1996; Harfe and Jinks-Robertson 2000). Msh2p forms a heterodimer with Msh6p (MutSα) to repair base–base mismatches and small insertions and/or deletions and with Msh3p (MutSβ) to repair large insertions and/or deletions (reviewed in Jiricny 2006). Mlh1p forms heterodimers with Pms1p, Mlh2p, and Mlh3p to coordinate the removal of these mismatches (Prolla et al. 1994; Wang et al. 1999). Mlh1p/Pms1p (MutLα) are involved in the repair of all types of mismatches in combination with MutSα and MutSβ, and in the absence of either protein a mutator phenotype is observed (Habraken et al. 1997, 1998). Mlh1p/Mlh2p (MutLβ) and Mlh1p/Mlh3p (MutLγ) are involved in the MutSβ pathway only, which repairs frameshift mutations caused by insertions or deletions. Consequently mlh3Δ mutants only exhibit a weak mutator phenotype, due to a lesser involvement in mismatch repair and a partial overlap in function with Pms1p (Flores-Rozas and Kolodner 1998; Harfe et al. 2000).Although the MutL homologs interact primarily through their C-terminal domains (Pang et al. 1997; Ban and Yang 1998), it is thought that the N-terminal domains must also interact for the complex to be fully functional (Ban and Yang 1998). Binding of ATP causes the proteins to undergo conformational changes, which are essential for the interaction between the N termini (Ban et al. 1999; Tran and Liskay 2000; Sacho et al. 2008). ATP hydrolysis and subsequent release of ADP is required to allow the protein complex to return to its initial state, completing the cycle so that the subunits are ready to bind ATP again if required. Using mutants of MLH1 and PMS1 that are presumed to be defective for ATP binding and/or ATP hydrolysis, it has been shown that both of these functions are essential for fully effective mismatch repair (Tran and Liskay 2000). However, the ATP binding and ATP hydrolysis mutants of PMS1 exhibited lower mitotic mutation rates than the corresponding MLH1 ATPase mutants, suggesting that there is functional asymmetry within the Mlh1p/Pms1p heterodimer (Tran and Liskay 2000; Hall et al. 2002). Another example of the asymmetry in the contributions of these subunits to function can be seen in assays that measure recombination between diverged sequences (homeologous recombination). The Mlh1p ATPase activity has been shown to be more important for the suppression of homeologous recombination than Pms1p ATPase activity (Welz-Voegele et al. 2002). This functional asymmetry is supported by in vitro biochemical analysis that demonstrated Pms1p has a lower ATP binding affinity than Mlh1p (Hall et al. 2002).As mentioned above, Mlh1p/Mlh3p function in the Msh4p/Msh5p pathway for meiotic recombination (Hunter and Borts 1997; Santucci-Darmanin et al. 2000). The Msh4p/Msh5p complex is thought to act in the stabilization of Holliday junction intermediates to allow their resolution in a crossover configuration (Snowden et al. 2004). The Mlh1p/Mlh3p complex has been suggested to act in the resolution of these structures, either directly or indirectly. Human Pms2 and its yeast homolog, Pms1p, have been shown to possess a latent endonuclease activity, conferred by a motif that is conserved among some of the MutL homologs, including Mlh3p (Kadyrov et al. 2006, 2007). Mutations in the DHQA(X)₂E(X)₄E motif in yeast MLH3 cause defects in both mismatch repair and meiotic recombination equivalent to mlh3Δ, suggesting that Mlh3p may also possess an endonuclease activity that is important for the generation of crossovers (Nishant et al. 2008).ATP binding by Mlh1p has been shown to be important for both of its meiotic functions (crossing over and repair of heteroduplex DNA) (Pang et al. 1997; Tran and Liskay 2000; Hoffmann et al. 2003). In contrast, the ATP hydrolysis mutant mlh1-E31A/mlh1-E31A appears to have no effect on meiotic recombination (Tran and Liskay 2000; Hoffmann et al. 2003). This may partly be explained by in vitro studies demonstrating that this mutant exhibits a low level of ATPase activity (Hall et al. 2002).The meiotic functions of MLH1 can be functionally separated as shown by mutating the same residue, G98, to different amino acids (Hoffmann et al. 2003). The residue G98 is situated in the ATPase motif in the GFRGEAL box (GYRGDAL in Mlh3p), which forms the lid of the ATP binding pocket. Mutations in this motif are predicted to affect ATP binding and/or heterodimerization with Pms1p (Ban and Yang 1998; Ban et al. 1999). Mutating the residue G98 in the ATP binding lid to alanine resulted in defective repair of heteroduplex DNA while crossing over was unaffected, but when the same residue was mutated to valine both mismatch repair and crossover functions were defective (Hoffmann et al. 2003). The mlh1-G98V mutant disrupts the interaction of Mlh1p with Pms1p, while mlh1-G98A does not (Pang et al. 1997). This may contribute to the difference observed in the effect on crossing over as Mlh1p is thought to interact with Pms1p and Mlh3p through the same residues (Wang et al. 1999; Kondo et al. 2001). Consequently if the interaction with Pms1p is affected then it is likely that the interaction with Mlh3p is also disrupted.We constructed mlh3 mutants corresponding to the ATP binding and ATP hydrolysis mutants of mlh1 to explore the role of Mlh3p in meiotic recombination. We also constructed mlh3-G97A and mlh3-G97V mutants, equivalent to the mlh1-G98A/V pair that has been shown to differentially affect the mitotic and meiotic functions of Mlh1p. All mutants were assayed for mitotic mismatch repair, meiotic heteroduplex repair, crossing over, and chromosome segregation.     

Using RNA Interference to Identify Specific Modifiers of a Temperature-Sensitive,Embryonic-Lethal Mutation in the Caenorhabditis elegans Ubiquitin-Like Nedd8 Protein Modification Pathway E1-Activating Gene rfl-1

The essential Caenorhabditis elegans gene rfl-1 encodes one subunit of a heterodimeric E1-activating enzyme in the Nedd8 ubiquitin-like protein conjugation pathway. This pathway modifies the Cullin scaffolds of E3 ubiquitin ligases with a single Nedd8 moiety to promote ligase function. To identify genes that influence neddylation, we used a synthetic screen to identify genes that, when depleted with RNAi, enhance or suppress the embryonic lethality caused by or198ts, a temperature-sensitive (ts) mutation in rfl-1. We identified reproducible suppressor and enhancer genes and employed a systematic specificity analysis for each modifier using four unrelated ts embryonic lethal mutants. Results of this analysis highlight the importance of specificity controls in identifying genetic interactions relevant to a particular biological process because 8/14 enhancers and 7/21 suppressors modified lethality in other mutants. Depletion of the strongest specific suppressors rescued the early embryonic cell division defects in rfl-1(or198ts) mutants. RNAi knockdown of some specific suppressors partially restored Cullin neddylation in rfl-1(or198ts) mutants, consistent with their gene products normally opposing neddylation, and GFP fusions to several suppressors were detected in the cytoplasm or the nucleus, similar in pattern to Nedd8 conjugation pathway components in early embryonic cells. In contrast, depletion of the two strongest specific enhancers did not affect the early embryonic cell division defects observed in rfl-1(or198ts) mutants, suggesting that they may act at later times in other essential processes. Many of the specific modifiers are conserved in other organisms, and most are nonessential. Thus, when controlled properly for specificity, modifier screens using conditionally lethal C. elegans mutants can identify roles for nonessential but conserved genes in essential processes.UBIQUITIN-mediated proteolysis regulates many biological processes (Nandi et al. 2006). In the early Caenorhabditis elegans embryo, these include oocyte maturation, cell cycle progression, cell polarization, and cell fate patterning, all of which require the timely destruction of maternally expressed proteins (Bowerman and Kurz 2006; Greenstein and Lee 2006). One C. elegans protein targeted for proteolysis early in embryogenesis is MEI-1, the AAA-ATPase subunit of the microtubule-severing complex called katanin (Mains et al. 1990; Dow and Mains 1998; Srayko et al. 2000; Kurz et al. 2002; Pintard et al. 2003a; Xu et al. 2003). Katanin is a heterodimer of two subunits called p60 and p80 in vertebrates and MEI-1 and MEI-2 in C. elegans. Katanin in C. elegans is required for proper assembly and function of the small, barrel-shaped meiotic spindles (Albertson and Thomson 1993; McNally et al. 2006) and must be degraded after meiotic divisions to permit assembly of the much larger first mitotic spindle in the one-cell zygote. In mutants that fail to degrade katanin after the completion of meiosis, the first mitotic spindle is fragmented and mis-oriented, cytokinesis is defective, and the embryos die without hatching (Dow and Mains 1998; Srayko et al. 2000; Kurz et al. 2002).The katanin subunit MEI-1 is targeted for poly-ubiquitylation and proteolytic destruction by a Cullin-based E3 ligase (Kurz et al. 2002). This complex includes the Cullin scaffolding protein CUL-3 and a substrate-specific adaptor called MEL-26 that binds to CUL-3 through a BTB domain and to MEI-1 through a MATH domain (Pintard et al. 2003b). Cullin 3-based E3 ligases in mammals also utilize substrate-specific adaptor proteins that, like MEL-26, have both a Cullin-binding BTB/POZ domain and another protein–protein interaction domain that binds to the substrate (Geyer et al. 2003; Cullinan et al. 2004; Angers et al. 2006). While MEI-1/Katanin downregulation by the CUL-3/MEL-26 E3 ligase is essential at most growth temperatures, a mel-26 null mutation is viable at the low growth temperature of 15° (Lu and Mains 2007). This bypass of mel-26 at 15° depends at least in part on the anaphase-promoting complex and its targeting of MEI-1 for proteolytic degradation (Lu and Mains 2007). Phosphorylation by the kinase MBK-2 primes MEI-1 for proteolysis (Quintin et al. 2003; Stitzel et al. 2007) and also promotes the downregulation of MEI-1 by the anaphase-promoting complex (Lu and Mains 2007).CUL-3 is the only C. elegans Cullin thus far identified that requires modification by the ubiquitin-like protein Nedd8 (Bowerman and Kurz 2006). In contrast, C. elegans CUL-2 is required for progression through meiosis and for the localized degradation of cell fate determinants in one-cell-stage embryos (Liu et al. 2004; Sonneville and Gonczy 2004), but neddylation-defective mutants do not exhibit these early defects (Bowerman and Kurz 2006). Cullin neddylation is mediated by the Nedd8 protein conjugation pathway, which begins with a heterodimeric E1-activating enzyme consisting of ULA-1 and RFL-1 (Uba3p in budding yeast) and also includes the E2-conjugating enzyme UBC-12 (Jones and Candido 2000; Srayko et al. 2000; Kurz et al. 2002) and the E3 ligase DCN-1 (Kurz et al. 2005).The downregulation of MEI-1/katanin by the CUL-3/MEL-26 E3 ligase requires a balance of both CUL-3 neddylation, which is mediated by the Nedd8 conjugation pathway, and deneddylation, which is mediated by the conserved COP-9 Signalosome (Pintard et al. 2003a). Other Cullin-based E3 ubiquitin ligases also require a balance of neddylation and deneddylation (Lyapina et al. 2001; Schwechheimer et al. 2001; Bornstein et al. 2006; Hetfeld et al. 2008). Deneddylation may modulate activation of the E3 ligase and thereby prevent the premature degradation of substrate adaptor proteins that also can become poly-ubiquitylated and degraded as a result of E3 ligase function.To identify additional factors that influence neddylation, and the downregulation of MEI-1/katanin after the completion of meiosis in C. elegans, we report here our use of RNA interference (RNAi) to reduce gene functions in a temperature-sensitive (ts) neddylation-defective mutant, rfl-1(or198ts). The discovery of RNAi and its systemic properties in C. elegans have made it possible to systematically target C. elegans genes for depletion by feeding worms bacterial strains that express double-strand RNAs corresponding to C. elegans gene sequences (Fire et al. 1998; Timmons et al. 2001; Feinberg and Hunter 2003; Baugh et al. 2005; Lehner et al. 2006; van Haaften et al. 2006). Furthermore, chemical mutagenesis screens have identified temperature-sensitive mutations in many essential C. elegans genes, which can be used for synthetic screens by choosing intermediate-growth temperatures that sensitize the genetic background and also optimize visual scoring of embryonic viability. Recently, genomewide RNAi screens have been used to identify C. elegans genes that, when reduced in function, restore viability to temperature-sensitive, embryonic-lethal mutants (Labbe et al. 2006; O''Rourke et al. 2007). Because a loss of suppressor function restores mutant viability, the suppressors may negatively regulate either the wild-type gene product or the process that requires the wild-type gene product.Here we report our identification of C. elegans genes that, when reduced in function by feeding RNAi, reproducibly suppressed or enhanced rfl-1(or198ts) embryonic lethality. Most suppressors were specific for rfl-1(or198ts), while specific enhancement was less common. Many of the rfl-1-specific suppressors and enhancers are conserved but appear nonessential. GFP fusions to several specific suppressors exhibit localization patterns that resemble those known for neddylation pathway components, and depletion of some of these partially restored CUL-3 neddylation in rfl-1(or198ts) mutants. In addition to identifying possible roles for conserved genes in cullin neddylation, we report the first quantitative analysis of specificity for both the enhancement and the suppression of a conditionally lethal mutant in C. elegans. Our results highlight the importance of testing genetic modifiers of conditionally lethal mutants for locus specificity.