Genomic approaches to analyzing natural variation in Arabidopsis thaliana

The genomics tools available for studying Arabidopsis thaliana are a great resource for researchers trying to characterize and understand the genetic basis of natural variation. Abundant polymorphic markers aid quantitative trait locus (QTL) mapping, the fully sequenced genome provides rapid identification of candidate loci, and extensive knockout collections allow those candidate loci to be tested. Combining QTL mapping of classic phenotypic traits with biochemical or expression analysis is providing mechanistic insight into the traits of interest. Conversely, natural variation studies are now being done on genomic traits such as methylation or chiasma frequency.     

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.     

The glucosinolate-myrosinase system in an ecological and evolutionary context

Functional analysis of natural variation in the model species Arabidopsis thaliana has enabled the cloning of many glucosinolate biosynthesis and hydrolysis genes. Variation in these genes is central to understanding the ecological role of the glucosinolate-myrosinase defense system, and allows us to dissect the evolutionary and ecological forces that shape polymorphism at underlying loci. These same genes are also variable in other crucifer species, suggesting the presence of recurring selection, possibly mediated by insects. By utilizing the genomic tools available in A. thaliana to investigate these loci fully, it might be possible to generate detailed evolutionary or ecological models to apply to other species.     

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.     

Birth, death and subfunctionalization in the Arabidopsis genome

Arabidopsis thaliana is now a model system, not just for plant biology but also for comparative genomics. The completion of the sequences of two closely related species, Arabidopsis lyrata and Brassica rapa, is complemented by genomic comparisons among A. thaliana accessions and mutation accumulation lines. Together these genomic data document the birth of new genes via gene duplication, transposon exaptation and de novo formation of new genes from noncoding sequence. Most novel loci exhibit low expression, and are undergoing pseudogenization or subfunctionalization. Comparatively, A. thaliana has lost large amounts of sequence through deletion, particularly of transposable elements. Intraspecific genomic variation indicates high rates of deletion mutations and deletion polymorphisms across accessions, shedding light on the history of Arabidopsis genome architecture.     

Natural variation in MAM within and between populations of Arabidopsis lyrata determines glucosinolate phenotype

The genetic variation that underlies the glucosinolate phenotype of Arabidopsis lyrata ssp. petraea was investigated between and within populations. A candidate glucosinolate biosynthetic locus (MAM, containing methylthioalkylmalate synthase genes) was mapped in A. lyrata to a location on linkage group 6 corresponding to the homologous location for MAM in A. thaliana. In A. thaliana MAM is responsible for side chain elongation in aliphatic glucosinolates, and the MAM phenotype can be characterized by the ratios of long- to short-chain glucosinolates. A quantitative trait loci (QTL) analysis of glucosinolate ratios in an A. lyrata interpopulation cross found one QTL at MAM. Additional QTL were identified for total indolic glucosinolates and for the ratio of aliphatic to indolic glucosinolates. MAM was then used as the candidate gene for a within-population cosegregation analysis in a natural A. lyrata population from Germany. Extensive variation in microsatellite markers at MAM was found and this variation cosegregated with the same glucosinolate ratios as in the QTL study. The combined results indicate that both between- and within-population genetic variation in the MAM region determines phenotypic variation in glucosinolate side chains in A. lyrata.     

Comparative mapping in Arabidopsis and Brassica, fine scale genome collinearity and congruence of genes controlling flowering time

The model dicotyledonous plant, Arabidopsis thaliana , is closely related to Brassica crop species. It is intended that information concerning the genetic control of basic biological processes in Arabidopsis will be transferable to other species. Genome collinearity and its potential to facilitate the identification of candidate genes in Arabidopsis homologous to genes controlling important agronomic traits in Brassica was investigated. Genetic mapping in B. nigra identified two loci influencing flowering time (FT), with loci on linkage groups 2 and 8 explaining 53% and 12% of the total variation in FT, respectively. The CO gene exerts an important control over FT in A. thaliana , and B. nigra homologues of CO probably also play an important role in regulating FT. B. nigra homologues of CO were identified on linkage groups 2 and 8, the homologue on group 2 was coincident with the major locus controlling FT while the homologue on group 8 was within the 90% confidence interval of the weaker FT gene. The CO homologue on group 2 exhibits abundant allelic variation suggesting that it naturally controls a wide range of flowering times. Fine-scale A. thaliana/B. nigra comparative mapping demonstrated short-range collinearity between the genomes of Arabidopsis and Brassica . Eleven DNA fragments spaced over a 1.5 Mb contig in A. thaliana were used as RFLP probes in B. nigra . Three collinear representations of the A. thaliana contig were identified in B. nigra , with one interrupted by a large chromosomal inversion. Collinearity over this range will allow the resources generated by the Arabidopsis genome project to facilitate map-based cloning in Brassica crops.     

Novel loci control variation in reproductive timing in Arabidopsis thaliana in natural environments

Molecular biologists are rapidly characterizing the genetic basis of flowering in model species such as Arabidopsis thaliana. However, it is not clear how the developmental pathways identified in controlled environments contribute to variation in reproductive timing in natural ecological settings. Here we report the first study of quantitative trait loci (QTL) for date of bolting (the transition from vegetative to reproductive growth) in A. thaliana in natural seasonal field environments and compare the results with those obtained under typical growth-chamber conditions. Two QTL specific to long days in the chamber were expressed only in spring-germinating cohorts in the field, and two loci specific to short days in the chamber were expressed only in fall-germinating cohorts, suggesting differential involvement of the photoperiod pathway in different seasonal environments. However, several other photoperiod-specific QTL with large effects in controlled conditions were undetectable in natural environments, indicating that expression of allelic variation at these loci was overridden by environmental factors specific to the field. Moreover, a substantial number of QTL with major effects on bolting date in one or more field environments were undetectable under controlled environment conditions. These novel loci suggest the involvement of additional genes in the transition to flowering under ecologically relevant conditions.     

Genome-wide comparison of nucleotide-binding site-leucine-rich repeat-encoding genes in Arabidopsis

Plants, like animals, use several lines of defense against pathogen attack. Prominent among genes that confer disease resistance are those encoding nucleotide-binding site-leucine-rich repeat (NB-LRR) proteins. Likely due to selection pressures caused by pathogens, NB-LRR genes are the most variable gene family in plants, but there appear to be species-specific limits to the number of NB-LRR genes in a genome. Allelic diversity within an individual is also increased by obligatory outcrossing, which leads to genome-wide heterozygosity. In this study, we compared the NB-LRR gene complement of the selfer Arabidopsis thaliana and its outcrossing close relative Arabidopsis lyrata. We then complemented and contrasted the interspecific patterns with studies of NB-LRR diversity within A. thaliana. Three important insights are as follows: (1) that both species have similar numbers of NB-LRR genes; (2) that loci with single NB-LRR genes are less variable than tandem arrays; and (3) that presence-absence polymorphisms within A. thaliana are not strongly correlated with the presence or absence of orthologs in A. lyrata. Although A. thaliana individuals are mostly homozygous and thus potentially less likely to suffer from aberrant interaction of NB-LRR proteins with newly introduced alleles, the number of NB-LRR genes is similar to that in A. lyrata. In intraspecific and interspecific comparisons, NB-LRR genes are also more variable than receptor-like protein genes. Finally, in contrast to Drosophila, there is a clearly positive relationship between interspecific divergence and intraspecific polymorphisms.     

Molecular population genetics of the Arabidopsis CLAVATA2 region. The genomic scale of variation and selection in a selfing species

The Arabidopsis thaliana CLAVATA2 (CLV2) gene encodes a leucine-rich repeat protein that regulates the development of the shoot meristem. The levels and patterns of nucleotide variation were assessed for CLV2 and 10 flanking genes that together span a 40-kb region of chromosome I. A total of 296 out of 7959 sequenced nucleotide sites were polymorphic. The mean levels of sequence diversity of the contiguous genes in this region are approximately twofold higher than those of other typical Arabidopsis nuclear loci. There is, however, wide variation in the levels and patterns of sequence variation among the 11 linked genes in this region, and adjacent genes appear to be subject to contrasting evolutionary forces. CLV2 has the highest levels of nucleotide variation in this region, a significant excess of intermediate frequency polymorphisms, and significant levels of intragenic linkage disequilibrium. Most alleles at CLV2 are found in one of three haplotype groups of moderate (>15%) frequency. These features suggest that CLV2 may harbor a balanced polymorphism.     

Plant genome organisation and diversity: the year of the junk!

Having gained a thorough understanding of the structure and organization of model plant genomes, such as those of Arabidopsis thaliana and rice, we have now started to investigate the most interesting aspect of genome structure - its variations. Variation in DNA sequence is responsible for the genetic component of phenotypic variation (i.e. the component upon which both natural and artificial selection act). Recent studies have started to shed light on sequence variation outside of the genic regions, owing mainly to large insertion/deletion (indel) polymorphisms caused by the presence or absence of transposable elements of different classes. In addition to long terminal repeat retrotransposons, DNA transposons have been shown to be responsible for these polymorphisms. These comprise Helitrons, CACTA and Mu-like elements that are capable of acquiring and piecing together fragments of plant genes and are often expressed. Future analyses of the functional roles of intergenic sequence variation will tell us if we will need to pay more attention not only to genes, but also to the 'junk' DNA surrounding them.     

Role of FRIGIDA and FLOWERING LOCUS C in determining variation in flowering time of Arabidopsis

Arabidopsis (Arabidopsis thaliana) accessions provide an excellent resource to dissect the molecular basis of adaptation. We have selected 192 Arabidopsis accessions collected to represent worldwide and local variation and analyzed two adaptively important traits, flowering time and vernalization response. There was huge variation in the flowering habit of the different accessions, with no simple relationship to latitude of collection site and considerable diversity occurring within local regions. We explored the contribution to this variation from the two genes FRIGIDA (FRI) and FLOWERING LOCUS C (FLC), previously shown to be important determinants in natural variation of flowering time. A correlation of FLC expression with flowering time and vernalization was observed, but it was not as strong as anticipated due to many late-flowering/vernalization-requiring accessions being associated with low FLC expression and early-flowering accessions with high FLC expression. Sequence analysis of FRI revealed which accessions were likely to carry functional alleles, and, from comparison of flowering time with allelic type, we estimate that approximately 70% of flowering time variation can be accounted for by allelic variation of FRI. The maintenance and propagation of 20 independent nonfunctional FRI haplotypes suggest that the loss-of-function mutations can confer a strong selective advantage. Accessions with a common FRI haplotype were, in some cases, associated with very different FLC levels and wide variation in flowering time, suggesting additional variation at FLC itself or other genes regulating FLC. These data reveal how useful these Arabidopsis accessions will be in dissecting the complex molecular variation that has led to the adaptive phenotypic variation in flowering time.