The pattern of polymorphism in Arabidopsis thaliana

Temporal dynamics and linkage disequilibrium in natural Caenorhabditis elegans populations

Caenorhabditis elegans is a major laboratory model system yet a newcomer to the field of population genetics, and relatively little is known of its biology in the wild. Recent studies of natural populations at a single time point revealed strong spatial population structure and suggested that these populations may be very dynamic. We have therefore studied several natural C. elegans populations over time and genotyped them at polymorphic microsatellite loci. While some populations appear to be genetically stable over the course of observation, others seem to go extinct, with full replacement of multilocus genotypes upon regrowth. The frequency of heterozygotes indicates that outcrossing occurs at a mean frequency of 1.7% and is variable between populations. However, in genetically stable populations, linkage disequilibrium between different chromosomes can be maintained over several years at a level much higher than expected from the heterozygote frequency. C. elegans seems to follow metapopulation dynamics, and the maintenance of linkage disequilibrium despite a low yet significant level of outcrossing suggests that selection may act against the progeny of outcrossings.

Natural genetic variation in Arabidopsis: tools, traits and prospects for evolutionary ecology

BACKGROUND

The model plant Arabidopsis thaliana (Arabidopsis) shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, the potential of Arabidopsis for molecular genetic analysis of this natural variation has increased dramatically in recent years.

SCOPE

Advanced genomics has accelerated molecular phylogenetic analysis and gene identification by quantitative trait loci (QTL) mapping and/or association mapping in Arabidopsis. In particular, QTL mapping utilizing natural accessions is now becoming a major strategy of gene isolation, offering an alternative to artificial mutant lines. Furthermore, the genomic information is used by researchers to uncover the signature of natural selection acting on the genes that contribute to phenotypic variation. The evolutionary significance of such genes has been evaluated in traits such as disease resistance and flowering time. However, although molecular hallmarks of selection have been found for the genes in question, a corresponding ecological scenario of adaptive evolution has been difficult to prove. Ecological strategies, including reciprocal transplant experiments and competition experiments, and utilizing near-isogenic lines of alleles of interest will be a powerful tool to measure the relative fitness of phenotypic and/or allelic variants.

CONCLUSIONS

As the plant model organism, Arabidopsis provides a wealth of molecular background information for evolutionary genetics. Because genetic diversity between and within Arabidopsis populations is much higher than anticipated, combining this background information with ecological approaches might well establish Arabidopsis as a model organism for plant evolutionary ecology.

A soybean transcript map: gene distribution, haplotype and single-nucleotide polymorphism analysis

The first genetic transcript map of the soybean genome was created by mapping one SNP in each of 1141 genes in one or more of three recombinant inbred line mapping populations, thus providing a picture of the distribution of genic sequences across the mapped portion of the genome. Single-nucleotide polymorphisms (SNPs) were discovered via the resequencing of sequence-tagged sites (STSs) developed from expressed sequence tag (EST) sequence. From an initial set of 9459 polymerase chain reaction primer sets designed to a diverse set of genes, 4240 STSs were amplified and sequenced in each of six diverse soybean genotypes. In the resulting 2.44 Mbp of aligned sequence, a total of 5551 SNPs were discovered, including 4712 single-base changes and 839 indels for an average nucleotide diversity of Theta= 0.000997. The analysis of the observed genetic distances between adjacent genes vs. the theoretical distribution based upon the assumption of a random distribution of genes across the 20 soybean linkage groups clearly indicated that genes were clustered. Of the 1141 genes, 291 mapped to 72 of the 112 gaps of 5-10 cM in the preexisting simple sequence repeat (SSR)-based map, while 111 genes mapped in 19 of the 26 gaps >10 cM. The addition of 1141 sequence-based genic markers to the soybean genome map will provide an important resource to soybean geneticists for quantitative trait locus discovery and map-based cloning, as well as to soybean breeders who increasingly depend upon marker-assisted selection in cultivar improvement.

Genetic basis for dosage sensitivity in Arabidopsis thaliana

Aneuploidy, the relative excess or deficiency of specific chromosome types, results in gene dosage imbalance. Plants can produce viable and fertile aneuploid individuals, while most animal aneuploids are inviable or developmentally abnormal. The swarms of aneuploid progeny produced by Arabidopsis triploids constitute an excellent model to investigate the mechanisms governing dosage sensitivity and aneuploid syndromes. Indeed, genotype alters the frequency of aneuploid types within these swarms. Recombinant inbred lines that were derived from a triploid hybrid segregated into diploid and tetraploid individuals. In these recombinant inbred lines, a single locus, which we call SENSITIVE TO DOSAGE IMBALANCE (SDI), exhibited segregation distortion in the tetraploid subpopulation only. Recent progress in quantitative genotyping now allows molecular karyotyping and genetic analysis of aneuploid populations. In this study, we investigated the causes of the ploidy-specific distortion at SDI. Allele frequency was distorted in the aneuploid swarms produced by the triploid hybrid. We developed a simple quantitative measure for aneuploidy lethality and using this measure demonstrated that distortion was greatest in the aneuploids facing the strongest viability selection. When triploids were crossed to euploids, the progeny, which lack severe aneuploids, exhibited no distortion at SDI. Genetic characterization of SDI in the aneuploid swarm identified a mechanism governing aneuploid survival, perhaps by buffering the effects of dosage imbalance. As such, SDI could increase the likelihood of retaining genomic rearrangements such as segmental duplications. Additionally, in species where triploids are fertile, aneuploid survival would facilitate gene flow between diploid and tetraploid populations via a triploid bridge and prevent polyploid speciation. Our results demonstrate that positional cloning of loci affecting traits in populations containing ploidy and chromosome number variants is now feasible using quantitative genotyping approaches.

Transcription-related mutations and GC content drive variation in nucleotide substitution rates across the genomes of Arabidopsis thaliana and Arabidopsis lyrata

Background

There has been remarkably little study of nucleotide substitution rate variation among plant nuclear genes, in part because orthology is difficult to establish. Orthology is even more problematic for intergenic regions of plant nuclear genomes, because plant genomes generally harbor a wealth of repetitive DNA. In theory orthologous intergenic data is valuable for studying rate variation because nucleotide substitutions in these regions should be under little selective constraint compared to coding regions. As a result, evolutionary rates in intergenic regions may more accurately reflect genomic features, like recombination and GC content, that contribute to nucleotide substitution.

Results

We generated a set of 66 intergenic sequences in Arabidopsis lyrata, a close relative of Arabidopsis thaliana. The intergenic regions included transposable element (TE) remnants and regions flanking the TEs. We verified orthology of these amplified regions both by comparison of existing A. lyrata – A. thaliana genetic maps and by using molecular features. We compared substitution rates among the 66 intergenic loci, which exhibit ~5-fold rate variation, and compared intergenic rates to a set of 64 orthologous coding sequences. Our chief observations were that the average rate of nucleotide substitution is slower in intergenic regions than in synonymous sites, that rate variation in both intergenic and coding regions correlate with GC content, that GC content alone is not sufficient to explain differences in rates between intergenic and coding regions, and that rates of evolution in intergenic regions correlate negatively with gene density.

Conclusion

Our observations indicated that mutation rates vary among genomics regions as a function of base composition, suggesting that previous observations of "selective constraint" on non-coding regions could more accurately be attributed to a GC effect instead of selection. The negative correlation between nucleotide substitution rate and gene density provides a potential neutral explanation for a previously documented correlation between gene density and polymorphism levels within A. thaliana. Finally, we discuss potential forces that could contribute to rapid synonymous rates, and provide evidence to suggest that transcription-related mutation contributes to rate differences between intergenic and synonymous sites.

The genetic architecture of shoot branching in Arabidopsis thaliana: a comparative assessment of candidate gene associations vs. quantitative trait locus mapping

Association mapping focused on 36 genes involved in branch development was used to identify candidate genes for variation in shoot branching in Arabidopsis thaliana. The associations between four branching traits and moderate-frequency haplogroups at the studied genes were tested in a panel of 96 accessions from a restricted geographic range in Central Europe. Using a mixed-model association-mapping method, we identified three loci--MORE AXILLARY GROWTH 2 (MAX2), MORE AXILLARY GROWTH 3 (MAX3), and SUPERSHOOT 1 (SPS1)--that were significantly associated with branching variation. On the basis of a more extensive examination of the MAX2 and MAX3 genomic regions, we find that linkage disequilibrium in these regions decays within approximately 10 kb and trait associations localize to the candidate genes in these regions. When the significant associations are compared to relevant quantitative trait loci (QTL) from previous Ler x Col and Cvi x Ler recombinant inbred line (RIL) mapping studies, no additive QTL overlapping these candidate genes are observed, although epistatic QTL for branching, including one that spans the SPS1, are found. These results suggest that epistasis is prevalent in determining branching variation in A. thaliana and may need to be considered in linkage disequilibrium mapping studies of genetically diverse accessions.

Highly variable patterns of linkage disequilibrium in multiple soybean populations

Prospects for utilizing whole-genome association analysis in autogamous plant populations appear promising due to the reported high levels of linkage disequilibrium (LD). To determine the optimal strategies for implementing association analysis in soybean (Glycine max L. Merr.), we analyzed the structure of LD in three regions of the genome varying in length from 336 to 574 kb. This analysis was conducted in four distinct groups of soybean germplasm: 26 accessions of the wild ancestor of soybean (Glycine soja Seib. et Zucc.); 52 Asian G. max Landraces, the immediate results of domestication from G. soja; 17 Asian Landrace introductions that became the ancestors of North American (N. Am.) cultivars, and 25 Elite Cultivars from N. Am. In G. soja, LD did not extend past 100 kb; however, in the three cultivated G. max groups, LD extended from 90 to 574 kb, likely due to the impacts of domestication and increased self-fertilization. The three genomic regions were highly variable relative to the extent of LD within the three cultivated soybean populations. G. soja appears to be ideal for fine mapping of genes, but due to the highly variable levels of LD in the Landraces and the Elite Cultivars, whole-genome association analysis in soybean may be more difficult than first anticipated.

Global eQTL mapping reveals the complex genetic architecture of transcript-level variation in Arabidopsis

The genetic architecture of transcript-level variation is largely unknown. The genetic determinants of transcript-level variation were characterized in a recombinant inbred line (RIL) population (n = 211) of Arabidopsis thaliana using whole-genome microarray analysis and expression quantitative trait loci (eQTL) mapping of transcript levels as expression traits (e-traits). Genetic control of transcription was highly complex: one-third of the quantitatively controlled transcripts/e-traits were regulated by cis-eQTL, and many trans-eQTL mapped to hotspots that regulated hundreds to thousands of e-traits. Several thousand eQTL of large phenotypic effect were detected, but almost all (93%) of the 36,871 eQTL were associated with small phenotypic effects (R(2) < 0.3). Many transcripts/e-traits were controlled by multiple eQTL with opposite allelic effects and exhibited higher heritability in the RILs than their parents, suggesting nonadditive genetic variation. To our knowledge, this is the first large-scale global eQTL study in a relatively large plant mapping population. It reveals that the genetic control of transcript level is highly variable and multifaceted and that this complexity may be a general characteristic of eukaryotes.

Background selection in single genes may explain patterns of codon bias

Background selection involves the reduction in effective population size caused by the removal of recurrent deleterious mutations from a population. Previous work has examined this process for large genomic regions. Here we focus on the level of a single gene or small group of genes and investigate how the effects of background selection caused by nonsynonymous mutations are influenced by the lengths of coding sequences, the number and length of introns, intergenic distances, neighboring genes, mutation rate, and recombination rate. We generate our predictions from estimates of the distribution of the fitness effects of nonsynonymous mutations, obtained from DNA sequence diversity data in Drosophila. Results for genes in regions with typical frequencies of crossing over in Drosophila melanogaster suggest that background selection may influence the effective population sizes of different regions of the same gene, consistent with observed differences in codon usage bias along genes. It may also help to cause the observed effects of gene length and introns on codon usage. Gene conversion plays a crucial role in determining the sizes of these effects. The model overpredicts the effects of background selection with large groups of nonrecombining genes, because it ignores Hill-Robertson interference among the mutations involved.

Maize Sh2 gene is constrained by natural selection but escaped domestication

In Zea mays L., we studied the molecular evolution of Shrunken2 (Sh2), a gene that encodes the large subunits of a major enzyme in endosperm starch biosynthesis, ADP-glucose pyrophosphorylase. We compared 4669 bp of the Sh2 coding region on 50 accessions of maize and teosinte. Very few nucleotide polymorphisms were found when compared with other genes in Z. mays, revealing an effect of purifying selection in the whole species that predates domestication. Additionally, the comparison of Sh2 sequences in all Z. mays subspecies and outgroups Z. diploperennis and Tripsacum dactyloides suggests the occurrence of an ancient selective sweep in the Sh2 3' region. The amount and nature of nucleotide diversity are similar in both maize and teosinte, confirming previous results that suggested that Sh2 has not been involved in maize domestication. The very low level of nucleotide diversity as well as the highly conserved protein sequence suggest that natural selection retained effective Sh2 allele(s) long before agriculture started, making human selection inefficient on this gene.

A molecular isolation mechanism associated with high intra-specific diversity in rice

High levels of inter-specific diversity are expected due to genetic isolation, the reproductive or geographical barriers, which lead to the accumulation of nucleotide variation. However, high levels of genetic variation are repeatedly observed even within species, notably at loci of the human major histocompatability complex and of plant resistance genes. Are molecular isolations responsible for the high intra-specific variation? To address this issue, we performed a genome-wide survey of the relationship between the possible factors that could cause genetic isolation, and the level of polymorphism, based on two rice genome comparisons. Here, we show that the levels of polymorphism in rice genes are positively correlated with the proportions of non-alignable flanking sequences, and that the correlation is observed even in single-copy genes. The physical locations of the genes were also investigated, and a strong association between the asymmetric architecture of genomes and the levels of polymorphism was revealed. These results suggest that the flank heterogeneity and the asymmetric architecture between genomes serve as isolation mechanisms at the molecular level that result in accumulation of higher genetic variation. This mechanism is of fundamental importance to understand natural genetic variation within species.

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits

A central challenge in evolutionary biology is to identify genes underlying ecologically important traits and describe the fitness consequences of naturally occurring variation at these loci. To address this goal, several novel approaches have been developed, including 'population genomics,' where a large number of molecular markers are scored in individuals from different environments with the goal of identifying markers showing unusual patterns of variation, potentially due to selection at linked sites. Such approaches are appealing because of (1) the increasing ease of generating large numbers of genetic markers, (2) the ability to scan the genome without measuring phenotypes and (3) the simplicity of sampling individuals without knowledge of their breeding history. Although such approaches are inherently applicable to non-model systems, to date these studies have been limited in their ability to uncover functionally relevant genes. By contrast, quantitative genetics has a rich history, and more recently, quantitative trait locus (QTL) mapping has had some success in identifying genes underlying ecologically relevant variation even in novel systems. QTL mapping, however, requires (1) genetic markers that specifically differentiate parental forms, (2) a focus on a particular measurable phenotype and (3) controlled breeding and maintenance of large numbers of progeny. Here we present current advances and suggest future directions that take advantage of population genomics and quantitative genetic approaches - in both model and non-model systems. Specifically, we discuss advantages and limitations of each method and argue that a combination of the two provides a powerful approach to uncovering the molecular mechanisms responsible for adaptation.

Genetic evidence for a second domestication of barley (Hordeum vulgare) east of the Fertile Crescent

Cereal agriculture originated with the domestication of barley and early forms of wheat in the Fertile Crescent. There has long been speculation that barley was domesticated more than once. We use differences in haplotype frequency among geographic regions at multiple loci to infer at least two domestications of barley; one within the Fertile Crescent and a second 1,500-3,000 km farther east. The Fertile Crescent domestication contributed the majority of diversity in European and American cultivars, whereas the second domestication contributed most of the diversity in barley from Central Asia to the Far East.

Single-nucleotide polymorphism frequency in a set of selected lines of bread wheat (Triticum aestivum L.)

Information on single-nucleotide polymorphisms (SNPs) in hexaploid bread wheat is still scarce. The goal of this study was to detect SNPs in wheat and examine their frequency. Twenty-six bread wheat lines from different origins worldwide were used. Specific PCR-products were obtained from 21 genes and directly sequenced. SNPs were discovered from the alignment of these sequences. The overall sequence polymorphism observed in this sample appears to be low; 64 single-base polymorphisms were detected in approximately 21.5 kb (i.e., 1 SNP every 335 bp). The level of polymorphism is highly variable among the different genes studied. Fifty percent of the genes studied contained no sequence polymorphism, whereas most SNPs detected were located in only 2 genes. As expected, taking into account a synthetic line created with a wild Triticum tauschii parent increases the level of polymorphism (101 SNPs; 1 SNP every 212 bp). The detected SNPs are available at http://urgi.versailles.inra.fr/GnpSNP">http://urgi.versailles.inra.fr/GnpSNP. Data on linkage disequilibrium (LD) are still preliminary. They showed a significant level of LD in the 2 most polymorphic genes. To conclude, the genome size of hexaploid wheat and its low level of polymorphism complicate SNP discovery in this species.

Unique pattern of R-gene variation within populations in Arabidopsis

An understanding of the variation pattern in disease resistance (R) genes is essential for its use in breeding programs aimed at neutralizing the threat of pathogens. Although the variation between populations is well known, there is little research about R-gene variation patterns within populations. Here, we investigate the polymorphism at three R-gene loci of 39 individual plants from nine populations of Arabidopsis thaliana. Our data suggest that alleles of each locus from individuals within a local population were either nearly identical, or highly diverse as ones between populations. The vast majority (92.5%) of within-population variation was shared globally, with high levels of allelic diversity (up to 11.7%) and abundant diverse-alleles. This unique pattern of within-population variation at R-loci suggests that individual plants within a population had the great potential to maintain a high level of globally-shared polymorphisms, and that the diversifying selection was the major force maintaining such polymorphisms. Consequently, the shared-polymorphism became recyclable for new R-genes, as the corresponding avirulence re-emerges in pathogen populations.

Effects of introgression and recombination on haplotype structure and linkage disequilibrium surrounding a locus encoding Bymovirus resistance in barley

We present a detailed analysis of linkage disequilibrium (LD) in the physical and genetic context of the barley gene Hv-eIF4E, which confers resistance to the barley yellow mosaic virus (BYMV) complex. Eighty-three SNPs distributed over 132 kb of Hv-eIF4E and six additional fragments genetically mapped to its flanking region were used to derive haplotypes from 131 accessions. Three haplogroups were recognized, discriminating between the alleles rym4 and rym5, which each encode for a spectrum of resistance to BYMV. With increasing map distance, haplotypes of susceptible genotypes displayed diverse patterns driven mainly by recombination, whereas haplotype diversity within the subgroups of resistant genotypes was limited. We conclude that the breakdown of LD within 1 cM of the resistance gene was generated mainly by susceptible genotypes. Despite the LD decay, a significant association between haplotype and resistance to BYMV was detected up to a distance of 5.5 cM from the resistance gene. The LD pattern and the haplotype structure of the target chromosomal region are the result of interplay between low recombination and recent breeding history.

Large-scale DNA polymorphism study of Oryza sativa and O. rufipogon reveals the origin and divergence of Asian rice

Polymorphism over approximately 26 kb of DNA sequence spanning 22 loci and one region distributed on chromosomes 1, 2, 3 and 4 was studied in 30 accessions of cultivated rice, Oryza sativa, and its wild relatives. Phylogenetic analysis using all the DNA sequences suggested that O. sativa ssp. indica and ssp. japonica were independently domesticated from a wild species O. rufipogon. O. sativa ssp. indica contained substantial genetic diversity (pi = 0.0024), whereas ssp. japonica exhibited extremely low nucleotide diversity (pi = 0.0001) suggesting the origin of the latter from a small number of founders. O. sativa ssp. japonica contained a larger number of derived and fixed non-synonymous substitutions as compared to ssp. indica. Nucleotide diversity and genealogical history substantially varied across the 22 loci. A locus, RLD15 on chromosome 2, showed a distinct genealogy with ssp. japonica sequences distantly separated from those of O. rufipogon and O. sativa ssp. indica. Linkage disequilibrium (LD) was analyzed in two different regions. LD in O. rufipogon decays within 5 kb, whereas it extends to approximately 50 kb in O. sativa ssp. indica.

Purdue ionomics information management system. An integrated functional genomics platform

The advent of high-throughput phenotyping technologies has created a deluge of information that is difficult to deal with without the appropriate data management tools. These data management tools should integrate defined workflow controls for genomic-scale data acquisition and validation, data storage and retrieval, and data analysis, indexed around the genomic information of the organism of interest. To maximize the impact of these large datasets, it is critical that they are rapidly disseminated to the broader research community, allowing open access for data mining and discovery. We describe here a system that incorporates such functionalities developed around the Purdue University high-throughput ionomics phenotyping platform. The Purdue Ionomics Information Management System (PiiMS) provides integrated workflow control, data storage, and analysis to facilitate high-throughput data acquisition, along with integrated tools for data search, retrieval, and visualization for hypothesis development. PiiMS is deployed as a World Wide Web-enabled system, allowing for integration of distributed workflow processes and open access to raw data for analysis by numerous laboratories. PiiMS currently contains data on shoot concentrations of P, Ca, K, Mg, Cu, Fe, Zn, Mn, Co, Ni, B, Se, Mo, Na, As, and Cd in over 60,000 shoot tissue samples of Arabidopsis (Arabidopsis thaliana), including ethyl methanesulfonate, fast-neutron and defined T-DNA mutants, and natural accession and populations of recombinant inbred lines from over 800 separate experiments, representing over 1,000,000 fully quantitative elemental concentrations. PiiMS is accessible at www.purdue.edu/dp/ionomics.

Variation and constraint in plant evolution and development

The goal of this short review is to consider the interrelated phenomena of phenotypic variation and genetic constraint with respect to plant diversity. The unique aspects of plants, including sessile habit, modular growth and diverse developmental programs expressed at the phytomer level, merit a specific examination of the genetic basis of their phenotypic variation, and how they experience and escape genetic constraint. Numerous QTL studies with wild and domesticated plants reveal that most phenotypic traits are polygenic but vary in the number and effect of the loci contributing, from a few loci of large effects to many with small effects. Further, somatic mutations, developmental plasticity and epigenetic variation, especially gene methylation, can contribute to increases in phenotypic variation. The flip side of these processes, genetic constraint, can similarly be the result of many factors, including pleiotropy, canalization and genetic redundancy. Genetic constraint is not only a mechanism to prevent change, however, it can also serve to direct evolution along certain paths. Ultimately, genetic constraint often comes full circle and is released through events such as hybridization, genome duplication and epigenetic remodeling. We are just beginning to understand how these processes can operate simultaneously during the evolution of ecologically important traits in plants.

VIM1, a methylcytosine-binding protein required for centromeric heterochromatinization

Epigenetic regulation in eukaryotes is executed by a complex set of signaling interactions among small RNA species and chromatin marks, including histone modification and DNA methylation. We identified vim1 (VARIANT IN METHYLATION 1), an Arabidopsis mutation causing cytosine hypomethylation and decondensation of centromeres in interphase. VIM1 is a member of a small gene family, encoding proteins containing PHD, RING, and SRA (SET- and RING-associated) domains, which are found together in mammalian proteins implicated in regulation of chromatin modification, transcription, and the cell cycle. VIM1 is an unconventional methylcytosine-binding protein that interacts in vitro with 5mCpG- and 5mCpHpG-modified DNA (via its SRA domain), as well as recombinant histones (H2B, H3, H4, and HTR12) in plant extracts. VIM1 associates with methylated genomic loci in vivo and is enriched in chromocenters. Our findings suggest that VIM1 acts at the DNA methylation-histone interface to maintain centromeric heterochromatin.

Regulatory network construction in Arabidopsis by using genome-wide gene expression quantitative trait loci

Accessions of a plant species can show considerable genetic differences that are analyzed effectively by using recombinant inbred line (RIL) populations. Here we describe the results of genome-wide expression variation analysis in an RIL population of Arabidopsis thaliana. For many genes, variation in expression could be explained by expression quantitative trait loci (eQTLs). The nature and consequences of this variation are discussed based on additional genetic parameters, such as heritability and transgression and by examining the genomic position of eQTLs versus gene position, polymorphism frequency, and gene ontology. Furthermore, we developed an approach for genetic regulatory network construction by combining eQTL mapping and regulator candidate gene selection. The power of our method was shown in a case study of genes associated with flowering time, a well studied regulatory network in Arabidopsis. Results that revealed clusters of coregulated genes and their most likely regulators were in agreement with published data, and unknown relationships could be predicted.

Levels of linkage disequilibrium in a wild bird population

Population-based mapping approaches are attractive for tracing the genetic background to phenotypic traits in wild species, given that it is often difficult to gather extensive and well-defined pedigrees needed for quantitative trait locus analysis. However, the feasibility of association or hitch-hiking mapping is dependent on the degree of linkage disequilibrium (LD) in the population, on which there is yet limited information for wild species. Here we use single nucleotide polymorphism (SNP) markers from 23 genes in a recently established linkage map of the Z chromosome of the collared flycatcher, to study the extent of LD in a natural bird population. In most but not all cases we find SNPs within the same intron (less than 500 bp) to be in perfect LD. However, LD then decays to background level at a distance 1cM or 400-500 kb. Although LD seems more extensive than in other species, if the observed pattern is representative for other regions of the genome and turns out to be a general feature of natural bird populations, dense marker maps might be needed for genome scans aimed at identifying association between marker and trait loci.

The spatial scale of genetic differentiation in a model organism: the wild yeast Saccharomyces paradoxus

Little information is presently available on the factors promoting genetic divergence in eukaryotic microbes. We studied the spatial distribution of genetic variation in Saccharomyces paradoxus, the wild relative of Saccharomyces cerevisiae, from the scale of a few centimetres on individual oak trees to thousands of kilometers across different continents. Genealogical analysis of six loci shows that isolates from Europe form a single recombining population, and within this population genetic differentiation increases with physical distance. Between different continents, strains are more divergent and genealogically independent, indicating well-differentiated lineages that may be in the process of speciation. Such replicated populations will be useful for studies in population genomics.

Fitness consequences of infection of Arabidopsis thaliana with its natural bacterial pathogen Pseudomonas viridiflava

Variation in plant resistance to pathogen infection is commonly observed in interactions between wild plants and their foliar pathogens. Models of host-pathogen interactions indicate that a large cost of infection is generally necessary to maintain this variation, yet there is limited evidence that foliar pathogens cause detectable fitness reductions in wild host plants. Most published work has focused on fungal pathogens. Pseudomonas viridiflava, a common bacterial pathogen of the annual weed Arabidopsis thaliana across its range, comprises two distinct genetic clades that cause disease symptoms of different severity. Here we measured the extent of infection of wild A. thaliana populations in the Midwest, USA, and examined the effect on seed production, in field and growth-chamber experiments, of experimental inoculation with isolates from the two clades. We found infection with P. viridiflava varied from 0 to 56% in Midwest A. thaliana populations, with the possibility of several leaves per plant infected later in the growing season. In the growth chambers, experimental inoculation reduced seed set by averages of 15 and 11% for clades A and B, respectively. In the field experiment, only clade A affected plant fitness significantly, reducing seed set by an average of 38%. Underlying these average effects we observed both negative and positive effects of infection, and variation in both fitness among plant genotypes and sensitivity to environmental conditions.

Microsatellite diversity and broad scale geographic structure in a model legume: building a set of nested core collection for studying naturally occurring variation in Medicago truncatula

BACKGROUND

Exploiting genetic diversity requires previous knowledge of the extent and structure of the variation occurring in a species. Such knowledge can in turn be used to build a core-collection, i.e. a subset of accessions that aim at representing the genetic diversity of this species with a minimum of repetitiveness. We investigate the patterns of genetic diversity and population structure in a collection of 346 inbred lines representing the breadth of naturally occurring diversity in the Legume plant model Medicago truncatula using 13 microsatellite loci distributed throughout the genome.

RESULTS

We confirm the uniqueness of all these genotypes and reveal a large amount of genetic diversity and allelic variation within this autogamous species. Spatial genetic correlation was found only for individuals originating from the same population and between neighbouring populations. Using a model-based clustering algorithm, we identified four main genetic clusters in the set of individuals analyzed. This stratification matches broad geographic regions. We also identified a set of "admixed" individuals that do not fit with this population structure scheme.

CONCLUSION

The stratification inferred is discussed considering potential historical events like expansion, refuge history and admixture between neighbouring groups. Information on the allelic richness and the inferred population structure are used to build a nested core-collection. The set of inbred lines and the core collections are publicly available and will help coordinating efforts for the study of naturally occurring variation in the growing Medicago truncatula community.

Selection, recombination and demographic history in Drosophila miranda

Selection, recombination, and the demographic history of a species can all have profound effects on genomewide patterns of variability. To assess the impact of these forces in the genome of Drosophila miranda, we examine polymorphism and divergence patterns at 62 loci scattered across the genome. In accordance with recent findings in D. melanogaster, we find that noncoding DNA generally evolves more slowly than synonymous sites, that the distribution of polymorphism frequencies in noncoding DNA is significantly skewed toward rare variants relative to synonymous sites, and that long introns evolve significantly slower than short introns or synonymous sites. These observations suggest that most noncoding DNA is functionally constrained and evolving under purifying selection. However, in contrast to findings in the D. melanogaster species group, we find little evidence of adaptive evolution acting on either coding or noncoding sequences in D. miranda. Levels of linkage disequilibrium (LD) in D. miranda are comparable to those observed in D. melanogaster, but vary considerably among chromosomes. These patterns suggest a significantly lower rate of recombination on autosomes, possibly due to the presence of polymorphic autosomal inversions and/or differences in chromosome sizes. All chromosomes show significant departures from the standard neutral model, including too much heterogeneity in synonymous site polymorphism relative to divergence among loci and a general excess of rare synonymous polymorphisms. These departures from neutral equilibrium expectations are discussed in the context of nonequilibrium models of demography and selection.

An Arabidopsis example of association mapping in structured samples

A potentially serious disadvantage of association mapping is the fact that marker-trait associations may arise from confounding population structure as well as from linkage to causative polymorphisms. Using genome-wide marker data, we have previously demonstrated that the problem can be severe in a global sample of 95 Arabidopsis thaliana accessions, and that established methods for controlling for population structure are generally insufficient. Here, we use the same sample together with a number of flowering-related phenotypes and data-perturbation simulations to evaluate a wider range of methods for controlling for population structure. We find that, in terms of reducing the false-positive rate while maintaining statistical power, a recently introduced mixed-model approach that takes genome-wide differences in relatedness into account via estimated pairwise kinship coefficients generally performs best. By combining the association results with results from linkage mapping in F2 crosses, we identify one previously known true positive and several promising new associations, but also demonstrate the existence of both false positives and false negatives. Our results illustrate the potential of genome-wide association scans as a tool for dissecting the genetics of natural variation, while at the same time highlighting the pitfalls. The importance of study design is clear; our study is severely under-powered both in terms of sample size and marker density. Our results also provide a striking demonstration of confounding by population structure. While statistical methods can be used to ameliorate this problem, they cannot always be effective and are certainly not a substitute for independent evidence, such as that obtained via crosses or transgenic experiments. Ultimately, association mapping is a powerful tool for identifying a list of candidates that is short enough to permit further genetic study.

There is currently tremendous interest in using association mapping to find the genes responsible for natural variation, particularly for human disease. In association mapping, researchers seek to identify regions of the genome where individuals who are phenotypically similar (e.g., they all have the same disease) are also unusually closely related. A potentially serious problem is that spurious correlations may arise if the population is structured so that members of a subgroup tend to be much more closely related. We have previously demonstrated that this problem can be severe in Arabidopsis thaliana, and that established statistical methods for controlling for population structure are insufficient. Here, we evaluate a broader range of methods. We find that a recently introduced mixed-model approach generally performs best. By combining the association results with results from linkage mapping in F2 crosses, we identify one previously known true positive and several promising new associations, but also demonstrate the existence of both false positives and false negatives. Our results illustrate the potential of genome-wide association scans as a tool for dissecting the genetics of natural variation, while at the same time highlighting the pitfalls.

Modern computational approaches for analysing molecular genetic variation data

An explosive growth is occurring in the quantity, quality and complexity of molecular variation data that are being collected. Historically, such data have been analysed by using model-based methods. Models are useful for sharpening intuition, for explanation and for prediction: they add to our understanding of how the data were formed, and they can provide quantitative answers to questions of interest. We outline some of these model-based approaches, including the coalescent, and discuss the applicability of the computational methods that are necessary given the highly complex nature of current and future data sets.

Natural variants of AtHKT1 enhance Na+ accumulation in two wild populations of Arabidopsis

Plants are sessile and therefore have developed mechanisms to adapt to their environment, including the soil mineral nutrient composition. Ionomics is a developing functional genomic strategy designed to rapidly identify the genes and gene networks involved in regulating how plants acquire and accumulate these mineral nutrients from the soil. Here, we report on the coupling of high-throughput elemental profiling of shoot tissue from various Arabidopsis accessions with DNA microarray-based bulk segregant analysis and reverse genetics, for the rapid identification of genes from wild populations of Arabidopsis that are involved in regulating how plants acquire and accumulate Na⁺ from the soil. Elemental profiling of shoot tissue from 12 different Arabidopsis accessions revealed that two coastal populations of Arabidopsis collected from Tossa del Mar, Spain, and Tsu, Japan (Ts-1 and Tsu-1, respectively), accumulate higher shoot levels of Na⁺ than do Col-0 and other accessions. We identify AtHKT1, known to encode a Na⁺ transporter, as being the causal locus driving elevated shoot Na⁺ in both Ts-1 and Tsu-1. Furthermore, we establish that a deletion in a tandem repeat sequence approximately 5 kb upstream of AtHKT1 is responsible for the reduced root expression of AtHKT1 observed in these accessions. Reciprocal grafting experiments establish that this loss of AtHKT1 expression in roots is responsible for elevated shoot Na⁺. Interestingly, and in contrast to the hkt1–1 null mutant, under NaCl stress conditions, this novel AtHKT1 allele not only does not confer NaCl sensitivity but also cosegregates with elevated NaCl tolerance. We also present all our elemental profiling data in a new open access ionomics database, the Purdue Ionomics Information Management System (PiiMS; http://www.purdue.edu/dp/ionomics). Using DNA microarray-based genotyping has allowed us to rapidly identify AtHKT1 as the casual locus driving the natural variation in shoot Na⁺ accumulation we observed in Ts-1 and Tsu-1. Such an approach overcomes the limitations imposed by a lack of established genetic markers in most Arabidopsis accessions and opens up a vast and tractable source of natural variation for the identification of gene function not only in ionomics but also in many other biological processes.

Unlike most animals, plants are sessile and cannot leave a poor-quality environment after germinating. They therefore need to tolerate the particular conditions they encounter to survive. This makes plants an ideal system for the study of adaptive variation, and this is particularly true of Arabidopsis thaliana (Arabidopsis), which shows substantial natural variation and for which numerous genetic tools exist. Using a combination of analytical chemistry, genetics, and genomics, the authors were able to identify the specific genetic alteration that drive the natural variation in shoot sodium (Na⁺) accumulation capacity observed in Arabidopsis populations from coastal regions of Spain and Japan (Tossa del Mar and Tsu, respectively). They observed that a deletion in the DNA responsible for regulating the expression of HKT1, a gene known to encode for a Na⁺ transporter, causes reduced expression of AtHKT1 in roots of both the Spanish and Japanese populations. Such altered expression results in the elevated shoot Na⁺ observed in these two populations. Interestingly, this novel version of the HKT1 genes is also associated genetically with the enhanced NaCl resistance they observe in the Japanese population.

Ecotype-specific and chromosome-specific expansion of variant centromeric satellites in Arabidopsis thaliana

Despite the conserved roles and conserved protein machineries of centromeres, their nucleotide sequences can be highly diverse even among related species. The diversity reflects rapid evolution, but the underlying mechanism is largely unknown. One approach to monitor rapid evolution is examination of intra-specific variation. Here we report variant centromeric satellites of Arabidopsis thaliana found through survey of 103 natural accessions (ecotypes). Among them, a cluster of variant centromeric satellites was detected in one ecotype, Cape Verde Islands (Cvi). Recombinant inbred mapping revealed that the variant satellites are distributed in centromeric region of the chromosome 5 (CEN5) of this ecotype. This apparently recent variant accumulation is associated with large deletion of a pericentromeric region and the expansion of satellite region. The variant satellite was bound to HTR12 (centromeric variant histone H3), although expansion of the satellite was not associated with comparable increase in the HTR12 binding. The results suggest that variant satellites with centromere function can rapidly accumulate in one centromere, supporting the model that the satellite repeats in the array are homogenized by occasional unequal crossing-over, which has a potential to generate an expansion of local sequence variants within a centromere cluster.

Population genetic structure of Arabidopsis lyrata in Europe

Population genetic theory predicts that the self-incompatible and perennial herb, Arabidopsis lyrata, will have a genetic structure that differs from the self-fertilizing, annual Arabidopsis thaliana. We quantified the genetic structure for eight populations of A. lyrata ssp. petraea in historically nonglaciated regions of central Europe. Analysis of 20 microsatellite loci for 344 individuals demonstrated that, in accordance with predictions, diploid populations had high genome-wide heterozygosity (H(O) = 0.48; H(E) = 0.52), high within-population diversity (83% of total) compatible with mutation-drift equilibrium, and moderate differentiation among populations (F(ST) = 0.17). Within a single population, the vast majority of genetic variability (92%) was found at the smallest spatial scale (< 3 m). Although there was no evidence of biparental inbreeding or clonal propagation at this scale (F(IS) = 0.003), significant fine-scale spatial autocorrelation indicated localized gene flow presumably due to gravity dispersed seeds (Sp = 0.018). Limited gene flow between isolated population clusters (regions) separated by hundreds of kilometres has given rise to an isolation by distance pattern of diversification, with low, but significant, differentiation among regions (F(ST) = 0.05). The maintenance of geographically widespread polymorphisms and uniformly high diversity throughout central Europe is consistent with periglacial survival of A. lyrata ssp. petraea north of the Alps in steppe-tundra habitats during the last glacial maximum. As expected of northern and previously glaciated localities, A. lyrata in Iceland was genetically less diverse and highly differentiated from central Europe (H(E) = 0.37; F(ST) = 0.27).

Meiotic recombination hotspots in plants

Many studies have demonstrated that the distribution of meiotic crossover events along chromosomes is non-random in plants and other species with sexual reproduction. Large differences in recombination frequencies appear at several scales. On a large scale, regions of high and low rates of crossover have been found to alternate along the chromosomes in all plant species studied. High crossover rates have been reported to be correlated with several chromosome features (e.g. gene density and distance to the centromeres). However, most of these correlations cannot be extended to all plant species. Only a few plant species have been studied on a finer scale. Hotspots of meiotic recombination (i.e. DNA fragments of a few kilobases in length with a higher rate of recombination than the surrounding DNA) have been identified in maize and rice. Most of these hotspots are intragenic. In Arabidopsis thaliana, we have identified several DNA fragments (less than 5 kb in size) with genetic recombination rates at least 5 times higher than the whole-chromosome average [4.6 cM (centimorgan)/Mb], which are therefore probable hotspots for meiotic recombination. Most crossover breakpoints lie in intergenic or non-coding regions. Major efforts should be devoted to characterizing meiotic recombination at the molecular level, which should help to clarify the role of this process in genome evolution.

Genome-wide DNA polymorphism analyses using VariScan

Background

DNA sequence polymorphisms analysis can provide valuable information on the evolutionary forces shaping nucleotide variation, and provides an insight into the functional significance of genomic regions. The recent ongoing genome projects will radically improve our capabilities to detect specific genomic regions shaped by natural selection. Current available methods and software, however, are unsatisfactory for such genome-wide analysis.

Results

We have developed methods for the analysis of DNA sequence polymorphisms at the genome-wide scale. These methods, which have been tested on a coalescent-simulated and actual data files from mouse and human, have been implemented in the VariScan software package version 2.0. Additionally, we have also incorporated a graphical-user interface. The main features of this software are: i) exhaustive population-genetic analyses including those based on the coalescent theory; ii) analysis adapted to the shallow data generated by the high-throughput genome projects; iii) use of genome annotations to conduct a comprehensive analyses separately for different functional regions; iv) identification of relevant genomic regions by the sliding-window and wavelet-multiresolution approaches; v) visualization of the results integrated with current genome annotations in commonly available genome browsers.

Conclusion

VariScan is a powerful and flexible suite of software for the analysis of DNA polymorphisms. The current version implements new algorithms, methods, and capabilities, providing an important tool for an exhaustive exploratory analysis of genome-wide DNA polymorphism data.

Resolving the genetic basis of invasiveness and predicting invasions

Considerable effort has been invested in determining traits underlying invasiveness. Yet, identifying a set of traits that commonly confers invasiveness in a range of species has proven elusive, and almost nothing is known about genetic loci affecting invasive success. Incorporating genetic model organisms into ecologically relevant studies is one promising avenue to begin dissecting the genetic underpinnings of invasiveness. Molecular biologists are rapidly characterizing genes mediating developmental responses to diverse environmental cues, i.e., genes for plasticity, as well as to environmental factors likely to impose strong selection on invading species, e.g., resistance to herbivores and competitors, coordination of life-history events with seasonal changes, and physiological tolerance of heat, drought, or cold. Here, we give an overview of molecular genetic tools increasingly used to characterize the genetic basis of adaptation and that may be used to begin identifying genetic mechanisms of invasiveness. Given the divergent traits that affect invasiveness, "invasiveness genes" common to many clades are unlikely, but the combination of developmental genetic advances with further evolutionary studies and modeling may provide a framework for identifying genes that account for invasiveness in related species.

Natural genetic variation of freezing tolerance in Arabidopsis

Low temperature is a primary determinant of plant growth and survival. Using accessions of Arabidopsis (Arabidopsis thaliana) originating from Scandinavia to the Cape Verde Islands, we show that freezing tolerance of natural accessions correlates with habitat winter temperatures, identifying low temperature as an important selective pressure for Arabidopsis. Combined metabolite and transcript profiling show that during cold exposure, global changes of transcripts, but not of metabolites, correlate with the ability of Arabidopsis to cold acclimate. There are, however, metabolites and transcripts, including several transcription factors, that correlate with freezing tolerance, indicating regulatory pathways that may be of primary importance for this trait. These data identify that enhanced freezing tolerance is associated with the down-regulation of photosynthesis and hormonal responses and the induction of flavonoid metabolism, provide evidence for naturally increased nonacclimated freezing tolerance due to the constitutive activation of the C-repeat binding factors pathway, and identify candidate transcriptional regulators that correlate with freezing tolerance.

Testing for effects of recombination rate on nucleotide diversity in natural populations of Arabidopsis lyrata

We investigated DNA sequence diversity for loci on chromosomes 1 and 2 in six natural populations of Arabidopsis lyrata and tested for the role of natural selection in structuring genomewide patterns of variability, specifically examining the effects of recombination rate on levels of silent polymorphism. In contrast with theoretical predictions from models of genetic hitchhiking, maximum-likelihood-based analyses of diversity and divergence do not suggest reduction of diversity in the region of suppressed recombination near the centromere of chromosome 1, except in a single population from Russia, in which the pericentromeric region may have undergone a local selective sweep or demographic process that reduced variability. We discuss various possibilities that might explain why nucleotide diversity in most A. lyrata populations is not related to recombination rate, including genic recombination hotspots, and low gene density in the low recombination rate region.

Genetic diversity contribution to errors in short oligonucleotide microarray analysis

DNA arrays based on short oligonucleotide (< or = 25-mer) probes are being developed for many species, and are being applied to quantify transcript abundance variation in species with high genetic diversity. To define the parameters necessary to design short oligo arrays for maize (Zea mays L.), a species with particularly high nucleotide (single nucleotide polymorphism, SNP) and insertion-deletion (indel) polymorphism frequencies, we analysed gene expression estimates generated for four maize inbred lines using a custom Affymetrix DNA array, and identified biases associated with high levels of polymorphism between lines. Statistically significant interactions between probes and maize inbreds were detected, affecting five or more probes (out of 30 probes per transcript) in the majority of cases. SNPs and indels were identified by re-sequencing; they are the primary source of probe-by-line interactions, affecting probeset level estimates and reducing the power of detecting transcript level variation between maize inbreds. This analysis identified 36,196 probes in 5118 probesets containing markers that may be used for genotyping in natural and segregating populations for association gene analysis and genetic mapping.

High nucleotide polymorphism and rapid decay of linkage disequilibrium in wild populations of Caenorhabditis remanei

The common ancestor of the self-fertilizing nematodes Caenorhabditis elegans and C. briggsae must have reproduced by obligate outcrossing, like most species in this genus. However, we have only a limited understanding about how genetic variation is patterned in such male-female (gonochoristic) Caenorhabditis species. Here, we report results from surveying nucleotide variation of six nuclear loci in a broad geographic sample of wild isolates of the gonochoristic C. remanei. We find high levels of diversity in this species, with silent-site diversity averaging 4.7%, implying an effective population size close to 1 million. Additionally, the pattern of polymorphisms reveals little evidence for population structure or deviation from neutral expectations, suggesting that the sampled C. remanei populations approximate panmixis and demographic equilibrium. Combined with the observation that linkage disequilibrium between pairs of polymorphic sites decays rapidly with distance, this suggests that C. remanei will provide an excellent system for identifying the genetic targets of natural selection from deviant patterns of polymorphism and linkage disequilibrium. The patterns revealed in this obligately outcrossing species may provide a useful model of the evolutionary circumstances in C. elegans' gonochoristic progenitor. This will be especially important if self-fertilization evolved recently in C. elegans history, because most of the evolutionary time separating C. elegans from its known relatives would have occurred in a state of obligate outcrossing.

The molecular population genetics of shoot development in Arabidopsis thaliana

Studies in Arabidopsis thaliana have provided us with a wealth of information about the genetic pathways that regulate plant morphogenesis. This developmental genetic treasure trove represents a fantastic resource for researchers interested in the microevolution of development. Several laboratories have begun using molecular population genetic analyses to investigate the evolutionary forces that act upon loci that regulate shoot morphogenesis. Much of this work has focused on coding sequence variation in transcription factors; however, recent studies have explored sequence variation in other types of proteins and in promoter regions. Several genes that regulate shoot development contain signatures of selective sweeps associated with positive selection or harbor putative balanced polymorphisms in coding and noncoding sequences. Other regulatory genes appear to be evolving neutrally, but have accumulated potentially deleterious replacement polymorphisms.

A genome-wide survey of R gene polymorphisms in Arabidopsis

We used polymorphism analysis to study the evolutionary dynamics of 27 disease resistance (R) genes by resequencing the leucine-rich repeat (LRR) region in 96 Arabidopsis thaliana accessions. We compared single nucleotide polymorphisms (SNPs) in these R genes to an empirical distribution of SNP in the same sample based on 876 fragments selected to sample the entire genome. LRR regions are highly polymorphic for protein variants but not for synonymous changes, suggesting that they generate many alleles maintained for short time periods. Recombination is also relatively common and important for generating protein variants. Although none of the genes is nearly as polymorphic as RPP13, a locus previously shown to have strong signatures of balancing selection, seven genes show weaker indications of balancing selection. Five R genes are relatively invariant, indicating young alleles, but all contain segregating protein variants. Polymorphism analysis in neighboring fragments yielded inconclusive evidence for recent selective sweeps at these loci. In addition, few alleles are candidates for rapid increases in frequency expected under directional selection. Haplotype sharing analysis revealed significant underrepresentation of R gene alleles with extended haplotypes compared with 1102 random genomic fragments. Lack of convincing evidence for directional selection or selective sweeps argues against an arms race driving R gene evolution. Instead, the data support transient or frequency-dependent selection maintaining protein variants at a locus for variable time periods.

Balancing selection and its effects on sequences in nearby genome regions

Our understanding of balancing selection is currently becoming greatly clarified by new sequence data being gathered from genes in which polymorphisms are known to be maintained by selection. The data can be interpreted in conjunction with results from population genetics models that include recombination between selected sites and nearby neutral marker variants. This understanding is making possible tests for balancing selection using molecular evolutionary approaches. Such tests do not necessarily require knowledge of the functional types of the different alleles at a locus, but such information, as well as information about the geographic distribution of alleles and markers near the genes, can potentially help towards understanding what form of balancing selection is acting, and how long alleles have been maintained.

A high-resolution map of Arabidopsis recombinant inbred lines by whole-genome exon array hybridization

Recombinant populations were the basis for Mendel's first genetic experiments and continue to be key to the study of genes, heredity, and genetic variation today. Genotyping several hundred thousand loci in a single assay by hybridizing genomic DNA to oligonucleotide arrays provides a powerful technique to improve precision linkage mapping. The genotypes of two accessions of Arabidopsis were compared by using a 400,000 feature exon-specific oligonucleotide array. Around 16,000 single feature polymorphisms (SFPs) were detected in ~8,000 of the ~26,000 genes represented on the array. Allelic variation at these loci was measured in a recombinant inbred line population, which defined the location of 815 recombination breakpoints. The genetic linkage map had a total length of 422.5 cM, with 676 informative SFP markers representing intervals of ~0.6 cM. One hundred fifteen single gene intervals were identified. Recombination rate, SFP distribution, and segregation in this population are not uniform. Many genomic regions show a clustering of recombination events including significant hot spots. The precise haplotype structure of the recombinant population was defined with unprecedented accuracy and resolution. The resulting linkage map allows further refinement of the hundreds of quantitative trait loci identified in this well-studied population. Highly variable recombination rates along each chromosome and extensive segregation distortion were observed in the population.

A goal of many genetic studies is to discover the underlying genetic condition (the genotype) of a specific physical manifestation in an organism (the phenotype), such as diabetes in humans or leaf rust in cultivated wheat. A limitation to making such discoveries is the ability to resolve genotype. Gene arrays carry representations of the genome, called features, at high-density on a surface the size of a thumbnail. In this study, microarrays designed to measure gene expression were used to detect DNA sequence polymorphisms. DNA from two different Arabidopsis strains was hybridized to arrays representing nearly the entire coding region of the genome. Differences in hybridization intensity indicated differences in DNA sequence. The sequence differences, termed single feature polymorphisms, were then assayed in a population of 100 plants derived through inbreeding the progeny from the two parental strains. The precise location of the genetic recombination breakpoints was defined for each line. As a result, Singer et al. were able to generate one of the first very high-resolution genotyping data sets in a multicellular organism that allowed the construction of a high-resolution genetic map of Arabidopsis. This map will greatly facilitate attempts to make definitive associations between genotypes and phenotypes.

The molecular genetic basis of plant adaptation

How natural selection on adaptive traits is filtered to the genetic level remains largely unknown. Theory and quantitative trait locus (QTL) mapping have provided insights into the number and effect of genes underlying adaptations, but these results have been hampered by questions of applicability to real biological systems and poor resolution, respectively. Advances in molecular technologies have expedited the cloning of adaptive genes through both forward and reverse genetic approaches. Forward approaches start with adaptive traits and attempt to characterize their underlying genetic architectures through linkage disequilibrium mapping, QTL mapping, and other methods. Reverse screens search large sequence data sets for genes that possess the signature of selection. Though both approaches have been successful in identifying adaptive genes in plants, very few, if any, of these adaptations' molecular bases have been fully resolved. The continued isolation of plant adaptive genes will lead to a more comprehensive understanding of natural selection's effect on genes and genomes.

Estimating the contribution of mutation, recombination and gene conversion in the generation of haplotypic diversity

Recombination occurs through both homologous crossing over and homologous gene conversion during meiosis. The contribution of recombination relative to mutation is expected to be dramatically reduced in inbreeding organisms. We report coalescent-based estimates of the recombination parameter (rho) relative to estimates of the mutation parameter (theta) for 18 genes from the highly self-fertilizing grass, wild barley, Hordeum vulgare ssp. spontaneum. Estimates of rho/theta are much greater than expected, with a mean rho/theta approximately 1.5, similar to estimates from outcrossing species. We also estimate rho with and without the contribution of gene conversion. Genotyping errors can mimic the effect of gene conversion, upwardly biasing estimates of the role of conversion. Thus we report a novel method for identifying genotyping errors in nucleotide sequence data sets. We show that there is evidence for gene conversion in many large nucleotide sequence data sets including our data that have been purged of all detectable sequencing errors and in data sets from Drosophila melanogaster, D. simulans, and Zea mays. In total, 13 of 27 loci show evidence of gene conversion. For these loci, gene conversion is estimated to contribute an average of twice as much as crossing over to total recombination.

Distribution of genetic variation within and among local populations of Arabidopsis thaliana over its species range

A detailed description of local population structure in Arabidopsis thaliana is presented, including an assessment of the genetic relatedness of individuals collected from the same field. A hierarchical sample of four individuals from 37 local populations, including North America, England, Eastern and Western Europe, and Asia, and a selection of ecotypes, were analysed for variation in Adh, ChiA, FAH1, F3H, Rpm1, Rps5, and five microsatellite loci. Twenty-eight of the 37 population samples contained individuals with identical multilocus haplotypes, 12 of which were fixed for a single haplotype. These monomorphic populations were evenly distributed over the species range. Only in North America did we find a single multilocus haplotype shared among different populations, perhaps indicating a continental founder event. Despite the occurrence of local inbreeding, a considerable amount of genetic variation was found segregating within and among local populations. A novel analysis of haplotype differences reveals that genetic differentiation occurs at every geographic scale in A. thaliana, where we find a surprising under-representation of recent migrants between local populations. This leads us to hypothesize that most dispersal between A. thaliana populations is by pollen rather than seed. Based on the structure of A. thaliana populations, it appears that regional groups of local populations may provide the most appropriate genetic material for linkage disequilibrium mapping of adaptive traits.

Moving forward in reverse: genetic technologies to enable genome-wide phenomic screens in Arabidopsis

Genome sequencing, in combination with various computational and empirical approaches to sequence annotation, has made possible the identification of more than 30,000 genes in Arabidopsis thaliana. Increasingly sophisticated genetic tools are being developed with the long-term goal of understanding how the coordinated activity of these genes gives rise to a complex organism. The combination of classical forward genetics with recently developed genome-wide, gene-indexed mutant collections is beginning to revolutionize the way in which gene functions are studied in plants. High-throughput screens using these mutant populations should provide a means to analyse plant gene functions--the phenome--on a genomic scale.

Association mapping with single-feature polymorphisms

We develop methods for exploiting "single-feature polymorphism" data, generated by hybridizing genomic DNA to oligonucleotide expression arrays. Our methods enable the use of such data, which can be regarded as very high density, but imperfect, polymorphism data, for genomewide association or linkage disequilibrium mapping. We use a simulation-based power study to conclude that our methods should have good power for organisms like Arabidopsis thaliana, in which linkage disequilibrium is extensive, the reason being that the noisiness of single-feature polymorphism data is more than compensated for by their great number. Finally, we show how power depends on the accuracy with which single-feature polymorphisms are called.

Challenges of detecting directional selection after a bottleneck: lessons from Sorghum bicolor

Multilocus surveys of sequence variation can be used to identify targets of directional selection, which are expected to have reduced levels of variation. Following a population bottleneck, the signal of directional selection may be hard to detect because many loci may have low variation by chance and the frequency spectrum of variation may be perturbed in ways that resemble the effects of selection. Cultivated Sorghum bicolor contains a subset of the genetic diversity found in its wild ancestor(s) due to the combined effects of a domestication bottleneck and human selection on traits associated with agriculture. As a framework for distinguishing between the effects of demography and selection, we sequenced 204 loci in a diverse panel of 17 cultivated S. bicolor accessions. Genomewide patterns of diversity depart strongly from equilibrium expectations with regard to the variance of the number of segregating sites, the site frequency spectrum, and haplotype configuration. Furthermore, gene genealogies of most loci with an excess of low frequency variants and/or an excess of segregating sites do not show the characteristic signatures of directional and diversifying selection, respectively. A simple bottleneck model provides an improved but inadequate fit to the data, suggesting the action of other population-level factors, such as population structure and migration. Despite a known history of recent selection, we find little evidence for directional selection, likely due to low statistical power and lack of an appropriate null model.

Genetic mechanisms and evolutionary significance of natural variation in Arabidopsis

Genomic studies of natural variation in model organisms provide a bridge between molecular analyses of gene function and evolutionary investigations of adaptation and natural selection. In the model plant species Arabidopsis thaliana, recent studies of natural variation have led to the identification of genes underlying ecologically important complex traits, and provided new insights about the processes of genome evolution, geographic population structure, and the selective mechanisms shaping complex trait variation in natural populations. These advances illustrate the potential for a new synthesis to elucidate mechanisms for the adaptive evolution of complex traits from nucleotide sequences to real-world environments.