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.     

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.     

C-Terminal Flap Endonuclease (rad27) Mutations: Lethal Interactions With a DNA Ligase I Mutation (cdc9-p) and Suppression by Proliferating Cell Nuclear Antigen (POL30) in Saccharomyces cerevisiae

During lagging-strand DNA replication in eukaryotic cells primers are removed from Okazaki fragments by the flap endonuclease and DNA ligase I joins nascent fragments. Both enzymes are brought to the replication fork by the sliding clamp proliferating cell nuclear antigen (PCNA). To understand the relationship among these three components, we have carried out a synthetic lethal screen with cdc9-p, a DNA ligase mutation with two substitutions (F43A/F44A) in its PCNA interaction domain. We recovered the flap endonuclease mutation rad27-K325* with a stop codon at residue 325. We created two additional rad27 alleles, rad27-A358* with a stop codon at residue 358 and rad27-pX8 with substitutions of all eight residues of the PCNA interaction domain. rad27-pX8 is temperature lethal and rad27-A358* grows slowly in combination with cdc9-p. Tests of mutation avoidance, DNA repair, and compatibility with DNA repair mutations showed that rad27-K325* confers severe phenotypes similar to rad27Δ, rad27-A358* confers mild phenotypes, and rad27-pX8 confers phenotypes intermediate between the other two alleles. High-copy expression of POL30 (PCNA) suppresses the canavanine mutation rate of all the rad27 alleles, including rad27Δ. These studies show the importance of the C terminus of the flap endonuclease in DNA replication and repair and, by virtue of the initial screen, show that this portion of the enzyme helps coordinate the entry of DNA ligase during Okazaki fragment maturation.CELLULAR maintenance of genomic integrity is essential for the continued viability of all organisms. The fidelity of DNA replication has to be maintained and DNA insults have to be repaired to ensure that deleterious mutations are not passed on to progeny or cause cancerous growth. A number of cellular proteins have multiple roles in DNA replication, mutation avoidance, and repair. In Saccharomyces cerevisiae, the flap endonuclease, proliferating cell nuclear antigen (PCNA), and DNA ligase I encoded by RAD27, POL30, and CDC9, respectively, are all required for proper replication and also function to avoid mutation and to facilitate repair.The flap endonuclease, FEN-1 in humans, is a highly conserved structure-specific nuclease that has both endonuclease and 5′–3′ exonuclease activity. During lagging-strand replication these activities function to remove primers from Okazaki fragments, either by endonucleolytic cleavage of a flap made by strand displacement (Liu et al. 2004) or by sequential exonucleolytic removal of single nucleotides at the 5′ end of the primer (Murante et al. 1994).While deletion of RAD27 is not lethal to yeast cells, the rad27Δ mutant exhibits temperature-sensitive growth, is a mutator, and undergoes genomic instability (Johnson et al. 1995; Reagan et al. 1995; Tishkoff et al. 1997b; Chen and Kolodner 1999). In addition, its sensitivity to low doses of the methylating agent methylmethane sulfonate (MMS) implicates the participation of the enzyme in base excision repair (BER) (Reagan et al. 1995; Wu and Wang 1999). rad27Δ mutants have been reported to be either mildly sensitive to UV light or not sensitive to UV light (Reagan et al. 1995; Sommers et al. 1995). In the strain background that the mutant is mildly sensitive, its combination with rad2Δ yields a double mutant more sensitive than each single mutant, implying that the enzyme does not participate in RAD2-mediated nucleotide excision repair (NER) (Reagan et al. 1995). The flap endonuclease has also been implicated in double-strand break (DSB) repair by virtue of the incompatibility of rad27Δ with mutations of the DSB repair pathways (Tishkoff et al. 1997b; Symington 1998). In addition, either the yeast enzyme or its human ortholog has been shown to participate in reactions of homologous recombination, nonhomologous end joining, and telomere maintenance (Parenteau and Wellinger 1999, 2002; Wu et al. 1999; Wang et al. 2004; Kikuchi et al. 2005). Curiously, the rad27Δ mutant is not sensitive to gamma radiation but is sensitive to high doses of MMS that are thought to act as a radiomimetic agent (Reagan et al. 1995; Sommers et al. 1995).PCNA is the replicative clamp that acts as a scaffold to facilitate the loading of DNA replication and repair proteins, including DNA ligase I and the flap endonuclease to DNA (Warbrick 2000, 2006; Maga and Hubscher 2003). PCNA (POL30) is essential for cell viability, which is indicative of its central role in DNA metabolism. Biochemical characterization of its effect on the flap endonuclease shows that it stimulates its activity ∼50-fold, evidencing the productive nature of the interaction (Gomes and Burgers 2000; Tom et al. 2000; Frank et al. 2001; Stucki et al. 2001). The ability of DNA ligase to efficiently catalyze the formation of phosphodiester bonds in the DNA backbone may also be facilitated by its binding to PCNA. Tom et al. (2001) showed that, in vitro, PCNA enhances the ligation reaction 5-fold and that the stable association of DNA ligase with nicked duplex DNA requires PCNA.Both DNA ligase and the flap endonuclease bind to PCNA via their respective PCNA interactive peptide domains (PIP box). The PIP box is a conserved sequence motif of the amino acids QXXLXXFF. The PIP box fits into the interdomain connector loop (IDCL) of PCNA to provide a protein–protein interaction surface (Gomes and Burgers 2000; Chapados et al. 2004; Sakurai et al. 2005; Pascal et al. 2006). Mutations in the PIP box or the IDCL that impair the interaction of DNA ligase and the flap endonuclease to PCNA lead to genomic instability (Amin and Holm 1996; Eissenberg et al. 1997; Gary et al. 1999; Refsland and Livingston 2005; Subramanian et al. 2005). We have reported that the double mutants made by combinations of cdc9-p, rad27-p, and pol30-90—mutations with alterations of the PIP box or the IDCL in the respective proteins—have synergistic phenotypes with respect to MMS sensitivity and to trinucleotide repeat instability (Refsland and Livingston 2005). These results suggest that the two enzymes function in a concerted manner that is facilitated by PCNA.The precise nature of how PCNA coordinates the entry of the flap endonuclease and DNA ligase into the replication fork is not well understood. Biochemical and structural studies have begun to elucidate a possible ordering of these PCNA-mediated interactions. The possibility of such an ordering is underscored by the observation that DNA ligase adopts a toroidal conformation by completely encircling duplex DNA while interacting with PCNA (Pascal et al. 2004). Moreover, both PCNA and DNA ligase may be loaded onto the DNA in a mechanism utilizing the replication clamp loader replication factor C (RFC) (Levin et al. 2004; Vijayakumar et al. 2009), again suggesting a complete encirclement of the DNA by DNA ligase as well as by PCNA. PCNA and DNA ligase are similar in size and their interaction is likely to extend along the face of PCNA in a manner that would prevent other proteins such as the flap endonuclease from binding to the IDCL (Pascal et al. 2004, 2006). A biochemical study with purified yeast proteins showed that the two enzymes cannot bind simultaneously to PCNA (Subramanian et al. 2005). These studies suggest that a coordinated sequential interaction among PCNA, DNA ligase, and the flap endonuclease is important for replication and repair.Alternatively, both the flap endonuclease and DNA ligase may bind to the same molecule of PCNA. Since PCNA is a homotrimer, DNA ligase can potentially bind to one monomer while the flap endonuclease binds to another, using its extended C-terminal tail in a conformation allowing it to be tethered to PCNA concurrently with DNA ligase (Gomes and Burgers 2000; Sakurai et al. 2005). DNA ligase could also bind to PCNA in an extended conformation while the flap endonuclease cleaves the DNA. Sulfolobus solfataricus DNA ligase has been shown to have an open, extended conformation while binding to PCNA (Pascal et al. 2006). Presumably, once the flap endonuclease has removed the 5′ flap, DNA ligase acquires a closed, ring-shaped conformation to catalyze the joining of Okazaki fragments (Pascal et al. 2006).Exactly how the interaction of these enzymes with PCNA is coordinated in vivo, whether singly or concurrently, is not well understood. To further elucidate how the interaction of DNA ligase with PCNA is ordered, we performed a genetic screen to identify mutations that are synthetically lethal with cdc9-p (F44A/F35A), an allele of DNA ligase that has impaired binding to PCNA (Refsland and Livingston 2005; Subramanian et al. 2005). We postulated that genes recovered from this screen would function in DNA repair, replication, and recombination or would be involved in ordering the DNA ligase–PCNA interaction. From the screen we recovered a truncated allele of RAD27, rad27-K325*. This allele encodes a protein that lacks the PIP box and the entire C-terminal domain of the enzyme but retains the N terminus containing the nuclease activities. We have characterized this allele and compared it to two other rad27 alleles in which we have created different alterations of the C-terminal end of the flap endonuclease.     

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.     

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?     

Sexual Isolation in Acinetobacter baylyi Is Locus-Specific and Varies 10,000-Fold Over the Genome

Naturally transformable bacteria acquire chromosomal DNA from related species at lower frequencies than from cognate DNA sources. To determine how genome location affects heterogamic transformation in bacteria, we inserted an nptI marker into random chromosome locations in 19 different strains of the Acinetobacter genus (>24% divergent at the mutS/trpE loci). DNA from a total of 95 nptI-tagged isolates was used to transform the recipient Acinetobacter baylyi strain ADP1. A total of >1300 transformation assays revealed that at least one nptI-tagged isolate for each of the strains/species tested resulted in detectable integration of the nptI marker into the ADP1 genome. Transformation frequencies varied up to ∼10,000-fold among independent nptI insertions within a strain. The location and local sequence divergence of the nptI flanking regions were determined in the transformants. Heterogamic transformation depended on RecA and was hampered by DNA mismatch repair. Our studies suggest that single-locus-based studies, and inference of transfer frequencies from general estimates of genomic sequence divergence, is insufficient to predict the recombination potential of chromosomal DNA fragments between more divergent genomes. Interspecies differences in overall gene content, and conflicts in local gene organization and synteny are likely important determinants of the genomewide variation in recombination rates between bacterial species.HORIZONTAL gene transfer (HGT) contributes to bacterial evolution by providing access to DNA evolved and retained in separate species or strains (Cohan 1994a,b; Bergstrom et al. 2000; Ochman et al. 2000; Feil et al. 2001; Koonin 2003; Lawrence and Hendrickson 2003; Fraser et al. 2007). Multilocus sequence typing (MLST) has provided strong evidence for frequent transfer and recombination of chromosomal DNA between related bacterial strains within the same species (Maiden et al. 1998; Enright et al. 2002). HGT occurring by natural transformation allows bacteria to exploit the presence of nucleic acids in their environment for the purposes of nutrition, DNA repair, reacquisition of lost genes, and/or acquisition of novel genetic diversity (Redfield 1993; Mehr and Seifert 1998; Dubnau 1999; Claverys et al. 2000; Szöllösi et al. 2006; Johnsen et al. 2009). It can be inferred from observations of the presence of extracellular DNA in most environments that bacteria are constantly exposed to DNA from a variety of sources, without such exposure necessarily producing observable changes in the genetic compositions of bacterial populations over evolutionary time (Thomas and Nielsen 2005; Nielsen et al. 2007a,b).The absence of sequence similarity between the donor DNA and the DNA of the recipient bacterium is the strongest barrier to the horizontal acquisition of chromosomal genes in bacteria (Matic et al. 1996; Vulic et al. 1997; Majewski 2001; Townsend et al. 2003) as illegitimate recombination occurs only at extremely low frequencies in bacteria (Hülter and Wackernagel 2008a). Single-locus transfer models have been extensively applied and have demonstrated a log-linear decrease in recombination frequencies with increasing sequence divergence for Bacillus subtilis (Roberts and Cohan, 1993; Zawadzki et al. 1995), Acinetobacter baylyi (Young and Ornston 2001), Escherichia coli (Shen and Huang 1986; Vulic et al. 1997), and Streptococcus pneumoniae (Majewski et al. 2000). For instance, heterogamic transformation between nonmutator isolates at the rpoB locus of B. mojavensis is undetectable at sequence divergences >16.7% (Zawadzki et al. 1995) and between S. pneumoniae isolates with sequence divergences >18% (Majewski et al. 2000). In A. baylyi, the nonmutator sequence divergence limit for detectable transformation at the pcaH locus of strain ADP1 was found to be 20% (Young and Ornston 2001), and up to 24% overall divergence yielded transformants at 16S rRNA loci in strain DSM587 (Strätz et al. 1996).Several recent studies also show that short stretches (<200 bp) of DNA sequence identity can facilitate additive or substitutive integration of longer stretches (>1000 bp) of heterologous DNA in bacteria (Prudhomme et al. 1991, 2002; de Vries and Wackernagel 2002; Hülter and Wackernagel 2008a). Thus, the uptake of DNA in bacteria can facilitate larger substitutions within gene sequences and the integration of additional DNA material on the basis of recombination initiated in flanking DNA stretches (either at one or both ends) with high sequence similarity (Nielsen et al. 2000). On the other hand, segments of heterologous DNA interrupting the synteny of homologous DNA have also been shown to be a barrier in intraspecies transformation in S. pneumoniae (Pasta and Sicard 1996, 1999).The various studies of the interspecies transfer potential of single genes demonstrate that the immediate local sequence divergence of the transferred locus is of high importance in determining recombination frequencies in hosts up to 20% divergent (at the housekeeping gene level). However, it can be hypothesized that the broader structural, organizational, and biochemical properties of the genome region surrounding a particular locus will determine its transfer potential to more divergent host species (Cohan 2001; Lawrence 2002). The interspecies transfer potential of various genome regions/loci between more diverged species (>20% at the housekeeping gene level) may therefore differ substantially from a log-linear model (determined experimentally for more closely related species) as local gene organization becomes less conserved with evolutionary time. The barriers to gene exchange between divergent bacterial species is likely a combination of inefficient recombination due to both mismatched base pairs (the main determinator in the log-linear model) and conflicting gene order and organization across the local recombining DNA regions. In addition, selective barriers due to negative effects on host fitness of the transferred DNA regions may become increasingly important for the removal of recombination events from the bacterial population. Recent bioinformatics-based genome analysis of E. coli and Salmonella genomes suggests various parts of the bacterial genome may have different suceptibilities to undergo evolutionarily successful recombination leading to temporal fragmentation of speciation (Lawrence 2002; Retchless and Lawrence 2007). Nevertheless, few studies have experimentally tested the effect of variable species and chromosome locations of genes on their transfer potential between bacteria (Ravin and Chen 1967; Ravin and Chakrabarti 1975; Siddiqui and Goldberg 1975; Cohan et al. 1991; Huang et al. 1991; Fall et al. 2007).Here, we determine to what extent genome location contributes to sexual isolation between the recipient A. baylyi strain ADP1 and 19 sequence divergent (24–27% divergent at the mutS/trpE loci) donor Acinetobacter strains and species (carrying a selectable nptI gene in a total of 95 random genome locations).