SUMMARYThe majority of agronomically important crop traits are quantitative, meaning that they are controlled by multiple genes each with a small effect (quantitative trait loci, QTLs). Mapping and isolation of QTLs is important for efficient crop breeding by marker-assisted selection (MAS) and for a better understanding of the molecular mechanisms underlying the traits. However, since it requires the development and selection of DNA markers for linkage analysis, QTL analysis has been time-consuming and labor-intensive. Here we report the rapid identification of plant QTLs by whole-genome resequencing of DNAs from two populations each composed of 20-50 individuals showing extreme opposite trait values for a given phenotype in a segregating progeny. We propose to name this approach QTL-seq as applied to plant species. We applied QTL-seq to rice recombinant inbred lines and F 2 populations and successfully identified QTLs for important agronomic traits, such as partial resistance to the fungal rice blast disease and seedling vigor. Simulation study showed that QTL-seq is able to detect QTLs over wide ranges of experimental variables, and the method can be generally applied in population genomics studies to rapidly identify genomic regions that underwent artificial or natural selective sweeps.

3 1 l e t t e r sThe shift from outcrossing to selfing is common in flowering plants 1,2 , but the genomic consequences and the speed at which they emerge remain poorly understood. An excellent model for understanding the evolution of self fertilization is provided by Capsella rubella, which became self compatible <200,000 years ago. We report a C. rubella reference genome sequence and compare RNA expression and polymorphism patterns between C. rubella and its outcrossing progenitor Capsella grandiflora. We found a clear shift in the expression of genes associated with flowering phenotypes, similar to that seen in Arabidopsis, in which self fertilization evolved about 1 million years ago. Comparisons of the two Capsella species showed evidence of rapid genome-wide relaxation of purifying selection in C. rubella without a concomitant change in transposable element abundance. Overall we document that the transition to selfing may be typified by parallel shifts in gene expression, along with a measurable reduction of purifying selection.

DNA methylation of coding regions, known as gene body methylation, is conserved across eukaryotic lineages. The function of body methylation is not known, but it may either prevent aberrant expression from intragenic promoters or enhance the accuracy of splicing. Given these putative functions, we hypothesized that body-methylated genes would be both longer and more functionally important than unmethylated genes. To test these hypotheses, we reanalyzed single-base resolution bisulfite sequence data from Arabidopsis thaliana to differentiate body-methylated genes from unmethylated genes using a probabilistic approach. Contrasting genic characteristics between the two groups, we found that body-methylated genes tend to be longer and to be more functionally important, as measured by phenotypic effects of insertional mutants and by gene expression, than unmethylated genes. We also found that methylated genes evolve more slowly than unmethylated genes, despite the potential for increased mutation rates in methylated CpG dinucleotides. We propose that slower rates in body-methylated genes are a function of higher selective constraint, lower nucleosome occupancy, and a lower proportion of CpG dinucleotides.

DNA methylation is a common feature of eukaryotic genomes and is especially common in noncoding regions of plants. Protein coding regions of plants are often methylated also, but the extent, function, and evolutionary consequences of gene body methylation remain unclear. Here we investigate gene body methylation using an explicit comparative evolutionary approach. We generated bisulfite sequencing data from two tissues of Brachypodium distachyon and compared genic methylation patterns to those of rice (Oryza sativa ssp. japonica). Gene body methylation was strongly conserved between orthologs of the two species and affected a biased subset of long, slowly evolving genes. Because gene body methylation is conserved over evolutionary time, it shapes important features of plant genome evolution, such as the bimodality of G+C content among grass genes. Our results superficially contradict previous observations of high cytosine methylation polymorphism within Arabidopsis thaliana genes, but reanalyses of these data are consistent with conservation of methylation within gene regions. Overall, our results indicate that the methylation level is a long-term property of individual genes and therefore of evolutionary consequence.epigenetics | methylome | Poaceae | molecular evolution

Little is known about patterns of genic DNA methylation across the plant kingdom or about the evolutionary processes that shape them. To characterize gene-body methylation (gbM) within exons, we have gathered single-base resolution methylome data that span the phylogenetic breadth of land plants. We find that a basal land plant, Marchantia polymorpha, lacks any evident signal of gbM within exons, but conifers have high levels of both CG and CHG (where H is A, C or T) methylation in expressed genes. To begin to understand the evolutionary forces that shape gbM, we first tested for correlations in methylation levels across orthologues(1,2). Genic CG methylation levels, but not CHG or CHH levels, are correlated across orthologues for species as distantly related as ferns and angiosperms. Hence, relative levels of CG methylation are a consistent property across genes, even for species that diverged ∼400 million years ago(3,4). In contrast, genic CHG methylation correlates with genome size, suggesting that the host epigenetic response to transposable elements also affects genes. Altogether, our data indicate that the evolutionary forces acting on DNA methylation vary substantially across species, genes and methylation contexts.

Early studies of protein evolution established that evolutionary rates vary among sites, genes, and lineages. Here we review rate variation within each these three components, with an emphasis on plant nuclear genes. Evolutionary rates vary among nucleotide sites as a consequence of selection and mutational biases. Mutation varies among sites as a function of the base, its neighboring bases, and its proximity to insertion-deletion polymorphisms. Nonsynonymous rate variation among genes correlates most strongly with gene expression, perhaps owing to selection for translational robustness. Among lineages, perennial plants evolve more slowly than annuals, but the mechanism driving this effect remains unclear. There are also interactions among these three main components; we discuss the patterns of interaction and their significance for understanding the evolutionary forces that shape nucleotide substitution rates in plants. Finally, we highlight gaps in our knowledge and future opportunities to answer fundamental questions in plant molecular evolution.

The semidwarf phenotype has been extensively selected during modern crop breeding as an agronomically important trait. Introduction of the semidwarf gene, semi-dwarf1 (sd1), which encodes a gibberellin biosynthesis enzyme, made significant contributions to the "green revolution" in rice (Oryza sativa L.). Here we report that SD1 was involved not only in modern breeding including the green revolution, but also in early steps of rice domestication. We identified two SNPs in O. sativa subspecies (ssp.) japonica SD1 as functional nucleotide polymorphisms (FNPs) responsible for shorter culm length and low gibberellin biosynthetic activity. Genetic diversity analysis among O. sativa ssp. japonica and indica, along with their wild ancestor O. rufipogon Griff, revealed that these FNPs clearly differentiate the japonica landrace and O. rufipogon. We also found a dramatic reduction in nucleotide diversity around SD1 only in the japonica landrace, not in the indica landrace or O. rufipogon. These findings indicate that SD1 has been subjected to artificial selection in rice evolution and that the FNPs participated in japonica domestication, suggesting that ancient humans already used the green revolution gene.

Convergent evolution is the independent evolution of similar traits in different species or lineages of the same species; this often is a result of adaptation to similar environments, a process referred to as convergent adaptation. We investigate here the molecular basis of convergent adaptation in maize to highland climates in Mesoamerica and South America, using genome-wide SNP data. Taking advantage of archaeological data on the arrival of maize to the highlands, we infer demographic models for both populations, identifying evidence of a strong bottleneck and rapid expansion in South America. We use these models to then identify loci showing an excess of differentiation as a means of identifying putative targets of natural selection and compare our results to expectations from recently developed theory on convergent adaptation. Consistent with predictions across a wide parameter space, we see limited evidence for convergent evolution at the nucleotide level in spite of strong similarities in overall phenotypes. Instead, we show that selection appears to have predominantly acted on standing genetic variation and that introgression from wild teosinte populations appears to have played a role in highland adaptation in Mexican maize.KEYWORDS convergent evolution; highland adaptation; maize; population genomics C ONVERGENT evolution occurs when multiple species or populations exhibit similar phenotypic adaptations to comparable environmental challenges (Wood et al. 2005;Arendt and Reznick 2008;Elmer and Meyer 2011). Evolutionary genetic analysis of a wide range of species has provided evidence for multiple pathways that lead to convergent evolution. One such route occurs when identical mutations arise independently and fix via natural selection in multiple populations. In humans, for example, malaria resistance due to mutations from Glu to Val at the sixth codon of the b-globin gene has arisen independently on multiple unique haplotypes (Currat et al. 2002;Kwiatkowski 2005). Convergent evolution can also be achieved when different mutations arise within the same locus yet produce similar phenotypic effects. Grain fragrance in rice appears to have evolved along these lines, as populations across East Asia have similar fragrances resulting from at least eight distinct loss-of-function alleles in the BADH2 gene (Kovach et al. 2009). Finally, convergent evolution may arise from natural selection acting on standing genetic variation in an ancestral population. In the three-spined stickleback, natural selection has repeatedly acted to reduce armor plating in independent colonizations of freshwater environments. Adaptation in these populations occurred both from new mutations and from standing variation at the Eda locus in marine populations (Colosimo et al. 2005).Not all convergent phenotypic evolution is the result of convergent evolution at the molecular level, however. Recent studies of adaptation to high elevation in humans, for example, reveal that the genes involved in highland adaptation are largely dist...