Haplotypic divergence coupled with lack of diversity at the Arabidopsis thaliana alcohol dehydrogenase locus: roles for both balancing and directional selection?

Insertions/Deletions-Associated Nucleotide Polymorphism in Arabidopsis thaliana

Although high levels of within-species variation are commonly observed, a general mechanism for the origin of such variation is still lacking. Insertions and deletions (indels) are a widespread feature of genomes and we hypothesize that there might be an association between indels and patterns of nucleotide polymorphism. Here, we investigate flanking sequences around 18 indels (>100 bp) among a large number of accessions of the plant, Arabidopsis thaliana. We found two distinct haplotypes, i.e., a nucleotide dimorphism, present around each of these indels and dimorphic haplotypes always corresponded to the indel-present/-absent patterns. In addition, the peaks of nucleotide diversity between the two divergent alleles were closely associated with these indels. Thus, there exists a close association between indels and dimorphisms. Further analysis suggests that indel-associated substitutions could be an important component of genetic variation shaping nucleotide polymorphism in Arabidopsis. Finally, we suggest a mechanism by which indels might generate these highly divergent haplotypes. This study provides evidence that nucleotide dimorphisms, which are frequently regarded as evidence of frequency-dependent selection, could be explained simply by structural variation in the genome.

1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana

Arabidopsis thaliana serves as a model organism for the study of fundamental physiological, cellular, and molecular processes. It has also greatly advanced our understanding of intraspecific genome variation. We present a detailed map of variation in 1,135 high-quality re-sequenced natural inbred lines representing the native Eurasian and North African range and recently colonized North America. We identify relict populations that continue to inhabit ancestral habitats, primarily in the Iberian Peninsula. They have mixed with a lineage that has spread to northern latitudes from an unknown glacial refugium and is now found in a much broader spectrum of habitats. Insights into the history of the species and the fine-scale distribution of genetic diversity provide the basis for full exploitation of A. thaliana natural variation through integration of genomes and epigenomes with molecular and non-molecular phenotypes.

Genetic diversity and population structure of Arabidopsis thaliana along an altitudinal gradient

The natural genetic variation within a plant species is primarily a consequence of its phylogeography and evolutionary history. This variation largely determines its present-day population structure. Arabidopsis thaliana, as a model plant, has been studied in great detail including its probable origin, local as well as global genetic diversity pattern, population structure, adaptation, etc. However, no such studies have so far been reported from the Indian Himalayan region. Here, we describe a comprehensive study on the genetic diversity and population structure of A. thaliana from an altitudinal range of 700-3400 m above mean sea level the highest altitudinal range reported so far. We also compare these populations with previously reported worldwide populations. A total of 48 accessions representing six populations were analysed using 19 microsatellites and 11 chloroplast markers. Genetic diversity analysis indicated populations to be highly diverse and comparable with worldwide populations. STRUCTURE, principal coordinate and isolation by distance (IBD) analyses showed that genetic variation in different populations is structured at geographical and altitudinal level. Further analyses indicate that these populations are genetically distinct from the rest of the world populations. Different parameters of the demographic expansion model support a rapid expansion. Based on mismatch distribution, the initial time of expansion of west Himalayan populations was found to be about 130 000 years. Bayesian analysis of divergence time indicated that these populations have a long evolutionary history in this region. Based on the results of genetic diversity parameters, demographic expansion and divergence time estimation, it appears that west Himalayan populations may be the source of the west-east expansion model.

A haploid genetics toolbox for Arabidopsis thaliana

Genetic analysis in haploids provides unconventional yet powerful advantages not available in diploid organisms. In Arabidopsis thaliana, haploids can be generated through seeds by crossing a wild-type strain to a transgenic strain with altered centromeres. Here we report the development of an improved haploid inducer (HI) strain, SeedGFP-HI, that aids selection of haploid seeds prior to germination. We also show that haploids can be used as a tool to accelerate a variety of genetic analyses, specifically pyramiding multiple mutant combinations, forward mutagenesis screens, scaling down a tetraploid to lower ploidy levels and swapping of nuclear and cytoplasmic genomes. Furthermore, the A. thaliana HI can be used to produce haploids from a related species A. suecica and generate homozygous mutant plants from strong maternal gametophyte lethal alleles, which is not possible via conventional diploid genetics. Taken together, our results demonstrate the utility and power of haploid genetics in A. thaliana.

A phylogenetic analysis of the genus Fragaria (strawberry) using intron-containing sequence from the ADH-1 gene

The genus Fragaria encompasses species at ploidy levels ranging from diploid to decaploid. The cultivated strawberry, Fragaria×ananassa, and its two immediate progenitors, F. chiloensis and F. virginiana, are octoploids. To elucidate the ancestries of these octoploid species, we performed a phylogenetic analysis using intron-containing sequences of the nuclear ADH-1 gene from 39 germplasm accessions representing nineteen Fragaria species and one outgroup species, Dasiphora fruticosa. All trees from Maximum Parsimony and Maximum Likelihood analyses showed two major clades, Clade A and Clade B. Each of the sampled octoploids contributed alleles to both major clades. All octoploid-derived alleles in Clade A clustered with alleles of diploid F. vesca, with the exception of one octoploid allele that clustered with the alleles of diploid F. mandshurica. All octoploid-derived alleles in clade B clustered with the alleles of only one diploid species, F. iinumae. When gaps encoded as binary characters were included in the Maximum Parsimony analysis, tree resolution was improved with the addition of six nodes, and the bootstrap support was generally higher, rising above the 50% threshold for an additional nine branches. These results, coupled with the congruence of the sequence data and the coded gap data, validate and encourage the employment of sequence sets containing gaps for phylogenetic analysis. Our phylogenetic conclusions, based upon sequence data from the ADH-1 gene located on F. vesca linkage group II, complement and generally agree with those obtained from analyses of protein-encoding genes GBSSI-2 and DHAR located on F. vesca linkage groups V and VII, respectively, but differ from a previous study that utilized rDNA sequences and did not detect the ancestral role of F. iinumae.

The BOY NAMED SUE quantitative trait locus confers increased meiotic stability to an adapted natural allopolyploid of Arabidopsis

Whole-genome duplication resulting from polyploidy is ubiquitous in the evolutionary history of plant species. Yet, polyploids must overcome the meiotic challenge of pairing, recombining, and segregating more than two sets of chromosomes. Using genomic sequencing of synthetic and natural allopolyploids of Arabidopsis thaliana and Arabidopsis arenosa, we determined that dosage variation and chromosomal translocations consistent with homoeologous pairing were more frequent in the synthetic allopolyploids. To test the role of structural chromosomal differentiation versus genetic regulation of meiotic pairing, we performed sequenced-based, high-density genetic mapping in F2 hybrids between synthetic and natural lines. This F2 population displayed frequent dosage variation and deleterious homoeologous recombination. The genetic map derived from this population provided no indication of structural evolution of the genome of the natural allopolyploid Arabidopsis suecica, compared with its predicted parents. The F2 population displayed variation in meiotic regularity and pollen viability that correlated with a single quantitative trait locus, which we named BOY NAMED SUE, and whose beneficial allele was contributed by A. suecica. This demonstrates that an additive, gain-of-function allele contributes to meiotic stability and fertility in a recently established allopolyploid and provides an Arabidopsis system to decipher evolutionary and molecular mechanisms of meiotic regularity in polyploids.

Genetic and epigenetic changes in a genomic region containing MIR172 in Arabidopsis allopolyploids and their progenitors

Combination of divergent genomes in allopolyploids creates genome-wide gene expression changes including many miRNA targets in Arabidopsis allotetraploids relative to the parents Arabidopsis thaliana and Arabidopsis arenosa. Here we report expression and epigenetic changes in a chromosomal region containing the MIR172b locus in the allotetraploids. Although mature miRNA sequences are conserved, A. thaliana and A. arenosa miRNA loci diverge rapidly in sequence and expression. Among four MIR172 loci in Arabidopsis, the level of nucleotide sequence divergence between A. thaliana and A. arenosa MIR172 loci is 15-25%, which is higher than that of protein-coding genes (∼5%). MIR172b locus and its flanking genes in A. arenosa were expressed at low levels relative to that in A. thaliana, which is associated with hypermethylation of this region in the allotetraploids. Consistently with this notion, pri-miR172 transcripts in the allotetraploids were primarily derived from the A. thaliana MIR172b locus. Expression of homoeologous alleles in miR172 target loci is associated with allelic loss, allelic changes in outcrossing A. arenosa or repression of A. thaliana alleles. These data suggest that gene expression changes in this homoeologous region are associated with genetic diversity and epigenetic variation of miRNA genes and their targets in allopolyploids.

Complex interplay of evolutionary forces in the ladybird homeobox genes of Drosophila melanogaster

Tandemly arranged paralogous genes lbe and lbl are members of the Drosophila NK homeobox family. We analyzed population samples of Drosophila melanogaster from Africa, Europe, North and South America, and single strains of D. sechellia, D. simulans, and D. yakuba within two linked regions encompassing partial sequences of lbe and lbl. The evolution of lbe and lbl is highly constrained due to their important regulatory functions. Despite this, a variety of forces have shaped the patterns of variation in lb genes: recombination, intragenic gene conversion and natural selection strongly influence background variation created by linkage disequilibrium and dimorphic haplotype structure. The two genes exhibited similar levels of nucleotide diversity and positive selection was detected in the noncoding regions of both genes. However, synonymous variability was significantly higher for lbe: no nonsynonymous changes were observed in this gene. We argue that balancing selection impacts some synonymous sites of the lbe gene. Stability of mRNA secondary structure was significantly different between the lbe (but not lbl) haplotype groups and may represent a driving force of balancing selection in epistatically interacting synonymous sites. Balancing selection on synonymous sites may be the first, or one of a few such observations, in Drosophila. In contrast, recurrent positive selection on lbl at the protein level influenced evolution at three codon sites. Transcription factor binding-site profiles were different for lbe and lbl, suggesting that their developmental functions are not redundant. Combined with our previous results on nucleotide variation in esterase and other homeobox genes, these results suggest that interplay of balancing and directional selection may be a general feature of molecular evolution in Drosophila and other eukaryote genomes.

The origin of populations of Arabidopsis thaliana in China, based on the chloroplast DNA sequences

BACKGROUND

In the studies incorporating worldwide sampling of A. thaliana populations, the samples from East Asia, especially from China, were very scattered; and the studies focused on global patterns of cpDNA genetic variation among accessions of A. thaliana are very few. In this study, chloroplast DNA sequence variability was used to infer phylogenetic relationships among Arabidopsis thaliana accessions from around the world, with the emphasis on samples from China.

RESULTS

A data set comprising 77 accessions of A. thaliana, including 19 field-collected Chinese accessions together with three related species (A. arenosa, A. suecica, and Olimarabidopsis cabulica) as the out-group, was compiled. The analysis of the nucleotide sequences showed that the 77 accessions of A. thaliana were partitioned into two major differentiated haplotype classes (MDHCs). The estimated divergence time of the two MDHCs was about 0.39 mya. Forty-nine haplotypes were detected among the 77 accessions, which exhibited nucleotide diversity (pi) of 0.00169. The Chinese populations along the Yangtze River were characterized by five haplotypes, and the two accessions collected from the middle range of the Altai Mountains in China shared six specific variable sites.

CONCLUSIONS

The dimorphism in the chloroplast DNA could be due to founder effects during late Pleistocene glaciations and interglacial periods, although introgression cannot be ruled out. The Chinese populations along the Yangtze River may have dispersed eastwards to their present-day locations from the Himalayas. These populations originated from a common ancestor, and a rapid demographic expansion began approximately 90,000 years ago. Two accessions collected from the middle range of the Altai Mountains in China may have survived in a local refugium during late Pleistocene glaciations. The natural populations from China with specific genetic characteristics enriched the gene pools of global A. thaliana collections.

Grouped nucleotide polymorphism: a major contributor to genetic variation in Arabidopsis

While information for single nucleotide polymorphism is accumulating in many organisms, little is known about the magnitude and the occurrence of nucleotide dimorphism or higher order structural polymorphisms (grouped nucleotide polymorphisms or GNPs). To address these questions, we systematically investigated the genetic variations of 996 loci in 96 Arabidopsis accessions. Our data suggest that GNP loci are highly frequent in the genomes (40.4% of 996 loci), and that 66.7% of genetic variation in 996 loci is attributed to GNPs. The frequency distribution and the linkage analysis of GNP loci demonstrate that GNPs occurred randomly, locally, ceaselessly and independently in general, although they are particularly abundant in the region near centromeres. The age distribution of GNP loci shows that the majority of these loci are under a transient phase of neutral evolution. The characteristics of GNPs imply that a molecular isolation exists in GNP loci, and most likely the indels caused isolation, which could explain how GNPs are generated and maintained.

Nucleotide polymorphism in the drought responsive gene Asr2 in wild populations of tomato

The Asr gene family (named after abscicic acid [ABA], stress, ripening), exclusively present in plant genomes, is involved in transcriptional regulation. Its members are up-regulated in roots and leaves of water- or salt-stressed plants. In previous work, evidence of adaptive evolution (as inferred from synonymous and nonsynonymous divergence rates) has been reported for Asr2 in Solanum chilense and S. arcanum, two species dwelling in habitats with different precipitation regimes. In this paper we investigate patterns of intraspecific nucleotide variation in Asr2 and the unlinked locus CT114 in S. chilense and S. arcanum. The extent of nucleotide diversity in Asr2 differed between species in more than one order of magnitude. In both species we detected evidence of non-neutral evolution, which may be ascribed to different selective regimes, potentially associated to unique climatic features, or, alternatively, to demographic events. The results are discussed in the light of demographic and selective hypotheses.

RNAi of met1 reduces DNA methylation and induces genome-specific changes in gene expression and centromeric small RNA accumulation in Arabidopsis allopolyploids

Changes in genome structure and gene expression have been documented in both resynthesized and natural allopolyploids that contain two or more divergent genomes. The underlying mechanisms for rapid and stochastic changes in gene expression are unknown. Arabidopsis suecica is a natural allotetraploid derived from the extant A. thaliana and A. arenosa genomes that are homeologous in the allotetraploid. Here we report that RNAi of met1 reduced DNA methylation and altered the expression of approximately 200 genes, many of which encode transposons, predicted proteins, and centromeric and heterochromatic RNAs. Reduced DNA methylation occurred frequently in promoter regions of the upregulated genes, and an En/Spm-like transposon was reactivated in met1-RNAi A. suecica lines. Derepression of transposons, heterochromatic repeats, and centromeric small RNAs was primarily derived from the A. thaliana genome, and A. arenosa homeologous loci were less affected by methylation defects. A high level of A. thaliana centromeric small RNA accumulation was correlated with hypermethylation of A. thaliana centromeres. The greater effects of reduced DNA methylation on transposons and centromeric repeats in A. thaliana than in A. arenosa are consistent with the repression of many genes that are expressed at higher levels in A. thaliana than in A. arenosa in the resynthesized allotetraploids. Moreover, non-CG (CC) methylation in the promoter region of A. thaliana At2g23810 remained in the resynthesized allotetraploids, and the methylation spread within the promoter region in natural A. suecica, leading to silencing of At2g23810. At2g23810 was demethylated and reactivated in met1-RNAi A. suecica lines. We suggest that many A. thaliana genes are transcriptionally repressed in resynthesized allotetraploids, and a subset of A. thaliana loci including transposons and centromeric repeats are heavily methylated and subjected to homeologous genome-specific RNA-mediated DNA methylation in natural allopolyploids.

Nucleotide sequence variability of the Adh gene of the coastal plant Calystegia soldanella (Convolvulaceae) in Japan

Calystegia soldanella (Convolvulaceae) is a self-incompatible perennial herb distributed on sandy seashores throughout the temperate zone of the world. In Japan, the species occasionally grows on the sandy shores of Lake Biwa. To clarify the genetic differentiation among local populations, we investigated the nucleotide sequence variability of the Adh gene. In a 1625-bp sequence between exon 2 and the 3' noncoding region of the Adh gene, a total of 44 polymorphic sites were found among 91 individuals from 19 populations. The nucleotide diversity for the entire sample was 0.00212. Similar values were determined for geographical groups of populations. No genetic differentiation among the groups of populations was found. The complete lack of genetic differentiation between the sea coastal populations and the inland populations could not be attributed to gene flow. Although the inland populations are geographically isolated from the sea coastal populations, the time since separation might be insufficient to establish significant genetic differentiation.

Unique nucleotide polymorphism of ankyrin gene cluster in Arabidopsis

The ankyrin (ANK) gene cluster is a part of a multigene family encoding ANK transmembrane proteins in Arabidopsis thaliana, and plays an important role in protein-protein interactions and in signal pathways. In contrast to other regions of a genome, the ANK gene cluster exhibits an extremely high level of DNA polymorphism in an approximately 5-kb region, without apparent decay. Phylogenetic analysis detects two clear, deeply differentiated haplotypes (dimorphism). The divergence between haplotypes of accession Col-0 and Ler-0 (Hap-C and Hap-L) is estimated to be 10.7%, approximately equal to the 10.5% average divergence between A. thaliana and A. lyrata. Sequence comparisons for the ANK gene cluster homologues in Col-0 indicate that the members evolve independently, and that the similarity among paralogues is lower than between alleles. Very little intralocus recombination or gene conversion is detected in ANK regions. All these characteristics of the ANK gene cluster are consistent with a tandem gene duplication and birth-and-death process. The possible mechanisms for and implications of this elevated nucleotide variation are also discussed, including the suggestion of balancing selection.

Genetic diversity of the natural populations of Arabidopsis thaliana in China

Although extensive studies have been conducted on the genetic structure of Arabidopsis thaliana (A. thaliana) populations worldwide, the populations from China have never been studied. In this study, we collected 560 individuals from 19 natural populations of A. thaliana distributed in East China along the lower reaches of the Yangtze River, and two populations from northwest China (Xinjiang Province). We adopted two kinds of molecular marker, inter-simple sequence repeats (ISSRs) and random amplified polymorphic DNA (RAPDs) to investigate the genetic diversity within and among populations, and the correlation between the genetic and geographic distances. Thirteen ISSR primers produced 165 polymorphic bands (PPB) (96%) and 11 RAPD primers produced 162 polymorphic bands (98%) in about 560 individuals. The two marker systems generated similar patterns of genetic diversity in these natural populations. The AMOVA analysis indicated about 42-45% of the total genetic variation existed within populations, and found possible geographic structure. The Mantel test revealed a significant correlation between the geographic distance and the genetic distance of these populations in general. A close genetic relationship was found among four populations in the Jiangxi Province, and these always appeared clustered together as a monophyletic group in unweighted pair-group method with arithmetic averages dendrograms based on both ISSR and RAPD data sets. Based on the observation of recolonization and extinction of naturally distributed populations of A. thaliana, and the pattern of their genetic differentiation, the distribution of this species in China might be a result of natural dispersal under the strong influence of human activity.

Sequence diversity and haplotype associations with phenotypic responses to crowding: GIGANTEA affects fruit set in Arabidopsis thaliana

Identifying the molecular genetic basis of intraspecific variation in quantitative traits promises to provide novel insight into their evolutionary history as well as genetic mechanisms of adaptation. In an attempt to identify genes responsible for natural variation in competitive responses in Arabidopsis thaliana, we examined DNA sequence diversity at seven loci previously identified as members of the phytochrome B signalling network. For one gene, GIGANTEA (GI), we detected significant haplotype structure. To test for GI haplogroup-phenotype associations, we genotyped 161 A. thaliana accessions at GI and censused the same accessions for total fruit set and the expression of three phenotypic traits (days to flowering, petiole length, and inflorescence height) in a greenhouse experiment where plants were grown in crowded and uncrowded environments. We detected a significant association between GI and total fruit set that resulted in a 14% difference in average fruit set among GI haplogroups. Given that fruit set is an important component of fitness in this species and given the magnitude of the effect, the question arises as to how variation at this locus is maintained. Our observation of frequent and significant epistasis between GI and background single nucleotide polymorphisms (SNP), where the fitness ranking of the GI allele either reverses or does not differ depending on the allele at the interacting SNP, suggests that epistatic selection may actively maintain or at least slow the loss of variation at GI. This result is particularly noteworthy in the light of the ongoing debate regarding the genetic underpinnings of phenotypic evolution and recent observations that epistasis for phenotypic traits and components of fitness is common in A. thaliana.

Studying genetics of adaptive variation in model organisms: flowering time variation in Arabidopsis lyrata

Arabidopsis thaliana has emerged as a model organism for plant developmental genetics, but it is also now being widely used for population genetic studies. Outcrossing relatives of A. thaliana are likely to provide suitable additional or alternative species for studies of evolutionary and population genetics. We have examined patterns of adaptive flowering time variation in the outcrossing, perennial A. lyrata. In addition, we examine the distribution of variation at marker genes in populations form North America and Europe. The probability of flowering in this species differs between southern and northern populations. Northern populations are much less likely to flower in short than in long days. A significant daylength by region interaction shows that the northern and southern populations respond differently to the daylength. The timing of flowering also differs between populations, and is made shorter by long days, and in some populations, by vernalization. North American and European populations show consistent genetic differentiation over microsatellite and isozyme loci and alcohol dehydrogenase sequences. Thus, the patterns of variation are quite different from those in A. thaliana, where flowering time differences show little relationship to latitude of origin and the genealogical trees of accessions vary depending on the genomic region studied. The genetic architecture of adaptation can be compared in these species with different life histories.

The pattern of polymorphism in Arabidopsis thaliana

We resequenced 876 short fragments in a sample of 96 individuals of Arabidopsis thaliana that included stock center accessions as well as a hierarchical sample from natural populations. Although A. thaliana is a selfing weed, the pattern of polymorphism in general agrees with what is expected for a widely distributed, sexually reproducing species. Linkage disequilibrium decays rapidly, within 50 kb. Variation is shared worldwide, although population structure and isolation by distance are evident. The data fail to fit standard neutral models in several ways. There is a genome-wide excess of rare alleles, at least partially due to selection. There is too much variation between genomic regions in the level of polymorphism. The local level of polymorphism is negatively correlated with gene density and positively correlated with segmental duplications. Because the data do not fit theoretical null distributions, attempts to infer natural selection from polymorphism data will require genome-wide surveys of polymorphism in order to identify anomalous regions. Despite this, our data support the utility of A. thaliana as a model for evolutionary functional genomics.

A systematic global survey of genomic DNA sequence polymorphism in Arabidopsis thaliana reveals that standard genetic tests for selection do not apply to this species but supports its status as a model organism.

Linkage disequilibrium and association studies in higher plants: present status and future prospects

During the last two decades, DNA-based molecular markers have been extensively utilized for a variety of studies in both plant and animal systems. One of the major uses of these markers is the construction of genome-wide molecular maps and the genetic analysis of simple and complex traits. However, these studies are generally based on linkage analysis in mapping populations, thus placing serious limitations in using molecular markers for genetic analysis in a variety of plant systems. Therefore, alternative approaches have been suggested, and one of these approaches makes use of linkage disequilibrium (LD)-based association analysis. Although this approach of association analysis has already been used for studies on genetics of complex traits (including different diseases) in humans, its use in plants has just started. In the present review, we first define and distinguish between LD and association mapping, and then briefly describe various measures of LD and the two methods of its depiction. We then give a list of different factors that affect LD without discussing them, and also discuss the current issues of LD research in plants. Later, we also describe the various uses of LD in plant genomics research and summarize the present status of LD research in different plant genomes. In the end, we discuss briefly the future prospects of LD research in plants, and give a list of softwares that are useful in LD research, which is available as electronic supplementary material (ESM).

A multilocus sequence survey in Arabidopsis thaliana reveals a genome-wide departure from a neutral model of DNA sequence polymorphism

The simultaneous analysis of multiple genomic loci is a powerful approach to studying the effects of population history and natural selection on patterns of genetic variation of a species. By surveying nucleotide sequence polymorphism at 334 randomly distributed genomic regions in 12 accessions of Arabidopsis thaliana, we examined whether a standard neutral model of nucleotide sequence polymorphism is consistent with observed data. The average nucleotide diversity was 0.0071 for total sites and 0.0083 for silent sites. Although levels of diversity are variable among loci, no correlation with local recombination rate was observed, but polymorphism levels were correlated for physically linked loci (<250 kb). We found that observed distributions of Tajima's D- and D/D(min)- and of Fu and Li's D-, D*- and F-, F*-statistics differed significantly from the expected distributions under a standard neutral model due to an excess of rare polymorphisms and high variances. Observed and expected distributions of Fay and Wu's H were not different, suggesting that demographic processes and not selection at multiple loci are responsible for the deviation from a neutral model. Maximum-likelihood comparisons of alternative demographic models like logistic population growth, glacial refugia, or past bottlenecks did not produce parameter estimates that were more consistent with observed patterns. However, exclusion of highly polymorphic "outlier loci" resulted in a fit to the logistic growth model. Various tests of neutrality revealed a set of candidate loci that may evolve under selection.

Nucleotide variation in the tinman and bagpipe homeobox genes of Drosophila melanogaster

The tinman (tin) and bagpipe (bap) genes are members of the NK homeobox gene family of Drosophila, so that tin occupies a higher position than bap in the regulatory hierarchy. Little is known about the level and pattern of genetic polymorphism in homeobox genes. We have analyzed nucleotide polymorphism in 27 strains of Drosophila melanogaster and one each of D. simulans and D. sechellia, within two closely linked regions encompassing a partial sequence of tin and the complete sequence of bap. The two genes exhibit different levels and patterns of nucleotide diversity. Two sets of sharply divergent sequence types are detected for tin. The haplotype structure of bap is more complex: about half of the sequences are identical (or virtually so), while the rest are fairly heterogeneous. The level of silent nucleotide variability is 0.0063 for tin but significantly higher, 0.0141, for bap, a level of polymorphism comparable to the most polymorphic structural genes of D. melanogaster. Recombination rate and gene conversion are also higher for bap than for tin. There is strong linkage disequilibrium, with the highest values in the introns of both genes and exon II of bap. The patterns of polymorphism in tin and bap are not compatible with an equilibrium model of selective neutrality. We suggest that negative selection and demographic history are the major factors shaping the pattern of nucleotide polymorphism in the tin and bap genes; moreover, there are clear indications of positive selection in the bap gene.

The maintenance of extreme amino acid diversity at the disease resistance gene, RPP13, in Arabidopsis thaliana

We have used the naturally occurring plant-parasite system of Arabidopsis thaliana and its common parasite Peronospora parasitica (downy mildew) to study the evolution of resistance specificity in the host population. DNA sequence of the resistance gene, RPP13, from 24 accessions, including 20 from the United Kingdom, revealed amino acid sequence diversity higher than that of any protein coding gene reported so far in A. thaliana. A significant excess of amino acid polymorphism segregating within this species is localized within the leucine-rich repeat (LRR) domain of RPP13. These results indicate that single alleles of the gene have not swept through the population, but instead, a diverse collection of alleles have been maintained. Transgenic complementation experiments demonstrate functional differences among alleles in their resistance to various pathogen isolates, suggesting that the extreme amino acid polymorphism in RPP13 is maintained through continual reciprocal selection between host and pathogen.

Functional divergence in tandemly duplicated Arabidopsis thaliana trypsin inhibitor genes

In multigene families, variation among loci and alleles can contribute to trait evolution. We explored patterns of functional and genetic variation in six duplicated Arabidopsis thaliana trypsin inhibitor (ATTI) loci. We demonstrate significant variation in constitutive and herbivore-induced transcription among ATTI loci that show, on average, 65% sequence divergence. Significant variation in ATTI expression was also found between two molecularly defined haplotype classes. Population genetic analyses for 17 accessions of A. thaliana showed that six ATTI loci arranged in tandem within 10 kb varied 10-fold in nucleotide diversity, from 0.0009 to 0.0110, and identified a minimum of six recombination events throughout the tandem array. We observed a significant peak in nucleotide and indel polymorphism spanning ATTI loci in the interior of the array, due primarily to divergence between the two haplotype classes. Significant deviation from the neutral equilibrium model for individual genes was interpreted within the context of intergene linkage disequilibrium and correlated patterns of functional differentiation. In contrast to the outcrosser Arabidopsis lyrata for which recombination is observed even within ATTI loci, our data suggest that response to selection was slowed in the inbreeding, annual A. thaliana because of interference among functionally divergent ATTI loci.

Molecular Population Genetics and the Search for Adaptive Evolution in Plants

The first papers on plant molecular population genetics were published approximately 10 years ago. Since that time, well over 50 additional studies of plant nucleotide polymorphism have been published, and many of these studies focused on detecting the signature of balancing or positive selection at a locus. In this review, we discuss some of the theoretical and statistical issues surrounding the detection of selection, with focus on plant populations, and we also summarize the empirical plant molecular population genetics literature. At face value, the literature suggests that a history of balancing or positive selection in plant genes is rampant. In two well-studied taxa (maize and Arabidopsis) over 20% of studied genes have been interpreted as containing the signature of selection. We argue that this is probably an overstatement of the prevalence of natural selection in plant genomes, for two reasons. First, demographic effects are difficult to incorporate and have generally not been well integrated into the plant population genetics literature. Second, the genes studied to date are not a random sample, so selected genes may be overrepresented. The next generation of studies in plant molecular population genetics requires additional sampling of local populations, explicit comparisons among loci, and improved theoretical methods to control for demography. Eventually, candidate loci should be confirmed by explicit consideration of phenotypic effects.

Neutral genetic variation among wild North American populations of the weedy plant Arabidopsis thaliana is not geographically structured

We investigated neutral genetic variation within and among 53 wild-collected populations of the weedy annual plant, Arabidopsis thaliana, in North America, using amplified fragment length polymorphism (AFLP) markers. A. thaliana is thought to have been introduced to North America from Eurasia by humans; such an introduction might be expected to leave a clear geographical signal in the genetic data. To detect such patterns, we sampled populations at several hierarchical geographical levels. We collected individuals from populations in two areas of the Southeast and one in the Midwest, as well as individuals from populations in the Pacific Northwest and Northeast. To estimate within-population variation, we sampled eight individuals from each of six populations in the Southeast and Midwest. Among all 95 individuals analysed, we detected 131 polymorphic AFLP fragments. We found no evidence for continental or regional diversification. Individuals sampled from Midwestern and Southeastern populations intermingled in a neighbour-joining tree, and Mantel tests conducted within the Midwestern and Southeastern regions as well as the full data set failed to detect any significant relationship between geographical and genetic distance. These results mirror those found for most global surveys of neutral genetic variability in A. thaliana. Surprisingly, we detected substantial amounts of neutral genetic variability within populations. The levels of genetic variation within populations, coupled with the nongeographical nature of divergence among populations, are consistent with contemporary gene flow and point to a complex and dynamic population history of A. thaliana in North America.

The development of an Arabidopsis model system for genome-wide analysis of polyploidy effects

Arabidopsis is a model system not only for studying numerous aspects of plant biology, but also for understanding mechanisms of the rapid evolutionary process associated with genome duplication and polyploidization. Although in animals interspecific hybrids are often sterile and aneuploids are related to disease syndromes, both Arabidopsis autopolyploids and allopolyploids occur in nature and can be readily formed in the laboratory, providing an attractive system for comparing changes in gene expression and genome structure among relatively 'young' and 'established' or 'ancient' polyploids. Powerful reverse and forward genetics in Arabidopsis offer an exceptional means by which regulatory mechanisms of gene and genome duplication may be revealed. Moreover, the Arabidopsis genome is completely sequenced; both coding and non-coding sequences are available. We have developed spotted oligo-gene and chromosome microarrays using the complete Arabidopsis genome sequence. The oligo-gene microarray consists of ~26 000 70-mer oligonucleotides that are designed from all annotated genes in Arabidopsis, and the chromosome microarray contains 1 kb genomic tiling fragments amplified from a chromosomal region or the complete sequence of chromosome 4. We have demonstrated the utility of microarrays for genome-wide analysis of changes in gene expression, genome organization and chromatin structure in Arabidopsis polyploids and related species.

Nested core collections maximizing genetic diversity in Arabidopsis thaliana

The successful exploitation of natural genetic diversity requires a basic knowledge of the extent of the variation present in a species. To study natural variation in Arabidopsis thaliana, we defined nested core collections maximizing the diversity present among a worldwide set of 265 accessions. The core collections were generated based on DNA sequence data from a limited number of fragments evenly distributed in the genome and were shown to successfully capture the molecular diversity in other loci as well as the morphological diversity. The core collections are available to the scientific community and thus provide an important resource for the study of genetic variation and its functional consequences in Arabidopsis. Moreover, this strategy can be used in other species to provide a rational framework for undertaking diversity surveys, including single nucleotide polymorphism (SNP) discovery and phenotyping, allowing the utilization of genetic variation for the study of complex traits.

Nucleotide Variation in the tinman and bagpipe Homeobox Genes of Drosophila melanogaster

The tinman (tin) and bagpipe (bap) genes are members of the NK homeobox gene family of Drosophila, so that tin occupies a higher position than bap in the regulatory hierarchy. Little is known about the level and pattern of genetic polymorphism in homeobox genes. We have analyzed nucleotide polymorphism in 27 strains of Drosophila melanogaster and one each of D. simulans and D. sechellia, within two closely linked regions encompassing a partial sequence of tin and the complete sequence of bap. The two genes exhibit different levels and patterns of nucleotide diversity. Two sets of sharply divergent sequence types are detected for tin. The haplotype structure of bap is more complex: about half of the sequences are identical (or virtually so), while the rest are fairly heterogeneous. The level of silent nucleotide variability is 0.0063 for tin but significantly higher, 0.0141, for bap, a level of polymorphism comparable to the most polymorphic structural genes of D. melanogaster. Recombination rate and gene conversion are also higher for bap than for tin. There is strong linkage disequilibrium, with the highest values in the introns of both genes and exon II of bap. The patterns of polymorphism in tin and bap are not compatible with an equilibrium model of selective neutrality. We suggest that negative selection and demographic history are the major factors shaping the pattern of nucleotide polymorphism in the tin and bap genes; moreover, there are clear indications of positive selection in the bap gene.

Geographic distribution and recombination of genomic fragments on the short arm of chromosome 2 of Arabidopsis thaliana

Range expansion from Pleistocene refugia and anthropogenic influences contribute to the present distribution pattern of Arabidopsis thaliana. We scored a genome-wide set of CAPSs and found two markers with an east-west geographic distribution across the Eurasian range of the species. Regions around the two SNPs were sequenced in 98 accessions, including newly collected plants from Middle Asia and Western Siberia. These regions correspond to a gene ( approximately 1500 bp) and a non-coding region ( approximately 500 bp) 300 kbp apart on chromosome 2. Nucleotide diversities, pi, of the two sequenced fragments were 0.0032 and 0.0130. The haplotypes of both sequences belonged to one of two groups: a rather uniform "Asian" and a more variable "European" haplotype group, on the basis of non-disjunct clusters of SNPs. Recombination between "Asian" and "European" haplotypes occurs where they meet. Especially in the "European" haplotype, many rare SNP variants representing independent mutations are scattered among the shared haplotype-specific SNPs. This agrees with previous suggestions of two large haplotype groups in A. thaliana and the post-glacial colonization of central Europe from the east and the west. A clear correlation between climatic factors and the haplotype distribution may reflect the dispersal history rather than local climate adaptation. The pattern of SNP variation within the contiguous sequences explains why only a minority of SNPs selected across the genome show evidence of this geographic pattern.

Sensitivity of 70-mer oligonucleotides and cDNAs for microarray analysis of gene expression in Arabidopsis and its related species

Synthetic oligonucleotides (oligos) represent an attractive alternative to cDNA amplicons for spotted microarray analysis in a number of model organisms, including Arabidopsis, C. elegans, Drosophila, human, mouse and yeast. However, little is known about the relative effectiveness of 60-70-mer oligos and cDNAs for detecting gene expression changes. Using 192 pairs of Arabidopsis thaliana cDNAs and corresponding 70-mer oligos, we performed three sets of dye-swap experiments and used analysis of variance (anova) to compare sources of variation and sensitivities for detecting gene expression changes in A. thaliana, A. arenosa and Brassica oleracea. Our major findings were: (1) variation among different RNA preparations from the same tissue was small, but large variation among dye-labellings and slides indicates the need to replicate these factors; (2) sources of variation were similar for experiments with all three species, suggesting these feature types are effective for analysing gene expression in related species; (3) oligo and cDNA features had similar sensitivities for detecting expression changes and they identified a common subset of significant genes, but results from quantitative RT-PCR did not support the use of one over the other. These findings indicate that spotted oligos are at least as effective as cDNAs for microarray analyses of gene expression. We are using oligos designed from approximately 26,000 annotated genes of A. thaliana to study gene expression changes in Arabidopsis and Brassica polyploids.

Natural variation in nucleolar dominance reveals the relationship between nucleolus organizer chromatin topology and rRNA gene transcription in Arabidopsis

In genetic hybrids, nucleolus formation on chromosomes inherited from only one parent is the epigenetic phenomenon, nucleolar dominance. By using Arabidopsis suecica, the allotetraploid hybrid of Arabidopsis thaliana and Arabidopsis arenosa, natural variation in nucleolar dominance was found to occur, providing a unique opportunity to examine homologous nucleolus organizer regions (NORs) in their active and inactive states. In A. suecica strain LC1, NORs derived from A. arenosa are active, whereas A. thaliana-derived NORs are silenced. In A. suecica strain 9502, NORs of both parental species are active. When active, NORs are partially, but not fully, decondensed. Both active and inactive LC1 NORs colocalize with the nucleolus, contradicting the long-standing assumption that rRNA gene transcription drives nucleolus association. Collectively, these observations clarify the relationships among NOR chromatin topology, rRNA gene transcription, and NOR-nucleolus associations. A. suecica strains LC1 and 9502 have each lost one pair of A. thaliana NORs during evolution, and amplified fragment-length polymorphism analysis further indicates that these strains are genetically very similar. These data suggest that nucleolar dominance can result from subtle genetic or epigenetic variation but is not a trait fundamental to a given interspecies hybrid combination.

Single-nucleotide polymorphisms in soybean

Single-nucleotide polymorphisms (SNPs) provide an abundant source of DNA polymorphisms in a number of eukaryotic species. Information on the frequency, nature, and distribution of SNPs in plant genomes is limited. Thus, our objectives were (1) to determine SNP frequency in coding and noncoding soybean (Glycine max L. Merr.) DNA sequence amplified from genomic DNA using PCR primers designed to complete genes, cDNAs, and random genomic sequence; (2) to characterize haplotype variation in these sequences; and (3) to provide initial estimates of linkage disequilibrium (LD) in soybean. Approximately 28.7 kbp of coding sequence, 37.9 kbp of noncoding perigenic DNA, and 9.7 kbp of random noncoding genomic DNA were sequenced in each of 25 diverse soybean genotypes. Over the >76 kbp, mean nucleotide diversity expressed as Watterson's theta was 0.00097. Nucleotide diversity was 0.00053 and 0.00111 in coding and in noncoding perigenic DNA, respectively, lower than estimates in the autogamous model species Arabidopsis thaliana. Haplotype analysis of SNP-containing fragments revealed a deficiency of haplotypes vs. the number that would be anticipated at linkage equilibrium. In 49 fragments with three or more SNPs, five haplotypes were present in one fragment while four or less were present in the remaining 48, thereby supporting the suggestion of relatively limited genetic variation in cultivated soybean. Squared allele-frequency correlations (r(2)) among haplotypes at 54 loci with two or more SNPs indicated low genome-wide LD. The low level of LD and the limited haplotype diversity suggested that the genome of any given soybean accession is a mosaic of three or four haplotypes. To facilitate SNP discovery and the development of a transcript map, subsets of four to six diverse genotypes, whose sequence analysis would permit the discovery of at least 75% of all SNPs present in the 25 genotypes as well as 90% of the common (frequency >0.10) SNPs, were identified.

DNA sequence variation in BpMADS2 gene in two populations of Betula pendula

The PISTILLATA (PI) homologue, BpMADS2, was isolated from silver birch (Betula pendula Roth) and used to study nucleotide polymorphism. Two regions (together about 2450 bp) comprising mainly untranslated sequences were sequenced from 10 individuals from each of two populations in Finland. The nucleotide polymorphism was low in the BpMADS2 locus, especially in the coding region. The synonymous site overall nucleotide diversity (pis) was 0.0043 and the nonsynonymous nucleotide diversity (pia) was only 0.000052. For the whole region, the pi values for the two populations were 0.0039 and 0.0045, and for the coding regions, the pi values were only 0 and 0.00066 (for the corresponding coding regions of Arabidopsis thaliana PI world-wide pi was 0.0021). Estimates of pi or theta did not differ significantly between the two populations, and the two populations were not diverged from each other. Two classes of BpMADS2 alleles were present in both populations, suggesting that this gene exhibits allelic dimorphism. In addition to the nucleotide site variation, two microsatellites were also associated within the haplotypes. This allelic dimorphism might be the result of postglacial re-colonization partly from northwestern, partly from southeastern/eastern refugia. The sequence comparison detected five recombination events in the regions studied. The large number of microsatellites in all of the three introns studied suggests that BpMADS2 is a hotspot for microsatellite formation.

DNA polymorphism at the FRIGIDA gene in Arabidopsis thaliana: extensive nonsynonymous variation is consistent with local selection for flowering time

FRIGIDA (FRI) is a major gene involved in the regulation of flowering time in Arabidopsis thaliana. Nucleotide variation at this gene was investigated by sequencing 25 field ecotypes collected from western Europe. Genetic diversity at FRI was characterized by a high number of haplotypes and an excess of low-frequency polymorphisms. A large excess of intraspecific nonsynonymous variation associated with low synonymous variation was detected along the first exon in the FRI gene. In contrast, no excess of nonsynonymous divergence was detected between A. thaliana and A. lyrata. The Tajima and McDonald and Kreitman tests, however, suggested that this gene has evolved in a nonneutral fashion. Nonsynonymous variation included eight loss-of-function mutations that have probably arisen recently and independently in several locations. A phenotypic evaluation of the sequenced ecotypes confirmed that these loss-of-function mutations were associated with an early-flowering phenotype. Taken together, our results suggest that DNA polymorphism at the FRI gene in A. thaliana from western Europe has been shaped by recent positive selection for earliness in a set of isolated populations.

Recombination and gene conversion in a 170-kb genomic region of Arabidopsis thaliana

Arabidopsis thaliana is a highly selfing plant that nevertheless appears to undergo substantial recombination. To reconcile its selfing habit with the observations of recombination, we have sampled the genetic diversity of A. thaliana at 14 loci of approximately 500 bp each, spread across 170 kb of genomic sequence centered on a QTL for resistance to herbivory. A total of 170 of the 6321 nucleotides surveyed were polymorphic, with 169 being biallelic. The mean silent genetic diversity (pi(s)) varied between 0.001 and 0.03. Pairwise linkage disequilibria between the polymorphisms were negatively correlated with distance, although this effect vanished when only pairs of polymorphisms with four haplotypes were included in the analysis. The absence of a consistent negative correlation between distance and linkage disequilibrium indicated that gene conversion might have played an important role in distributing genetic diversity throughout the region. We tested this by coalescent simulations and estimate that up to 90% of recombination is due to gene conversion.

DNA variation in the acidic chitinase locus (ChiA) region in Arabis gemmifera and its related species

We analyzed DNA variation at the acidic chitinase (ChiA) locus of Arabis gemmifera and among its eight related species. Nucleotide diversity (pi) of the entire locus in A. gemmifera was 0.0032, which was one third that of A. thaliana. In A. gemmifera, an excess of unique polymorphisms yielded significantly negative results with the tests of Tajima and Fu and Li. The McDonald and Kreitman test revealed that the ratio of nucleotide replacement to synonymous changes in A. gemmifera was significantly greater than those between A. gemmifera and A. glabra, A. gemmifera and A. griffithiana, A. gemmifera and A. korshinskyi, A. gemmifera and A. wallichii, and A. gemmifera and A. himalaica. These results indicated that the neutrality assumption, the equilibrium population assumption, or both, could not be applied to the ChiA locus of A. gemmifera. The small size and relative isolation of local subpopulations of A. gemmifera could explain the excess of unique polymorphisms and the high proportion of replacement changes. The specific sampling scheme of this study, where one strain each was sampled from each local subpopulation might also have led to an excess of singletons. Interspecific comparison among Arabidopsis, Arabis and Cardaminopsis species showed that Ka was always lower than Ks, providing evidence against the adaptive evolution of ChiA. However, Ka/Ks was greater between closely related species than between more distant related species. ChiA had a higher level of replacement divergence and a lower level of synonymous divergence compared than did Adh. We suggest that both the mutation rate at the nucleotide level and the selective constraints at the protein level are lower in ChiA than in Adh.

Sequence variation and haplotype structure surrounding the flowering time locus FRI in Arabidopsis thaliana

Linkage disequilibrium in highly selfing organisms is expected to extend well beyond the scale of individual genes. The pattern of polymorphism in such species must thus be studied over a larger scale. We sequenced 14 short (0.5-1 kb) fragments from a 400-kb region surrounding the flowering time locus FRI in a sample of 20 accessions of Arabidopsis thaliana. The distribution of allele frequencies, as quantified by Tajima's D, varies considerably over the region and is incompatible with a standard neutral model. The region is characterized by extensive haplotype structure, with linkage disequilibrium decaying over 250 kb. In particular, recombination is evident within 35 kb of FRI in a haplotype associated with a functionally important allele. This suggests that A. thaliana may be highly suitable for linkage disequilibrium mapping.

Centromeric localization and adaptive evolution of an Arabidopsis histone H3 variant

Centromeric H3-like histones, which replace histone H3 in the centromeric chromatin of animals and fungi, have not been reported in plants. We identified a histone H3 variant from Arabidopsis thaliana that encodes a centromere-identifying protein designated HTR12. By immunological detection, HTR12 localized at centromeres in both mitotic and meiotic cells. HTR12 signal revealed tissue- and stage-specific differences in centromere morphology, including a distended bead-like structure in interphase root tip cells. The anti-HTR12 antibody also detected spherical organelles in meiotic cells. Although the antibody does not label centromeres in the closely related species Arabidopsis arenosa, HTR12 signal was found on all centromeres in allopolyploids of these two species. Comparison of the HTR12 genes of A. thaliana and A. arenosa revealed striking adaptive evolution in the N-terminal tail of the protein, similar to the pattern seen in its counterpart in Drosophila. This finding suggests that the same evolutionary forces shape centromeric chromatin in both animals and plants.

Molecular characterization of a second copy of holocarboxylase synthetase gene (hcs2) in Arabidopsis thaliana

Holocarboxylase synthetase (HCS), catalyzing the covalent attachment of biotin, is ubiquitously represented in living organisms. Indeed, the biotinylation is a post-translational modification that allows the transformation of inactive biotin-dependent carboxylases, which are committed in fundamental metabolisms such as fatty acid synthesis, into their active holo form. Among other living organisms, plants present a peculiarly complex situation. In pea, HCS activity has been detected in three subcellular compartments and the systematic sequencing of the Arabidopsis genome revealed the occurrence of two hcs genes (hcs1 and hcs2). Hcs1 gene product had been previously characterized at molecular and biochemical levels. Here, by PCR amplification, we cloned an hcs2 cDNA from Arabidopsis thaliana (Ws ecotype) mRNA. We observed the occurrence of multiple cDNA forms which resulted from the alternative splicing of hcs2 mRNA. Furthermore, we evidenced a nucleotide polymorphism at the hcs2 gene within the Ws ecotype, which affected splicing of hcs2 mRNA. This contrasted sharply with the situation at hcs1 locus. However, this polymorphism had no apparent effect on total HCS activity in planta. Finally, hcs2 mRNAs were found 4-fold less abundant than hcs1 mRNA and the most abundant hcs2 mRNA spliced variant should code for a truncated protein. We discuss the possible role of such a multiplicity of putative HCS proteins in plants and discuss the involvement of each of hcs genes in the correct realization of biotinylation.

Genetic diversity in Arabidopsis thaliana L. Heynh. investigated by cleaved amplified polymorphic sequence (CAPS) and inter-simple sequence repeat (ISSR) markers

In this study, we investigated genetic diversity among 37 accessions in Arabidopsis thaliana from Eurasia, North Africa and North America using morphological traits and two polymerase chain reaction (PCR)-based marker systems: cleaved amplified polymorphic sequences (CAPS) and inter-simple sequence repeats (ISSR). Cluster analysis based on genetic similarities calculated from CAPS data grouped the accessions roughly according to their geographical origin: one large group contained accessions from Western, Northern and Southern Europe as well as North Africa, a second group consisted of Eastern European and Asian continental accessions. North American accessions were interspersed into these groups. Contrary to the CAPS analysis, the dendrogram obtained from the ISSR data did not reflect the geographical origin of the accessions, and the calculated genetic distances did not match the CAPS results. This could be attributable to an uneven genomic distribution of ISSR markers as substantiated by a database search for ISSR binding sites in A. thaliana genomic DNA sequence files, or to the ISSR's different mode of evolution. We recommend CAPS markers for diversity analysis in A. thaliana because a careful selection of markers can ascertain an even representation of the entire genome.

The extent of linkage disequilibrium in Arabidopsis thaliana

Linkage disequilibrium (LD), the nonrandom occurrence of alleles in haplotypes, has long been of interest to population geneticists. Recently, the rapidly increasing availability of genomic polymorphism data has fueled interest in LD as a tool for fine-scale mapping, in particular for human disease loci. The chromosomal extent of LD is crucial in this context, because it determines how dense a map must be for associations to be detected and, conversely, limits how finely loci may be mapped. Arabidopsis thaliana is expected to harbor unusually extensive LD because of its high degree of selfing. Several polymorphism studies have found very strong LD within individual loci, but also evidence of some recombination. Here we investigate the pattern of LD on a genomic scale and show that in global samples, LD decays within approximately 1 cM, or 250 kb. We also show that LD in local populations may be much stronger than that of global populations, presumably as a result of founder events. The combination of a relatively high level of polymorphism and extensive haplotype structure bodes well for developing a genome-wide LD map in A. thaliana.

The oxylipin pathway in Arabidopsis

Oxylipins are acyclic or cyclic oxidation products derived from the catabolism of fatty acids which regulate many defense and developmental pathways in plants. The dramatic increase in the volume of publications and reviews on these compounds since 1997 documents the increasing interest in this compound and its role in plants. Research on this topic has solidified our understanding of the chemistry and biosynthetic pathways for oxylipin production. However, more information is still needed on how free fatty acids are produced and the role of beta-oxidation in the biosynthetic pathway for oxylipins. It is also becoming apparent that oxylipin content and composition changes during growth and development and during pathogen or insect attack. Oxylipins such as jasmonic acid (JA) or 12-oxo-phytodienoic acid modulate the expression of numerous genes and influence specific aspects of plant growth, development and responses to abiotic and biotic stresses. Although oxylipins are believed to act alone, several examples were presented to illustrate that JA-induced responses are modulated by the type and the nature of crosstalk with other signaling molecules such as ethylene and salicylic acid. How oxylipins cause changes in gene expression and instigate a physiological response is becoming understood with the isolation of mutations in both positive and negative regulators in the jasmonate signaling pathway and the use of cDNA microarrays.

Trichome distribution in Arabidopsis thaliana and its close relative Arabidopsis lyrata: molecular analysis of the candidate gene GLABROUS1

GLABROUS1 (GL1) belongs to the large family of MYB transcription factors and is known to play a central role in trichome initiation. We studied trichome distribution and the molecular variation of GL1 in 28 A. thaliana accessions. Trichome density on rosette leaves was highly variable among those accessions. On the molecular level, we detected substantial sequence variation in a 3-kb fragment which included the complete coding region of the GL1 locus (pi = 0.01). Phylogenetic analysis of GL1 indicates the presence of two diverged clades among 28 accessions. Using ANOVA, we show that the phenotypic variation in trichome density cannot be explained by the sequence divergence between the two phylogenetic lineages. Sequence analysis of wild-type Arabidopsis thaliana and Arabidopsis lyrata accessions indicates that all amino acid substitutions are located outside of the conserved helix-turn-helix DNA-binding domains R2 and R3. Using plants of A. thaliana and A. lyrata with either naturally occurring or ethyl methane sulfonate--induced glabrous phenotypes, we demonstrate that the last 14 C-terminal amino acids of the GL1 gene have no major impact on the initiation of trichomes.

DNA variation in the 5' upstream region of the Adh locus of the wild plants Arabidopsis thaliana and Arabis gemmifera

To investigate the level and pattern of DNA polymorphism in the noncoding regulatory region in the plant nuclear genome, 2.4 kb of nucleotide sequence of the 5' upstream region of ADH: was determined for 14 ecotypes of Arabidopsis thaliana and five accessions of Arabis gemmifera. Using this data set and previously determined ADH: sequence data, DNA variation was analyzed in a 4.4-kb region of the locus. Two divergent sequence types detected in the transcriptional unit of ADH: were not present in the 5' region of the ADH: gene in A. thaliana. Nucleotide diversity of the entire 5' region was estimated to be 0.0040, which is lower than that in the transcriptional unit. The level of variation was not uniform. There were peaks of variations in a approximately 400-bp region where cis-regulatory elements for ADH: expression were clustered and in exon 4. In interspecific comparison with A. gemmifera, lower divergence was observed in the 5' flanking region than in the exons. High peaks of divergence in the 400-bp regulatory region and exon 4 were also detected, although there were many other peaks. These results indicate that regions of functional importance have a high level of polymorphism and divergence in the ADH: locus of these genera. The possibility of balancing selection in the ADH: gene of these plants is discussed.

Heterogeneous geographic patterns of nucleotide sequence diversity between two alcohol dehydrogenase genes in wild barley (Hordeum vulgare subspecies spontaneum)

Patterns of nucleotide sequence diversity in the predominantly self-fertilizing species Hordeum vulgare subspecies spontaneum (wild barley) are compared between the putative alcohol dehydrogenase 3 locus (denoted "adh3") and alcohol dehydrogenase 1 (adh1), two related but unlinked loci. The data consist of a sequence sample of 1,873 bp of "adh3" drawn from 25 accessions that span the species range. There were 104 polymorphic sites in the sequenced region of "adh3." The data reveal a strong geographic pattern of diversity at "adh3" despite geographic uniformity at adh1. Moreover, levels of nucleotide sequence diversity differ by nearly an order of magnitude between the two loci. Genealogical analysis resolved two distinct clusters of "adh3" alleles (dimorphic sequence types) that coalesce roughly 3 million years ago. One type consists of accessions from the Middle East, and the other consists of accessions predominantly from the Near East. The two "adh3" sequence types are characterized by a high level of differentiation between clusters ( approximately 2.2%), which induces an overall excess of intermediate frequency variants in the pooled sample. Finally, there is evidence of intralocus recombination in the "adh3" data, despite the high level of self-fertilization characteristic of wild barley.

Heterogeneous geographic patterns of nucleotide sequence diversity between two alcohol dehydrogenase genes in wild barley (Hordeum vulgare subspecies spontaneum)

Patterns of nucleotide sequence diversity in the predominantly self-fertilizing species Hordeum vulgare subspecies spontaneum (wild barley) are compared between the putative alcohol dehydrogenase 3 locus (denoted "adh3") and alcohol dehydrogenase 1 (adh1), two related but unlinked loci. The data consist of a sequence sample of 1,873 bp of "adh3" drawn from 25 accessions that span the species range. There were 104 polymorphic sites in the sequenced region of "adh3." The data reveal a strong geographic pattern of diversity at "adh3" despite geographic uniformity at adh1. Moreover, levels of nucleotide sequence diversity differ by nearly an order of magnitude between the two loci. Genealogical analysis resolved two distinct clusters of "adh3" alleles (dimorphic sequence types) that coalesce roughly 3 million years ago. One type consists of accessions from the Middle East, and the other consists of accessions predominantly from the Near East. The two "adh3" sequence types are characterized by a high level of differentiation between clusters ( approximately 2.2%), which induces an overall excess of intermediate frequency variants in the pooled sample. Finally, there is evidence of intralocus recombination in the "adh3" data, despite the high level of self-fertilization characteristic of wild barley.

Nucleotide sequence variation at two genes of the phenylpropanoid pathway, the FAH1 and F3H genes, in Arabidopsis thaliana

The FAH1 and F3H genes encode ferulate-5-hydroxylase and flavanone-3-hydroxylase, which are enzymes in the pathways leading to the synthesis of sinapic acid esters and flavonoids, respectively. Nucleotide variation at these genes was surveyed by sequencing a sample of 20 worldwide Arabidopsis thaliana ecotypes and one Arabidopsis lyrata spp. petraea stock. In contrast with most previously studied genes, the percentage of singletons was rather low in both the FAH1 and the F3H gene regions. There was, therefore, no footprint of a recent species expansion in the pattern of nucleotide variation in these regions. In both FAH1 and F3H, nucleotide variation was structured into two major highly differentiated haplotypes. In both genes, there was a peak of silent polymorphism in the 5' part of the coding region without a parallel increase in silent divergence. In FAH1, the peak was centered at the beginning of the second exon. In F3H, nucleotide diversity was highest at the beginning of the gene. The observed pattern of variation in both FAH1 and F3H, although suggestive of balancing selection, was compatible with a neutral model with no recombination.

Genetic isolation by distance in Arabidopsis thaliana: biogeography and postglacial colonization of Europe

Arabidopsis thaliana provides a useful model system for functional, evolutionary and ecological studies in plant biology. We have analysed natural genetic variation in A. thaliana in order to infer its biogeographical and historical distribution across Eurasia. We analysed 79 amplified fragment length polymorphism (AFLP) markers in 142 accessions from the species' native range, and found highly significant genetic isolation by distance among A. thaliana accessions from Eurasia and southern Europe. These spatial patterns of genetic variation suggest that A. thaliana colonized central and northern Europe from Asia and from Mediterranean Pleistocene refugia, a trend which has been identified in other species. Statistically significant levels of multilocus linkage disequilibrium suggest intermediate levels of disequilibrium among subsets of loci, and analysis of genetic relationships among accessions reveal a star or bush-like dendrogram with low bootstrap support. Taken together, it appears that there has been sufficient historical recombination in the A. thaliana genome such that accessions do not conform to a tree-like, bifurcating pattern of evolution - there is no 'ecotype phylogeny.' Nonetheless, significant isolation by distance provides a framework upon which studies of natural variation in A. thaliana may be designed and interpreted.