Genome-level evolution of resistance genes in Arabidopsis thaliana

Pathogen resistance genes represent some of the most abundant and diverse gene families found within plant genomes. However, evolutionary mechanisms generating resistance gene diversity at the genome level are not well understood. We used the complete Arabidopsis thaliana genome sequence to show that most duplication of individual NBS-LRR sequences occurs at close physical proximity to the parent sequence and generates clusters of closely related NBS-LRR sequences. Deploying the statistical strength of phylogeographic approaches and using chromosomal location as a proxy for spatial location, we show that apparent duplication of NBS-LRR genes to ectopic chromosomal locations is largely the consequence of segmental chromosome duplication and rearrangement, rather than the independent duplication of individual sequences. Although accounting for a smaller fraction of NBS-LRR gene duplications, segmental chromosome duplication and rearrangement events have a large impact on the evolution of this multigene family. Intergenic exchange is dramatically lower between NBS-LRR sequences located in different chromosome regions as compared to exchange between sequences within the same chromosome region. Consequently, once translocated to new chromosome locations, NBS-LRR gene copies have a greater likelihood of escaping intergenic exchange and adopting new functions than do gene copies located within the same chromosomal region. We propose an evolutionary model that relates processes of genome evolution to mechanisms of evolution for the large, diverse, NBS-LRR gene family.     

Metallochaperone-like genes in Arabidopsis thaliana

A complete inventory of metallochaperone-like proteins containing a predicted HMA domain in Arabidopsis revealed a large family of 67 proteins. 45 proteins, the HIPPs, have a predicted isoprenylation site while 22 proteins, the HPPs, do not. Sequence comparisons divided the proteins into seven major clusters (I-VII). Cluster IV is notable for the presence of a conserved Asp residue before the CysXXCys, metal binding motif, analogous to the Zn binding motif in E. coli ZntA. HIPP20, HIPP21, HIPP22, HIPP26 and HIPP27 in Cluster IV were studied in more detail. All but HIPP21 could rescue the Cd-sensitive, ycf1 yeast mutant but failed to rescue the growth of zrt1zrt2, zrc1cot1 and atx1 mutants. In Arabidopsis, single and double mutants did not show a phenotype but the hipp20/21/22 triple mutant was more sensitive to Cd and accumulated less Cd than the wild-type suggesting the HIPPs can have a role in Cd-detoxification, possibly by binding Cd. Promoter-GUS reporter expression studies indicated variable expression of these HIPPs. For example, in roots, HIPP22 and HIPP26 are only expressed in lateral root tips while HIPP20 and HIPP25 show strong expression in the root vasculature.     

Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A. thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the ‘core’ proteome of green plants; (iii) 70% of A. thaliana genes could be traced back to 450 million years ago; and (iv) intriguingly, A. thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome‐wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A. thaliana.     

Identifying essential genes in Arabidopsis thaliana

Eight years after publication of the Arabidopsis genome sequence and two years before completing the first phase of an international effort to characterize the function of every Arabidopsis gene, plant biologists remain unable to provide a definitive answer to the following basic question: what is the minimal gene set required for normal growth and development? The purpose of this review is to summarize different strategies employed to identify essential genes in Arabidopsis, an important component of the minimal gene set in plants, to present an overview of the datasets and specific genes identified to date, and to discuss the prospects for future saturation of this important class of genes. The long-term goal of this collaborative effort is to facilitate basic research in plant biology and complement ongoing research with other model organisms.     

Targeted disruption of the TGA3 locus in Arabidopsis thaliana

A major drawback to study gene functions in plant systems is the lack of an effective gene knockout strategy. With a large number of plant genes isolated and the accelerating pace by which this collection is growing, the need for their functional analyses at the whole plant level has become increasingly urgent. Here evidence is reported for the first successful disruption of a non-selectable gene in Arabidopsis thaliana by creating a mutant of the TGA3 locus via targeted insertion of the bacterial neo gene conferring kanamycin (Km) resistance. A β-glucuronidase (GUS) expression unit outside the region of homology was used as a screenable marker to distinguish homologous recombination events from those of ectopic insertions. PCR amplification coupled with Southern blot screening identified two putative homologous recombination events among 2580 Km^r calli. One callus line was subsequently isolated and the structure of the targeted TGA3 allele confirmed by Southern blot analyses. This study demonstrates the feasibility of targeting a non-selectable locus in Arabidopsis. Combined with future improvements in negative selection strategies and efficient, transformation methodologies, gene replacement studies in plants could become a routine technique.     

Sequence and organization of 5S ribosomal RNA-encoding genes of Arabidopsis thaliana.

We have isolated a genomic clone containing Arabidopsis thaliana 5S ribosomal RNA (rRNA)-encoding genes (rDNA) by screening an A. thaliana library with a 5S rDNA probe from flax. The clone isolated contains seven repeat units of 497 bp, plus 11 kb of flanking genomic sequence at one border. Sequencing of individual subcloned repeat units shows that the sequence of the 5S rRNA coding region is very similar to that reported for other flowering plants. Four A. thaliana ecotypes were found to contain approx. 1000 copies of 5S rDNA per haploid genome. Southern-blot analysis of genomic DNA indicates that 5S rDNA occurs in long tandem arrays, and shows the presence of numerous restriction-site polymorphisms among the six ecotypes studied.     

Mode of amplification and reorganization of resistance genes during recent Arabidopsis thaliana evolution.

The NBS-LRR (nucleotide-binding site plus leucine-rich repeat) genes represent the major class of disease resistance genes in flowering plants and comprise 166 genes in the ecotype Col-0 of Arabidopsis thaliana. NBS-LRR genes are organized in single-gene loci, clusters, and superclusters. Phylogenetic analysis reveals nine monophyletic clades and a few phylogenetic orphans. Most clusters contain only genes from the same phylogenetic lineage, reflecting their origin from the exchange of sequence blocks as a result of intralocus recombination. Multiple duplications increased the number of NBS-LRR genes in the progenitors of Arabidopsis, suggesting that the present complexity in Col-0 may derive from as few as 17 progenitors. The combination of physical and phylogenetic analyses of the NBS-LRR genes makes it possible to detect relatively recent gene rearrangements, which increased the number of NBS-LRR genes by about 50, but which are almost never associated with large segmental duplications. The identification of 10 heterogeneous clusters containing members from different clades demonstrates that sequence sampling between different resistance gene loci and clades has occurred. Such events may have taken place early during flowering plant evolution, but they generated modules that have been duplicated and remobilized also more recently.     

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.     

Genetic structure and evolution of RAC-GTPases in Arabidopsis thaliana

Rho GTPases regulate a number of important cellular functions in eukaryotes, such as organization of the cytoskeleton, stress-induced signal transduction, cell death, cell growth, and differentiation. We have conducted an extensive screening, characterization, and analysis of genes belonging to the Ras superfamily of GTPases in land plants (embryophyta) and found that the Rho family is composed mainly of proteins with homology to RAC-like proteins in terrestrial plants. Here we present the genomic and cDNA sequences of the RAC gene family from the plant Arabidopsis thaliana. On the basis of amino acid alignments and genomic structure comparison of the corresponding genes, the 11 encoded AtRAC proteins can be divided into two distinct groups of which one group apparently has evolved only in vascular plants. Our phylogenetic analysis suggests that the plant RAC genes underwent a rapid evolution and diversification prior to the emergence of the embryophyta, creating a group that is distinct from rac/cdc42 genes in other eukaryotes. In embryophyta, RAC genes have later undergone an expansion through numerous large gene duplications. Five of these RAC duplications in Arabidopsis thaliana are reported here. We also present an hypothesis suggesting that the characteristic RAC proteins in higher plants have evolved to compensate the loss of RAS proteins.     

Receptor-like protein kinase genes of Arabidopsis thaliana

The isolation of a maize cDNA clone that encodes a membrane spanning protein kinase related to the self-incompatibility glycoproteins (SLG) of Brassica and structurally similar to the growth factor receptor tyrosine kinases has recently been reported. Three distinct receptor-like protein kinase (RLK) cDNA clones from Arabidopsis thaliana have now been identified. Two of the Arabidopsis RLK genes encode SLG-related protein kinases but have different patterns of expression: one is expressed predominantly in rosettes while the other is expressed primarily in roots. The third RLK gene contains an extracellular domain that consists of 21 leucine-rich repeats that are analogous to the leucine-rich repeats found in proteins from humans, flies and yeast. The Arabidopsis leucine-rich gene is expressed at equivalent levels in roots and rosettes. These results show that there are several genes in higher plants that encode members of the receptor protein kinase superfamily. The structural diversity and differential expression of these genes suggest that each plays a distinct and possibly important role in cellular signaling in plants.     

Distribution of 5-methylcytosine residues in 5S rRNA genes in Arabidopsis thaliana and Secale cereale

Bisulfite genomic sequencing was used to localise 5-methylcytosine residues (mC) in 5S rRNA genes of Arabidopsis thaliana and Secale cereale. The maps of mC distribution were compared with the previously published map of the corresponding region in Nicotiana tabacum. In all three species, the level of methylation of 5S rRNA genes was generally higher than the average for the entire genome. The ratio of 5S rDNA methylation to average overall methylation was 44%/30-33% for N. tabacum, 27%/4-6% for A. thaliana and 24%/20-22% for S. cereale. With the exception of one clone from S. cereale, no methylation-free 5S rDNA was detected. The level of methylation at different sequence motifs in 5S rDNA was calculated for N. tabacum/A. thaliana/ S. cereale, and this analysis yielded the following values (expressed as a percentage of total C): mCG 90%/78%/85%, mCWG 89%/41%/53%, mCmCG 72%/32%/16%, mCCG 4%/2%/0%, mCHH 15%/6%/1%, where W=A or T, and H=A or C or T. Non-symmetrical methylation was almost negligible in the large genome of S. cereale but relatively frequent in N. tabacum and A. thaliana, suggesting that the strict correlation between genome size and cytosine methylation might be violated for this type of methylation. Among non-symmetrical motifs the mCWA triplets were significantly over-represented in Arabidopsis, while in tobacco this preference was not as pronounced. The differences in methylation levels in different sequence contexts might be of phylogenetic significance, but further species in related and different taxa need to be studied before firm conclusions can be drawn.     

Leaf-variegated mutations and their responsible genes in Arabidopsis thaliana

Leaf variegation has long been known as a recessive genetic trait in higher plants. Unlike albino mutants, leaf-variegated mutants are non-lethal and thus enable us to study a novel mechanism of plastid development and maintenance. Variegation results from a defect that makes chloroplast development unstable, since at least part of the tissues gives rise to normal chloroplasts. Despite the fact that leaf-variegated mutants have contributed to the findings of maternal inheritance or have been used as genetic markers, these mutations and the responsible loci have been poorly understood at the molecular level. A comprehensive study of the leaf-variegated mutants is possible in Arabidopsis, since such mutants have been known and the cloning can be at relative ease as a model plant. Here I summarize recent progress on characterization of the Arabidopsis leaf-variegated mutants. Detailed analysis of the responsible loci revealed that variegation is caused by a defect in various metabolic pathways related to organelle functions. Thus, studies on these genes provide us with novel redundant mechanisms by which heteroplasmic organelles such as plastids and mitochondria can survive from an environmental stress.     

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.     

Fine mapping a self-fertility locus in perennial ryegrass

Key message

A self-fertility locus was fine mapped to a 1.6 cM region on linkage group 5 in a perennial ryegrass population. This locus was the main determinant of pollen self-compatibility.

Abstract

In grasses, self-incompatibility (SI) is characterized by a two-loci gametophytic (S and Z) mechanism acting together in the recognition and inhibition of self-pollen. Mutations affecting the expression of SI have been reported in a few grass species. In perennial ryegrass (Lolium perenne L.), a mutation independent from S and Z, and mapping on linkage group 5 (LG 5), was previously reported to produce self-fertile plants. Here, we describe fine mapping of the self-fertility (SF) gene in a perennial ryegrass population and determine whether there is any effect of other genomic regions on the pollen compatibility. The phenotypic segregation of SF showed a bimodal distribution with one mean at 49% pollen compatibility and the other at 91%. Marker-trait association analysis showed that only markers on LG 5 were significantly associated with the trait. A single gene model explained 82% of the observed variability and no effects of the other regions were detected. Using segregation and linkage analysis, the SF locus was located to a 1.6 cM region on LG 5. The flanking marker sequences were aligned to rice and Brachypodium distachyon reference genomes to estimate the physical distance. We provide markers tightly linked to SF that can be used for introgression of this trait into advanced breeding germplasm. Moreover, our results represent a further step towards the identification of the SF gene in LG 5.