Rate and polarity of gene fusion and fission in Oryza sativa and Arabidopsis thaliana

Eukaryotic gene fusion and fission events are mechanistically more complicated than in prokaryotes, and their quantitative contributions to genome evolution are still poorly understood. We have identified all differentially composite or split genes in 2 fully sequenced plant genomes, Oryza sativa and Arabidopsis thaliana. Out of 10,172 orthologous gene pairs, 60 (0.6% of the total) revealed a verified fusion or fission event in either lineage after the divergence of O. sativa and A. thaliana. Polarizing these events by outgroup comparison revealed differences in the rate of gene fission but not of gene fusion in the rice and Arabidopsis lineages. Gene fission occurred at a higher rate than gene fusion in the O. sativa lineage and was furthermore more common in rice than in Arabidopsis. Nucleotide insertion bias has promoted gene fission in the O. sativa lineage, consistent with its generally longer nucleotide sequences than A. thaliana in selectively neutral regions, and with the abundance of transposable elements in rice. The divergence time of monocots and dicots (140-200 Myr) indicates that gene fusion/fission events occur at an average rate of 1x10(-11) to 2x10(-11) events per gene per year, approximately 100-fold slower than the average per site nuclear nucleotide substitution rate in these lineages. Gene fusion and fission are thus rare and slow processes in higher plant genomes; they should be of utility to address deeper evolutionary relationships among plants--and the relationship of plants to other eukaryotic lineages--where sequence-based phylogenies provide equivocal or conflicting results.     

Isolation and characterization of conserved non-coding sequences among rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) paralogous regions

Segmental duplication is particularly frequent within plant genomes and the ability of the original single-copy gene to gain a new function for the change of regulatory elements is one of the prominent consequences of duplication. Thus, it is important to study the pattern of conserved non-coding sequence (CNS) between paralogous genes. We report the result of a survey of CNSs among paralogous regions in rice (Oryza sativa L.), as well as the comparison of CNS dataset between rice and Arabidopsis thaliana. Some common properties, such as the change of A + T content near the CNS boundaries and CNS are enriched in regulatory genes, were observed. However, the content of CNSs differs between rice and Arabidopsis, and it is interesting that the rice metabolic network includes both CNS-poor and CNS-rich genes, which indicated a fine-tuned metabolic network presents in rice.     

Exploring the plant transcriptome through phylogenetic profiling

Publicly available protein sequences represent only a small fraction of the full catalog of genes encoded by the genomes of different plants, such as green algae, mosses, gymnosperms, and angiosperms. By contrast, an enormous amount of expressed sequence tags (ESTs) exists for a wide variety of plant species, representing a substantial part of all transcribed plant genes. Integrating protein and EST sequences in comparative and evolutionary analyses is not straightforward because of the heterogeneous nature of both types of sequence data. By combining information from publicly available EST and protein sequences for 32 different plant species, we identified more than 250,000 plant proteins organized in more than 12,000 gene families. Approximately 60% of the proteins are absent from current sequence databases but provide important new information about plant gene families. Analysis of the distribution of gene families over different plant species through phylogenetic profiling reveals interesting insights into plant gene evolution, and identifies species- and lineage-specific gene families, orphan genes, and conserved core genes across the green plant lineage. We counted a similar number of approximately 9,500 gene families in monocotyledonous and eudicotyledonous plants and found strong evidence for the existence of at least 33,700 genes in rice (Oryza sativa). Interestingly, the larger number of genes in rice compared to Arabidopsis (Arabidopsis thaliana) can partially be explained by a larger amount of species-specific single-copy genes and species-specific gene families. In addition, a majority of large gene families, typically containing more than 50 genes, are bigger in rice than Arabidopsis, whereas the opposite seems true for small gene families.     

Mutational bias affects protein evolution in flowering plants

Amino acid sequences from several thousand homologous gene pairs were compared for two plant genomes, Oryza sativa and Arabidopsis thaliana. The Arabidopsis genes all have similar G+C (guanine plus cytosine) contents, whereas their homologs in rice span a wide range of G+C levels. The results show that those rice genes that display increased divergence in their nucleotide composition (specifically, increased G+C content) showed a corresponding, predictable change in the amino acid compositions of the encoded proteins relative to their Arabidopsis homologs. This trend was not seen in a "control" set of rice genes that had nucleotide contents closer to their Arabidopsis homologs. In addition to showing an overall difference in the amino acid composition of the homologous proteins, we were also able to investigate the biased patterns of amino acid substitution since the divergence of these two species. We found that the amino acid exchange matrix was highly asymmetric when comparing the High G+C rice genes with their Arabidopsis homologs. Finally, we investigated the possible causes of this biased pattern of sequence evolution. Our results indicate that the biased pattern of protein evolution is the consequence, rather than the cause, of the corresponding changes in nucleotide content. In fact, there is an even more marked asymmetry in the patterns of substitution at synonymous nucleotide sites. Surprisingly, there is a very strong negative correlation between the level of nucleotide bias and the length of the coding sequences within the rice genome. This difference in gene length may provide important clues about the underlying mechanisms.     

Ultraconserved elements between the genomes of the plants Arabidopsis thaliana and rice

Ultraconserved elements, sequences with 100% identity with no insertions or deletions between genomes, have been found in both vertebrate and invertebrate genomes; whether plant genomes contain ultraconserved elements, however, is unknown. We consequently compared the genomes of Arabidopsis thaliana and rice, which diverged about 200 million years ago, and identified 25 ultraconserved elements that are longer than 100 bp. Similar to those previously found, ultraconserved elements in plants tend to occur in clusters and locate at noncoding regions; nevertheless, they have many distinct features. For instance, the longest ultraconserved element between the 2 plant genomes is 1491 bp, much longer than the longest one (779 bp) between the human and rodent genomes. Some biological implications are discussed, but the functions of these plant ultraconserved elements and the reasons why they are practically frozen during the evolution of millions of years remain a mystery.     

Genome-wide analysis of plant glutaredoxin systems

The recent release of the first tree genome (Populus trichocarpa) has allowed a comparison to be made of the multigenic glutaredoxin (Grx) and glutathione reductase (GR) families of this tree with those of other sequenced organisms and especially of the two other fully sequenced plant species, Arabidopsis thaliana and Oryza sativa. Grxs are small proteins involved in disulphide bridge or protein-glutathione adduct reduction, and they are maintained in a reduced form using glutathione and an NADPH-dependent GR. While the P. trichocarpa and O. sativa genomes are nearly five times larger than that of A. thaliana, they contain approximately 45 000 and 37 500 genes compared with the 25 500 genes of A. thaliana. On the one hand, the GR gene composition varies little between species and the gene structures are relatively conserved. On the other hand, the Grx gene family can be divided into three subgroups and the gene content is larger in P. trichocarpa (36 genes) compared with A. thaliana and O. sativa (31 and 27 genes, respectively). This could be partly explained by the occurrence of more duplication events, and this is especially true for one of the three identified Grx subgroups (subgroup III). The expression of most of these genes was confirmed by analysing expressed sequence tags present in various databases. In addition, the expression of Grx of subgroups I and II was examined by RT-PCR in various poplar organs. A complete classification based essentially on gene structure and sequence identity is proposed.     

Genome-wide comparative analysis of putative bidirectional promoters from rice, Arabidopsis and Populus

A bidirectional promoter can regulate the expression of two flanking genes arranged in a divergent manner. Although reports pertaining to bidirectional promoters on a genomic scale exist in mammals, little progress has been made in plants. In the present study, we performed a computational analysis of this unique class of promoters to identify overrepresented cis-regulatory motifs from three sequenced plant genomes: rice (Oryza sativa), Arabidopsis thaliana, and Populus trichocarpa using the Plant Cis-acting Regulatory DNA Elements (PLACE) and PLANT CARE databases. We describe these overrepresented elements and their possible regulatory mechanisms. We also discuss similarities and differences with human bidirectional promoters. Furthermore, we describe in detail a few coexpressed and evolutionarily conserved divergent gene pairs and their bidirectional promoters. This study provides insights into bidirectional promoters in three plant species, thereby laying a foundation for their experimental analysis.     

MUSTANG is a novel family of domesticated transposase genes found in diverse angiosperms

While transposons have traditionally been viewed as genomic parasites or "junk DNA," the discovery of transposon-derived host genes has fueled an ongoing debate over the evolutionary role of transposons. In particular, while mobility-related open reading frames have been known to acquire host functions, the contribution of these types of events to the evolution of genes is not well understood. Here we report that genome-wide searches for Mutator transposase-derived host genes in Arabidopsis thaliana (Columbia-0) and Oryza sativa ssp. japonica (cv. Nipponbare) (domesticated rice) identified 121 sequences, including the taxonomically conserved MUSTANG1. Syntenic MUSTANG1 orthologs in such varied plant species as rice, poplar, Arabidopsis, and Medicago truncatula appear to be under purifying selection. However, despite the evidence of this pathway of gene evolution, MUSTANG1 belongs to one of only two Mutator-like gene families with members in both monocotyledonous and dicotyledonous plants, suggesting that Mutator-like elements seldom evolve into taxonomically widespread host genes.     

Genome evolution in diploid and tetraploid Coffea species as revealed by comparative analysis of orthologous genome segments

Sequence comparison of orthologous regions enables estimation of the divergence between genomes, analysis of their evolution and detection of particular features of the genomes, such as sequence rearrangements and transposable elements. Despite the economic importance of Coffea species, little genomic information is currently available. Coffea is a relatively young genus that includes more than one hundred diploid species and a single tetraploid species. Three Coffea orthologous regions of 470-900 kb were analyzed and compared: both subgenomes of allotetraploid Coffea arabica (contributed by the diploid species Coffea eugenioides and Coffea canephora) and the genome of diploid C. canephora. Sequence divergence was calculated on global alignments or on coding and non-coding sequences separately. A search for transposable elements detected 43 retrotransposons and 198 transposons in the sequences analyzed. Comparative insertion analysis made it possible to locate 165 TE insertions in the phylogenetic tree of the three genomes/subgenomes. In the tetraploid C. arabica, a homoeologous non-reciprocal transposition (HNRT) was detected and characterized: a 50 kb region of the C. eugenioides derived subgenome replaced the C. canephora derived counterpart. Comparative sequence analysis on three Coffea genomes/subgenomes revealed almost perfect gene synteny, low sequence divergence and a high number of shared transposable elements. Compared to the results of similar analysis in other genera (Aegilops/Triticum and Oryza), Coffea genomes/subgenomes appeared to be dramatically less diverged, which is consistent with the relatively recent radiation of the Coffea genus. Based on nucleotide substitution frequency, the HNRT was dated at 10,000-50,000 years BP, which is also the most recent estimation of the origin of C. arabica.     

Characterization of T-DNA insertion sites in Arabidopsis thaliana and the implications for saturation mutagenesis

A key component of a sound functional genomics infrastructure is the availability of a knockout mutant for every gene in the genome. A fruitful approach to systematically knockingout genes in the plant Arabidopsis thaliana has been the use of transferred-DNA (T-DNA) from Agrobacterium tumefaciens as an insertional mutagen. One of the assumptions underlying the use of T-DNA as a mutagen is that the insertion of these DNA elements into the Arabidopsis genome occurs at randomly selected locations. We have directly investigated the distribution of T-DNA insertions sites in populations of transformed Arabidopsis using two different approaches. To begin with, we utilized a polymerase chain reaction (PCR) procedure to systematically catalog the precise locations of all the T-DNA elements inserted within a 65 kb segment of chromosome IV. Of the 47 T-DNA insertions identified, 30% were found within the coding regions of genes. We also documented the insertion of T-DNA elements within the centromeric region of chromosome IV. In addition to these targeted T-DNA screens, we also mapped the genomic locations of 583 randomly chosen T-DNA elements by sequencing the genomic DNA flanking the insertion sites from individual T-DNA-transformed lines. 35% of these randomly chosen T-DNA insertions were located within the coding regions of genes. For comparison, coding sequences account for 44% of the Arabidopsis genome. Our results demonstrate that there is a small bias towards recovering T-DNA insertions within intergenic regions. However, this bias does not limit the utility of T-DNA as an effective insertional mutagen for use in reverse-genetic strategies.     

The evolution of RNA editing and pentatricopeptide repeat genes

The pentatricopeptide repeat (PPR) is a degenerate 35-amino-acid structural motif identified from analysis of the sequenced genome of the model plant Arabidopsis thaliana. From the wealth of sequence information now available from plant genomes, the PPR protein family is now known to be one of the largest families in angiosperm species, as most genomes encode 400-600 members. As the number of PPR genes is generally only c. 10-20 in other eukaryotic organisms, including green algae, the family has obviously greatly expanded during land plant evolution. This provides a rare opportunity to study selection pressures driving a 50-fold expansion of a single gene family. PPR proteins are sequence-specific RNA-binding proteins involved in many aspects of RNA processing in organelles. In this review, we will summarize our current knowledge about the evolution of PPR genes, and will discuss the relevance of the dramatic expansion in the family to the functional diversification of plant organelles, focusing primarily on RNA editing.     

Comparative analysis of divergent and convergent gene pairs and their expression patterns in rice, Arabidopsis, and populus

Comparative analysis of the organization and expression patterns of divergent and convergent gene pairs in multiple plant genomes can identify patterns that are shared by more than one species or are unique to a particular species. Here, we study the coexpression and interspecies conservation of divergent and convergent gene pairs in three plant species: rice (Oryza sativa), Arabidopsis (Arabidopsis thaliana), and black cottonwood (Populus trichocarpa). Strongly correlated expression levels between divergent and convergent genes were found to be quite common in all three species, and the frequency of strong correlation appears to be independent of intergenic distance. Conservation of divergent or convergent arrangement among these species appears to be quite rare. However, conserved arrangement is significantly more frequent when the genes display strongly correlated expression levels or have one or more Gene Ontology (GO) classes in common. A correlation between intergenic distance in divergent and convergent gene pairs and shared GO classes was observed, in varying degrees, in rice and Populus but not in Arabidopsis. Furthermore, multiple GO classes were either overrepresented or underrepresented in Arabidopsis and Populus gene pairs, while only two GO classes were underrepresented in rice divergent gene pairs. Three cis-regulatory elements common to both Arabidopsis and rice were overrepresented in the intergenic regions of strongly correlated divergent gene pairs compared to those of noncorrelated pairs. Our results suggest that shared as well as unique mechanisms operate in shaping the organization and function of divergent and convergent gene pairs in different plant species.     

Implementing a rational and consistent nomenclature for serine/arginine-rich protein splicing factors (SR proteins) in plants

Growing interest in alternative splicing in plants and the extensive sequencing of new plant genomes necessitate more precise definition and classification of genes coding for splicing factors. SR proteins are a family of RNA binding proteins, which function as essential factors for constitutive and alternative splicing. We propose a unified nomenclature for plant SR proteins, taking into account the newly revised nomenclature of the mammalian SR proteins and a number of plant-specific properties of the plant proteins. We identify six subfamilies of SR proteins in Arabidopsis thaliana and rice (Oryza sativa), three of which are plant specific. The proposed subdivision of plant SR proteins into different subfamilies will allow grouping of paralogous proteins and simple assignment of newly discovered SR orthologs from other plant species and will promote functional comparisons in diverse plant species.     

Species-specific double-strand break repair and genome evolution in plants

Even closely related eukaryotic species may differ drastically in genome size. While insertion of retroelements represents a major source of genome enlargement, the mechanism mediating species- specific deletions is fairly obscure. We analyzed the formation of deletions during double-strand break (DSB) repair in Arabidopsis thaliana and tobacco, two dicotyledonous plant species differing >20-fold in genome size. DSBs were induced by the rare cutting restriction endonuclease I-SCE:I and deletions were identified by loss of function of a negative selectable marker gene containing an I-SCE:I site. Whereas the partial use of micro-homologies in junction formation was similar in both species, in tobacco 40% of the deletions were accompanied by insertions. No insertions could be detected in Arabidopsis , where larger deletions were more frequent, indicating a putative inverse correlation between genome size and the average length of deletions. Such a correlation has been postulated before by a theoretical study on the evolution of related insect genomes and our study now identifies a possible molecular cause for the phenomenon, indicating that species-specific differences in DSB repair might indeed influence genome evolution.     

EMF基因与植物营养生长

在高等植物中,外源和内源因素共同调控着植物从营养生长到生殖生长的转换。拟南芥EMF1和EMF2基因缺失的突变体不经过任何营养生长,种子萌发后便开花,这说明EMF基因是植物花发育的抑制基因。目前已从水稻、玉米、拟南芥等植物中克隆得到EMF同源基因,但其功能研究大多停留在拟南芥上。研究表明,EMF基因决定着植物营养生长阶段的发育,抑制植物开花。因此,开展EMF基因的分离、克隆和功能研究,有利于阐述植物营养生长过程阶段的抑花机制。对EMF基因的研究进展进行了综述,并提出EMF基因表达调控的闸门模型,以对EMF基因功能的进一步分析提供参考。