A gain-of-function polymorphism controlling complex traits and fitness in nature.

Identification of the causal genes that control complex trait variation remains challenging, limiting our appreciation of the evolutionary processes that influence polymorphisms in nature. We cloned a quantitative trait locus that controls plant defensive chemistry, damage by insect herbivores, survival, and reproduction in the natural environments where this polymorphism evolved. These ecological effects are driven by duplications in the BCMA (branched-chain methionine allocation) loci controlling this variation and by two selectively favored amino acid changes in the glucosinolate-biosynthetic cytochrome P450 proteins that they encode. These changes cause a gain of novel enzyme function, modulated by allelic differences in catalytic rate and gene copy number. Ecological interactions in diverse environments likely contribute to the widespread polymorphism of this biochemical function.

Genome wide analyses reveal little evidence for adaptive evolution in many plant species.

The relative contribution of advantageous and neutral mutations to the evolutionary process is a central problem in evolutionary biology. Current estimates suggest that whereas Drosophila, mice, and bacteria have undergone extensive adaptive evolution, hominids show little or no evidence of adaptive evolution in protein-coding sequences. This may be a consequence of differences in effective population size. To study the matter further, we have investigated whether plants show evidence of adaptive evolution using an extension of the McDonald-Kreitman test that explicitly models slightly deleterious mutations by estimating the distribution of fitness effects of new mutations. We apply this method to data from nine pairs of species. Altogether more than 2,400 loci with an average length of approximately 280 nucleotides were analyzed. We observe very similar results in all species; we find little evidence of adaptive amino acid substitution in any comparison except sunflowers. This may be because many plant species have modest effective population sizes.

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

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

Comparative genetic mapping in Boechera stricta, a close relative of Arabidopsis.

The angiosperm family Brassicaceae contains both the research model Arabidopsis (Arabidopsis thaliana) and the agricultural genus Brassica. Comparative genomics in the Brassicaceae has largely focused on direct comparisons between Arabidopsis and the species of interest. However, the reduced genome size and chromosome number (n = 5) of Arabidopsis complicates comparisons. Arabidopsis shows extensive genome and chromosome reshuffling compared to its close relatives Arabidopsis lyrata and Capsella rubella, both with n = 8. To facilitate comparative genomics across the Brassicaceae we recently outlined a system of 24 conserved chromosomal blocks based on their positions in an ancestral karyotype of n = 8, rather than by their position in Arabidopsis. In this report we use this system as a tool to understand genome structure and evolution in Boechera stricta (n = 7). B. stricta is a diploid, sexual, and highly self-fertilizing species occurring in mostly montane regions of western North America. We have created an F(2) genetic map of B. stricta based on 192 individuals scored at 196 microsatellite and candidate gene loci. Single-nucleotide polymorphism genotyping of 94 of the loci was done simultaneously using an Illumina bead array. The total map length is 725.8 cM, with an average marker spacing of 3.9 cM. There are no gaps greater than 19.3 cM. The chromosomal reduction from n = 8 to n = 7 and other genomic changes in B. stricta likely involved a pericentric inversion, a chromosomal fusion, and two reciprocal translocations that are easily visualized using the genomic blocks. Our genetic map will facilitate the analysis of ecologically relevant quantitative variation in Boechera. Sequence data from this article can be found in the GenBank/EMBL data libraries under accession numbers DU 667459 to DU 708532.

Comparative genomics in the Brassicaceae: a family-wide perspective.

Comparative genomics of Arabidopsis relatives has great potential to improve our understanding of molecular function and evolutionary processes. Recent studies of phylogenetic relationships within Brassicaceae and the publication of a new tribal classification scheme provide an important framework for comparative genomics research. Comparative linkage mapping and chromosome painting in the close relatives of Arabidopsis have inferred an ancestral karyotype of these species. In addition, comparative mapping to Brassica has identified genomic blocks that have been maintained since the divergence of the Arabidopsis and Brassica lineages. Several analyses of conserved non-coding regions have identified putative cis-regulatory sequences, and have highlighted the need for comparative sequencing at greater evolutionary distances. The development of new model species with novel physiological and ecological traits allows analysis of phenotypes that are not available in A. thaliana. Looking towards the future, we suggest a prioritized research agenda for comparative genomics in the Brassicaceae.

Positive selection driving diversification in plant secondary metabolism.

In Arabidopsis thaliana and related plants, glucosinolates are a major component in the blend of secondary metabolites and contribute to resistance against herbivorous insects. Methylthioalkylmalate synthases (MAM) encoded at the MAM gene cluster control an early step in the biosynthesis of glucosinolates and, therefore, are central to the diversification of glucosinolate metabolism. We sequenced bacterial artificial chromosomes containing the MAM cluster from several Arabidopsis relatives, conducted enzyme assays with heterologously expressed MAM genes, and analyzed MAM nucleotide variation patterns. Our results show that gene duplication, neofunctionalization, and positive selection provide the mechanism for biochemical adaptation in plant defense. These processes occur repeatedly in the history of the MAM gene family, indicating their fundamental importance for the evolution of plant metabolic diversity both within and among species.

Partial shotgun sequencing of the Boechera stricta genome reveals extensive microsynteny and promoter conservation with Arabidopsis.

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value < or = 10(-30)) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5' to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5' noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.

Comparative genomics as a tool for gene discovery.

With the increasing availability of data from multiple eukaryotic genome sequencing projects, attention has focused on interspecific comparisons to discover novel genes and transcribed genomic sequences. Generally, these extrinsic strategies combine ab initio gene prediction with expression and/or homology data to identify conserved gene candidates between two or more genomes. Interspecific sequence analyses have proven invaluable for the improvement of existing annotations, automation of annotation, and identification of novel coding regions and splice variants. Further, comparative genomic approaches hold the promise of improved prediction of terminal or small exons, microRNA precursors, and small peptide-encoding open reading frames--sequence elements that are difficult to identify through purely intrinsic methodologies in the absence of experimental data.

Geographic and evolutionary diversification of glucosinolates among near relatives of Arabidopsis thaliana (Brassicaceae).

Glucosinolates are biologically active secondary metabolites that display both intra- and interspecific variation in the order Brassicales. Glucosinolate profiles have not been interpreted within a phylogenic framework and little is known regarding the processes that influence the evolution of glucosinolate diversity at a macroevolutionary scale. We have analyzed leaf glucosinolate profiles from members of the Brassicaceae that have diverged from Arabidopsis thaliana within the last 15 million years and interpreted our findings relative to the phylogeny of this group. We identified several interspecific polymorphisms in glucosinolate composition. A majority of these polymorphisms are lineage-specific secondary losses of glucosinolate characters, but a gain-of-character polymorphism was also detected. The genetic basis of most observed polymorphisms appears to be regulatory. In the case of A. lyrata, geographic distribution is also shown to contribute to glucosinolate metabolic diversity. Further, we observed evidence of gene-flow between sympatric species, parallel evolution, and the existence of genetic constraints on the evolution of glucosinolates within the Brassicaceae.