Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution

Independent ancient polyploidy events in the sister families Brassicaceae and Cleomaceae

Recent studies have elucidated the ancient polyploid history of the Arabidopsis thaliana (Brassicaceae) genome. The studies concur that there was at least one polyploidy event occurring some 14.5 to 86 million years ago (Mya), possibly near the divergence of the Brassicaceae from its sister family, Cleomaceae. Using a comparative genomics approach, we asked whether this polyploidy event was unique to members of the Brassicaceae, shared with the Cleomaceae, or an independent polyploidy event in each lineage. We isolated and sequenced three genomic regions from diploid Cleome spinosa (Cleomaceae) that are each homoeologous to a duplicated region shared between At3 and At5, centered on the paralogs of SEPALLATA (SEP) and CONSTANS (CO). Phylogenetic reconstructions and analysis of synonymous substitution rates support the hypothesis that a genomic triplication in Cleome occurred independently of and more recently than the duplication event in the Brassicaceae. There is a strong correlation in the copy number (single versus duplicate) of individual genes, suggesting functionally consistent influences operating on gene copy number in these two independently evolving lineages. However, the amount of gene loss in Cleome is greater than in Arabidopsis. The genome of C. spinosa is only 1.9 times the size of A. thaliana, enabling comparative genome analysis of separate but related polyploidy events.

Independent ancient polyploidy events in the sister families Brassicaceae and Cleomaceae

Recent studies have elucidated the ancient polyploid history of the Arabidopsis thaliana (Brassicaceae) genome. The studies concur that there was at least one polyploidy event occurring some 14.5 to 86 million years ago (Mya), possibly near the divergence of the Brassicaceae from its sister family, Cleomaceae. Using a comparative genomics approach, we asked whether this polyploidy event was unique to members of the Brassicaceae, shared with the Cleomaceae, or an independent polyploidy event in each lineage. We isolated and sequenced three genomic regions from diploid Cleome spinosa (Cleomaceae) that are each homoeologous to a duplicated region shared between At3 and At5, centered on the paralogs of SEPALLATA (SEP) and CONSTANS (CO). Phylogenetic reconstructions and analysis of synonymous substitution rates support the hypothesis that a genomic triplication in Cleome occurred independently of and more recently than the duplication event in the Brassicaceae. There is a strong correlation in the copy number (single versus duplicate) of individual genes, suggesting functionally consistent influences operating on gene copy number in these two independently evolving lineages. However, the amount of gene loss in Cleome is greater than in Arabidopsis. The genome of C. spinosa is only 1.9 times the size of A. thaliana, enabling comparative genome analysis of separate but related polyploidy events.

Multiple mechanisms promote the retained expression of gene duplicates in the tetraploid frog Xenopus laevis

Gene duplication provides a window of opportunity for biological variants to persist under the protection of a co-expressed copy with similar or redundant function. Duplication catalyzes innovation (neofunctionalization), subfunction degeneration (subfunctionalization), and genetic buffering (redundancy), and the genetic survival of each paralog is triggered by mechanisms that add, compromise, or do not alter protein function. We tested the applicability of three types of mechanisms for promoting the retained expression of duplicated genes in 290 expressed paralogs of the tetraploid clawed frog, Xenopus laevis. Tests were based on explicit expectations concerning the ka/ks ratio, and the number and location of nonsynonymous substitutions after duplication. Functional constraints on the majority of paralogs are not significantly different from a singleton ortholog. However, we recover strong support that some of them have an asymmetric rate of nonsynonymous substitution: 6% match predictions of the neofunctionalization hypothesis in that (1) each paralog accumulated nonsynonymous substitutions at a significantly different rate and (2) the one that evolves faster has a higher ka/ks ratio than the other paralog and than a singleton ortholog. Fewer paralogs (3%) exhibit a complementary pattern of substitution at the protein level that is predicted by enhancement or degradation of different functional domains, and the remaining 13% have a higher average ka/ks ratio in both paralogs that is consistent with altered functional constraints, diversifying selection, or activity-reducing mutations after duplication. We estimate that these paralogs have been retained since they originated by genome duplication between 21 and 41 million years ago. Multiple mechanisms operate to promote the retained expression of duplicates in the same genome, in genes in the same functional class, over the same period of time following duplication, and sometimes in the same pair of paralogs. None of these paralogs are superfluous; degradation or enhancement of different protein subfunctions and neofunctionalization are plausible hypotheses for the retained expression of some of them. Evolution of most X. laevis paralogs, however, is consistent with retained expression via mechanisms that do not radically alter functional constraints, such as selection to preserve post-duplication stoichiometry or temporal, quantitative, or spatial subfunctionalization.

Gene duplication plays a fundamental role in biological innovation but it is not clear how both copies of a duplicated gene manage to circumvent degradation by mutation if neither is unique. This study explores genetic mechanisms that could make each copy of a duplicate gene different, and therefore distinguishable and potentially preserved by natural selection. It is based on DNA sequences of the protein-coding region of 290 expressed duplicated genes in a frog, Xenopus laevis, that underwent complete duplication of its entire genome. Results provide evidence for multiple mechanisms acting within the same genome, within the same functional classes of genes, within the same period of time following duplication, and even on the same set of duplicated genes. Each copy of a duplicate gene may be subject to distinct evolutionary constraints, and this could be associated with degradation or enhancement of function. Functional constraints of most of these duplicates, however, are not substantially different from a single copy gene; their persistence in the first dozens of millions of years after duplication may more frequently be explained by mechanisms acting on their expression rather than their function.

The evolution and diversification of Dicers in plants

Most multicellular organisms regulate developmental transitions by microRNAs, which are generated by an enzyme, Dicer. Insects and fungi have two Dicer-like genes, and many animals have only one, yet the plant, Arabidopsis, has four. Examining the poplar and rice genomes revealed that they contain five and six Dicer-like genes, respectively. Analysis of these genes suggests that plants require a basic set of four Dicer types which were present before the divergence of mono- and dicotyledonous plants ( approximately 200 million years ago), but after the divergence of plants from green algae. A fifth type of Dicer seems to have evolved in monocots.

Partial shotgun sequencing of the Boechera stricta genome reveals extensive microsynteny and promoter conservation with Arabidopsis

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value < or = 10(-30)) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5' to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5' noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.

Evolution of Arabidopsis microRNA families through duplication events

Recently there has been a great interest in the identification of microRNAs and their targets as well as understanding the spatial and temporal regulation of microRNA genes. To understand how microRNA genes evolve, we looked at several rapidly evolving families in Arabidopsis thaliana, and found that they arose from a process of genome-wide duplication, tandem duplication, and segmental duplication followed by dispersal and diversification, similar to the processes that drive the evolution of protein gene families. Using multiple expression data sets to examine the transcription patterns of different members of the microRNA families, we find the sequence diversification of duplicated microRNA genes to be accompanied by a change in spatial and temporal expression patterns, suggesting that duplicated copies acquire new functionality as they evolve.

Ancient signals: comparative genomics of plant MAPK and MAPKK gene families

MAPK signal transduction modules play crucial roles in regulating many biological processes in plants, and their components are encoded by highly conserved genes. The recent availability of genome sequences for rice and poplar now makes it possible to examine how well the previously described Arabidopsis MAPK and MAPKK gene family structures represent the broader evolutionary situation in plants, and analysis of gene expression data for MPK and MKK genes in all three species allows further refinement of those families, based on functionality. The Arabidopsis MAPK nomenclature appears sufficiently robust to allow it to be usefully extended to other well-characterized plant systems.

Paleopolyploidy and gene duplication in soybean and other legumes

Two of the most important observations from whole-genome sequences have been the high rate of gene birth and death and the prevalence of large-scale duplication events, including polyploidy. There is also a growing appreciation that polyploidy is more than the sum of the gene duplications it creates, in part because polyploidy duplicates the members of entire regulatory networks. Thus, it may be important to distinguish paralogs that are produced by individual gene duplications from the homoeologous sequences produced by (allo)polyploidy. This is not a simple task, for several reasons, including the chromosomally cryptic nature of many duplications and the variable rates of gene evolution. Recent progress has been made in understanding patterns of gene and genome duplication in the legume family, specifically in soybean.

Transgressive segregation reveals two Arabidopsis TIR-NB-LRR resistance genes effective against Leptosphaeria maculans, causal agent of blackleg disease

In a cross between the two resistant accessions Col-0 and Ler-0, a 15:1 segregation was found in F2, suggesting the presence of unlinked resistance loci to Leptosphaeria maculans. One hundred Col-4 x Ler-0, and 50 Ler-2 x Cvi-1 recombinant inbred lines, and seven susceptible Ler-0 x Ws-0 F2 progenies were examined to identify the two loci. Resistance in Col-4, Ws-0 and Cvi-1 (RLM1) was mapped to the marker m305 on chromosome 1. Col-4 x Ler-0 and Ler-2 x Cvi-1 mapping populations located RLM2(Ler) on the same arm of chromosome 4. A tight physical location of RLM2 was established through near-isogenic lines. This region was found to correspond to an ancient duplication event between the RLM1 and RLM2 loci. Two independent T-DNA mutants in a TIR-NB-LRR R gene (At1g64070) displayed susceptibility, and L. maculans susceptible mutant phenotypes were confirmed to be allelic for rlm1 in F1 after crosses with susceptible rlm1(Ler)rlm2(Col) plants. Complementation of rlm1(Ler)rlm2(Col) with the genomic Col-0 sequence of At1g64070 conferred resistance. In addition, two T-DNA mutants in a neighbouring homologous TIR-NB-LRR gene (At1g63880) displayed moderate susceptibility to L. maculans. Sequence analysis revealed that At1g64070 was truncated by a premature stop codon, and that At1g63880 was absent in Ler-0. RNA interference confirmed that Ler-0 resistance is dependent on genes structurally related to RLM1. Camalexin was identified as a quantitative co-dominant resistance factor of Col-0 origin, but independent of RLM1. RLM1/RLM2 resistance was, however, found to require RAR1 and partially HSP90.1.

Promoter Analysis of MADS-Box Genes in Eudicots Through Phylogenetic Footprinting

The MIKC MADS-box gene family has been shaped by extensive gene duplications giving rise to subfamilies of genes with distinct functions and expression patterns. However, within these subfamilies the functional assignment is not that clear-cut, and considerable functional redundancy exists. One way to investigate the diversity in regulation present in these subfamilies is promoter sequence analysis. With the advent of genome sequencing projects, we are now able to exert a comparative analysis of Arabidopsis and poplar promoters of MADS-box genes belonging to the same subfamily. Based on the principle of phylogenetic footprinting, sequences conserved between the promoters of homologous genes are thought to be functional. Here, we have investigated the evolution of MADS-box genes at the promoter level and show that many genes have diverged in their regulatory sequences after duplication and/or speciation. Furthermore, using phylogenetic footprinting, a distinction can be made between redundancy, neo/nonfunctionalization, and subfunctionalization.

A simplified explanation for the frameshift mutation that created a novel C-terminal motif in the APETALA3 gene lineage

BACKGROUND

The evolution of type II MADS box genes has been extensively studied in angiosperms. One of the best-understood subfamilies is that of the Arabidopsis gene APETALA3 (AP3). Previous work has demonstrated that the ancestral paleoAP3 lineage was duplicated at some point within the basal eudicots to give rise to the paralogous TM6 and euAP3 lineages. This event was followed in euAP3 orthologs by the replacement of the C-terminal paleoAP3 motif with the derived euAP3 motif. It has been suggested that the new motif was created by an eight-nucleotide insertion that produced a translational frameshift.

RESULTS

The addition of 25 eudicot AP3 homologs to the existing dataset has allowed us to clarify the process by which the euAP3 motif evolved. Phylogenetic analysis indicates that the euAP3/TM6 duplication maps very close to the base of the core eudicots, associated with the families Trochodendraceae and Buxaceae. We demonstrate that although the transformation of paleoAP3 into euAP3 was due to a frameshift mutation, this was the result of a single nucleotide deletion. The use of ancestral character state reconstructions has allowed us to demonstrate that the frameshift was accompanied by few other nucleotide changes. We further confirm that the sequence is evolving as coding region.

CONCLUSION

This study demonstrates that the simplest of genetic changes can result in the remodeling of protein sequence to produce a kind of molecular 'hopeful monster.' Moreover, such a novel protein motif can become conserved almost immediately on the basis of what appears to be a rapidly generated new function. Given that the existing data on the function of such C-terminal motifs are somewhat disparate and contradictory, we have sought to synthesize previous findings within the context of the current analysis and thereby highlight specific hypotheses that require further investigation before the significance of the euAP3 frameshift event can be fully understood.

Evolution of Arabidopsis microRNA families through duplication events

Recently there has been a great interest in the identification of microRNAs and their targets as well as understanding the spatial and temporal regulation of microRNA genes. To understand how microRNA genes evolve, we looked at several rapidly evolving families in Arabidopsis thaliana, and found that they arose from a process of genome-wide duplication, tandem duplication, and segmental duplication followed by dispersal and diversification, similar to the processes that drive the evolution of protein gene families. Using multiple expression data sets to examine the transcription patterns of different members of the microRNA families, we find the sequence diversification of duplicated microRNA genes to be accompanied by a change in spatial and temporal expression patterns, suggesting that duplicated copies acquire new functionality as they evolve.

Diagnosing duplications – can it be done?

New genes arise through duplication and modification of DNA sequences on a range of scales: single gene duplication, duplication of large chromosomal fragments and whole-genome duplication. Each duplication mechanism has specific characteristics that influence the fate of the resulting duplicates, such as the size of the duplicated fragment, the potential for dosage imbalance, the preservation or disruption of regulatory control and genomic context. The ability to diagnose or identify the mechanism that produced a pair of paralogs has the potential to increase our ability to reconstruct evolutionary history, to understand the processes that govern genome evolution and to make functional predictions based on paralogy. The recent availability of large amounts of whole-genome sequence, often from several closely related species, has stimulated a wealth of new computational methods to diagnose gene duplications.

Nonrandom divergence of gene expression following gene and genome duplications in the flowering plant Arabidopsis thaliana

Analysis of expression data of duplicated genes in Arabidopsis thaliana shows that the mode of duplication, the time since duplication and the function of the duplicated genes play a role in the divergence of their expression.

Background

Genome analyses have revealed that gene duplication in plants is rampant. Furthermore, many of the duplicated genes seem to have been created through ancient genome-wide duplication events. Recently, we have shown that gene loss is strikingly different for large- and small-scale duplication events and highly biased towards the functional class to which a gene belongs. Here, we study the expression divergence of genes that were created during large- and small-scale gene duplication events by means of microarray data and investigate both the influence of the origin (mode of duplication) and the function of the duplicated genes on expression divergence.

Results

Duplicates that have been created by large-scale duplication events and that can still be found in duplicated segments have expression patterns that are more correlated than those that were created by small-scale duplications or those that no longer lie in duplicated segments. Moreover, the former tend to have highly redundant or overlapping expression patterns and are mostly expressed in the same tissues, while the latter show asymmetric divergence. In addition, a strong bias in divergence of gene expression was observed towards gene function and the biological process genes are involved in.

Conclusion

By using microarray expression data for Arabidopsis thaliana, we show that the mode of duplication, the function of the genes involved, and the time since duplication play important roles in the divergence of gene expression and, therefore, in the functional divergence of genes after duplication.

Molecular characterization and differential expression of cytokinin-responsive type-A response regulators in rice (Oryza sativa)

BACKGROUND

The response regulators represent the elements of bacterial two-component system and have been characterized from dicot plants like Arabidopsis but little information is available on the monocots, including the cereal crops. The aim of this study was to characterize type-A response regulator genes from rice, and to investigate their expression in various organs as well as in response to different hormones, including cytokinin, and environmental stimuli.

RESULTS

By analysis of the whole genome sequence of rice, we have identified ten genes encoding type-A response regulators based upon their high sequence identity within the receiver domain. The exon-intron organization, intron-phasing as well as chromosomal location of all the RT-PCR amplified rice (Oryza sativa) response regulator (OsRR) genes have been analyzed. The transcripts of OsRR genes could be detected by real-time PCR in all organs of the light- and dark-grown rice seedlings/plants, although there were quantitative differences. The steady-state transcript levels of most of the OsRR genes increased rapidly (within 15 min) on exogenous cytokinin application even in the presence of cycloheximide. Moreover, the expression of the OsRR6 gene was enhanced in rice seedlings exposed to salinity, dehydration and low temperature stress.

CONCLUSION

Ten type-A response regulator genes identified in rice, the model monocot plant, show overlapping/differential expression patterns in various organs and in response to light. The induction of OsRR genes by cytokinin even in the absence of de novo protein synthesis qualifies them to be primary cytokinin response genes. The induction of OsRR6 in response to different environmental stimuli indicates its role in cross-talk between abiotic stress and cytokinin signaling. These results provide a foundation for further investigations on specific as well as overlapping cellular functions of type-A response regulators in rice.

Buffering of crucial functions by paleologous duplicated genes may contribute cyclicality to angiosperm genome duplication

Genome duplication followed by massive gene loss has permanently shaped the genomes of many higher eukaryotes, particularly angiosperms. It has long been believed that a primary advantage of genome duplication is the opportunity for the evolution of genes with new functions by modification of duplicated genes. If so, then patterns of genetic diversity among strains within taxa might reveal footprints of selection that are consistent with this advantage. Contrary to classical predictions that duplicated genes may be relatively free to acquire unique functionality, we find among both Arabidopsis ecotypes and Oryza subspecies that SNPs encode less radical amino acid changes in genes for which there exists a duplicated copy at a "paleologous" locus than in "singleton" genes. Preferential retention of duplicated genes encoding long complex proteins and their unexpectedly slow divergence (perhaps because of homogenization) suggest that a primary advantage of retaining duplicated paleologs may be the buffering of crucial functions. Functional buffering and functional divergence may represent extremes in the spectrum of duplicated gene fates. Functional buffering may be especially important during "genomic turmoil" immediately after genome duplication but continues to act approximately 60 million years later, and its gradual deterioration may contribute cyclicality to genome duplication in some lineages.

Genes involved in biosynthesis and signalisation of ethylene in Brassica oleracea and Arabidopsis thaliana: identification and genome comparative mapping of specific gene homologues

The study reported was aimed at the identification and determination of the chromosomal organisation of genes involved in the ethylene biosynthesis and signalling pathways in Brassica oleracea, on the basis of the Arabidopsis thaliana DNA probes and in silico genome analysis. Because of its polyploidal origin, the B. oleracea genome is characterised by extensive gene redundancy. Therefore, an important aspect of gene expression in B. oleracea response to environmental stimuli is to identify the specific gene copy involved. This aspect should also be taken into consideration while studying the genetic basis of biosynthesis and signal transduction in relation to basic phytohormones. Our present work concerns the identification of homologue genes involved in ethylene biosynthesis such as SAM, ACS and ACO, as well as those involved in the ethylene signalling pathway, mainly ETR1, CTR1, MKK4, MKK5, EIN2, EIN3, EREBP, ERF5 and ERF7 on the basis of the restriction fragment length polymorphism (RFLP) and PCR mapping. In the case of ACC synthases, (ACSs) the in silico analysis of gene variants in the genome of A. thaliana was followed by the identification of homologues to ACS2, ACS6 and ACS7 in the B. oleracea database. In total, 22 loci with sequence homology to the genes under analysis were included in the existing B. oleracea RFLP chromosomal map. Based on the stress responsiveness of most of the A. thaliana genes analysed in this study, we performed initial functional analysis of some gene homologues mapped. With the use of the RT-PCR approach the conservation of differential transcriptional induction of ACS homologues in the B. oleracea and A. thaliana was demonstrated during ozone stress.

Genome-wide analysis of the ERF gene family in Arabidopsis and rice

Genes in the ERF family encode transcriptional regulators with a variety of functions involved in the developmental and physiological processes in plants. In this study, a comprehensive computational analysis identified 122 and 139 ERF family genes in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa L. subsp. japonica), respectively. A complete overview of this gene family in Arabidopsis is presented, including the gene structures, phylogeny, chromosome locations, and conserved motifs. In addition, a comparative analysis between these genes in Arabidopsis and rice was performed. As a result of these analyses, the ERF families in Arabidopsis and rice were divided into 12 and 15 groups, respectively, and several of these groups were further divided into subgroups. Based on the observation that 11 of these groups were present in both Arabidopsis and rice, it was concluded that the major functional diversification within the ERF family predated the monocot/dicot divergence. In contrast, some groups/subgroups are species specific. We discuss the relationship between the structure and function of the ERF family proteins based on these results and published information. It was further concluded that the expansion of the ERF family in plants might have been due to chromosomal/segmental duplication and tandem duplication, as well as more ancient transposition and homing. These results will be useful for future functional analyses of the ERF family genes.

Developmental and gene expression analyses of a cotton naked seed mutant

Cotton fiber development is a fundamental biological phenomenon, yet the molecular basis of fiber cell initiation is poorly understood. We examined molecular and cellular events of fiber cell development in the naked seed mutant (N1N1) and its isogenic line of cotton (Gossypium hirsutum L. cv. Texas Marker-1, TM-1). The dominant mutation not only delayed the process of fiber cell formation and elongation but also reduced the total number of fiber cells, resulting in sparsely distributed short fibers. Gene expression changes in TM-1 and N1N1 mutant lines among four tissues were analyzed using spotted cotton oligo-gene microarrays. Using the Arabidopsis genes, we selected and designed approximately 1,334 70-mer oligos from a subset of cotton fiber ESTs. Statistical analysis of the microarray data indicates that the number of significantly differentially expressed genes was 856 in the leaves compared to the ovules (3 days post-anthesis, DPA), 632 in the petals relative to the ovules (3 DPA), and 91 in the ovules at 0 DPA compared to 3 DPA, all in TM-1. Moreover, 117 and 30 genes were expressed significantly different in the ovules at three and 0 DPA, respectively, between TM-1 and N1N1. Quantitative RT-PCR analysis of 23 fiber-associated genes in seven tissues including ovules, fiber-bearing ovules, fibers, and non-fiber tissues in TM-1 and N1N1 indicates a mode of temporal regulation of the genes involved in transcriptional and translational regulation, signal transduction, and cell differentiation during early stages of fiber development. Suppression of the fiber-associated genes in the mutant may suggest that the N1N1 mutation disrupts temporal regulation of gene expression, leading to a defective process of fiber cell elongation and development.

Genome-wide analysis of the ERF gene family in Arabidopsis and rice

Genes in the ERF family encode transcriptional regulators with a variety of functions involved in the developmental and physiological processes in plants. In this study, a comprehensive computational analysis identified 122 and 139 ERF family genes in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa L. subsp. japonica), respectively. A complete overview of this gene family in Arabidopsis is presented, including the gene structures, phylogeny, chromosome locations, and conserved motifs. In addition, a comparative analysis between these genes in Arabidopsis and rice was performed. As a result of these analyses, the ERF families in Arabidopsis and rice were divided into 12 and 15 groups, respectively, and several of these groups were further divided into subgroups. Based on the observation that 11 of these groups were present in both Arabidopsis and rice, it was concluded that the major functional diversification within the ERF family predated the monocot/dicot divergence. In contrast, some groups/subgroups are species specific. We discuss the relationship between the structure and function of the ERF family proteins based on these results and published information. It was further concluded that the expansion of the ERF family in plants might have been due to chromosomal/segmental duplication and tandem duplication, as well as more ancient transposition and homing. These results will be useful for future functional analyses of the ERF family genes.

Genetic and epigenetic consequences of recent hybridization and polyploidy in Spartina (Poaceae)

To study the consequences of hybridization and genome duplication on polyploid genome evolution and adaptation, we used independently formed hybrids (Spartina x townsendii and Spartina x neyrautii) that originated from natural crosses between Spartina alterniflora, an American introduced species, and the European native Spartina maritima. The hybrid from England, S. x townsendii, gave rise to the invasive allopolyploid, salt-marsh species, Spartina anglica. Recent studies indicated that allopolyploid speciation may be associated with rapid genetic and epigenetic changes. To assess this in Spartina, we performed AFLP (amplified fragment length polymorphism) and MSAP (methylation sensitive amplification polymorphism) on young hybrids and the allopolyploid. By comparing the subgenomes in the hybrids and the allopolyploid to the parental species, we inferred structural changes that arose repeatedly in the two independently formed hybrids. Surprisingly, 30% of the parental methylation patterns are altered in the hybrids and the allopolyploid. This high level of epigenetic regulation might explain the morphological plasticity of Spartina anglica and its larger ecological amplitude. Hybridization rather than genome doubling seems to have triggered most of the methylation changes observed in Spartina anglica.

Varying signals of the effects of natural selection during teleost growth hormone gene evolution

The growth hormone (GH) gene of teleost fish exhibits a higher degree of variability compared with other vertebrate groups. However, the different selective constraints at the sequence level are not well understood. In this study, maximum-likelihood (ML) models of codon substitutions were used to investigate Darwinian adaptive evolution of the GH gene in teleost fishes. Complete GH gene sequences of 54 fish species were classified into 4 orders, and the variable nature of GH was examined by determining the dNand dSrate variation and the rates of molecular evolution for each teleost order. The results indicate that although the overall evolution rate for teleost GH is high ((1.15 ± 0.01) × 10–9substitutions/(aa site·y)) compared with the "slow phases" in mammals ((0.21 to 0.28 ± 0.05) × 10–9), the vital structure of this gene has been retained. While the majority of the amino acid changes appear to be due to relaxation of purifying selection, some positively selected sites were detected in regions with no specifically identified role in protein function. The positively selected regions observed in salmoniformes lineage suggests a possible role for positive selection driving functional divergence in paralogous forms of the GH gene after whole-genome duplication in this lineage.Key words: teleost fish, growth hormone, positive selection, synonymous substitution, non-synonymous substitution, molecular evolution.

Genome histories clarify evolution of the expansin superfamily: new insights from the poplar genome and pine ESTs

Expansins comprise a superfamily of plant cell wall-loosening proteins that has been divided into four distinct families, EXPA, EXPB, EXLA and EXLB. In a recent analysis of Arabidopsis thaliana and Oryza sativa expansins, we proposed a further subdivision of the families into 17 clades, representing independent lineages in the last common ancestor of monocots and eudicots. This division was based on both traditional sequence-based phylogenetic trees and on position-based trees, in which genomic locations and dated segmental duplications were used to reconstruct gene phylogeny. In this article we review recent work concerning the patterns of expansin evolution in angiosperms and include additional insights gained from the genome of a second eudicot species, Populus trichocarpa, which includes at least 36 expansin genes. All of the previously proposed monocot-eudicot orthologous groups, but no additional ones, are represented in this species. The results also confirm that all of these clades are truly independent lineages. Furthermore, we have used position-based phylogeny to clarify the history of clades EXPA-II and EXPA-IV. Most of the growth of the expansin superfamily in the poplar lineage is likely due to a recent polyploidy event. Finally, some monocot-eudicot clades are shown to have diverged before the separation of the angiosperm and gymnosperm lineages.

The photorespiratory Arabidopsis shm1 mutant is deficient in SHM1

Mitochondrial serine hydroxymethyltransferase (SHMT), combined with glycine decarboxylase, catalyzes an essential sequence of the photorespiratory C2 cycle, namely, the conversion of two molecules of glycine into one molecule each of CO2, NH4+, and serine. The Arabidopsis (Arabidopsis thaliana) mutant shm (now designated shm1-1) is defective in mitochondrial SHMT activity and displays a lethal photorespiratory phenotype when grown at ambient CO2, but is virtually unaffected at elevated CO2. The Arabidopsis genome harbors seven putative SHM genes, two of which (SHM1 and SHM2) feature predicted mitochondrial targeting signals. We have mapped shm1-1 to the position of the SHM1 gene (At4g37930). The mutation is due to a G --> A transition at the 5' splice site of intron 6 of SHM1, causing aberrant splicing and a premature termination of translation. A T-DNA insertion allele of SHM1, shm1-2, and the F1 progeny of a genetic cross between shm1-1 and shm1-2 displayed the same conditional lethal phenotype as shm1-1. Expression of wild-type SHM1 under the control of either the cauliflower mosaic virus 35S or the SHM1 promoter in shm1-1 abrogated the photorespiratory phenotype of the shm mutant, whereas overexpression of SHM2 or expression of SHM1 under the control of the SHM2 promoter did not rescue the mutant phenotype. Promoter-beta-glucuronidase analyses revealed that SHM1 is predominantly expressed in leaves, whereas SHM2 is mainly transcribed in the shoot apical meristem and roots. Our findings establish SHM1 as the defective gene in the Arabidopsis shm1-1 mutant.

Global studies of cell type-specific gene expression in plants

Technological advances in expression profiling and in the ability to collect minute quantities of tissues have come together to allow a growing number of global transcriptional studies at the cell level in plants. Microarray technology, with a choice of cDNA or oligo-based slides, is now well established, with commercial full-genome platforms for rice and Arabidopsis and extensive expressed sequence tag (EST)-based designs for many other species. Microdissection and cell sorting are two established methodologies that have been used in conjunction with microarrays to provide an early glimpse of the transcriptional landscape at the level of individual cell types. The results indicate that much of the transcriptome is compartmentalized. A minor but consistent percentage of transcripts appear to be unique to specific cell types. Functional analyses of cell-specific patterns of gene expression are providing important clues to cell-specific functions. The spatial dissection of the transcriptome has also yielded insights into the localized mediators of hormone inputs and promises to provide detail on cell-specific effects of microRNAs.

Genomewide nonadditive gene regulation in Arabidopsis allotetraploids

Polyploidy has occurred throughout the evolutionary history of all eukaryotes and is extremely common in plants. Reunification of the evolutionarily divergent genomes in allopolyploids creates regulatory incompatibilities that must be reconciled. Here we report genomewide gene expression analysis of Arabidopsis synthetic allotetraploids, using spotted 70-mer oligo-gene microarrays. We detected >15% transcriptome divergence between the progenitors, and 2105 and 1818 genes were highly expressed in Arabidopsis thaliana and A. arenosa, respectively. Approximately 5.2% (1362) and 5.6% (1469) genes displayed expression divergence from the midparent value (MPV) in two independently derived synthetic allotetraploids, suggesting nonadditive gene regulation following interspecific hybridization. Remarkably, the majority of nonadditively expressed genes in the allotetraploids also display expression changes between the parents, indicating that transcriptome divergence is reconciled during allopolyploid formation. Moreover, >65% of the nonadditively expressed genes in the allotetraploids are repressed, and >94% of the repressed genes in the allotetraploids match the genes that are expressed at higher levels in A. thaliana than in A. arenosa, consistent with the silencing of A. thaliana rRNA genes subjected to nucleolar dominance and with overall suppression of the A. thaliana phenotype in the synthetic allotetraploids and natural A. suecica. The nonadditive gene regulation is involved in various biological pathways, and the changes in gene expression are developmentally regulated. In contrast to the small effects of genome doubling on gene regulation in autotetraploids, the combination of two divergent genomes in allotetraploids by interspecific hybridization induces genomewide nonadditive gene regulation, providing a molecular basis for de novo variation and allopolyploid evolution.

Genome-wide comparative analysis of the IQD gene families in Arabidopsis thaliana and Oryza sativa

BACKGROUND

Calcium signaling plays a prominent role in plants for coordinating a wide range of developmental processes and responses to environmental cues. Stimulus-specific generation of intracellular calcium transients, decoding of calcium signatures, and transformation of the signal into cellular responses are integral modules of the transduction process. Several hundred proteins with functions in calcium signaling circuits have been identified, and the number of downstream targets of calcium sensors is expected to increase. We previously identified a novel, calmodulin-binding nuclear protein, IQD1, which stimulates glucosinolate accumulation and plant defense in Arabidopsis thaliana. Here, we present a comparative genome-wide analysis of a new class of putative calmodulin target proteins in Arabidopsis and rice.

RESULTS

We identified and analyzed 33 and 29 IQD1-like genes in Arabidopsis thaliana and Oryza sativa, respectively. The encoded IQD proteins contain a plant-specific domain of 67 conserved amino acid residues, referred to as the IQ67 domain, which is characterized by a unique and repetitive arrangement of three different calmodulin recruitment motifs, known as the IQ, 1-5-10, and 1-8-14 motifs. We demonstrated calmodulin binding for IQD20, the smallest IQD protein in Arabidopsis, which consists of a C-terminal IQ67 domain and a short N-terminal extension. A striking feature of IQD proteins is the high isoelectric point (approximately 10.3) and frequency of serine residues (approximately 11%). We compared the Arabidopsis and rice IQD gene families in terms of gene structure, chromosome location, predicted protein properties and motifs, phylogenetic relationships, and evolutionary history. The existence of an IQD-like gene in bryophytes suggests that IQD proteins are an ancient family of calmodulin-binding proteins and arose during the early evolution of land plants.

CONCLUSION

Comparative phylogenetic analyses indicate that the major IQD gene lineages originated before the monocot-eudicot divergence. The extant IQD loci in Arabidopsis primarily resulted from segmental duplication and reflect preferential retention of paralogous genes, which is characteristic for proteins with regulatory functions. Interaction of IQD1 and IQD20 with calmodulin and the presence of predicted calmodulin binding sites in all IQD family members suggest that IQD proteins are a new class of calmodulin targets. The basic isoelectric point of IQD proteins and their frequently predicted nuclear localization suggest that IQD proteins link calcium signaling pathways to the regulation of gene expression. Our comparative genomics analysis of IQD genes and encoded proteins in two model plant species provides the first step towards the functional dissection of this emerging family of putative calmodulin targets.

Cross-species amplification of microsatellite loci within the dioecious, polyploid genus Actinidia (Actinidiaceae)

Microsatellite marker transfer across species in the dioecious genus Actinidia (kiwifruit) could offer an efficient and time-effective technique for use during trait transfer for vine and fruit improvement in breeding programmes. We evaluated the cross-species amplification of 20 EST-derived microsatellite markers that were fully informative in an Actinidia chinensis mapping family. We tested all 20 markers on 120 genotypes belonging to 21 species, 5 with varieties and/or chromosome races. These 26 taxa included 16 diploids, 7 tetraploids, 2 hexaploids and 1 octaploid, and represented all four taxonomic sections in the genus. All 20 markers showed some level of cross-species amplification. The most successful marker amplified in all genotypes from all species from all sections of the genus, the least successful amplified fragments only in A. chinensis and A. deliciosa. One species, A. glaucophylla, failed to amplify with all but 2 markers. PIC (Polymorphism information content) values were high, with 14 of 17 markers recording values of 0.90 and above. Sequence data demonstrated the presence of the microsatellite in all the amplified products. Sequence homology was less 5' of the microsatellite and increased toward the start codon of the translated region of the EST from which the marker was derived. The data confirm that EST-derived microsatellite markers from Actinidia species show cross-species amplification with high levels of polymorphism which could make them useful markers in breeding programmes.

Thioredoxins in Arabidopsis and other plants

Regulation of disulfide dithiol exchange has become increasingly important in our knowledge of plant life. Initially discovered as regulators of light-dependent malate biosynthesis in the chloroplast, plant thioredoxins are now implicated in a large panel of reactions related to metabolism, defense and development. In this review we describe the numerous thioredoxin types encoded by the Arabidopsis genome, and provide evidence that they are present in all higher plants. Some results suggest cross-talk between thioredoxins and glutaredoxins, the second family of disulfide dithiol reductase. The development of proteomics in plants revealed an unexpectedly large number of putative target proteins for thioredoxins and glutaredoxins. Nevertheless, we are far from a clear understanding of the actual function of each thioredoxin in planta. Although hampered by functional redundancies between genes, genetic approaches are probably unavoidable to define which thioredoxin interacts with which target protein and evaluate the physiological consequences.

Functional analysis of four members of the PsbP family in photosystem II in Nicotiana tabacum using differential RNA interference

Gene redundancy is frequently found in higher plants and complicates genetic analysis. In this study, a method referred to as 'differential RNA interference (dRNAi)' was used to investigate the psbP gene family in Nicotiana tabacum. PsbP is a membrane-extrinsic subunit of PSII and plays important roles in the water splitting reaction. N. tabacum has four psbP isogenes and the function of each isogene has not yet been characterized in vivo. To obtain transgenic tobacco plants with various amounts and compositions of PsbP members, the psbP isogenes were differentially silenced by RNA interference (RNAi) using the 3'-untranslated region (UTR) as a silencing trigger (dRNAi). In addition, the extra psbP genes without the 3'-UTR were complementarily transformed into the above silenced plants, which accumulated PsbP originating from the exogenous gene while differential silencing of the endogenous target was maintained. By using dRNAi and subsequent complementation (substitution) in dRNAi, we clearly demonstrated that, regardless of the of PsbP members that were accumulated, PSII activity was linearly correlated with the total amount of PsbP. Therefore, we concluded that the protein functions of the PsbP members in N. tabacum are equivalent in vivo, whereas full expression of the four isogenes is required for optimum PSII activity. These results demonstrate that the use of dRNAi and subsequent complementation/substitution in dRNAi would provide a new experimental approach for studying the function of multigene families in plants.

The Arabidopsis genome: A foundation for plant research

The sequence of the first plant genome was completed and published at the end of 2000. This spawned a series of large-scale projects aimed at discovering the functions of the 25,000+ genes identified in Arabidopsis thaliana (Arabidopsis). This review summarizes progress made in the past five years and speculates about future developments in Arabidopsis research and its implications for crop science. The provision of large populations of gene disruption lines to the research community has greatly accelerated the impact of genomics on many areas of plant science. The tools and community organization required for plant integrative and systems biology approaches are now ready to accomplish the next big step in plant biology--the integration of knowledge and modeling of biological processes. In the future, plant science will continue to be enriched by the alignment of high-quality basic research (generally conducted in Arabidopsis), with strategic objectives in crop plants. The sequence and analysis of an increasing number of crop plant genomes enhance this alignment and provide new insights into genome evolution and crop plant domestication.

Expression Pattern Shifts Following Duplication Indicative of Subfunctionalization and Neofunctionalization in Regulatory Genes of Arabidopsis

Gene duplication plays an important role in the evolution of diversity and novel function and is especially prevalent in the nuclear genomes of flowering plants. Duplicate genes may be maintained through subfunctionalization and neofunctionalization at the level of expression or coding sequence. In order to test the hypothesis that duplicated regulatory genes will be differentially expressed in a specific manner indicative of regulatory subfunctionalization and/or neofunctionalization, we examined expression pattern shifts in duplicated regulatory genes in Arabidopsis. A two-way analysis of variance was performed on expression data for 280 phylogenetically identified paralogous pairs. Expression data were extracted from global expression profiles for wild-type root, stem, leaf, developing inflorescence, nearly mature flower buds, and seedpod. Gene, organ, and gene by organ interaction (G x O) effects were examined. Results indicate that 85% of the paralogous pairs exhibited a significant G x O effect indicative of regulatory subfunctionalization and/or neofunctionalization. A significant G x O effect was associated with complementary expression patterns in 45% of pairwise comparisons. No association was detected between a G x O effect and a relaxed evolutionary constraint as detected by the ratio of nonsynonymous to synonymous substitutions. Ancestral gene expression patterns inferred across a Type II MADS-box gene phylogeny suggest several cases of regulatory neofunctionalization and organ-specific nonfunctionalization. Complete linkage clustering of gene expression levels across organs suggests that regulatory modules for each organ are independent or ancestral genes had limited expression. We propose a new classification, regulatory hypofunctionalization, for an overall decrease in expression level in one member of a paralogous pair while still having a significant G x O effect. We conclude that expression divergence specifically indicative of subfunctionalization and/or neofunctionalization contributes to the maintenance of most if not all duplicated regulatory genes in Arabidopsis and hypothesize that this results in increasing expression diversity or specificity of regulatory genes after each round of duplication.

Use of genomic history to improve phylogeny and understanding of births and deaths in a gene family

Polyploidy events have played an important role in the evolution of angiosperm genomes. Here, we demonstrate how genomic histories can increase phylogenetic resolution in a gene family, specifically the expansin superfamily of cell wall proteins. There are 36 expansins in Arabidopsis and 58 in rice. Traditional sequence-based phylogenetic trees yield poor resolution below the family level. To improve upon these analyses, we searched for gene colinearity (microsynteny) between Arabidopsis and rice genomic segments containing expansin genes. Multiple rounds of genome duplication and extensive gene loss have obscured synteny. However, by simultaneously aligning groups of up to 10 potentially orthologous segments from the two species, we traced the history of 49 out of 63 expansin-containing segments back to the ancestor of monocots and eudicots. Our results indicate that this ancestor had 15-17 expansin genes, each ancestral to an extant clade. Some clades have strikingly different growth patterns in the rice and Arabidopsis lineages, with more than half of all rice expansins arising from two ancestral genes. Segmental duplications, most of them part of polyploidy events, account for 12 out of 21 new expansin genes in Arabidopsis and 16 out of 44 in rice. Tandem duplications explain most of the rest. We were also able to estimate a minimum of 28 gene deaths in the Arabidopsis lineage and nine in rice. This analysis greatly clarifies expansin evolution since the last common ancestor of monocots and eudicots and the method should be broadly applicable to many other gene families.

Genetic mapping of expressed sequences in onion and in silico comparisons with rice show scant colinearity

The Poales (which include the grasses) and Asparagales [which include onion (Allium cepa L.) and other Allium species] are the two most economically important monocot orders. Enormous genomic resources have been developed for the grasses; however, their applicability to other major monocot groups, such as the Asparagales, is unclear. Expressed sequence tags (ESTs) from onion that showed significant similarities (80% similarity over at least 70% of the sequence) to single positions in the rice genome were selected. One hundred new genetic markers developed from these ESTs were added to the intraspecific map derived from the BYG15-23xAC43 segregating family, producing 14 linkage groups encompassing 1,907 cM at LOD 4. Onion linkage groups were assigned to chromosomes using alien addition lines of Allium fistulosum L. carrying single onion chromosomes. Visual comparisons of genetic linkage in onion with physical linkage in rice revealed scant colinearity; however, short regions of colinearity could be identified. Our results demonstrate that the grasses may not be appropriate genomic models for other major monocot groups such as the Asparagales; this will make it necessary to develop genomic resources for these important plants.

Organization and evolution of resistance gene analogs in peanut

The scarcity of genetic polymorphism in Arachis hypogaea (peanut), as in other monophyletic polyploid species, makes it especially vulnerable to nematode, bacterial, fungal, and viral pathogens. Although no disease resistance genes have been cloned from peanut itself, the conserved motifs in cloned resistance genes from other plant species provide a means to isolate and analyze similar genes from peanut. To survey the number, diversity, evolutionary history, and genomic organization of resistance gene-like sequences in peanut, we isolated 234 resistance gene analogs (RGAs) by using primers designed from conserved regions of different classes of resistance genes including NBS-LRR, and LRR-TM classes. Phylogenetic and sequence analyses were performed to explore evolutionary relationships both among peanut RGAs and with orthologous genes from other plant taxa. Fifty-six overgos designed from the RGA sequences on the basis of their phyletic association were applied to a peanut BAC library; 736 hybridizing BAC clones were fingerprinted and contigs were formed in order to gain insights into the genomic organization of these genes. All the fingerprinting gels were blotted and screened with the respective overgos in order to verify the authenticity of the hits from initial screens, and to explore the physical organization of these genes in terms of both copy number and distribution in the genome. As a result, we identified 250 putative resistance gene loci. A correlation was found between the phyletic positions of the sequences and their physical locations. The BACs isolated here will serve as a valuable resource for future applications, such as map-based cloning, and will help improve our understanding of the evolution and organization of these genes in the peanut genome.

Epigenetics and plant evolution

A fundamental precept of evolutionary biology is that natural selection acts on phenotypes determined by DNA sequence variation within natural populations. Recent advances in our understanding of gene regulation, however, have elucidated a spectrum of epigenetic molecular phenomena capable of altering the temporal, spatial, and abundance patterns of gene expression. These modifications may have morphological, physiological, and ecological consequences, and are heritable across generations, suggesting they are important in evolution. A corollary is that genetic variation per se is not always a prerequisite to evolutionary change. Here, we provide an introduction to epigenetic mechanisms in plants, and highlight some of the empirical studies illustrative of the possible connections between evolution and epigenetically mediated alterations in gene expression and morphology.

Structure and expression analysis of early auxin-responsive Aux/IAA gene family in rice (Oryza sativa)

Auxin exerts pleiotropic effects on plant growth and development by regulating the expression of early auxin-responsive genes of auxin/indoleacetic acid (Aux/IAA), small auxin-up RNA, and GH3 classes. These genes have been studied extensively in dicots like soybean and Arabidopsis. We had earlier characterized a cDNA of the first monocot member of Aux/IAA family from rice. The achievement of the large scale rice genome sequencing combined with the availability of full-length cDNA sequences from Knowledge-based Oryza Molecular Biological Encyclopedia provided us the opportunity to draw up the first comprehensive list of Aux/IAA genes in a monocot. By screening the available databases, we have identified 31 Aux/IAA genes having high sequence identity within the conserved domains I, II, III, and IV. The genomic organization as well as chromosomal location of all the Oryza sativa indoleacetic acid (OsIAA) genes is reported. The rice Aux/IAA proteins can be classified in two groups (A and B) on the basis of their phylogenetic relationship with Arabidopsis Aux/IAA proteins. An evolutionary pattern of the rice Aux/IAA genes has been discussed by analyzing their structure (exon/intron organization) and duplications. Interestingly, the duplication of rice Aux/IAA genes was found to be associated with chromosomal block duplication events in rice. The in-silico analysis has been complemented with real-time polymerase chain reaction analysis to quantify transcript levels of all Aux/IAA family members. OsIAA genes showed differential and overlapping organ-specific expression patterns in light- and dark-grown seedlings/plants. Although auxin enhanced the transcript abundance of most of the OsIAA genes, the effect was more pronounced on OsIAA9, 14, 19, 20, 24, and 31. These results provide a foundation for future studies on elucidating the precise role of rice Aux/IAA genes in early steps of auxin signal transduction.

Gene order in plants: a slow but sure shuffle

Comparative mapping studies have revealed a great deal about the patterns of gene order and gene content evolution in plants. These findings have practical importance for leveraging genomic information from model to nonmodel plant species. However, there is much to be learned about the processes by which gene order and content evolve. The role of gene duplication and loss in the evolution of plant gene order, in particular, appears to be more important than commonly appreciated. An exciting area of current research is the study of gene order and content polymorphism within species. Some recent findings suggest that there may be a functional, and adaptive, relationship between gene order and phenotype that is mediated by the effects of gene order on transcriptional regulation.

Flower development and evolution: gene duplication, diversification and redeployment

Gene duplication and diversification can provide the raw material for the evolution of new morphologies. In plants, the numbers of MADS-box genes have multiplied considerably, resulting in a plethora of these transcriptional regulators in the angiosperms (flowering plants). MADS-box genes have been implicated in the regulation of a variety of flower developmental processes; therefore, understanding the functional consequences of duplication and diversification in this gene family can shed light on the evolution of different floral forms. Recent functional analyses of MADS-box gene lineages have demonstrated that in various instances these genes have swapped roles, acquired novel roles, or retained ancestral roles. These studies underscore the idea that gene function cannot be extrapolated from structural orthology.

Expression divergence between duplicate genes

A general picture of the role of expression divergence in the evolution of duplicate genes is emerging, thanks to the availability of completely sequenced genomes and functional genomic data, such as microarray data. It is now clear that expression divergence, regulatory-motif divergence and coding-sequence divergence all increase with the age of duplicate genes, although their exact interrelationships remain to be determined. It is also clear that gene duplication increases expression diversity and enables tissue or developmental specialization to evolve. However, the relative roles of subfunctionalization and neofunctionalization in the retention of duplicate genes remain to be clarified, especially for higher eukaryotes. In addition, the relationship between gene duplication and evolution of transcriptional regulatory networks is largely unexplored.

Molecular evidence for asymmetric evolution of sister duplicated blocks after cereal polyploidy

Polyploidy (genome duplication) is thought to have contributed to the evolution of the eukaryotic genome, but complex genome structures and massive gene loss during evolution has complicated detection of these ancestral duplication events. The major factors determining the fate of duplicated genes are currently unclear, as are the processes by which duplicated genes evolve after polyploidy. Fine-scale analysis between homologous regions may allow us to better understand post-polyploidy evolution. Here, using gene-by-gene and gene-by-genome strategies, we identified the S5 region and four homologous regions within the japonica genome. Additional phylogenomic analyses of the comparable duplicated blocks indicate that four successive duplication events gave rise to these five regions, allowing us to propose a model for this local chromosomal evolution. According to this model, gene loss may play a major role in post-duplication genetic evolution at the segmental level. Moreover, we found molecular evidence that one of the sister duplicated blocks experienced more gene loss and a more rapid evolution subsequent to two recent duplication events. Given that these two recent duplication events were likely involved in polyploidy, this asymmetric evolution (gene loss and gene divergence) may be one possible mechanism accounting for the diploidization at the segmental level.

Homeodomain leucine zipper class I genes in Arabidopsis. Expression patterns and phylogenetic relationships

Members of the homeodomain leucine zipper (HDZip) family of transcription factors are present in a wide range of plants, from mosses to higher plants, but not in other eukaryotes. The HDZip genes act in developmental processes, including vascular tissue and trichome development, and several of them have been suggested to be involved in the mediation of external signals to regulate plant growth. The Arabidopsis (Arabidopsis thaliana) genome contains 47 HDZip genes, which, based on sequence criteria, have been grouped into four different classes: HDZip I to IV. In this article, we present an overview of the class I HDZip genes in Arabidopsis. We describe their expression patterns, transcriptional regulation properties, duplication history, and phylogeny. The phylogeny of HDZip class I genes is supported by data on the duplication history of the genes, as well as the intron/exon patterning of the HDZip-encoding motifs. The HDZip class I genes were found to be widely expressed and partly to have overlapping expression patterns at the organ level. Further, abscisic acid or water deficit treatments and different light conditions affected the transcript levels of a majority of the HDZip I genes. Within the gene family, our data show examples of closely related HDZip genes with similarities in the function of the gene product, but a divergence in expression pattern. In addition, six HDZip class I proteins tested were found to be activators of gene expression. In conclusion, several HDZip I genes appear to regulate similar cellular processes, although in different organs or tissues and in response to different environmental signals.

Transcription factor families have much higher expansion rates in plants than in animals

Transcription factors (TFs), which are central to the regulation of gene expression, are usually members of multigene families. In plants, they are involved in diverse processes such as developmental control and elicitation of defense and stress responses. To investigate if differences exist in the expansion patterns of TF gene families between plants and other eukaryotes, we first used Arabidopsis (Arabidopsis thaliana) TFs to identify TF DNA-binding domains. These DNA-binding domains were then used to identify related sequences in 25 other eukaryotic genomes. Interestingly, among 19 families that are shared between animals and plants, more than 14 are larger in plants than in animals. After examining the lineage-specific expansion of TF families in two plants, eight animals, and two fungi, we found that TF families shared among these organisms have undergone much more dramatic expansion in plants than in other eukaryotes. Moreover, this elevated expansion rate of plant TF is not simply due to higher duplication rates of plant genomes but also to a higher degree of expansion compared to other plant genes. Further, in many Arabidopsis-rice (Oryza sativa) TF orthologous groups, the degree of lineage-specific expansion in Arabidopsis is correlated with that in rice. This pattern of parallel expansion is much more pronounced than the whole-genome trend in rice and Arabidopsis. The high rate of expansion among plant TF genes and their propensity for parallel expansion suggest frequent adaptive responses to selection pressure common among higher plants.

Duplication of floral regulatory genes in the Lamiales

Duplication of some floral regulatory genes has occurred repeatedly in angiosperms, whereas others are thought to be single-copy in most lineages. We selected three genes that interact in a pathway regulating floral development conserved among higher tricolpates (LFY/FLO, UFO/FIM, and AP3/DEF) and screened for copy number among families of Lamiales that are closely related to the model species Antirrhinum majus. We show that two of three genes have duplicated at least twice in the Lamiales. Phylogenetic analyses of paralogs suggest that an ancient whole genome duplication shared among many families of Lamiales occurred after the ancestor of these families diverged from the lineage leading to Veronicaceae (including the single-copy species A. majus). Duplication is consistent with previous patterns among angiosperm lineages for AP3/DEF, but this is the first report of functional duplicate copies of LFY/FLO outside of tetraploid species. We propose Lamiales taxa will be good models for understanding mechanisms of duplicate gene preservation and how floral regulatory genes may contribute to morphological diversity.

It's all GO for plant scientists

The Gene Ontology project (http://www.geneontology.org/) produces structured, controlled vocabularies and gene product annotations. Gene products are classified according to the cellular locations and biological process in which they act, and the molecular functions that they carry out. We annotate gene products from a broad range of model species and provide support for those groups that wish to contribute annotation of further model species. The Gene Ontology facilitates the exchange of information between groups of scientists studying similar processes in different model organisms, and so provides a broad range of opportunities for plant scientists.

Constitutively high expression of the histidine biosynthetic pathway contributes to nickel tolerance in hyperaccumulator plants

Plants that hyperaccumulate Ni exhibit an exceptional degree of Ni tolerance and the ability to translocate Ni in large amounts from root to shoot. In hyperaccumulator plants in the genus Alyssum, free His is an important Ni binding ligand that increases in the xylem proportionately to root Ni uptake. To determine the molecular basis of the His response and its contribution to Ni tolerance, transcripts representing seven of the eight enzymes involved in His biosynthesis were investigated in the hyperaccumulator species Alyssum lesbiacum by RNA gel blot analysis. None of the transcripts changed in abundance in either root or shoot tissue when plants were exposed to Ni, but transcript levels were constitutively higher in A. lesbiacum than in the congeneric nonaccumulator A. montanum, especially for the first enzyme in the biosynthetic pathway, ATP-phosphoribosyltransferase (ATP-PRT). Comparison with the weak hyperaccumulator A. serpyllifolium revealed a close correlation between Ni tolerance, root His concentration, and ATP-PRT transcript abundance. Overexpression of an A. lesbiacum ATP-PRT cDNA in transgenic Arabidopsis thaliana increased the pool of free His up to 15-fold in shoot tissue, without affecting the concentration of any other amino acid. His-overproducing lines also displayed elevated tolerance to Ni but did not exhibit increased Ni concentrations in either xylem sap or shoot tissue, suggesting that additional factors are necessary to recapitulate the complete hyperaccumulator phenotype. These results suggest that ATP-PRT expression plays a major role in regulating the pool of free His and contributes to the exceptional Ni tolerance of hyperaccumulator Alyssum species.

Computational identification of 69 retroposons in Arabidopsis

Retroposition is a shot-gun strategy of the genome to achieve evolutionary diversities by mixing and matching coding sequences with novel regulatory elements. We have identified 69 retroposons in the Arabidopsis (Arabidopsis thaliana) genome by a computational approach. Most of them were derivatives of mature mRNAs, and 20 genes contained relics of the reverse transcription process, such as truncations, deletions, and extra sequence additions. Of them, 22 are processed pseudogenes, and 52 genes are likely to be actively transcribed, especially in tissues from apical meristems (roots and flowers). Functional compositions of these retroposon parental genes imply that not the mRNA itself but its expression in gamete cells defines a suitable template for retroposition. The presence/absence patterns of retroposons can be used as cladistic markers for biogeographic research. Effects of human and the Mediterranean Pleistocene refugia in Arabidopsis biogeographic distributions were revealed based on two recent retroposons (At1g61410 and At5g52090). An evolutionary rate of new gene creation by retroposition was calculated as 0.6 genes per million years. Retroposons can also be used as molecular fossils of the parental gene expressions in ancient time. Extensions of 3' untranslated regions for those expressed parental genes are revealed as a possible trend of plant transcriptome evolution. In addition, we reported the first plant functional chimeric gene that adapts to intercompartmental transport by capturing two additional exons after retroposition.

Local coexpression domains of two to four genes in the genome of Arabidopsis

Expression of genes in eukaryotic genomes is known to cluster, but cluster size is generally loosely defined and highly variable. We have here taken a very strict definition of cluster as sets of physically adjacent genes that are highly coexpressed and form so-called local coexpression domains. The Arabidopsis (Arabidopsis thaliana) genome was analyzed for the presence of such local coexpression domains to elucidate its functional characteristics. We used expression data sets that cover different experimental conditions, organs, tissues, and cells from the Massively Parallel Signature Sequencing repository and microarray data (Affymetrix) from a detailed root analysis. With these expression data, we identified 689 and 1,481 local coexpression domains, respectively, consisting of two to four genes with a pairwise Pearson's correlation coefficient larger than 0.7. This number is approximately 1- to 5-fold higher than the numbers expected by chance. A small (5%-10%) yet significant fraction of genes in the Arabidopsis genome is therefore organized into local coexpression domains. These local coexpression domains were distributed over the genome. Genes in such local domains were for the major part not categorized in the same functional category (GOslim). Neither tandemly duplicated genes nor shared promoter sequence nor gene distance explained the occurrence of coexpression of genes in such chromosomal domains. This indicates that other parameters in genes or gene positions are important to establish coexpression in local domains of Arabidopsis chromosomes.