The evolution of plant nuclear genes

Coevolution in Hybrid Genomes: Nuclear-Encoded Rubisco Small Subunits and Their Plastid-Targeting Translocons Accompanying Sequential Allopolyploidy Events in Triticum

The Triticum/Aegilops complex includes hybrid species resulting from homoploid hybrid speciation and allopolyploid speciation. Sequential allotetra- and allohexaploidy events presumably result in two challenges for the hybrids, which involve 1) cytonuclear stoichiometric disruptions caused by combining two diverged nuclear genomes with the maternal inheritance of the cytoplasmic organellar donor; and 2) incompatibility of chimeric protein complexes with diverged subunits from nuclear and cytoplasmic genomes. Here, we describe coevolution of nuclear rbcS genes encoding the small subunits of Rubisco (ribulose 1,5-bisphosphate carboxylase/oxygenase) and nuclear genes encoding plastid translocons, which mediate recognition and translocation of nuclear-encoded proteins into plastids, in allopolyploid wheat species. We demonstrate that intergenomic paternal-to-maternal gene conversion specifically occurred in the genic region of the homoeologous rbcS3 gene from the D-genome progenitor of wheat (abbreviated as rbcS3D) such that it encodes a maternal-like or B-subgenome-like SSU3D transit peptide in allohexaploid wheat but not in allotetraploid wheat. Divergent and limited interaction between SSU3D and the D-subgenomic TOC90D translocon subunit is implicated to underpin SSU3D targeting into the chloroplast of hexaploid wheat. This implicates early selection favoring individuals harboring optimal maternal-like organellar SSU3D targeting in hexaploid wheat. These data represent a novel dimension of cytonuclear evolution mediated by organellar targeting and transportation of nuclear proteins.

Sequencing and variation of terpene synthase gene (TPS2) as the major gene in biosynthesis of thymol in different Thymus species

Thyme (Thymus spp.) is a valuable genus of Lamiaceae family with different pharmaceutical and food properties. Thymol has also been considered as the major essential oil compound in most of the studied Thymus species. In this research, the gene encoding γ-terpinene synthase (Ttps2) was sequenced in T. vulgaris and in eight Iranian thymes including T. carmanicus, T. daenensis, T. fedtschenkoi, T. kotschyanus, T. migricus, T. pubescens, T. serpyllum, and T. trautvetteri. Genetic relationships based on terpene synthase genes were also determined among the studied species. Rapid Amplification of cDNA Ends (RACE) PCR was done to complete the sequence of all species. The cDNA of the studied species possessed an open reading frame ranging from 1788 to 1794 bp that encode for a protein of 596-598 amino acids, presenting all the conserved motifs characteristics of monoterpene synthases. The taxonomic status of Thymus species was determined based on eight reported sections. The species were classified in three major groups. The first and second group comprised species of Micantes and Mastichina sections. The third cluster included the species belonging to Serpyllum and Pseudothymbra sections. Overall, phylogenetic analysis according to whole sequence of Ttps2 gene can help providing insights in respect to its evolutionary process. Finally, clustering based on the amount of main essential oils components (thymol and carvacrol) was compared with that based on Ttps2 gene classification in the studied Thymus species, showing that clustering is not always in accordance.

Cytonuclear Variation of Rubisco in Synthesized Rice Hybrids and Allotetraploids

The allopolyploid speciation process faces the genomic challenge of stoichiometric disruption caused by merging biparental nuclear genomes with only one (usually maternal) of the two sets of progenitor cytoplasmic genomes. The photosynthetic protein 1,5-bisphosphate carboxylase/oxygenase (Rubisco) is composed of nuclear-encoded small subunits (SSUs) and plastome-encoded large subunits (LSUs), making it an ideal enzyme to explore the evolution process of cytonuclear accommodation. We investigated the variation of SSUs and their encoding genes in synthetic nascent rice ( L.) allotetraploid lineages, formed from the parental subspecies and of Asian rice. The evolution of genes in rice subspecies involves both mutation and concerted homogenization. Within reciprocal rice hybrids and allopolyploids, there was no consistent pattern of biased expression of alleles or homeologs, nor was there biased gene conversion favoring the maternal gene copies. Instead, we observed an apparently stochastic pattern of intergenomic gene conversions and biased expression of homeologs. We conclude that in young rice allopolyploids, cytonuclear coordination either is not selectively favored because of high parental sequence similarity or because there has been insufficient time for subtle selective effects to become observable.

Origin and Evolution of Allopolyploid Wheatgrass Elymus fibrosus (Schrenk) Tzvelev (Poaceae: Triticeae) Reveals the Effect of Its Origination on Genetic Diversity

Origin and evolution of tetraploid Elymus fibrosus (Schrenk) Tzvelev were characterized using low-copy nuclear gene Rpb2 (the second largest subunit of RNA polymerase II), and chloroplast region trnL-trnF (spacer between the tRNA Leu (UAA) gene and the tRNA-Phe (GAA) gene). Ten accessions of E. fibrosus along with 19 Elymus species with StH genomic constitution and diploid species in the tribe Triticeae were analyzed. Chloroplast trnL-trnF sequence data suggested that Pseudoroegneria (St genome) was the maternal donor of E. fibrosus. Rpb2 data confirmed the presence of StH genomes in E. fibrosus, and suggested that St and H genomes in E. fibrosus each is more likely originated from single gene pool. Single origin of E. fibrosus might be one of the reasons causing genetic diversity in E. fibrosus lower than those in E. caninus and E. trachycaulus, which have similar ecological preferences and breeding systems with E. fibrosus, and each was originated from multiple sources. Convergent evolution of St and H copy Rpb2 sequences in some accessions of E. fibrosus might have occurred during the evolutionary history of this allotetraploid.

Diversification of R2R3-MYB Transcription Factors in the Tomato Family Solanaceae

MYB transcription factors play an important role in regulating key plant developmental processes involving defense, cell shape, pigmentation, and root formation. Within this gene family, sequences containing an R2R3 MYB domain are the most abundant type and exhibit a wide diversity of functions. In this study, we identify 559 R2R3 MYB genes using whole genome data from four species of Solanaceae and reconstruct their evolutionary relationships. We compare the Solanaceae R2R3 MYBs to the well-characterized Arabidopsis thaliana sequences to estimate functional diversity and to identify gains and losses of MYB clades in the Solanaceae. We identify numerous R2R3 MYBs that do not appear closely related to Arabidopsis MYBs, and thus may represent clades of genes that have been lost along the Arabidopsis lineage or gained after the divergence of Rosid and Asterid lineages. Despite differences in the distribution of R2R3 MYBs across functional subgroups and species, the overall size of the R2R3 subfamily has changed relatively little over the roughly 50 million-year history of Solanaceae. We added our information regarding R2R3 MYBs in Solanaceae to other data and performed a meta-analysis to trace the evolution of subfamily size across land plants. The results reveal many shifts in the number of R2R3 genes, including a 54 % increase along the angiosperm stem lineage. The variation in R2R3 subfamily size across land plants is weakly positively correlated with genome size and strongly positively correlated with total number of genes. The retention of such a large number of R2R3 copies over long evolutionary time periods suggests that they have acquired new functions and been maintained by selection. Discovering the nature of this functional diversity will require integrating forward and reverse genetic approaches on an -omics scale.

De novo sequencing transcriptome of endemic Gentiana straminea (Gentianaceae) to identify genes involved in the biosynthesis of active ingredients

Gentiana straminea is a popular Tibetan medicine that has been used for thousands of years in China to treat various diseases and conditions. Although it has multiple pharmaceutical purposes and important economic plant resource in China, transcriptome and molecular base still known limited. In flowering season, samples were collected from different tissues, using the NGS Illumina. Solexa platform, about 58.85 million sequencing reads were generated and assembled de novo, yielding 78,764 high quality unigenes with an average length of 1090bp. Gene Ontology (GO), KEGG pathway mapping showed that 49,033 of these were identified as putative homologs of annotated sequences in the protein databases. Among them, candidate genes associated with iridoid, flavonoid and anthocyanin were identified. Further the key enzymes involved to iridoid and flavonoid synthesis pathway were analyzed by quantitative real-time polymerase chain reaction (qRT-PCR) on different tissues, the flower and root had the higher expression than leaves. In addition, 7591 SSR markers were identified from the unigenes of the G. straminea transcriptome. The foundation of G. straminea provided the important resource for facilitating to study molecular and functional genomics of it and related this species on the Qinghai-Tibet Plateau.

Origin and Reticulate Evolutionary Process of Wheatgrass Elymus trachycaulus (Triticeae: Poaceae)

To study origin and evolutionary dynamics of tetraploid Elymus trachycaulus that has been cytologically defined as containing StH genomes, thirteen accessions of E. trachycaulus were analyzed using two low-copy nuclear gene Pepc (phosphoenolpyruvate carboxylase) and Rpb2 (the second largest subunit of RNA polymerase II), and one chloroplast region trnL–trnF (spacer between the tRNA Leu (UAA) gene and the tRNA-Phe (GAA) gene). Our chloroplast data indicated that Pseudoroegneria (St genome) was the maternal donor of E. trachycaulus. Rpb2 data indicated that the St genome in E. trachycaulus was originated from either P. strigosa, P. stipifolia, P. spicata or P. geniculate. The Hordeum (H genome)-like sequences of E. trachycaulus are polyphyletic in the Pepc tree, suggesting that the H genome in E. trachycaulus was contributed by multiple sources, whether due to multiple origins or introgression resulting from subsequent hybridization. Failure to recovering St copy of Pepc sequence in most accessions of E. trachycaulus might be caused by genome convergent evolution in allopolyploids. Multiple copies of H-like Pepc sequence from each accession with relative large deletions and insertions might be caused by either instability of Pepc sequence in H- genome or incomplete concerted evolution. Our results highlighted complex evolutionary history of E. trachycaulus.

Cytonuclear evolution of rubisco in four allopolyploid lineages

Allopolyploidization in plants entails the merger of two divergent nuclear genomes, typically with only one set (usually maternal) of parental plastidial and mitochondrial genomes and with an altered cytonuclear stoichiometry. Thus, we might expect cytonuclear coevolution to be an important dimension of allopolyploid evolution. Here, we investigate cytonuclear coordination for the key chloroplast protein rubisco (ribulose 1,5-bisphosphate carboxylase/oxygenase), which is composed of nuclear-encoded, small subunits (SSUs) and plastid-encoded, large subunits. By studying gene composition and diversity as well as gene expression in four model allopolyploid lineages, Arabidopsis, Arachis, Brassica, and Nicotiana, we demonstrate that paralogous nuclear-encoded rbcS genes within diploids are subject to homogenization via gene conversion and that such concerted evolution via gene conversion characterizes duplicated genes (homoeologs) at the polyploid level. Many gene conversions in the polyploids are intergenomic with respect to the diploid progenitor genomes, occur in functional domains of the homoeologous SSUs, and are directionally biased, such that the maternal amino acid states are favored. This consistent preferential maternal-to-paternal gene conversion is mirrored at the transcriptional level, with a uniform transcriptional bias of the maternal-like rbcS homoeologs. These data, repeated among multiple diverse angiosperm genera for an important photosynthetic enzyme, suggest that cytonuclear coevolution may be mediated by intergenomic gene conversion and altered transcription of duplicated, now homoeologous nuclear genes.

A phylogenetic analysis of the genus Fragaria (strawberry) using intron-containing sequence from the ADH-1 gene

The genus Fragaria encompasses species at ploidy levels ranging from diploid to decaploid. The cultivated strawberry, Fragaria×ananassa, and its two immediate progenitors, F. chiloensis and F. virginiana, are octoploids. To elucidate the ancestries of these octoploid species, we performed a phylogenetic analysis using intron-containing sequences of the nuclear ADH-1 gene from 39 germplasm accessions representing nineteen Fragaria species and one outgroup species, Dasiphora fruticosa. All trees from Maximum Parsimony and Maximum Likelihood analyses showed two major clades, Clade A and Clade B. Each of the sampled octoploids contributed alleles to both major clades. All octoploid-derived alleles in Clade A clustered with alleles of diploid F. vesca, with the exception of one octoploid allele that clustered with the alleles of diploid F. mandshurica. All octoploid-derived alleles in clade B clustered with the alleles of only one diploid species, F. iinumae. When gaps encoded as binary characters were included in the Maximum Parsimony analysis, tree resolution was improved with the addition of six nodes, and the bootstrap support was generally higher, rising above the 50% threshold for an additional nine branches. These results, coupled with the congruence of the sequence data and the coded gap data, validate and encourage the employment of sequence sets containing gaps for phylogenetic analysis. Our phylogenetic conclusions, based upon sequence data from the ADH-1 gene located on F. vesca linkage group II, complement and generally agree with those obtained from analyses of protein-encoding genes GBSSI-2 and DHAR located on F. vesca linkage groups V and VII, respectively, but differ from a previous study that utilized rDNA sequences and did not detect the ancestral role of F. iinumae.

Identification and characterization of a second isogene encoding γ-terpinene synthase in Thymus caespititius

Thymus caespititius Brot. is an Iberian endemic species, whose essential oils possess high polymorphism. They consist mostly of mono- and sesquiterpene, some of them with interest for the pharmaceutical and food industries. The search for terpene synthase genes was performed in three in vitro T. caespititius genotypes. For these plants, the expression of a previously described γ-terpinene synthase gene, Tctps2, was confirmed, occurring concomitantly with a new gene encoding an enzyme with similar activity, named Thymus caespititius terpene synthase 4 (Tctps4). The two isogenes were isolated and functionally characterized in the three plant genotypes. Alignment of the two Tctps revealed a transit peptide much shorter in Tctps4 than in Tctps2 (3-4 amino acids instead of 47). The Tctps4 open reading frame is shorter than Tctps2 (1665 bp versus 1794 bp). The amino acid sequence of both γ-terpinene synthases shared an 88% pairwise identity. The fact that T. caespititius carries two isogenes for γ-terpinene synthases, suggests gene duplication along the evolutionary process, followed by mutations leading to the differentiation of both genes. These mutations didn't compromise protein activity. A high accumulation of transcripts from both genes was found in shoots of in vitro plantlets, while in roots they could not be detected. Still, γ-terpinene levels in aerial parts were reduced, probably due to fast conversion into carvacrol and thymol, the main components from T. caespititius essential oils. This study is a contribution to the identification of terpene synthase genes in Lamiaceae.

Purification and enzymatic characterization of alcohol dehydrogenase from Arabidopsis thaliana

Alcohol dehydrogenases (ADH) catalyze the interconversion between alcohols and aldehydes with the reduction of nicotinamide adenine dinucleotide (NAD(+)) to NADH. In this study, for the first time we report an over-expression and purification strategy for the Arabidosis thaliana ADH (AtADH), and characterize its enzymatic properties. AtADH was expressed in an Escherichia coli system, the polyhistidine-tag was removed after the recombinant AtADH protein was purified by metal chelating affinity chromatography. Activity assays demonstrated that AtADH has distinct enzymatic properties when compared with many well-known ADHs. It held peak activity at pH 10.5 and showed broad substrate selectivity for primary and secondary alcohols. The kinetic Km parameters for both ethanol and coenzyme were in the order of mM. This relative low affinity may reflect the need of the plant to maintain a supply of NAD(+) in nature. Different from yeast ADH, AtADH showed almost the same activity for short straight chain alcohols and reduced activity for secondary alcohols. This broad spectrum in alcohol selection and the observed higher catalytic activity (high Vmax (EtOH)) may result from the requirement of the single enzyme to accommodate many substrates.

The evolution of phenylpropanoid metabolism in the green lineage

Phenolic secondary metabolites are only produced by plants wherein they play important roles in both biotic and abiotic defense in seed plants as well as being potentially important bioactive compounds with both nutritional and medicinal benefits reported for animals and humans as a consequence of their potent antioxidant activity. During the long evolutionary period in which plants have adapted to the environmental niches in which they exist (and especially during the evolution of land plants from their aquatic algal ancestors), several strategies such as gene duplication and convergent evolution have contributed to the evolution of this pathway. In this respect, diversity and redundancy of several key genes of phenolic secondary metabolism such as polyketide synthases, cytochrome P450s, Fe(2+)/2-oxoglutarate-dependent dioxygenases and UDP-glycosyltransferases have played an essential role. Recent technical developments allowing affordable whole genome sequencing as well as a better inventory of species-by-species chemical diversity have resulted in a dramatic increase in the number of tools we have to assess how these pathways evolved. In parallel, reverse genetics combined with detailed molecular phenotyping is allowing us to elucidate the functional importance of individual genes and metabolites and by this means to provide further mechanistic insight into their biological roles. In this review, phenolic metabolite-related gene sequences (for a total of 65 gene families including shikimate biosynthetic genes) are compared across 23 independent species, and the phenolic metabolic complement of various plant species are compared with one another, in attempt to better understand the evolution of diversity in this crucial pathway.

Expansion of the gamma-gliadin gene family in Aegilops and Triticum

BACKGROUND

The gamma-gliadins are considered to be the oldest of the gliadin family of storage proteins in Aegilops/Triticum. However, the expansion of this multigene family has not been studied in an evolutionary perspective.

RESULTS

We have cloned 59 gamma-gliadin genes from Aegilops and Triticum species (Aegilops caudata L., Aegilops comosa Sm. in Sibth. & Sm., Aegilops mutica Boiss., Aegilops speltoides Tausch, Aegilops tauschii Coss., Aegilops umbellulata Zhuk., Aegilops uniaristata Vis., and Triticum monococcum L.) representing eight different genomes: Am, B/S, C, D, M, N, T and U. Overall, 15% of the sequences contained internal stop codons resulting in pseudogenes, but this percentage was variable among genomes, up to over 50% in Ae. umbellulata. The most common length of the deduced protein, including the signal peptide, was 302 amino acids, but the length varied from 215 to 362 amino acids, both obtained from Ae. speltoides. Most genes encoded proteins with eight cysteines. However, all Aegilops species had genes that encoded a gamma-gliadin protein of 302 amino acids with an additional cysteine. These conserved nine-cysteine gamma-gliadins may perform a specific function, possibly as chain terminators in gluten network formation in protein bodies during endosperm development. A phylogenetic analysis of gamma-gliadins derived from Aegilops and Triticum species and the related genera Lophopyrum, Crithopsis, and Dasypyrum showed six groups of genes. Most Aegilops species contained gamma-gliadin genes from several of these groups, which also included sequences from the genera Lophopyrum, Crithopsis, and Dasypyrum. Hordein and secalin sequences formed separate groups.

CONCLUSIONS

We present a model for the evolution of the gamma-gliadins from which we deduce that the most recent common ancestor (MRCA) of Aegilops/Triticum-Dasypyrum-Lophopyrum-Crithopsis already had four groups of gamma-gliadin sequences, presumably the result of two rounds of duplication of the locus.

The cytonuclear dimension of allopolyploid evolution: an example from cotton using rubisco

During allopolyploid speciation, two divergent nuclear genomes merge, yet only one (usually the maternal) of the two sets of progenitor organellar genomes is maintained. Rubisco (1,5-bisphosphate carboxylase/oxygenase) is composed of nuclear-encoded small subunits (SSUs) and plastome-encoded large subunits (LSUs), providing an ideal system to explore the evolutionary process of cytonuclear accommodation. Here, we take initial steps in this direction, using Gossypium allopolyploids as our model. SSU copies from divergent (5-10 My) progenitor diploids ("A" and "D" genomes) were combined at the time of polyploid formation 1-2 Ma, with the LSU encoded by the maternal A-genome parent. LSU genes from A- and D-genome diploids and AD-genome allopolyploids were sequenced, revealing several nonsynonymous substitutions and suggesting the possibility of differential selection on the nuclear-encoded rbcS partner following allopolyploid formation. Sequence data for the rbcS gene family revealed nonreciprocal homoeologous recombination between A- and D-rbcS homoeologs in all polyploid species but not in a synthetic intergenomic F1 hybrid, demonstrating "gene conversion" during allopolyploid evolution. All progenitor rbcS genes are retained and expressed in the five extant allopolyploid species, but analysis of the leaf transcriptome showed that A-homoeologs are preferentially expressed in both the allopolyploid and hybrid, consistent with the maternal origin of rbcL. Although rbcS genes from both progenitor genomes are expressed, some appear to have experienced mutations that may represent cytonuclear coevolution.

Xanthomonas oryzae pv oryzae triggers immediate transcriptomic modulations in rice

BACKGROUND

Xanthomonas oryzae pv oryzae is a devastating pathogen of rice and has been extensively studied as a model pathogen of monocotyledons. Expressional studies in both the contenders have been undertaken in past to understand the molecular mechanism underlying the compatible and incompatible interactions in the pathosystem. Continuous update on database and gene annotations necessitates constant updating on the roles of the new entities as well as reinterpretation of regulations of the previous ones. Moreover the past endeavors have addressed the middle or late defense responses of the rice plant whereas in the present study an attempt has been made to investigate the early defense responses taking place immediately after inoculation.

RESULTS

Microarray was used to study the transcriptional modulations in eighteen days old rice seedling leaves of both susceptible and resistant genotypes one hour after inoculation. In resistant plants as compared to susceptible ones 274 genes were found to be differentially expressed. Annotations could be assigned to 112 up- and 73 down-regulated transcripts and gene interaction maps were generated for 86 transcripts. Expressional data and interaction maps were used to develop a hypothetical scheme of the molecular events taking place during early defense response. Network analysis with the differential transcripts showed up-regulation of major clusters of cell signaling proteins and transcription factors while growth and basal metabolic components were largely found to be down-regulated.

CONCLUSIONS

This study provides an understanding of the early defense signaling in rice cells. Components of the calcium and lipid signaling as well as MAPK cascade were modulated, by signals from surface receptors and cytosolic R-proteins, to arouse jasmonic acid and ethylene signaling and suppress auxin signaling through various transcription factors. Abscisic acid modulation was also evident through the expression regulation of transcription factors involved with its functions. Moreover adjustments in expression levels of components of primary as well as secondary metabolism, protein trafficking and turnout were apparent, highlighting the complexity of defense response.

The evolution of the plastid chromosome in land plants: gene content, gene order, gene function

This review bridges functional and evolutionary aspects of plastid chromosome architecture in land plants and their putative ancestors. We provide an overview on the structure and composition of the plastid genome of land plants as well as the functions of its genes in an explicit phylogenetic and evolutionary context. We will discuss the architecture of land plant plastid chromosomes, including gene content and synteny across land plants. Moreover, we will explore the functions and roles of plastid encoded genes in metabolism and their evolutionary importance regarding gene retention and conservation. We suggest that the slow mode at which the plastome typically evolves is likely to be influenced by a combination of different molecular mechanisms. These include the organization of plastid genes in operons, the usually uniparental mode of plastid inheritance, the activity of highly effective repair mechanisms as well as the rarity of plastid fusion. Nevertheless, structurally rearranged plastomes can be found in several unrelated lineages (e.g. ferns, Pinaceae, multiple angiosperm families). Rearrangements and gene losses seem to correlate with an unusual mode of plastid transmission, abundance of repeats, or a heterotrophic lifestyle (parasites or myco-heterotrophs). While only a few functional gene gains and more frequent gene losses have been inferred for land plants, the plastid Ndh complex is one example of multiple independent gene losses and will be discussed in detail. Patterns of ndh-gene loss and functional analyses indicate that these losses are usually found in plant groups with a certain degree of heterotrophy, might rendering plastid encoded Ndh1 subunits dispensable.

The plant ADH gene family

The structures, evolution and functions of alcohol dehydrogenase gene families and their products have been scrutinized for half a century. Our understanding of the enzyme structure and catalytic activity of plant alcohol dehydrogenase (ADH-P) is based on the vast amount of information available for its animal counterpart. The probable origins of the enzyme from a simple β-coil and eventual emergence from a glutathione-dependent formaldehyde dehydrogenase have been well described. There is compelling evidence that the small ADH gene families found in plants today are the survivors of multiple rounds of gene expansion and contraction. To the probable original function of their products in the terminal reaction of anaerobic fermentation have been added roles in yeast-like aerobic fermentation and the production of characteristic scents that act to attract animals that serve as pollinators or agents of seed dispersal and to protect against herbivores.

Adaptive evolution of digestive RNASE1 genes in leaf-eating monkeys revisited: new insights from ten additional colobines

Pancreatic RNase genes implicated in the adaptation of the colobine monkeys to leaf eating have long intrigued evolutionary biologists since the identification of a duplicated RNASE1 gene with enhanced digestive efficiencies in Pygathrix nemaeus. The recent emergence of two contrasting hypotheses, that is, independent duplication and one-duplication event hypotheses, make it into focus again. Current understanding of Colobine RNASE1 gene evolution of colobine monkeys largely depends on the analyses of few colobine species. The present study with more intensive taxonomic and character sampling not only provides a clearer picture of Colobine RNASE1 gene evolution but also allows to have a more thorough understanding about the molecular basis underlying the adaptation of Colobinae to the unique leaf-feeding lifestyle. The present broader and detailed phylogenetic analyses yielded two important findings: 1) All trees based on the analyses of coding, noncoding, and both regions provided consistent evidence, indicating RNASE1 duplication occurred after Asian and African colobines speciation, that is, independent duplication hypothesis; 2) No obvious evidence of gene conversion in RNASE1 gene was found, favoring independent evolution of Colobine RNASE1 gene duplicates. The conclusion drawn from previous studies that gene conversion has played a significant role in the evolution of Colobine RNASE1 was not supported. Our selective constraint analyses also provided interesting insights, with significant evidence of positive selection detected on ancestor lineages leading to duplicated gene copies. The identification of a handful of new adaptive sites and amino acid changes that have not been characterized previously also provide a necessary foundation for further experimental investigations of RNASE1 functional evolution in Colobinae.

Characterization and differential expression analysis of artichoke phenylalanine ammonia-lyase-coding sequences

Sequences encoding phenylalanine ammonia-lyase were isolated from artichoke, by using a sequence homology strategy, by screening a genomic library and by 3'-rapid amplification of cDNA end (RACE) technology. These analyses and Southern blots suggested that, in artichoke, phenylalanine ammonia-lyase (PAL) is encoded by a small gene family. The sequences isolated from genomic DNA possess two exons and one intron at the conserved position as in most plant pal characterized to date. The 3'-RACE analysis also indicated that each member of the artichoke pal gene family was present as a pool of transcripts, different in the length of 3'-untranslated region. The deduced amino acid sequences were highly similar to those of PAL from lettuce and sunflower. One of the artichoke pal genes was completely sequenced, and its 5' upstream region contained TATA, CAAT box and cis regulatory elements identified in other phenylpropanoid pathway genes as playing a role in UV and elicitor induction. The expression of three of the identified artichoke pal sequences was evaluated in different plant parts, in developmental stages and after wounding, using gene-specific primers/probe combinations in real-time polymerase chain reaction assays. The three putative genes were differentially expressed in the plant parts analysed and were developmentally regulated. Moreover, after leaf mechanical injury, all of them were differentially regulated. The possible involvement of the single pal genes in different physiological processes is discussed.

Phylogenetic and molecular analysis of the ribulose-1,5-bisphosphate carboxylase small subunit gene family in banana

Despite being the number one fruit crop in the world, very little is known about the phylogeny and molecular biology of banana (Musa spp.). Six banana rbcS gene families encoding the small subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase from six different Musa spp. are presented. For a comprehensive phylogenetic study using Musa rbcS genes, a total of 57 distinct rbcS sequences was isolated from six accessions that contained different combinations of the A and B ancestral/parental genomes. As a result, five of the six members of the rbcS gene family could be affiliated with the A and/or B Musa genomes and at least three of the six gene families most likely existed before Musa A and B genomes separated. By combining sequence data with quantitative real-time PCR it was determined that the different Musa rbcS gene family members are also often multiply represented in each genome, with the highest copy numbers in the B genome. Expression of some of the rbcS genes varied in intensity and in different tissues indicating differences in regulation. To analyse and compare regulatory sequences of Musa rbcS genes, promoter and terminator regions were cloned for three Musa rbcS genes. Transient transformation assays using promoter-reporter-terminator constructs in maize, wheat, and sugarcane demonstrated that the rbcS-Ma1, rbcS-Ma3, and rbcS-Ma5 promoters could be useful for transgene expression in heterologous expression systems. Furthermore, the rbcS-Ma1 terminator resulted in a 2-fold increase of transgene expression when directly compared with the widely used Nos terminator.

Phytochrome-mediated development in land plants: red light sensing evolves to meet the challenges of changing light environments

Phytochromes are photoreceptors that provide plants with circadian, seasonal, and positional information critical for the control of germination, seedling development, shade avoidance, reproduction, dormancy, and sleep movements. Phytochromes are unique among photoreceptors in their capacity to interconvert between a red-absorbing form (absorption maximum of approximately 660 nm) and a far-red absorbing form (absorption maximum of approximately 730 nm), which occur in a dynamic equilibrium within plant cells, corresponding to the proportions of red and far-red energy in ambient light. Because pigments in stems and leaves absorb wavelengths below about 700 nm, this provides plants with an elegant system for detecting their position relative to other plants, with which the plants compete for light. Certain aspects of phytochrome-mediated development outside of flowering plants are strikingly similar to those that have been characterized in Arabidopsis thaliana and other angiosperms. However, early diverging land plants have fewer distinct phytochrome gene lineages, suggesting that both diversification and subfunctionalization have been important in the evolution of the phytochrome gene family. There is evidence that subfunctionalization proceeded by the partitioning among paralogues of photosensory specificity, physiological response modes, and light-regulated gene expression and protein stability. Parallel events of duplication and functional divergence may have coincided with the evolution of canopy shade and the increasing complexity of the light environment. Within angiosperms, patterns of functional divergence are clade-specific and the roles of phytochromes in A. thaliana change across environments, attesting to the evolutionary flexibility and contemporaneous plasticity of phytochrome signalling in the control of development.

Experimental reconstruction of functional gene transfer from the tobacco plastid genome to the nucleus

Eukaryotic cells arose through the uptake of free-living bacteria by endosymbiosis and their gradual conversion into organelles (plastids and mitochondria). Capture of the endosymbionts was followed by massive translocation of their genes to the genome of the host cell. How genes were transferred from the (prokaryotic) organellar genome to the (eukaryotic) nuclear genome and how the genes became functional in their new eukaryotic genetic environment is largely unknown. Here, we report the successful experimental reconstruction of functional gene transfer between an organelle and the nucleus, a process that normally occurs only on large evolutionary timescales. In consecutive genetic screens, we first transferred a chloroplast genome segment to the nucleus and then selected for gene activation in the nuclear genome. We show that DNA-mediated gene transfer can give rise to functional nuclear genes if followed by suitable rearrangements in the nuclear genome. Acquisition of gene function involves (1) transcriptional activation by capture of the promoter of an upstream nuclear gene and (2) utilization of AT-rich noncoding sequences downstream of the plastid gene as RNA cleavage and polyadenylation sites. Our results reveal the molecular mechanisms of how organellar DNA transferred to the nucleus gives rise to functional genes and reproduce in the laboratory a key process in the evolution of eukaryotic cells.

Functional analyses of two tomato APETALA3 genes demonstrate diversification in their roles in regulating floral development

The floral homeotic APETALA3 (AP3) gene in Arabidopsis thaliana encodes a MADS box transcription factor required for specifying petal and stamen identities. AP3 is a member of the euAP3 lineage, which arose by gene duplication coincident with radiation of the core eudicots. Although Arabidopsis lacks genes in the paralogous Tomato MADS box gene 6 (TM6) lineage, tomato (Solanum lycopersicum) possesses both euAP3 and TM6 genes, which have functionally diversified. A loss-of-function mutation in Tomato AP3 (TAP3) resulted in homeotic transformations of both petals and stamens, whereas RNA interference-induced reduction in TM6 function resulted in flowers with homeotic defects primarily in stamens. The functional differences between these genes can be ascribed partly to different expression domains. When overexpressed in an equivalent domain, both genes can partially rescue the tap3 mutant, indicating that relative levels as well as spatial patterns of expression contribute to functional differences. Our results also indicate that the two proteins have differing biochemical capabilities. Together, these results suggest that TM6 and TAP3 play qualitatively different roles in floral development; they also support the ideas that the ancestral role of AP3 lineage genes was in specifying stamen development and that duplication and divergence in the AP3 lineage allowed for the acquisition of a role in petal specification in the core eudicots.

Genome evolution of allopolyploids: a process of cytological and genetic diploidization

Allopolyploidy is a prominent mode of speciation in higher plants. Due to the coexistence of closely related genomes, a successful allopolyploid must have the ability to invoke and maintain diploid-like behavior, both cytologically and genetically. Recent studies on natural and synthetic allopolyploids have raised many discrepancies. Most species have displayed non-Mendelian behavior in the allopolyploids, but others have not. Some species have demonstrated rapid genome changes following allopolyploid formation, while others have conserved progenitor genomes. Some have displayed directed, non-random genome changes, whereas others have shown random changes. Some of the genomic changes have appeared in the F1 hybrids, which have been attributed to the union of gametes from different progenitors, while other changes have occurred during or after genome doubling. Although these observations provide significant novel insights into the evolution of allopolyploids, the overall mechanisms of the event are still elusive. It appears that both genetic and epigenetic operations are involved in the diploidization process of allopolyploids. Overall, genetic and epigenetic variations are often associated with the activities of repetitive sequences and transposon elements. Specifically, genomic sequence elimination and chromosome rearrangement are probably the major forces guiding cytological diploidization. Gene non-functionalization, sub-functionalization, neo-functionalization, as well as other kinds of epigenetic modifications, are likely the leading factors promoting genetic diploidization.

Gene structure and molecular analysis of the laccase-like multicopper oxidase (LMCO) gene family in Arabidopsis thaliana

Completed genome sequences have made it clear that multicopper oxidases related to laccase are widely distributed as multigene families in higher plants. Laccase-like multicopper oxidase (LMCO) sequences culled from GenBank and the Arabidopsis thaliana genome, as well as those from several newly cloned genes, were used to construct a gene phylogeny that clearly divided plant LMCOs into six distinct classes, at least three of which predate the evolutionary divergence of angiosperms and gymnosperms. Alignments of the predicted amino acid sequences highlighted regions of variable sequence flanked by the highly conserved copper-binding domains that characterize members of this enzyme family. All of the predicted proteins contained apparent signal sequences. The expression of 13 of the 17 LMCO genes in A. thaliana was assessed in different tissues at various stages of development using RT-PCR. A diversity of expression patterns was demonstrated with some genes being expressed in a constitutive fashion, while others were only expressed in specific tissues at a particular stage of development. Only a few of the LMCO genes were expressed in a pattern that could be considered consistent with a major role for these enzymes in lignin deposition. These results are discussed in the context of other potential physiological functions for plant LMCOs, such as iron metabolism and wound healing.

Molecular phylogeny and evolution of alcohol dehydrogenase (Adh) genes in legumes

BACKGROUND

Nuclear genes determine the vast range of phenotypes that are responsible for the adaptive abilities of organisms in nature. Nevertheless, the evolutionary processes that generate the structures and functions of nuclear genes are only now be coming understood. The aim of our study is to isolate the alcohol dehydrogenase (Adh) genes in two distantly related legumes, and use these sequences to examine the molecular evolutionary history of this nuclear gene.

RESULTS

We isolated the expressed Adh genes from two species of legumes, Sophora flavescens Ait. and Wisteria floribunda DC., by a RT-PCR based approach and found a new Adh locus in addition to homologues of the Adh genes found previously in legumes. To examine the evolution of these genes, we compared the species and gene trees and found gene duplication of the Adh loci in the legumes occurred as an ancient event.

CONCLUSION

This is the first report revealing that some legume species have at least two Adh gene loci belonging to separate clades. Phylogenetic analyses suggest that these genes resulted from relatively ancient duplication events.

Phylogenetic analysis of the nuclear alcohol dehydrogenase (Adh) gene family in Carex section Acrocystis (Cyperaceae) and combined analyses of Adh and nuclear ribosomal ITS and ETS sequences for inferring species relationships

We analyzed sequence variation for the alcohol dehydrogenase (Adh) gene family in Carex section Acrocystis (Cyperaceae) to reconstruct Adh gene trees for Acrocystis species and to characterize the structure of the Adh gene family in Carex. Two Adh loci were included with ITS and ETS sequences in a combined Bayesian inference analysis of Carex section Acrocystis to gain a better understanding of species relationships in the section. In addition, we comment on how the results presented here contribute to our knowledge of the birth-death process of the Adh gene family in angiosperms. It appears that the structure of the Adh gene family in Carex is complex with possibly six loci present in the gene family. Additionally, variation among Acrocystis species within loci is quite low, and there is little phylogenetic resolution in the individual datasets. Bayesian inference analysis of the combined ITS, ETS, Adh1, and Adh2 datasets resulted in a moderately well-supported phylogenetic hypothesis of relationships in the section which is discussed in relation to previous hypotheses of relationships.

Comparative population genetics of the panicoid grasses: sequence polymorphism, linkage disequilibrium and selection in a diverse sample of sorghum bicolor

Levels of genetic variation and linkage disequilibrium (LD) are critical factors in association mapping methods as well as in identification of loci that have been targets of selection. Maize, an outcrosser, has a high level of sequence variation and a limited extent of LD. Sorghum, a closely related but largely self-pollinating panicoid grass, is expected to have higher levels of LD. As a first step in estimation of population genetic parameters in sorghum, we surveyed 27 diverse S. bicolor accessions for sequence variation at a total of 29,186 bp in 95 short regions derived from genetically mapped RFLPs located throughout the genome. Consistent with its higher level of inbreeding, the extent of LD is at least severalfold greater in sorghum than in maize. Total sequence variation in sorghum is about fourfold lower than that in maize, while synonymous variation is fivefold lower, suggesting a smaller effective population size in sorghum. Because we surveyed a species-wide sample, the mating system, which primarily affects population-level diversity, may not be primarily responsible for this difference. Comparisons of polymorphism and divergence suggest that both directional and diversifying selection have played important roles in shaping variation in the sorghum genome.

Evolution of a RNA Polymerase Gene Family in Silene (Caryophyllaceae)—Incomplete Concerted Evolution and Topological Congruence Among Paralogues

Four low-copy nuclear DNA intron regions from the second largest subunits of the RNA polymerase gene family (RPA2, RPB2, RPD2a, and RPD2b), the internal transcribed spacers (ITSs) from the nuclear ribosomal regions, and the rps16 intron from the chloroplast were sequenced and used in a phylogenetic analysis of 29 species from the tribe Sileneae (Caryophyllaceae). We used a low stringency nested polymerase chain reaction (PCR) approach to overcome the difficulties of constructing specific primers for amplification of the low copy nuclear DNA regions. Maximum parsimony analyses resulted in largely congruent phylogenetic trees for all regions. We tested overall model congruence in a likelihood context using the software PLATO and found that ITSs, RPA2, and RPB2 deviated from the maximum likelihood model for the combined data. The topology parameter was then isolated and topological congruence assessed by nonparametric bootstrapping. No strong topological incongruence was found. The analysis of the combined data sets resolves previously poorly known major relationships within Sileneae. Two paralogues of RPD2 were found, and several independent losses and incomplete concerted evolution were inferred. The among-site rate variation was significantly lower in the RNA polymerase introns than in the rps16 intron and ITSs, a property that is attractive in phylogenetic analyses.

Unintended effects and their detection in genetically modified crops

The commercialisation of GM crops in Europe is practically non-existent at the present time. The European Commission has instigated changes to the regulatory process to address the concerns of consumers and member states and to pave the way for removing the current moratorium. With regard to the safety of GM crops and products, the current risk assessment process pays particular attention to potential adverse effects on human and animal health and the environment. This document deals with the concept of unintended effects in GM crops and products, i.e. effects that go beyond that of the original modification and that might impact primarily on health. The document first deals with the potential for unintended effects caused by the processes of transgene insertion (DNA rearrangements) and makes comparisons with genetic recombination events and DNA rearrangements in traditional breeding. The document then focuses on the potential value of evolving "profiling" or "omics" technologies as non-targeted, unbiased approaches, to detect unintended effects. These technologies include metabolomics (parallel analysis of a range of primary and secondary metabolites), proteomics (analysis of polypeptide complement) and transcriptomics (parallel analysis of gene expression). The technologies are described, together with their current limitations. Importantly, the significance of unintended effects on consumer health are discussed and conclusions and recommendations presented on the various approaches outlined.

Conservation and divergence of ASK1 and ASK2 gene functions during male meiosis in Arabidopsis thaliana

Selective proteolysis of regulatory proteins mediated by the ubiquitin pathway is an important mechanism for controlling many biological events. The SCF (Skpl-Cullin-F-box protein) class of E3 ubiquitin ligases controls the ubiquitination of a wide variety of substrates, thereby mediating their degradation by the 26S proteasome. The Arabidopsis genome contains 21 genes encoding Skp1-like proteins that are named as ASKs (Arabidopsis Skp1-like). So far, only the ASK1 gene has been characterized genetically, and is known to be required for male meiosis, flower development, and auxin response. The ASK2 gene is most similar to ASK1 in terms of both the amino acid sequence and expression pattern. To compare ASK2 with ASK1 functionally in male meiosis, different transgenic lines over-expressing ASK1 and ASK2 were tested for their ability to complement the male meiosis defect of the ask1-1 mutant. The genomic ASK1 rescued the ask1-1 mutant defects. The 35S::ASK1 transgene restored male fertility to the ask1-1 mutant, although the percentages of normal pollen grains and tetrads were reduced. 35S::ASK2 lines in the ask1-1 background exhibited partial fertility with even fewer normal pollen grains and tetrads than those of the 35S::ASK1 lines. Detailed analysis of chromosome behavior during male meiosis demonstrated that 35S::ASK1 and 35S::ASK2 lines had different fractions of pollen mother cells undergoing normal meiosis. Our results suggest that ASK2 partially substitutes for ASK1 if expressed at higher than normal levels.

Evolutionary implications of infra- and interspecific molecular variability of pathogenesis-related proteins

We have examined phylogenetic relationships in seven pathogenesis-related (PR) protein families. Within-family comparisons involved 79 species, 166 amino acid sequences, and 1,791 sites. For 37 species, 124 different PR isoforms were identified (an average of 3.3 per species). Thirty-one of the 37 species investigated tended to cluster together (84%). Of the 17 clusters distinguished in the seven phylogenetic trees, 10 (59%) were in agreement with their taxonomic status, ascertained at the family level. The strong similarities among the intraspecific forms, as compared to interspecific differences, argue for some kind of gene conversion, but the rare occurrence of widely different isoforms also suggests diversifying selection. PRs 1, 6, and 4 seem to be less differentiated than PRs 3, 2, 10, and 5.

Role of the small subunit in ribulose-1,5-bisphosphate carboxylase/oxygenase

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) catalyzes the rate-limiting step of CO2 fixation in photosynthesis, but O2 competes with CO2 for substrate ribulose 1,5-bisphosphate, leading to the loss of fixed carbon. Interest in genetically engineering improvements in carboxylation catalytic efficiency and CO2/O2 specificity has focused on the chloroplast-encoded large subunit because it contains the active site. However, there is another type of subunit in the holoenzyme of plants, which, like the large subunit, is present in eight copies. The role of these nuclear-encoded small subunits in Rubisco structure and function is poorly understood. Small subunits may have originated during evolution to concentrate large-subunit active sites, but the extensive divergence of structures among prokaryotes, algae, and land plants seems to indicate that small subunits have more-specialized functions. Furthermore, plants and green algae contain families of differentially expressed small subunits, raising the possibility that these subunits may regulate the structure or function of Rubisco. Studies of interspecific hybrid enzymes have indicated that small subunits are required for maximal catalysis and, in several cases, contribute to CO2/O2 specificity. Although small-subunit genetic engineering remains difficult in land plants, directed mutagenesis of cyanobacterial and green-algal genes has identified specific structural regions that influence catalytic efficiency and CO2/O2 specificity. It is thus apparent that small subunits will need to be taken into account as strategies are developed for creating better Rubisco enzymes.

Duplication and adaptive evolution of the chalcone synthase genes of Dendranthema (Asteraceae)

Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoids, which are important for the pigmentation of flowers and act as attractants to the pollinators. Genes encoding CHS constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. Plants of the Dendranthema genus have white, yellow, and pink flowers, exhibiting considerable variation in flower color. In this article, 18 CHS genes from six Dendranthema species were sequenced. Two of them were found to be pseudogenes. The functional Dendranthema CHS genes formed three well-supported subfamilies: SF1, SF2, and SF3. The inferred phylogeny of the CHS genes of Dendranthema and Gerbera suggests that those genes originated as a result of duplications before divergence of these two genera, and the function of Dendranthema CHS genes have diverged in a similar fashion to the Gerbera CHS genes; i.e., the genes of SF1 and SF3 code for typical CHS enzymes expressed during different stages of development, whereas the genes of SF2 code for another enzyme that is different from CHS in substrate specificity and reaction. Relative rate tests revealed that the Dendranthema CHS genes significantly deviated from clocklike evolution at nonsynonymous sites. Maximum likelihood analysis showed that the nonsynonymous-synonymous (omega = d(N)/d(S)) rate ratio for the lineage ancestral to SF2 was much higher than for other lineages, with some sites having a ratio well above one. Positive selective pressure appears to have driven the divergence of SF2 from SF1 and SF3.

Co-induction of glutathione-S-transferases and multidrug resistance associated protein by xenobiotics in wheat

Herbicide safeners are known to protect monocotyledonous crops from herbicide injury by accelerating the metabolism of herbicides. We have investigated the effects of the safener cloquintocetmexyl, which protects small-grain cereals against the graminicidal herbicide, clodinafop-propargyl. Subtractive suppression hybridisation was used to identify wheat genes which are up-regulated by treatment not only with cloquintocet-mexyl but also with phenobarbital, which is known to stimulate xenobiotic metabolism in animals and plants. DNA sequences of five glutathione transferases (GSTs) belonging to three different classes and a multidrug resistance associated protein (MRP) homologue were identified in the screen. The chemical inducibility of these clones was confirmed by Northern analysis. The MRP protein was shown to be induced by treatments with cloquintocet-mexyl and phenobarbital and to be localised to the tonoplast. Since clodinafop-propargyl is not known to be metabolised by glutathionylation, the significance of GST induction is interpreted in terms of a generalised response to chemical stress, particularly the generation of active oxygen species. This work establishes herbicide safeners as useful tools for the identification of genes encoding herbicide-metabolising enzymes.

Patterns of nucleotide substitution among simultaneously duplicated gene pairs in Arabidopsis thaliana

We characterized rates and patterns of synonymous and nonsynonymous substitution in 242 duplicated gene pairs on chromosomes 2 and 4 of Arabidopsis thaliana. Based on their collinear order along the two chromosomes, the gene pairs were likely duplicated contemporaneously, and therefore comparison of genetic distances among gene pairs provides insights into the distribution of nucleotide substitution rates among plant nuclear genes. Rates of synonymous substitution varied 13.8-fold among the duplicated gene pairs, but 90% of gene pairs differed by less than 2.6-fold. Average nonsynonymous rates were approximately fivefold lower than average synonymous rates; this rate difference is lower than that of previously studied nonplant lineages. The coefficient of variation of rates among genes was 0.65 for nonsynonymous rates and 0.44 for synonymous rates, indicating that synonymous and nonsynonymous rates vary among genes to roughly the same extent. The causes underlying rate variation were explored. Our analyses tentatively suggest an effect of physical location on synonymous substitution rates but no similar effect on nonsynonymous rates. Nonsynonymous substitution rates were negatively correlated with GC content at synonymous third codon positions, and synonymous substitution rates were negatively correlated with codon bias, as observed in other systems. Finally, the 242 gene pairs permitted investigation of the processes underlying divergence between paralogs. We found no evidence of positive selection, little evidence that paralogs evolve at different rates, and no evidence of differential codon usage or third position GC content.

Differential evolutionary dynamics of duplicated paralogous Adh loci in allotetraploid cotton (Gossypium)

Levels and patterns of nucleotide diversity vary widely among lineages. Because allopolyploid species contain duplicated (homoeologous) genes, studies of nucleotide diversity at homoeologous loci may facilitate insight into the evolutionary dynamics of duplicated loci. In this study, we describe patterns of sequence diversity from an alcohol dehydrogenase homoeologous locus pair (AdhC) in allotetraploid cotton (Gossypium, Malvaceae). These data are compared with equivalent information from another homoeologous alcohol dehydrogenase gene pair (AdhA, Small, Ryburn, and Wendel 1999. Mol. Biol. Evol. 16:491-501) which has an overall slower evolutionary rate than AdhC. As expected from the predicted correlation between nucleotide diversity and evolutionary rate, nucleotide diversity was higher for AdhC than for AdhA. In addition, nucleotide diversity is higher in the D-subgenome of allotetraploid cotton for AdhC, confirming earlier observations for AdhA. These observations indicate that for these two pairs of Adh loci, the null hypothesis of equivalent evolutionary dynamics for duplicated genes in allotetraploid cotton is rejected.

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

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

Rubisco: structure, regulatory interactions, and possibilities for a better enzyme

Ribulose-1,5-bisphosphate (RuBP) carboxylase/oxygenase (Rubisco) catalyzes the first step in net photosynthetic CO2 assimilation and photorespiratory carbon oxidation. The enzyme is notoriously inefficient as a catalyst for the carboxylation of RuBP and is subject to competitive inhibition by O2, inactivation by loss of carbamylation, and dead-end inhibition by RuBP. These inadequacies make Rubisco rate limiting for photosynthesis and an obvious target for increasing agricultural productivity. Resolution of X-ray crystal structures and detailed analysis of divergent, mutant, and hybrid enzymes have increased our insight into the structure/function relationships of Rubisco. The interactions and associations relatively far from the Rubisco active site, including regulatory interactions with Rubisco activase, may present new approaches and strategies for understanding and ultimately improving this complex enzyme.

Evolution of intron/exon structure of DEAD helicase family genes in Arabidopsis, Caenorhabditis, and Drosophila

The DEAD box RNA helicase (RH) proteins are homologs involved in diverse cellular functions in all of the organisms from prokaryotes to eukaryotes. Nevertheless, there is a lack of conservation in the splicing pattern in the 53 Arabidopsis thaliana (AtRHs), the 32 Caenorhabditis elegans (CeRHs) and the 29 Drosophila melanogaster (DmRHs) genes. Of the 153 different observed intron positions, 4 are conserved between AtRHs, CeRHs, and DmRHs, and one position is also found in RHs from yeast and human. Of the 27 different AtRH structures with introns, 20 have at least one predicted ancient intron in the regions coding for the catalytic domain. In all of the organisms examined, we found at least one gene with most of its intron predicted to be ancient. In A. thaliana, the large diversity in RH structures suggests that duplications of the ancestral RH were followed by a high number of intron deletions and additions. The very high bias toward phase 0 introns is in favor of intron addition, preferentially in phase 0. Results from this comparative study of the same gene family in a plant and in two animals are discussed in terms of the general mechanisms of gene family evolution.

The Arabidopsis thaliana genome contains at least 29 active genes encoding SET domain proteins that can be assigned to four evolutionarily conserved classes

SET domains are conserved amino acid motifs present in chromosomal proteins that function in epigenetic control of gene expression. These proteins can be divided into four classes as typified by their Drosophila members E(Z), TRX, ASH1 and SU(VAR)3-9. Homologs of all four classes have been identified in yeast and mammals, but not in plants. A BLASTP screening of the Arabidopsis genome identified 37 genes: three E(z) homologs, five trx homologs, four ash1 homologs and 15 genes similar to Su(var)3-9. Seven genes were assigned as trx-related and three as ash1-related. Only four genes have been described previously. Our classification is based on the characteristics of the SET domains, cysteine-rich regions and additional conserved domains, including a novel YGD domain. RT-PCR analysis, cDNA cloning and matching ESTs show that at least 29 of the genes are active in diverse tissues. The high number of SET domain genes, possibly involved in epigenetic control of gene activity during plant development, can partly be explained by extensive genome duplication in Arabidopsis. Additionally, the lack of introns in the coding region of eight SU(VAR)3-9 class genes indicates evolution of new genes by retrotransposition. The identification of putative nuclear localization signals and AT-hooks in many of the proteins supports an anticipated nuclear localization, which was demonstrated for selected proteins.

Low rates of silent substitution in nuclear genes of two distantly related Scrophulariaceae (Antirrhinum and Verbascum)

Low levels of genetic diversity and divergence at nuclear loci have previously been observed for cycloidea and fil1-like genes within and between several Antirrhinum species, and divergence at these loci is also low between species in genera at different levels of relatedness in the former family Scrophulariaceae (Digitalis and Verbascum). The low divergence values are surprising, because (based on the sequences of chloroplast loci) the Scrophulariaceae are thought to be polyphyletic, with two anciently diverged clades, and the species we compared belonged to the two different clades. Here, we extend our studies of sequence divergence to more nuclear genes: fil2, far, globosa, and ADH: Detailed studies revealed that in Antirrhinum these genes belong to gene families. Low levels of divergence between Antirrhinum and Verbascum were observed for four of the loci studied, fil2-1, fil2-2, far-L, and globosa, similar to our previous observations. We discuss hypotheses to explain these low synonymous divergence values. For Adh, no cases of very similar sequences were found, but, rather, our sequences from the three different genera (Antirrhinum, Digitalis, and Verbascum) were all very diverged. Repeated gene duplication and loss of elements in the Adh gene family is likely in these lineages, making it impossible to determine orthology of the Adh genes.

Genomic organization of plant terpene synthases and molecular evolutionary implications

Terpenoids are the largest, most diverse class of plant natural products and they play numerous functional roles in primary metabolism and in ecological interactions. The first committed step in the formation of the various terpenoid classes is the transformation of the prenyl diphosphate precursors, geranyl diphosphate, farnesyl diphosphate, and geranylgeranyl diphosphate, to the parent structures of each type catalyzed by the respective monoterpene (C(10)), sesquiterpene (C(15)), and diterpene synthases (C(20)). Over 30 cDNAs encoding plant terpenoid synthases involved in primary and secondary metabolism have been cloned and characterized. Here we describe the isolation and analysis of six genomic clones encoding terpene synthases of conifers, [(-)-pinene (C(10)), (-)-limonene (C(10)), (E)-alpha-bisabolene (C(15)), delta-selinene (C(15)), and abietadiene synthase (C(20)) from Abies grandis and taxadiene synthase (C(20)) from Taxus brevifolia], all of which are involved in natural products biosynthesis. Genome organization (intron number, size, placement and phase, and exon size) of these gymnosperm terpene synthases was compared to eight previously characterized angiosperm terpene synthase genes and to six putative terpene synthase genomic sequences from Arabidopsis thaliana. Three distinct classes of terpene synthase genes were discerned, from which assumed patterns of sequential intron loss and the loss of an unusual internal sequence element suggest that the ancestral terpenoid synthase gene resembled a contemporary conifer diterpene synthase gene in containing at least 12 introns and 13 exons of conserved size. A model presented for the evolutionary history of plant terpene synthases suggests that this superfamily of genes responsible for natural products biosynthesis derived from terpene synthase genes involved in primary metabolism by duplication and divergence in structural and functional specialization. This novel molecular evolutionary approach focused on genes of secondary metabolism may have broad implications for the origins of natural products and for plant phylogenetics in general.