Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

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

Katie E Weeks1, Nadia A Chuzhanova2, Iain S Donnison3 and Ian M Scott4*

Author Affiliations

1 Cardiff School of Computer Science, Cardiff University, Queen's Buildings, 5 The Parade, Roath, Cardiff, CF24 3AA, UK

2 Department of Biological Sciences, University of Central Lancashire, Preston, PR1 2HE, UK

3 Institute of Grassland and Environmental Research, Aberystwyth, Ceredigion, SY23 3EB, UK

4 Institute of Biological Sciences, University of Wales, Aberystwyth, Ceredigion, SY23 3DA, UK

For all author emails, please log on.

BMC Evolutionary Biology 2007, 7:51  doi:10.1186/1471-2148-7-51

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/7/51


Received:2 October 2006
Accepted:2 April 2007
Published:2 April 2007

© 2007 Weeks et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

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 (CLZ), 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 CLZ 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 CLZ 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.

Background

Promoter sequences have been described as a vast and largely uncharted territory for evolutionary biologists [1]. One impediment to exploration is the difficulty of motif prediction in noncoding sequences (NCS): motif-discovery tools achieved detection rates of only 22–35% for transcription factor (TF) binding sites in recent benchmark studies [2,3]. Although it has long been recognized in principle [4] that evidence for motifs can be enhanced by comparing sequences of common ancestry, 'phylogenetic footprinting' of higher eukaryotes is still in a development and evaluation phase [5-8]. There are also perceived challenges in the use of sequence alignment for phylogenetic analysis of NCS [9], as complex mutational processes (slipped-strand mispairing, stem-loop secondary structure excision/repair, minute inversions, intramolecular recombination) are prevalent. In practice, however, Bremer et al. [10] found chloroplast NCS to be of similar utility to coding sequences in phylogenetic tree construction for asterids. This result confirmed that plant NCS contain evolutionary signal, which might be hypothesized to reside in the conserved motifs sought in phylogenetic footprinting. The present study sought to explore the extent to which phylogenetic footprints in plant 5'-NCS could be subjected to evolutionary analysis and interpretation. For this objective, we needed to conduct sufficiently comprehensive phylogenetic footprinting for meaningful evolutionary analysis of conserved sequence blocks.

We employed a greater taxonomic range than other phylogenetic footprinting studies of plant NCS, which have been confined to single families [6,8,11,12] or to a couple of species [13]. Much of the interest in promoter evolution lies in comparisons of paralogous genes (i.e. genes that diverged after a duplication event). In consequence, it must be noted, our dataset included several multigene families, and therefore was not optimized to investigate taxonomic phylogenies in the manner of Bremer et al. [10].

Recognizing limitations in individual motif discovery tools [2,3,7], we sought to maximize detection of conservation by combining distinct methodologies. Analysis of generalized Lempel-Ziv complexity (CLZ), played several roles in our study. CLZ measures the complexity of a text as the minimal number of steps in a defined procedure of its synthesis with the parsing rule: the next phrase is the longest seen previously. Many text compression algorithms are based on Lempel-Ziv parsing [14]. Computation of CLZ thus involves a decomposition of the text into repeated blocks, and an application to the discovery of structural regularities in genetic 'texts' was realized by Gusev et al. [15]. This method has identified arrays of conserved sequence blocks in NCS of vertebrates from fish to humans [16,17]. CLZ analysis has also been used to study human mutagenic mechanisms [18,19] and genomic architecture [20].

Our second tool was MOTIF SAMPLER, in which the probability of finding a particular motif is estimated using Gibbs sampling and modelling of the background sequence with a Markov model [21].

We complemented these tools with sequence alignment, including the DIALIGN and ALIGN-M algorithms designed for highly divergent sequences with only localized similarities, as seen in 5'-NCS. DIALIGN is based on a segment-to-segment comparison [22,23], while ALIGN-M uses a non-progressive local approach to guide alignments [24].

We focused on 5'-NCS of ribulose-1,5-bisphosphate carboxylase small subunit (rbcS) genes because of the exceptional corpus of knowledge against which analytical outcomes could be benchmarked. As the earliest nuclear protein-coding sequences in plants to be cloned, rbcS genes became paradigms for functional studies of plant promoters [25], and several classes of cis-elements were originally defined in rbcS promoters. Thus, the prototype of trihelix TFs was the nuclear protein GT-1, which binds to the 14-bp Box II and related motifs in light-responsive regions of the pea rbcS-3A promoter [26]. Box II versions featured in the earliest rbcS promoter alignments [25,27], and occur in other light-responsive genes [28], where they may be targets of calcium/calmodulin phototransduction [29].

Two further cis-elements discovered in rbcS promoters, the G-box and I-box, are common features in light-responsive promoters [28,30], and have been functionally characterized as dual components of a minimal light-responsive unit [31]. G-box binding factors (GBFs), identified using tomato rbcS-3A upstream sequences [32,33], are basic leucine zipper TFs interacting with the G-box core, CACGTG [34]. Dicot rbcS G-boxes interact with the HY5 GBF, which mediates phytochrome and cryptochrome signals in concert with COP and DET regulators [31,35].

The I-box, core motif GATAAGR, was also defined in rbcS promoters [27,33,36]. Its reverse, YCTTATC, was highlighted in rbcS and other light-regulated promoters by early motif searches [37,38]. Binding factors for the I-box are still being clarified. Functional interactions occurred in yeast between I-box sequences and recombinant zinc-finger GATA TFs from Arabidopsis [39]. I-box binding nuclear proteins reported in several species [40,41] may therefore include GATA TFs, though the first cloned I-box binding protein was a tomato Myb-like TF [42]. While the above rbcS cis-elements are the most studied, there is evidence for numerous further elements and DNA-protein interactions in rbcS promoters [30,32,43-48].

There is a particularly extensive history of characterization of rbcS promoters from pea, Petunia and tomato [25,49]. We analyzed these along with other studied rbcS 5'-NCS such as those of Arabidopsis [50] to provide a gradation of taxonomic relations and evolutionary distances. Conserved features shared by the plant families analyzed would have persisted since the Cretaceous, to which can be dated the divergence of eurosids I (represented by the Fabaceae) from eurosids II (Brassicaeae), and both from asterids (Solanaceae) [51].

Results

Phylogenetic footprinting

5'-NCS of up to 400 bp including proximal promoter and 5'-untranslated regions (5'-UTRs) were analyzed for 27 dicot rbcS genes. The rosid complement comprised all four Arabidopsis genes (three from a tandem locus), plus genes from Brassica and the legumes Phaseolus, Medicago and Pisum (pea). The Lycopersicon (tomato), Solanum (potato), Petunia and Nicotiana genes included representatives of all three solanaceous rbcS loci, which are distinguished by features including an extra (third) intron in 'locus 2' genes, and tandem duplicates at 'locus 3' [25]. Phylogenetic footprinting analyses were performed on the entire dataset, and separately on various subgroups, e.g. rosid, brassica, legume, or solanaceous genes, or genes of each solanaceous locus. Three methodologies were employed:

(1) CLZ analysis was used as proposed by Gusev et al. [15] to search for recurrent sequence blocks in the rbcS 5'-NCS. The CLZ measure is based on representation of a sequence by fragments that have been encountered before (in the same or other sequences). Let S = s1 ... sL be a nucleotide sequence of length L. Denote by S [i:j] the substring of S that starts at position i and ends at position j. A Lempel-Ziv decomposition of S is a partition of S into m consecutive fragments, S = S [1:i1] S [i1+1:i2]...S [im-1:L], such that the k-th component S [ik-1+1:ik] is the longest fragment downstream of position ik-1 for which an exact repeat has been encountered somewhere upstream of position ik-1+1. The number of fragments in the decomposition, CLZ (S) = m, is called the complexity of S with respect to direct repeats. For example, if S = TCGATCGAGAT, then the decomposition of S with respect to direct repeats is T-C-G-A-TCGA-GAT. Fragments 1, 2, 3 and 4 in this decomposition are of length one since respective nucleotides T, C, G and A occur for the first time. Exact copies of fragments 5 and 6 occur in positions 1 and 3 respectively. The CLZ of the sequence with respect to direct repeats equals 6. To find fragments repeated in different rbcS 5'-NCS, we concatenated multiple sequences for CLZ analysis.

(2) Overrepresented motifs were sought with MOTIF SAMPLER, using a range of program options for prior probabilities, lengths, numbers and overlaps of motifs. MOTIF SAMPLER can also vary the background Markov model order (i.e. dependency on a given number of preceding sequence positions). Thijs et al. [52] found higher order models improved robustness of motif recovery in Arabidopsis data. We found that the optimal Markov model order differed for different motifs: in 40 repeat runs, optimal model orders were zero for detection of blocks 06, 22 and 29, first for 10, 25 and 30, second for 23 and 28, and third for 08 and 20. (Blocks are defined in Table 1.)

Table 1. Conserved Blocks in rbcS 5'-NCS

(3) Sites of local congruence were sought in multiple sequence alignments produced by CLUSTALW, ALIGN-M and DIALIGN, with various gap penalty options for the first two. Collation of methodologies by mapping output from CLZ and Gibbs sampling analyses onto alignments yielded useful synergies. In particular, the alignments revealed arrays of blocks that occurred in several sequences in the same order, which increased confidence in less conserved block versions that occurred in the appropriate position relative to other blocks.

Our initial CLZ procedure specified blocks ≥ 8 bp with up to two mismatches, which identified 218 instances of 34 conserved blocks (on average 90% identical with their definitions in Table 1). Relaxation of the mismatch criterion for DIALIGN-aligned versions of these 34 CLZ-defined blocks exposed an additional 109 instances (of average 76% identity with definitions).

Conversely, mapping blocks from other tools clarified often complex alignments. When the full dataset was aligned by DIALIGN, 67% of aligned blocks split into an average 3.5 fragments, and 86% of blocks were co-aligned on average with 1.7 others. Nonetheless, with support from CLZ and MOTIF SAMPLER, 323 instances of 35 blocks were identified within alignments. MOTIF SAMPLER used independently found 291 instances of 35 blocks.

The complementarity of our different phylogenetic footprinting methods was demonstrated by the benchmarking exercise in Figure 1. In this exercise, each tool independently analyzed the full set of 27 5'-NCS, to test performance (versus the methodological consensus) in scoring each instance of the 12 most frequent blocks. Performance parameters, following Tompa et al. [3], were:

thumbnailFigure 1. Comparison of phylogenetic footprinting tools in predicting the 12 most frequent blocks in the 27 dicot rbcS 5'-NCS. See text for Sensitivity (equation 1) and PPV (equation 2) performance parameters. MOTIF SAMPLER was run 8 times each for background model orders 0–3, and with the prior probability of motif (p) at 0.3, the empirical value from the analytical consensus. Other option settings were s 0, M 1, n 3, w 11, x 0, r 5. Gap penalties for CLUSTALW (1) were: opening 15.0, extension 6.66; and for ALIGN-M and CLUSTALW (2): opening 8.0, extension 0.5. Mean performance parameters shown with standard deviation bars. Values sharing alphabet labels were not significantly different (Mann-Whitney U test, P > 0.05).

Sensitivity = nTP/(nTP+nFN)(1)

Positive Predictive Value (PPV) = nTP/(nTP+nFP)(2)

where nTP = number of 'true' positives (identified blocks found also by other tools), nFN = 'false' negatives (blocks not found though supported by other tools), and nFP = 'false' positives (blocks found but not supported by other tools). (Since every block instance had not been verified as a cis-element, the 'true' and 'false' concepts in these equations reflected sequence analysis performance rather than functionality prediction.)

CLZ analysis and MOTIF SAMPLER showed greater PPV in block prediction, but weaker sensitivities, than the best alignment methods (Figure 1). MOTIF SAMPLER's sensitivity for individual blocks correlated (r = 0.85, P < 0.001) with its log-likelihood statistic [21] that is optimized during Gibbs sampling. Among the alignment tools, DIALIGN and ALIGN-M, designed for highly divergent sequences with localized alignments, outperformed the CLUSTALW global alignment algorithm (Figure 1). The performance of CLUSTALW was significantly improved by reducing the gap penalties, though the PPV of DIALIGN and ALIGN-M remained superior (Figure 1).

MOTIF SAMPLER outputs statistical data, which helped estimate the significance of our phylogenetic footprinting results. Ten dummy datasets with different randomizations of every sequence were analyzed by MOTIF SAMPLER using background model orders 0–3. Randomization caused MOTIF SAMPLER to find on average 6.5-fold fewer pseudo-motif instances. Log-likelihood scores [21] for pseudo-motifs in the 10 dummy datasets were much lower (mean = 49.6, standard deviation = 17.3) than those of the original sequence motifs (mean = 188.7, standard deviation = 60.2), which differed from random at significance levels of P < 0.0001 (Kruskal-Wallis tests).

In summation, phylogenetic footprinting defined 36 conserved blocks, representing contiguous nucleotide sequences occurring in two or more rbcS 5'-NCS, and being of sufficient length, sequence fidelity and positional similarity to make their common evolutionary origin probable. A total of 338 instances of these blocks were identified in the dataset. A large majority (275 instances of 33 blocks) were supported by all three methodologies. Of these 33 blocks, another 37 instances were supported only by CLZ and alignments, and 5 more only by CLZ and MOTIF SAMPLER. Two other blocks (11 instances) were defined using only MOTIF SAMPLER and alignments, and a single block (10 instances) only by CLZ analysis.

Conserved blocks

All block instances are mapped for the rosid (brassica and legume) 5'-NCS in Figure 2, and solanaceous 5'-NCS in Figure 3. An average of 12.5 blocks were found in each gene. The blocks occurred in arrays whose relative order was absolutely conserved, so that the number-codes detailed in Table 1 consistently reflect relative block positions from 5' to 3' in all sequences. We therefore confirmed observations of Argüello-Astorga and Herrera-Estrella [28] on the existence in light-responsive plant promoters of 'conserved modular arrays' (CMAs), which they defined as 'short promoter regions, including at least two different DNA stretches larger than 6 bp (putative individual factor binding sites or phylogenetic footprints), in which nucleotidic sequence, spacing, and position relative to the transcription start site are conserved in a phylogenetic series'.

thumbnailFigure 2. Block structures of rosid rbcS 5'-NCS. Blocks individually coloured, and numbered (as Table 1) on first appearance from top. Horizontal dimension = block length (bp), vertical dimension proportional to identity with definitions in Table 1 (range: 40–100%). Complete [06-08-10-11] CMAs in dotted-line boxes. Blocks common to brassica and legume 5'-NCS joined by lines with block numbers in boxes. Blocks also found in solanaceous 5'-NCS indicated by unfilled arrowheads on Arabidopsis ats1A and Phaseolus rbcS-2. Red arrowheads show experimentally determined transcription start sites [46, 50, 86, 87].

thumbnailFigure 3. Block structures of solanaceous rbcS 5'-NCS. Sequences grouped as the 3 loci of Dean et al. [25]. Blocks individually coloured, and numbered (as Table 1) on first appearance from top. Horizontal dimension = block length (bp), vertical dimension proportional to identity with Table 1 definitions (range: 47–100%). Complete [06-08-10-11] CMAs in dotted-line boxes. Blocks also in rosid 5'-NCS indicated by unfilled arrowheads on Nicotiana rbcS-8B. Red arrowheads show experimentally determined transcription start sites [53, 54, 88, 89]. Stowaway-Le2 transposable element is mapped in tomato rbcS-1 [56].

Over a third of blocks were conserved in two or more plant families, but the remainder were distinctive to single families, or, in the case of the solanaceous genes, to particular orthologous loci identified by Dean et al. [25].

Blocks are listed in Table 1 with 'definitions' as typical instances, since for variable blocks a consensus would be dominated by ambiguous IUPAC code. The degree of conservation of each instance relative to the 'definition' is indicated by the vertical block dimensions in Figures 2 and 3; the 'definitions' were chosen to maximize these dimensions and do not necessarily represent importance in functional terms. Full sequences and locations of all block instances are in Additional File 1.

Additional file 1. Conserved blocks in rbcS 5'-NCS. Alignments and locations of conserved blocks in all sequences.

Format: PDF Size: 54KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

The 18 blocks asterisked in Table 1 have been recognized in past rbcS research. Of these, the motif most represented was the I-box (blocks 08, 09, 12, 18, 29). The reverse-strand I-box (block 29) immediately upstream of the TATA-box (block 30) was found by Grob and Stüber [37], who termed it the light-responsive element (LRE).

The I-box block 08 functions in a light-responsive dual unit with the G-box block 10 [31]. The I-G boxes unit represented by blocks [08–10] was found to be common in light-responsive promoters, and termed rbcS-CMA5 by Argüello-Astorga and Herrera-Estrella [28]. In rosid NCS, a second I-box downstream (block 12) occurred in an I-G-I boxes array postulated as ancestral by these authors. The TG-rich block 11, between the G-box (10) and second I-box (12), formed part of rbcS-CMA4 of Argüello-Astorga and Herrera-Estrella [28]. Block 11 usually corresponded to Motif 4 of Manzara and Gruissem [27] (but see later on Box II).

In the present study, the largest CMA found in all three plant families comprised blocks [06-08-10-11], in dotted-line boxes in Figures 2 and 3. Block 06 is a previously overlooked motif, but we found identical versions in similar relative locations in the caryophyllid genes Mesembryanthemum crystallinum rbcS-1 [EMBL L10212, -241 bp] and Spinacia oleracea rbcS-1 [EMBL X73236, -363 bp]. In pea rbcS-3A, block 06 overlapped the 5' flank of the box III* inverted GT-1 site [32], and so might be a site for a factor like 3AF5, a light-regulated phosphoprotein that binds the 5' flank of the similar downstream Box III [43]. The pea rbcS-3A 3AF5 and Box III sites themselves corresponded to legume-specific blocks 19 and 21, which with block 17 are equivalent to rbcS-CMA3 of Argüello-Astorga and Herrera-Estrella [28].

Block 17 coincided with the pea rbcS-3A Box II element, which is the prototype of GT-1 trihelix TF binding sites and a target of the calcium phototransduction pathway [26,29]. The variability of Box II-like motifs [28] was reflected in the low MOTIF SAMPLER consensus score [21] for block 17 (1.04), but this block was recognized by MOTIF SAMPLER with 85% sensitivity, and aligned in all dicot NCS by DIALIGN and ALIGN-M on its relatively conserved TGTGG sub-fragment. The Box II motifs of earlier alignments [25,27] corresponded to block 17 for most sequences, but to block 11 for tomato rbcS-2 and rbcS-3A. (Local alignments of sequences not available to the earlier authors confirm our assignments.)

The solanaceous 5'-NCS (Figure 2) yielded further previously identified motifs, whose functions generally remain uncertain. Blocks 22, 23 and 25 were components of rbcS-CMA2 [28] and identified by Manzara et al. [49] (Table 1). Likewise, the blocks 20 and 28 flanking rbcS-CMA2 were found by Manzara et al. [49] (Table 1).

On average, 10.2% of the length of those sequences with known transcription start sites was occupied by 5'-UTRs, though these were highly variable in extent (Figures 2, 3). Blocks 32–34 occurred in the proximities of transcription start sites. Only two blocks, 35 and 36, were located fully within 5'-UTRs, but each featured in multiple sequences in several species (Figures 2, 3).

Precisely half the blocks in Table 1 were newly identified in this study. These novel blocks were confined to single plant families, apart from the brassica blocks 02, 24 and 35 also found in a legume species. In most of the novel blocks, potential cis-elements could be speculatively identified using promoter databases (Table 1).

Protein-DNA interactions in tomato rbcS 5'-NCS have been extensively mapped by Gruissem and colleagues, using DNase I footprinting of promoter fragments in nuclear extracts from different organs [49,53,54]. As shown for locus 3 (Figure 4), over 90% of our conserved blocks overlapped with DNase-protected regions in the 5'-NCS where these authors had defined DNase footprints for both DNA strands. DNase-protected regions also included blocks 31, 34, 36, which have not been defined in past studies.

thumbnailFigure 4. Protein binding relative to sequence structure in rbcS 5'-NCS of tomato locus 3. Round-cornered rectangles with alphabet labels correspond to mapped regions of DNase protection, colour-coded by organ [53, 54]. Square-cornered rectangles are conserved blocks as Figure 3. Line plots show local CLZ bp-1 with respect to the [AT][GC] alphabet in 6-bp sliding-window profiles in 2-bp steps.

On the other hand, one-third of DNase-protected regions did not overlap with well defined blocks (Figure 4). These additional DNase-protected sequences tended to be very variable between genes and dominated by particular nucleotides (e.g. AT-rich regions). The latter feature can be formally translated as low complexity, as shown by the sliding-window profiles of CLZ [55] with respect to the [AT] [GC] alphabet in Figure 4. The association of DNase-protection with CLZ troughs implied functional roles for low-complexity regions.

One characterized mechanism for the introduction of AT-richness (and thus low CLZ) into an rbcS 5'-NCS sequence is the Stowaway-Le2 inverted repeat element in the tomato rbcS-1 sequence (Figure 3) [56,57]. Sliding-window CLZ profiles confirmed the Stowaway-Le2 sequence as one of the main low-complexity regions of the tomato rbcS-1 5'-NCS (not shown). DNase-protected regions do occur within the Stowaway-Le2 sequence [53,56,57].

Evolutionary analysis

The absolutely conserved relative order of blocks indicated common ancestry of all the dicot 5'-NCS studied (Figures 2, 3). This provided basic confirmation of the potential for evolutionary analysis of phylogenetic footprints, as these must share the evolutionary history of the plant taxa or gene loci with which they are associated. Minimum ages for blocks found in different species were estimated by reference to molecular clock dates for relevant taxon divergences (Table 2). For blocks common to paralogous loci, further evidence on minimum ages was available from recent studies of ancestral genome duplications (Table 2). Blanc et al. [58] produced a database of 45 duplicate chromosome segment pairs in the Arabidopsis genome, one of which (Figure 5) encompassed the ats1A and B genes. Comparisons of synonymous substitutions (Ks) in the duplicate genes indicated the relevant polyploidy event was roughly twice as ancient as the Brassica-Arabidopsis divergence [58]. Bowers et al. [59] similarly identified a duplication event, prior to the Brassica-Arabidopsis split, that generated 34 chromosome segment pairs, of which their segment α 25 encompassed the ats1A and B genes.

Table 2. Minimum Age Estimates for rbcS 5'-NCS Blocks

thumbnailFigure 5. Duplicated segment pairs on Arabidopsis chromosomes 1 and 5 encompassing ats1A and B family genes. Duplicate genes linked by inter-chromosome lines, with dotted lines for tandem arrays. Gene labels are for clarity, and may refer merely to putative functions. Based on block 0105451100840, PARALOGONS IN ARABIDOPSIS THALIANA database [73].

Within the duplicate segments containing the ats1A and B genes, the latter were among several examples of tandem arrays, others including LRK10L receptor-like Ser/Thr kinases [60]. Such tandem arrays, presumed due to unequal crossing over, account for up to 17% of all Arabidopsis genes [61], but their age range is currently uncertain [62].

In the Solanaceae, a large-scale genome duplication was dated to 18–23 million years ago (mya) from Ks distributions of duplicate tomato and potato genes [62]. Ks values for inter-locus comparisons of tomato rbcS coding sequences were consistent with formation of the 3 loci in this event (Figure 6). This must have occurred in a common ancestor, as Ks values for tomato and potato rbcS orthologues (Figure 6) were consistent with the much more recent speciation date estimated at 1.6–3.3 mya by Blanc and Wolfe [62].

thumbnailFigure 6. Levels of synonymous substitutions (Ks) in solanaceous rbcS coding sequences. Mean Ks (with standard deviation bars) are shown for comparisons of: all gene pairs from the 3 loci within tomato (white bars) or potato (black bars); or all paired tomato-potato (T-P) orthologues (grey bar). Brackets indicate Ks distribution peaks attributed by Blanc and Wolfe [62] to genome duplication or speciation events.

Estimated dates for major lineage divergences implied that 13 taxonomically widespread blocks were of Cretaceous antiquity at least (Table 2). These included the [06-08-10-11] CMA, block 12 (the second I-box in rosids), and blocks 17 (Box II), 23 (CAAT-box), 29 (LRE) and 30 (TATA-box). Another, block 20, remains poorly characterized in functional terms, but does bind protein (Figure 4), and was further noted in genes from the Amaranthaceae (Spinacia oleracea rbcS-1 [EMBL X73236, -191 bp]), and Malvaceae (Gossypium hirsutum rbcS [EMBL X54091, -186 bp]). Other Cretaceous blocks were three rosid blocks (02, 24 and 35) discovered in the present study. The remaining blocks were found only in single families but could be dated by clade divergence or gene duplication events to 18–54 mya (Table 2).

The occurrence of particular phylogenetic footprints at different levels in the taxonomic hierarchy (Table 2) indicated that the 5'-NCS might be amenable to phylogenetic analysis. Opinions differ, however, about phylogenetic analysis of NCS, particularly at higher taxonomic levels. NCS are seen as problematic for alignment and phylogenetic analysis because of their structural constraints, non-randomness of evolution, and mutational changes such as slipped-strand mispairing, stem-loop secondary structure excision/repair, minute inversions, and intramolecular recombination [9]. In practice, however, Bremer et al. [10] found chloroplast NCS to be of similar utility to coding sequences for asterid phylogenetics.

In view of the technical uncertainties and limited precedents for exploring evolutionary relations between 5'-NCS [9,10], we compared several distinct methodologies. First, given the role of CLZ analysis in our phylogenetic footprinting, comparison of the 5'-NCS based on this methodology was pertinent. A set of N sequences can be described in terms of their pairwise complexities, in the form of N vectors each containing N components. The (i,j) component is the pairwise CLZ with respect to direct repeats between sequences i and j. To some extent, pairwise CLZ measures an evolutionary distance between sequences by the number of steps required to produce sequence j from sequence i using it as a source of building blocks. Hierarchical cluster analysis of 5'-NCS in this format produced the dendrogram in Figure 7A. (As CLZ depends on sequence length, 5'-NCS shorter than the maximum length of 400 bp had the potential to yield anomalous results. Short sequences were therefore analyzed only if overall topology was robust to their inclusion; only Petunia SSU11A was omitted in consequence.)

thumbnailFigure 7. Dendrograms of rbcS 5'-NCS relations constructed by 3 methods. Groupings highlighted for rosid (brassica and legume) genes, and the 3 solanaceous loci. Non-unique gene symbols prefixed with binomial species initials. (A) Hierarchical cluster analysis, with each sequence defined as vector of CLZ values from pairwise decomposition by each of the others. Numerals indicate nodes with multiscale bootstrap resampling values ≥ 50% obtained by PVCLUST. (B) Parsimony analysis by PAUP* of DIALIGN alignments. 50% majority-rule consensus of 234 most-parsimonious trees shown with bootstrap values ≥ 50%. (C) Parsimony analysis by PAUP* of sequences defined by block characters. 50% majority-rule consensus of 882 most-parsimonious trees shown with bootstrap values ≥ 50%.

Secondly, more conventional analyses based on DNA parsimony or distance were applied to 5'-NCS aligned using DIALIGN, or CLUSTALW with the gap penalties found to be most effective in phylogenetic footprinting (Figure 1). (ALIGN-M was not usable as it does not produce complete alignments where sequence tracts are too divergent.) Figure 7B shows a consensus of most-parsimonious trees of DIALIGN-aligned 5'-NCS. (The short Petunia SSU11A sequence was also omitted from this tree for reasons discussed for Figure 7A.)

Our third method (Figure 7C) was a cladistic analysis of character-states defined as presence or absence of conserved blocks. All blocks in Figures 2, 3 were included: of these 96.9% had ≥ 50% identity with the definitions in Table 1. The remainder averaged 45.6% identity, and all but one had been found by three phylogenetic footprinting methods. Close inspections of aligned locations scored as absences were often suggestive of degenerate residues of blocks.

Several points of congruence between the dendrograms produced by these diverse analyses were identifiable, though bootstrap support for nodes was often moderate or weak (Figure 7). Themes included the clustering of the 5'-NCS by gene loci rather than by species. Thus, 5'-NCS of the Arabidopsis atsB tandem locus showed more affinity with the Brassica sequence than with Arabidopsis ats1A. This accorded with the conclusion of Bowers et al. [59] that the ancestral α duplication event occurred prior to the Brassica-Arabidopsis split, because 49–64% of relevant Brassica genes were more similar to one Arabidopsis gene than was the Arabidopsis duplicate.

Another theme was the segregation of the solanaceous 5'-NCS as the three loci deduced by Dean et al. [25] (Figure 7). Pairings of tomato and potato orthologues received particularly strong bootstrap support, consistent with a recent speciation [62]. In contrast, the coding sequences of tomato and potato locus3 instead segregated by species (Figure 8). Similar discrepancies between noncoding- and coding-sequence trees in several organisms have been attributed to gene conversion processes that have a greater effect on coding sequences [63]. Also consistent with gene conversion in the locus3 coding sequences were very low intralocus Ks values that would imply tandem duplication near the tomato-potato speciation time (Figure 6), which would be hard to reconcile with the more ancient relationships of their 5'-NCS to Petunia orthologues (Figure 7). Gene conversion in the Petunia locus 3 genes themselves was suggested by Dean et al. [64].

thumbnailFigure 8. Parsimony analysis of solanaceous rbcS coding sequences aligned by CLUSTALW. 50% majority-rule consensus of 182 most-parsimonious trees from branch-and-bound search shown with pea rbcS-3C as outgroup. Numerals indicate bootstrap values of nodes. Gene labels as Figure 7.

The locus 3 5'-NCS presented a consistent picture in that the tomato and potato A genes were resolved as basal members of a monophyletic group (Figure 7). In fact, tomato rbcS-3A was the only gene retaining the ancestral [06-08-10-11] CMA in the analyzed region (Figure 3). The most likely counterpart among the Petunia 5'-NCS analyzed was SSU11A (Figure 7C). Petunia SSU112, SSU491 and SSU911 grouped with the more derived 5'-NCS of the tomato and potato B and C genes.

The remaining solanaceous 5'-NCS grouped into loci 1 and 2 (Figure 7). In evolutionary trees based on CLUSTALW rather than DIALIGN alignments, the Petunia locus 1 gene SSU611 formed the outgroup to locus 2, while the tomato and potato locus 1 genes grouped with the A genes of locus 3 (not shown). The DIALIGN trees were preferred as they were supported by the alternative dendrograms, and because we rated the alignments from this algorithm most highly (Figure 1). Moreover, CLUSTALW alignments of transit peptides supported affinity of SSU611 with the other solanaceous locus 1 genes (not shown).

The basal solanaceous locus could not be confidently identified, as the basal position of locus 2 in two of the dendrograms in Figure 7 had only moderate bootstrap support. Clear guidance was not forthcoming from the coding sequence Ks values (Figure 6), or from parsimony analyses (Figure 8), in which outgroup choice influenced topology with respect to these two loci.

Discussion

Conserved blocks revealed by phylogenetic footprinting in dicot rbcS 5'-NCS formed an evolutionary hierarchy, from those common to plant families that diverged in the Cretaceous, to family-specific blocks with minimum estimated ages of only about 20 million years. Similar heterogeneity in longevity and clade-specificity of promoter motifs has been found in other organisms of ancient divergence. Among homologous human and rodent TF binding sites, for example, Dermitzakis and Clark [65] found 33 with shared functions, while 14 were human-specific and 17 rodent-specific.

The most ancient conserved blocks we found included those recognized earliest in rbcS research on the basis of functional importance (I-boxes, G-box, Box II, CAAT-box, TATA-box) [25,27,32], though several relatively unknown ones also fell in this category. Furthermore, we were able to extend CMAs postulated in previous studies [28]. Younger blocks were generally of less widely recognized function, and presumably had acquired roles in the more recent clades in which they had evolved. Simulations by Stone and Wray [66] of the acquisition by point mutation of novel TF binding sites, and their subsequent fixation within populations, indicated the evolution of new sites must be virtually inevitable over millions of years. In a theoretical population of 106 Arabidopsis plants with two generations per year, the fixation time for two 6-bp binding sites in a 200-bp region was only 270,000 years.

Evolutionary information in the 5'-NCS was sufficient for several formal computational methods to produce dendrograms in accordance with the existing classification of solanaceous rbcS loci based descriptively on sequence similarities, intron features and linkage relations [25,27,67]. The solanaceous locus 2 genes are distinguished as the only land plant rbcS genes with introns at three positions, while locus 3 is a distinctive tandem array of three 2-intron rbcS genes in tomato and potato, and probably six in Petunia. Gene duplications appear to have provided additional impetus in functional evolution of rbcS genes. For instance, the strongly expressed locus 3 tomato genes rbcS-3B and rbcS-3C [68] represented the most derived members of tandem arrays according to our dendrograms. It has been suggested that gene duplicates are conserved and subfunctionalized by regulatory mutations, because each duplicate must survive to complement lost expression for essential subfunctions in the other [69]. Duplicate gene preservation by such a process could be < 4 million years for a gene with ≥ 5 regulatory elements and a mutation rate of 10-7 per year [69]. Such a rapid preservation of duplicates may need to be invoked for locus 3, because of coincident estimates (18–23 mya) for the major ancestral genome duplication event [62] and for divergence of the Petunia clade [51]. In our dendrograms, segregation of Petunia SSU112, SSU491 and SSU911 with the more derived tomato and potato genes of locus 3 indicated that tandem duplications at locus 3 had occurred prior to the Petunia divergence, and had undergone relatively little subsequent sequence evolution.

Point mutations do not appear to have been the only evolutionary processes governing protein interactions in the rbcS 5'-NCS. Mechanisms such as slipped-strand mispairing [9] probably generated the relatively extensive and variable low-complexity tracts that coincided with known DNase footprints in the locus 3 tomato genes. Another example of the gross mutational processes that can occur in 5'-NCS was the Stowaway-Le2 transposable element in the tomato rbcS-1 sequence (Figure 3). The absence of this transposable element from the potato sequence [57] implies a recent insertion event in tomato.

A primary factor that facilitated our study was a suite of phylogenetic footprinting tools that complemented and cross-validated each other. The least known member of our toolkit was probably CLZ analysis, whose use deserves to increase with its availability as an internet tool [55]. Its intuitive process of sequence decomposition by repeated fragments proved useful not only for identification of conserved motifs, but also for highlighting low-complexity regions such as AT-rich tracts, and as a similarity measure for global sequence comparisons and hence dendogram construction. Otu and Sayood [70] formally examined CLZ as a new sequence distance measure for phylogenetic tree construction, and demonstrated that its lack of dependence on alignments or evolutionary models was particularly suited for sequences subject to segment-based modifications, including whole mitochondrial genomes of eutherian mammals. Promising alternative alignment-independent methods of sequence comparison have also been proposed using the general information theoretical concept of Kolmogorov complexity [71,72], of which CLZ is one explicitly computable implementation.

The dendrograms we produced using CLZ, and those obtained by parsimony analysis of DIALIGN alignments or block characters, were of sufficient consistency to confirm the presence of evolutionary information in plant 5'-NCS. The dataset was not designed to investigate taxonomic phylogenies, as it included several multigene families. Moreover, we would not claim that the dendrograms rival in quality those produced using coding sequences, as bootstrap support for nodes was often moderate or weak, and there were points of variance between the dendrograms. Further investigation is needed to establish the extent to which NCS might contribute to molecular phylogenetics. We do, however, conclude that current computational methods provide the potential for analysis of the evolution of gene expression in terms of promoter structure.

Conclusion

Comprehensive phylogenetic footprinting of dicot 5'-NCS revealed conserved modular arrays of recurrent sequence blocks. Transcriptional functionality was confirmed as an evolutionary basis for this conservation by coincidence of recurrent blocks with cis-elements and protein-binding sites. Evolutionary hierarchies were discernible within the assemblage of blocks, such that taxonomically widespread, and hence ancient, blocks could be distinguished from taxon-specific, more recent, ones.

Methods

Database information

Noncoding sequences (NCS) up to 400 bp including and immediately 5' to the ATG codon were obtained for the following genes [accession numbers, bp analyzed]: Arabidopsis thaliana ats1A [EMBL:X13611, 400], ats1B [EMBL:X14564, 400], ats2B [EMBL:X14564, 400], ats3B [EMBL:X14564, 400]; Brassica napus rbcS [EMBL:X61097, 400]; Phaseolus vulgaris rbcS-2 [EMBL:AF028707, 400]; Pisum sativum rbcS-E9 [EMBL:X00806, 400], rbcS-3A [EMBL:M21356, 400], rbcS-3C [EMBL:X04334, 331]; Medicago sativa rbcSK-1A [EMBL:X96847, 400]; Lycopersicon esculentum rbcS-1 [EMBL:X05982, 338], rbcS-2 [EMBL:X05983, 400], rbcS-3A [EMBL:X05984, 380], rbcS-3B [EMBL:X05985, 283], rbcS-3C [EMBL:X05986, 300]; Petunia × hybrida SSU112 [EMBL:X12990, 351], SSU11A [EMBL:X03821, 281], SSU301 [EMBL:X12986, 400], SSU491 [EMBL:X12988, 400], SSU611 [EMBL:X12987, 400], SSU911 [EMBL:X12989, 400]; Nicotiana plumbaginifolia rbcS-8B [EMBL:X13711, 400]; Solanum tuberosum rbcS-1 [EMBL:X69759, 400], rbcS-2A [EMBL:X69760, 400], rbcS-2B [EMBL:X69761, 400], rbcS-2C [EMBL:X69762, 400], rbcS-3 [EMBL:X69763, 382]; Zea mays rbcSZm1 [EMBL:S42508, 400]. Coding sequences were from the same accessions, except P. sativum rbcS-3A [EMBL:X04333] and Zea mays ZmrbcS [EMBL:Y00322].

Duplicated ancestral chromosome segments encompassing the Arabidopsis ats genes were identified (as block 0105451100840) in the PARALOGONS IN ARABIDOPSIS THALIANA database [58,73]. Potential cis-elements in the 5'-NCS were identified using the PLACE [74,75] database.

Sequence analysis

Recurrent sequence blocks were identified in rbcS 5'-NCS by Lempel-Ziv complexity (CLZ) decomposition. Lempel and Ziv [76] suggested measurement of sequence complexity by the number of steps required for the iterative generation (recovery) of a given sequence S from scratch, using two possible 'recovery' operations per iteration: either copy a fragment that has already been encountered in the recovered part of the sequence; or add (generate) a new symbol not encountered before. This iterative process, called a decomposition, represents a sequence S as a concatenation of m consecutive fragments, H(S) = S [1:i1]S [i1+1: i2]...S [im-1+1: im = N], where S [ik-1+1: ik] is a fragment copied or generated at k-th step, N is the length of the sequence and m = mH(S) is the number of steps in decomposition process. Among all possible decompositions the one with the minimum number of steps defines sequence complexity, i.e. <a onClick="popup('http://www.biomedcentral.com/1471-2148/7/51/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/7/51/mathml/M1">View MathML</a>. The minimum is ensured by copying at each step of the decomposition process the longest fragment that has been encountered before. Similarly, one can define a pair-wise complexity of sequences S and Q, C(S|Q), as the number of steps needed to recover Q from S (or S from Q). In this case, each fragment in the decomposition of Q is the longest one whose copy occurs anywhere in sequence S. Gusev et al. [15] proposed a linear algorithm for sequence decomposition and computation of CLZ with respect to various types of repeat (including direct and inverted repeats or any combination of them). Its implementation is available online at LZCOMPOSER [55,77].

The full algorithm used in this study follows. Step 1. For N 5'-NCS denoted as S1,..., SN (with corresponding lengths |S1|,..., |SN|), a new, concatenated sequence ŋ = S1#...# SN of length L = Σ |Si|, i = 1,..., N was defined. (The arbitrary symbol # separated the concatenated sequences.) Step 2. A Lempel-Ziv decomposition of ŋ into m consecutive fragments, ŋ [1:i1]ŋ [i1+1:i2] ... ŋ [im-1: L], was computed, such that ŋ [ik-1+1:ik] was the longest fragment downstream of position ik-1 for which a direct repeat occurred starting from position j(k) somewhere upstream of position ik-1+1, and ŋ [ik-1+1:ik] did not contain #. Pointers j(k) were expressed as pairs (sequence number, position within the sequence). Step 3. Fragments ≥ 8 bp that were common for at least two sequences were included in a vocabulary of 'blocks'. Only exact matches were considered in the decomposition process. However, when two or more consecutive fragments of a decomposition were identical to the respective substrings in another sequence, and when these fragments were separated by a similar number of nucleotides (± 1) then they were merged into a single block. All remaining sequences from the given dataset were scanned for the occurrence of these blocks. For each block, the origin in the decomposition and the entire track of occurrences in different sequences were traced, ensuring that the fragments found were independent of the sequence order in ŋ. Step 4. Fragments defined as the same block were aligned, including an extra 10 bp either side to check for possible block extension, and their consensus sequence was defined allowing for a given number of mismatches (initially two). Step 5. All sequences were then scanned for each block defined by its consensus. No steps involved a priori knowledge of cis-elements. A similar algorithm was used to search for inverted repeats [16], but we found too few of these for detailed analysis. The decomposition process is available online at LZCOMPOSER [55,77].

Matrices (N × N) of pairwise CLZ values for N sequences were produced on LZCOMPOSER (using the symmetrized matrix output with diagonals adjusted to 0). Sliding-window profiles of local CLZ along single sequences were also generated on LZCOMPOSER.

Overrepresented motifs in the 5'-NCS were also sought using MOTIF SAMPLER v3.1 [21]. The 5'-NCS were analyzed in 19 combinations, with program options in the following ranges: search (s), single stranded; prior probability of motif (p), 0.3–0.8; length of motif (w), 9–25 bp; number of different motifs (n), 3–20; number of instances of each motif per sequence (M), 1, 2 or undefined; allowed overlap (x), 1–9 bp; program repeat runs (r), 0–99. Background models of order 0–3 were used in the analysis.

Multiple alignments of 5'-NCS were performed with three algorithms: CLUSTALW v1.83 [78]; DIALIGN 2 [22] in the QALIGN v1.10T software of Sammeth et al. [79]; and ALIGN-M v2.3 [24]. Unless stated, gap penalties in both CLUSTALW and the S2P step of ALIGN-M were: opening 8.0; extension0.5. In the search process for conserved blocks, a total of 26 different sequence combinations were aligned with DIALIGN and/or ALIGN-M, and the blocks from CLZ analysis and MOTIF SAMPLER were mapped in the alignments.

Levels of synonymous substitutions (Ks) were obtained by multiple alignment of all the tomato and potato rbcS coding sequences by CLUSTALW (default gap penalties), followed by estimation of the matrix of pairwise Ks values by the method of Li [80] implemented in the R package SEQINR [81].

Dendrograms

Three strategies were used to produce dendrograms of the 27 dicot sequences, with the 5'-NCS of Zea rbcSZm1 included as outgroup. (1) Hierarchical cluster analysis was performed on matrices of pairwise CLZ values (see above), using Euclidean distance to measure similarity of the different rows. Dendrograms were produced by the unweighted pair group method with arithmetic mean (UPGMA), using PAST v1.34 [82] and the R package PVCLUST [81]. Statistical support was assessed using PVCLUST to calculate the approximately unbiased (AU) values of Shimodaira [83] by multiscale bootstrap resampling of 1000 pseudoreplications. (2)Evolutionary trees were produced from multiple sequence alignments created with DIALIGN or CLUSTALW. Trees were obtained, using PAUP* v4.0b10 [84] and PHYLIP v3.64 [85], by DNA parsimony or, by the neighbour-joining, UPGMA or Fitch-Margoliash methods, from DNA distance matrices produced with the Jukes-Cantor substitution model. (3)Cladistic analysis, using PAUP* v4.0b10 and PAST v1.34, was performed on the conserved blocks identified by sequence analyses. A character-state matrix of absence (0) or presence (1) of each block was created. Characters were assigned equal weight and Dollo status (i.e. a block could evolve only once, but could disappear at several points on the tree). The tree-bisection-reconnection heuristic was used to search for the most parsimonious topologies. For methods (2) and (3), the Zea sequence rbcSZm1 was specified as outgroup, and nodal support was estimated from 100 tenfold-replicated bootstrap pseudoreplicates.

Evolutionary trees of coding sequences were obtained by bootstrapped parsimony analysis in PAUP*v4.0b10 of sequences aligned by CLUSTALW (default gap penalties) or DIALIGN.

List of abbreviations

CLZ, Lempel-Ziv complexity; CMA, conserved modular array; GBF, G-box binding factor; LRE, light-responsive element; Ks, level of synonymous substitutions; mya, million years ago; PPV, Positive Predictive Value; NCS, noncoding sequences; TF, transcription factor; UPGMA, unweighted pair group method with arithmetic mean; UTR, untranslated region.

Authors' contributions

Bioinformatic and dendrogram analyses were carried out by KW, NC and IS. The study was designed and coordinated by NC, ID and IS. The manuscript was drafted by IS with contributions and approval by all authors.

Acknowledgements

This research was funded by the BBSRC Genes and Developmental Biology Committee. We are also grateful to Prof. N.J. Fiddian of Cardiff University School of Computer Science for his support.

References

  1. Wray GA, Hahn MW, Abouheif E, Balhoff JP, Pizer M, Rockman MV, Romano LA: The evolution of transcriptional regulation in eukaryotes.

    Mol Biol Evol 2003, 20:1377-1419. PubMed Abstract | Publisher Full Text OpenURL

  2. Hu J, Li B, Kihara D: Limitations and potentials of current motif discovery algorithms.

    Nucleic Acids Res 2005, 33:4889-4913. OpenURL

  3. Tompa M, Li N, Bailey TL, Church GM, De Moor B, Eskin E, Favorov AA, Frith MC, Fu Y, Kent WJ, Makeev VJ, Mironov AA, Noble WS, Pavesi G, Pesole G, Régnier M, Simonis N, Sinha S, Thijs G, Van Helden J, Vandenbogaert M, Weng Z, Workman C, Ye C, Zhu Z: Assessing computational tools for the discovery of transcription factor binding sites.

    Nature Biotech 2005, 23:137-144. Publisher Full Text OpenURL

  4. Karlin S, Ghandour G, Foulser DE: DNA sequence comparisons of the human, mouse, and rabbit immunoglobulin kappa gene.

    Mol Biol Evol 1985, 2:35-52. PubMed Abstract | Publisher Full Text OpenURL

  5. Clark AG: The search for meaning in noncoding DNA.

    Genome Res 2001, 11:1319-1320. PubMed Abstract | Publisher Full Text OpenURL

  6. Guo H, Moose SP: Conserved noncoding sequences among cultivated cereal genomes identify candidate regulatory sequence elements and patterns of promoter evolution.

    Plant Cell 2003, 15:1143-1158. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Van Hellemont R, Monsieurs P, Thijs G, De Moor B, de Peer YV, Marchal K: A novel approach to identifying regulatory motifs in distantly related genomes.

    Genome Biol 2005, 6(13):R113. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  8. Colinas J, Birnbaum K, Benfey PN: Using cauliflower to find conserved non-coding regions in Arabidopsis.

    Plant Physiol 2002, 129:451-454. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Kelchner SA: The evolution of non-coding chloroplast DNA and its application in plant systematics.

    Ann Mo Bot Gard 2000, 87:482-498. Publisher Full Text OpenURL

  10. Bremer B, Bremer K, Heidari N, Erixon P, Olmstead RG, Anderberg AA, Källersjö M, Barkhordian E: Phylogenetics of asterids based on 3 coding and 3 non-coding chloroplast DNA markers and the utility of non-coding DNA at higher taxonomic levels.

    Mol Phylogenet Evol 2002, 24:274-301. PubMed Abstract | Publisher Full Text OpenURL

  11. Koch MA, Weisshaar B, Kroymann J, Haubold B, Mitchell-Olds T: Comparative genomics and regulatory evolution: conservation and function of the Chs and Apetala3 promoters.

    Mol Biol Evol 2001, 18:1882-1891. PubMed Abstract | Publisher Full Text OpenURL

  12. Kaplinsky NJ, Braun DM, Penterman J, Goff SA, Freeling M: Utility and distribution of conserved noncoding sequences in the grasses.

    Proc Natl Acad Sci USA 2002, 99:6147-6151. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. De Bodt S, Theissen G, Van de Peer Y: Promoter analysis of MADS-box genes in eudicots through phylogenetic footprinting.

    Mol Biol Evol 2006, 23:1293-1303. PubMed Abstract | Publisher Full Text OpenURL

  14. Verdú S: Fifty years of Shannon theory.

    IEEE T Inform Theory 1998, 44:2057-2078. Publisher Full Text OpenURL

  15. Gusev VD, Nemytikova LA, Chuzhanova NA: On the complexity measures of genetic sequences.

    Bioinformatics 1999, 15:994-999. PubMed Abstract | Publisher Full Text OpenURL

  16. Chuzhanova NA, Krawczak M, Nemytikova LA, Gusev VD, Cooper DN: Promoter shuffling has occurred during the evolution of the vertebrate growth hormone gene.

    Gene 2000, 254:9-18. PubMed Abstract | Publisher Full Text OpenURL

  17. Chuzhanova NA, Krawczak M, Thomas N, Nemytikova LA, Gusev VD, Cooper DN: The evolution of the vertebrate beta-globin gene promoter.

    Evolution 2002, 56:224-232. PubMed Abstract | Publisher Full Text OpenURL

  18. Bacolla A, Jaworski A, Larson JE, Jakupciak JP, Chuzhanova N, Abeysinghe SS, O'Connell CD, Cooper DN, Wells RD: Breakpoints of gross deletions coincide with non-B DNA conformations.

    Proc Natl Acad Sci USA 2004, 101:14162-14167. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Chuzhanova NA, Anassis EJ, Ball EV, Krawczak M, Cooper DN: Meta-analysis of indels causing human genetic disease: Mechanisms of mutagenesis and the role of local DNA sequence complexity.

    Hum Mutat 2003, 21:28-44. PubMed Abstract | Publisher Full Text OpenURL

  20. Kehrer-Sawatzki H, Sandig C, Chuzhanova N, Goidts V, Szamalek JM, Tanzer S, Muller S, Platzer M, Cooper DN, Hameister H: Breakpoint analysis of the pericentric inversion distinguishing human chromosome 4 from the homologous chromosome in the chimpanzee (Pan troglodytes).

    Hum Mutat 2005, 25:45-55. PubMed Abstract | Publisher Full Text OpenURL

  21. Thijs G, Marchal K, Lescot M, Rombauts S, de Moor B, Rouzé P, Moreau Y: A Gibbs sampling method to detect overrepresented motifs in the upstream regions of coexpressed genes.

    J Comput Biol 2002, 9:447-464. PubMed Abstract | Publisher Full Text OpenURL

  22. Morgenstern B: DIALIGN 2: improvement of the segment-to-segment approach to multiple sequence alignment.

    Bioinformatics 1999, 15:211-218. PubMed Abstract | Publisher Full Text OpenURL

  23. Morgenstern B, Frech K, Dress A, Werner T: DIALIGN: Finding local similarities by multiple sequence alignment.

    Bioinformatics 1998, 14:290-294. PubMed Abstract | Publisher Full Text OpenURL

  24. Van Walle I, Lasters I, Wyns L: Align-m - a new algorithm for multiple alignment of highly divergent sequences.

    Bioinformatics 2004, 20:1428-1435. PubMed Abstract | Publisher Full Text OpenURL

  25. Dean C, Pichersky E, Dunsmuir P: Structure, evolution and regulation of RbcS genes in higher plants.

    Annu Rev Plant Physiol Plant Mol Biol 1989, 40:415-439. Publisher Full Text OpenURL

  26. Gilmartin PM, Memelink J, Hiratsuka K, Kay SA, Chua NH: Characterization of a gene encoding a DNA binding protein with specificity for a light-responsive element.

    Plant Cell 1992, 4:839-849. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Manzara T, Gruissem W: Organization and expression of the genes encoding ribulose-1,5-bisphosphate carboxylase in higher plants.

    Photosynth Res 1988, 16:117-139. Publisher Full Text OpenURL

  28. Argüello-Astorga G, Herrera-Estrella L: Ancestral multipartite units in light-responsive plant promoters have structural features correlating with specific phototransduction pathways.

    Plant Physiol 1996, 112:1151-1166. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Wu Y, Hiratsuka K, Neuhaus G, Chua NH: Calcium and cGMP target distinct phytochrome-responsive elements.

    Plant J 1996, 10:1149-1154. PubMed Abstract | Publisher Full Text OpenURL

  30. Argüello-Astorga G, Herrera-Estrella L: Evolution of light-regulated plant promoters.

    Annu Rev Plant Physiol Plant Mol Biol 1998, 49:525-555. PubMed Abstract | Publisher Full Text OpenURL

  31. Martínez-Hernández A, López-Ochoa L, Argüello-Astorga G, Herrera-Estrella L: Functional properties and regulatory complexity of a minimal RBCS light-responsive unit activated by phytochrome, cryptochrome, and plastid signals.

    Plant Physiol 2002, 128:1223-1233. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Gilmartin PM, Sarokin L, Memelink J, Chua NH: Molecular light switches for plant genes.

    Plant Cell 1990, 2:369-378. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Guiliano G, Pichersky E, Malik VS, Timko MP, Scolnik PA, Cashmore A: An evolutionarily conserved protein binding sequence upstream of a plant light-regulated gene.

    Proc Natl Acad Sci USA 1988, 85:7089-7093. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Sibéril Y, Doireau P, Gantet P: Plant bZIP G-box binding factors.

    Eur J Biochem 2001, 268:5655-5666. PubMed Abstract | Publisher Full Text OpenURL

  35. Chattopadhyay S, Ang LH, Puente P, Deng XW, Nei W: Arabidopsis bZIP protein HY5 directly interacts with light-responsive promoters in mediating light control of gene expression.

    Plant Cell 1998, 10:673-683. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Donald RGK, Cashmore A: Mutation of either G box or I box sequences profoundly affects expression from the Arabidopsis rbcS1A promoter.

    EMBO J 1990, 9:1717-1726. PubMed Abstract | PubMed Central Full Text OpenURL

  37. Grob U, Stüber K: Discrimination of phytochrome dependent light inducible from non-light inducible plant genes. Prediction of a common light-responsive element (LRE) in phytochrome dependent light inducible plant genes.

    Nucleic Acids Res 1987, 15:9957-9973. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Dean C, Vandenelzen P, Tamaki S, Dunsmuir P, Bedbrook J: Differential expression of the 8 genes of the petunia ribulose bisphosphate carboxylase small subunit multi-gene family.

    EMBO J 1985, 4:3055-3061. PubMed Abstract | PubMed Central Full Text OpenURL

  39. Teakle GR, Manfield IW, Graham JF, Gilmartin PM: Arabidopsis thaliana GATA factors: organisation, expression and DNA-binding characteristics.

    Plant Mol Biol 2002, 50:43-57. PubMed Abstract | Publisher Full Text OpenURL

  40. Borello U, Ceccarelli E, Guiliano G: Constitutive, light-responsive and circadian clock-responsive factors compete for the different I box elements in plant light-regulated promoters.

    Plant J 1993, 4:611-619. PubMed Abstract | Publisher Full Text OpenURL

  41. Buzby JS, Yamada T, Tobin EM: A light-regulated DNA-binding activity interacts with a conserved region of a Lemna gibba rbcS promoter.

    Plant Cell 1990, 2:805-814. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Rose A, Meier I, Wienand U: The tomato I-box binding factor LeMYB1 is a member of a novel class of Myb-like proteins.

    Plant J 1999, 20:641-652. PubMed Abstract | Publisher Full Text OpenURL

  43. Sarokin L, Chua NH: Binding sites for two novel phosphoproteins, 3AF5 and 3AF3, are required for rbcS-3A expression.

    Plant Cell 1992, 4:473-483. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Lam E, Kano-Murakami Y, Gilmartin P, Niner B, Chua NH: A metal-dependent DNA-binding protein interacts with a constitutive element of a light-responsive promoter.

    Plant Cell 1990, 2:857-866. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Meier I, Callan KL, Fleming AJ, Gruissem W: Organ-specific differential regulation of a promoter subfamily for the ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit genes in tomato.

    Plant Physiol 1995, 107:1105-1118. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Urwin NAR, Jenkins GI: A sucrose repression element in the Phaseolus vulgaris rbcS2 gene promoter resembles elements responsible for sugar stimulation of plant and mammalian genes.

    Plant Mol Biol 1997, 35:929-942. PubMed Abstract | Publisher Full Text OpenURL

  47. Ueda T, Pichersky E, Malik VS, Cashmore A: Level of expression of the tomato rbcS-3A gene is modulated by a far upstream promoter element in a developmentally regulated manner.

    Plant Cell 1989, 1:217-227. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Datta N, Cashmore A: Binding of a pea nuclear protein to promoters of certain photoregulated genes is modulated by phosphorylation.

    Plant Cell 1989, 1:1069-1077. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Manzara T, Carrasco P, Gruissem W: Developmental and organ-specific changes in promoter DNA-protein interactions in the tomato rbcS gene family.

    Plant Cell 1991, 3:1305-1316. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Krebbers E, Seurinck J, Herdies L, Cashmore AR, Timko MP: Four genes in two diverged subfamilies encode the ribulose-1,5-bisphosphate carboxylase small subunit polypeptides of Arabidopsis thaliana.

    Plant Mol Biol 1988, 11:745-759. Publisher Full Text OpenURL

  51. Wikström N, Savolainen V, Chase MW: Evolution of the angiosperms: calibrating the family tree.

    Proc R Soc Lond B 2001, 268:2211-2220. Publisher Full Text OpenURL

  52. Thijs G, Lescot M, Marchal K, Rombauts S, De Moor B, Rouze P, Moreau Y: A higher-order background model improves the detection of promoter regulatory elements by Gibbs sampling.

    Bioinformatics 2001, 17:1113-1122. PubMed Abstract | Publisher Full Text OpenURL

  53. Manzara T, Carrasco P, Gruissem W: Developmental and organ-specific changes in DNA-protein interactions in the tomato rbcS1, rbcS2 and rbcS3A promoter regions.

    Plant Mol Biol 1993, 21:69-88. PubMed Abstract | Publisher Full Text OpenURL

  54. Carrasco P, Manzara T, Gruissem W: Developmental and organ-specific changes in DNA-protein interactions in the tomato rbcS3B and rbcS3C promoter regions.

    Plant Mol Biol 1993, 21:1-15. PubMed Abstract | Publisher Full Text OpenURL

  55. Orlov YL, Potapov VN: Complexity: an internet resource for analysis of DNA sequence complexity.

    Nucleic Acids Res 2004, 32:W628-W633. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Blume B, Barry CS, Hamilton AJ, Bouzayen M, Grierson D: Identification of transposon-like elements in non-coding regions of tomato ACC oxidase genes.

    Mol Gen Genet 1997, 254:297-303. PubMed Abstract | Publisher Full Text OpenURL

  57. Bureau TE, Wessler SR: Stowaway: a new family of inverted repeat elements associated with the genes of both monocotyledonous and dicotyledonous plants.

    Plant Cell 1994, 6:907-916. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Blanc G, Hokamp K, Wolfe K: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome.

    Genome Res 2003, 13:137-144. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Bowers JE, Chapman BA, Rong J, Paterson AH: Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events.

    Nature 2003, 422:433-438. PubMed Abstract | Publisher Full Text OpenURL

  60. Shiu SH, Bleecker AB: Expansion of the receptor-like kinase/Pelle gene family and receptor-like proteins in Arabidopsis.

    Plant Physiol 2003, 132:530-543. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Zhang L, Gaut BS: Does recombination shape the distribution and evolution of tandemly arrayed genes (TAGs) in the Arabidopsis thaliana genome?

    Genome Res 2003, 13:2533-2540. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Blanc G, Wolfe KH: Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes.

    Plant Cell 2004, 16:1667-1678. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Angers B, Gharbi K, Estoup A: Evidence of gene conversion events between paralogous sequences produced by tetraploidization in Salmoninae fish.

    J Mol Evol 2002, 54:501-510. PubMed Abstract | Publisher Full Text OpenURL

  64. Dean C, van den Elzen P, Tamaki S, Black M, Dunsmuir P, Bedbrook J: Molecular characterization of the rbcS multi-gene family of Petunia (Mitchell).

    Mol Gen Genet 1987, 206:465-474. Publisher Full Text OpenURL

  65. Dermitzakis ET, Clark AG: Evolution of transcription factor binding sites in mammalian gene regulatory regions: conservation and turnover.

    Mol Biol Evol 2002, 19:1114-1121. PubMed Abstract | Publisher Full Text OpenURL

  66. Stone JR, Wray GA: Rapid evolution of cis-regulatory sequences via local point mutations.

    Mol Biol Evol 2001, 18:1764-1770. PubMed Abstract | Publisher Full Text OpenURL

  67. Fritz CC, Wolter FP, Schenkemeyer V, Herget T, Schreier PH: The gene family encoding the ribulose-(1,5)-bisphosphate carboxylase/oxygenase (Rubisco) small subunit of potato.

    Gene 1993, 137:271-274. PubMed Abstract | Publisher Full Text OpenURL

  68. Sugita M, Gruissem W: Developmental, organ-specific, and light-dependent expression of the tomato ribulose-1,5-bisphosphate carboxylase small subunit gene family.

    Proc Natl Acad Sci USA 1987, 84:7104-7108. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  69. Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations.

    Genetics 1999, 151:1531-1545. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  70. Otu HH, Sayood K: A new sequence distance measure for phylogenetic tree construction.

    Bioinformatics 2003, 19:2122-2130. PubMed Abstract | Publisher Full Text OpenURL

  71. Varré JS, Delahaye JP, Rivals E: Transformation distances: a family of dissimilarity measures based on movements of segments.

    Bioinformatics 1999, 15:194-202. PubMed Abstract | Publisher Full Text OpenURL

  72. Li M, Badger JH, Chen X, Kwong S, Kearney P, Zhang H: An information-based sequence distance and its application to whole mitochondrial genome phylogeny.

    Bioinformatics 2001, 17:149-154. PubMed Abstract | Publisher Full Text OpenURL

  73. Paralogons in Arabidopsis thaliana [http://wolfe.gen.tcd.ie/athal/dup] webcite

  74. Higo K, Ugawa Y, Iwamoto M, Korenaga T: Plant cis-acting regulatory DNA elements (PLACE) database: 1999.

    Nucleic Acids Res 1999, 27:297-300. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  75. Plant cis-acting regulatory DNA elements (PLACE) database [http://www.dna.affrc.go.jp/PLACE/] webcite

  76. Lempel A, Ziv J: Complexity of finite sequences.

    IEEE T Inform Theory 1976, 22(1):75-81. Publisher Full Text OpenURL

  77. LZcomposer [http://wwwmgs.bionet.nsc.ru/mgs/programs/lzcomposer/complete.html] webcite

  78. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position specific gap penalties and weight matrix choice.

    Nucleic Acids Res 1994, 22:4673-4680. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  79. Sammeth M, Rothgänger J, Esser W, Albert J, Stoye J, Harmsen D: QAlign: quality-based multiple alignments with dynamic phylogenetic analysis.

    Bioinformatics 2003, 19:1592–1593. OpenURL

  80. Li WH: Unbiased estimation of the rates of synonymous and nonsynonymous substitution.

    J Mol Evol 1993, 36:96-99. PubMed Abstract | Publisher Full Text OpenURL

  81. The R Project for Statistical Computing [http://www.r-project.org/] webcite

  82. Hammer O, Harper DAT, Ryan PD: PAST: paleontological statistics software package for education and data analysis.

    Palaeontologia Electronica 2001, 4(1):http://palaeo-electronica.org. OpenURL

  83. Shimodaira H: An approximately unbiased test of phylogenetic tree selection.

    Syst Biol 2002, 51:492-508. PubMed Abstract | Publisher Full Text OpenURL

  84. Swofford DL: PAUP*: Phylogenetic Analysis using Parsimony (* and Other Methods). Version 4. Sunderland, Massachusetts , Sinauer Associates; 1998. PubMed Abstract OpenURL

  85. Felsenstein J: PHYLIP (Phylogeny Inference Package). Version 3.64 edition. Seattle , University of Washington; 2004.

  86. Fiebig C, Link G: 5'-upstream cis-elements and binding factor(s) potentially involved in light-regulated expression of a Brassica napus rbcS gene.

    Curr Genet 1992, 21:161-168. PubMed Abstract | Publisher Full Text OpenURL

  87. Fluhr R, Moses P, Morelli G, Coruzzi G, Chua NH: Expression dynamics of the pea rbcS multigene family and organ distribution of the transcripts.

    EMBO J 1986, 5:2063-2071. PubMed Abstract | PubMed Central Full Text OpenURL

  88. Dean C, Favreau M, Dunsmuir P, Bedbrook J: Confirmation of the relative expression levels of the Petunia (Mitchell) rbcS genes.

    Nucleic Acids Res 1987, 15:4655-4668. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  89. Poulsen C, Chua NH: Dissection of 5' upstream sequences for selective expression of the Nicotiana plumbaginifolia rbcS-8B gene.

    Mol Gen Genet 1988, 214:16-23. PubMed Abstract | Publisher Full Text OpenURL

  90. Kuhlemeier C, Cuozzo M, Green PJ, Goyvaerts E, Ward K, Chua NH: Localization and conditional redundancy of regulatory elements in rbcS-3A, a pea gene encoding the small subunit of ribulose-bisphosphate carboxylase.

    Proc Natl Acad Sci USA 1988, 85:4662-4666. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  91. Lavin M, Herendeen PS, Wojciechowski MF: Evolutionary rates analysis of Leguminosae implicates a rapid diversification of lineages during the Tertiary.

    Syst Biol 2005, 54:530–549. Publisher Full Text OpenURL

  92. Koch MA, Haubold B, Mitchell-Olds T: Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae).

    Mol Biol Evol 2000, 17:1483-1498. PubMed Abstract | Publisher Full Text OpenURL