Conserved Non-Coding Regulatory Signatures in Arabidopsis Co-Expressed Gene Modules

Complex traits and other polygenic processes require coordinated gene expression. Co-expression networks model mRNA co-expression: the product of gene regulatory networks. To identify regulatory mechanisms underlying coordinated gene expression in a tissue-enriched context, ten Arabidopsis thaliana co-expression networks were constructed after manually sorting 4,566 RNA profiling datasets into aerial, flower, leaf, root, rosette, seedling, seed, shoot, whole plant, and global (all samples combined) groups. Collectively, the ten networks contained 30% of the measurable genes of Arabidopsis and were circumscribed into 5,491 modules. Modules were scrutinized for cis regulatory mechanisms putatively encoded in conserved non-coding sequences (CNSs) previously identified as remnants of a whole genome duplication event. We determined the non-random association of 1,361 unique CNSs to 1,904 co-expression network gene modules. Furthermore, the CNS elements were placed in the context of known gene regulatory networks (GRNs) by connecting 250 CNS motifs with known GRN cis elements. Our results provide support for a regulatory role of some CNS elements and suggest the functional consequences of CNS activation of co-expression in specific gene sets dispersed throughout the genome.     

Plant conserved non-coding sequences and paralogue evolution

Genome duplication is a powerful evolutionary force and is arguably most prominent in plants, where several ancient whole-genome duplication events have been documented. Models of gene evolution predict that functional divergence between duplicates (subfunctionalization) is caused by the loss of regulatory elements. Studies of conserved non-coding sequences (CNSs), which are putative regulatory elements, indicate that plants have far fewer CNSs per gene than mammals, suggesting that plants have less complex regulatory mechanisms. Furthermore, a recent study of a duplicated gene pair in maize suggests that CNSs are lost in a complementary fashion, perhaps driving subfunctionalization. If subfunctionalization is common, one expects duplicate genes to diverge in expression; recent microarray analyses in Arabidopsis thalinia suggest that this is the case. Plant genomes are relatively complex on a genomic level because of the prevalence of whole-genome duplication and, paradoxically, subfunctionalization after duplication can lead to relatively simple regulatory regions on a per gene basis.     

Evolution of Conserved Non-Coding Sequences Within the Vertebrate Hox Clusters Through the Two-Round Whole Genome Duplications Revealed by Phylogenetic Footprinting Analysis

As a result of two-round whole genome duplications, four or more paralogous Hox clusters exist in vertebrate genomes. The paralogous genes in the Hox clusters show similar expression patterns, implying shared regulatory mechanisms for expression of these genes. Previous studies partly revealed the expression mechanisms of Hox genes. However, cis-regulatory elements that control these paralogous gene expression are still poorly understood. Toward solving this problem, the authors searched conserved non-coding sequences (CNSs), which are candidates of cis-regulatory elements. When comparing orthologous Hox clusters of 19 vertebrate species, 208 intergenic conserved regions were found. The authors then searched for CNSs that were conserved not only between orthologous clusters but also among the four paralogous Hox clusters. The authors found three regions that are conserved among all the four clusters and eight regions that are conserved between intergenic regions of two paralogous Hox clusters. In total, 28 CNSs were identified in the paralogous Hox clusters, and nine of them were newly found in this study. One of these novel regions bears a RARE motif. These CNSs are candidates for gene expression regulatory regions among paralogous Hox clusters. The authors also compared vertebrate CNSs with amphioxus CNSs within the Hox cluster, and found that two CNSs in the HoxA and HoxB clusters retain homology with amphioxus CNSs through the two-round whole genome duplications.     

Expression pattern shifts following duplication indicative of subfunctionalization and neofunctionalization in regulatory genes of Arabidopsis

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

Investigating ancient duplication events in the Arabidopsis genome

The complete genomic analysis of Arabidopsis thaliana has shown that a major fraction of the genome consists of paralogous genes that probably originated through one or more ancient large-scale gene or genome duplication events. However, the number and timing of these duplications still remains unclear, and several different hypotheses have been put forward recently. Here, we reanalyzed duplicated blocks found in the Arabidopsis genome described previously and determined their date of divergence based on silent substitution estimations between the paralogous genes and, where possible, by phylogenetic reconstruction. We show that methods based on averaging protein distances of heterogeneous classes of duplicated genes lead to unreliable conclusions and that a large fraction of blocks duplicated much more recently than assumed previously. We found clear evidence for one large-scale gene or even complete genome duplication event somewhere between 70 to 90 million years ago. Traces pointing to a much older (probably more than 200 million years) large-scale gene duplication event could be detected. However, for now it is impossible to conclude whether these old duplicates are the result of one or more large-scale gene duplication events. abbreviations dA, fraction of amino acid substitutions; Kn, number of nonsynonymous substitutions per nonsynonymous site; Ks, number of synonymous substitutions per synonymous site; MYA, million years ago     

Sequence-level analysis of the diploidization process in the triplicated FLOWERING LOCUS C region of Brassica rapa

Strong evidence exists for polyploidy having occurred during the evolution of the tribe Brassiceae. We show evidence for the dynamic and ongoing diploidization process by comparative analysis of the sequences of four paralogous Brassica rapa BAC clones and the homologous 124-kb segment of Arabidopsis thaliana chromosome 5. We estimated the times since divergence of the paralogous and homologous lineages. The three paralogous subgenomes of B. rapa triplicated 13 to 17 million years ago (MYA), very soon after the Arabidopsis and Brassica divergence occurred at 17 to 18 MYA. In addition, a pair of BACs represents a more recent segmental duplication, which occurred approximately 0.8 MYA, and provides an exception to the general expectation of three paralogous segments within the B. rapa genome. The Brassica genome segments show extensive interspersed gene loss relative to the inferred structure of the ancestral genome, whereas the Arabidopsis genome segment appears little changed. Representatives of all 32 genes in the Arabidopsis genome segment are represented in Brassica, but the hexaploid complement of 96 has been reduced to 54 in the three subgenomes, with compression of the genomic region lengths they occupy to between 52 and 110 kb. The gene content of the recently duplicated B. rapa genome segments is identical, but intergenic sequences differ.     

Computational identification and analysis of immune-associated nucleotide gene family in Arabidopsis thaliana

GTP-binding proteins represent a ubiquitous regulatory mechanism in controlling growth and development in eukaryotes under normal and stress conditions. The IAN/GIMAP proteins belong to a novel family of functionally uncharacterized GTP-binding proteins expressed in both plant and vertebrate cells during anti-pathogenic responses. To gain novel insights into their roles in plants, we did genome-wide analysis of the IAN/GIMAP gene family. We identified 13 Arabidopsis IAN/GIMAP genes, which share similar gene structures and mostly reside in a tandem cluster on chromosomes. Sequence comparison reveals that these genes encode 26–52 kDa proteins with one GTP-binding domain and a conserved box unique to the family. Phylogenetic analysis suggests that the IAN/GIMAP genes of angiosperms and vertebrates may have evolved by independent gene duplication events. GENEVESTIGATOR sources were mined for comprehensive and comparative Arabidopsis IAN/GIMAP gene family expression analysis. These data reveal that IAN/GIMAPs exhibit diverse expression patterns during development and in response to external stimuli, indicating that these paralogous genes are likely involved in complex biological processes in Arabidopsis. Our present findings provide a basis for elucidating the novel GTPase family protein-mediated regulatory mechanisms in the future.     

Patterns of molecular evolution among paralogous floral homeotic genes.

The plant MADS-box regulatory gene family includes several loci that control different aspects of inflorescence and floral development. Orthologs to the Arabidopsis thaliana MADS-box floral meristem genes APETALA1 and CAULIFLOWER and the floral organ identity genes APETALA3 and PISTILLATA were isolated from the congeneric species Arabidopsis lyrata. Analysis of these loci between these two Arabidopsis species, as well as three other more distantly related taxa, reveal contrasting dynamics of molecular evolution between these paralogous floral regulatory genes. Among the four loci, the CAL locus evolves at a significantly faster rate, which may be associated with the evolution of genetic redundancy between CAL and AP1. Moreover, there are significant differences in the distribution of replacement and synonymous substitutions between the functional gene domains of different floral homeotic loci. These results indicate that divergence in developmental function among paralogous members of regulatory gene families is accompanied by changes in rate and pattern of sequence evolution among loci.     

Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes

It is often anticipated that many of today's diploid plant species are in fact paleopolyploids. Given that an ancient large-scale duplication will result in an excess of relatively old duplicated genes with similar ages, we analyzed the timing of duplication of pairs of paralogous genes in 14 model plant species. Using EST contigs (unigenes), we identified pairs of paralogous genes in each species and used the level of synonymous nucleotide substitution to estimate the relative ages of gene duplication. For nine of the investigated species (wheat [Triticum aestivum], maize [Zea mays], tetraploid cotton [Gossypium hirsutum], diploid cotton [G. arboretum], tomato [Lycopersicon esculentum], potato [Solanum tuberosum], soybean [Glycine max], barrel medic [Medicago truncatula], and Arabidopsis thaliana), the age distributions of duplicated genes contain peaks corresponding to short evolutionary periods during which large numbers of duplicated genes were accumulated. Large-scale duplications (polyploidy or aneuploidy) are strongly suspected to be the cause of these temporal peaks of gene duplication. However, the unusual age profile of tandem gene duplications in Arabidopsis indicates that other scenarios, such as variation in the rate at which duplicated genes are deleted, must also be considered.     

Structure and expression analysis of rice paleo duplications

Having a well-known history of genome duplication, rice is a good model for studying structural and functional evolution of paleo duplications. Improved sequence alignment criteria were used to characterize 10 major chromosome-to-chromosome duplication relationships associated with 1440 paralogous pairs, covering 47.8% of the rice genome, with 12.6% of genes that are conserved within sister blocks. Using a micro-array experiment, a genome-wide expression map has been produced, in which 2382 genes show significant differences of expression in root, leaf and grain. By integrating both structural (1440 paralogous pairs) and functional information (2382 differentially expressed genes), we identified 115 paralogous gene pairs for which at least one copy is differentially expressed in one of the three tissues. A vast majority of the 115 paralogous gene pairs have been neofunctionalized or subfunctionalized as 88%, 89% and 96% of duplicates, respectively, expressed in grain, leaf and root show distinct expression patterns. On the basis of a Gene Ontology analysis, we have identified and characterized the gene families that have been structurally and functionally preferentially retained in the duplication showing that the vast majority (>85%) of duplicated have been either lost or have been subfunctionalized or neofunctionalized during 50–70 million years of evolution.     

Functional diversification of the two C-class MADS box genes OSMADS3 and OSMADS58 in Oryza sativa

The C-class MADS box gene AGAMOUS (AG) plays crucial roles in Arabidopsis thaliana development by regulating the organ identity of stamens and carpels, the repression of A-class genes, and floral meristem determinacy. To examine the conservation and diversification of C-class gene function in monocots, we analyzed two C-class genes in rice (Oryza sativa), OSMADS3 and OSMADS58, which may have arisen by gene duplication before divergence of rice and maize (Zea mays). A knockout line of OSMADS3, in which the gene is disrupted by T-DNA insertion, shows homeotic transformation of stamens into lodicules and ectopic development of lodicules in the second whorl near the palea where lodicules do not form in the wild type but carpels develop almost normally. By contrast, RNA-silenced lines of OSMADS58 develop astonishing flowers that reiterate a set of floral organs, including lodicules, stamens, and carpel-like organs, suggesting that determinacy of the floral meristem is severely affected. These results suggest that the two C-class genes have been partially subfunctionalized during rice evolution (i.e., the functions regulated by AG have been partially partitioned into two paralogous genes, OSMADS3 and OSMADS58, which were produced by a recent gene duplication event in plant evolution).     

Organisation and structural evolution of the rice glutathione <Emphasis Type="Italic">S</Emphasis>-transferase gene family

Glutathione S-transferases (GSTs) comprise a large family of key defence enzymes against xenobiotic toxicity. Here we describe the comprehensive characterisation of this important multigene family in the model monocot species rice [Oryza sativa (L.)]. Furthermore, we investigate the molecular evolution of the family based on the analysis of (1) the patterns of within-genome duplication, and (2) the phylogenetic relationships and evolutionary divergence among rice, Arabidopsis, maize and soybean GSTs. By in-silico screening of the EST and genome divisions of the Genbank/EMBL/DDBJ database we have isolated 59 putative genes and two pseudogenes, making this the largest plant GST family characterised to date. Of these, 38 (62%) are represented by genomic and EST sequences and 23 (38%) are known only from their genomic sequences. A preliminary survey of EST collections shows a large degree of variability in gene expression between different tissues and environmental conditions, with a small number of genes (13) accounting for 80% of all ESTs. Rice GSTs are organised in four main phylogenetic classes, with 91% of all rice genes belonging to the two plant-specific classes Tau (40 genes) and Phi (16 genes). Pairwise identity scores range between 17 and 98% for proteins of the same class, and 7 and 21% for interclass comparisons. Rapid evolution by gene duplication is suggested by the discovery of two large clusters of 7 and 23 closely related genes on chromosomes 1 and 10, respectively. A comparison of the complete GST families in two monocot and two dicot species suggests a monophyletic origin for all Theta and Zeta GSTs, and no more than three common ancestors for all Phi and Tau genes.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by M.-A. Grandbastien     

Diversification of non-TIR class NB-LRR genes in relation to whole-genome duplication events in Arabidopsis

Arabidopsis thaliana is believed to have experienced at least two and possibly three whole-genome duplication events in its evolutionary history. In order to investigate the evolutionary relationships between these duplication events and diversification of disease resistance (R) genes, segmental-duplication events containing R genes belonging to the nucleotide binding-leucine rich repeat (NB-LRR) class were identified. Of 153 segmental-duplication events containing NB-LRR genes, only 22 contained NB-LRR genes in both members of the duplication pair, indicating a high frequency of NB-LRR gene loss after whole-genome duplication. The relative age of the duplication events was estimated based on the average synonymous substitution rate of the duplicated gene pairs in the segments. These data were combined with phylogenetic analyses. NB-LRR genes present in segment pairs derived from the most recent whole-genome duplication event, estimated to have occurred only 20 to 40 million years ago, occupy very distant branches of the NB-LRR phylogenetic tree. These data suggest that when NB-LRR clusters are duplicated as part of a whole-genome duplication, homoeologous NB-LRR genes are preferentially lost, either by eliminating one copy of the cluster or by eliminating individual genes such that only paralogous NB-LRR genes are maintained.     

Modeling the evolution of a classic genetic switch

Background  

The regulatory network underlying the yeast galactose-use pathway has emerged as a model system for the study of regulatory network evolution. Evidence has recently been provided for adaptive evolution in this network following a whole genome duplication event. An ancestral gene encoding a bi-functional galactokinase and co-inducer protein molecule has become subfunctionalized as paralogous genes (GAL1 and GAL3) in Saccharomyces cerevisiae, with most fitness gains being attributable to changes in cis-regulatory elements. However, the quantitative functional implications of the evolutionary changes in this regulatory network remain unexplored.