Genetic Modifiers of dFMR1 Encode RNA Granule Components in Drosophila

Mechanisms of neuronal mRNA localization and translation are of considerable biological interest. Spatially regulated mRNA translation contributes to cell-fate decisions and axon guidance during development, as well as to long-term synaptic plasticity in adulthood. The Fragile-X Mental Retardation protein (FMRP/dFMR1) is one of the best-studied neuronal translational control molecules and here we describe the identification and early characterization of proteins likely to function in the dFMR1 pathway. Induction of the dFMR1 in sevenless-expressing cells of the Drosophila eye causes a disorganized (rough) eye through a mechanism that requires residues necessary for dFMR1/FMRP''s translational repressor function. Several mutations in dco, orb2, pAbp, rm62, and smD3 genes dominantly suppress the sev-dfmr1 rough-eye phenotype, suggesting that they are required for dFMR1-mediated processes. The encoded proteins localize to dFMR1-containing neuronal mRNPs in neurites of cultured neurons, and/or have an effect on dendritic branching predicted for bona fide neuronal translational repressors. Genetic mosaic analyses indicate that dco, orb2, rm62, smD3, and dfmr1 are dispensable for translational repression of hid, a microRNA target gene, known to be repressed in wing discs by the bantam miRNA. Thus, the encoded proteins may function as miRNA- and/or mRNA-specific translational regulators in vivo.THE subcellular localization and regulated translation of stored mRNAs contributes to cellular asymmetry and subcellular specialization (Lecuyer et al. 2007; Martin and Ephrussi 2009). In mature neurons, local protein synthesis at active synapses may contribute to synapse-specific plasticity that underlies persistent forms of memory (Casadio et al. 1999; Ashraf et al. 2006; Sutton and Schuman 2006; Richter and Klann 2009). During this process, synaptic activity causes local translation of mRNAs normally stored in translationally repressed synaptic mRNPs (Sutton and Schuman 2006; Richter and Klann 2009). While specific neuronal translational repressors and microRNAs have been implicated in this process, their involvement in local translation that underlies memory, as well as the underlying mechanisms, are generally not well understood (Schratt et al. 2006; Keleman et al. 2007; Kwak et al. 2008; Li et al. 2008; Richter and Klann 2009). Furthermore, it remains possible that there are neuron-specific, mRNA-specific, and stimulus-pattern specific pathways for neuronal translational control (Raab-Graham et al. 2006; Giorgi et al. 2007).The Fragile-X Mental Retardation protein (FMRP) is among the best studied of neuronal translational repressors, in part due to its association with human neurodevelopmental disease (Pieretti et al. 1991; Mazroui et al. 2002; Gao 2008). Consistent with function in synaptic translation required for memory formation, mutations in FMRP are associated with increased synaptic translation, enhanced LTD, increased synapse growth, and also with enhanced long-term memory (Zhang et al. 2001; Huber et al. 2002; Bolduc et al. 2008; Dictenberg et al. 2008).FMRP co-immunoprecipitates with components of the RNAi and miRNA machinery and appears to be required for aspects of miRNA function in neurons (Caudy et al. 2002; Ishizuka et al. 2002; Jin et al. 2004b; Gao 2008). In addition, FMRP associates with neuronal polyribosomes as well as with Staufen-containing ribonucleoprotein (mRNP) granules easily observed in neurites of cultured neurons (Feng et al. 1997; Krichevsky and Kosik 2001; Mazroui et al. 2002; Kanai et al. 2004; Barbee et al. 2006; Bramham and Wells 2007; Bassell and Warren 2008; Dictenberg et al. 2008). FMRP-containing neuronal mRNPs contain not only several ubiquitous translational control molecules, but also CaMKII and Arc mRNAs, whose translation is locally controlled at synapses (Rook et al. 2000; Krichevsky and Kosik 2001; Kanai et al. 2004; Barbee et al. 2006). Thus, FMRP-containing RNA particles are probably translationally repressed and transported along microtubules from the neuronal cell body to synaptic sites in dendrites where local synaptic activity can induce their translation (Kiebler and Bassell 2006; Dictenberg et al. 2008).The functions of FMRP/dFMR1 in mRNA localization as well as miRNA-dependent and independent forms of translational control is likely to require several other regulatory proteins. To identify such proteins, we used a previously designed and validated genetic screen (Wan et al. 2000; Jin et al. 2004a; Zarnescu et al. 2005). The overexpression of dFMR1 in the fly eye causes a “rough-eye” phenotype through a mechanism that requires (a) key residues in dFMR1 that mediate translational repression in vitro; (b) Ago1, a known components of the miRNA pathway; and (c) a DEAD-box helicase called Me31B, which is a highly conserved protein from yeast (Dhh1p) to humans (Rck54/DDX6) functioning in translational repression and present on neuritic mRNPs (Wan et al. 2000; Laggerbauer et al. 2001; Jin et al. 2004a; Coller and Parker 2005; Barbee et al. 2006; Chu and Rana 2006). To identify other Me31B-like translational repressors and neuronal granule components, we screened mutations in 43 candidate proteins for their ability to modify dFMR1 induced rough-eye phenotype. We describe the results of this genetic screen and follow up experiments to address the potential cellular functions of five genes identified as suppressors of sev-dfmr1.     

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.     

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.     

Profiling Sex-Specific piRNAs in Zebrafish

Piwi proteins and their partner small RNAs play an essential role in fertility, germ-line stem cell development, and the basic control and evolution of animal genomes. However, little knowledge exists regarding piRNA biogenesis. Utilizing microfluidic chip analysis, we present a quantitative profile of zebrafish piRNAs expressed differentially between testis and ovary. The sex-specific piRNAs are derived from separate loci of repeat elements in the genome. Ovarian piRNAs can be categorized into groups that reach up to 92 members, indicating a sex-specific arrangement of piRNA genes in the genome. Furthermore, precursor piRNAs preferentially form a hairpin structure at the 3′end, which seem to favor the generation of mature sex-specific piRNAs. In addition, the mature piRNAs from both the testis and the ovary are 2′-O-methylated at their 3′ ends.SMALL RNAs, ranging from 19 to 30 nucleotides (nt) in length, constitute a large family of regulatory molecules with diverse functions in invertebrates, vertebrates, plants, and fungi (Bartel 2004; Nakayashiki 2005). Two major classes of small RNAs are microRNAs (miRNAs) and small interfering RNAs (siRNAs). The functions of small RNAs have been conserved through evolution; they have been shown to inhibit gene expression at the levels of mRNA degradation, translational repression, chromatin modification, heterochromatin formation, and DNA elimination (Mochizuki et al. 2002; Bartel 2004; Kim et al. 2005; Brodersen and Voinnet 2006; Lee and Collins 2006; Vaucheret 2006).Over the past few years, focus on the genetics of small RNAs has helped clarify the mechanisms behind the regulation of these molecules. While hundreds of small RNAs have been identified from mammalian somatic tissues, relatively little is known about small RNAs in germ cells. A recent breakthrough has been the identification of small RNAs that associate with Piwi proteins (piRNAs) from Drosophila and mammalian gonads (Aravin et al. 2001, 2006; Girard et al. 2006; Grivna et al. 2006; Vagin et al. 2006; Watanabe et al. 2006). piRNAs and their interacting proteins Ziwi/Zili have also been identified in zebrafish (Houwing et al. 2007, 2008). Increasing evidence indicates that piRNAs play roles mainly in germ cell differentiation and genomic stability (Carthew 2006; Lau et al. 2006; Vagin et al. 2006; Brennecke et al. 2007; Chambeyron et al. 2008; Klattenhoff and Theurkauf 2008; Kuramochi-Miyagawa et al. 2008; Kim et al. 2009; Lim et al. 2009; Unhavaithaya et al. 2009). Moreover, although piRNAs are mostly expressed in germ line cells, recent studies showed piRNA expression in nongerm cells, for example, T-cell lines (Jurkat cells and MT4) (Azuma-Mukai et al. 2008; Yeung et al. 2009), indicating other functions such as in the immune system. piRNAs do not appear to be derived from double-stranded RNA precursors, and their biogenesis mechanisms, although unclear, may be distinct from those of siRNA and miRNA. Recently, two distinct piRNA production pathways were further proposed: the “ping-pong” model (Brennecke et al. 2007; Gunawardane et al. 2007) and the Ago3-independent piRNA pathway centered on Piwi in somatic cells (Li et al. 2009; Malone et al. 2009). However, the mechanistic pathways of piRNA activity and their biogenesis are still largely unknown.Teleost fishes comprise >24,000 species, accounting for more than half of extant vertebrate species, displaying remarkable variation in morphological and physiological adaptations (see review in Zhou et al. 2001). Recently, Houwing et al. (2007, 2008) reported findings on Ziwi/Zili and associated piRNAs, implicating roles in germ cell differentiation, meiosis, and transposon silencing in the germline of the zebrafish. However, some of the identified zebrafish piRNAs are nonrepetitive and nontransposon-related piRNAs, suggesting that piRNAs may have additional unknown roles. In this study, we show that for males and females, piRNAs are specifically derived from separate loci of the repeat elements, and that ovarian piRNAs are far more often associated in groups. Genomic analysis of piRNAs indicates a tendency to folding at the 3′ end of the piRNA precursor, which may favor cleavage of the piRNA precursor to generate mature sex-specific piRNAs. Furthermore, methylation modification occurs at the 2′-O-hydroxyl group on the ribose of the final 3′ nucleotide in both the testis and the ovary.     

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.     

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.     

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?     

Estimating Selection Intensity on Synonymous Codon Usage in a Nonequilibrium Population

Codon usage bias is the nonrandom use of synonymous codons for the same amino acid. Most population genetic models of codon usage evolution assume that the population is at mutation–selection–drift equilibrium. Natural populations, however, frequently deviate from equilibrium, often because of recent demographic changes. Here, we construct a matrix model that includes the effects of a recent change in population size on estimates of selection on preferred vs. unpreferred codons. Our results suggest that patterns of synonymous polymorphisms affecting codon usage can be quite erratic after such a change; statistical methods that fail to take demographic effects into account can then give incorrect estimates of important parameters. We propose a new method that can accurately estimate both demographic and codon usage parameters. The method also provides a simple way of testing for the effects of covariates such as gene length and level of gene expression on the intensity of selection, which we apply to a large Drosophila melanogaster polymorphism data set. Our analyses of twofold degenerate codons reveal that (i) selection acts in favor of preferred codons, (ii) there is mutational bias in favor of unpreferred codons, (iii) shorter genes and genes with higher expression levels are under stronger selection, and (iv) there is little evidence for a recent change in population size in the Zimbabwe population of D. melanogaster.CODONS specifying the same amino acid are called synonymous codons. These are often used nonrandomly, with some codons appearing more frequently than others. This biased usage of synonymous codons has been found in many organisms such as Drosophila, yeast, and bacteria (Ikemura 1985; Duret and Mouchiroud 1999; Hershberg and Petrov 2008). Conventionally, synonymous codons for a given amino acid are divided into two classes: preferred and unpreferred codons (Ikemura 1985; Akashi 1994; Duret and Mouchiroud 1999). Several observations indicate that codon usage is affected by natural selection. First, in species with codon usage bias, preferred codons generally correspond to the most abundant tRNA species (Ikemura 1981). Second, highly expressed genes usually have higher codon usage bias than genes with low expression (Sharp and Li 1986; Duret and Mouchiroud 1999; Hey and Kliman 2002). Third, the synonymous substitution rate of a gene has been shown to be negatively correlated with its degree of codon usage bias (Sharp and Li 1986; Bierne and Eyre-Walker 2006). The most commonly cited explanations of the apparent fitness differences between preferred and unpreferred codons are selection for translation efficiency, translational accuracy, and mRNA stability (Ikemura 1985; Eyre-Walker and Bulmer 1993; Akashi 1994; Drummond et al. 2005). Recently, it has been proposed that exon splicing also affects codon usage bias (Warnecke and Hurst 2007).From a population genetics perspective, the extent of codon usage bias is ultimately a product of the joint effects of mutation, selection, genetic drift, recombination, and demographic history. The Li–Bulmer model of drift, selection, and reversible mutation between preferred and unpreferred codons at a site is the most widely used model (Li 1987; Bulmer 1991; McVean and Charlesworth 1999). Applications of this model generally assume that the population is at mutation–selection–drift equilibrium. However, empirical studies have suggested that changes in the strengths of various driving forces may not be unusual. For example, in Drosophila melanogaster, there is evidence that the population size (Li and Stephan 2006; Thornton and Andolfatto 2006; Keightley and Eyre-Walker 2007; Stephan and Li 2007), recombinational landscape (Takano-Shimizu 1999), and mutational process (Takano-Shimizu 2001; Kern and Begun 2005) may have changed significantly over the species'' evolutionary history.Such changes cause departures from equilibrium. Theoretical models show that it takes a very long time, proportional to the reciprocal of the mutation rate, for the population to approach a new equilibrium state (Tachida 2000; Comeron and Kreitman 2002). Before reaching equilibrium, the population often shows counterintuitive patterns of evolution (Eyre-Walker 1997; Takano-Shimizu 1999, 2001; Comeron and Kreitman 2002; Comeron and Guthrie 2005; Charlesworth and Eyre-Walker 2007). Despite these theoretical results, details of the patterns of polymorphism and substitution rates following a recent change in population size, and their effects on estimates of strength of selection, have not been determined.The above findings point to the importance of incorporating nonequilibrium factors into the study of codon usage bias. To this end, we extend the Li–Bulmer model to allow population size to vary over time, by representing the evolutionary process by a transition matrix. By analyzing this matrix model, we show that a recent change in population size can result in erratic patterns of codon usage and that methods failing to take into account these demographic effects can give false estimates of the intensity of selection.To solve these problems, we propose a new method, which does not require polarizing ancestral vs. derived states using outgroup data (cf. Cutter and Charlesworth 2006), but requires only knowledge of preferred vs. unpreferred states defined by patterns of codon usage. We use information on both polymorphic and fixed sites, which enables both mutational bias and the strength of selection to be estimated, in contrast to previous methods that use information on polymorphisms alone. Simulations indicate that this method can accurately estimate both demographic and codon usage parameters and can distinguish between selection and demography. We use the new method to analyze a large D. melanogaster polymorphism data set (Shapiro et al. 2007) and find evidence for natural selection on synonymous codons. We use our approach to show that genes with shorter coding sequences and higher levels of expression are under significantly stronger selection than longer genes with lower expression.