Comparative genomics and regulatory evolution: conservation and function of the Chs and Apetala3 promoters

Transcriptomics of temperature-sensitive R gene-mediated resistance identifies a WAKL10 protein interaction network

Understanding temperature-sensitivity of R gene-mediated resistance against apoplastic pathogens is important for sustainable food production in the face of global warming. Here, we show that resistance of Brassica napus cotyledons against Leptosphaeria maculans was temperature-sensitive in introgression line Topas-Rlm7 but temperature-resilient in Topas-Rlm4. A set of 1,646 host genes was differentially expressed in Topas-Rlm4 and Topas-Rlm7 in response to temperature. Amongst these were three WAKL10 genes, including BnaA07g20220D, representing the temperature-sensitive Rlm7-1 allele and Rlm4. Network analysis identified a WAKL10 protein interaction cluster specifically for Topas-Rlm7 at 25 °C. Diffusion analysis of the Topas-Rlm4 network identified WRKY22 as a putative regulatory target of the ESCRT-III complex-associated protein VPS60.1, which belongs to the WAKL10 protein interaction community. Combined enrichment analysis of gene ontology terms considering gene expression and network data linked vesicle-mediated transport to defence. Thus, dysregulation of effector-triggered defence in Topas-Rlm7 disrupts vesicle-associated resistance against the apoplastic pathogen L. maculans.

Chalcone synthases (CHSs): the symbolic type III polyketide synthases

Present review provides a thorough insight on some significant aspects of CHSs over a period of about past three decades with a better outlook for future studies toward comprehending the structural and mechanistic intricacy of this symbolic enzyme. Polyketide synthases (PKSs) form a large family of iteratively acting multifunctional proteins that are involved in the biosynthesis of spectrum of natural products. They exhibit remarkable versatility in the structural configuration and functional organization with an incredible ability to generate different classes of compounds other than the characteristic secondary metabolite constituents. Architecturally, chalcone synthase (CHS) is considered to be the simplest representative of Type III PKSs. The enzyme is pivotal for phenylpropanoid biosynthesis and is also well known for catalyzing the initial step of the flavonoid/isoflavonoid pathway. Being the first Type III enzyme to be discovered, CHS has been subjected to ample investigations which, to a greater extent, have tried to understand its structural complexity and promiscuous functional behavior. In this context, we vehemently tried to collect the fragmented information entirely focussed on this symbolic enzyme from about past three-four decades. The aim of this review is to selectively summarize data on some of the fundamental aspects of CHSs viz, its history and distribution, localization, structure and analogs in non-plant hosts, promoter analyses, and role in defense, with an emphasis on mechanistic studies in different species and vis-à-vis mutation-led changes, and evolutionary significance which has been discussed in detail. The present review gives an insight with a better perspective for the scientific community for future studies devoted towards delimiting the mechanistic and structural basis of polyketide biosynthetic machinery vis-à-vis CHS.

A quartet pollen phenotype identified in a population of Brassica interspecific hybrids shows incomplete penetrance and variable response to temperature

Quartet pollen, where pollen grains remain attached to each other post-meiosis, is useful for tetrad analysis, crossover assessment and centromere mapping. We observed the quartet pollen phenotype for the first time in the agriculturally significant Brassica genus, in an experimental population of allohexaploid Brassica hybrids derived from the cross (Brassica napus × B. carinata) × B. juncea followed by two self-pollination generations. Quartet pollen production was assessed in 144 genotypes under glasshouse conditions, following which a set of 16 genotypes were selected to further investigate the effect of environment (warm: 25 °C and cold: 10 °C temperatures) on quartet pollen production in growth cabinets. Under glasshouse phenotyping conditions, only 92 out of 144 genotypes produced enough pollen to score: of these, 30 did not produce any observable quartet pollen, while 62 genotypes produced quartet pollen at varying frequencies. Quartet pollen production appeared quantitative and did not clearly fall into phenotypic or qualitative categories indicative of major gene expression. No consistent effect of temperature on quartet pollen production was identified, with some genotypes producing more and some producing less quartet pollen under different temperature treatments. The genetic heterogeneity and frequent pollen infertility of this population prevents strong conclusions being made. However, it is clear that the quartet phenotype in this Brassica population does not show complete penetrance and shows variable (likely genotype-specific) response to temperature stress. In future, identification of quartet phenotypes in Brassica would perhaps best be carried out via screening of diploid (e.g. B. rapa) TILLING populations.

Comparative genomic and transcriptomic analysis of selected fatty acid biosynthesis genes and CNL disease resistance genes in oil palm

Comparative genomics and transcriptomic analyses were performed on two agronomically important groups of genes from oil palm versus other major crop species and the model organism, Arabidopsis thaliana. The first analysis was of two gene families with key roles in regulation of oil quality and in particular the accumulation of oleic acid, namely stearoyl ACP desaturases (SAD) and acyl-acyl carrier protein (ACP) thioesterases (FAT). In both cases, these were found to be large gene families with complex expression profiles across a wide range of tissue types and developmental stages. The detailed classification of the oil palm SAD and FAT genes has enabled the updating of the latest version of the oil palm gene model. The second analysis focused on disease resistance (R) genes in order to elucidate possible candidates for breeding of pathogen tolerance/resistance. Ortholog analysis showed that 141 out of the 210 putative oil palm R genes had homologs in banana and rice. These genes formed 37 clusters with 634 orthologous genes. Classification of the 141 oil palm R genes showed that the genes belong to the Kinase (7), CNL (95), MLO-like (8), RLK (3) and Others (28) categories. The CNL R genes formed eight clusters. Expression data for selected R genes also identified potential candidates for breeding of disease resistance traits. Furthermore, these findings can provide information about the species evolution as well as the identification of agronomically important genes in oil palm and other major crops.

Regulation of Fig (Ficus carica L.) Fruit Color: Metabolomic and Transcriptomic Analyses of the Flavonoid Biosynthetic Pathway

Combined metabolomic and transcriptomic analyses were carried out with fig cultivar Green Peel and its color mutant "Purple Peel." Five and twenty-two metabolites were identified as having significantly different contents between fruit peels of the two cultivars at young and mature stages, respectively. Cyanidin O-malonylhexoside demonstrated a 3,992-fold increase in the mature purple peel, the first identification of a major cyanidin in fig fruit; cyanidin 3-O-glucoside, cyanidin O-malonylhexoside O-hexoside and cyanidin-3,5-O-diglucoside were upregulated 100-fold, revealing the anthocyanins underlying the purple mutation. Beyond the visible differences, there was very significant accumulation of the colorless flavonoids procyanidin B1, luteolin-3',7-di-O-glucoside, epicatechin and quercetin-3-O-rhamnoside in the mature "Purple Peel" compared to "Green Peel." At the young stage, only cyanidin O-malonylhexoside, cyanidin O-malonylhexoside O-hexoside and esculetin were upregulated a few fold in the mutant. Transcriptome analysis revealed a downregulated expression trend of genes encoding phenylpropanoid and flavonoid biosynthetic pathway enzyme in the young "Purple Peel" compared to the young "Green Peel," whereas significant and simultaneous upregulation was revealed in almost all of the flavonoid and anthocyanin pathway components and relevant transcription factors in the mature-stage mutant. The role of R2R3-MYB transcription factors in the color morph mutation and its possible relation to the activity of retrotransposons are discussed. Moreover, large-scale upregulation of small heat-shock protein genes was found in the mature mutant. This is the first work to reveal comprehensive metabolome and transcriptome network changes underlying a fig mutation in a single horticultural attribute, and its profound effects on fruit nutrition and quality.

Analysis of cis-acting regulatory elements of Respiratory burst oxidase homolog (Rboh) gene families in Arabidopsis and rice provides clues for their diverse functions

NADPH oxidase (NOX) is a critical enzyme in the production of reactive oxygen species (ROS). It catalyzes the production of apoplastic superoxide (O2(-)), that regulates a wide array of biological functions in different organisms. Plant Noxes are homologs of catalytic subunit of mammalian NADPH oxidase and are well-known as Respiratory burst oxidase homologs (Rbohs). In recent years, there has been growing interest to study plant Noxes due to their versatile roles in plant systems. In the present work, comprehensive analysis on upstream regions from 10 Rbohs from Arabidopsis thaliana and 9 from Oryza sativa japonica was conducted. The distribution of various cis-elements, CpG islands and tandem repeats were analyzed to uncover the 5' regulatory region in wide array of functions from Rbohs. Information retrieved from cis-elements analysis was also correlated with the microarray data. Present study which involves uncovering transcription regulatory elements provided vital clues for diverse functions of plant Rbohs.

Genome-wide conserved non-coding microsatellite (CNMS) marker-based integrative genetical genomics for quantitative dissection of seed weight in chickpea

Phylogenetic footprinting identified 666 genome-wide paralogous and orthologous CNMS (conserved non-coding microsatellite) markers from 5'-untranslated and regulatory regions (URRs) of 603 protein-coding chickpea genes. The (CT)n and (GA)n CNMS carrying CTRMCAMV35S and GAGA8BKN3 regulatory elements, respectively, are abundant in the chickpea genome. The mapped genic CNMS markers with robust amplification efficiencies (94.7%) detected higher intraspecific polymorphic potential (37.6%) among genotypes, implying their immense utility in chickpea breeding and genetic analyses. Seventeen differentially expressed CNMS marker-associated genes showing strong preferential and seed tissue/developmental stage-specific expression in contrasting genotypes were selected to narrow down the gene targets underlying seed weight quantitative trait loci (QTLs)/eQTLs (expression QTLs) through integrative genetical genomics. The integration of transcript profiling with seed weight QTL/eQTL mapping, molecular haplotyping, and association analyses identified potential molecular tags (GAGA8BKN3 and RAV1AAT regulatory elements and alleles/haplotypes) in the LOB-domain-containing protein- and KANADI protein-encoding transcription factor genes controlling the cis-regulated expression for seed weight in the chickpea. This emphasizes the potential of CNMS marker-based integrative genetical genomics for the quantitative genetic dissection of complex seed weight in chickpea.

Evolutionary conservation of cold-induced antisense RNAs of FLOWERING LOCUS C in Arabidopsis thaliana perennial relatives

Antisense RNA (asRNA) COOLAIR is expressed at A. thaliana FLOWERING LOCUS C (FLC) in response to winter temperatures. Its contribution to cold-induced silencing of FLC was proposed but its functional and evolutionary significance remain unclear. Here we identify a highly conserved block containing the COOLAIR first exon and core promoter at the 3′ end of several FLC orthologues. Furthermore, asRNAs related to COOLAIR are expressed at FLC loci in the perennials A. alpina and A. lyrata, although some splicing variants differ from A. thaliana. Study of the A. alpina orthologue, PERPETUAL FLOWERING 1 (PEP1), demonstrates that AaCOOLAIR is induced each winter of the perennial life cycle. Introduction of PEP1 into A. thaliana reveals that AaCOOLAIR cis-elements confer cold-inducibility in this heterologous species while the difference between PEP1 and FLC mRNA patterns depends on both cis-elements and species-specific trans-acting factors. Thus, expression of COOLAIR is highly conserved, supporting its importance in FLC regulation.

FLOWERING LOCUS C (FLC) is thought to control the flowering time of A. thaliana in response to winter temperatures, in a process known as vernalization. Here, the authors suggest that the COOLAIR antisense RNA, which is conserved across plant species, acts to repress the expression of FLC during vernalization.

The AtCathB3 gene, encoding a cathepsin B-like protease, is expressed during germination of Arabidopsis thaliana and transcriptionally repressed by the basic leucine zipper protein GBF1

Protein hydrolysis plays an important role during seed germination and post-germination seedling establishment. In Arabidopsis thaliana, cathepsin B-like proteases are encoded by a gene family of three members, but only the AtCathB3 gene is highly induced upon seed germination and at the early post-germination stage. Seeds of a homozygous T-DNA insertion mutant in the AtCathB3 gene have, besides a reduced cathepsin B activity, a slower germination than the wild type. To explore the transcriptional regulation of this gene, we used a combined phylogenetic shadowing approach together with a yeast one-hybrid screening of an arrayed library of approximately 1200 transcription factor open reading frames from Arabidopsis thaliana. We identified a conserved CathB3-element in the promoters of orthologous CathB3 genes within the Brassicaceae species analysed, and, as its DNA-interacting protein, the G-Box Binding Factor1 (GBF1). Transient overexpression of GBF1 together with a PAtCathB3::uidA (β-glucuronidase) construct in tobacco plants revealed a negative effect of GBF1 on expression driven by the AtCathB3 promoter. In stable P35S::GBF1 lines, not only was the expression of the AtCathB3 gene drastically reduced, but a significant slower germination was also observed. In the homozygous knockout mutant for the GBF1 gene, the opposite effect was found. These data indicate that GBF1 is a transcriptional repressor of the AtCathB3 gene and affects the germination kinetics of Arabidopsis thaliana seeds. As AtCathB3 is also expressed during post-germination in the cotyledons, a role for the AtCathB3-like protease in reserve mobilization is also inferred.

Analysis of the formation of flower shapes in wild species and cultivars of tree peony using the MADS-box subfamily gene

Tree peony (Paeonia suffricotisa) cultivars have a unique character compared with wild species; the stamen petalody results in increased whorls of petals and generates different flower forms, which are one of the most important traits for cultivar classification. In order to investigate how petaloid stamens are formed, we obtained the coding sequence (666 bp) and genomic DNA sequence of the PsTM6 genes (belongs to B subfamily of MADS-box gene family) from 23 tree peony samples, Five introns and six exons consisted of the genomic DNA sequence. The analysis of cis-acting regulatory elements in the third and fourth intron indicated that they were highly conserved in all samples. Partial putative amino acids were analyzed and the results suggested that functional differentiation of PsTM6 paralogs apparently affected stamen petalody and flower shape formation due to due to amino acid substitution caused by differences in polarity and electronic charge. Sliding window analysis indicated that the different regions of PsTM6 were subjected to different selection forces, especially in the K domain. This is the first attempt to investigate genetic control of the stamen petalody based on the PsTM6 sequence. This will provide a basis for understanding the evolution of PsTM6 and its the function of in determining stamen morphology of tree peony.

Boechera, a model system for ecological genomics

The selection and development of a study system for evolutionary and ecological functional genomics (EEFG) depend on a variety of factors. Here, we present the genus Boechera as an exemplary system with which to address ecological and evolutionary questions. Our focus on Boechera is based on several characteristics as follows: (i) native populations in undisturbed habitats where current environments reflect historical conditions over several thousand years; (ii) functional genomics benefitting from its close relationship to Arabidopsis thaliana; (iii) inbreeding tolerance enabling development of recombinant inbred lines, near-isogenic lines and positional cloning; (iv) interspecific crosses permitting mapping for genetic analysis of speciation; (v) apomixis (asexual reproduction by seeds) in a genetically tractable diploid; and (vi) broad geographic distribution in North America, permitting ecological genetics for a large research community. These characteristics, along with the current sequencing of three Boechera species by the Joint Genome Institute, position Boechera as a rapidly advancing system for EEFG studies.

Robustness and evolvability in the B-system of flower development

BACKGROUND

Gene duplication has often been invoked as a key mechanism responsible for evolution of new morphologies. The floral homeotic B-group gene family has undergone a number of gene duplication events, and yet the functions of these genes appear to be largely conserved. However, detailed comparative analysis has indicated that such duplicate genes have considerable cryptic variability in their functions. In the Solanaceae, two duplicate B-group gene lineages have been retained in three subfamilies. Comparisons of orthologous genes across members of the Solanaceae have demonstrated that the combined function of all four B-gene members is to establish petal and stamen identity, but that this function was partitioned differently in each species. These observations emphasize both the robustness and the evolvability of the B-system.

SCOPE

We provide an overview of how the B-function genes can robustly specify petal and stamen identity and at the same time evolve through changes in protein-protein interaction, gene expression patterns, copy number variation or alterations in the downstream genes they control. By using mathematical models we explore regulatory differences between species and how these impose constraints on downstream gene regulation.

CONCLUSIONS

Evolvability of the B-genes can be understood through the multiple ways in which the B-system can be robust. Quantitative approaches should allow for the incorporation of more biological realism in the representations of these regulatory systems and this should contribute to understanding the constraints under which different B-systems can function and evolve. This, in turn, can provide a better understanding of the ways in which B-genes have contributed to flower diversity.

Arctic mustard flower color polymorphism controlled by petal-specific downregulation at the threshold of the anthocyanin biosynthetic pathway

Intra- and interspecific variation in flower color is a hallmark of angiosperm diversity. The evolutionary forces underlying the variety of flower colors can be nearly as diverse as the colors themselves. In addition to pollinator preferences, non-pollinator agents of selection can have a major influence on the evolution of flower color polymorphisms, especially when the pigments in question are also expressed in vegetative tissues. In such cases, identifying the target(s) of selection starts with determining the biochemical and molecular basis for the flower color variation and examining any pleiotropic effects manifested in vegetative tissues. Herein, we describe a widespread purple-white flower color polymorphism in the mustard Parrya nudicaulis spanning Alaska. The frequency of white-flowered individuals increases with increasing growing-season temperature, consistent with the role of anthocyanin pigments in stress tolerance. White petals fail to produce the stress responsive flavonoid intermediates in the anthocyanin biosynthetic pathway (ABP), suggesting an early pathway blockage. Petal cDNA sequences did not reveal blockages in any of the eight enzyme-coding genes in white-flowered individuals, nor any color differentiating SNPs. A qRT-PCR analysis of white petals identified a 24-fold reduction in chalcone synthase (CHS) at the threshold of the ABP, but no change in CHS expression in leaves and sepals. This arctic species has avoided the deleterious effects associated with the loss of flavonoid intermediates in vegetative tissues by decoupling CHS expression in petals and leaves, yet the correlation of flower color and climate suggests that the loss of flavonoids in the petals alone may affect the tolerance of white-flowered individuals to colder environments.

Isolation and phylogenetic footprinting analysis of the 5'-regulatory region of the floral homeotic gene OrcPI from Orchis italica (Orchidaceae)

The nucleotide sequences of regulatory elements from homologous genes can be strongly divergent. Phylogenetic footprinting, a comparative analysis of noncoding regions, can detect putative transcription factor binding sites (TFBSs) shared among the regulatory regions of 2 or more homologous genes. These conserved motifs have the potential to serve the same regulatory function in distantly related taxa. We isolated the 5'-noncoding region of the OrcPI gene, a MADS-box transcription factor involved in flower development in Orchis italica, using the thermal asymmetric interlaced polymerase chain reaction technique. This region (comprising 1352 bp) induced transient beta-glucuronidase expression in the petal tissue of white Rosa hybrida flowers and represents the 5'-regulatory sequence of the OrcPI gene. Phylogenetic footprinting analysis detected conserved regions within the 5'-regulatory sequence of OrcPI and the homologous regions of Oryza sativa, Lilium regale, and Arabidopsis thaliana. Some of these sequences are known TFBSs described in databases of plant regulatory elements. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the following accession numbers: AF198055 promoter region of the PISTILLATA (PI) gene of A. thaliana; AB094985 cDNA of OrcPI (PI/GLOBOSA [PI/GLO] homologue) of O. italica; AB378089 5'-regulatory region of the OrcPI gene of O. italica; AP008211 putative promoter region of OSMADS2 (PI/GLO homologue) of O. sativa; AP008207 putative promoter region of OSMADS4 (PI/GLO homologue) of O. sativa; and AB158292 putative promoter region of the PI/GLO homologue of L. regale.

Genomics of fungal disease resistance in tomato

Tomato (Solanum lycopersicum) is an important vegetable crop worldwide. Often times, its production is hindered by fungal diseases. Important fungal diseases limiting tomato production are late blight, caused by Phytophthora infestans, early blight, caused by Alternaria solanii, and septoria leaf spot, caused by Septoria lycopersici, fusarium wilt caused by Fusarium oxysporium fsp. oxysporium, and verticilium wilt caused by Verticilium dahlea. The Phytophthora infestans is the same fungus that caused the devastating loss of potato in Europe in 1845. A similar magnitude of crop loss in tomato has not occurred but Phytophthora infestans has caused the complete loss of tomato crops around the world on a small scale. Several attempts have been made through conventional breeding and the molecular biological approaches to understand the biology of host-pathogen interaction so that the disease can be managed and crop loss prevented. In this review, we present a comprehensive analysis of information produced by molecular genetic and genomic experiments on host-pathogen interactions of late blight, early blight, septoria leaf spot, verticilim wilt and fusarium wilt in tomato. Furthermore, approaches adopted to manage these diseases in tomato including genetic transformation are presented. Attempts made to link molecular markers with putative genes and their use in crop improvement are discussed.

Phylogenetic and transcriptional analysis of a strictosidine synthase-like gene family in Arabidopsis thaliana reveals involvement in plant defence responses

Protein domains with similarity to plant strictosidine synthase-like (SSL) sequences have been uncovered in the genomes of all multicellular organisms sequenced so far and are known to play a role in animal immune responses. Among several distinct groups of Arabidopsis thaliana SSL sequences, four genes (AtSSL4-AtSSL7) arranged in tandem on chromosome 3 show more similarity to SSL genes from Drosophila melanogaster and Caenorhabditis elegans than to other Arabidopsis SSL genes. To examine whether any of the four AtSSL genes are immune-inducible, we analysed the expression of each of the four AtSSL genes after exposure to microbial pathogens, wounding and plant defence elicitors using real-time quantitative RT-PCR, Northern blot hybridisation and Western blot analysis with antibodies raised against recombinant AtSSL proteins. While the AtSSL4 gene was constitutively expressed and not significantly induced by any treatment, the other three AtSSL genes were induced to various degrees by plant defence signalling compounds, such as salicylic acid, methyl jasmonate and ethylene, as well as by wounding and exposure to the plant pathogens Alternaria brassicicola and cucumber mosaic virus. Our data demonstrate that the four SSL-coding genes are regulated individually, suggesting specific roles in basal (SSL4) and inducible (SSL5-7) plant defence mechanisms.

Phylogenetic footprints in fern spore- and seed-specific gene promoters

Spermatophyte seed-storage proteins have descended from a group of proteins involved in cellular desiccation/hydration processes. Conserved protein structures are found across all plant phyla and in the fungi and Archaea. We investigated whether conservation in the coding region sequence is paralleled by common gene regulatory processes. Seed- and spore-specific gene promoters of three phylogenetically diverse plants were analysed by transient and transgenic expression in Arabidopsis thaliana and tobacco. The transcription factors FUS3 and ABI3, which are central regulators of seed maturation processes, interact with cis-motifs of seed-specific promoters from distantly related plants. The promoter of a fern spore-specific gene encoding a seed-storage globulin-like protein exhibits strong seed-specific activity in both Arabidopsis and tobacco. The existence of phylogenetic footprints indicates good conservation of regulatory pathways controlling gene expression in fern spores and in gymnosperm and angiosperm seeds, reflecting the concerted evolution of coding and regulatory regions.

Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera

As a model system, Arabidopsis thaliana and its wild relatives have played an important role in the study of genomics and evolution in plants. In this study, we examined the genetic diversity of the chalcone synthase (Chs) gene, which encodes a key enzyme of the flavonoid pathway and is located on chromosome five, as well as two Chs-like genes on the first and fourth chromosomes of Arabidopsis. The objectives of the study are to determine if natural selection operates differentially on the paralogs of the Chs gene family in A. thaliana and Arabidopsis halleri ssp. gemmifera. The mode of selection was inferred from Tajima's D values from noncoding and coding regions, as well as from the ratio of nonsynonymous to synonymous substitutions. Both McDonald-Kreitman and HKA tests revealed the effects of selection on the allelic distribution, except for the chromosome 1 paralog in ssp. gemmifera. The Chs gene on chromosome 5 was under purifying selection in both species. Significant, negative Tajima's D values at synonymous sites and positive Fay and Wu's H values within coding region, plus reduced genetic variability in introns, indicated effects of background selection in shaping the evolution of this gene region in A. thaliana. The Chs paralog on chromosome 1 was under positive selection in A. thaliana, while interspecific introgression and balancing selection determined the fates of the paralog and resulted in high heterogeneity in ssp. gemmifera. Local adaptation differentiated populations of Japan and China at the locus. In contrast, the other Chs-paralog of chromosome 4 was shaped by purifying selection in A. thaliana, while under positive selection in ssp. gemmifera, as indicated by dn/ds>1. Moreover, these contrasting patterns of selection have likely resulted in functional divergence in Arabidopsis, as indicated by radical amino acid substitutions at the chalcone synthase/stilbene synthase motif of the Chs genes. Unlike previous studies of the evolutionary history of A. thaliana, the high levels of genetic diversity in most gene regions of Chs paralogs and nonsignificant Tajima's D in the intron sequences of the Chs gene family in A. thaliana did not reflect the effects of a recent demographic expansion.

Evolutionary hierarchies of conserved blocks in 5'-noncoding sequences of dicot rbcS genes

Background

Evolutionary processes in gene regulatory regions are major determinants of organismal evolution, but exceptionally challenging to study. We explored the possibilities of evolutionary analysis of phylogenetic footprints in 5'-noncoding sequences (NCS) from 27 ribulose-1,5-bisphosphate carboxylase small subunit (rbcS) genes, from three dicot families (Brassicaceae, Fabaceae and Solanaceae).

Results

Sequences of up to 400 bp encompassing proximal promoter and 5'-untranslated regions were analyzed. We conducted phylogenetic footprinting by several alternative methods: generalized Lempel-Ziv complexity (C_LZ), multiple alignments with DIALIGN and ALIGN-M, and the MOTIF SAMPLER Gibbs sampling algorithm. These tools collectively defined 36 conserved blocks of mean length 12.8 bp. On average, 12.5 blocks were found in each 5'-NCS. The blocks occurred in arrays whose relative order was absolutely conserved, confirming the existence of 'conserved modular arrays' in promoters. Identities of half of the blocks confirmed past rbcS research, including versions of the I-box, G-box, and GT-1 sites such as Box II. Over 90% of blocks overlapped DNase-protected regions in tomato 5'-NCS. Regions characterized by low C_LZ in sliding-window analyses were also frequently associated with DNase-protection. Blocks could be assigned to evolutionary hierarchies based on taxonomic distribution and estimated age. Lineage divergence dates implied that 13 blocks found in all three plant families were of Cretaceous antiquity, while other family-specific blocks were much younger. Blocks were also dated by formation of multigene families, using genome and coding sequence information. Dendrograms of evolutionary relations of the 5'-NCS were produced by several methods, including: cluster analysis using pairwise C_LZ values; evolutionary trees of DIALIGN sequence alignments; and cladistic analysis of conserved blocks.

Conclusion

Dicot 5'-NCS contain conserved modular arrays of recurrent sequence blocks, which are coincident with functional elements. These blocks are amenable to evolutionary interpretation as hierarchies in which ancient, taxonomically widespread blocks can be distinguished from more recent, taxon-specific ones.

Cis-regulatory evolution of chalcone-synthase expression in the genus Arabidopsis

The contribution of cis-regulation to adaptive evolutionary change is believed to be essential, yet little is known about the evolutionary rules that govern regulatory sequences. Here, we characterize the short-term evolutionary dynamics of a cis-regulatory region within and among two closely related species, A. lyrata and A. halleri, and compare our findings to A. thaliana. We focused on the cis-regulatory region of chalcone synthase (CHS), a key enzyme involved in the synthesis of plant secondary metabolites. We observed patterns of nucleotide diversity that differ among species but do not depart from neutral expectations. Using intra- and interspecific F1 progeny, we have evaluated functional cis-regulatory variation in response to light and herbivory, environmental cues, which are known to induce CHS expression. We find that substantial cis-regulatory variation segregates within and among populations as well as between species, some of which results from interspecific genetic introgression. We further demonstrate that, in A. thaliana, CHS cis-regulation in response to herbivory is greater than in A. lyrata or A. halleri. Our work indicates that the evolutionary dynamics of a cis-regulatory region is characterized by pervasive functional variation, achieved mostly by modification of response modules to one but not all environmental cues. Our study did not detect the footprint of selection on this variation.

Toward a Global Phylogeny of the Brassicaceae

The Brassicaceae is a large plant family (338 genera and 3,700 species) of major scientific and economic importance. The taxonomy of this group has been plagued by convergent evolution in nearly every morphological feature used to define tribes and genera. Phylogenetic analysis of 746 nrDNA internal transcribed spacer (ITS) sequences, representing 24 of the 25 currently recognized tribes, 146 genera, and 461 species of Brassicaceae, produced the most comprehensive, single-locus-based phylogenetic analysis of the family published to date. Novel approaches to nrDNA ITS analysis and extensive taxonomic sampling offered a test of monophyly for a large complement of the currently recognized tribes and genera of Brassicaceae. In the most comprehensive analysis, tribes Alysseae, Anchonieae plus Hesperideae, Boechereae, Cardamineae, Eutremeae, Halimolobeae, Iberideae, Noccaeeae, Physarieae, Schizopetaleae, Smelowskieae, and Thlaspideae were all monophyletic. Several broadly defined genera (e.g., Draba and Smelowskia) were supported as monophyletic, whereas others (e.g., Sisymbrium and Alyssum) were clearly polyphyletic. Analyses of ITS data identified several problematic sequences attributable to errors in sample identification or database submission. Results from parsimony ratchet and Bayesian analyses recovered little support for the backbone of the phylogeny, suggesting that many lineages of Brassicaceae have undergone rapid radiations that may ultimately be difficult to resolve with any single locus. However, the development of a preliminary supermatrix including the combination of 10 loci for 65 species provides an initial estimate of intertribal relations and suggests that broad application of such a method will provide greater understanding of relationships in the family.

Analysis of the Petunia TM6 MADS box gene reveals functional divergence within the DEF/AP3 lineage

Antirrhinum majus DEFICIENS (DEF) and Arabidopsis thaliana APETALA3 (AP3) MADS box proteins are required to specify petal and stamen identity. Sampling of DEF/AP3 homologs revealed two types of DEF/AP3 proteins, euAP3 and TOMATO MADS BOX GENE6 (TM6), within core eudicots, and we show functional divergence in Petunia hybrida euAP3 and TM6 proteins. Petunia DEF (also known as GREEN PETALS [GP]) is expressed mainly in whorls 2 and 3, and its expression pattern remains unchanged in a blind (bl) mutant background, in which the cadastral C-repression function in the perianth is impaired. Petunia TM6 functions as a B-class organ identity protein only in the determination of stamen identity. Atypically, Petunia TM6 is regulated like a C-class rather than a B-class gene, is expressed mainly in whorls 3 and 4, and is repressed by BL in the perianth, thereby preventing involvement in petal development. A promoter comparison between DEF and TM6 indicates an important change in regulatory elements during or after the duplication that resulted in euAP3- and TM6-type genes. Surprisingly, although TM6 normally is not involved in petal development, 35S-driven TM6 expression can restore petal development in a def (gp) mutant background. Finally, we isolated both euAP3 and TM6 genes from seven solanaceous species, suggesting that a dual euAP3/TM6 B-function system might be the rule in the Solanaceae.

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.

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.

Allopolyploid origin of Cardamine asarifolia (Brassicaceae): incongruence between plastid and nuclear ribosomal DNA sequences solved by a single-copy nuclear gene

Interspecific hybridization and polyploidization have played central roles in plant diversification. However, technical difficulties in the analyses of low-copy genes have limited the study of the origins of hybrid and polyploid plants. Here, we present a phylogenetic analysis of the hexaploid Cardamine asarifolia, distributed in the southern European Alps and northern Apennines. Our study included all relevant taxa of the genus found in Europe. A marked discrepancy was revealed between the trnL-trnF region of cpDNA and internal transcribed spacer (nrDNA ITS) sequences. To solve the incongruence, we sequenced a single-copy nuclear CHS gene (chalcone synthase) using a novel method to design homoeologue-specific PCR primers to bypass artefacts caused by artificial recombination of homoeologues during PCR and/or cloning. Three homoeologues were isolated from C. asarifolia, providing evidence for its allopolyploid origin. One homoeologue, showing the same phylogenetic position as the ITS sequences, most likely originated from an extinct parent. Furthermore, we documented recurrent polytopic hybridizations between C. asarifolia and diploid C. amara. The allohexaploidization and the following hybridization with a diploid species exemplify the ongoing dynamic processes of speciation in the genus Cardamine.

Variation in conserved non-coding sequences on chromosome 5q and susceptibility to asthma and atopy

BACKGROUND

Evolutionarily conserved sequences likely have biological function.

METHODS

To determine whether variation in conserved sequences in non-coding DNA contributes to risk for human disease, we studied six conserved non-coding elements in the Th2 cytokine cluster on human chromosome 5q31 in a large Hutterite pedigree and in samples of outbred European American and African American asthma cases and controls.

RESULTS

Among six conserved non-coding elements (> 100 bp, > 70% identity; human-mouse comparison), we identified one single nucleotide polymorphism (SNP) in each of two conserved elements and six SNPs in the flanking regions of three conserved elements. We genotyped our samples for four of these SNPs and an additional three SNPs each in the IL13 and IL4 genes. While there was only modest evidence for association with single SNPs in the Hutterite and European American samples (P < 0.05), there were highly significant associations in European Americans between asthma and haplotypes comprised of SNPs in the IL4 gene (P < 0.001), including a SNP in a conserved non-coding element. Furthermore, variation in the IL13 gene was strongly associated with total IgE (P = 0.00022) and allergic sensitization to mold allergens (P = 0.00076) in the Hutterites, and more modestly associated with sensitization to molds in the European Americans and African Americans (P < 0.01).

CONCLUSION

These results indicate that there is overall little variation in the conserved non-coding elements on 5q31, but variation in IL4 and IL13, including possibly one SNP in a conserved element, influence asthma and atopic phenotypes in diverse populations.

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.

Allele-specific assay reveals functional variation in the chalcone synthase promoter of Arabidopsis thaliana that is compatible with neutral evolution

Promoters are thought to play a major role in adaptive evolution, yet little is known about the regulatory diversity within species, where microevolutionary processes take place. To investigate the potential for evolutionary change in the promoter of a gene, we examined nucleotide and functional variation of the Chalcone Synthase (CHS) cis-regulatory region in Arabidopsis thaliana. CHS is the branch point enzyme of a biosynthetic pathway that leads to the production of secondary metabolites influencing the interaction between the plant and its environment. We found that nucleotide diversity in the intergenic region encompassing the CHS promoter (pi=0.003) is compatible with neutral expectations. To quantify functional variation specifically as a result of cis-regulation of CHS mRNA levels, we developed an assay using F1 individuals in which distinct promoter alleles are compared within a common trans-regulatory background. We examined functional cis-regulatory variation in response to different stimuli representing a variety of CHS transcriptional environments (dark, light, and insect feeding). We observed extensive functional variation, some of which appeared to be independent of the trans-regulatory background. Comparison of functional and nucleotide diversity suggested a candidate point mutation that may explain cis-regulatory differences in light response. Our results indicate that functional changes in promoters can arise from a few mutations, pointing to promoter regions as a fundamental determinant of functional genetic variation.

Molecular evolution of FLORICAULA/LEAFY orthologs in the Andropogoneae (Poaceae)

Members of the grass family (Poaceae) exhibit a broad range of inflorescence structures and other morphologies, making the grasses an interesting model system for studying the evolution of development. Here we present an analysis of the molecular evolution of FLORICAULA/LEAFY-like genes, which are important developmental regulatory loci known to affect inflorescence development in a wide range of flowering plant species. We have focused on sequences from the Andropogoneae, a tribe within the grass family that includes maize (Zea mays ssp. mays) and Sorghum (Sorghum bicolor). The FLORICAULA/LEAFY gene phylogeny we generated largely agrees with previously published phylogenies for the Andropogoneae using other nuclear genes but is unique in that it includes both members of one of the many duplicate gene sets present in maize. The placement of these sequences in the phylogeny suggests that the duplication of the maize FLORICAULA/LEAFY orthologs, zfl1 and zfl2, is a consequence of a proposed tetraploidy event that occurred in the common ancestor of Zea and a closely related genus, Tripsacum. Our data are consistent with the hypothesis that the transcribed regions of the FLORICAULA/LEAFY-like genes in the Andropogoneae are functionally constrained at both nonsynonymous and synonymous sites and show no evidence of directional selection. We also examined conservation of short noncoding sequences in the first intron, which may play a role in gene regulation. Finally, we investigated the genetic diversity of one of the two maize FLORICAULA/LEAFY orthologs, zfl2, in maize and its wild ancestor, teosinte (Z. mays ssp. parviglumis), and found no evidence for selection pressure resulting from maize domestication within the zfl2-coding region.

Differential combinatorial interactions of cis-acting elements recognized by R2R3-MYB, BZIP, and BHLH factors control light-responsive and tissue-specific activation of phenylpropanoid biosynthesis genes

Chalcone synthase (CHS), chalcone flavanone isomerase (CFI), flavanone 3-hydroxylase (F3H) and flavonol synthase (FLS) catalyze successive steps in the biosynthetic pathway leading to the production of flavonols. We show that in Arabidopsis thaliana all four corresponding genes are coordinately expressed in response to light, and are spatially coexpressed in siliques, flowers and leaves. Light regulatory units (LRUs) sufficient for light responsiveness were identified in all four promoters. Each unit consists of two necessary elements, namely a MYB-recognition element (MRE) and an ACGT-containing element (ACE). C1 and Sn, a R2R3-MYB and a BHLH factor, respectively, known to control tissue specific anthocyanin biosynthesis in Z. mays, were together able to activate the AtCHS promoter. This activation of the CHS promoter required an intact MRE and a newly identified sequence designated R response element (RREAtCHS) containing the BHLH factor consensus binding site CANNTG. The RRE was dispensable for light responsiveness, and the ACE was not necessary for activation by C1/Sn. These data suggest that a BHLH and a R2R3-MYB factor cooperate in directing tissue-specific production of flavonoids, while an ACE-binding factor, potentially a BZIP, and a R2R3-MYB factor work together in conferring light responsiveness.

Activation of CRABS CLAW in the Nectaries and Carpels of Arabidopsis

CRABS CLAW (CRC), a member of the YABBY gene family, is required for nectary and carpel development. To further understand CRC regulation in Arabidopsis thaliana, we performed phylogenetic footprinting analyses of 5' upstream regions of CRC orthologs from three Brassicaceae species, including Arabidopsis. Phylogenetic footprinting efficiently identified functionally important regulatory regions (modules), indicating that CRC expression is regulated by a combination of positive and negative regulatory elements in the modules. Within the conserved modules, we identified putative binding sites of LEAFY and MADS box proteins, and functional in vivo analyses revealed their importance for CRC expression. Both expression and genetic studies demonstrate that potential binding sites for MADS box proteins within the conserved regions are functionally significant for the transcriptional regulation of CRC in nectaries. We propose that in wild-type flowers, a combination of floral homeotic gene activities, specifically the B class genes APETALA3 and PISTILLATA and the C class gene AGAMOUS act redundantly with each other and in combination with SEPALLATA genes to activate CRC in the nectaries and carpels. In the absence of B and C class gene activities, other genes such as SHATTERPROOF1/2 can substitute if they are ectopically expressed, as in an A class mutant background (apetala2). These MADS box proteins may provide general floral factors that must work in conjunction with specific factors in the activation of CRC in the nectaries and carpels.

Characterization of the two arginine decarboxylase (polyamine biosynthesis) paralogues of the endemic subantarctic cruciferous species Pringlea antiscorbutica and analysis of their differential expression during development and response to environmental stress

Arginine decarboxylase (ADC) is a key enzyme involved in the synthesis of polyamines, which have been implicated in developmental processes and stress responses in higher plants. An ancestral ADC gene appears to have been duplicated at the origin of the Brassicaceae family, thus yielding two paralogues in the derived taxa. ADC gene structure was investigated in Pringlea antiscorbutica R. Br., a geographically isolated Brassicaceae species that is endemic from the subantarctic region. P. antiscorbutica exhibits several biochemical and physiological adaptations related to this cold and harsh environment, including high levels of polyamines, which is unusual in higher plants, and especially high levels of agmatine, the product of the ADC-catalysed reaction. Various ADC clones were obtained from P. antiscorbutica. Sequence and phylogenetic analysis showed that all these clones fitted with the presence of two paralogues, PaADC1 and PaADC2, in P. antiscorbutica. Expression of these ADC paralogues was analyzed in P. antiscorbutica during vegetative development and response to stress. Whereas PaADC2 was expressed at both seedling and mature stages, PaADC1 transcripts were hardly detected during early development and were significantly expressed in mature plants. Moreover, PaADC2, but not PaADC1, expression was up-regulated in response to chilling and salt stress at seedling stage. Analysis of 5' regulatory regions of these ADC genes revealed several differences in putative cis-regulatory elements, which could be associated with specific expression patterns. These results were compared to ADC paralogue expression in Arabidopsis thaliana and are discussed in the evolutionary context of genetic diversity resulting from gene duplication.

Phylogenetic analysis of 5'-noncoding regions from the ABA-responsive rab16/17 gene family of sorghum, maize and rice provides insight into the composition, organization and function of cis-regulatory modules

Phylogenetic analysis of sequences from gene families and homologous genes from species of varying divergence can be used to identify conserved noncoding regulatory elements. In this study, phylogenetic analysis of 5'-noncoding sequences was optimized using rab17, a well-characterized ABA-responsive gene from maize, and five additional rab16/17 homologs from sorghum and rice. Conserved 5'-noncoding sequences among the maize, sorghum, and rice rab16/17 homologs were identified with the aid of the software program FootPrinter and by screening for known transcription-factor-binding sites. Searches for 7 of 8 (7/8)bp sequence matches within aligned 5'-noncoding segments of the rab genes identified many of the cis-elements previously characterized by biochemical analysis in maize rab17 plus several additional putative regulatory elements. Differences in the composition of conserved noncoding sequences among rab16/17 genes were related to variation in rab gene mRNA levels in different tissues and to response to ABA treatment using qRT-PCR. Absence of a GRA-like element in the promoter of sorghum dhn2 relative to maize rab17 was correlated with an approximately 85-fold reduction of dhn2 RNA in sorghum shoots. Overall, we conclude that phylogenetic analysis of gene families among rice, sorghum, and maize will help identify regulatory sequences in the noncoding regions of genes and contribute to our understanding of grass gene regulatory networks.

cis-Regulatory and protein evolution in orthologous and duplicate genes

The relationship between protein and regulatory sequence evolution is a central question in molecular evolution. It is currently not known to what extent changes in gene expression are coupled with the evolution of protein coding sequences, or whether these changes differ among orthologs (species homologs) and paralogs (duplicate genes). Here, we develop a method to measure the extent of functionally relevant cis-regulatory sequence change in homologous genes, and validate it using microarray data and experimentally verified regulatory elements in different eukaryotic species. By comparing the genomes of Caenorhabditis elegans and C. briggsae, we found that protein and regulatory evolution is weakly coupled in orthologs but not paralogs, suggesting that selective pressure on gene expression and protein evolution is quite similar and persists for a significant amount of time following speciation but not gene duplication. Additionally, duplicates of both species exhibit a dramatic acceleration of both regulatory and protein evolution compared to orthologs, suggesting increased directional selection and/or relaxed selection on both gene expression patterns and protein function in duplicate genes.

Comparative genomics: methods and applications

Interpreting the functional content of a given genomic sequence is one of the central challenges of biology today. Perhaps the most promising approach to this problem is based on the comparative method of classic biology in the modern guise of sequence comparison. For instance, protein-coding regions tend to be conserved between species. Hence, a simple method for distinguishing a functional exon from the chance absence of stop codons is to investigate its homologue from closely related species. Predicting regulatory elements is even more difficult than exon prediction, but again, comparisons pinpointing conserved sequence motifs upstream of translation start sites are helping to unravel gene regulatory networks. In addition to interspecific studies, intraspecific sequence comparison yields insights into the evolutionary forces that have acted on a species in the past. Of particular interest here is the identification of selection events such as selective sweeps. Both intra- and interspecific sequence comparisons are based on a variety of computational methods, including alignment, phylogenetic reconstruction, and coalescent theory. This article surveys the biology and the central computational ideas applied in recent comparative genomics projects. We argue that the most fruitful method of understanding the functional content of genomes is to study them in the context of related genomic sequences. In particular, such a study may reveal selection, a fundamental pointer to biological relevance.

Characterization and structural features of a chalcone synthase mutation in a white-flowering line of Matthiola incana R. Br. (Brassicaceae)

For Matthiola incana (Brassicaceae), used as a model system to study biochemical and genetical aspects of anthocyanin biosynthesis, several nearly isogenic colored wild type lines and white-flowering mutant lines are available, each with a specific defect in the genes responsible for anthocyanin production (genes e, f, and g). For gene f supposed to code for chalcone synthase (CHS; EC 2.3.1.74), the key enzyme of the flavonoid/anthocyanin biosynthesis pathway belonging to the group of type III polyketide synthases (PKS), the wild type genomic sequence of M. incana line 04 was determined in comparison to the white-flowering CHS mutant line 18. The type of mutation in the chs gene was characterized as a single nucleotide substitution in a triplet AGG coding for an evolutionary conserved arginine into AGT coding for serine (R72S). Northern blots and RT-PCR demonstrated that the mutated gene is expressed in flower petals. Heterologous expression of the wild type and mutated CHS cDNA in E. Scherichia coli, verified by Western blotting and enzyme assays with various starter molecules, revealed that the mutant protein had no detectable activity, indicating that the strictly conserved arginine residue is essential for the enzymatic reaction. This mutation, which previously was not detected by mutagenic screening, is discussed in the light of structural and functional information on alfalfa CHS and related type III PKS enzymes.

Promoter analysis and transcription profiling: Integration of genetic data enhances understanding of gene expression

It is increasingly evident that transcription control might be conserved among organisms. For this reason, genome sequencing and gene expression profiling methods, which have yielded a plethora of data in different organisms, may be applied in species where genomic sequence is limited to mostly expression array and EST data. The identification of transcription factors and promoters associated with gene expression profiles and ESTs could therefore contribute to elucidate and predict complex regulatory events in plants.

Functional and phylogenetic analyses of a conserved regulatory program in the phloem of minor veins

The minor-vein phloem of mature leaves is developmentally and physiologically distinct from the phloem in the rest of the vascular system. Phloem loading of transport sugars occurs in the minor veins, and consistent with this, galactinol synthase is expressed in the minor veins of melon (Cucumis melo) as part of the symplastic-loading mechanism that operates in this species. A galactinol synthase promoter from melon drives gene expression in the minor-vein companion cells of both transgenic tobacco (Nicotiana tabacum) and Arabidopsis. Neither of these plants use galactinol in the phloem-loading process, implying that the promoter responds to a minor-vein-specific regulatory cascade that is highly conserved across a broad range of eudicotyledons. Detailed analysis of this promoter by truncation and mutagenesis identified three closely coupled sequences that unambiguously modulate tissue specificity. These sequences cooperate in a combinatorial fashion: two promote expression throughout the vascular system of the plant, whereas the third functions to repress expression in the larger bundles. In a complementary approach, phylogenetic footprinting was used to obtain single-nucleotide resolution of conserved sites in orthologous promoters from diverse members of the Cucurbitaceae. This comparative analysis confirmed the importance of the closely coupled sites but also revealed other highly conserved sequences that may modulate promoter strength or contribute to expression patterns outside of the phloem. The conservation of this regulatory design among species that phloem load by different mechanisms supports a model for organismal development in which tissues and cell types are controlled by relatively ancient and conserved paradigms but expression of genes influencing final form and function are relatively plastic.

And then there were many: MADS goes genomic

During the past decade, MADS-box genes have become known as key regulators in both reproductive and vegetative plant development. Traditional genetics and functional genomics tools are now available to elucidate the expression and function of this complex gene family on a much larger scale. Moreover, comparative analysis of the MADS-box genes in diverse flowering and non-flowering plants, boosted by bioinformatics, contributes to our understanding of how this important gene family has expanded during the evolution of land plants. Therefore, the recent advances in comparative and functional genomics should enable researchers to identify the full range of MADS-box gene functions, which should help us significantly in developing a better understanding of plant development and evolution.

Construction of gene expression system in hop (Humulus lupulus) lupulin gland using valerophenone synthase promoter

The promoter region of the valerophenone synthase (VPS) gene was isolated from hop (Humulus lupulus). VPS, a member of the chalcone synthase (CHS) super-family, catalyzes the biosynthesis reaction of the hop resin that significantly accumulates in the cone's secretory gland called the "lupulin gland". The typical H-box and G-box sequences, which exist in many plants' CHS promoters and act as cis-elements for tissue specificity, UV-light induction, etc., were not found in the isolated VPS promoter, although the H-box-like sequence (CCTTACC, CCTAACC) and the core sequence (ACGT) of the G-box were observed. The transformation experiment using the VPS promoter-UIDA gene fusion revealed that the promoter acts not only in the lupulin gland but also in the glands of leaf and stem. On the other hand, the VPS promoter activity was not induced by UV-irradiation.

Computational approaches to identify promoters and cis-regulatory elements in plant genomes

The identification of promoters and their regulatory elements is one of the major challenges in bioinformatics and integrates comparative, structural, and functional genomics. Many different approaches have been developed to detect conserved motifs in a set of genes that are either coregulated or orthologous. However, although recent approaches seem promising, in general, unambiguous identification of regulatory elements is not straightforward. The delineation of promoters is even harder, due to its complex nature, and in silico promoter prediction is still in its infancy. Here, we review the different approaches that have been developed for identifying promoters and their regulatory elements. We discuss the detection of cis-acting regulatory elements using word-counting or probabilistic methods (so-called "search by signal" methods) and the delineation of promoters by considering both sequence content and structural features ("search by content" methods). As an example of search by content, we explored in greater detail the association of promoters with CpG islands. However, due to differences in sequence content, the parameters used to detect CpG islands in humans and other vertebrates cannot be used for plants. Therefore, a preliminary attempt was made to define parameters that could possibly define CpG and CpNpG islands in Arabidopsis, by exploring the compositional landscape around the transcriptional start site. To this end, a data set of more than 5,000 gene sequences was built, including the promoter region, the 5'-untranslated region, and the first introns and coding exons. Preliminary analysis shows that promoter location based on the detection of potential CpG/CpNpG islands in the Arabidopsis genome is not straightforward. Nevertheless, because the landscape of CpG/CpNpG islands differs considerably between promoters and introns on the one side and exons (whether coding or not) on the other, more sophisticated approaches can probably be developed for the successful detection of "putative" CpG and CpNpG islands in plants.

Effects of inbreeding on the genetic diversity of populations

The study of variability within species is important to all biologists who use genetic markers. Since the discovery of molecular variability among normal individuals, data have been collected from a wide range of organisms, and it is important to understand the major factors affecting diversity levels and patterns. Comparisons of inbreeding and outcrossing populations can contribute to this understanding, and therefore studying plant populations is important, because related species often have different breeding systems. DNA sequence data are now starting to become available from suitable plant and animal populations, to measure and compare variability levels and test predictions.

Regulatory elements of the floral homeotic gene AGAMOUS identified by phylogenetic footprinting and shadowing

In Arabidopsis thaliana, cis-regulatory sequences of the floral homeotic gene AGAMOUS (AG) are located in the second intron. This 3-kb intron contains binding sites for two direct activators of AG, LEAFY (LFY) and WUSCHEL (WUS), along with other putative regulatory elements. We have used phylogenetic footprinting and the related technique of phylogenetic shadowing to identify putative cis-regulatory elements in this intron. Among 29 Brassicaceae species, several other motifs, but not the LFY and WUS binding sites identified previously, are largely invariant. Using reporter gene analyses, we tested six of these motifs and found that they are all functionally important for the activity of AG regulatory sequences in A. thaliana. Although there is little obvious sequence similarity outside the Brassicaceae, the intron from cucumber AG has at least partial activity in A. thaliana. Our studies underscore the value of the comparative approach as a tool that complements gene-by-gene promoter dissection but also demonstrate that sequence-based studies alone are insufficient for a complete identification of cis-regulatory sites.

Conserved noncoding sequences among cultivated cereal genomes identify candidate regulatory sequence elements and patterns of promoter evolution

Surveys for conserved noncoding sequences (CNS) among genes from monocot cereal species were conducted to assess the general properties of CNS in grass genomes and their correlation with known promoter regulatory elements. Initial comparisons of 11 orthologous maize-rice gene pairs found that previously defined regulatory motifs could be identified within short CNS but could not be distinguished reliably from random sequence matches. Among the different phylogenetic footprinting algorithms tested, the VISTA tool yielded the most informative alignments of noncoding sequence. VISTA was used to survey for CNS among all publicly available genomic sequences from maize, rice, wheat, barley, and sorghum, representing >300 gene comparisons. Comparisons of orthologous maize-rice and maize-sorghum gene pairs identified 20 bp as a minimal length criterion for a significant CNS among grass genes, with few such CNS found to be conserved across rice, maize, sorghum, and barley. The frequency and length of cereal CNS as well as nucleotide substitution rates within CNS were consistent with the known phylogenetic distances among the species compared. The implications of these findings for the evolution of cereal gene promoter sequences and the utility of using the nearly completed rice genome sequence to predict candidate regulatory elements in other cereal genes by phylogenetic footprinting are discussed.

Transcriptional regulation: a genomic overview

The availability of the Arabidopsis thaliana genome sequence allows a comprehensive analysis of transcriptional regulation in plants using novel genomic approaches and methodologies. Such a genomic view of transcription first necessitates the compilation of lists of elements. Transcription factors are the most numerous of the different types of proteins involved in transcription in eukaryotes, and the Arabidopsis genome codes for more than 1,500 of them, or approximately 6% of its total number of genes. A genome-wide comparison of transcription factors across the three eukaryotic kingdoms reveals the evolutionary generation of diversity in the components of the regulatory machinery of transcription. However, as illustrated by Arabidopsis, transcription in plants follows similar basic principles and logic to those in animals and fungi. A global view and understanding of transcription at a cellular and organismal level requires the characterization of the Arabidopsis transcriptome and promoterome, as well as of the interactome, the localizome, and the phenome of the proteins involved in transcription.