The Evolution of Sex is Tempered by Costly Hybridization in Boechera (Rock Cress)

Despite decades of research, the evolution of sex remains an enigma in evolutionary biology. Typically, research addresses the costs of sex and asexuality to characterize the circumstances favoring one reproductive mode. Surprisingly few studies address the influence of common traits that are, in many organisms, obligately correlated with asexuality, including hybridization and polyploidy. These characteristics have substantial impacts on traits under selection. In particular, the fitness consequences of hybridization (i.e., reduced fitness due to interspecific reproductive isolation) will influence the evolution of sex. This may comprise a cost of either sex or asexuality due to the link between hybridity and asexuality. We examined reproductive isolation in the formation of de novo hybrid lineages between 2 widespread species in the ecological model system Boechera. Seventeen percent of 664 crosses produced F1 fruits, and only 10% of these were viable, suggesting that postmating prezygotic and postzygotic barriers inhibit hybrid success in this system. The postmating prezygotic barrier was asymmetrical, with 110 of 115 total F1 fruits produced when Boechera stricta acted as maternal parent. This asymmetry was confirmed in wild-collected lineages, using a chloroplast phylogeny of wild-collected B. stricta, Boechera retrofracta, and hybrids. We next compared fitness of F2 hybrids and selfed parental B. stricta lines, finding that F2 fitness was reduced by substantial hybrid sterility. Multiple reproductively isolating barriers influence the formation and fitness of hybrid lineages in the wild, and the costs of hybridization likely have profound impacts on the evolution of sex in the natural environment.

Ecological differentiation facilitates fine-scale coexistence of sexual and asexual Boechera

PREMISE OF THE STUDY

Ecological differentiation (ED) between sexual and asexual organisms may permit the maintenance of reproductive polymorphism. Several studies of sexual/asexual ED in plants have shown that the geographic ranges of asexuals extend beyond those of sexuals, often in areas of higher latitude or elevation. But very little is known about ED at fine scales, wherein coexistence of sexuals and asexuals may be permitted by differential niche occupation.

METHODS

We used 149 populations of sexual and apomictic lineages in the genus Boechera (rock cress) collected across a portion of this mustard's vast range. We characterized reproductive mode, ploidy, and species identity or hybrid parentage of each individual, and then used a multipronged statistical approach to (1) identify ED between sexuals and asexuals; (2) investigate the impacts of two confounding factors, polyploidy and hybridization, on ED; and (3) determine the environmental variables underlying ED.

KEY RESULTS

We found that sexuals and asexuals are significantly ecologically differentiated across the landscape, despite fine-scale interdigitation of these two reproductive forms. Asexual reproduction was strongly associated with greater disturbance, reduced slope, and greater environmental variability. Although ploidy had little effect on the patterns observed, hybridization has a unique impact on the relationships between asexual reproduction and specific environmental variables.

CONCLUSIONS

Ecological differentiation along the axes of disturbance, slope, and climatic variability, as well as the effects of heterozygosity, may contribute to the maintenance of sexuality and asexuality across the landscape, ultimately impacting the establishment and spread of asexual lineages.

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.

Phenotypic plasticity and adaptive evolution contribute to advancing flowering phenology in response to climate change

Anthropogenic climate change has already altered the timing of major life-history transitions, such as the initiation of reproduction. Both phenotypic plasticity and adaptive evolution can underlie rapid phenological shifts in response to climate change, but their relative contributions are poorly understood. Here, we combine a continuous 38 year field survey with quantitative genetic field experiments to assess adaptation in the context of climate change. We focused on Boechera stricta (Brassicaeae), a mustard native to the US Rocky Mountains. Flowering phenology advanced significantly from 1973 to 2011, and was strongly associated with warmer temperatures and earlier snowmelt dates. Strong directional selection favoured earlier flowering in contemporary environments (2010-2011). Climate change could drive this directional selection, and promote even earlier flowering as temperatures continue to increase. Our quantitative genetic analyses predict a response to selection of 0.2 to 0.5 days acceleration in flowering per generation, which could account for more than 20 per cent of the phenological change observed in the long-term dataset. However, the strength of directional selection and the predicted evolutionary response are likely much greater now than even 30 years ago because of rapidly changing climatic conditions. We predict that adaptation will likely be necessary for long-term in situ persistence in the context of climate change.

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 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. Boechera 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₂ lines that have previously been used for genetic map construction. We also performed herbivory assays on six F₃ replicates per F₂ 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 ∼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.

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.

Comparing the linkage maps of the close relatives Arabidopsis lyrata and A. thaliana

We have constructed a genetic map of Arabidopsis lyrata, a self-incompatible relative of the plant model species A. thaliana. A. lyrata is a diploid (n = 8) species that diverged from A. thaliana (n = 5) approximately 5 MYA. Mapping was conducted in a full-sib progeny of two unrelated F(1) hybrids between two European populations of A. lyrata ssp. petraea. We used the least-squares method of the Joinmap program for map construction. The gross chromosomal differences between the two species were most parsimoniously explained with three fusions, two reciprocal translocations, and one inversion. The total map length was 515 cM, and the distances were 12% larger than those between corresponding markers in the linkage map of A. thaliana. The 72 markers, consisting of microsatellites and gene-based markers, were spaced on average every 8 cM. Transmission ratio distortion was extensive, and most distortions were specific to each reciprocal cross, suggesting cytoplasmic interactions. We estimate locations and most probable genotype frequencies of transmission ratio distorting loci (TRDL) with a Bayesian method and discuss the possible reasons for the observed distortions.

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.

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.

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.

Large-scale identification and analysis of genome-wide single-nucleotide polymorphisms for mapping in Arabidopsis thaliana

Genetic markers such as single nucleotide polymorphisms (SNPs) are essential tools for positional cloning, association, or quantitative trait locus mapping and the determination of genetic relationships between individuals. We identified and characterized a genome-wide set of SNP markers by generating 10,706 expressed sequence tags (ESTs) from cDNA libraries derived from 6 different accessions, and by analysis of 606 sequence tagged sites (STS) from up to 12 accessions of the model flowering plant Arabidopsis thaliana. The cDNA libraries for EST sequencing were made from individuals that were stressed by various means to enrich for transcripts from genes expressed under such conditions. SNPs discovered in these sequences may be useful markers for mapping genes involved in interactions with the biotic and abiotic environment. The STS loci are distributed randomly over the genome. By comparison with the Col-0 genome sequence, we identified a total of 8051 SNPs and 637 insertion/deletion polymorphisms (InDel). Analysis of STS-derived SNPs shows that most SNPs are rare, but that it is possible to identify intermediate frequency framework markers that can be used for genetic mapping in many different combinations of accessions. A substantial proportion of SNPs located in ORFs caused a change of the encoded amino acid. A comparison of the density of our SNP markers among accessions in both the EST and STS datasets, revealed that Cvi-0 is the most divergent accession from Col-0 among the 12 accessions studied. All of these markers are freely available via the internet.

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 <or= 0.0006) than observed previously at other functional loci. A significant excess of high frequency derived alleles in the signal sequence and 5' UTR of ATTI2 was consistent with a model of positive selection. However, demographic processes such as population subdivision and expansion, both likely to have occurred in A. thaliana during the last 10 000 years, can also give rise to similar deviations from neutrality. A single population of A. lyrata subspp. petraea in Germany had up to an order of magnitude more standing genetic variation at ATTI loci than the species-wide sample of A. thaliana. Although the level of variability for ATTI1 and ATTI2 within this single population was similar to, or even greater than, observed species-wide diversity for other loci in A. lyrata, there was little evidence to reject an equilibrium neutral model. A spatially explicit sample of 87 A. lyrata subspp. petraea individuals detected outbreeding (FIS = -0.16; FIT = -0.15) but little population subdivision (FST = 0.006) in this self-incompatible perennial herb. Genetic differences between Arabidopsis species were consistent with, but not fully explained by, divergence in ecology and life history. Diversification appears to have occurred in different functional domains for the tandemly duplicated ATTI1 and ATTI2 genes; the majority of fixed replacements in ATTI1 surround the enzyme binding site of the mature protein, whereas in ATTI2 most functional evolutionary change is located in the signal peptide. This pattern is consistent with a hypothesis of subfunctionalization in trypsin inhibitory function.

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.

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.

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.

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.

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.

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.

Local adaptation across a climatic gradient despite small effective population size in the rare sapphire rockcress

When assigning conservation priorities in endangered species, two common management strategies seek to protect remnant populations that (i) are the most genetically divergent or (ii) possess the highest diversity at neutral genetic markers. These two approaches assume that variation in molecular markers reflects variation in ecologically important traits and ignore the possibility of local adaptation among populations that show little divergence or variation at marker loci. Using common garden experiments, we demonstrate that populations of the rare endemic plant Arabis fecunda are physiologically adapted to the local microclimate. Local adaptation occurs despite (i) the absence of divergence at almost all marker loci and (ii) very small effective population sizes, as evidenced by extremely low levels of allozyme and DNA sequence polymorphism. Our results provide empirical evidence that setting conservation priorities based exclusively on molecular marker diversity may lead to the loss of locally adapted populations.

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.

Genetic control of natural variation in Arabidopsis glucosinolate accumulation

Glucosinolates are biologically active secondary metabolites of the Brassicaceae and related plant families that influence plant/insect interactions. Specific glucosinolates can act as feeding deterrents or stimulants, depending upon the insect species. Hence, natural selection might favor the presence of diverse glucosinolate profiles within a given species. We determined quantitative and qualitative variation in glucosinolates in the leaves and seeds of 39 Arabidopsis ecotypes. We identified 34 different glucosinolates, of which the majority are chain-elongated compounds derived from methionine. Polymorphism at only five loci was sufficient to generate 14 qualitatitvely different leaf glucosinolate profiles. Thus, there appears to be a modular genetic system regulating glucosinolate profiles in Arabidopsis. This system allows the rapid generation of new glucosinolate combinations in response to changing herbivory or other selective pressures. In addition to the qualitative variation in glucosinolate profiles, we found a nearly 20-fold difference in the quantity of total aliphatic glucosinolates and were able to identify a single locus that controls nearly three-quarters of this variation.

Molecular systematics of the Brassicaceae: evidence from coding plastidic matK and nuclear Chs sequences

Phylogenetic relationships were inferred using nucleotide sequence variation of the nuclear-encoded chalcone synthase gene (Chs) and the chloroplast gene matK for members of five tribes from the family Brassicaceae to analyze tribal and subtribal structures. Phylogenetic trees from individual data sets are mostly in congruence with the results from a combined matK-Chs analysis with a total of 2721 base pairs, but with greater resolution and higher statistical support for deeper branching patterns. The analysis indicates that tribes Lepidieae, Arabideae, and Sisymbrieae are not monophyletic. Among taxa under study four different lineages each were detected in tribes Arabideae and Lepidieae, interspersed with taxa from tribes Sisymbrieae, Hesperideae, and Brassiceae. It is concluded that tribe Brassiceae might be the only monophyletic group of the traditional tribes. From our data we estimated several divergence times for different lineages among cruciferous plants: 5.8 mya (million years ago) for the Arabidopsis-Cardaminopsis split, 20 mya for the Brassica-Arabidopsis split, and ∼40 mya for the age of the deepest split between the most basal crucifer Aethionema and remaining cruciferous taxa.

Gene duplication in the diversification of secondary metabolism: tandem 2-oxoglutarate-dependent dioxygenases control glucosinolate biosynthesis in Arabidopsis

Secondary metabolites are a diverse set of plant compounds believed to have numerous functions in plant-environment interactions. The large chemical diversity of secondary metabolites undoubtedly arises from an equally diverse set of enzymes responsible for their biosynthesis. However, little is known about the evolution of enzymes involved in secondary metabolism. We are studying the biosynthesis of glucosinolates, a large group of secondary metabolites, in Arabidopsis to investigate the evolution of enzymes involved in secondary metabolism. Arabidopsis contains natural variations in the presence of methylsulfinylalkyl, alkenyl, and hydroxyalkyl glucosinolates. In this article, we report the identification of genes encoding two 2-oxoglutarate--dependent dioxygenases that are responsible for this variation. These genes, AOP2 and AOP3, which map to the same position on chromosome IV, result from an apparent gene duplication and control the conversion of methylsulfinylalkyl glucosinolate to either the alkenyl or the hydroxyalkyl form. By heterologous expression in Escherichia and the correlation of gene expression patterns to the glucosinolate phenotype, we show that AOP2 catalyzes the conversion of methylsulfinylalkyl glucosinolates to alkenyl glucosinolates. Conversely, AOP3 directs the formation of hydroxyalkyl glucosinolates from methylsulfinylalkyl glucosinolates. No ecotype coexpressed both genes. Furthermore, the absence of functional AOP2 and AOP3 leads to the accumulation of the precursor methylsulfinylalkyl glucosinolates. A third member of this gene family, AOP1, is present in at least two forms and found in all ecotypes examined. However, its catalytic role is still uncertain.

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.

Induced plant defense responses against chewing insects. Ethylene signaling reduces resistance of Arabidopsis against Egyptian cotton worm but not diamondback moth

The induction of plant defenses by insect feeding is regulated via multiple signaling cascades. One of them, ethylene signaling, increases susceptibility of Arabidopsis to the generalist herbivore Egyptian cotton worm (Spodoptera littoralis; Lepidoptera: Noctuidae). The hookless1 mutation, which affects a downstream component of ethylene signaling, conferred resistance to Egyptian cotton worm as compared with wild-type plants. Likewise, ein2, a mutant in a central component of the ethylene signaling pathway, caused enhanced resistance to Egyptian cotton worm that was similar in magnitude to hookless1. Moreover, pretreatment of plants with ethephon (2-chloroethanephosphonic acid), a chemical that releases ethylene, elevated plant susceptibility to Egyptian cotton worm. By contrast, these mutations in the ethylene-signaling pathway had no detectable effects on diamondback moth (Plutella xylostella) feeding. It is surprising that this is not due to nonactivation of defense signaling, because diamondback moth does induce genes that relate to wound-response pathways. Of these wound-related genes, jasmonic acid regulates a novel beta-glucosidase 1 (BGL1), whereas ethylene controls a putative calcium-binding elongation factor hand protein. These results suggest that a specialist insect herbivore triggers general wound-response pathways in Arabidopsis but, unlike a generalist herbivore, does not react to ethylene-mediated physiological changes.

Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae)

We analyzed sequence variation for chalcone synthase (Chs) and alcohol dehydrogenase (Adh) loci in 28 species in the genera Arabidopsis and Arabis and related taxa from tribe Arabideae. Chs was single-copy in nearly all taxa examined, while Adh duplications were found in several species. Phylogenies constructed from both loci confirmed that the closest relatives of Arabidopsis thaliana include Arabidopsis lyrata, Arabidopsis petraea, and Arabidopsis halleri (formerly in the genus Cardaminopsis). Slightly more distant are the North American n = 7 Arabis (Boechera) species. The genus Arabis is polyphyletic-some unrelated species appear within this taxonomic classification, which has little phylogenetic meaning. Fossil pollen data were used to compute a synonymous substitution rate of 1.5 x 10 substitutions per site per year for both Chs and Adh. Arabidopsis thaliana diverged from its nearest relatives about 5 MYA, and from Brassica roughly 24 MYA. Independent molecular and fossil data from several sources all provide similar estimates of evolutionary timescale in the Brassicaceae.

Rapid evolution in plant chitinases: molecular targets of selection in plant-pathogen coevolution

Many pathogen recognition genes, such as plant R-genes, undergo rapid adaptive evolution, providing evidence that these genes play a critical role in plant-pathogen coevolution. Surprisingly, whether rapid adaptive evolution also occurs in genes encoding other kinds of plant defense proteins is unknown. Unlike recognition proteins, plant chitinases attack pathogens directly, conferring disease resistance by degrading chitin, a component of fungal cell walls. Here, we show that nonsynonymous substitution rates in plant class I chitinase often exceed synonymous rates in the plant genus Arabis (Cruciferae) and in other dicots, indicating a succession of adaptively driven amino acid replacements. We identify individual residues that are likely subject to positive selection by using codon substitution models and determine the location of these residues on the three-dimensional structure of class I chitinase. In contrast to primate lysozymes and plant class III chitinases, structural and functional relatives of class I chitinase, the adaptive replacements of class I chitinase occur disproportionately in the active site cleft. This highly unusual pattern of replacements suggests that fungi directly defend against chitinolytic activity through enzymatic inhibition or other forms of chemical resistance and identifies target residues for manipulating chitinolytic activity. These data also provide empirical evidence that plant defense proteins not involved in pathogen recognition also evolve in a manner consistent with rapid coevolutionary interactions.

Genetics and evolution of insect resistance in Arabidopsis

The genetic and molecular tools available in Arabidopsis allow identification of insect resistance genes. Many functional aspects of pest recognition and signal transduction are conserved in the defensive physiology of a broad range of plant species. Therefore, studies of insect resistance in Arabidopsis may be extended to functional genomics studies in many plant species of agricultural and ecological importance. Because of public concerns for field release of genetically modified organisms, naturally occurring genetic variation for resistance to insect herbivores will be valuable in plant breeding. Combined studies employing QTL mapping and candidate resistance genes are necessary to find and understand the genes responsible for variation in resistance. We review experiments showing that plant populations contain high levels of genetic variation for defensive physiology and disease and insect resistance, and that this variation can be manipulated to alter resistance and its components in a predictable fashion. In Arabidopsis, we can map the genes controlling physiological variation, and estimate the importance of regulatory or enzyme-encoding loci. Finally, we review functional genomics approaches for identification of insect resistance genes in Arabidopsis.

Plant-insect interactions

Recent research shows partially overlapping signal transduction pathways controlling responses to wounding, insects, and pathogens. Chemical and behavioral assays show that plants release herbivore-specific volatiles, and that parasitic wasps can distinguish between these emission patterns. QTL mapping and candidate gene studies are beginning to identify polymorphic resistance genes, and ecological analyses provide information on the physiological and fitness costs of resistance. Such multidisciplinary approaches can elucidate the physiological causes and ecological consequences of plant-herbivore interactions.

The molecular basis of quantitative genetic variation in central and secondary metabolism in Arabidopsis

To find the genes controlling quantitative variation, we need model systems where functional information on physiology, development, and gene regulation can guide evolutionary inferences. We mapped quantitative trait loci (QTLs) influencing quantitative levels of enzyme activity in primary and secondary metabolism in Arabidopsis. All 10 enzymes showed highly significant quantitative genetic variation. Strong positive genetic correlations were found among activity levels of 5 glycolytic enzymes, PGI, PGM, GPD, FBP, and G6P, suggesting that enzymes with closely related metabolic functions are coregulated. Significant QTLs were found influencing activity of most enzymes. Some enzyme activity QTLs mapped very close to known enzyme-encoding loci (e.g., hexokinase, PGI, and PGM). A hexokinase QTL is attributable to cis-acting regulatory variation at the AtHXK1 locus or a closely linked regulatory locus, rather than polypeptide sequence differences. We also found a QTL on chromosome IV that may be a joint regulator of GPD, PGI, and G6P activity. In addition, a QTL affecting PGM activity maps within 700 kb of the PGM-encoding locus. This QTL is predicted to alter starch biosynthesis by 3.4%, corresponding with theoretical models, suggesting that QTLs reflect pleiotropic effects of mutant alleles.

Genetics of Brassica rapa (syn. campestris). 2. Multiple disease resistance to three fungal pathogens: Peronospora parasitica, Albugo candida and Leptosphaeria maculans

Although the genetic basis of multiple disease resistance (MDR) is poorly understood, it is of great value for understanding the evolution of disease resistance in natural plant populations and for increasing crop yields in agriculture. In Brassica rapa, we studied genetic correlations among levels of disease resistance to three fungal pathogens: Peronospora parasitica, Albugo candida and Leptosphaeria maculans. A large, replicated quantitative genetics experiment used artificial selection on resistance to individual pathogens, and examined correlated responses to selection for resistance to other, unselected pathogens. Data from 9518 plants, each measured simultaneously for resistance to three fungal pathogens, showed heritable genetic variation for resistance to each pathogen and a positive genetic correlation between resistance to P. parasitica and L. maculans. This indicates that some resistance genes provide defence against fundamental characteristics common to two taxonomic orders of fungal pathogens. Conceivably, such MDR could contribute to a durable defence that might not be easily circumvented by rapidly evolving fungal pathogens.

Interval mapping of viability loci causing heterosis in Arabidopsis

The genetic basis of heterosis has implications for many problems in genetics and evolution. Heterosis and inbreeding depression affect human genetic diseases, maintenance of genetic variation, evolution of breeding systems, agricultural productivity, and conservation biology. Despite decades of theoretical and empirical studies, the genetic basis of heterosis has remained unclear. I mapped viability loci contributing to heterosis in Arabidopsis. An overdominant factor with large effects on viability mapped to a short interval on chromosome I. Homozygotes had 50% lower viability than heterozygotes in this chromosomal region. Statistical analysis of viability data in this cross indicates that observed viability heterosis is better explained by functional overdominance than by pseudo-overdominance. Overdominance sometimes may be an important cause of hybrid vigor, especially in habitually inbreeding species. Finally, I developed a maximum likelihood interval mapping procedure that can be used to examine chromosomal regions showing segregation distortion or viability selection.

Statistical genetics of an annual plant, Impatiens capensis. I. Genetic basis of quantitative variation

Analysis of quantitative genetics in natural populations has been hindered by computational and methodological problems in statistical analysis. We developed and validated a jackknife procedure to test for existence of broad sense heritabilities and dominance or maternal effects influencing quantitative characters in Impatiens capensis. Early life cycle characters showed evidence of dominance and/or maternal effects, while later characters exhibited predominantly environmental variation. Monte Carlo simulations demonstrate that these jackknife tests of variance components are extremely robust to heterogeneous error variances. Statistical methods from human genetics provide evidence for either a major locus influencing germination date, or genes that affect phenotypic variability per se. We urge explicit consideration of statistical behavior of estimation and testing procedures for proper biological interpretation of statistical results.