Update on sumoylation: defining core components of the plant SUMO conjugation system by phylogenetic comparison

The conjugation of the small ubiquitin-related modifier, SUMO, to substrate proteins is a reversible and dynamic process, and an important response of plants to environmental challenges. Nevertheless, reliable data have so far been restricted largely to the model plant Arabidopsis thaliana. The increasing availability of genome information for other plant species offers the possibility to identify a core set of indispensable components, and to discover species-specific features of the sumoylation pathway. We analyzed the enzymes responsible for the conjugation of SUMO to substrates for their conservation between dicots and monocots. We thus assembled gene sets that relate the Arabidopsis SUMO conjugation system to that of the dicot species tomato, grapevine and poplar, and to four plant species from the monocot class: rice, Brachypodium distachyon, Sorghum bicolor and maize. We found that a core set of genes with clear assignment in Arabidopsis had highly conserved homologs in all tested plants. However, we also observed a variation in the copy number of homologous genes, and sequence variations that suggested monocot-specific variants. Generally, SUMO ligases and proteases showed the most pronounced differences. Finally, we identified potential SUMO chain-binding ubiquitin ligases, pointing to an in vivo function of SUMO chains as degradation signals in plants.     

Dichotomy in the NRT gene families of dicots and grass species

A large proportion of the nitrate (NO(3)(-)) acquired by plants from soil is actively transported via members of the NRT families of NO(3)(-) transporters. In Arabidopsis, the NRT1 family has eight functionally characterised members and predominantly comprises low-affinity transporters; the NRT2 family contains seven members which appear to be high-affinity transporters; and there are two NRT3 (NAR2) family members which are known to participate in high-affinity transport. A modified reciprocal best hit (RBH) approach was used to identify putative orthologues of the Arabidopsis NRT genes in the four fully sequenced grass genomes (maize, rice, sorghum, Brachypodium). We also included the poplar genome in our analysis to establish whether differences between Arabidopsis and the grasses may be generally applicable to monocots and dicots. Our analysis reveals fundamental differences between Arabidopsis and the grass species in the gene number and family structure of all three families of NRT transporters. All grass species possessed additional NRT1.1 orthologues and appear to lack NRT1.6/NRT1.7 orthologues. There is significant separation in the NRT2 phylogenetic tree between NRT2 genes from dicots and grass species. This indicates that determination of function of NRT2 genes in grass species will not be possible in cereals based simply on sequence homology to functionally characterised Arabidopsis NRT2 genes and that proper functional analysis will be required. Arabidopsis has a unique NRT3.2 gene which may be a fusion of the NRT3.1 and NRT3.2 genes present in all other species examined here. This work provides a framework for future analysis of NO(3)(-) transporters and NO(3)(-) transport in grass crop species.     

A T-DNA mutation in the RNA helicase eIF4A confers a dose-dependent dwarfing phenotype in Brachypodium distachyon

In a survey of the BrachyTAG mutant population of Brachypodium distachyon, we identified a line carrying a T-DNA insertion in one of the two eukaryotic initiation factor 4A (eIF4A) genes present in the nuclear genome. The eif4a homozygous mutant plants were slow-growing, and exhibited reduced final plant stature due to a decrease in both cell number and cell size, consistent with roles for eIF4A in both cell division and cell growth. Hemizygous plants displayed a semi-dwarfing phenotype, in which stem length was reduced but leaf length was normal. Linkage between the insertion site and phenotype was confirmed, and we show that the level of eIF4A protein is strongly reduced in the mutant. Transformation of the Brachypodium homozygous mutant with a genomic copy of the Arabidopsis eIF4A-1 gene partially complemented the growth phenotype, indicating that gene function is conserved between mono- and dicotyledonous species. This study identifies eIF4A as a novel dose-dependent regulator of stem elongation, and demonstrates the utility of Brachypodium as a model for grass and cereals research.     

The phytochrome gene family in grasses (Poaceae): a phylogeny and evidence that grasses have a subset of the loci found in dicot angiosperms

The phytochrome nuclear gene family encodes photoreceptor proteins that mediate developmental responses to red and far red light throughout the life of the plant. From studies of the dicot flowering plant Arabidopsis, the family has been modeled as comprising five loci, PHYA- PHYE. However, it has been shown recently that the Arabidopsis model may not completely represent some flowering plant groups because additional PHY loci related to PHYA and PHYB of Arabidopsis apparently have evolved independently several times in dicots, and monocot flowering plants may lack orthologs of PHYD and PHYE of Arabidopsis. Nonetheless, the phytochrome nucleotide data were informative in a study of organismal evolution because the loci occur as single copy sequences and appear to be evolving independently. We have continued our investigation of the phytochrome gene family in flowering plants by sampling extensively in the grass family. The phytochrome nuclear DNA data were cladistically analyzed to address the following questions: (1) Are the data consistent with a pattern of differential distribution of phytochrome genes among monocots and higher dicots, with homologs of PHYA, B, C, D, and E present in higher dicots, but of just PHYA, B, and C in monocots, and (2) what phylogenetic pattern within Poaceae do they reveal? Results of these analyses, and of Southern blot experiments, are consistent with the observation that the phytochrome gene family in grasses comprises the same subset of loci detected in other monocots. Furthermore, for studies of organismal phylogeny in the grass family, the data are shown to provide significant support for relationships that are just weakly resolved by other data sets.      

Complete structure of the chloroplast genome of Arabidopsis thaliana.

The complete nucleotide sequence of the chloroplast genome of Arabidopsis thaliana has been determined. The genome as a circular DNA composed of 154,478 bp containing a pair of inverted repeats of 26,264 bp, which are separated by small and large single copy regions of 17,780 bp and 84,170 bp, respectively. A total of 87 potential protein-coding genes including 8 genes duplicated in the inverted repeat regions, 4 ribosomal RNA genes and 37 tRNA genes (30 gene species) representing 20 amino acid species were assigned to the genome on the basis of similarity to the chloroplast genes previously reported for other species. The translated amino acid sequences from respective potential protein-coding genes showed 63.9% to 100% sequence similarity to those of the corresponding genes in the chloroplast genome of Nicotiana tabacum, indicating the occurrence of significant diversity in the chloroplast genes between two dicot plants. The sequence data and gene information are available on the World Wide Web database KAOS (Kazusa Arabidopsis data Opening Site) at http://www.kazusa.or.jp/arabi/.     

Comparative Genomics of Flowering Time Pathways Using Brachypodium distachyon as a Model for the Temperate Grasses

Brachypodium distachyon (Brachypodium) is a model for the temperate grasses which include important cereals such as barley, wheat and oats. Comparison of the Brachypodium genome (accession Bd21) with those of the model dicot Arabidopsis thaliana and the tropical cereal rice (Oryza sativa) provides an opportunity to compare and contrast genetic pathways controlling important traits. We analysed the homologies of genes controlling the induction of flowering using pathways curated in Arabidopsis Reactome as a starting point. Pathways include those detecting and responding to the environmental cues of day length (photoperiod) and extended periods of low temperature (vernalization). Variation in these responses has been selected during cereal domestication, providing an interesting comparison with the wild genome of Brachypodium. Brachypodium Bd21 has well conserved homologues of circadian clock, photoperiod pathway and autonomous pathway genes defined in Arabidopsis and homologues of vernalization pathway genes defined in cereals with the exception of VRN2 which was absent. Bd21 also lacked a member of the CO family (CO3). In both cases flanking genes were conserved showing that these genes are deleted in at least this accession. Segmental duplication explains the presence of two CO-like genes in temperate cereals, of which one (Hd1) is retained in rice, and explains many differences in gene family structure between grasses and Arabidopsis. The conserved fine structure of duplications shows that they largely evolved to their present structure before the divergence of the rice and Brachypodium. Of four flowering-time genes found in rice but absent in Arabidopsis, two were found in Bd21 (Id1, OsMADS51) and two were absent (Ghd7, Ehd1). Overall, results suggest that an ancient core photoperiod pathway promoting flowering via the induction of FT has been modified by the recruitment of additional lineage specific pathways that promote or repress FT expression.     

Colinearity and gene density in grass genomes

Grasses are the single most important plant family in agriculture. In the past years, comparative genetic mapping has revealed conserved gene order (colinearity) among many grass species. Recently, the first studies at gene level have demonstrated that microcolinearity of genes is less conserved: small scale rearrangements and deletions complicate the microcolinearity between closely related species, such as sorghum and maize, but also between rice and other crop plants. In spite of these problems, rice remains the model plant for grasses as there is limited useful colinearity between Arabidopsis and grasses. However, studies in rice have to be complemented by more intensive genetic work on grass species with large genomes (maize, Triticeae). Gene-rich chromosomal regions in species with large genomes, such as wheat, have a high gene density and are ideal targets for partial genome sequencing.     

Adaptive evolution of Xa21 homologs in Gramineae

The XA21 protein has broad spectrum resistance against Xanthomonas oryzae pv. oryzae. Although Xa21-mediated immunity is well characterized, little is known about the origin and evolutionary history of this gene in grasses. Therefore, we analyzed all Xa21 gene homologs in eight whole-genome sequenced rice lines, as well as in four gramineous genomes, rice, Brachypodium, sorghum and maize; using Arabidopsis Xa21 homologs as outgroups, 17, 7, 7 and 3 Xa21 homologs were detected in these four grasses, respectively. Synteny and phylogenetic analysis showed that frequent gene translocation, duplication and/or loss, have occurred at Xa21 homologous loci, suggesting that they have undergone or are undergoing rapid generation of copy number variations. Within the rice species, the high level of nucleotide diversity between Xa21-like orthologs showed a strong association with the presence/absence haplotypes, suggesting that the genetic structure of rice lines plays an important role in the variations between these Xa21-like orthologs. Strongly positive selection was detected in the core region of the leucine-rich repeat domains of the Xa21 subclade among the rice lines, indicating that the rapid gene diversification of Xa21 homologs may be a strategy for a given species to adapt to the changing spectrum of species-specific pathogens.     

Phylogenetic analysis of the acetyl-CoA carboxylase and 3-phosphoglycerate kinase loci in wheat and other grasses

We have applied a two-gene system based on the sequences of nuclear genes encoding multi-domain plastid acetyl-CoA carboxylase (ACCase) and plastid 3-phosphoglycerate kinase (PGK) to study grass evolution. Our analysis revealed that these genes are single-copy in most of the grass species studied, allowing the establishment of orthologous relationships between them. These relationships are consistent with the known facts of their evolution: the eukaryotic origin of the plastid ACCase, created by duplication of a gene encoding the cytosolic multi-domain ACCase gene early in grass evolution, and the prokaryotic (endosymbiont) origin of the plastid PGK. The major phylogenetic relationships among grasses deduced from the nucleotide sequence comparisons of ACCase and PGK genes are consistent with each other and with the milestones of grass evolution revealed by other methods. Nucleotide substitution rates were calculated based on multiple pairwise sequence comparisons. On a relative basis, with the divergence of the Pooideae and Panicoideae subfamilies set at 60 million years ago (MYA), events leading to the Triticum/Aegilops complex occurred at the following intervals: divergence of Lolium (Lolium rigidum) at 35 MYA, divergence of Hordeum (Hordeum vulgare) at 11 MYA and divergence of Secale (Secale cereale) at 7 MYA. On the same scale, gene duplication leading to the multi-domain plastid ACCase in grasses occurred at 129 MYA, divergence of grass and dicot plastid PGK genes at 137 MYA, and divergence of grass and dicot cytosolic PGK genes at 155 MYA. The ACCase and PGK genes provide a well-understood two-locus system to study grass phylogeny, evolution and systematics.     

Micro-colinearity between rice, <Emphasis Type="Italic">Brachypodium</Emphasis>, and <Emphasis Type="Italic">Triticum monococcum</Emphasis> at the wheat domestication locus <Emphasis Type="Italic">Q</Emphasis>

Brachypodium, a wild temperate grass with a small genome, was recently proposed as a new model organism for the large-genome grasses. In this study, we evaluated gene content and microcolinearity between diploid wheat (Triticum monococcum), Brachypodium sylvaticum, and rice at a local genomic region harboring the major wheat domestication gene Q. Gene density was much lower in T. monococcum (one per 41 kb) because of gene duplication and an abundance of transposable elements within intergenic regions as compared to B. sylvaticum (one per 14 kb) and rice (one per 10 kb). For the Q gene region, microcolinearity was more conserved between wheat and rice than between wheat and Brachypodium because B. sylvaticum contained two genes apparently not present within the orthologous regions of T. monococcum and rice. However, phylogenetic analysis of Q and leukotriene A-4 hydrolase-like gene orthologs, which were colinear among the three species, showed that Brachypodium is more closely related to wheat than rice, which agrees with previous studies. We conclude that Brachypodium will be a useful tool for gene discovery, comparative genomics, and the study of evolutionary relationships among the grasses but will not preclude the need to conduct large-scale genomics experiments in the Triticeae.     

Analysis of sequence, map position, and gene expression reveals conserved essential genes for iron uptake in Arabidopsis and tomato

Arabidopsis (Arabidopsis thaliana) and tomato (Lycopersicon esculentum) show similar physiological responses to iron deficiency, suggesting that homologous genes are involved. Essential gene functions are generally considered to be carried out by orthologs that have remained conserved in sequence and map position in evolutionarily related species. This assumption has not yet been proven for plant genomes that underwent large genome rearrangements. We addressed this question in an attempt to deduce functional gene pairs for iron reduction, iron transport, and iron regulation between Arabidopsis and tomato. Iron uptake processes are essential for plant growth. We investigated iron uptake gene pairs from tomato and Arabidopsis, namely sequence, conserved gene content of the regions containing iron uptake homologs based on conserved orthologous set marker analysis, gene expression patterns, and, in two cases, genetic data. Compared to tomato, the Arabidopsis genome revealed more and larger gene families coding for the iron uptake functions. The number of possible homologous pairs was reduced if functional expression data were taken into account in addition to sequence and map position. We predict novel homologous as well as partially redundant functions of ferric reductase-like and iron-regulated transporter-like genes in Arabidopsis and tomato. Arabidopsis nicotianamine synthase genes encode a partially redundant family. In this study, Arabidopsis gene redundancy generally reflected the presumed genome duplication structure. In some cases, statistical analysis of conserved gene regions between tomato and Arabidopsis suggested a common evolutionary origin. Although involvement of conserved genes in iron uptake was found, these essential genes seem to be of paralogous rather than orthologous origin in tomato and Arabidopsis.     

Comparison of a high-density genetic linkage map to genome features in the model grass <Emphasis Type="Italic">Brachypodium distachyon</Emphasis>

The small annual grass Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to study questions unique to the grasses. Many Brachypodium resources have been developed including a whole genome sequence, highly efficient transformation and a large germplasm collection. We developed a genetic linkage map of Brachypodium using single nucleotide polymorphism (SNP) markers and an F₂ mapping population of 476 individuals. SNPs were identified by targeted resequencing of single copy genomic sequences. Using the Illumina GoldenGate Genotyping platform we placed 558 markers into five linkage groups corresponding to the five chromosomes of Brachypodium. The unusually long total genetic map length, 1,598 centiMorgans (cM), indicates that the Brachypodium mapping population has a high recombination rate. By comparing the genetic map to genome features we found that the recombination rate was positively correlated with gene density and negatively correlated with repetitive regions and sites of ancestral chromosome fusions that retained centromeric repeat sequences. A comparison of adjacent genome regions with high versus low recombination rates revealed a positive correlation between interspecific synteny and recombination rate.     

Lr34 multi-pathogen resistance ABC transporter: molecular analysis of homoeologous and orthologous genes in hexaploid wheat and other grass species

The Triticum aestivum (bread wheat) disease resistance gene Lr34 confers durable, race non-specific protection against three fungal pathogens, and has been a highly relevant gene for wheat breeding since the green revolution. Lr34, located on chromosome 7D, encodes an ATP-binding cassette (ABC) transporter. Both wheat cultivars with and without Lr34-based resistance encode a putatively functional protein that differ by only two amino acid polymorphisms. In this study, we focused on the identification and characterization of homoeologous and orthologous Lr34 genes in hexaploid wheat and other grasses. In hexaploid wheat we found an expressed and putatively functional Lr34 homoeolog located on chromosome 4A, designated Lr34-B. Another homoeologous Lr34 copy, located on chromosome 7A, was disrupted by the insertion of repetitive elements. Protein sequences of LR34-B and LR34 were 97% identical. Orthologous Lr34 genes were detected in the genomes of Oryza sativa (rice) and Sorghum bicolor (sorghum). Zea mays (maize), Brachypodium distachyon and Hordeum vulgare (barley) lacked Lr34 orthologs, indicating independent deletion of this particular ABC transporter. Lr34 was part of a gene-rich island on the wheat D genome. We found gene colinearity on the homoeologous A and B genomes of hexaploid wheat, but little microcolinearity in other grasses. The homoeologous LR34-B protein and the orthologs from rice and sorghum have the susceptible haplotype for the two critical polymorphisms distinguishing the LR34 proteins from susceptible and resistant wheat cultivars. We conclude that the particular Lr34-haplotype found in resistant wheat cultivars is unique. It probably resulted from functional gene diversification that occurred after the polyploidization event that was at the origin of cultivated bread wheat.     

An essential role of caffeoyl shikimate esterase in monolignol biosynthesis in Medicago truncatula

Biochemical and genetic analyses have previously identified caffeoyl shikimate esterase (CSE) as an enzyme in the monolignol biosynthesis pathway in Arabidopsis thaliana, although the generality of this finding has been questioned. Here we show the presence of CSE genes and associated enzyme activity in barrel medic (Medicago truncatula, dicot, Leguminosae), poplar (Populus deltoides, dicot, Salicaceae), and switchgrass (Panicum virgatum, monocot, Poaceae). Loss of function of CSE in transposon insertion lines of M. truncatula results in severe dwarfing, altered development, reduction in lignin content, and preferential accumulation of hydroxyphenyl units in lignin, indicating that the CSE enzyme is critical for normal lignification in this species. However, the model grass Brachypodium distachyon and corn (Zea mays) do not possess orthologs of the currently characterized CSE genes, and crude protein extracts from stems of these species exhibit only a weak esterase activity with caffeoyl shikimate. Our results suggest that the reaction catalyzed by CSE may not be essential for lignification in all plant species.