A Microarray Based Genomic Hybridization Method for Identification of New Genes in Plants: Case Analyses of Arabidopsis and Oryza

To systematically estimate the gene duplication events in closely related species, we have to use comparative genomic approaches, either through genomic sequence comparison or comparative genomic hybridization (CGH). Given the scarcity of complete genomic sequences of plant species, in the present study we adopted an array based CGH to investigate gene duplications in the genus Arabidopsis. Fragment genomic DNA from four species, namely Arabidopsis thaliana, A. lyrata subsp. lyrata, A. lyrata subsp. petraea, and A. halleri, was hybridized to Affymetrix (Santa Clara, CA, USA) tiling arrays that are designed from the genomic sequences of A. thaliana. Pairwise comparisons of signal intensity were made to infer the potential duplicated candidates along each phylo-genetic branch. Ninety-four potential candidates of gene duplication along the genus were identified. Among them, the majority (69 of 94) were A. thaliana lineage specific. This result indicates that the array based CGH approach may be used to identify candidates of duplication in other plant genera containing closely related species, such as Oryza, particularly for the AA genome species. We compared the degree of gene duplication through retrotransposon between O. sativa and A. thaliana and found a strikingly higher number of chimera retroposed genes in rice. The higher rate of gene duplication through retroposition and other mechanisms may indicate that the grass species is able to adapt to more diverse environments.     

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.     

Does Speciation between Arabidopsis halleri and Arabidopsis lyrata Coincide with Major Changes in a Molecular Target of Adaptation?

Ever since Darwin proposed natural selection as the driving force for the origin of species, the role of adaptive processes in speciation has remained controversial. In particular, a largely unsolved issue is whether key divergent ecological adaptations are associated with speciation events or evolve secondarily within sister species after the split. The plant Arabidopsis halleri is one of the few species able to colonize soils highly enriched in zinc and cadmium. Recent advances in the molecular genetics of adaptation show that the physiology of this derived ecological trait involves copy number expansions of the AhHMA4 gene, for which orthologs are found in single copy in the closely related A. lyrata and the outgroup A. thaliana. To gain insight into the speciation process, we ask whether adaptive molecular changes at this candidate gene were contemporary with important stages of the speciation process. We first inferred the scenario and timescale of speciation by comparing patterns of variation across the genomic backgrounds of A. halleri and A. lyrata. Then, we estimated the timing of the first duplication of AhHMA4 in A. halleri. Our analysis suggests that the historical split between the two species closely coincides with major changes in this molecular target of adaptation in the A. halleri lineage. These results clearly indicate that these changes evolved in A. halleri well before industrial activities fostered the spread of Zn- and Cd-polluted areas, and suggest that adaptive processes related to heavy-metal homeostasis played a major role in the speciation process.     

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.     

Ecotype-specific and chromosome-specific expansion of variant centromeric satellites in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Despite the conserved roles and conserved protein machineries of centromeres, their nucleotide sequences can be highly diverse even among related species. The diversity reflects rapid evolution, but the underlying mechanism is largely unknown. One approach to monitor rapid evolution is examination of intra-specific variation. Here we report variant centromeric satellites of Arabidopsis thaliana found through survey of 103 natural accessions (ecotypes). Among them, a cluster of variant centromeric satellites was detected in one ecotype, Cape Verde Islands (Cvi). Recombinant inbred mapping revealed that the variant satellites are distributed in centromeric region of the chromosome 5 (CEN5) of this ecotype. This apparently recent variant accumulation is associated with large deletion of a pericentromeric region and the expansion of satellite region. The variant satellite was bound to HTR12 (centromeric variant histone H3), although expansion of the satellite was not associated with comparable increase in the HTR12 binding. The results suggest that variant satellites with centromere function can rapidly accumulate in one centromere, supporting the model that the satellite repeats in the array are homogenized by occasional unequal crossing-over, which has a potential to generate an expansion of local sequence variants within a centromere cluster. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

The transition to self-compatibility in Arabidopsis thaliana and evolution within S-haplotypes over 10 Myr

A recent investigation found evidence that the transition of Arabidopsis thaliana from ancestral self-incompatibility (SI) to full self-compatibility occurred very recently and suggested that this occurred through a selective fixation of a nonfunctional allele (PsiSCR1) at the SCR gene, which determines pollen specificity in the incompatibility response. The main evidence is the lack of polymorphism at the SCR locus in A. thaliana. However, the nearby SRK gene, which determines stigma specificity in self-incompatible Brassicaceae species, has extremely high sequence diversity, with 3 very divergent SRK haplotypes, 2 of them present in multiple strains. Such high diversity is extremely unusual in this species, and it suggests the possibility that multiple, different SRK haplotypes may have been preserved from A. thaliana's self-incompatible ancestor. To study the evolution of S-haplotypes in the A. thaliana lineage, we searched the 2 most closely related Arabidopsis species Arabidopsis lyrata and Arabidopsis halleri, in which most populations have retained SI, and found SRK sequences corresponding to all 3 A. thaliana haplogroup sequences. Our molecular evolutionary analyses of these 3 S-haplotypes provide an independent estimate of the timing of the breakdown of SI and again exclude an ancient transition to selfing in A. thaliana. Comparing sequences of each of the 3 haplogroups between species, we find that 2 of the 3 SRK sequences (haplogroups A and B) are similar throughout their length, suggesting that little or no recombination with other SRK alleles has occurred since these species diverged. The diversity difference between the SCR and SRK loci in A. thaliana, however, suggests crossing-over, either within SRK or between the SCR and SRK loci. If the loss of SI involved fixation of the PsiSCR1 sequence, the exchange must have occurred during its fixation. Divergence between the species is much lower at the S-locus, compared with reference loci, and we discuss two contributory possibilities. Introgression may have occurred between A. lyrata and A. halleri and between their ancestral lineage and A. thaliana, at least for some period after their split. In addition, the coalescence times of sequences of individual S-haplogroups are expected to be less than those of alleles at non-S-loci.     

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.     

The allopolyploid Arabidopsis kamchatica originated from multiple individuals of Arabidopsis lyrata and Arabidopsis halleri

Polyploidization, or genome duplication, has played a critical role in the diversification of animals, fungi and plants. Little is known about the population structure and multiple origins of polyploid species because of the difficulty in identifying multiple homeologous nuclear genes. The allotetraploid species Arabidopsis kamchatica is closely related to the model species Arabidopsis thaliana and is distributed in a broader climatic niche than its parental species. Here, we performed direct sequencing of homeologous pairs of the low-copy nuclear genes WER and CHS by designing homeolog-specific primers, and obtained also chloroplast and ribosomal internal transcribed spacer sequences. Phylogenetic analysis showed that 50 individuals covering the distribution range including North America are allopolyploids derived from Arabidopsis lyrata and Arabidopsis halleri . Three major clusters within A. kamchatica were detected using Bayesian clustering. One cluster has widespread distribution. The other two are restricted to the southern part of the distribution range including Japan, where the parent A. lyrata is not currently distributed. This suggests that the mountains in Central Honshu and surrounding areas in Japan served as refugia during glacial–interglacial cycles and retained this diversity. We also found that multiple haplotypes of nuclear and chloroplast sequences of A. kamchatica are identical to those of their parental species. This indicates that multiple diploid individuals contributed to the origin of A. kamchatica . The haplotypes of low-copy nuclear genes in Japan suggest independent polyploidization events rather than introgression. Our findings suggest that self-compatibility and gene silencing occurred independently in different origins.     

Centromeric localization and adaptive evolution of an Arabidopsis histone H3 variant

Centromeric H3-like histones, which replace histone H3 in the centromeric chromatin of animals and fungi, have not been reported in plants. We identified a histone H3 variant from Arabidopsis thaliana that encodes a centromere-identifying protein designated HTR12. By immunological detection, HTR12 localized at centromeres in both mitotic and meiotic cells. HTR12 signal revealed tissue- and stage-specific differences in centromere morphology, including a distended bead-like structure in interphase root tip cells. The anti-HTR12 antibody also detected spherical organelles in meiotic cells. Although the antibody does not label centromeres in the closely related species Arabidopsis arenosa, HTR12 signal was found on all centromeres in allopolyploids of these two species. Comparison of the HTR12 genes of A. thaliana and A. arenosa revealed striking adaptive evolution in the N-terminal tail of the protein, similar to the pattern seen in its counterpart in Drosophila. This finding suggests that the same evolutionary forces shape centromeric chromatin in both animals and plants.     

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.     

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.     

The genetic basis of zinc tolerance in the metallophyte Arabidopsis halleri ssp. halleri (Brassicaceae): an analysis of quantitative trait loci

The species Arabidopsis halleri, an emerging model for the study of heavy metal tolerance and accumulation in plants, has evolved a high level of constitutive zinc tolerance. Mapping of quantitative trait loci (QTL) was used to investigate the genetic architecture of zinc tolerance in this species. A first-generation backcross progeny of A. halleri ssp. halleri from a highly contaminated industrial site and its nontolerant relative A. lyrata ssp. petraea was produced and used for QTL mapping of zinc tolerance. A genetic map covering most of the A. halleri genome was constructed using 85 markers. Among these markers, 65 were anchored in A. thaliana and revealed high synteny with other Arabidopsis genomes. Three QTL of comparable magnitude on three different linkage groups were identified. At all QTL positions zinc tolerance was enhanced by A. halleri alleles, indicating directional selection for higher zinc tolerance in this species. The two-LOD support intervals associated with these QTL cover 24, 4, and 13 cM. The importance of each of these three regions is emphasized by their colocalization with HMA4, MTP1-A, and MTP1-B, respectively, three genes well known to be involved in metal homeostasis and tolerance in plants.     

Linkage disequilibrium between incompatibility locus region genes in the plant Arabidopsis lyrata

We have studied diversity in Arabidopsis lyrata of sequences orthologous to the ARK3 gene of A. thaliana. Our main goal was to test for recombination in the S-locus region. In A. thaliana, the single-copy ARK3 gene is closely linked to the non-functional copies of the self-incompatibility loci, and the ortholog in A. lyrata (a self-incompatible species) is in the homologous genome region and is known as Aly8. It is thus of interest to test whether Aly8 sequence diversity is elevated due to close linkage to the highly polymorphic incompatibility locus, as is theoretically predicted. However, Aly8 is not a single-copy gene, and the presence of paralogs could also lead to the appearance of elevated diversity. We established a typing approach based on different lengths of Aly8 PCR products and show that most A. lyrata haplotypes have a single copy, but some have two gene copies, both closely linked to the incompatibility locus, one being a pseudogene. We determined the phase of multiple haplotypes in families of plants from Icelandic and other populations. Different Aly8 sequence types are associated with different SRK alleles, while haplotypes with the same SRK sequences tend to have the same Aly8 sequence. There is evidence of some exchange of sequences between different Aly8 sequences, making it difficult to determine which ones are allelic or to estimate the diversity. However, the homogeneity of the Aly8 sequences of each S-haplotype suggests that recombination between the loci has been very infrequent over the evolutionary history of these populations. Overall, the results suggest that recombination rarely occurs in the interval between the S-loci and Aly8 and that linkage to the S-loci can probably account for the observed high Aly8 diversity.     

Zinc tolerance and accumulation in metallicolous and nonmetallicolous populations of Arabidopsis halleri (Brassicaceae)

Zinc tolerance was investigated in five populations of Arabidopsis halleri (syn.: Cardaminopsis halleri ) raised from seeds collected from contaminated and uncontaminated sites. Tolerance was measured by determining the concentration which inhibited root growth (EC₁₀₀). A. halleri populations from contaminated and uncontaminated sites were found to be Zn-tolerant compared with the Zn-nontolerant species Arabidopsis thaliana and A. lyrata subsp. petraea . At very high Zn concentrations, populations of A. halleri from uncontaminated sites were slightly less Zn-tolerant than those from contaminated sites. These observations support the hypothesis that in A. halleri , Zn tolerance is largely a constitutive property. One population from an uncontaminated site and one population from a contaminated site were studied for Zn uptake. Zinc content was measured in shoots and roots using a colorimetric test under laboratory conditions. The results showed that whatever their origin, individuals from both populations are Zn accumulators compared with the nonaccumulator species A. thaliana . Moreover, the population from the uncontaminated area accumulated Zn in its shoots and roots more quickly than the population from the contaminated site. These results suggest that, in A. halleri , Zn accumulation to very high concentration is a constitutive property.     

Genome wide gene expression in artificially synthesized amphidiploids of Arabidopsis

The merging of two different genomes occurs during the formation of amphidiploids, and the merged regulatory networks have the potential to generate a new gene expression pattern. We examined the genome-wide gene expression of two newly synthesized amphidiploids between Arabidopsis thaliana and the related species Arabidopsis lyrata subsp. lyrata and Arabidopsis halleri subsp. gemmifera. 1,137 (4.7%) and 1,316 (5.4%) of probesets showed differential gene expression in A. thaliana-A. halleri and A. thaliana-A. lyrata hybrids respectively, compared to the mid parent value and of these, 489 were in common. Genes that differed in expression between the parental lines tended to have an expression level in both hybrids differing from the mid parent value. In contrast to protein coding genes, there is little differential expression of transposons. Genes in the categories of chloroplast-targeted and response to stress were overrepresented in the non-additively expressed genes in both amphidiploids. As these genes have the potential to contribute directly to the plant phenotype, we suggest that rapid changes of gene expression in amphidiploids might be important for producing greater biomass.     

High diversity due to balancing selection in the promoter region of the Medea gene in Arabidopsis lyrata

Molecular imprinting is the differential expression and/or silencing of alleles according to their parent of origin [1, 2]. Conflicts between parents, or parents and offspring, should cause "arms races," with accelerated evolution of the genes involved in imprinting. This should be detectable in the evolution of imprinting genes' protein sequences and in the promoter regions of imprinted genes. Previous studies, however, found no evidence of more amino acid substitutions in imprinting genes [1, 3]. We have analyzed sequence diversity of the Arabidopsis lyrata Medea (MEA) gene and divergence from the A. thaliana sequence, including the first study of the promoter region. In A. thaliana, MEA is imprinted, with paternal alleles silenced in endosperm cells [4, 5], and also functions in the imprinting machinery [4, 6]; MEA protein binding at the MEA promoter region indicates self-regulated imprinting [7-9]. We find the same paternal MEA allele silencing in A. lyrata endosperm but no evidence for adaptive evolution in the coding region, whereas the 5' flanking region displays high diversity, with distinct haplotypes, suggesting balancing selection in the promoter region.     

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.