Subdivision and haplotype structure in natural populations of Arabidopsis lyrata

High levels of inbreeding are expected to cause a strong reduction in levels of genetic variability, effective recombination rates and in adaptation compared with related outcrossing populations. We examined patterns of DNA polymorphism at five nuclear loci and one chloroplast locus within and between four populations of the outcrossing plant Arabidopsis lyrata, a close relative of the highly self-fertilizing model species A. thaliana. The observed patterns are compared with species-wide polymorphism at orthologous loci, as well as within- and between-population patterns at other studied loci in A. thaliana. In addition to evidence for much higher average within-population diversity, species-wide levels of silent polymorphism are generally higher in A. lyrata than in A. thaliana, unlike the results from a previous study of the ADH locus. However, polymorphism is also low in the North American A. lyrata subspecies lyrata compared with the European subspecies petraea, suggesting either a population bottleneck in North American populations or recent admixture involving diverged European populations. Differentiation between the two subspecies is strong, although there are few fixed differences, suggesting that their isolation is recent. Estimates of intralocus recombination rates and analysis of haplotype structure in European A. lyrata populations indicate lower recombination than predicted based on the variability together with physical recombination rates estimated from A. thaliana. This may be due to strong population subdivision, or to recent departures from demographic equilibrium such as a bottleneck or population admixture. Alternatively, there may be consistently lower recombination rates in the outcrossing species. In contrast, estimates of recombination rates from species-wide samples of A. thaliana are close to the values expected assuming a high rate of self-fertilization. Complex population histories in both A. thaliana and A. lyrata complicate theoretical predictions and empirical tests of the effects of inbreeding on polymorphism and molecular evolution.     

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.     

Rates of deleterious mutation and the evolution of sex in Caenorhabditis

A variety of models propose that the accumulation of deleterious mutations plays an important role in the evolution of breeding systems. These models make predictions regarding the relative rates of protein evolution and deleterious mutation in taxa with contrasting modes of reproduction. Here we compare available coding sequences from one obligately outcrossing and two primarily selfing species of Caenorhabditis to explore the potential for mutational models to explain the evolution of breeding system in this clade. If deleterious mutations interact synergistically, the mutational deterministic hypothesis predicts that a high genomic deleterious mutation rate (U) will offset the reproductive disadvantage of outcrossing relative to asexual or selfing reproduction. Therefore, C. elegans and C. briggsae (both largely selfing) should both exhibit lower rates of deleterious mutation than the obligately outcrossing relative C. remanei. Using a comparative approach, we estimate U to be equivalent (and < 1) among all three related species. Stochastic mutational models, Muller's ratchet and Hill-Robertson interference, are expected to cause reductions in the effective population size in species that rarely outcross, thereby allowing deleterious mutations to accumulate at an elevated rate. We find only limited support for more rapid molecular evolution in selfing lineages. Overall, our analyses indicate that the evolution of breeding system in this group is unlikely to be explained solely by available mutational models.     

Higher intron loss rate in Arabidopsis thaliana than A. lyrata is consistent with stronger selection for a smaller genome

The number of introns varies considerably among different organisms. This can be explained by the differences in the rates of intron gain and loss. Two factors that are likely to influence these rates are selection for or against introns and the mutation rate that generates the novel intron or the intronless copy. Although it has been speculated that stronger selection for a compact genome might result in a higher rate of intron loss and a lower rate of intron gain, clear evidence is lacking, and the role of selection in determining these rates has not been established. Here, we studied the gain and loss of introns in the two closely related species Arabidopsis thaliana and A. lyrata as it was recently shown that A. thaliana has been undergoing a faster genome reduction driven by selection. We found that A. thaliana has lost six times more introns than A. lyrata since the divergence of the two species but gained very few introns. We suggest that stronger selection for genome reduction probably resulted in the much higher intron loss rate in A. thaliana, although further analysis is required as we could not find evidence that the loss rate increased in A. thaliana as opposed to having decreased in A. lyrata compared with the rate in the common ancestor. We also examined the pattern of the intron gains and losses to better understand the mechanisms by which they occur. Microsimilarity was detected between the splice sites of several gained and lost introns, suggesting that nonhomologous end joining repair of double-strand breaks might be a common pathway not only for intron gain but also for intron loss.     

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.     

Genome evolution among cruciferous plants: a lecture from the comparison of the genetic maps of three diploid species--Capsella rubella, Arabidopsis lyrata subsp. petraea, and A. thaliana

Comparative mapping in cruciferous plants is ongoing, and recently two additional genetic maps of diploid Capsella and Arabidopsis lyrata subsp. petraea have been presented. We compared both maps with each other using the sequence map and genomic data resources from Arabidopsis thaliana as a reference. The ancestors of the species pair Capsella-Arabidopsis diverged from one another approximately 10-14 million years ago (mya), whereas Arabidopsis thaliana and Arabidopsis lyrata have been separated since roughly 5-6 mya. Our analysis indicated that among diploid Capsella and Arabidopsis lyrata all eight genetic linkage groups are totally colinear to each other, with only two inversions significantly differentiating these two species.By minimizing the number of chromosomal rearrangements during genome evolution, we presented a model of chromosome evolution involving all three species. From this scenario, it is obvious that Arabidopsis thaliana underwent a dramatic genome reconstruction, with a base chromosome number reduction from five to eight and with approximately 1.3 chromosomal rearrangements per million years. In contrast, the terminal lineage leading to Capsella has only undergone less than 0.09 rearrangements per million years. This is the same rate as calculated for Arabidopsis lyrata since its separation from the Capsella lineage 10-14 mya. These results are in strong contrast to all overestimated rates calculated from comparisons of the systems Arabidopsis thaliana and Brassica, and our data demonstrate the problematic nature of both model systems.     

Population dynamics of an Ac-like transposable element in self- and cross-pollinating arabidopsis.

Theoretical models predict that the mating system should be an important factor driving the dynamics of transposable elements in natural populations due to differences in selective pressure on both element and host. We used a PCR-based approach to examine the abundance and levels of insertion polymorphism of Ac-III, a recently identified Ac-like transposon family, in natural populations of the selfing plant Arabidopsis thaliana and its close outcrossing relative, Arabidopsis lyrata. Although several insertions appeared to be ancient and shared between species, there is strong evidence for recent activity of this element family in both species. Sequences of the regions flanking insertions indicate that all Ac-III transposons segregating in natural populations are in noncoding regions and provide no evidence for local transposition events. Transposon display analysis suggests the presence of slightly higher numbers of insertion sites per individual but fewer total polymorphic insertions in the self-pollinating A. thaliana than A. lyrata. Element insertions appear to be segregating at significantly lower frequencies in A. lyrata than A. thaliana, which is consistent with a reduction in transposition rate, reduction in effective population size, or reduced efficacy of natural selection against element insertions in selfing populations.     

Patterns of molecular evolution in Caenorhabditis preclude ancient origins of selfing

The evolution of self-fertilization can mediate pronounced changes in genomes as a by-product of a drastic reduction in effective population size and the concomitant accumulation of slightly deleterious mutations by genetic drift. In the nematode genus Caenorhabditis, a highly selfing lifestyle has evolved twice independently, thus permitting an opportunity to test for the effects of mode of reproduction on patterns of molecular evolution on a genomic scale. Here we contrast rates of nucleotide substitution and codon usage bias among thousands of orthologous groups of genes in six species of Caenorhabditis, including the classic model organism Caenorhabditis elegans. Despite evidence that weak selection on synonymous codon usage is pervasive in the history of all species in this genus, we find little difference among species in the patterns of codon usage bias and in replacement-site substitution. Applying a model of relaxed selection on codon usage to the C. elegans and C. briggsae lineages suggests that self-fertilization is unlikely to have evolved more than approximately 4 million years ago, which is less than a quarter of the time since they shared a common ancestor with outcrossing species. We conclude that the profound changes in mating behavior, physiology, and developmental mechanisms that accompanied the transition from an obligately outcrossing to a primarily selfing mode of reproduction evolved in the not-too-distant past.     

Duplication of centromeric histone H3 (HTR12) gene in Arabidopsis halleri and A. lyrata, plant species with multiple centromeric satellite sequences

Arabidopsis halleri and lyrata have three different major centromeric satellite sequences, a unique finding for a diploid Arabidopsis species. Since centromeric histones coevolve with centromeric satellites, these proteins would be predicted to show signs of selection when new centromere satellites have recently arisen. We isolated centromeric protein genes from A. halleri and lyrata and found that one of them, HTR12 (CENP-A), is duplicated, while CENP-C is not. Phylogenetic analysis indicates that the HTR12 duplication occurred after these species diverged from A. thaliana. Genetic mapping shows that HTR12 copy B has the same genomic location as the A. thaliana gene; the other copy (A, at the other end of the same chromosome) is probably the new copy. To test for selection since the duplication, we surveyed diversity at both HTR12 loci within A. lyrata. Overall, there is no strong evidence for an "evolutionary arms race" causing multiple replacement substitutions. The A. lyrata HTR12B sequences fall into three classes of haplotypes, apparently maintained for a long time, but they all encode the same amino acid sequence. In contrast, HTR12A has low diversity, but many variants are amino acid replacements, possibly due to independent selective sweeps within populations of the species.     

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.     

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.      

Poorly known relatives of Arabidopsis thaliana

Non-model Arabidopsis species have been widely used as outgroup taxa in studies of molecular evolution. In Arabidopsis lyrata, Arabidopsis halleri and Arabidopsis arenosa, traits pertaining to self-incompatibility, heavy metal tolerance and inter-specific hybridization have been subjected to detailed genetic analysis. However, the full potential for exploring the causes and consequences of natural variation in complex traits within the genus Arabidopsis has not been widely appreciated or realized. Here, we draw on broadly dispersed information to characterize the basic biology, ecology, population genetics and molecular evolution for these three non-model Arabidopsis species. We illustrate how the wealth of functional and genomic tools pioneered in A. thaliana can be applied to gain insights into adaptive evolution of ecologically important traits and genome-wide processes, such as polyploidy, speciation and reticulate evolution, within and among Arabidopsis species.     

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.     

Intra- and interspecific DNA variation and codon bias of the alcohol dehydrogenase (Adh) locus in Arabis and Arabidopsis species

Sequence variation at the alcohol dehydrogenase (Adh) locus was analyzed for six species each of the genera Arabis and Arabidopsis. Phylogenetic analysis showed that investigated species were grouped into three clusters, and the generic classification did not correspond to the clusterings. The results indicated that the genera could not be distinguished on the basis of the Adh variation. A significant difference in the ratio of silent to replacement sites was detected by MK test in two comparisons, with Arabidopsis thaliana polymorphism due to excess silent divergence. Silent changes were predominant in the evolution of the Adh locus in Arabis and Arabidopsis. To infer evolutionary significance of silent substitutions, codon bias was studied. The degree of codon bias of the Adh region was relatively constant over Arabis and Arabidopsis species. "Preferred" codons of A. thaliana were determined. No evidence of natural selection on codon change was detected in the Adh regions of A. thaliana and Arabis gemmifera.     

Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera

As a model system, Arabidopsis thaliana and its wild relatives have played an important role in the study of genomics and evolution in plants. In this study, we examined the genetic diversity of the chalcone synthase (Chs) gene, which encodes a key enzyme of the flavonoid pathway and is located on chromosome five, as well as two Chs-like genes on the first and fourth chromosomes of Arabidopsis. The objectives of the study are to determine if natural selection operates differentially on the paralogs of the Chs gene family in A. thaliana and Arabidopsis halleri ssp. gemmifera. The mode of selection was inferred from Tajima's D values from noncoding and coding regions, as well as from the ratio of nonsynonymous to synonymous substitutions. Both McDonald-Kreitman and HKA tests revealed the effects of selection on the allelic distribution, except for the chromosome 1 paralog in ssp. gemmifera. The Chs gene on chromosome 5 was under purifying selection in both species. Significant, negative Tajima's D values at synonymous sites and positive Fay and Wu's H values within coding region, plus reduced genetic variability in introns, indicated effects of background selection in shaping the evolution of this gene region in A. thaliana. The Chs paralog on chromosome 1 was under positive selection in A. thaliana, while interspecific introgression and balancing selection determined the fates of the paralog and resulted in high heterogeneity in ssp. gemmifera. Local adaptation differentiated populations of Japan and China at the locus. In contrast, the other Chs-paralog of chromosome 4 was shaped by purifying selection in A. thaliana, while under positive selection in ssp. gemmifera, as indicated by dn/ds>1. Moreover, these contrasting patterns of selection have likely resulted in functional divergence in Arabidopsis, as indicated by radical amino acid substitutions at the chalcone synthase/stilbene synthase motif of the Chs genes. Unlike previous studies of the evolutionary history of A. thaliana, the high levels of genetic diversity in most gene regions of Chs paralogs and nonsignificant Tajima's D in the intron sequences of the Chs gene family in A. thaliana did not reflect the effects of a recent demographic expansion.     

Genome-wide comparison of nucleotide-binding site-leucine-rich repeat-encoding genes in Arabidopsis

Plants, like animals, use several lines of defense against pathogen attack. Prominent among genes that confer disease resistance are those encoding nucleotide-binding site-leucine-rich repeat (NB-LRR) proteins. Likely due to selection pressures caused by pathogens, NB-LRR genes are the most variable gene family in plants, but there appear to be species-specific limits to the number of NB-LRR genes in a genome. Allelic diversity within an individual is also increased by obligatory outcrossing, which leads to genome-wide heterozygosity. In this study, we compared the NB-LRR gene complement of the selfer Arabidopsis thaliana and its outcrossing close relative Arabidopsis lyrata. We then complemented and contrasted the interspecific patterns with studies of NB-LRR diversity within A. thaliana. Three important insights are as follows: (1) that both species have similar numbers of NB-LRR genes; (2) that loci with single NB-LRR genes are less variable than tandem arrays; and (3) that presence-absence polymorphisms within A. thaliana are not strongly correlated with the presence or absence of orthologs in A. lyrata. Although A. thaliana individuals are mostly homozygous and thus potentially less likely to suffer from aberrant interaction of NB-LRR proteins with newly introduced alleles, the number of NB-LRR genes is similar to that in A. lyrata. In intraspecific and interspecific comparisons, NB-LRR genes are also more variable than receptor-like protein genes. Finally, in contrast to Drosophila, there is a clearly positive relationship between interspecific divergence and intraspecific polymorphisms.     

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.     

Lineage‐specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the ‘edges’ of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage‐specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage‐specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage‐specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.