Evolution of microsatellites in Arabis petraea and Arabis lyrata, outcrossing relatives of Arabidopsis thaliana

We examined microsatellite variation in two diploid, outcrossing relatives of Arabidopsis thaliana, Arabis petraea and Arabis lyrata. The primer sequences were derived from A. thaliana. About 50% (14 loci) of the A. thaliana primers could successfully amplify microsatellites in the related species. Analysis of microsatellite structure in the related species showed that there had been large changes in the microsatellites: there were large differences in repeat numbers and many of the A. thaliana simple repeats were shorter in the related species. For the loci we compared, the related species had a much lower level of variability at the microsatellites than Japanese wild populations of A. thaliana. This is presumably related to the different microsatellite structures, because allozyme data showed that the outcrossing relatives were highly polymorphic compared to other outcrossing herbaceous species. Use of microsatellites in assessing variability or phylogenetic relationships between different species requires caution, because changes in microsatellite structure may alter evolutionary rates.      

Effects of gene expression on molecular evolution in Arabidopsis thaliana and Arabidopsis lyrata

We analyzed the complete genome sequence of Arabidopsis thaliana and sequence data from 83 genes in the outcrossing A. lyrata, to better understand the role of gene expression on the strength of natural selection on synonymous and replacement sites in Arabidopsis. From data on tRNA gene abundance, we find a good concordance between codon preferences and the relative abundance of isoaccepting tRNAs in the complete A. thaliana genome, consistent with models of translational selection. Both EST-based and new quantitative measures of gene expression (MPSS) suggest that codon preferences derived from information on tRNA abundance are more strongly associated with gene expression than those obtained from multivariate analysis, which provides further support for the hypothesis that codon bias in Arabidopsis is under selection mediated by tRNA abundance. Consistent with previous results, analysis of protein evolution reveals a significant correlation between gene expression level and amino acid substitution rate. Analysis by MPSS estimates of gene expression suggests that this effect is primarily the result of a correlation between the number of tissues in which a gene is expressed and the rate of amino acid substitution, which indicates that the degree of tissue specialization may be an important determinant of the rate of protein evolution in Arabidopsis.     

Evolution of Stress-Regulated Gene Expression in Duplicate Genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >~0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments.     

Centromere locations and associated chromosome rearrangements in Arabidopsis lyrata and A. thaliana

We analyzed linkage and chromosomal positions of genes in A. lyrata ssp. petraea that are located near the centromere (CEN) regions of A. thaliana, using at least two genes from the short and long arms of each chromosome. In our map, genes from all 10 A. thaliana chromosome arms are also tightly linked in A. lyrata. Genes from the regions on the two sides of CEN5 have distant map localizations in A. lyrata (genes on the A. thaliana short-arm genes are on linkage group AL6, and long-arm genes are on AL7), but genes from the other four A. thaliana centromere regions remain closely linked in A. lyrata. The observation of complete linkage between short- and long-arm centromere genes, but not between genes in other genome regions that are separated by similar physical distances, suggests that crossing-over frequencies near the A. lyrata ssp. petraea centromere regions are low, as in A. thaliana. Thus, the centromere positions appear to be conserved between A. thaliana and A. lyrata, even though three centromeres have been lost in A. thaliana, and the core satellite sequences in the two species are very different. We can now definitively identify the three centromeres that were eliminated in the fusions that formed the A. thaliana chromosomes. However, we cannot tell whether genes were lost along with these centromeres, because such genes are absent from the A. thaliana genome, which is the sole source of markers for our mapping.     

Gene targeting in Arabidopsis thaliana.

Summary Gene targeting of a chromosomally integrated transgene in Arabidopsis thaliana is reported. A chimeric gene consisting of the promoter of the 35S RNA of CaMV, the polyadenylation signal of the octopine synthase gene and the coding region of the bacterial hygromycin phosphotransferase gene (hpt), which was rendered non-functional by deletion of 19 bp, was introduced into the genome of A. thaliana using Agrobacterium-mediated gene transfer. A total of 3.46 x 10⁸ protoplasts isolated from 17 independent transgenic Arabidopsis lines harbouring the defective chimeric hpt gene were transformed via direct gene transfer using various DNA forms containing only the intact coding region of the hpt gene. Out of 150 hygromycin-resistant colonies appearing in the course of these experiments, four were the result of targeted recombination of the incoming DNA with the defective chromosomal locus as revealed by PCR and Southern blot analysis. Comparison with the number of transformants obtained when an hpt gene controlled by a promoter and terminator from the nopaline synthase gene was employed results in a maximal ratio of homologous to non-homologous transformation in A. thaliana of 1 x 10^–4.     

Chromosome arrangement and nuclear architecture but not centromeric sequences are conserved between Arabidopsis thaliana and Arabidopsis lyrata

In contrast to the situation described for mammals and Drosophila, chromosome territory (CT) arrangement and somatic homologous pairing in interphase nuclei of Arabidopsis thaliana (n = 5) are predominantly random except for a more frequent association of the chromosomes bearing a homologous nucleolus organizer region. To find out whether this chromosome arrangement is also characteristic for other species of the genus Arabidopsis, we investigated Arabidopsis lyrata ssp. lyrata (n = 8), one of the closest relatives of A. thaliana. First, we determined the size of each chromosome and chromosome arm, the sequence type of centromeric repeats and their distribution between individual centromeres and the position of the 5S/45S rDNA arrays in A. lyrata. Then we demonstrated that CT arrangement, homologous pairing and sister chromatid alignment of distinct euchromatic and/or heterochromatic regions within A. lyrata interphase nuclei are similar to that in A. thaliana nuclei. Thus, the arrangement of interphase chromosomes appears to be conserved between both taxa that diverged about 5 million years ago. Since the chromosomes of A. lyrata resemble those of the presumed ancestral karyotype, a similar arrangement of interphase chromosomes is also to be expected for other closely related diploid species of the Brassicaceae family.     

Linkage maps for Arabidopsis lyrata subsp. lyrata and Arabidopsis lyrata subsp. petraea combining anonymous and Arabidopsis thaliana-derived markers.

Arabidopsis lyrata, a close relative of the model plant Arabidopsis thaliana, is 1 of a few plant species for which the genome is to be entirely sequenced, which promises to yield important insights into genome evolution. Only 2 sparse linkage maps have been published, and these were based solely on markers derived from the A. thaliana genome. Because the genome of A. lyrata is practically twice as large as that of A. thaliana, the extent of map coverage of the A. lyrata genome remains uncertain. In this study, a 2-way pseudo-testcross strategy was used to construct genetic linkage maps of A. lyrata subsp. petraea and A. lyrata subsp. lyrata, using simple sequence repeat (SSR) and cleaved amplified polymorphic sequence (CAPS) markers from the A. thaliana genome, and anonymous amplified fragment length polymorphism (AFLP) markers that could potentially uncover regions unique to the A. lyrata genome. The SSR and CAPS markers largely confirmed the relationships between linkage groups in A. lyrata and A. thaliana. AFLP markers slightly increased the coverage of the A. lyrata maps, but mostly increased marker density on the linkage groups. We noted a much lower level of polymorphism and a greater segregation distortion in A. lyrata subsp. lyrata markers. The implications of these findings for the sequencing of the A. lyrata genome are discussed.     

Gene prediction and gene classes in Arabidopsis thaliana

Gene prediction methods for eukaryotic genomes still are not fully satisfying. One way to improve gene prediction accuracy, proven to be relevant for prokaryotes, is to consider more than one model of genes. Thus, we used our classification of Arabidopsis thaliana genes in two classes (CU(1) and CU(2)), previously delineated according to statistical features, in the GeneMark gene identification program. For each gene class, as well as for the two classes combined, a Markov model was developed (respectively, GM-CU(1), GM-CU(2) and GM-all) and then used on a test set of 168 genes to compare their respective efficiency. We concluded from this analysis that GM-CU(1) is more sensitive than GM-CU(2) which seems to be more specific to a gene type. Besides, GM-all does not give better results than GM-CU(1) and combining results from GM-CU(1) and GM-CU(2) greatly improve prediction efficiency in comparison with predictions made with GM-all only. Thus, this work confirms the necessity to consider more than one gene model for gene prediction in eukaryotic genomes, and to look for gene classes in order to build these models.     

Herbicide Safener-Inducible Gene Expression in Arabidopsis thaliana

The potential use of a new chemical-inducibie gene expressionsystem in Arabidopsis thaliana has been examined. The systemis based on the maize In2-2 promoter which is activated by benzenesulfonamideherbicide safe-ners. Plants transformed with the ß-glucuronidase(gus) reporter gene under the control of the In2-2 promoterwere grown in the presence of different safeners and the inducedGUS activity pattern was studied histochemically. In the absenceof safeners, the In2-2 promoter was not active. Applicationof different safeners induced distinct gus expression patterns,including expression in the root, hy-dathodes, and the shootapical meristem. Plants maintained continuously on inducingconcentrations of the safeners were retarded in growth. Thegrowth inhibition effects of the Sa5 safener could be overcomein a sul-fonylurea-resistant background. In2-2 promoter activitycould also be induced by the sulfonylurea herbicide chlor-sulfuron.In the sulfonylurea-resistant background, which derives fromherbicide-resistant acetolactate synthase activity, inductionof the In2-2 promoter by chlorsulfuron was lower. Furthermore,branched-chain amino acids, known to inhibit acetolactate synthaseactivity, also induced In2-2 promoter activity. Our data suggesta strong correlation between In2-2 expression and inhibitionof the acetolactate synthase activity. (Received November 12, 1996; Accepted February 21, 1997)     

Cesium Toxicity Alters MicroRNA Processing and AGO1 Expressions in Arabidopsis thaliana

MicroRNAs (miRNAs) are short RNA fragments that play important roles in controlled gene silencing, thus regulating many biological processes in plants. Recent studies have indicated that plants modulate miRNAs to sustain their survival in response to a variety of environmental stimuli, such as biotic stresses, cold, drought, nutritional starvation, and toxic heavy metals. Cesium and radio-cesium contaminations have arisen as serious problems that both impede plant growth and enter the food chain through contaminated plants. Many studies have been performed to define plant responses against cesium intoxication. However, the complete profile of miRNAs in plants during cesium intoxication has not been established. Here we show the differential expression of the miRNAs that are mostly down-regulated during cesium intoxication. Furthermore, we found that cesium toxicity disrupts both the processing of pri-miRNAs and AGONOUTE 1 (AGO1)-mediated gene silencing. AGO 1 seems to be especially destabilized by cesium toxicity, possibly through a proteolytic regulatory pathway. Our study presents a comprehensive profile of cesium-responsive miRNAs, which is distinct from that of potassium, and suggests two possible mechanisms underlying the cesium toxicity on miRNA metabolism.     

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.     

Evolution of chloroplast mononucleotide microsatellites in Arabidopsis thaliana

The level of variation and the mutation rate were investigated in an empirical study of 244 chloroplast microsatellites in 15 accessions of Arabidopsis thaliana. In contrast to SNP variation, microsatellite variation in the chloroplast was found to be common, although less common than microsatellite variation in the nucleus. No microsatellite variation was found in coding regions of the chloroplast. To evaluate different models of microsatellite evolution as possible explanations for the observed pattern of variation, the length distribution of microsatellites in the published DNA sequence of the A. thaliana chloroplast was subsequently used. By combining information from these two analyses we found that the mode of evolution of the chloroplast mononucleotide microsatellites was best described by a linear relation between repeat length and mutation rate, when the repeat lengths exceeded about 7 bp. This model can readily predict the variation observed in non-coding chloroplast DNA. It was found that the number of uninterrupted repeat units had a large impact on the level of chloroplast microsatellite variation. No other factors investigated—such as the position of a locus within the chromosome, or imperfect repeats—appeared to affect the variability of chloroplast microsatellites. By fitting the slippage models to the Genbank sequence of chromosome 1, we show that the difference between microsatellite variation in the nucleus and the chloroplast is largely due to differences in slippage rate. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

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.     

Selection-Driven Divergence After Gene Duplication in Arabidopsis thaliana

Gene duplications are one of the most important mechanisms for the origin of evolutionary novelties. Even though various models of the fate of duplicated genes have been established, current knowledge about the role of divergent selection after gene duplication is rather limited. In this study, we analyzed sequence divergence in response to neo- and subfunctionalization of segmentally duplicated genes in the genome of Arabidopsis thaliana. We compared the genomes of A. thaliana and the poplar Populus trichocarpa to identify orthologous pairs of genes and their corresponding inparalogs. Maximum-likelihood analyses of the nonsynonymous and synonymous substitution rate ratio [Formula: see text] of pairs of A. thaliana inparalogs were used to detect differences in the evolutionary rates of protein coding sequences. We analyzed 1,924 A. thaliana paralogous pairs and our results indicate that around 6.9% show divergent ω values between the lineages for a fraction of sites. We observe an enrichment of regulatory sequences, a reduced level of co-expression and an increased number of substitutions that can be attributed to positive selection based on an McDonald-Kreitman type of analysis. Taken together, these results show that selection after duplication contributes substantially to gene novelties and hence functional divergence in plants.