Jonathan M. Flowers

Molecular evidence for a single evolutionary origin of domesticated rice

Asian rice, Oryza sativa, is one of worlds oldest and most important crop species. Rice is believed to have been domesticated ∼9,000 y ago, although debate on its origin remains contentious. A single-origin model suggests that two main subspecies of Asian rice, indica and japonica, were domesticated from the wild rice O. rufipogon. In contrast, the multiple independent domestication model proposes that these two major rice types were domesticated separately and in different parts of the species range of wild rice. This latter view has gained much support from the observation of strong genetic differentiation between indica and japonica as well as several phylogenetic studies of rice domestication. We reexamine the evolutionary history of domesticated rice by resequencing 630 gene fragments on chromosomes 8, 10, and 12 from a diverse set of wild and domesticated rice accessions. Using patterns of SNPs, we identify 20 putative selective sweeps on these chromosomes in cultivated rice. Demographic modeling based on these SNP data and a diffusion-based approach provide the strongest support for a single domestication origin of rice. Bayesian phylogenetic analyses implementing the multispecies coalescent and using previously published phylogenetic sequence datasets also point to a single origin of Asian domesticated rice. Finally, we date the origin of domestication at ∼8,200–13,500 y ago, depending on the molecular clock estimate that is used, which is consistent with known archaeological data that suggests rice was first cultivated at around this time in the Yangtze Valley of China.

A Conserved Mutation in an Ethylene Biosynthesis Enzyme Leads to Andromonoecy in Melons

Andromonoecy is a widespread sexual system in angiosperms characterized by plants carrying both male and bisexual flowers. In melon, this sexual form is controlled by the identity of the alleles at the andromonoecious (a) locus. Cloning of the a gene reveals that andromonoecy results from a mutation in the active site of 1-aminocyclopropane-1-carboxylic acid synthase. Expression of the active enzyme inhibits the development of the male organs and is not required for carpel development. A causal single-nucleotide polymorphism associated with andromonoecy was identified, which suggests that the a allele has been under recent positive selection and may be linked to the evolution of this sexual system.

Population Genomics of the Arabidopsis thaliana Flowering Time Gene Network

The time to flowering is a key component of the life-history strategy of the model plant Arabidopsis thaliana that varies quantitatively among genotypes. A significant problem for evolutionary and ecological genetics is to understand how natural selection may operate on this ecologically significant trait. Here, we conduct a population genomic study of resequencing data from 52 genes in the flowering time network. McDonald-Kreitman tests of neutrality suggested a strong excess of amino acid polymorphism when pooling across loci. This excess of replacement polymorphism across the flowering time network and a skewed derived frequency spectrum toward rare alleles for both replacement and noncoding polymorphisms relative to synonymous changes is consistent with a large class of deleterious polymorphisms segregating in these genes. Assuming selective neutrality of synonymous changes, we estimate that approximately 30% of amino acid polymorphisms are deleterious. Evidence of adaptive substitution is less prominent in our analysis. The photoperiod regulatory gene, CO, and a gibberellic acid transcription factor, AtMYB33, show evidence of adaptive fixation of amino acid mutations. A test for extended haplotypes revealed no examples of flowering time alleles with haplotypes comparable in length to those associated with the null fri(Col) allele reported previously. This suggests that the FRI gene likely has a uniquely intense or recent history of selection among the flowering time genes considered here. Although there is some evidence for adaptive evolution in these life-history genes, it appears that slightly deleterious polymorphisms are a major component of natural molecular variation in the flowering time network of A. thaliana.

Domestication history and geographical adaptation inferred from a SNP map of African rice

African rice (Oryza glaberrima Steud.) is a cereal crop species closely related to Asian rice (Oryza sativa L.) but was independently domesticated in West Africa ∼3,000 years ago. African rice is rarely grown outside sub-Saharan Africa but is of global interest because of its tolerance to abiotic stresses. Here we describe a map of 2.32 million SNPs of African rice from whole-genome resequencing of 93 landraces. Population genomic analysis shows a population bottleneck in this species that began ∼13,000–15,000 years ago with effective population size reaching its minimum value ∼3,500 years ago, suggesting a protracted period of population size reduction likely commencing with predomestication management and/or cultivation. Genome-wide association studies (GWAS) for six salt tolerance traits identify 11 significant loci, 4 of which are within ∼300 kb of genomic regions that possess signatures of positive selection, suggesting adaptive geographical divergence for salt tolerance in this species.

Natural selection in gene dense regions shapes the genomic pattern of polymorphism in wild and domesticated rice

Levels of nucleotide variability are frequently positively correlated with recombination rate and negatively associated with gene density due to the effects of selection on linked variation. These relationships are determined by properties that frequently differ among species, including the mating system, and aspects of genome organization such as how genes are distributed along chromosomes. In rice, genes are found at highest density in regions with frequent crossing-over. This association between gene density and recombination rate provides an opportunity to evaluate the effects of selection in a genomic context that differs from other model organisms. Using single-nucleotide polymorphism data from Asian domesticated rice Oryza sativa ssp. japonica and ssp. indica and their progenitor species O. rufipogon, we observe a significant negative association between levels of polymorphism and both gene and coding site density, but either no association, or a negative correlation, between nucleotide variability and recombination rate. We establish that these patterns are unlikely to be explained by neutral mutation rate biases and demonstrate that a model of background selection with variable rates of deleterious mutation is sufficient to account for the gene density effect in O. rufipogon. In O. sativa ssp. japonica, we report a strong negative correlation between polymorphism and recombination rate and greater losses of variation during domestication in the euchromatic chromosome arms than heterochromatin. This is consistent with Hill-Robertson interference in low-recombination regions, which may limit the efficacy of selection for domestication traits. Our results suggest that the physical distribution of selected mutations is a primary factor that determines the genomic pattern of polymorphism in wild and domesticated rice species.

Variation, Sex, and Social Cooperation: Molecular Population Genetics of the Social Amoeba Dictyostelium discoideum

Dictyostelium discoideum is a eukaryotic microbial model system for multicellular development, cell–cell signaling, and social behavior. Key models of social evolution require an understanding of genetic relationships between individuals across the genome or possibly at specific genes, but the nature of variation within D. discoideum is largely unknown. We re-sequenced 137 gene fragments in wild North American strains of D. discoideum and examined the levels and patterns of nucleotide variation in this social microbial species. We observe surprisingly low levels of nucleotide variation in D. discoideum across these strains, with a mean nucleotide diversity (π) of 0.08%, and no strong population stratification among North American strains. We also do not find any clear relationship between nucleotide divergence between strains and levels of social dominance and kin discrimination. Kin discrimination experiments, however, show that strains collected from the same location show greater ability to distinguish self from non-self than do strains from different geographic areas. This suggests that a greater ability to recognize self versus non-self may arise among strains that are more likely to encounter each other in nature, which would lead to preferential formation of fruiting bodies with clonemates and may prevent the evolution of cheating behaviors within D. discoideum populations. Finally, despite the fact that sex has rarely been observed in this species, we document a rapid decay of linkage disequilibrium between SNPs, the presence of recombinant genotypes among natural strains, and high estimates of the population recombination parameter ρ. The SNP data indicate that recombination is widespread within D. discoideum and that sex as a form of social interaction is likely to be an important aspect of the life cycle.

Phylogeography of Asian wild rice, Oryza rufipogon: a genome-wide view.

Asian wild rice (Oryza rufipogon) that ranges widely across the eastern and southern part of Asia is recognized as the direct ancestor of cultivated Asian rice (O. sativa). Studies of the geographic structure of O. rufipogon, based on chloroplast and low‐copy nuclear markers, reveal a possible phylogeographic signal of subdivision in O. rufipogon. However, this signal of geographic differentiation is not consistently observed among different markers and studies, with often conflicting results. To more precisely characterize the phylogeography of O. rufipogon populations, a genome‐wide survey of unlinked markers, intensively sampled from across the entire range of O. rufipogon is critical. In this study, we surveyed sequence variation at 42 genome‐wide sequence tagged sites (STS) in 108 O. rufipogon accessions from throughout the native range of the species. Using Bayesian clustering, principal component analysis and amova, we conclude that there are two genetically distinct O. rufipogon groups, Ruf‐I and Ruf‐II. The two groups exhibit a clinal variation pattern generally from north‐east to south‐west. Different from many earlier studies, Ruf‐I, which is found mainly in China and the Indochinese Peninsula, shows genetic similarity with one major cultivated rice variety, O. satvia indica, whereas Ruf‐II, mainly from South Asia and the Indochinese Peninsula, is not found to be closely related to cultivated rice varieties. The other major cultivated rice variety, O. sativa japonica, is not found to be similar to either O. rufipogon groups. Our results support the hypothesis of a single origin of the domesticated O. sativa in China. The possible role of palaeoclimate, introgression and migration–drift balance in creating this clinal variation pattern is also discussed.

Variation in Arabidopsis flowering time associated with cis-regulatory variation in CONSTANS.

The onset of flowering, the change from vegetative to reproductive development, is a major life history transition in flowering plants. Recent work suggests that mutations in cis-regulatory mutations should play critical roles in the evolution of this (as well as other) important adaptive traits, but thus far there has been little evidence that directly links regulatory mutations to evolutionary change at the species level. While several genes have previously been shown to affect natural variation in flowering time in Arabidopsis thaliana, most either show protein-coding changes and/or are found at low frequency (<5%). Here we identify and characterize natural variation in the cis-regulatory sequence in the transcription factor CONSTANS that underlies flowering time diversity in Arabidopsis. Mutation in this regulatory motif evolved recently and has spread to high frequency in Arabidopsis natural accessions, suggesting a role for these cis-regulatory changes in adaptive variation of flowering time.

Whole genome re-sequencing of date palms yields insights into diversification of a fruit tree crop

Date palms (Phoenix dactylifera) are the most significant perennial crop in arid regions of the Middle East and North Africa. Here, we present a comprehensive catalogue of approximately seven million single nucleotide polymorphisms in date palms based on whole genome re-sequencing of a collection of 62 cultivars. Population structure analysis indicates a major genetic divide between North Africa and the Middle East/South Asian date palms, with evidence of admixture in cultivars from Egypt and Sudan. Genome-wide scans for selection suggest at least 56 genomic regions associated with selective sweeps that may underlie geographic adaptation. We report candidate mutations for trait variation, including nonsense polymorphisms and presence/absence variation in gene content in pathways for key agronomic traits. We also identify a copia-like retrotransposon insertion polymorphism in the R2R3 myb-like orthologue of the oil palm virescens gene associated with fruit colour variation. This analysis documents patterns of post-domestication diversification and provides a genomic resource for this economically important perennial tree crop.

The evolution of plant genomes — scaling up from a population perspective

Plant genomes exhibit tremendous diversity in both their size and structure, with genome sizes across land plants ranging over two to three orders of magnitude and significant variation in structural organization was observed across species (EA Kellogg, JL Bennetzen, The evolution of nuclear genome structure in seed plants, Am J Bot 2004, 91:1709-1725). Five plant genomes are now either completely sequenced or in the draft stage; the grape (O Jaillon et al., The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla, Nature 2007, 449:463-467) and papaya (R Ming et al., The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus), Nature 2008, 452:991-997) whole genome sequences were reported most recently. Moreover, sequencing of 41 additional genomes is in progress. There is now an emerging consensus that understanding genome evolution requires consideration of the population genetics of genome diversification, and that description of evolutionary forces at the level of populations and within species can help identify the features that led to plant genome diversity (M Lynch, JS Conery, The origins of genome complexity, Science 2003, 302:1401-1404). In this review we focus on advances in our understanding of the mechanisms that drive the diversification of genomes. In particular, we look at the extent to which demographic features such as effective population size changes within species can drive genome evolution, discuss population genetic models of genome diversification associated with transposable element (TE) mobilization, and describe recent studies on the evolution of gene families.