The evolution of quantitative traits in complex environments.

Species inhabit complex environments and respond to selection imposed by numerous abiotic and biotic conditions that vary in both space and time. Environmental heterogeneity strongly influences trait evolution and patterns of adaptive population differentiation. For example, heterogeneity can favor local adaptation, or can promote the evolution of plastic genotypes that alter their phenotypes based on the conditions they encounter. Different abiotic and biotic agents of selection can act synergistically to either accelerate or constrain trait evolution. The environmental context has profound effects on quantitative genetic parameters. For instance, heritabilities measured in controlled conditions often exceed those measured in the field; thus, laboratory experiments could overestimate the potential for a population to respond to selection. Nevertheless, most studies of the genetic basis of ecologically relevant traits are conducted in simplified laboratory environments, which do not reflect the complexity of nature. Here, we advocate for manipulative field experiments in the native ranges of plant species that differ in mating system, life-history strategy and growth form. Field studies are vital to evaluate the roles of disparate agents of selection, to elucidate the targets of selection and to develop a nuanced perspective on the evolution of quantitative traits. Quantitative genetics field studies will also shed light on the potential for natural populations to adapt to novel climates in highly fragmented landscapes. Drawing from our experience with the ecological model system Boechera (Brassicaceae), we discuss advancements possible through dedicated field studies, highlight future research directions and examine the challenges associated with field studies.

The calmodulin-binding transcription factor SIGNAL RESPONSIVE1 is a novel regulator of glucosinolate metabolism and herbivory tolerance in Arabidopsis.

The Arabidopsis Ca(2+)/calmodulin (CaM)-binding transcription factor SIGNAL RESPONSIVE1 (AtSR1/CAMTA3) was previously identified as a key negative regulator of plant immune responses. Here, we report a new role for AtSR1 as a critical component of plant defense against insect herbivory. Loss of AtSR1 function impairs tolerance to feeding by the generalist herbivore Trichoplusia ni as well as wound-induced jasmonate accumulation. The susceptibility of the atsr1 mutant is associated with decreased total glucosinolate (GS) levels. The two key herbivory deterrents, indol-3-ylmethyl (I3M) and 4-methylsulfinylbutyl (4MSOB), showed the most significant reductions in atsr1 plants. Further, changes in AtSR1 transcript levels led to altered expression of several genes involved in GS metabolism including IQD1, MYB51 and AtST5a. Overall, our results establish AtSR1 as an important component of plant resistance to insect herbivory as well as one of only three described proteins involved in Ca(2+)/CaM-dependent signaling to function in the regulation of GS metabolism, providing a novel avenue for future investigations of plant-insect interactions.

Ecological genomics of Boechera stricta: identification of a QTL controlling the allocation of methionine- vs branched-chain amino acid-derived glucosinolates and levels of insect herbivory.

In the Brassicaceae, glucosinolates influence the feeding, reproduction and development of many insect herbivores. Glucosinolate production and effects on herbivore feeding have been extensively studied in the model species, Arabidopsis thaliana and Brassica crops, both of which constitutively produce leaf glucosinolates mostly derived from the amino acid, methionine. Much less is known about the regulation or role in defense of glucosinolates derived from other aliphatic amino acids, such as the branched-chain amino acids (BCAA), valine and isoleucine. We have identified a glucosinolate polymorphism in Boechera stricta controlling the allocation to BCAA- vs methionine-derived glucosinolates in both leaves and seeds. B. stricta is a perennial species that grows in mostly undisturbed habitats of western North America. We have measured glucosinolate profiles and concentrations in 192 F(2) lines that have earlier been used for genetic map construction. We also performed herbivory assays on six F(3) replicates per F(2) line using the generalist lepidopteran, Trichoplusia ni. Quantitative trait locus (QTL) analysis identified a single locus controlling both glucosinolate profile and levels of herbivory, the branched chain-methionine allocation or BCMA QTL. We have delimited this QTL to a small genomic region with a 1.0 LOD confidence interval just 1.9 cm wide, which, in A. thaliana, contains approximately 100 genes. We also found that methionine-derived glucosinolates provided significantly greater defense than the BCAA-derived glucosinolates against feeding by this generalist insect herbivore. The future positional cloning of this locus will allow for testing various adaptive explanations.

Cis-regulatory evolution of chalcone-synthase expression in the genus Arabidopsis.

The contribution of cis-regulation to adaptive evolutionary change is believed to be essential, yet little is known about the evolutionary rules that govern regulatory sequences. Here, we characterize the short-term evolutionary dynamics of a cis-regulatory region within and among two closely related species, A. lyrata and A. halleri, and compare our findings to A. thaliana. We focused on the cis-regulatory region of chalcone synthase (CHS), a key enzyme involved in the synthesis of plant secondary metabolites. We observed patterns of nucleotide diversity that differ among species but do not depart from neutral expectations. Using intra- and interspecific F1 progeny, we have evaluated functional cis-regulatory variation in response to light and herbivory, environmental cues, which are known to induce CHS expression. We find that substantial cis-regulatory variation segregates within and among populations as well as between species, some of which results from interspecific genetic introgression. We further demonstrate that, in A. thaliana, CHS cis-regulation in response to herbivory is greater than in A. lyrata or A. halleri. Our work indicates that the evolutionary dynamics of a cis-regulatory region is characterized by pervasive functional variation, achieved mostly by modification of response modules to one but not all environmental cues. Our study did not detect the footprint of selection on this variation.

Biogeographic distribution of polyploidy and B chromosomes in the apomictic Boechera holboellii complex.

The Boechera holboellii complex comprises B. holboellii and B. drummondii, both of which can reproduce through sex or apomixis. Sexuality is associated with diploid individuals, whereas apomictic individuals are diploid or triploid and may additionally have B chromosomes. Using flow cytometry and karyotype analysis, we have shown that B chromosomes (a) occur in both diploid and triploid apomictic B. holboellii, (b) may occur in triploid B. drummondii, and (c) are dispensable for the plant. Both diploid and triploid karyotypes are found in multiple chloroplast haplotypes of both species, suggesting that triploid forms have originated multiple times during the evolution of this complex. B chromosome carriers are found in geographically and genetically distinct populations, but it is unknown whether the extra chromosomes are shared by common descent (single origin) or have originated via introgressive hybridization and repeated transitions from diploidy to triploidy. Diploid plants containing the Bs reproduce apomictically, suggesting that the supernumerary elements are associated with apomixis. Finally, our analyses of pollen size and viability suggest that irregular chromosome segregation in some triploid lineages may lead to the generation of diploid individuals which carry the B chromosomes.

Is the aneuploid chromosome in an apomictic Boechera holboellii a genuine B chromosome?

The Boechera holboellii complex comprises B. holboellii and B. drummondii, both of which can reproduce through sex or apomixis. Sexuality is associated with diploidy, whereas apomictic individuals can either be diploid, aneuploid or triploid. Aneuploid individuals are found in geographically and genetically distinct populations and contain a single extra chromosome. It is unknown whether the supernumerary chromosomes are shared by common descent (single origin) or have originated via introgressive hybridizations associated with the repeated transition from diploidy to triploidy. Diploid plants containing the extra chromosome(s) reproduce apomictically, suggesting that the supernumerary elements are associated with apomixis. In this study we compared flow cytometry data, chromosome morphology, and DNA sequences of sexual diploid and apomictic aneuploids in order to establish whether the extra chromosome fits the classical concept of a B chromosome. Karyotype analyses revealed that the supernumerary chromosome in the metaphase complement is heterochromatic and often smaller than the A chromosomes, and differs in length between apomictic plants from different populations. DNA sequence analyses furthermore demonstrated elevated levels of non-synonymous substitutions in one of the alleles, likely that on the aneuploid chromosome. Although the extra chromosome in apomictic Boechera does not go through normal reductional meiosis, in which it may get eliminated or accumulated by a B-chromosome-specific process, its variable size and heterochromatic nature does meet the remaining criteria for a genuine B chromosome in other species. Its prevalence and conserved genetic composition nonetheless implies that this chromosome, if truly a B, may be atypical with respect to its influence on its carriers.

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.

Extensive chloroplast haplotype variation indicates Pleistocene hybridization and radiation of North American Arabis drummondii, A. x divaricarpa, and A. holboellii (Brassicaceae).

Arabis drummondii, A. holboellii and their hybrid A. x divaricarpa are widespread perennials of open habitats in North America. A phylogenetic analysis based on noncoding chloroplast DNA sequences (trnL intron and trnL/F intergenic spacer) resolved A. drummondii as a monophyletic taxon, but found A. holboellii to bear chloroplast haplotypes from highly diverged evolutionary lineages. This raised the question of a possible polyphyletic origin of A. holboellii. Arabis x divaricarpa was found to be of recent and polytopic origin, a result consistent with its presumed hybrid origin. One hundred and three chloroplast haplotypes were detected within 719 Arabis accessions investigated. The majority of chloroplast-types were estimated to have arisen prior to the Wisconsin glaciation. Phylogeographical analysis using nested clade analysis, suggested for A. holboellii (i). past fragmentation events, partitioning genetic variation in several instances between the Sierra Nevada, the Southern Rocky Mountains and the Colorado Plateau on the one hand and the Central to Northern Rockies of the United States and adjacent Cascades on the other; and for both parental species (ii). recolonization of major areas formerly covered by the Wisconsin glaciation by three haplotypes; and (iii). restricted gene flow indicating isolation by distance in areas south of the last glacial maximum. Arabis x divaricarpa is closely codistributed with its parental species and resampled their haplotypes. The highest genetic diversity was found in the Rocky Mountains from Idaho and Montana south to Utah and Colorado. This area was further hypothesized to have played a major role in the origin of both parental species and probably represented an important glacial refugium. However, evidence for glacial refugia was also found in arctic and boreal regions of Alaska and near the Great Lakes. In comparison to nuclear ribosomal internal transcribed spacer data, chloroplast DNA divergence was very high and evidently predated the origin of A. drummondii and possibly A. holboellii. Divergence of major chloroplast lineages dates back to the middle of the Pleistocene at least. Extensive hybridization is the most likely evolutionary factor working on A. holboellii to explain the revealed discrepancy in nuclear DNA and chloroplast DNA diversification.

Establishment of a high-efficiency SNP-based framework marker set for Arabidopsis.

The major goal of this project was the establishment of a tool for rapid mapping of new mutations and genotyping in Arabidopsis consisting of at least 100 evenly spaced framework markers. We assembled a single nucleotide polymorphism (SNP)-based marker set consisting of 112 polymorphic sites with average spacing of 1.15 Mbp derived from an SNP database that we recently developed. This information was used to set up efficient SNP detection reactions based on multiplexed primer extension assays. The 112 Columbia (Col-0)/C24 framework markers were used to assemble 18 multiplexed SNaPshot assays with which up to eight separate loci can be genotyped in a single-tube/single-capillary format. In addition, for 110 framework markers matrix-assisted laser desorption/ionization time of flight (MALDI-ToF) assays have been established for high throughput analyses. We demonstrated the usefulness and the robustness of both procedures of this tool by genotyping 48 BC3F1 individuals created between the accessions Col-0 and C24. Subsets of 10-62 of the established markers discriminate between various combinations of the accessions Col-0, C24, Landsberg erecta (Ler), Cape Verdi Islands (Cvi) and Niederzenz (Nd). Using a subset of 17 evenly distributed and established SNP markers that are also polymorphic between Ler and Col-0, we were able to rapidly map a mutant gene (tbr1) to an interval of 2.3 Mbp in an Ler (tbr1) x Col-0 cross.

Evolution of plant resistance at the molecular level: ecological context of species interactions.

Molecular data regarding the diversity of plant loci involved in resistance to herbivores or pathogens are becoming increasingly available. These genes demonstrate variable patterns of diversity, suggesting that they differ in their evolutionary history. In parallel, the study of natural variation for resistance, generally conducted at the phenotypic level, has shown that resistance does not evolve solely under selection pressures exerted by enemies. Metapopulation dynamics and other ecological characteristics of interacting species also appear to have a large impact on resistance evolution. Until now, studies of resistance at the molecular level have been conducted separately from ecological studies in extant populations. Future progress requires an evolutionary approach integrating both molecular and ecological aspects of resistance evolution. Such an approach will contribute greatly to our understanding of the evolution of molecular diversity at loci involved in biotic stress.

Genetics of drought adaptation in Arabidopsis thaliana: I. Pleiotropy contributes to genetic correlations among ecological traits.

We examined patterns of genetic variance and covariance in two traits (i) carbon stable isotope ratio delta13C (dehydration avoidance) and (ii) time to flowering (drought escape), both of which are putative adaptations to local water availability. Greenhouse screening of 39 genotypes of Arabidopsis thaliana native to habitats spanning a wide range of climatic conditions, revealed a highly significant positive genetic correlation between delta13C and flowering time. Studies in a range of C3 annuals have also reported large positive correlations, suggesting the presence of a genetically based trade-off between mechanisms of dehydration avoidance (delta13C) and drought escape (early flowering). We examined the contribution of pleiotropy by using a combination of mutant and near-isogenic lines to test for positive mutational covariance between delta13C and flowering time. Ecophysiological mutants generally showed variation in delta13C but not flowering time. However, flowering time mutants generally demonstrated pleiotropic effects consistent with natural variation. Mutations that caused later flowering also typically resulted in less negative delta13C and thus probably higher water use efficiency. We found strong evidence for pleiotropy using near-isogenic lines of Frigida and Flowering locus C, cloned loci known to be responsible for natural variation in flowering time. These data suggest the correlated evolution of delta13C and flowering time is explained in part by the fixation of pleiotropic alleles that alter both delta13C and time to flowering.

Population genetics of tandem trypsin inhibitor genes in Arabidopsis species with contrasting ecology and life history.

Duplicated genes are important in the evolution and ecology of plant-defences because herbivore and pathogen attack can be countered via functional diversification at two levels: among duplicated loci and within loci. We explore molecular sequence variation for two members of a defence-related gene family, Arabidopsis thaliana trypsin inhibitors (ATTI), in A. thaliana and a closely related species, A. lyrata subspp. petraea. A worldwide sample of the inbreeding annual A. thaliana had less genetic variation at two ATTI loci (piTOTAL

Cross-species microsatellite markers for elucidating population genetic structure in Arabidopsis and Arabis (Brassicaeae).

Species closely related to model organisms present the opportunity to efficiently apply molecular and functional tools developed by a large research community to taxa with different ecological and evolutionary histories. We complied 42 microsatellite loci that amplify under common conditions in four closely related Arabidopsis: A. thaliana; A. halleri; A. lyrata ssp. lyrata; and A. lyrata ssp. petraea, as well as in one more distantly related crucifer; Arabis drummondii. Variation at these loci is amenable to a diversity of applications including population genetics, phylogeographical analyses, mapping of inter and intraspecific crosses, and recombination mapping. Our analysis of microsatellite variation illustrates significant differences in population genetic parameters among three Arabidopsis species. A population of A. thaliana, an inbreeding annual plant associated with disturbed habitats, was highly monomorphic (P = 8% percent polymorphic loci) and only 0.2% heterozygous for 648 locus-by-individual combinations. A population of the self-incompatible perennial herb, A. halleri, was more genetically variable (P = 71%) and had an excess of heterozygosity that may reflect a recent population bottleneck associated with human-mediated founder events. A population of the self-incompatible perennial herb, A. lyrata ssp. petraea, was even more genetically variable (P = 86%) and appeared to be at mutation-drift equilibrium. Population structure estimated from neutrally evolving loci provides an empirical expectation against which hypotheses of adaptive evolution at functional loci can be tested.

The Arabidopsis epithiospecifier protein promotes the hydrolysis of glucosinolates to nitriles and influences Trichoplusia ni herbivory.

Glucosinolates are anionic thioglucosides that have become one of the most frequently studied groups of defensive metabolites in plants. When tissue damage occurs, the thioglucoside linkage is hydrolyzed by enzymes known as myrosinases, resulting in the formation of a variety of products that are active against herbivores and pathogens. In an effort to learn more about the molecular genetic and biochemical regulation of glucosinolate hydrolysis product formation, we analyzed leaf samples of 122 Arabidopsis ecotypes. A distinct polymorphism was observed with all ecotypes producing primarily isothiocyanates or primarily nitriles. The ecotypes Columbia (Col) and Landsberg erecta (Ler) differed in their hydrolysis products; therefore, the Col x Ler recombinant inbred lines were used for mapping the genes controlling this polymorphism. The major quantitative trait locus (QTL) affecting nitrile versus isothiocyanate formation was found very close to a gene encoding a homolog of a Brassica napus epithiospecifier protein (ESP), which causes the formation of epithionitriles instead of isothiocyanates during glucosinolate hydrolysis in the seeds of certain Brassicaceae. The heterologously expressed Arabidopsis ESP was able to convert glucosinolates both to epithionitriles and to simple nitriles in the presence of myrosinase, and thus it was more versatile than previously described ESPs. The role of ESP in plant defense is uncertain, because the generalist herbivore Trichoplusia ni (the cabbage looper) was found to feed more readily on nitrile-producing than on isothiocyanate-producing Arabidopsis. However, isothiocyanates are frequently used as recognition cues by specialist herbivores, and so the formation of nitriles instead of isothiocyanates may allow Arabidopsis to be less apparent to specialists.

A gene controlling variation in Arabidopsis glucosinolate composition is part of the methionine chain elongation pathway.

Arabidopsis and other Brassicaceae produce an enormous diversity of aliphatic glucosinolates, a group of methionine (Met)-derived plant secondary compounds containing a beta-thio-glucose moiety, a sulfonated oxime, and a variable side chain. We fine-scale mapped GSL-ELONG, a locus controlling variation in the side-chain length of aliphatic glucosinolates. Within this locus, a polymorphic gene was identified that determines whether Met is extended predominantly by either one or by two methylene groups to produce aliphatic glucosinolates with either three- or four-carbon side chains. Two allelic mutants deficient in four-carbon side-chain glucosinolates were shown to contain independent missense mutations within this gene. In cell-free enzyme assays, a heterologously expressed cDNA from this locus was capable of condensing 2-oxo-4-methylthiobutanoic acid with acetyl-coenzyme A, the initial reaction in Met chain elongation. The gene methylthioalkylmalate synthase1 (MAM1) is a member of a gene family sharing approximately 60% amino acid sequence similarity with 2-isopropylmalate synthase, an enzyme of leucine biosynthesis that condenses 2-oxo-3-methylbutanoate with acetyl-coenzyme A.

Recurrent polyploid origins and chloroplast phylogeography in the Arabis holboellii complex (Brassicaceae).

Arabis holboellii is a North American member of the Brassicaceae that can reproduce via sex or apomixis. Previous studies have shown sexual individuals to be diploid, whilst apomictic individuals can be diploid (and aneuploid) or polyploid (typically 3x). Apomictic individuals can furthermore be facultative (i.e. both sexual and apomictic seed production in a single individual). Using flow cytometry, ploidy variation in 245 accessions of A. holboellii and A. drummondii from western North America and Greenland has been examined. Additionally, the chloroplast trnL intron region from each accession was sequenced for phylogenetic analysis of ploidy variation. Based upon 17 informative single nucleotide and insertion-deletion polymorphisms, we identified seven and 14 chloroplast haplotypes for A. drummondii and A. holboellii, respectively. Six of the haplotypes were found in both species. Ten of the chloroplast haplotypes were characterized by diploid, aneuploid, and triploid individuals, and thus we conclude that polyploidy has repeatedly and independently arisen within the species complex. As triploid individuals, which undergo normal meiosis, can only reproduce through apomixis, this may imply that the phenotype apomixis has also arisen multiple times. Arabis holboellii thus appears to have some predisposition to evolve apomictic reproduction.

SGP-1: prediction and validation of homologous genes based on sequence alignments.

Conventional methods of gene prediction rely on the recognition of DNA-sequence signals, the coding potential or the comparison of a genomic sequence with a cDNA, EST, or protein database. Reasons for limited accuracy in many circumstances are species-specific training and the incompleteness of reference databases. Lately, comparative genome analysis has attracted increasing attention. Several analysis tools that are based on human/mouse comparisons are already available. Here, we present a program for the prediction of protein-coding genes, termed SGP-1 (Syntenic Gene Prediction), which is based on the similarity of homologous genomic sequences. In contrast to most existing tools, the accuracy of depends little on species-specific properties such as codon usage or the nucleotide distribution. may therefore be applied to nonstandard model organisms in vertebrates as well as in plants, without the need for extensive parameter training. In addition to predicting genes in large-scale genomic sequences, the program may be useful to validate gene structure annotations from databases. To this end, SGP-1 output also contains comparisons between predicted and annotated gene structures in HTML format. The program can be accessed via a Web server at http://soft.ice.mpg.de/sgp-1. The source code, written in ANSI C, is available on request from the authors.

Comparative quantitative trait loci mapping of aliphatic, indolic and benzylic glucosinolate production in Arabidopsis thaliana leaves and seeds.

Secondary metabolites are a diverse set of plant compounds believed to have numerous functions in plant-environment interactions. Despite this importance, little is known about the regulation of secondary metabolite accumulation. We are studying the regulation of glucosinolates, a large group of secondary metabolites, in Arabidopsis to investigate how secondary metabolism is controlled. We utilized Ler and Cvi, two ecotypes of Arabidopsis that have striking differences in both the types and amounts of glucosinolates that accumulate in the seeds and leaves. QTL analysis identified six loci determining total aliphatic glucosinolate accumulation, six loci controlling total indolic glucosinolate concentration, and three loci regulating benzylic glucosinolate levels. Our results show that two of the loci controlling total aliphatic glucosinolates map to biosynthetic loci that interact epistatically to regulate aliphatic glucosinolate accumulation. In addition to the six loci regulating total indolic glucosinolate concentration, mapping of QTL for the individual indolic glucosinolates identified five additional loci that were specific to subsets of the indolic glucosinolates. These data show that there are a large number of variable loci controlling glucosinolate accumulation in Arabidopsis thaliana.

Comparative genomics and regulatory evolution: conservation and function of the Chs and Apetala3 promoters.

DNA sequence variations of chalcone synthase (Chs) and Apetala3 gene promoters from 22 cruciferous plant species were analyzed to identify putative conserved regulatory elements. Our comparative approach confirmed the existence of numerous conserved sequences which may act as regulatory elements in both investigated promoters. To confirm the correct identification of a well-conserved UV-light-responsive promoter region, a subset of Chs promoter fragments were tested in Arabidopsis thaliana protoplasts. All promoters displayed similar light responsivenesses, indicating the general functional relevance of the conserved regulatory element. In addition to known regulatory elements, other highly conserved regions were detected which are likely to be of functional importance. Phylogenetic trees based on DNA sequences from both promoters (gene trees) were compared with the hypothesized phylogenetic relationships (species trees) of these taxa. The data derived from both promoter sequences were congruent with the phylogenies obtained from coding regions of other nuclear genes and from chloroplast DNA sequences. This indicates that promoter sequence evolution generally is reflective of species phylogeny. Our study also demonstrates the great value of comparative genomics and phylogenetics as a basis for functional analysis of promoter action and gene regulation.