Contrasting patterns of nucleotide polymorphism at the alcohol dehydrogenase locus in the outcrossing Arabidopsis lyrata and the selfing Arabidopsis thaliana

Nucleotide variation at the alcohol dehydrogenase locus (Adh) was studied in the outcrossing Arabidopsis lyrata, a close relative of the selfing Arabidopsis thaliana. Overall, estimated nucleotide diversity in the North American ssp. lyrata and two European ssp. petraea populations was 0.0038, lower than the corresponding specieswide estimate for A. thaliana at the same set of nucleotide sites. The distribution of segregating sites across the gene differed between the two species. Estimated sequence diversity within an A. lyrata population with a large sample size (0.0023) was much higher than has previously been observed for A. thaliana. This North American population has an excess of sites at intermediate frequencies compared with neutral expectation (Tajima's D = 2.3, P < 0.005), suggestive of linked balancing selection or a recent population bottleneck. In contrast, an excess of rare polymorphisms has been found in A. thaliana. Polymorphism within A. lyrata and divergence from A. thaliana appear to be correlated across the Adh gene sequence. The geographic distribution of polymorphism was quite different from that of A. thaliana, for which earlier studies of several genes found low within-population nucleotide site polymorphism and no overall continental differentiation of variation despite large differences in site frequencies between local populations. Differences between the outcrossing A. lyrata and the selfing A. thaliana reflect the impact of differences in mating system and the influence of bottlenecks in A. thaliana during rapid colonization on DNA sequence polymorphism. The influence of additional variability-reducing mechanisms, such as background selection or hitchhiking, may not be discernible.     

Intraspecific genetic variations, fitness cost and benefit of RPW8, a disease resistance locus in Arabidopsis thaliana

The RPW8 locus of Arabidopsis thaliana confers broad-spectrum resistance to powdery mildew pathogens. In many A. thaliana accessions, this locus contains two homologous genes, RPW8.1 and RPW8.2. In some susceptible accessions, however, these two genes are replaced by HR4, a homolog of RPW8.1. Here, we show that RPW8.2 from A. lyrata conferred powdery mildew resistance in A. thaliana, suggesting that RPW8.2 might have gained the resistance function before the speciation of A. thaliana and A. lyrata. To investigate how RPW8 has been maintained in A. thaliana, we examined the nucleotide sequence polymorphisms in RPW8 from 51 A. thaliana accessions, related disease reaction phenotypes to the evolutionary history of RPW8.1 and RPW8.2, and identified mutations that confer phenotypic variations. The average nucleotide diversities were high at RPW8.1 and RPW8.2, showing no sign of selective sweep. Moreover, we found that expression of RPW8 incurs fitness benefits and costs on A. thaliana in the presence and absence of the pathogens, respectively. Our results suggest that polymorphisms at the RPW8 locus in A. thaliana may have been maintained by complex selective forces, including those from the fitness benefits and costs both associated with RPW8.     

Selection on amino acid substitutions in Arabidopsis

Studies of nucleotide diversity have found an excess of low-frequency amino acid polymorphisms segregating in Arabidopsis thaliana, suggesting a predominance of weak purifying selection acting on amino acid polymorphism in this inbreeding species. Here, we investigate levels of diversity and divergence at synonymous and nonsynonymous sites in 6 circumpolar populations of the outbreeding Arabidopsis lyrata and compare these results with A. thaliana, to test for differences in mutation and selection parameters across genes, populations, and species. We find that A. lyrata shows an excess of low-frequency nonsynonymous polymorphisms both within populations and species wide, consistent with weak purifying selection similar to the patterns observed in A. thaliana. Furthermore, nonsynonymous polymorphisms tend to be more restricted in their population distribution in A. lyrata, consistent with purifying selection preventing their geographic spread. Highly expressed genes show a reduced ratio of amino acid to synonymous change for both polymorphism and fixed differences, suggesting a general pattern of stronger purifying selection on high-expression proteins.     

Multilocus analysis of variation using a large empirical data set: phenylpropanoid pathway genes in Arabidopsis thaliana

Detecting the signature of adaptation on nucleotide variation is often difficult in species that like Arabidopsis thaliana might have a complex demographic history. Recent re-sequencing surveys in this species provided genome-wide information that would mainly reflect its demographic history. We have used a large empirical data set (LED) as well as multilocus coalescent simulations to analyse sequence variation at loci involved in the phenylpropanoid pathway of this species. We surveyed and examined DNA sequence variation at nine of these loci (about 19.7 kb) in 23 accessions of A. thaliana and one accession of its closely related species Arabidopsis lyrata . Nucleotide variation was lower at nonsynonymous sites than at silent sites in all loci, indicating generalized functional constraint at the protein level. No association between variation and position in the metabolic pathway was detected. When the data were contrasted against the standard neutral model, significant deviations for silent variation were detected with Tajima's D , Fu's F_S and Fay and Wu's H multilocus test statistics. These deviations were in the same direction than in previous large-scale multilocus analyses, suggesting a genome-wide effect. When the nine-locus data set was contrasted against the large empirical data set, the level (Watterson's θ) and pattern of variation (Tajima's D ) detected in these loci did not deviate either at the single-locus or multilocus level from the corresponding empirical distributions. These results would support an important role of the demographic history of A. thaliana in shaping nucleotide variation at the nine studied phenylpropanoid loci. The potential and limitations of the empirical distribution approach are discussed.     

Nucleotide variation at the CHALCONE ISOMERASE locus in Arabidopsis thaliana

An approximately 1.9-kb region encompassing the CHI gene, which encodes chalcone isomerase, was sequenced in 24 worldwide ecotypes of Arabidopsis thaliana (L.) Heynh. and in 1 ecotype of A. lyrata ssp. petraea. There was no evidence for dimorphism at the CHI region. A minimum of three recombination events was inferred in the history of the sampled ecotypes of the highly selfing A. thaliana. The estimated nucleotide diversity theta(TOTAL) = 0.004, theta(SIL) = 0. 005 was on the lower part of the range of the corresponding estimates for other gene regions. The skewness of the frequency spectrum toward an excess of low-frequency polymorphisms, together with the bell-shaped distribution of pairwise nucleotide differences at CHI, suggests that A. thaliana has recently experienced a rapid population growth. Although this pattern could also be explained by a recent selective sweep at the studied region, results from the other studied loci and from an AFLP survey seem to support the expansion hypothesis. Comparison of silent polymorphism and divergence at the CHI region and at the Adh1 and ChiA revealed in some cases a significant deviation of the direct relationship predicted by the neutral theory, which would be compatible with balancing selection acting at the latter regions.     

Cell Specific,Cross-Species Expression of Myrosinases in Brassica Napus,Arabidopsis Thaliana and Nicotiana Tabacum

A prototypical characteristic of the Brassicaceae is the presence of the myrosinase-glucosinolate system. Myrosinase, the only known S-glycosidase in plants, degrades glucosinolates, thereby initiating the formation of isothiocyanates, nitriles and other reactive products with biological activities. We have used myrosinase gene promoters from Brassica napus and Arabidopsis thaliana fused to the beta -glucuronidase (GUS) reporter gene and introduced into Arabidopsis thaliana, Brassica napus and/or Nicotiana tabacum plants to compare and determine the cell types expressing the myrosinase genes and the GUS expression regulated by these promoters. The A. thaliana TGG1 promoter directs expression to guard cells and phloem myrosin cell idioblasts of transgenic A. thaliana plants. Expression from the same promoter construct in transgenic tobacco plants lacking the myrosinase enzyme system also directs expression to guard cells. The B. napus Myr1.Bn1 promoter directs a cell specific expression to idioblast myrosin cells of immature and mature seeds and myrosin cells of phloem of B. napus. In A. thaliana the B. napus promoter directs expression to guard cells similar to the expression pattern of TGG1. The Myr1.Bn1 signal peptide targets the gene product to the reticular myrosin grains of myrosin cells. Our results indicate that myrosinase gene promoters from Brassicaceae direct cell, organ and developmental specific expression in B. napus, A. thaliana and N. tabacum.     

Comparative gene mapping in Arabidopsis lyrata chromosomes 1 and 2 and the corresponding A. thaliana chromosome 1: recombination rates, rearrangements and centromere location

To add detail to the genetic map of Arabidopsis lyrata, and compare it with that of A. thaliana, we have developed many additional markers in the A. lyrata linkage groups, LG1 and LG2, corresponding to A. thaliana chromosome 1. We used a newly developed method for marker development for single nucleotide polymorphisms present in gene sequences, plus length differences, to map genes in an A. lyrata family, including variants in several genes close to the A. thaliana centromere 1, providing the first data on the location of an A. lyrata centromere; we discuss the implications for the evolution of chromosome 1 of A. thaliana. With our larger marker density, large rearrangements between the two Arabidopsis species are excluded, except for a large inversion on LG2. This was previously known in Capsella; its presence in A. lyrata suggests that, like most other rearrangements, it probably arose in the A. thaliana lineage. Knowing that marker orders are similar, we can now compare homologous, non-rearranged map distances to test the prediction of more frequent crossing-over in the more inbreeding species. Our results support the previous conclusion of similar distances in the two species for A. lyrata LG1 markers. For LG2 markers, the distances were consistently, but non-significantly, larger in A. lyrata. Given the two species' large DNA content difference, the similarity of map lengths, particularly for LG1, suggests that crossing-over is more frequent across comparable physical distances in the inbreeder, A. thaliana, as predicted.     

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 相似文献   

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

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

Birth, death and subfunctionalization in the Arabidopsis genome

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

Effective population size and tests of neutrality at cytoplasmic genes in Arabidopsis.

Cytoplasmic genomes typically lack recombination, implying that genetic hitch-hiking could be a predominant force structuring nucleotide polymorphism in the chloroplast and mitochondria. We test this hypothesis by analysing nucleotide polymorphism data at 28 loci across the chloroplast and mitochondria of the outcrossing plant Arabidopsis lyrata, and compare patterns with multiple nuclear loci, and the highly selfing Arabidopsis thaliana. The maximum likelihood estimate of the ratio of effective population size at cytoplasmic relative to nuclear genes in A. lyrata does not depart from the neutral expectation of 0.5. Similarly, the ratio of effective size in A. thaliana is close to unity, the neutral expectation for a highly selfing species. The results are thus consistent with neutral organelle polymorphism in these species or with comparable effects of hitch-hiking in both cytoplasmic and nuclear genes, in contrast to the results of recent studies on gynodioecious taxa. The four-gamete test and composite likelihood estimation provide evidence for very low levels of recombination in the organelles of A. lyrata, although permutation tests do not suggest that adjacent polymorphic sites are more closely linked than more distant sites across the two genomes, suggesting that mutation hotspots or very low rates of gene conversion could explain the data.     

Multilocus analysis of variation and speciation in the closely related species Arabidopsis halleri and A. lyrata

Nucleotide variation in eight effectively unlinked genes was surveyed in species-wide samples of the closely related outbreeding species Arabidopsis halleri and A. lyrata ssp. petraea and in three of these genes in A. lyrata ssp. lyrata and A. thaliana. Significant genetic differentiation was observed more frequently in A. l. petraea than in A. halleri. Average estimates of nucleotide variation were highest in A. l. petraea and lowest in A. l. lyrata, reflecting differences among species in effective population size. The low level of variation in A. l. lyrata is concordant with a bottleneck effect associated with its origin. The A. halleri/A. l. petraea speciation process was studied, considering the orthologous sequences of an outgroup species (A. thaliana). The high number of ancestral mutations relative to exclusive polymorphisms detected in A. halleri and A. l. petraea, the significant results of the multilocus Fay and Wu H tests, and haplotype sharing between the species indicate introgression subsequent to speciation. Average among-population variation in A. halleri and A. l. petraea was approximately 1.5- and 3-fold higher than that in the inbreeder A. thaliana. The detected reduction of variation in A. thaliana is less than that expected from differences in mating system alone, and therefore from selective processes related to differences in the effective recombination rate, but could be explained by differences in population structure.     

Guard cell- and phloem idioblast-specific expression of thioglucoside glucohydrolase 1 (myrosinase) in Arabidopsis

Thioglucoside glucohydrolase 1 (TGG1) is one of two known functional myrosinase enzymes in Arabidopsis. The enzyme catalyzes the hydrolysis of glucosinolates into compounds that are toxic to various microbes and herbivores. Transgenic Arabidopsis plants carrying beta-glucuronidase and green fluorescent protein reporter genes fused to 0.5 or 2.5 kb of the TGG1 promoter region were used to study spatial promoter activity. Promoter activity was found to be highly specific and restricted to guard cells and distinct cells of the phloem. No promoter activity was detected in the root or seed. All guard cells show promoter activity. Positive phloem cells are distributed in a discontinuous pattern and occur more frequent in young tissues. Immunocytochemical localization of myrosinase in transverse and longitudinal sections of embedded material show that the TGG1 promoter activity reflects the position of the myrosinase enzyme. In the flower stalk, the myrosinase-containing phloem cells are located between phloem sieve elements and glucosinolate-rich S cells. Our results suggest a cellular separation of myrosinase enzyme and glucosinolate substrate, and that myrosinase is contained in distinct cells. We discuss the potential advantages of locating defense and communication systems to only a few specific cell types.     

Genetic variation within and among populations of Arabidopsis thaliana.

We investigated levels of nucleotide polymorphism within and among populations of the highly self-fertilizing Brassicaceous species, Arabidopsis thaliana. Four-cutter RFLP data were collected at one mitochondrial and three nuclear loci from 115 isolines representing 11 worldwide population collections, as well as from seven commonly used ecotypes. The collections include multiple populations from North America and Eurasia, as well as two pairs of collections from locally proximate sites, and thus allow a hierarchical geographic analysis of polymorphism. We found no variation at the mitochondrial locus Nad5 and very low levels of intrapopulation nucleotide diversity at Adh, Dhs1, and Gpa1. Interpopulation nucleotide diversity was also consistently low among the loci, averaging 0.0014. gst, a measure of population differentiation, was estimated to be 0.643. Interestingly, we found no association between geographical distance between populations and genetic distance. Most haplotypes have a worldwide distribution, suggesting a recent expansion of the species or long-distance gene flow. The low level of polymorphism found in this study is consistent with theoretical models of neutral mutations and background selection in highly self-fertilizing species.     

Self-incompatibility in the genus Arabidopsis: characterization of the S locus in the outcrossing A. lyrata and its autogamous relative A. thaliana

As a starting point for a phylogenetic study of self-incompatibility (SI) in crucifers and to elucidate the genetic basis of transitions between outcrossing and self-fertilizing mating systems in this family, we investigated the SI system of Arabidopsis lyrata. A. lyrata is an outcrossing close relative of the self-fertile A. thaliana and is thought to have diverged from A. thaliana approximately 5 million years ago and from Brassica spp 15 to 20 million years ago. Analysis of two S (sterility) locus haplotypes demonstrates that the A. lyrata S locus contains tightly linked orthologs of the S locus receptor kinase (SRK) gene and the S locus cysteine-rich protein (SCR) gene, which are the determinants of SI specificity in stigma and pollen, respectively, but lacks an S locus glycoprotein gene. As described previously in Brassica, the S haplotypes of A. lyrata differ by the rearranged order of their genes and by their variable physical sizes. Comparative mapping of the A. lyrata and Brassica S loci indicates that the S locus of crucifers is a dynamic locus that has undergone several duplication events since the Arabidopsis--Brassica split and was translocated as a unit between two distant chromosomal locations during diversification of the two taxa. Furthermore, comparative analysis of the S locus region of A. lyrata and its homeolog in self-fertile A. thaliana identified orthologs of the SRK and SCR genes and demonstrated that self-compatibility in this species is associated with inactivation of SI specificity genes.     

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

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

Arabidopsis myrosinases TGG1 and TGG2 have redundant function in glucosinolate breakdown and insect defense

In Arabidopsis and other Brassicaceae, the enzyme myrosinase (beta-thioglucoside glucohydrolase, TGG) degrades glucosinolates to produce toxins that deter herbivory. A broadly applicable selection for meiotic recombination between tightly linked T-DNA insertions was developed to generate Arabidopsis tgg1tgg2 double mutants and study myrosinase function. Glucosinolate breakdown in crushed leaves of tgg1 or tgg2 single mutants was comparable to that of wild-type, indicating redundant enzyme function. In contrast, leaf extracts of tgg1tgg2 double mutants had undetectable myrosinase activity in vitro, and damage-induced breakdown of endogenous glucosinolates was apparently absent for aliphatic and greatly slowed for indole glucosinolates. Maturing leaves of myrosinase mutants had significantly increased glucosinolate levels. However, developmental decreases in glucosinolate content during senescence and germination were unaffected, showing that these processes occur independently of TGG1 and TGG2. Insect herbivores with different host plant preferences and feeding styles varied in their responses to myrosinase mutations. Weight gain of two Lepidoptera, the generalist Trichoplusia ni and the facultative Solanaceae-specialist Manduca sexta, was significantly increased on tgg1tgg2 double mutants. Two crucifer-specialist Lepidoptera had differing responses. Whereas Plutella xylostella was unaffected by myrosinase mutations, Pieris rapae performed better on wild-type, perhaps due to reduced feeding stimulants in tgg1tgg2 mutants. Reproduction of two Homoptera, Myzus persicae and Brevicoryne brassicae, was unaffected by myrosinase mutations.