Genome-wide patterns of nucleotide polymorphism in domesticated rice

Inferring the joint demographic history of multiple populations from multidimensional SNP frequency data

Demographic models built from genetic data play important roles in illuminating prehistorical events and serving as null models in genome scans for selection. We introduce an inference method based on the joint frequency spectrum of genetic variants within and between populations. For candidate models we numerically compute the expected spectrum using a diffusion approximation to the one-locus, two-allele Wright-Fisher process, involving up to three simultaneous populations. Our approach is a composite likelihood scheme, since linkage between neutral loci alters the variance but not the expectation of the frequency spectrum. We thus use bootstraps incorporating linkage to estimate uncertainties for parameters and significance values for hypothesis tests. Our method can also incorporate selection on single sites, predicting the joint distribution of selected alleles among populations experiencing a bevy of evolutionary forces, including expansions, contractions, migrations, and admixture. We model human expansion out of Africa and the settlement of the New World, using 5 Mb of noncoding DNA resequenced in 68 individuals from 4 populations (YRI, CHB, CEU, and MXL) by the Environmental Genome Project. We infer divergence between West African and Eurasian populations 140 thousand years ago (95% confidence interval: 40–270 kya). This is earlier than other genetic studies, in part because we incorporate migration. We estimate the European (CEU) and East Asian (CHB) divergence time to be 23 kya (95% c.i.: 17–43 kya), long after archeological evidence places modern humans in Europe. Finally, we estimate divergence between East Asians (CHB) and Mexican-Americans (MXL) of 22 kya (95% c.i.: 16.3–26.9 kya), and our analysis yields no evidence for subsequent migration. Furthermore, combining our demographic model with a previously estimated distribution of selective effects among newly arising amino acid mutations accurately predicts the frequency spectrum of nonsynonymous variants across three continental populations (YRI, CHB, CEU).

The demographic history of our species is reflected in patterns of genetic variation within and among populations. We developed an efficient method for calculating the expected distribution of genetic variation, given a demographic model including such events as population size changes, population splits and joins, and migration. We applied our approach to publicly available human sequencing data, searching for models that best reproduce the observed patterns. Our joint analysis of data from African, European, and Asian populations yielded new dates for when these populations diverged. In particular, we found that African and Eurasian populations diverged around 100,000 years ago. This is earlier than other genetic studies suggest, because our model includes the effects of migration, which we found to be important for reproducing observed patterns of variation in the data. We also analyzed data from European, Asian, and Mexican populations to model the peopling of the Americas. Here, we find no evidence for recurrent migration after East Asian and Native American populations diverged. Our methods are not limited to studying humans, and we hope that future sequencing projects will offer more insights into the history of both our own species and others.

Indel-associated mutation rate varies with mating system in flowering plants

A recently proposed mutational mechanism, indel-associated mutation (IDAM), posits that heterozygous insertions/deletions (indels) increase the point mutation rate at nearby nucleotides due to errors during meiosis. This mechanism could have especially dynamic consequences for the evolution of plant genomes, because the high degree of variation in the rate of self-fertilization among plant species causes differences in the heterozygosity of alleles, including indel alleles, segregating in plant species. In this study, we investigated the consequences of IDAM for species differing in mating system using both forward population genetic simulations and genomewide DNA resequencing data from Arabidopsis thaliana, Oryza sativa, and Oryza rufipogon. Simulations of different levels of selfing suggest that the effect of IDAM on surrounding nucleotide diversity should decrease with increasing selfing rate. Further simulations incorporating selfing rates and the time of onset of selfing suggest that the time since the switch to selfing also affects patterns of nucleotide diversity due to IDAM. Population genetic analyses of A. thaliana and Oryza DNA sequence data sets empirically confirmed our simulation results, revealing the strongest effect of IDAM in the outcrossing O. rufipogon, a weaker effect in the recently evolved selfer O. sativa, and the weakest effect in the relatively ancient selfer A. thaliana. These results support the novel idea that differences in life history, such as the level of selfing, can affect the per-individual mutation rate among species.

A population genomics study of the Arabidopsis core cell cycle genes shows the signature of natural selection

Large-scale comparison of sequence polymorphism and divergence at numerous genomic loci within and between closely related species can reveal signatures of natural selection. Here, we present a population genomics study based on direct sequencing of 61 mitotic cell cycle genes from 30 Arabidopsis thaliana accessions and comparison of the resulting data to the close relative Arabidopsis lyrata. We found that the Arabidopsis core cell cycle (CCC) machinery is not highly constrained but is subject to different modes of selection. We found patterns of purifying selection for the cyclin-dependent kinase (CDK), CDK subunit, retinoblastoma, and WEE1 gene families. Other CCC gene families often showed a mix of one or two constrained genes and relaxed purifying selection on the other genes. We found several large effect mutations in CDKB1;2 that segregate in the species. We found a strong signature of adaptive protein evolution in the Kip-related protein KRP6 and departures from equilibrium at CDKD;1 and CYCA3;3 consistent with the operation of selection in these gene regions. Our data suggest that within Arabidopsis, the genetic robustness of cell cycle-related processes is more due to functional redundancy than high selective constraint.

The origin and evolution of fragrance in rice (Oryza sativa L.)

Fragrance in the grain is one of the most highly valued grain quality traits in rice, yet the origin and evolution of the betaine aldehyde dehydrogenase gene (BADH2) underlying this trait remains unclear. In this study, we identify eight putatively nonfunctional alleles of the BADH2 gene and show that these alleles have distinct geographic and genetic origins. Despite multiple origins of the fragrance trait, a single allele, badh2.1, is the predominant allele in virtually all fragrant rice varieties today, including the widely recognized Basmati and Jasmine types. Haplotype analysis allowed us to establish a single origin of the badh2.1 allele within the Japonica varietal group and demonstrate the introgression of this allele from Japonica to Indica. Basmati-like accessions were nearly identical to the ancestral Japonica haplotype across a 5.3-Mb region flanking BADH2 regardless of their fragrance phenotype, demonstrating a close evolutionary relationship between Basmati varieties and the Japonica gene pool. These results clarify the relationships among fragrant rice varieties and challenge the traditional assumption that the fragrance trait arose in the Indica varietal group.

Evolutionary history of GS3, a gene conferring grain length in rice

Unlike maize and wheat, where artificial selection is associated with an almost uniform increase in seed or grain size, domesticated rice exhibits dramatic phenotypic diversity for grain size and shape. Here we clone and characterize GS3, an evolutionarily important gene controlling grain size in rice. We show that GS3 is highly expressed in young panicles in both short- and long-grained varieties but is not expressed in leaves or panicles after flowering, and we use genetic transformation to demonstrate that the dominant allele for short grain complements the long-grain phenotype. An association study revealed that a C to A mutation in the second exon of GS3 (A allele) was associated with enhanced grain length in Oryza sativa but was absent from other Oryza species. Linkage disequilibrium (LD) was elevated and there was a 95.7% reduction in nucleotide diversity (theta(pi)) across the gene in accessions carrying the A allele, suggesting positive selection for long grain. Haplotype analysis traced the origin of the long-grain allele to a Japonica-like ancestor and demonstrated introgression into the Indica gene pool. This study indicates a critical role for GS3 in defining the seed morphologies of modern subpopulations of O. sativa and enhances the potential for genetic manipulation of grain size in rice.

Evolutionary history of GS3, a gene conferring grain length in rice

Unlike maize and wheat, where artificial selection is associated with an almost uniform increase in seed or grain size, domesticated rice exhibits dramatic phenotypic diversity for grain size and shape. Here we clone and characterize GS3, an evolutionarily important gene controlling grain size in rice. We show that GS3 is highly expressed in young panicles in both short- and long-grained varieties but is not expressed in leaves or panicles after flowering, and we use genetic transformation to demonstrate that the dominant allele for short grain complements the long-grain phenotype. An association study revealed that a C to A mutation in the second exon of GS3 (A allele) was associated with enhanced grain length in Oryza sativa but was absent from other Oryza species. Linkage disequilibrium (LD) was elevated and there was a 95.7% reduction in nucleotide diversity (theta(pi)) across the gene in accessions carrying the A allele, suggesting positive selection for long grain. Haplotype analysis traced the origin of the long-grain allele to a Japonica-like ancestor and demonstrated introgression into the Indica gene pool. This study indicates a critical role for GS3 in defining the seed morphologies of modern subpopulations of O. sativa and enhances the potential for genetic manipulation of grain size in rice.

Patterns of DNA sequence variation at candidate gene loci in black poplar (Populus nigra L.) as revealed by single nucleotide polymorphisms

Black poplar (Populus nigra L.) is an economically and ecologically important tree species and an ideal organism for studies of genetic variation. In the present work, we use a candidate gene approach to infer the patterns of DNA variation in natural populations of this species. A total of 312 single nucleotide polymorphisms (SNPs) are found among 8,056 bp sequenced from nine drought-adaptation and photosynthesis-related gene loci. The median SNP frequency is one site per 26 bp. The average nucleotide diversity is calculated to be theta(W) = 0.01074 and pi(T) = 0.00702, higher values than those observed in P. tremula, P. trichocarpa and most conifer species. Tests of neutrality for each gene reveal a general excess of low-frequency mutations, a greater number of haplotypes than expected and an excess of high-frequency derived variants in P. nigra, which is consistent with previous findings that genetic hitchhiking has occurred in this species. Linkage disequilibrium is low, decaying rapidly from 0.45 to 0.20 or less within a distance of 300 bp, although the declines of r(2) are variable among different loci. This is similar to the rate of decay reported in most other tree species. Our dataset is expected to enhance understanding of how evolutionary forces shape genetic variation, and it will contribute to molecular breeding in black poplar.

Plant genomes do a balancing act

Balancing selection is one mechanism that may explain why diversity is maintained in wild populations. However, relatively few examples of genes showing evidence of balancing selection have been identified, particularly in plants. In this issue, Reininga et al. (2009) present three Arabidopsis loci that show strong evidence of balancing selection. The loci, discovered using a genome-scanning approach, encode proteins with diverse predicted functions: starch synthesis and control of gene expression. These three genes were identified by scanning only a small fraction of the Arabidopsis genome, suggesting that balancing selection may be more prevalent than previously known.

DNA changes tell us about rice domestication

Crop domestication can be considered a model system of plant evolution. Genome analyses of rice have revealed the fine population structure of this major crop associated with local origins of landraces. Recent cloning of rice domestication-related genes and identification of the responsible functional nucleotide polymorphisms in landraces, while taking into account their population structures, have revealed the existence of historical signatures of the DNA involved in the domestication process. These signatures imply the importance of multiple selection steps wherein natural variants were combined to improve crop performance during domestication. These analyses will provide new insights into the relationship between Darwinian selection for agronomical phenotypes and DNA changes in terms of plant evolution.

Genic non-coding microsatellites in the rice genome: characterization, marker design and use in assessing genetic and evolutionary relationships among domesticated groups

Background

Completely sequenced plant genomes provide scope for designing a large number of microsatellite markers, which are useful in various aspects of crop breeding and genetic analysis. With the objective of developing genic but non-coding microsatellite (GNMS) markers for the rice (Oryza sativa L.) genome, we characterized the frequency and relative distribution of microsatellite repeat-motifs in 18,935 predicted protein coding genes including 14,308 putative promoter sequences.

Results

We identified 19,555 perfect GNMS repeats with densities ranging from 306.7/Mb in chromosome 1 to 450/Mb in chromosome 12 with an average of 357.5 GNMS per Mb. The average microsatellite density was maximum in the 5' untranslated regions (UTRs) followed by those in introns, promoters, 3'UTRs and minimum in the coding sequences (CDS). Primers were designed for 17,966 (92%) GNMS repeats, including 4,288 (94%) hypervariable class I types, which were bin-mapped on the rice genome. The GNMS markers were most polymorphic in the intronic region (73.3%) followed by markers in the promoter region (53.3%) and least in the CDS (26.6%). The robust polymerase chain reaction (PCR) amplification efficiency and high polymorphic potential of GNMS markers over genic coding and random genomic microsatellite markers suggest their immediate use in efficient genotyping applications in rice. A set of these markers could assess genetic diversity and establish phylogenetic relationships among domesticated rice cultivar groups. We also demonstrated the usefulness of orthologous and paralogous conserved non-coding microsatellite (CNMS) markers, identified in the putative rice promoter sequences, for comparative physical mapping and understanding of evolutionary and gene regulatory complexities among rice and other members of the grass family. The divergence between long-grained aromatics and subspecies japonica was estimated to be more recent (0.004 Mya) compared to short-grained aromatics from japonica (0.006 Mya) and long-grained aromatics from subspecies indica (0.014 Mya).

Conclusion

Our analyses showed that GNMS markers with their high polymorphic potential would be preferred candidate functional markers in various marker-based applications in rice genetics, genomics and breeding. The CNMS markers provided encouraging implications for their use in comparative genome mapping and understanding of evolutionary complexities in rice and other members of grass family.

Methods for human demographic inference using haplotype patterns from genomewide single-nucleotide polymorphism data

We propose a novel approximate-likelihood method to fit demographic models to human genomewide single-nucleotide polymorphism (SNP) data. We divide the genome into windows of constant genetic map width and then tabulate the number of distinct haplotypes and the frequency of the most common haplotype for each window. We summarize the data by the genomewide joint distribution of these two statistics-termed the HCN statistic. Coalescent simulations are used to generate the expected HCN statistic for different demographic parameters. The HCN statistic provides additional information for disentangling complex demography beyond statistics based on single-SNP frequencies. Application of our method to simulated data shows it can reliably infer parameters from growth and bottleneck models, even in the presence of recombination hotspots when properly modeled. We also examined how practical problems with genomewide data sets, such as errors in the genetic map, haplotype phase uncertainty, and SNP ascertainment bias, affect our method. Several modifications of our method served to make it robust to these problems. We have applied our method to data collected by Perlegen Sciences and find evidence for a severe population size reduction in northwestern Europe starting 32,500-47,500 years ago.

SNP deserts of Asian cultivated rice: genomic regions under domestication

When performing a genome-wide comparison between indica (93-11) and japonica (Nipponbare), we find 8% of the genome, which have an extremely low SNP rate (< 1 SNP/kb). Inside these 'SNP deserts', experimentally confirmed genes show increased K(a)/K(s) that indicate adaptive selection. To further elucidate this connection, we survey the level and pattern of genetic variation in both cultivated and wild rice groups, using 155 noncoding regions located within SNP deserts. The results suggest that cultivated rice has greatly reduced genetic variation within SNP deserts as compared to either the nondesert or corresponding genomic regions in wild rice. Consistent with this reduction in genetic variation, we find a biased distribution of derived allele frequency in the cultivated group, indicative of positive selection. Furthermore, over half of the confirmed, domestication-related genes are found within SNP deserts, also suggesting that SNP deserts are strongly related to domestication, and might be the key sites in the process of domestication.

The nature of selection during plant domestication

Plant domestication is an outstanding example of plant-animal co-evolution and is a far richer model for studying evolution than is generally appreciated. There have been numerous studies to identify genes associated with domestication, and archaeological work has provided a clear understanding of the dynamics of human cultivation practices during the Neolithic period. Together, these have provided a better understanding of the selective pressures that accompany crop domestication, and they demonstrate that a synthesis from the twin vantage points of genetics and archaeology can expand our understanding of the nature of evolutionary selection that accompanies domestication.

The nature of selection during plant domestication

Plant domestication is an outstanding example of plant-animal co-evolution and is a far richer model for studying evolution than is generally appreciated. There have been numerous studies to identify genes associated with domestication, and archaeological work has provided a clear understanding of the dynamics of human cultivation practices during the Neolithic period. Together, these have provided a better understanding of the selective pressures that accompany crop domestication, and they demonstrate that a synthesis from the twin vantage points of genetics and archaeology can expand our understanding of the nature of evolutionary selection that accompanies domestication.

Linkage disequilibrium and demographic history of wild and domestic canids

Assessing the extent of linkage disequilibrium (LD) in natural populations of a nonmodel species has been difficult due to the lack of available genomic markers. However, with advances in genotyping and genome sequencing, genomic characterization of natural populations has become feasible. Using sequence data and SNP genotypes, we measured LD and modeled the demographic history of wild canid populations and domestic dog breeds. In 11 gray wolf populations and one coyote population, we find that the extent of LD as measured by the distance at which r2=0.2 extends <10 kb in outbred populations to >1.7 Mb in populations that have experienced significant founder events and bottlenecks. This large range in the extent of LD parallels that observed in 18 dog breeds where the r2 value varies from approximately 20 kb to >5 Mb. Furthermore, in modeling demographic history under a composite-likelihood framework, we find that two of five wild canid populations exhibit evidence of a historical population contraction. Five domestic dog breeds display evidence for a minor population contraction during domestication and a more severe contraction during breed formation. Only a 5% reduction in nucleotide diversity was observed as a result of domestication, whereas the loss of nucleotide diversity with breed formation averaged 35%.

Analysis of site frequency spectra from Arabidopsis with context-dependent corrections for ancestral misinference

Previous studies have shown that the pattern of single nucleotide polymorphism (SNP) in Arabidopsis (Arabidopsis thaliana) deviates from the distribution expected under a neutral model. Here, we test whether or not ancestral misinference could explain this deviation. We start by showing that there are significant and complex influences of context on mutation dynamics as inferred from SNP frequency, in Arabidopsis, and compare the results to observations about context dependency that have been made on a previous analysis of a maize (Zea mays) SNP dataset. The data concerning heterogeneity across sites are then used to make corrections for ancestral misinference in a context-dependent manner. Using Arabidopsis lyrata to infer the ancestral state for SNPs, we show that the resulting unfolded site frequency spectrum (SFS) in Arabidopsis is skewed toward sites with high frequency derived nucleotides. Sites are also partitioned into two general functional classes, second codon position and 4-fold degenerate sites. These two classes show different SFS; although both show an overrepresentation of high frequency derived sites, low frequency derived sites are vastly overrepresented at the second codon position, but significantly underrepresented at 4-fold degenerate sites. We find that these results are robust to corrections for ancestral misinference, even when context-dependent variation in mutation properties is taken into consideration. The data suggest that, in addition to purifying selection, complex demographic events and/or linked positive selection need to be invoked to explain the SFS, and they highlight the importance of sequence context in analyses of genome-wide variation.

A comparative view of the evolution of grasses under domestication

Crop grasses were among the first plants to be domesticated c. 12,000 yr ago, and they still represent the main staple crops for humans. During domestication, as did many other crops, grasses went through dramatic genetic and phenotypic changes. The recent massive increase in genomic data has provided new tools to investigate the genetic basis and consequences of domestication. Beyond the genetics of domestication, many aspects of grass biology, including their phylogeny and developmental biology, are also increasingly well studied, offering a unique opportunity to analyse the domestication process in a comparative way. Taking such a comparative point of view, we review the history of domesticated grasses and how domestication affected their phenotypic and genomic diversity. Considering recent theoretical developments and the accumulation of genetic data, we revisit more specifically the role of mating systems in the domestication process. We close by suggesting future directions for the study of domestication in grasses.

Multilocus patterns of nucleotide diversity, population structure and linkage disequilibrium in Boechera stricta, a wild relative of Arabidopsis

Information about polymorphism, population structure, and linkage disequilibrium (LD) is crucial for association studies of complex trait variation. However, most genomewide studies have focused on model systems, with very few analyses of undisturbed natural populations. Here, we sequenced 86 mapped nuclear loci for a sample of 46 genotypes of Boechera stricta and two individuals of B. holboellii, both wild relatives of Arabidopsis. Isolation by distance was significant across the species range of B. stricta, and three geographic groups were identified by structure analysis, principal coordinates analysis, and distance-based phylogeny analyses. The allele frequency spectrum indicated a genomewide deviation from an equilibrium neutral model, with silent nucleotide diversity averaging 0.004. LD decayed rapidly, declining to background levels in approximately 10 kb or less. For tightly linked SNPs separated by <1 kb, LD was dependent on the reference population. LD was lower in the specieswide sample than within populations, suggesting that low levels of LD found in inbreeding species such as B. stricta, Arabidopsis thaliana, and barley may result from broad geographic sampling that spans heterogeneous genetic groups. Finally, analyses also showed that inbreeding B. stricta and A. thaliana have approximately 45% higher recombination per kilobase than outcrossing A. lyrata.

A genomic scan for selection reveals candidates for genes involved in the evolution of cultivated sunflower (Helianthus annuus)

Genomic scans for selection are a useful tool for identifying genes underlying phenotypic transitions. In this article, we describe the results of a genome scan designed to identify candidates for genes targeted by selection during the evolution of cultivated sunflower. This work involved screening 492 loci derived from ESTs on a large panel of wild, primitive (i.e., landrace), and improved sunflower (Helianthus annuus) lines. This sampling strategy allowed us to identify candidates for selectively important genes and investigate the likely timing of selection. Thirty-six genes showed evidence of selection during either domestication or improvement based on multiple criteria, and a sequence-based test of selection on a subset of these loci confirmed this result. In view of what is known about the structure of linkage disequilibrium across the sunflower genome, these genes are themselves likely to have been targeted by selection, rather than being merely linked to the actual targets. While the selection candidates showed a broad range of putative functions, they were enriched for genes involved in amino acid synthesis and protein catabolism. Given that a similar pattern has been detected in maize (Zea mays), this finding suggests that selection on amino acid composition may be a general feature of the evolution of crop plants. In terms of genomic locations, the selection candidates were significantly clustered near quantitative trait loci (QTL) that contribute to phenotypic differences between wild and cultivated sunflower, and specific instances of QTL colocalization provide some clues as to the roles that these genes may have played during sunflower evolution.

Independent losses of function in a polyphenol oxidase in rice: differentiation in grain discoloration between subspecies and the role of positive selection under domestication

Asian rice (Oryza sativa) cultivars originated from wild rice and can be divided into two subspecies by several criteria, one of which is the phenol reaction (PHR) phenotype. Grains of indica cultivars turn brown in a phenol solution that accelerates a similar process that occurs during prolonged storage. By contrast, the grains of japonica do not discolor. This distinction may reflect the divergent domestication of these two subspecies. The PHR is controlled by a single gene, Phr1; here, we report the cloning of Phr1, which encodes a polyphenol oxidase. The Phr1 gene is indeed responsible for the PHR phenotype, as transformation with a functional Phr1 can complement a PHR negative cultivar. Phr1 is defective in all japonica lines but functional in nearly all indica and wild strains. Phylogenetic analysis showed that the defects in Phr1 arose independently three times. The multiple recent origins and rapid spread of phr1 in japonica suggest the action of positive selection, which is further supported by several population genetic tests. This case may hence represent an example of artificial selection driving the differentiation among domesticated varieties.

Patterns of molecular evolution associated with two selective sweeps in the Tb1-Dwarf8 region in maize

We focused on a region encompassing a major maize domestication locus, Tb1, and a locus involved in the flowering time variation, Dwarf8 (D8), to investigate the consequences of two closely linked selective sweeps on nucleotide variation and gain some insights into maize geographical diffusion, through climate adaptation. First, we physically mapped D8 at approximately 300 kb 3' of Tb1. Second, we analyzed patterns of nucleotide variation at Tb1, D8, and seven short regions (400-700 bp) located in the Tb1-D8 region sequenced on a 40 maize inbred lines panel encompassing early-flowering temperate and late-flowering tropical lines. The pattern of polymorphism along the region is characterized by two valleys of depleted polymorphism while the region in between exhibits an appreciable amount of diversity. Our results reveal that a region approximately 100 kb upstream of the D8 gene exhibits hallmarks of divergent selection between temperate and tropical lines and is likely closer than the D8 gene to the target of selection for climate adaptation. Selection in the tropical lines appears more recent than in the temperate lines, suggesting an initial domestication of early-flowering maize. Simulation results indicate that the polymorphism pattern is consistent with two interfering selective sweeps at Tb1 and D8.

Inference of the japonica rice domestication process from the distribution of six functional nucleotide polymorphisms of domestication-related genes in various landraces and modern cultivars

Crop domestication can serve as a model of plant evolutionary processes. It involves a series of selection events from standing natural variation and newly occurring mutations and combinations of mutations as a result of natural crossings in populations during local adaptation and propagation of plant lines to other cultivation areas. Our earlier identification of three functional nucleotide polymorphisms (FNPs) of distinct genes involved in the rice domestication process led us to propose a model of the japonica rice domestication process. Here, we examined three more FNPs in two domestication-related genes involved in pigment synthesis during the development of seed pericarp color (Rc and Rd) in 91 landraces (and some modern cultivars) of japonica rice collected from throughout the area of distribution of rice. These polymorphisms were assigned by using genome-wide patterns of restriction fragment length polymorphisms (RFLPs) and the local origins of the landraces. The results led us to infer the process of japonica rice domestication in more detail and propose a more refined model of the japonica domestication process. In this model, the critical role of the Rc FNP at an early step of the domestication process was highlighted. Independent artificial selections of two defective Rd alleles were found, suggesting a role for Rd other than in pigment synthesis during rice domestication.

The subtelomere of Oryza sativa chromosome 3 short arm as a hot bed of new gene origination in rice

Despite general observations of non-random genomic distribution of new genes, it is unclear whether or not new genes preferentially occur in certain genomic regions driven by related molecular mechanisms. Using 1.5 Mb of genomic sequences from short arms of chromosome 3 of Oryza glaberrima and O. punctata, we conducted a comparative genomic analysis with the reference O. sativa ssp. japonica genome. We identified a 60-kb segment located in the middle of the subtelomeric region of chromosome 3, which is unique to the species O. sativa. The region contained gene duplicates that occurred in Asian cultivated rice species that diverged from the ancestor of Asian and African cultivated rice one million years ago (MYA). For the 12 genes and one complete retrotransposon identified in this segment in O. sativa ssp. japonica, we searched for their parental genes. The high similarity between duplicated paralogs further supports the recent origination of these genes. We found that this segment was recently generated through multiple independent gene recombination and transposon insertion events. Among the 12 genes, we found that five had chimeric gene structures derived from multiple parental genes. Nine out of the 12 new genes seem to be functional, as suggested by Ka/Ks analysis and the presence of cDNA and/or MPSS data. Furthermore, for the eight transcribed genes, at least two genes could be classified as defense or stress response-related genes. Given these findings, and the fact that subtelomeres are associated with high rates of recombination and transcription, it is likely that subtelomeres may facilitate gene recombination and transposon insertions and serve as hot spots for new gene origination in rice genomes.

Multilocus patterns of nucleotide polymorphism and the demographic history of Populus tremula

I have studied nucleotide polymorphism and linkage disequilibrium using multilocus data from 77 fragments, with an average length of fragments of 550 bp, in the deciduous tree Populus tremula (Salicaceae). The frequency spectrum across loci showed a modest excess of mutations segregating at low frequency and a marked excess of high-frequency derived mutations at silent sites, relative to neutral expectations. These excesses were also seen at replacement sites, but were not so pronounced for high-frequency derived mutations. There was a marked excess of low-frequency mutations at replacement sites, likely indicating deleterious amino acid-changing mutations that segregate at low frequencies in P. tremula. I used approximate Bayesian computation (ABC) to evaluate a number of different demographic scenarios and to estimate parameters for the best-fitting model. The data were found to be consistent with a historical reduction in the effective population size of P. tremula through a bottleneck. The timing inferred for this bottleneck is largely consistent with geological data and with data from several other long-lived plant species. The results show that P. tremula harbors substantial levels of nucleotide polymorphism with the posterior mode of the scaled mutation rate, theta = 0.0177 across loci. The ABC analyses also provided an estimate of the scaled recombination rate that indicates that recombination rates in P. tremula are likely to be 2-10 times higher than the mutation rate. This study reinforces the notion that linkage disequilibrium is low and decays to negligible levels within a few hundred base pairs in P. tremula.

Molecular evolution of the Pi-ta gene resistant to rice blast in wild rice (Oryza rufipogon)

Rice blast disease resistance to the fungal pathogen Magnaporthe grisea is triggered by a physical interaction between the protein products of the host R (resistance) gene, Pi-ta, and the pathogen Avr (avirulence) gene, AVR-pita. The genotype variation and resistant/susceptible phenotype at the Pi-ta locus of wild rice (Oryza rufipogon), the ancestor of cultivated rice (O. sativa), was surveyed in 36 locations worldwide to study the molecular evolution and functional adaptation of the Pi-ta gene. The low nucleotide polymorphism of the Pi-ta gene of O. rufipogon was similar to that of O. sativa, but greatly differed from what has been reported for other O. rufipogon genes. The haplotypes can be subdivided into two divergent haplogroups named H1 and H2. H1 is derived from H2, with nearly no variation and at a low frequency. H2 is common and is the ancestral form. The leucine-rich repeat (LRR) domain has a high pi(non)/pi(syn) ratio, and the low polymorphism of the Pi-ta gene might have primarily been caused by recurrent selective sweep and constraint by other putative physiological functions. Meanwhile, we provide data to show that the amino acid Ala-918 of H1 in the LRR domain has a close relationship with the resistant phenotype. H1 might have recently arisen during rice domestication and may be associated with the scenario of a blast pathogen-host shift from Italian millet to rice.

TriMEDB: a database to integrate transcribed markers and facilitate genetic studies of the tribe Triticeae

BACKGROUND

The recent rapid accumulation of sequence resources of various crop species ensures an improvement in the genetics approach, including quantitative trait loci (QTL) analysis as well as the holistic population analysis and association mapping of natural variations. Because the tribe Triticeae includes important cereals such as wheat and barley, integration of information on the genetic markers in these crops should effectively accelerate map-based genetic studies on Triticeae species and lead to the discovery of key loci involved in plant productivity, which can contribute to sustainable food production. Therefore, informatics applications and a semantic knowledgebase of genome-wide markers are required for the integration of information on and further development of genetic markers in wheat and barley in order to advance conventional marker-assisted genetic analyses and population genomics of Triticeae species.

DESCRIPTION

The Triticeae mapped expressed sequence tag (EST) database (TriMEDB) provides information, along with various annotations, regarding mapped cDNA markers that are related to barley and their homologues in wheat. The current version of TriMEDB provides map-location data for barley and wheat ESTs that were retrieved from 3 published barley linkage maps (the barley single nucleotide polymorphism database of the Scottish Crop Research Institute, the barley transcript map of Leibniz Institute of Plant Genetics and Crop Plant Research, and HarvEST barley ver. 1.63) and 1 diploid wheat map. These data were imported to CMap to allow the visualization of the map positions of the ESTs and interrelationships of these ESTs with public gene models and representative cDNA sequences. The retrieved cDNA sequences corresponding to each EST marker were assigned to the rice genome to predict an exon-intron structure. Furthermore, to generate a unique set of EST markers in Triticeae plants among the public domain, 3472 markers were assembled to form 2737 unique marker groups as contigs. These contigs were applied for pairwise comparison among linkage maps obtained from different EST map resources.

CONCLUSION

TriMEDB provides information regarding transcribed genetic markers and functions as a semantic knowledgebase offering an informatics facility for the acceleration of QTL analysis and for population genetics studies of Triticeae.

Patterns of polymorphism and demographic history in natural populations of Arabidopsis lyrata

BACKGROUND

Many of the processes affecting genetic diversity act on local populations. However, studies of plant nucleotide diversity have largely ignored local sampling, making it difficult to infer the demographic history of populations and to assess the importance of local adaptation. Arabidopsis lyrata, a self-incompatible, perennial species with a circumpolar distribution, is an excellent model system in which to study the roles of demographic history and local adaptation in patterning genetic variation.

PRINCIPAL FINDINGS

We studied nucleotide diversity in six natural populations of Arabidopsis lyrata, using 77 loci sampled from 140 chromosomes. The six populations were highly differentiated, with a median FST of 0.52, and structure analysis revealed no evidence of admixed individuals. Average within-population diversity varied among populations, with the highest diversity found in a German population; this population harbors 3-fold higher levels of silent diversity than worldwide samples of A. thaliana. All A. lyrata populations also yielded positive values of Tajima's D. We estimated a demographic model for these populations, finding evidence of population divergence over the past 19,000 to 47,000 years involving non-equilibrium demographic events that reduced the effective size of most populations. Finally, we used the inferred demographic model to perform an initial test for local adaptation and identified several genes, including the flowering time gene FCA and a disease resistance locus, as candidates for local adaptation events.

CONCLUSIONS

Our results underscore the importance of population-specific, non-equilibrium demographic processes in patterning diversity within A. lyrata. Moreover, our extensive dataset provides an important resource for future molecular population genetic studies of local adaptation in A. lyrata.

Nonindependent domestication of the two rice subspecies, Oryza sativa ssp. indica and ssp. japonica, demonstrated by multilocus microsatellites

The origins of the Asian cultivated rice Oryza sativa from its wild ancestor O. rufipogon have been debated for decades. The question mainly concerns whether it originated monophyletically or polyphyletically. To shed light on the origins and demographic history of rice domestication, we genotyped a total of 92 individual plants from the two O. sativa subspecies and O. rufipogon for 60 microsatellites. An approximate Bayesian method was applied to estimate demographic parameters for O. rufipogon vs. O. sativa ssp. indica and O. rufipogon vs. O. sativa ssp. japonica. We showed that the japonica subspecies suffered a more severe bottleneck than the indica subspecies and thus a greater loss of genetic variation during its domestication. Across microsatellite loci there is a significant positive correlation in the reduction of genetic diversity between the two subspecies. The results suggest that completely independent domestication of indica and japonica subspecies may not explain our data and that there is at least partial sharing of their ancestral populations and/or recent gene flow between them.

Genomic consequences of selection on self-incompatibility genes

Frequency-dependent selection at plant self-incompatibility systems is inherent and well understood theoretically. A self-incompatibility locus leads to a strong peak of diversity in the genome, to a unique distribution of diversity across the species and possibly to increased introgression between closely related species. We review recent empirical studies demonstrating these features and relate the empirical findings to theoretical predictions. We show how these features are being exploited in searches for other genes under multi-allelic balancing selection and for inference on recent breakdown of self-incompatibility.

Leveraging natural diversity: back through the bottleneck

Plant breeders have long recognized the existence of useful genetic variation in the wild ancestors of our domesticated crop species. In cultivated rice (Oryza sativa), crosses between high-yielding elite cultivars and low-yielding wild accessions often give rise to superior offspring, with wild alleles conferring increased performance in the context of the elite cultivar genetic background. Because the breeding value of wild germplasm cannot be determined by examining the performance of wild accessions, a phylogenetic approach is recommended to determine which interspecific combinations are most likely to be useful in a breeding program. As we deepen our understanding of how genetic diversity is partitioned within and between cultivated and wild gene pools of Oryza, breeders will have increased power to make predictions about the most efficient strategies for utilizing wild germplasm for rice improvement.

Selection in the making: a worldwide survey of haplotypic diversity around a causative mutation in porcine IGF2

Domestic species allow us to study dramatic evolutionary changes at an accelerated rate due to the effectiveness of modern breeding techniques and the availability of breeds that have undergone distinct selection pressures. We present a worldwide survey of haplotype variability around a known causative mutation in porcine gene IGF2, which increases lean content. We genotyped 34 SNPs spanning 27 kb in 237 domestic pigs and 162 wild boars. Although the selective process had wiped out variability for at least 27 kb in the haplotypes carrying the mutation, there was no indication of an overall reduction in genetic variability of international vs. European local breeds; there was also no evidence of a reduction in variability caused by domestication. The haplotype structure and a plot of Tajima's D against the frequency of the causative mutation across breeds suggested a temporal pattern, where each breed corresponded to a different selective stage. This was observed comparing the haplotype neighbor-joining (NJ) trees of breeds that have undergone increasing selection pressures for leanness, e.g., European local breeds vs. Pietrain. These results anticipate that comparing current domestic breeds will decisively help to recover the genetic history of domestication and contemporary selective processes.

The extent of linkage disequilibrium in rice (Oryza sativa L.)

Despite its status as one of the world's major crops, linkage disequilibrium (LD) patterns have not been systematically characterized across the genome of Asian rice (Oryza sativa). Such information is critical to fully exploit the genome sequence for mapping complex traits using association techniques. Here we characterize LD in five 500-kb regions of the rice genome in three major cultivated rice varieties (indica, tropical japonica, and temperate japonica) and in the wild ancestor of Asian rice, Oryza rufipogon. Using unlinked SNPs to determine the amount of background linkage disequilibrium in each population, we find that the extent of LD is greatest in temperate japonica (probably >500 kb), followed by tropical japonica (approximately 150 kb) and indica (approximately 75 kb). LD extends over a shorter distance in O. rufipogon (<<40 kb) than in any of the O. sativa groups assayed here. The differences in the extent of LD among these groups are consistent with differences in outcrossing and recombination rate estimates. As well as heterogeneity between groups, our results suggest variation in LD patterns among genomic regions. We demonstrate the feasibility of genomewide association mapping in cultivated Asian rice using a modest number of SNPs.