Skip to main content

Evolution of plant senescence

Abstract

Background

Senescence is integral to the flowering plant life-cycle. Senescence-like processes occur also in non-angiosperm land plants, algae and photosynthetic prokaryotes. Increasing numbers of genes have been assigned functions in the regulation and execution of angiosperm senescence. At the same time there has been a large expansion in the number and taxonomic spread of plant sequences in the genome databases. The present paper uses these resources to make a study of the evolutionary origins of angiosperm senescence based on a survey of the distribution, across plant and microbial taxa, and expression of senescence-related genes.

Results

Phylogeny analyses were carried out on protein sequences corresponding to genes with demonstrated functions in angiosperm senescence. They include proteins involved in chlorophyll catabolism and its control, homeoprotein transcription factors, metabolite transporters, enzymes and regulators of carotenoid metabolism and of anthocyanin biosynthesis. Evolutionary timelines for the origins and functions of particular genes were inferred from the taxonomic distribution of sequences homologous to those of angiosperm senescence-related proteins. Turnover of the light energy transduction apparatus is the most ancient element in the senescence syndrome. By contrast, the association of phenylpropanoid metabolism with senescence, and integration of senescence with development and adaptation mediated by transcription factors, are relatively recent innovations of land plants. An extended range of senescence-related genes of Arabidopsis was profiled for coexpression patterns and developmental relationships and revealed a clear carotenoid metabolism grouping, coordinated expression of genes for anthocyanin and flavonoid enzymes and regulators and a cluster pattern of genes for chlorophyll catabolism consistent with functional and evolutionary features of the pathway.

Conclusion

The expression and phylogenetic characteristics of senescence-related genes allow a framework to be constructed of decisive events in the evolution of the senescence syndrome of modern land-plants. Combining phylogenetic, comparative sequence, gene expression and morphogenetic information leads to the conclusion that biochemical, cellular, integrative and adaptive systems were progressively added to the ancient primary core process of senescence as the evolving plant encountered new environmental and developmental contexts.

Background

We present evidence that the genes and metabolism representing the core of the process of senescence in modern land-plants can be traced back to primaeval unicellular photoautotrophs and that specific elaborations and regulatory mechanisms have been progressively added to the senescence program at decisive points during plant evolution. Senescence is a prominent characteristic of angiosperm morphogenesis and ecology. The term is applied at all levels of plant biology, from processes at the biome-wide scale that define the season of autumn in temperate climates [1, 2] through to the cellular dimension in the terminal development of individual tissues and organs [3, 4]. Plant senescence is commonly considered to be programmed [5] and several genes have been identified as essential for normal initiation and/or execution of the syndrome [6, 7]. Here we focus particularly on the metabolism of pigments in green tissues and the colorful lateral organs and structures of angiosperms [8, 9] and, based on comparative cell biology and protein sequence relationships across the evolutionary range of plants and their phototrophic ancestors, we reconstruct a progression of genetic events leading to the establishment of senescence as an essential element in the terrestrial plant life-cycle.

Evolutionary origin of photosynthetic pigments

A challenging question sometimes asked of biologists by children and physicists is: if leaves are solar panels, why are they green and not black [10–12]? This query in turn implies further questions about the nature of the earliest life-forms, how they used light as an energy-source, how the photosynthetic apparatus became organised and reorganised during evolution and (of particular relevance to the present paper) the evolutionary origin of pigment metabolism and its control in plant development. The light absorption profile of chlorophyll-based green membranes is characterised by a broad gap in the middle of the visible wavelength range. It has been proposed that the earliest phototrophs had light receptors that occupied the centre of the spectrum between the red and blue peaks of chlorophyll-based green membranes, and that chlorophyll evolved as a complementary pigment, filling in the missing wavelengths at the edges of the spectrum. According to this hypothesis the earliest photoautotrophs may have been similar to purple bacteria in using bacteriorhodopsin as the light-intercepting pigment-protein and bacteriochlorophyll at the reaction centre, driving a proton pump coupled to ATPase [10, 13, 14]. Phylogenetic analyses based on genes for photosynthesis as well as other markers such as rRNA and heat shock proteins are suggesting a number of alternative trajectories for the early evolution of pigment metabolism and organisation, but the widespread prevalence of horizontal gene transfer makes it difficult to agree the definitive scenario [13, 15]. We can be certain, however, that isoprenoid-derived pigments continued to carry out the function of intercepting light in the middle of the visible spectrum while chlorophyll adopted an increasingly prominent role in light harvesting as well as charge separation, leading to modern streptophytes including land plants.

Phylogenetic analysis of genes for Mg-tetrapyrrole biosynthesis indicates that anoxygenic photosynthetic organisms are ancestral to oxygen-evolving cyanobacteria and that the pathway in purple bacteria may be closest to that of the most ancient phototrophs [16]. By contrast to chlorophyll synthesis, the catabolic side of chlorophyll metabolism has not hitherto been examined in molecular phylogenetic terms, though a detailed survey of the enzymological properties of RCC reductase (RCCR) in a range of terrestrial species carried out by Hörtensteiner et al. [17] allowed some inferences to be drawn about the evolution of the pigment breakdown pathway in land plants. In the present paper we analyse the occurrence of genes for two key enzymes of chlorophyll catabolism (phaeophorbide a oxygenase – PaO – and RCCR; [18]), one post-translational regulator (Sgr; [19]) and two catabolite transporters (WBC23 and MRP2; [20]) in the lineages leading to modern angiosperms.

From photosynthetic eukaryotes to the invasion of land by plants

The significant early events in the evolution of the photosynthetic apparatus comprised diversification of reaction centre and antenna structure, based on chlorophyll and isoprenoids, and organisation and delineation of oxygenic and anoxygenic photosynthesis. According to Strother [21] 'the endosymbiotic origin of plastids...indicates that perhaps all of the original host cell lineages for the algae were phagotrophic...' Based on nutrient conditions and the evolution of food webs in the Palaeo- and Meso-proterozoic ocean [21, 22], it is reasonable to suppose that facultative autotrophy would have been the rule in early photosynthetic organisms. The formation of endosymbiotic relationships resulting in the evolution of plastid-containing eukaryotic cells probably occurred just once in the lineage, although there is evidence of several secondary events, including horizontal gene transfer, too [23–25]. It seems likely that the foundation endosymbiont was mixotrophic in character and that the earliest green eukaryotes were facultatively autotrophic [21, 26]. This is significant for the origin of plant senescence because aspects of the biochemistry and cell biology of plastids in terrestrial plants suggest that the transition from chloroplast to the plastids of senescing leaves (gerontoplasts) or fruit and floral parts (chromoplasts) is a change in status from auto- to hetero-trophy [27]. A strikingly senescence-like response during the change of trophic condition on removal of light and nitrogen source is seen in cultures of Auxenochlorella (formerly Chlorella) protothecoides. Cells turn yellow, accompanied by changes in protein and enzyme activities that resemble those of senescing mesophyll cells [28], and a red pigment is secreted into the medium [29–31]. Chemical analysis shows the pigment to be a bilin derived from chlorophyll by removal of phytol and Mg and by oxygenolysis of the methine bridge between pyrrole groups A and B [32]. This structure and its biochemical origin relate directly to the metabolic sequence for chlorophyll catabolism in angiosperms [33]. Similar production of chlorophyll catabolites has been reported in a Chlamydomonas mutant [34]. The luciferin of the bioluminescent dinoflagellate Pyrocystis lunula is also a bilin degradation product of chlorophyll, though in this case the chlorin macrocycle is opened between rings A and D [35]. A luciferin of similar structure is found in the krill organism Euphausia pacifica [35, 36], presumably originating in ingested phytoplankton. We conclude that at least part of the enzymic pathway of chlorophyll catabolism was present during the aquatic phase of plant evolution and this has led us to seek phylogenetic continuity between genes for trophic responses and pigment catabolism from unicellular algae to angiosperms.

Amongst the genes expressed during the auto- to hetero-trophic transition in Auxenochlorella is an amino acid permease, dee4 [28]. Transcripts of the Arabidopsis homologue of dee4, At1g44100, were identified in isolated nuclei of Arabidopsis phloem cells, suggesting a role in nitrogen translocation from sources to sinks [37]. Bleecker [38] considered nutrient salvage to be a driving force in the evolution of the plant senescence program. Nutrient recovery implies structural differentiation into source and destination tissues, with a transport system between them [39, 40]. To assess their relevance to the evolution of nutrient relations during senescence, we carried out protein sequence tree and expression analyses of dee4 and At1g44100 respectively. We also analysed a NAC gene with a regulatory function in organ development and senescence. Uauy et al. [41] isolated NAM-2, a NAC transcription factor, by positional cloning and demonstrated a regulatory role in nutrient redistribution from senescing leaves to developing grain in cereals. AtNAP, a NAM-like factor, expression of which is closely associated with the senescence process of Arabidopsis rosette leaves [42], was included in our analyses of expression and taxonomic distribution.

From endosymbiont origins of photosynthetic eukaryotes to angiosperms is a sequence of major events including differentiation of two- and three-dimensional body plans and emergence from the aquatic environment onto land [43]. Cell specialisation, intra- and inter-cellular signalling and elaboration of body plan are regulated by homeoprotein transcription factors and the integration of senescence into morphogenetic and adaptive programs is likely to be apparent in the phylogenetic record. Prominent amongst the regulators of signalling networks activated during plant senescence are transcription factors of the Wrky family. Sequence analysis of Wrky53, the most intensively studied of the senescence-related signal-transduction components [44], suggests the point during evolution at which such factors may have became available for networking with the core senescence program.

Appearance of coloured organs in land-plants driven by coevolution

As well as promoting diversification of structure, and adaptation to terrestrial ecological niches, achieving landfall eventually exposed plants to new biotic and abiotic interactions. The ancient role of carotenoids as quenchers of excess light energy intercepted by photosystems [45] acquired new significance in environments where mismatch between photon flux, temperature and the availability of water and nutrients is the rule and a strongly oxidising atmosphere supports the uncontrolled propagation of damage by reactive chemical species. Coevolution between plant pigments and the visual systems of insects and other animals introduced a new role for carotenoids as signalling molecules, advertising to pollinators, to dispersers or to herbivores and other predators [46, 47]. Investigations of the phylogenetic origin of carotenoids in relation to colour changes in senescence and ripening are frustrated by the ubiquity and antiquity of the corresponding genes and the inordinate degree of horizontal gene transfer apparent in the evolutionary record [48, 49]. Nevertheless, we have identified four genes associated with the subcellular organisation and developmental regulation of carotenoid metabolism in senescing and ripening tissues that provide insights into how this aspect of senescence may have evolved. Lipid bodies (plastoglobules) accumulate during the transition from chloroplast to gerontoplast or chromoplast [50]. They have been observed in the plastids of Auxenochlorella during the switch from auto- to hetero-trophy [51]. Plastoglobule-like lipid bodies are very widely distributed amongst prokaryotes [52]. Existing and newly-synthesised carotenoids concentrated in plastoglobules are responsible for the colours of fruits and senescing leaves in some species [8, 53, 54]. We examined the distribution of sequences similar to PAP-Fibrillin, a plastoglobule protein with functions in plastid development and stress responses [55].

Carotenoid metabolism in senescence and ripening is closely integrated with the transdifferentiation of chloroplasts into gerontoplasts and chromoplasts. Chromoplast development is regulated by Or (dnaJ-type) chaperone/transcription factors [56], two of which we included in our analyses. Carotenoid derivatives also have hormonal and cell signalling roles in plants. Genes of the Carotenoid Cleavage Dioxygenase (CCD) group are known to participate in regulation of a number of aspects of development, including leaf senescence [57]. Accordingly, we carried out an analysis of the protein sequences related to that of senescence-associated CCD8.

Another major class of pigments in land-plants comprises products of phenylpropanoid metabolism [58]. Anthocyanins are considered to have adaptive value both in stress responses and in signalling to beneficial or predatory animals, and there is an extensive literature on the significance of the range of possible anthocyanin functions in senescing leaves [8, 59]. By contrast with carotenoids, which are highly hydrophobic and confined to the lipophilic environment of plastid membranes and plastoglobules, anthocyanins are water-soluble and accumulate in the cell vacuole [60]. The behaviour of the two classes of pigment during senescence is indicative of the state of their respective subcellular compartments. The anthocyanin responsible for the red colour of autumn leaves is commonly cyanidin-3-glucoside. The final step in its synthesis is catalysed by anthocyanidin 3-O-glucosyltransferase, encoded by the Bronze1 (Bz1) gene [61]. A myb transcription factor C1 regulates expression of the anthocyanin pathway [62]. We included Bz1 and C1 in our evolutionary survey of senescence-associated genes. As a senescence-related marker for the evolution of specific tonoplast membrane functions we also selected MRP2, a chlorophyll catabolite transporter [63]. The proteome of the vegetative plant cell vacuole is consistent with a lytic role in defence and tissue differentiation [64]. Transposon insertional mutagenesis of See2, a gene of maize encoding a vacuolar processing endopeptidase (VPE) has a range of effects on plant development, including delay of senescence [65]. We carried out an evolutionary survey of See2 and confirmed senescence-related expression of the Arabidopsis homologue At4g32940.

Relationships between the ontogeny and phylogeny of plant senescence

As well as analysing the distribution across plant and microbial taxa of the selection of senescence-related protein sequences described above, we applied the clustering tools in Genevestigator to an extended group of Arabidopsis genes with likely functions in senescence, with the object of identifying common regulatory units. It is also of interest to seek evidence for homologies between terminal physiological changes in different organs, which in turn would begin to answer the question as to whether there is a common senescence program which operates, with variations, in all plant cells, tissues and organs.

Results

Gene expression profiles, protein structures and molecular phylogenies are discussed in relation to chlorophyll catabolism and its control, homeoprotein transcription factors with known functions in senescence, metabolite transporters, carotenoid metabolism and anthocyanin biosynthesis (Table 1). Many senescence-associated genes have been described but the physiological functions of relatively few of them have been established, and fewer still have been demonstrated to be necessary for initiation or execution of the syndrome [6, 7, 9]. Those listed in Table 1 have been shown, by mutation, transgenic manipulation or biochemical analysis, to have defined roles in senescence. The tools of phylogenetic analysis were applied to amino acid rather than nucleotide sequences to avoid difficulties with codon usage across the very wide taxonomic scale considered in this study, though where a particular taxon yielded no hit when a critical protein was BLASTed, we checked for the existence of corresponding pseudogenes or other untranslated nucleotide sequence. Since we are principally interested in the first detectable occurrence and subsequent evolutionary continuity of particular protein sequences, we have not explicitly distinguished between orthologues and paralogues, a known issue when establishing phylogenetic and functional relationships [66]. The relative paucity of sequence data for many of the species examined is an acknowledged limitation in this study but, with a few exceptions discussed in detail, the analyses reported here make a consistent and biologically plausible picture.

Table 1 Proteins and their corresponding Arabidopsis genes with known functions in senescence, selected for phylogenetic and expression analyses.

In most cases the trees are presented in unrooted radial form and are based on global alignments of protein sequences from one or two representative genera of each major taxon. Tree topologies for individual proteins often deviated from those predicted by the evolutionary relationships between the corresponding taxa. The non-congruence of gene and species trees is a well-described feature of molecular phylogenies and indeed is statistically to be expected [67]. Bootstrap values and branch lengths with their error estimations are presented in Additional Files 1 and 2 respectively. The range of senescence-related sequences (Table 1) is subsequently extended (Figure 1) for cluster analysis to reveal units of coexpression and developmental relationships amongst plant structures with respect to their senescence programs.

Figure 1
figure 1

Arabidopsis genes with functions and expression patterns related to senescence. The list includes the genes in Table 1 plus senescence-enhanced genes identified by a search of the TAIR database. Genes are colour-coded according to broad functional categories, where known or inferred.

Chlorophyll catabolism

The observation of bilin products of chlorophyll catabolism in unicellular algae focuses attention on the enzymes responsible for opening the chlorin macrocycle, namely PaO and RCCR. Figure 2 is a Genevestigator heat map representation of the relative expression of all the Arabidopsis genes whose protein trees are reported in the present paper, including RCCR and PaO (AGI accession numbers respectively At4g37000 and At3g44880). Expression of RCCR is not particularly strong in senescing leaves. It gives a noticeably weak signal in roots and other non-green structures like pollen. PaO is also largely undetectable in non-green tissues but exhibits clear differential expression during development, being most prominent in senescing leaf, cauline leaf and sepal.

Figure 2
figure 2

Genevestigator Anatomy expression profiles of 15 senescence-related Arabidopsis genes (L to R: RCCR, PaO, Sgr1, WBC23, Wrky53, NAP, dee4, Fibrillin, CCD8, OrI, OrII, Bronze1, myb C1, MRP2, See2). Heat maps represent expression level in Affimetrix® arrays as colour intensity computed in Genevestigator according to the MAS5 normalization algorithm [192].

Sequence similarities to the protein corresponding to RCCR gene At4g37000 (UniProt accession number for the barley sequence Q8LDU4; [68]) were analysed. Figure 3 shows the relationships between these proteins from representative genera in tree form. Three major groupings are apparent. Angiosperms form one clade, with monocots and dicots resolved within it. Non-angiosperm land plants are grouped together. Cyanobacteria form a third clade, characterised by relatively high E scores, although examination of the sequences for conserved domains identified these cyanobacterial proteins as members of the RCCR family. For example SMART analysis recognised the RCCR conserved domain in the Anabaena protein, recording an extremely low E score of 4.80e-159. Protein and DNA BLAST searches, both of general databases and of the JGI Chlamydomonas reinhardtii site, failed to record a single hit within the green algae. The significance of the discontinuity in the distribution of RCCR across phyla is considered in the Discussion section.

Figure 3
figure 3

Unrooted radial bootstrapped protein sequence tree for RCCR based on PHYLIP analyses using Arabidopsis RCCR (UniProt NP_195417) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1a). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2a gives branch length values and confidence limits from the PROML outfile.

The tree representation of proteins similar to the Arabidopsis PaO gene At3g44880 (UniProt accession number NP_190074; Figure 4) reveals that the land plants form a group of closely related sequences (E value of 0). The cyanobacteria, represented by Cyanothece and Lyngbya, make a distinct cluster, the green algae Chlamydomonas and Ostreococcus make another and the proteobacteria Phenylobacterium and Ralstonia are more remotely related. The protein sequences of eukaryotes and cyanobacterial species include conserved Rieske Fe-S and PaO motifs (Figure 5). The Ralstonia sequence has a Rieske Fe-S domain near the N terminus but no PaO motif (Figure 5c). This is generally true of bacterial sequences with similarity to PaO. The Ralstonia protein is recognised as a member of the HcaE phenylpropionate dioxygenase multidomain group, which may be close to the evolutionary source of the PaO gene of phototrophs. It is striking that the cyanobacterial branch has a higher support value and their sequences have a higher similarity to those of land plants than do those of the eukaryotic green algae (Figure 4).

Figure 4
figure 4

Unrooted radial bootstrapped protein sequence tree for PaO based on PHYLIP analyses using Arabidopsis RCCR (UniProt NP_190074) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1b). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2b gives branch length values and confidence limits from the PROML outfile.

Figure 5
figure 5

(a) Comparison of PaO structures in representative genera by BLASTP alignment. Shading density indicates % sequence similarity: black = 100%; dark grey = 80%; light grey = 60%. Conserved motifs are shown. Both the Rieske and PaO domains are present in the Arabidopsis sequence (b) and other species except that of the b-proteobacterium Ralstonia (c), in which the PaO motif is missing.

Stay-green (Sgr) is a post-translational regulator of chlorophyll-protein complex degradation [19, 69, 70]. The Genevestigator profile shows a tightly-specified expression pattern, with highest transcript abundance in senescent leaf and significant signal in floral parts (Figure 2). The sequence of the protein corresponding to Arabidopsis Sgr1 gene At4g22920 (UniProt accession number NP_567673) was aligned with similar proteins from other species and the results were used to generate the tree shown in Figure 6. Dicots, monocots and non-flowering land plants form three distinct groupings, with low E values that gradually increase with evolutionary distance from the reference species. Between the land plant and green algal clades there is a marked step up in E value. Grouped with the algae is a well-supported outlying clade, comprising a number of firmicutes (clostridia and bacilli). We further analysed the features of the Sgr-like genes in these species to determine whether Sgr may be indicative of a photosynthetic Clostridium ancestor, possibly something like the firmicute Heliobacterium. The rRNA [71, 72] and Sgr trees are similar (Figure 7), but not congruent, and only 14 out of a total of 283 genomes searched (BLASTP at the JGI IMG site) contained genes encoding Sgr homologues. These were restricted to the genera Clostridium (11 of 44 genomes) and Bacillus (3 of 41 genomes) and there was no obvious conservation of genome context of Sgr between these organisms. We therefore conclude that Sgr is either an ancestral feature, which has been lost in all but a few remaining lineages; or, more likely, a feature that has been acquired more than once via horizontal gene transfer (Figure 7).

Figure 6
figure 6

Unrooted radial bootstrapped tree for plant and bacterial Sgr-like proteins based on PHYLIP analyses using Arabidopsis Sgr1 (UniProt NP_567673) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1c). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2c gives branch length values and confidence limits from the PROML outfile.

Figure 7
figure 7

Structures and phylogenies of Clostridium Sgr proteins. (a) Comparison of Sgr structures in representative plant and Firmicute genera by BLASTP alignment. Shading density indicates % sequence similarity: black = 100%; dark grey = 80%; light grey = 60%. (b) Unrooted tree of phylogenetic relationships between Clostridium spp based on sequence alignments of ribosomal RNA [71, 72]. (c) Protein tree of Sgr-like sequences from clostridia.

Aubry [20] identified WBC23, a transporter in the plastid envelope responsible for removing the linear tetrapyrrole product of the PaO-RCCR reaction. The corresponding gene, At5g06530, is strongly expressed in senescing leaves (Figure 2). ATP-dependent membrane transporters are very widely distributed across all eukaryotes [73], as evidenced by hits, with low E values, on the cnidarian Nematostella and the choanoflagellate Monosiga when BLASTed against Arabidopsis WBC23 (UniProt code NP_850781; Figure 8). Land plants group together, with clear discrimination into monocots, dicots and mosses. Chlamydomonas and Ostreococcus are associated at some distance from the land plants and each other. At higher E values, but still significantly less than chance, BLAST also identifies alignments with proteins from representatives of various insect groups.

Figure 8
figure 8

Unrooted radial bootstrapped protein sequence tree for family of ATP-dependent bilin transporters based on PHYLIP analyses using Arabidopsis WBC23 (UniProt NP_850781) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1d). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2d gives branch length values and confidence limits from the PROML outfile.

Cell and tissue specialisation

In a study of genes expressed during Arabidopsis leaf senescence, Hinderhofer and Zentgraf [44] identified Wrky53 as a transcription factor with a major coordinating role in senescence regulation. Phylogenetic analysis of the Wrky family by Zhang and Wang [74] established that Wrky53 is a member of Group 3 that appeared before divergence of monocots and dicots but later than diversification of the bryophytes 160 – 330 million years ago (mya). Wrky53 is expressed exclusively in mature and senescing leaves (Figure 2), suggesting a role beginning in the earliest phase of senescence. Sequence similarities to the protein corresponding to Wrky53 gene At4g23810 (UniProt accession number NP_194112) were analysed and the resulting tree topology was in agreement with that of Zhang and Wang [74], with the exception of the appearance of a single hit in Chlamydomonas (Figure 9). The E score for this protein is 2e-11 and analysis of conserved domains reveals two Wrky superfamily sequences with E-values of 9.7e-33 and 1.26e-35 respectively (Figure 9). We conclude that this represents a true algal Wrky53 homologue and is evidence that the appearance of Group 3 Wrkys may have been significantly earlier than the 330 mya originally inferred.

Figure 9
figure 9

Unrooted radial bootstrapped protein sequence tree for Wrky53 based on PHYLIP analyses using Arabidopsis Wrky53 (UniProt NP_194112) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1e). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2e gives branch length values and confidence limits from the PROML outfile. Conserved domains in Chlamydomonas Wrky53-like protein are shown (low-complexity sequence indicated in pink).

During senescence metabolites are relocated within the cell, from tissue to tissue and from source to sink organs. Uauy et al. [41] isolated a NAC transcription factor that corresponded to a major QTL for nutrient mobilisation from vegetative tissues during cereal grain-fill. Guo and Gan [42] also identified a senescence-associated NAC factor (AtNAP, AGI At1g69490) in Arabidopsis. Figure 2 reiterates the expression patterns reported by Guo and Gan [42], confirming strong upregulation in senescing leaves and high levels in cauline leaf, sepal and petal. Alignment of the AtNAP protein (UniProt NP_564966) showed similarity to the cereal sequence. The tree based on protein sequence alignments is presented in Figure 10. Exhaustive protein and nucleotide BLAST searches support the conclusion that NAP- and NAM-like proteins are exclusively found in land plants.

Figure 10
figure 10

Unrooted radial bootstrapped tree for NAP-like proteins based on PHYLIP analyses using Arabidopsis AtNAP (UniProt NP_564966) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1f). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2f gives branch length values and confidence limits from the PROML outfile.

A gene encoding an amino acid permease (dee4, UniProt Q9XFY8) is upregulated when Auxenochlorella protothecoides is transferred from auto- to hetero-trophic culture conditions; the Arabidopsis homologue (At1g44100) is also upregulated in senescing leaves (Figure 2) and has been identified in phloem [37]. Similar proteins to dee4 occur in angiosperms, mosses and conifers but are not present in taxa that diverged before the appearance of green algae (Figure 11).

Figure 11
figure 11

Unrooted radial bootstrapped protein sequence tree for senescence-related amino acid permeases based on PHYLIP analyses using Auxenochlorella dee4 (UniProt Q9XFY8) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1g). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2g gives branch length values and confidence limits from the PROML outfile.

Gerontoplasts and chromoplasts

Plastids differentiate, dedifferentiate and morph from one type to another in a network of developmental transitions (Figure 12; [75–78]). Orange and red carotenoids are unmasked, concentrated or synthesised de novo in gerontoplasts and chromoplasts during senescence and ripening [8] and are sequestered in plastoglobules that become abundant as chloroplasts transdifferentiate [55, 79]. We examined the expression and sequence relationships of Arabidopsis PAP-Fibrillin (UniProt NP_192311, AGI At4g04020), a plastoglobule-associated protein [55]. Genevestigator profiling shows expression to be absent from roots and to be prominent in cotyledons, floral parts and leaves (Figure 2). The protein tree for Fibrillin (Figure 13) reveals dicot and monocot clusters and a conifer branch with a comparatively low bootstrap value. Eukaryotic unicellular green algae and cyanobacteria make a distinct separate group. The Ostreococcus sequence with a high E score appears to be a chimeric protein comprising two additional FHA (forkhead-associated, putative nuclear-signalling) domains, separated from each other and the Fibrillin domain by regions of low complexity (Figure 13).

Figure 12
figure 12

The plastid development network (based on [75, 78]).

Figure 13
figure 13

Unrooted radial bootstrapped tree for Fibrillin-like proteins based on PHYLIP analyses using Arabidopsis PAP-Fibrillin (UniProt NP_192311) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1h). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2h gives branch length values and confidence limits from the PROML outfile. The schematic diagram shows conserved domains within the Ostreococcus protein, including FHA (forkhead-associated) motifs, a coiled-coil sequence (green) and regions of low complexity (pink).

The gene Carotenoid Cleavage Dioxygenase8 (CCD8; UniProt NP_195007, AGI At4g32810) has a regulatory function in branch production, root growth, flower development and leaf senescence [57]. It is mostly expressed in shoots, particularly hypocotyl tissues, stem, node and senescent leaf (Figure 2). The protein is highly conserved from proteobacteria (eg Coxiella) to angiosperms (Figure 14) and a search for conserved domains establishes that CCD8 and its homologues are members of the RPE65 retinal pigment epithelial membrane protein superfamily which includes neoxanthin cleavage enzymes in plants and lignostilbene-α, β-dioxygenase in bacteria.

Figure 14
figure 14

Unrooted radial bootstrapped protein sequence tree for CCD8 based on PHYLIP analyses using Arabidopsis CCD8 (UniProt NP_195007) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1i). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2i gives branch length values and confidence limits from the PROML outfile.

The Or genes of Brassica encode proteins of the dnaJ chaperone class that play an important role in the differentiation of chromoplasts [56, 80]. At5g61670 is the Arabidopsis gene corresponding to the cauliflower Or protein ABH07405 (UniProt), and is referred to here as OrI. It is expressed in floral parts, cotyledons and in cauline and senescent leaves (Figure 2). At5g06130 and NP_851031 are the AGI and UniProt codes for a second Arabidopsis Or-like dnaJ gene and protein respectively, designated OrII. The expression profile of this gene is broadly similar to that of At5g61670 (Figure 2), with somewhat stronger expression in senescing leaf, petal, stamen and pedicel. Trees of sequence relationships for OrI and OrII are presented in Figure 15. Both proteins are highly conserved from prasinophyceans to angiosperms. Prokaryotes gave no significant BLAST hit. The occurrence of Or-like genes in plants that do not differentiate chromoplasts, and their upregulation in a range of tissues including floral parts and senescing leaves (Figure 2), suggest that they may have a more generic function in plastid transdifferentiation.

Figure 15
figure 15

Unrooted radial bootstrapped trees for two dnaJ chaperone-like proteins based on PHYLIP analyses using (a) OrI of Brassica (UniProt ABH07405) and (b) Arabidopsis OrII (UniProt NP_851031) as the reference sequences. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1j and 1k). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2j and 2k give branch length values and confidence limits from the PROML outfiles.

Vacuoles

Anthocyanin production and leaf yellowing are associated with senescence responses to nitrogen limitation, but are genetically independent [81]. We analysed the expression and protein sequence relationships of two genes with functions in anthocyanin biosynthesis. Bronze1 (Bz1) encodes the glycosyl transferase that synthesises the red vacuolar pigment cyanidin-3-glucoside. Zea Bz1 (UniProt P16167) was aligned with similar sequences; the homologous Arabidopsis gene is AGI At5g17050. The anthocyanin pathway is under the control of myb transcription factor C1 (UniProt accession number of Zea C1 P10290; AGI code of the homologous Arabidopsis gene At4g34990). Both genes are widely expressed in different organs and tissues (Figure 2). Bz1 is most prominent in senescing leaf. C1 is also expressed in senescing leaf but is most strongly represented in stigma and hypocotyl tissues. The Bz1 tree reveals the occurrence of homologues in angiosperms, conifers and mosses (Figure 16). At an E value of 1e-08, only a single BLAST hit was recorded for a non-tracheophyte. The Chlamydomonas protein has a UDP glycosyl transferase region towards the C terminal end but low similarity to Bz1 elsewhere in the polypeptide. A BLAST search against the Chlamydomonas protein identifies other plant UDP glycosyl transferases, but E scores are high, and the algal protein aligns more closely with a range of nucleoside glycosyl transferases from metazoa. We conclude that this algal protein has a distant structural relationship to tracheophyte Bz1-like sequences but functions as a transferase in a different metabolic pathway from that leading to anthocyanin synthesis. It follows that the anthocyanin pathway appeared relatively late in the evolution of senescence processes, around the time of cryptogam diversification. By contrast, genes related to C1, the myb transcription factor responsible for regulating expression of anthocyanin synthesis, are well represented throughout plant taxa, from Zea to Ostreococcus (Figure 17). BLAST also identifies similar genes in the actinomycete Frankia and other bacteria.

Figure 16
figure 16

Unrooted radial bootstrapped protein sequence tree for a family of anthocyanidin 3-O-glucosyltransferases based on PHYLIP analyses using Zea Bronze 1 (UniProt P16167) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1l). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2l gives branch length values and confidence limits from the PROML outfile.

Figure 17
figure 17

Unrooted radial bootstrapped tree for plant and bacterial myb C1-like proteins based on PHYLIP analyses using Zea C1 transcription factor (UniProt P10290) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1m). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2m gives branch length values and confidence limits from the PROML outfile.

The final products of chlorophyll catabolism are deposited in the vacuole. They are moved across the tonoplast membrane by ATP mediated transporters of the MRP type [63]. MRP2 (AGI At2g34660) is expressed in senescent and cauline leaves as well as sepal, stamen and root tissues (Figure 2). Multidrug-resistance transporter proteins like AtMRP2 are very widely distributed across plants and animals. The tree of the corresponding protein (UniProt NP_181013) identifies a group of closely-related structures within angiosperms, mosses and green algae. Interestingly the MRP2-like protein of the predatory bacterium Bdellovibrio aligns to Arabidopsis with a lower E score than the cyanobacterial homologue (Figure 18a). BLAST also identifies a more distantly-diverged, but still highly similar, branch consisting of protozoa (Placozoa, Choanocyte) and vertebrates. The algal and angiosperm MRP2-like sequences have similar domain structures (Figure 18b, c).

Figure 18
figure 18

(a) Unrooted radial bootstrapped protein sequence tree for ABC-type transporters based on PHYLIP analyses using to Arabidopsis AtMRP2 (UniProt NP_181013) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1n). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2n gives branch length values and confidence limits from the PROML outfile. Conserved domains in (b) Arabidopsis and (c) Ostreococcus MRP2 sequences. AAA, ATPase region; blue block, transmembrane segment; pink block, low complexity region.

As well as representing a destination for pigments and their catabolites, the vacuole has complex functions in intracellular turnover, lytic processes and defence against biotic invasion. See2 (UniProt CAB64545) is the product of a Zea gene originally identified amongst those upregulated during flag leaf senescence [82]. It is a representative of the VPE group of protein-processing legumains and may function as part of a proteolytic cascade, rather in the manner of the caspases involved in apoptosis [83]. The Arabidopsis homologue (At4g32940) is most strongly expressed in senescent and cauline leaf, with some significant signal in petal, sepal, xylem and cork (Figure 2). BLAST identified similar proteins in angiosperms, conifers, mosses and green algae (Figure 19). A distantly diverged group comprises marine myxobacteria, diverse vertebrate and invertebrate taxa.

Figure 19
figure 19

Unrooted radial bootstrapped tree for VPE-like proteins based on PHYLIP analyses using Zea See2 (UniProt CAB64545) as the reference sequence. Representative genera are shown, together with E values, indicating the numerical expectation that the alignment with the reference sequence has arisen by chance. Support values (in red) were determined by bootstrap analysis (Additional file 1o). Branch lengths are proportional to number of amino acid substitutions per position. Additional file 2o gives branch length values and confidence limits from the PROML outfile.

Cluster analysis of senescence-related genes

In the evolution of the senescence program, genes would be expected to show a degree of coordinate expression in groupings that reflect phylogenetic trends. Similarly, a consequence of the differentiation of plant tissues and organs from ancestral structures is the sharing of common senescence programs. We used the Genevestigator Clustering Analysis tool to look for relationships between a collection of Arabidopsis genes comprising those listed in Table 1 and a further group that includes the results of a search of the TAIR database for senescence-related genes. Figure 1 lists a total of 63 genes, each of which is assigned to a colour-coded broad functional grouping (no colour indicates an unknown, unclassifiable or hypothetical protein). The Genevestigator heat map shown in Figure 20 clusters these genes according to profile and tissue, in the order left to right in which they are listed in Figure 1. The distribution of colour codes indicates blocks of genes with common functions and coordinated expression patterns. Tissues with similar complements of expressed genes are hierarchically arranged close to each other. The gene and tissue associations presented in Figure 20 are suggestive of developmental and evolutionary relationships that are considered further in the Discussion section below.

Figure 20
figure 20

Genevestigator heat map of Arabidopsis genes with functions and transcription patterns related to senescence, clustered by tissue and expression profile according to transcript abundance. Heat maps represent expression level in Affimetrix® arrays as shading intensity computed in Genevestigator according to the MAS5 normalization algorithm [192]. Genes are colour-coded by broad functional categories and listed in order as shown in Figure 1.

Discussion

The present account of the evolutionary origins of plant senescence seeks to follow the major features of the syndrome in angiosperms through the lineage of organisms leading from the earliest phototrophs, as revealed by analysing sequence relationships among key proteins. Additional File 3 summarises innovations in the evolutionary time-line that are relevant to the construction of a senescence program integrated into plant development and adaptation, and represents a framework for the discussion that follows. Figure 21 relates this time-line to the phylogenetic origins of the senescence-related genes analysed in the present paper. The foundation of the argument developed here is that a core of primary metabolic processes concerned with building and maintaining the photoreceptors of the earliest prokaryotes has been invested with successive layers of adaptive metabolic and regulatory circuitry as factors in the evolution of complexity in trophic status, life-cycle, anatomy, morphology, adaptation to new environments and biotic and abiotic challenges.

Figure 21
figure 21

Proposed evolutionary timeline for senescence-related genes discussed in the present study.

Pre-endosymbiotic phase

The Archaeozoic and Mesoproterozoic eons lasted from 3800-1000 mya. Although the detailed phylogeny of this period is difficult to infer from molecular analysis of contemporary prokaryotes [84], it is probable that the fundamental design of the photosynthetic apparatus was established early, in the form of pigment-protein photoreceptors driving an ATP-generating proton pump [13]. The ubiquity of genes for isoprenoid synthesis [49] is evidence of their ancient and fundamental role in light interception and photochemical quenching [45, 85, 86]. Genes for (bacterio)chlorophyll biosynthesis have been used to create well-rooted trees leading back to the putative earliest photosynthetic organisms [15]. Assembling complexes requires translational regulation [87] and efficient facilitated shape-matching [88], to coordinate the supply of structural elements. In modern autotrophs a major factor in achieving stoichiometry is post-transcriptional culling of excess unmatched components [89–91]. Thus in mutants deficient in pigment or protein components of photosystems, most of the genes responsible for building the complexes are transcribed and translated normally, but because they cannot make stable assemblages the proteins are degraded by fidelity-testing peptide hydrolases. For example, Fu et al. [92] showed that the cyclophilin AtCYP38 is required to stabilise thylakoid complexes in Arabidopsis: mutants deficient in CYP38 make normal amounts of photosystem structural components but these fail to assemble correctly and are turned over. Significantly, the CYP38 gene (At3g01480) is strongly down-regulated in senescing compared with mature and juvenile Arabidopsis leaves (data not shown). Xiong [15] described how a functional photosynthetic apparatus could have evolved in the earliest prokaryotes by multi-staged recruitment of reaction centre and antenna proteins with previously diverse evolutionary histories and functions. We consider that disassembly of complexes is correspondingly hierarchical in origin. Phylogenetic analyses of plant proteases [93] and of proteasome components and organisation [94] shows convincingly that the catabolic machinery for protein turnover is as ancient as the Chlorobacteria and Proteobacteria at the base of the evolutionary tree.

Charge-separation is hazardous for pigment-protein reaction centres and here again there will be a requirement for turnover systems to maintain functioning structures [95]. Genes specifying cellular systems for turning over photoreceptor complexes may therefore be expected to be ancient and we hypothesise that the distribution of protein sequences implicated in some features of angiosperm senescence represents evidence for continuity back to the pre-endosymbiotic origins of green plants. Studies of their function in chloroplast senescence suggest that Sgr and PaO are components of a machine that destabilises thylakoid pigment-proteolipid complexes by removing chlorophyll and exposing the associated proteins to attack by peptide hydrolases [19, 65, 96]. Plastoglobules play a part in this process, receiving hydrophobic products released from dismantled membranes, including carotenoids, phytol from chlorophyll and products of lipolysis [55, 97]. The structural similarity between the PaO- and Fibrillin-like proteins of prokaryotic phototrophs and the corresponding proteins of angiosperms (Figures 4, 13) supports the proposition that these components of the machinery for unpacking pigment-protein complexes are ancient.

Endosymbiotic associations in the origin of eukaryotic plants

The origin of the chloroplast as a cyanobacterium-protozoan symbiosis occurred, according to most estimates, around 1500 mya or earlier [98], though there is an alternative view placing it as late as 600 mya [23]. The primary event soon gave rise to two lineages, the glaucophytes and the green plant/red algal group. Mixotrophy and opportunist organotrophy were probably well established in the earliest phototrophic prokaryotes and, in response to increasing competition for nutrients, vestiges of facultative trophic behaviour have persisted throughout the period from the enslavement of the endosymbiotic ancestors of eukaryote organelles to the plastid differentiation network of modern angiosperms (Figure 12; [24]). The metabolism of green cells of angiosperms changes from auto- to hetero-trophic as the photosynthetic apparatus is remodelled during the transdifferentiation of chloroplasts into chromoplasts or gerontoplasts. Light is a declining source of energy in these cells and respiratory processes become increasingly significant [99]. Fruit ripening in many species, and leaf senescence in a few, is accompanied by a respiratory burst (climacteric; [100]). The amino acid products of protein recycling are subjected to transamination, deamination and amide synthesis to provide carbon skeletons for respiration while salvaging nitrogen in a transportable form [101]. A gene encoding an amino acid permease is upregulated when Auxenochlorella protothecoides is transferred from auto- to hetero-trophic culture conditions; the Arabidopsis homologue is also upregulated in senescing leaves (Figure 2). The exact identity of the protease(s) responsible for releasing amino acids from photosynthetic proteins is not agreed [102]. A number of plastid-located proteases are known [90], amongst which at least four (ClpP3, ClpD, FtsH7 and FtsH8) are upregulated in senescence (Figure 20). BLAST searches show that highly homologous sequences to all four proteins are widely distributed across taxa including phototrophic and non-phototrophic bacteria. There is evidence that Sgr is essential for making chlorophyll-binding proteins available to proteases [9, 69, 89]. Sgr is not detectable in photosynthetic prokaryotes, although Sgr-like sequences are present in some firmicutes, possibly as the result of horizontal transfer (Figures 6, 7). We suggest that Sgr was recruited to the system for dismantling pigment-proteolipid complexes when subcellular topology became more complex with the appearance of the eukaryotic condition. This in turn suggests that protein mobilisation and colour change in tracheophyte senescence have their origins in the intraplastidic turnover of pigment-proteins related to regulated assembly of macromolecular complexes and trophic flexibility.

An alternative, or additional, mechanism for proteolysis in senescence involves exchange of protein substrates and/or proteases between plastids and vacuoles in as-yet uncharacterised vesicles [103–105]. This would assign a role in senescence to vacuolar proteases such as SAG12 ([106]; Figure 20) and to components such as See2 (Figures 2, 19) that may trigger proteolytic cascades [83, 107]. Martinez and Diaz [93] have studied the evolution in plants of C1A papain-like cysteine endopeptidases (SAG12 is one of these) and the C13 VPE-like legumains (includes See2) that may activate them. C1A proteases are present in angiosperms, mosses, Volvox, Chlamydomonas and Ostreococcus. Martinez and Diaz [93] report that Ostreococcus has no C13 VPE sequence, but the NCBI BLAST search reported here has produced two hits: one (UniProt Q00VF5) from the published study of O tauri by Derelle et al. [108], the other from the US DOE Joint Genome Institute O lucimarinus database (B Palenik et al. unpublished submission 2007). In both cases the E value for alignment with See2 is better than 1e-13 and for the C13 peptidase superfamily conserved domain E scores were of the order of 1e-80. We conclude that VPE-like proteases are present in Ostreococcus (Figure 19).

An ancient function of the cell vacuole is to act as an autolytic body in cell death. 'Metacaspases' (a name for vacuolar proteases that may function like the caspases of apoptotic animal cells but which have little or no homology to them) have been implicated in a cell death program in diatoms [109]. It is doubtful that autolysis plays a part in the initiation and execution of the senescence program in angiosperms [110] but it may be critical in the terminal phase of tissue development when senescence is ending – for example in vascular tissue where the cells are dead at maturity [111, 112] or as a cauterisation measure in senescent mesophyll [106, 110]. We speculate that the progression from Ostreococcus to Volvox to Chlamydomonas shown in Figure 19 may signify the appearance and development of a capacity for proteolytic processing by See2-like VPEs during the period of green alga evolution associated with a new role for the vacuole in organelle protein turnover.

The auto- to hetero-trophic switch induced in Auxenochlorella protothecoides by light and nitrogen limitation is associated with excretion of bilin products of chlorophyll catabolism [28]. It is not known whether such catabolites have antibiotic or allelopathic properties, though they are certainly highly photodynamic. Also unknown is the transporter that moves bilins across the cell membrane. ATP-dependent transporters that export chlorophyll catabolites from plastids to cytosol and from cytosol to vacuole have been identified in Arabidopsis and the corresponding genes (At5g06530 encoding the plastid envelope transporter, At2g34660 encoding one of the tonoplast MRP2 transporters) are upregulated during senescence (Figure 2). BLAST identifies algal sequences with similarity to AtMRP2, and comparisons of domain structure suggest common functions (Figure 18). However the sequence of the AtMRP2 protein has higher similarity to that of AtMRP4 than to its Ostreococcus homologue. Both AtMRP2 and AtMRP4 proteins have the conserved domain structure of full-molecule ABC transporters, but the former is located in the tonoplast whereas the latter is a plasmalemma component [73]. Such an MRP in algae, with similarities to both MRP2 and MRP4, would have the expected properties of the transporter responsible for ejecting chlorophyll catabolites from the cell. Subsequently, it is conceivable that tonoplast and plasmalemma transporter localisations and functions diverged as the participation of the vacuole in chlorophyll catabolism became established as part of a novel physiological module added downstream of bilin formation when multicellular plants emerged onto the land [113].

Development of multicellularity

The Ordovician-Silurian period (450-410 mya) saw the rise of the charophytes, which are considered to link the unicellular chlorophytes with the first land-plants. The lack of sequence data for members of the Charophyta is a major gap in the phylogenetic record, but it is possible to identify features within charophyte algae likely to be relevant to senescence evolution, notably the elaboration of multicellular morphology, branching, 3-D cell division and interconnection of cells by plasmodesmata [43]. A key genus is Coleochaete, which includes branching filaments and thallus-like parenchymatous forms. Molecular phylogenies based on chloroplast and ribosomal sequences put this group close to the evolutionary line between unicellular chlorophyceans and bryophytes [114].

The presence of lignin-like compounds in Coleochaete [115] is evidence of a capacity for phenylpropanoid metabolism which would have been the ancestral basis for the elaboration of lignified vascular and mechanical issues, and later the rise of flavonoid pigmentation, in the land flora. Long-distance symplastic transport was already well established in green algae, including charophytes, and is observed in bryophytes alongside the beginnings of vascular translocation structures [40]. Long-range nutrient redistribution, from green cells in which the balance of turnover favours net dismantling of macromolecules and cell structures, is thus anatomically feasible in the charophytes.

The development of multicellular two- and three-dimensional forms changes the relationship between the plant body and the external medium. Direct exchange of materials across the plasmalemma has to be subordinated to the structural and physiological requirements of the tissue. It seems likely that this was the impetus for the evolution of more complex roles for the cell vacuole, which took on the status of 'inner space', in a sense adopting the role of a remnant of the ancestral aquatic milieu inside each vegetative cell. Extensive studies of the vacuoles of Chara and Nitella show them to perform sophisticated functions in ionic regulation, hydraulic maintenance of tissue morphology and turgor-driven cell wall growth [116–118]. They are also highly lytic, sequestering a range of proteases and other acid hydrolases [119–121]. The vacuole of charophytes thus possesses, or is en route to developing, most of the characteristics of the subcellular compartment that functions in angiosperm growth and senescence [122].

Plants invade the land

Green plants made landfall in the Silurian-Jurassic period, more than 410 mya. In this new and hostile environment, desiccation was avoided and ionic homeostasis maintained by the early development of an all-enveloping epidermal and cuticular 'space suit' [123] and a water-transport system [40, 124]. The capacity for lysigeny and schizogeny was established in the first land plants, in connection with both the differentiation of vascular systems and the shedding of parts [124, 125]. Along with mechanisms for smoothing out extreme variations in water-supply, it was necessary to establish physiological and developmental adaptations to deal with wide fluctuations in the inputs to photosynthesis. Many of the adjustments to the primary events of light capture and CO2 assimilation within chloroplasts are regulated by systems present in Chlamydomonas and other unicellular chlorophyceans [45]. Atmospheric CO2 concentration 410 mya is believed to have been around 3000 to 4000 ppmv [126]. This, combined with exposure to light fluxes no longer attenuated by the water-column, has led to the suggestion that the move to the terrestrial environment stimulated carbon fixation and utilisation to the point of incontinence – in the words of Harper [127] 'the green plant may indeed be a pathological overproducer of carbohydrates'. It allowed plants to develop profligate, throw-away lifestyles based on lysigeny and schizogeny on a vast scale. For example, primitive trees such as Archaeopteris (359-349 mya) probably controlled canopy morphology through abscission of lower megaphylls [128]. Abscission may have also been important in dehiscence of sporangia even before it became significant in vegetative shaping of plants [128], and a deciduous Glossopteris flora was already well established in Gondwanaland during the Carboniferous era [129]. Preserved senescent plastids have distinctive features that enable them to be identified in fossilised material [130]. Thomas and Sadras [131] discussed senescence as a strategy evolved by early land plants to deal with the promiscuous productivity of a photosynthetic machinery optimised for the relatively stable ancestral aquatic environment.

By leaving the water, plants lost the option of expelling unwanted or aggressive metabolites directly from individual cells into the outside environment. To a great degree the cell vacuole has taken over as the destination for such compounds. Chlorophyll degradation during the transdifferentiation of chloroplasts to gerontoplasts (Figure 12) illustrates this change in the relationship between cells and their environment. The catabolic pathway of chlorophyll in land plants is organised so as to limit photodynamic damage by free pigments. Chlorophylls and their coloured derivatives must be moved around and interconverted within strongly quenching microenvironments inside the plastid until the risk of photodamage is finally removed by opening the tetrapyrrole ring [9]. Sgr, PaO and Fibrillin are important participants in this phase of the process. Algae like Auxenochlorella expel RCC, the product of the PaO reaction, directly into the medium, possibly via an MRP2/4-like transporter as discussed above. Land plants, on the other hand, have requisitioned a detoxification pathway to redirect the products of chlorophyll catabolism to the cell vacuole. The two activities that make the link between the evolutionarily ancient first section of the catabolic pathway and the cytosolic-vacuolar steps are RCCR and the envelope transporter WBC23. The presence of the latter in Chlamydomonas and Ostreococcus is identified by BLAST; BLAST analysis also records hits on similar choanoflagellate and cnidarian sequences (Figure 8), suggesting that an ABC transporter similar to that which moves straight-chain bilins across the plastid envelope was present at or before the evolution of the chloroplastic endosymbiont.

Using site-specific mutagenesis methods, Pružinská et al. [132] have shown unequivocally the in vivo role of RCCR in chlorophyll catabolism; but the enzyme appears to have some other function too, concerned with mitochondrial participation in programmed cell death [133]. The occurrence of RCCR-like sequences is consistent with the appearance of this enzyme at the point of land plant evolution [17], though there are BLAST hits in the cyanobacteria (Figure 3). Why is there a discontinuity in the distribution of RCCR? Although there is a considerable step up in E value between the land plants and cyanobacteria, the Nostoc and Anabaena proteins have all the structural features of RCCR when examined by the conserved domain tools of NCBI BLAST and Pfam. BLASTing against CYANOBASE (the Cyanobacteria genome database) revealed weak similarity with phycoerythrobilin ferredoxin oxidoreductases, which may conceivably have been the distant ancestral sources of plant RCCRs. BLASTing the Nostoc punctiforme genome database reveals two genes located side-by-side (Npun_R3279 and Npun_R3277). They are highly similar and may represent a duplication. Pfam recognises both proteins as possessing the conserved domains of RCCR. Global transcription data are available for Nostoc and show that 3277 is strongly expressed during hormogonia differentiation but not in other developmental stages. By contrast 3279 was down-regulated during hormogonia differentiation [134]. Nostoc is an extracellular symbiont or endophyte in the hornwort Anthoceros [135], but an intracellular symbiont in the angiosperm Gunnera [136]. The hormogonium phase is the one that colonises the host. The available Anabaena genome (A variabilis) has only one RCCR-homologous sequence. Anabaena is symbiotic with the water fern Azolla [137].

It is significant that the only cyanobacteria known to have RCCR homologues (functional at the transcript level in the case of Nostoc at least) are those that form symbioses with green plants in tissues which are exposed to light. In these associations (especially the intracellular symbiosis), the normal option of excreting RCC into an aqueous medium would not be available. We speculate that this is the reason for the anomalous occurrence of RCCR genes in symbiotic Nostoc and Anabaena. The source of these genes remains problematical. Did land plants acquire RCCR by gene transfer between a cyanobacterium and a hornwort-like early colonizer host; or (more likely, since there is no obvious trace of the gene in other cyanobacteria) did a cyanobacterial symbiont pick up the gene from its host? The possibility of gene transfer between hosts and symbionts as a factor in the adaptation of early land plants to the new environment merits further study.

RCCR converts the phototoxic chlorophyll catabolite RCC into pFCC, which the envelope transporter WBC23 exports to the cytosol. The subsequent reaction sequence – conjugation followed by sequestration in the vacuole – is the usual means by which intracellular waste products or xenobiotics are detoxified [138, 139]. The tonoplast transporter MRP2 was discussed above.

Beerling [140] has proposed that the high atmospheric CO2 concentrations at the time when plants invaded the land persisted for 40–50 my, until the late Devonian, when a steep fall in CO2 drove the evolutionary transition from microphyllous plants (psilopsids, lycopsids, bryophytes) to the macrophyllous and euphyllous ferns, conifers and angiosperms. Beerling points out that the genes for leaf morphogenesis were to a large extent already present in clubmosses, spikemosses and quillworts and proposes that high CO2 constrained the expression of this morphogenetic potential. The elaboration of megaphyllous laterals from the Carboniferous onward led directly to the evolution of the specialised reproductive and propagative structures of the angiosperms according to long-established ontogenetic principles [39, 141]. Although light, carbon and (with qualifications) water may have been superabundant, mineral nutrients are likely to have been a limiting factor in early terrestrial ecosystems. Nitrogen in particular is released in only very small amounts from the weathering of rocks, though it is possible that biological N2 fixation by (cyano)bacteria and archea and nitrogen oxide fertilization by electrical storms could have contributed to the formation of skeletal oligotrophic soils as far back as 1200 mya [142]. The requirement for parsimony in the internal nitrogen economy of early land plants would have favoured the development of a number of structural, developmental and physiological traits contributing to the adoption of an integrated senescence program [38, 142].

Efficient internal nitrogen recycling requires cell integrity and metabolic regulation in the senescing source tissue, a directional translocation system connecting the source to the developing sink, means of signalling supply and demand status between the different regions of the plant body and a morphogenetic context that maintains structural and functional fitness [40, 143, 144]. Development of the thalli of liverworts clearly includes the capacity for senescence-like terminal changes (see, for example, Koeberl and Maravolo; [145]). Eckstein and Karlsson [146] showed the recycling of nitrogen from senescing to young ramets in mosses lacking a specialised vascular system. This suggests that source-sink integration of senescence into development was well established before the evolution of tracheophytes [40]. A search of the JGI Physcomitrella patens online database at the time of writing yields more than 330 hits with high sequence similarity to senescence-related genes in Arabidopsis, rice and Medicago, and proteomics studies, such as the analysis of responses of the P patens phosphoproteome to cytokinin treatment described by Heintz et al. [147], also identify gene products that may have a function in the developmental integration of senescence in this species. Highly similar homologues of all the representative senescence-related proteins analysed in the present study have been identified within the bryophytes by sequence-based analysis (Figure 21). We conclude that the genetic potential for expressing and regulating the senescence syndrome and deploying it as an adaptation to oligotrophic conditions was well established in the earliest land plants.

The senescence program is integral to strategies for dealing with the biotic and abiotic stresses of the terrestrial environment. Wrky53, member of a family of regulatory genes that diversified greatly with the evolution of the tracheophytes ([74]; Figure 9), encodes a transcription factor that functions both in senescence pathways and in a network of stress responses [148]. Wrky53 interacts with more than 60 genes of various kinds, including those for other transcription factors and for the senescence-related cysteine endopeptidase SAG12 [149]. Wrky53 is induced by H2O2, and autoregulates its own synthesis by feedback inhibition. Salicylic acid and jasmonic acid, via the regulatory protein ESR, control Wrky53 functions in senescence and resistance responses to pathogens [150]. Acclimation and adaptation to environmental stresses is often mediated by redox sensing. The stress-activated transcriptional cofactor NPR1 exists in the cytoplasm in the form of disulfide-bonded oligomers. A change in redox conditions results in reduction of the intermolecular disulfide bonds and relocation of monomeric NPR1 to the nucleus where it activates genes [151]. Miao et al. [149] reported that At4g26120, a gene encoding an NPR1-like protein, is a target for Wrky53 binding. The corresponding protein sequence shows better than 99% identity with that encoded by At1g64280, transcription of which is enhanced in senescence (Figure 20). Other genes listed in Figure 1 and profiled in Figure 20 with which Wrky53 interacts include At4g33030 (sulpholipid synthase), At5g45890 (SAG12 protease) and At5g14930 (SAG 101 lipase). It may be significant that the latter gene clusters immediately adjacent to Wrky53 in the Genevestigator diagram (Figure 20).

Tracheophytes, angiosperms and co-evolution

Angiosperms appeared in the Jurassic period, 150 mya. Carpels, flowers and tectate pollen occur in the fossil record at around this time and within 30 my angiosperms had become the dominant component of global floras. Before evolution of the angiosperms the colour world of vegetation would have consisted largely of greens, yellows and browns, much as it does in modern conifer-dominated forest ecosystems. A rapid increase in the profusion of natural pigmentation accompanied the Cretaceous explosion in angiosperm evolution and was likely to have been an adaptive response to a combination of abiotic and biotic factors [152, 153]. Radiation of the first angiosperms into a wide range of new and often hostile environments was made possible by physiological and morphogenetic innovations or consolidations including, for example, the development of flexible photosynthetic mechanisms and the proliferation of architectures and life-forms [152, 154]. Climate change may also have been a driving force at the onset of angiosperm diversification, since a transition between icehouse and greenhouse conditions occurred at this time [155].

Photoinhibition and the damaging effect of excess light energy within assimilatory tissues is a common factor in the experience of abiotic stresses [95, 156]. The ancient role of carotenoids in the dissipation of excess light energy within the photosynthetic apparatus [45] assumed new significance in land plants as they moved out into ecological niches in which the abiotic environment was increasingly subject to wide fluctuations that differentially affected the capture and utilisation of light energy [157, 158]. Amongst the carotenoid-related stress responses shared with senescence is the photoprotection role of Fibrillin [159]. Upregulation of Fibrillin expression in senescence (Figure 2) is consistent with the requirement for the maintenance of photoprotection during chlorophyll catabolism and protein mobilisation within the transdifferentiating plastid [9].

Coevolution with pollinators and dispersers is an important factor in establishing carotenoid- and phenylpropanoid-based pigmentation of floral organs and propagating structures such as seeds and fruits [153, 160–163]. The interaction between the colours of plant parts and the animal visual system is an example of extreme convergence in evolution, in the sense that the isoprenoid-protein light receptor structures and metabolic systems in the earliest photoautotrophs (evolved to exploit the main central part of the visible spectrum) show remarkable gene and protein sequence conservation through to both the carotenoid biochemistry of angiosperm organs and the eyes of the animals that interact with them. For example the highly conserved carotenoid-metabolising enzyme CCD8 (Figure 14) is a member of the same superfamily (RPE65) as the membrane proteins that bind retinal pigments in the eye. In contrast with the extremely ancient biochemistry and functions of isoprenoids, molecular phylogeny identifies flavonoids, including the hydrophilic vacuolar pigments of leaves, flowers and fruits, as innovations of the land plants (Figure 16, Figure 21; [164]), and their diversification as a factor in angiosperm radiation [153, 162, 165]. The evolutionary origin and significance of red anthocyanin coloration in senescing leaves are subjects of much current debate [8, 59, 166]. Modelling, and some experimental evidence, supports the idea of coevolution with herbivores, in which foliar anthocyanins act as signals to insects and other predators [167]. Alternatively, or additionally, anthocyanin accumulation may be a plant response to, or insurance against the effects of, abiotic stresses such as extremes of, or imbalances in, light, temperature, water availability or nutrient status [168]. We shall return to the question of the evolutionary relationship between the pigments of senescing foliage and the colours of flowers and fruits.

According to Ruddiman [169] the present era is the Anthropocene (a term coined by Crutzen and Stoermer; [170]) and dates from the origins of agriculture in human evolution, more than 8000 years ago. Human activity has intruded on plant evolution during the Anthropocene not only by subjecting it to environmental and climatic perturbations but also, and more directly, through selection, crop improvement and biotechnology. The domestication syndrome comprises, amongst a range of traits, annual growth habit, hypertrophied grains and foliage, rapid seedling establishment and high harvest-index [171]: all of which are influenced by, and influence, the timing and execution of senescence programs [172]. The senescence-associated NAC transcription factors ([41]; Figure 10) illustrate the point. Ooka et al. [173] analysed the NAC family of Arabidopsis and rice and placed At1g69490 and three other closely related sequences (At1g61110, At3g15510, At1g52880) in the NAP subgroup of Group I. Genevestigator expression profiling shows that each of the four AtNAPs in this family has a distinctly different tissue specificity: highest abundance of At1g69490 transcripts is in senescing leaves, with petals and sepals also giving a strong signal (Figure 2). At1g52880 expression is most intense in seed and fruit tissues, with some signal in senescing leaves; At3g15510 transcripts are most abundant in cork and xylem; whereas At1g61110 expression is confined to flowers and stamens (data not shown). At1g69490 is the top hit when the Hordeum senescence-associated NAM protein sequence is BLASTed against the TAIR8 protein dataset. Sequence alignment also places rice ONAC010 in the NAP subgroup, and BLASTing the rice sequence hits both the Hordeum and Triticum senescence-related NAM proteins with E values of better than 3e-100. A transcriptomic study by Hirose et al. [174] showed that ONAC010 expression in rice leaves is suppressed by treatment with the anti-senescence hormone zeatin. Divergence of function and structure is evidence of the plasticity of genes of the NAP/NAM family under natural and, potentially, human-mediated selection. Functional alleles of the NAM-B1 gene, encoding a transcription factor that accelerates senescence and increases nutrient remobilization, occur in ancestral wild wheat. Modern wheat varieties have been selected for delayed senescence, a trait correlated with enhanced grain yield [172]. They carry a non-functional allele of NAM-B1 [41]. Many further examples of human intervention to alter the senescence and ripening patterns of domesticated plants (with consequences for pigmentation or productivity or both) are discussed in [175–177].

Relationship of ontogeny and phylogeny in the evolution of senescence

During the present discussion we have moved freely across the full evolutionary range of taxonomic groupings and between cells, tissues and organs and their genes, regulators, enzymes, transporters, pigments and metabolites. To justify treating processes and structures, as well as genes and proteins, as homologues, we propose here a simple unifying developmental scheme. Following the principles of morphological and allometric transformation first described in 1917 by d'Arcy Thompson (revised edition 1992; [178]), a linear branch- or microphyll-like structure may evolve by 'rubber sheet' distortion into a two-dimensional lamina (a leaf or floral part for example), which in turn can be expanded in a third dimension to produce a fleshy spheroidal structure – a fruit, say (Figure 22). By adding a developmental time dimension to this figure, the progress of senescence may be represented by a sequence of pigmentation changes, green through yellow to red and finally the post-senescence transition to cell death accompanied by non-physiological darkening or bleaching. The result (Figure 22) is a spatial-temporal grid across which the aerial structures of multicellular green plants may travel freely during development or evolution. The grid itself is plastic and deformable so that particular structural or temporal stages are curtailed, extended or even subject to reversion, thereby giving rise to the profusion of forms and senescence behaviours seen in the land flora.

Figure 22
figure 22

Ontogenetic and phylogenetic interrelationships of plant organ senescence. The diagram represents morphogenetic relationships between one-, two- and three-dimensional organ structures (for example, linear axes, laminae and fleshy fruits respectively) and the time-related progression through senescence/ripening to death signified by characteristic changes in pigmentation.

This schematic description of the homologous interrelationships of organ senescence has important implications when considered together with the molecular relationships described in the present study. It is clear that the competence to transdifferentiate chloroplasts into gerontoplasts was already well established in the bryophytes and ferns, whereas chromoplasts appeared later when coevolution with the visual systems of interacting animals began to drive diversification of the colours and morphologies of reproductive structures (Additional File 3, Figure 22). In this regard it is significant that genes for Fibrillin, OrI and OrII, which originated phylogenetically long before chromoplasts had evolved (Figures 13, 15, are close to each other when clustered by expression profile (Figure 20, tracks 32, 27, 28). We argue, therefore, that there is a direct relationship between the phylogenetic precedence of gerontoplasts over chromoplasts on the one hand and their corresponding ontogenies on the other [113]. This implies that leaves which lose their chlorophyll and reveal and/or accumulate non-green pigments (carotenoids in all groups, anthocyanins arising with the cryptogams and angiosperms – Figure 21) are the evolutionary progenitors of the highly coloured flowers, fruits and seeds of dicots and monocots. It means that, phylogenetically and ontogenetically, the pigmented sepals, petals and dispersal structures of angiosperms are essentially heterochronic senescing leaves [110].

Moreover, if the developmental sequence of the two-dimensional structure in Figure 22 represents the ontogeny of a leaf and a homologous perianth structure such as a petal, then it is formation and pigmentation of the floral organ rather than the subsequent fading and death of the structure that expresses the modified senescence program. This view is not widely accepted (see [3, 179]), but receives some support from the literature (for example [180–182]) as well as from genomics. When senescence-related transcription profiles are used to cluster tissues and organs according to similarity of expressed gene complement (Figure 20), senescent leaf, petal, sepal and other floral parts form a clear grouping. On the other hand xylem and cork, tissues in which developmental cell death is an essential part of the differentiation program, are relatively distant from the leaf and floral clades (Figure 20), which is consistent with the argument that foliar senescence and homologous processes such as perianth development should be classed as transdifferentiation phenomena distinct from terminal cell death events [110].

The clustering patterns of Figure 20 provide further insights into the possible nature and origins of plant senescence and its regulation. It is striking that 11 of the 14 genes for isoprenoid metabolism analysed form a group between tracks 21 and 34. A coordinated unit of genes determining carotenoid functions in transdifferentiating plastids is suggested. Similarly, of 9 phenylpropanoid genes, 8 cluster between columns 49 and 61, indicating an integrated regulatory grouping. PaO, chlorophyll b reductase and Sgr are immediately adjacent to each other (tracks 37–39) and near to WBC23 (track 43). The remoteness of chlorophyllase from this cluster (column 13) may reflect growing opinion that the enzyme product of the corresponding gene (At5g43860) may not function in chlorophyll degradation in vivo [183]. The status of phaeophorbidase (column 20) in physiological chlorophyll catabolism is also uncertain. Like RCCR, which is nearby (track 23), no algal homologue has been detected [33], suggesting it too may be a relatively recent recruit to the pathway. The isolation of MRP2 (track 1) from other components of the chlorophyll breakdown pathway may reflect its variety of functions and ancient evolutionary origins. The associations revealed by cluster analysis (Figure 20) and the phylogenetic diversity of senescence-related genes (Figure 21) are consistent with the proposition that senescence is more a timetabled network than a single executable program [5]. Senescence has its evolutionary origin in a core of metabolic processes concerned with building and rebuilding photosynthetic complexes. Other biochemical, cellular, integrative and adaptive systems became accreted on this armature as the evolving plant encountered new environmental and developmental contexts. Current systems biology approaches [184] are leading towards the idea of senescence as a complex of reticulated gene interactions with nodes and hubs and away from the older conception of a master switch from which subsequent events propagated [5]. The evolutionary perspective is entirely consistent with such a system network model of plant senescence and its control.

Conclusion

In recent times the quantity and quality of DNA sequence information from plants across the taxonomic and evolutionary range has increased greatly. At the same time our understanding of the molecular and cell biology of plant senescence has developed to the point where we can assign functions to some enzymes and regulatory proteins. In the present paper we have taken the first steps towards a synthesis that allows us to frame, and begin to test, hypotheses about the evolutionary origins of angiosperm senescence, its control and integration into plant development and adaptation. Because major gaps remain in the taxonomic spread of sequenced species, and the assignment of functions to senescence-associated genes, the arguments presented here include a number of conjectural interpolations and extrapolations. For example, it is something of an article of faith for molecular evolution studies to consider a modern member of an ancient phylum (the lycopod Selaginella, for example) as representative of the genomic status of its newly-evolved ancestral form. Many more sequences from many more species, particularly covering the span between the cyanobacteria and the angiosperms, are needed before this belief can be held with full confidence. For the specific case of plant senescence, the scarcity of sequences from ferns and charophytes is a particular limitation. Even so, the present study has allowed a feasible scheme to be proposed that accounts for much of the available information on senescence-related genes, their functions and interactions.

Methods

Databases accessed

The following database sources were used in this study.

TAIR – the Arabidopsis Information Resource

http://www.arabidopsis.org/index.jsp

TIGR Arabidopsis database

http://www.tigr.org/tdb/e2k1/ath1/

SwissProt and other protein databases accessed through NCBI Entrez

Genbank and other nucleotide databases accessed through NCBI Entrez

Conserved domains accessed through NCBI Entrez

http://www.ncbi.nlm.nih.gov/Database/

UniProtKB protein knowledgebase

http://www.uniprot.org/

JGI Integrated Microbial Genomes (IMG) system

http://img.jgi.doe.gov/cgi-bin/pub/main.cgi

JGI Chlamydomonas reinhardtii v 3.0 BLAST

http://genome.jgi-psf.org/cgi-bin/runAlignment?db=Chlre3&advanced=1

CYANOBASE

http://bacteria.kazusa.or.jp/cyanobase/

JGI Physcomitrella patens online database

http://genome.jgi-psf.org/Phypa1_1/Phypa1_1.home.html

Gramene data resource for comparative genome analysis in the grasses

http://www.gramene.org/

Tree of Life web project

http://www.tolweb.org/tree/

Sequence alignment and structural comparison tools used

BLAST, including protein BLAST (blastp), nucleotide BLAST (blastn, megablast), blastx, tblastn, tree viewer and BLINK precomputed BLAST

http://blast.ncbi.nlm.nih.gov/Blast.cgi

EBI FASTA fast protein comparison

http://www.ebi.ac.uk/Tools/fasta/index.html

EBI ClustalW2 sequence alignment tool

http://www.ebi.ac.uk/Tools/clustalw2/index.html

GENEDOC sequence alignment editor

http://www.nrbsc.org/gfx/genedoc/index.html

EMBL SMART domain architecture tool

http://smart.embl-heidelberg.de/smart/set_mode.cgi?NORMAL=1

Pfam protein domain database

http://pfam.sanger.ac.uk/

Creating and displaying protein sequence trees

All protein trees reported here were generated using, as appropriate, BLASTP and FASTA for local alignments, ClustalW2 for global alignment, and TREE-PUZZLE, PUZZLEBOOT and the PHYLIP Phylogeny Inference Package [185, 186] for phylogeny and bootstrap analysis.

Protein sequence trees were constructed from amino acid sequence data using PROML (Protein Maximum Likelihood Program) in PHYLIP 3.62 for inferring evolutionary trees. PHYLIP (current version 3.68, released in August 2008) is distributed by the author, J Felsenstein, Department of Genome Sciences, University of Washington, Seattle. http://evolution.genetics.washington.edu/phylip.html. Input files for PROML, consisting of multiple sequence alignments, were created using ClustalW2. The Jones-Taylor-Thornton model was used to compute probabilities of changes between amino acid sequences.

For all amino acid sequence data sets, PROML was run with three different settings for the R (rate variation among sites) parameter: Gamma distribution [187], Hidden Markov Model (HMM [188]) and constant rate of change. The Log Likelihood (Ln) values of these three approaches for each set of protein sequences studied are summarised in Additional File 4. Although there were some minor variations in Ln values using the different approaches for estimating the among-site rate heterogeneity, the final tree conformations were the same regardless of which approach was used. Thus the among-site rate variation did not have significant influence on the topology of trees produced in this study.

Bootstrapping analyses were carried out on all datasets to measure how consistently the data support given taxon bipartitions. The bootstrapping majority-rule consensus tree and support values were generated through the use of SEQBOOT, NEIGHBOR and CONSENSE programs in PHYLIP, and PUZZLEBOOT http://rogerlab.biochemistryandmolecularbiology.dal.ca/Software/Software.htm#puzzleboot. High bootstrapping values mean uniform support; a support value of less than 50 is treated as giving no support. Support values generated in this study are presented in Additional File 1 and indicated in each tree figure. Tree branch lengths, generated using PROML with Gamma distribution, are summarised in Additional File 2. Bayesian analysis using the TOPALI http://www.topali.org/ package was also performed on all datasets and the output of support values from this analysis is consistent with the support values presented in this study.

To estimate the degree of rate heterogeneity, TREE-PUZZLE 5.2 [189, 190]http://www.tree-puzzle.de/ was run, using alignments created in ClustalW2, to generate the parameters such as the alpha (Gamma rate heterogeneity parameter) values, number of rate categories, individual rates and their respective probabilities to specify, as appropriate, coefficient of variation for a Gamma distribution model and how many categories of substitution rate there would be in a HMM of rate variation.

Protein trees were displayed using the TreeView visualisation tool:

http://taxonomy.zoology.gla.ac.uk/rod/treeview.html

Genevestigator

We profiled the expression of senescence-related genes in Arabidopsis using the Genevestigator reference expression database and meta-analysis system for Affymetrix® expression arrays [191]. Figure 2 presents output in heat-map format from the Meta-profile analysis tool based on screening of the ATH1 22 k array source for strength of expression categorised by Anatomy. Use of the Clustering analysis tool is shown in Figure 20. Both Anatomy and Gene categories were subjected to Hierarchical Clustering. Genevestigator displays absolute values either as MAS5-normalized signal values, or as MAS5 normalized values relative to the expression potential of each gene. The MAS5 normalization algorithm implemented by Genevestigator is described by Gentleman et al. [192]. The expression potential of a given gene is defined as the average of the top 1% signal values of a probe set across all arrays in the database. Its value is represented by the darkest blue colour, against which other levels of expression are normalised.

https://www.genevestigator.ethz.ch/gv/index.jsp

References

  1. Keskitalo J, Bergquist G, Gardeström P, Jansson S: A cellular timetable of autumn senescence. Plant Physiology. 2005, 139: 1635-1648.

    PubMed Central  CAS  Google Scholar 

  2. Petit RJ, Hampe A: Some evolutionary consequences of being a tree. Annual Review of Ecology, Evolution, and Systematics. 2006, 37: 187-214.

    Google Scholar 

  3. Rogers HJ: Programmed cell death in floral organs: how and why do flowers die?. Annals of Botany. 2006, 97: 309-315.

    PubMed Central  CAS  Google Scholar 

  4. Hopkins M, Taylor C, Liu Z, Ma F, McNamara L, Wang T-W, Thompson JE: Regulation and execution of molecular disassembly and catabolism during senescence. New Phytologist. 2007, 175: 201-214.

    CAS  Google Scholar 

  5. Jansson S, Thomas H: Senescence: developmental program or timetable?. New Phytologist. 2008, 179: 575-579.

    Google Scholar 

  6. Guo Y, Gan S: Leaf senescence: signals, execution, and regulation. Current Topics in Developmental Biology. 2005, 71: 83-112.

    CAS  Google Scholar 

  7. Lim PO, Kim HJ, Nam HG: Leaf senescence. Annual Review of Plant Biology. 2006, 58: 115-136.

    Google Scholar 

  8. Ougham HJ, Morris P, Thomas H: The colors of autumn leaves as symptoms of cellular recycling and defenses against environmental stresses. Current Topics in Developmental Biology. 2005, 66: 135-160.

    CAS  Google Scholar 

  9. Ougham H, Hörtensteiner S, Armstead I, Donnison I, King I, Thomas H, Mur L: The control of chlorophyll catabolism and the status of yellowing as a biomarker of leaf senescence. Plant Biology. 2008, 10 (supplement 1): 4-14.

    CAS  Google Scholar 

  10. Goldsworthy A: Why did Nature select green plants?. Nature. 1987, 328: 207-208.

    Google Scholar 

  11. Goldsworthy A: Why trees are green. New Scientist. 1987, 116 (1590): 48-52.

    Google Scholar 

  12. Dyson F: Our Biotech Future. New York Review of Books. 2007, 54 (12): [http://www.nybooks.com/articles/20370]

    Google Scholar 

  13. Bryant DA, Frigaard N-U: Prokaryotic photosynthesis and phototrophy illuminated. Trends in Microbiology. 2006, 14: 488-496.

    CAS  Google Scholar 

  14. DasSarma S: Extreme halophiles are models for astrobiology. Microbe. 2006, 1: 120-126.

    Google Scholar 

  15. Xiong J: Photosynthesis: what color was its origin?. Genome Biology. 2006, 7: 245-

    PubMed Central  Google Scholar 

  16. Xiong J, Bauer CE: Complex evolution of photosynthesis. Annual Review of Plant Biology. 2002, 53: 503-21.

    CAS  Google Scholar 

  17. Hörtensteiner S, Rodoni S, Schellenberg M, Vicentini F, Nandi OI, Qiu Y-L, Matile P: Evolution of chlorophyll degradation: the significance of RCC reductase. Plant Biology. 2000, 2: 63-67.

    Google Scholar 

  18. Rodoni S, Mühlecker W, Anderl M, Kräutler B, Moser D, Thomas H, Matile P, Hörtensteiner S: Chlorophyll breakdown in senescent chloroplasts. Cleavage of pheophorbide a in two enzymic steps. Plant Physiology. 1997, 115: 669-676.

    PubMed Central  CAS  Google Scholar 

  19. Armstead I, Donnison I, Aubry S, Harper J, Hörtensteiner S, James C, Mani J, Moffet M, Ougham H, Roberts L, Thomas A, Weeden N, Thomas H, King I: From crop to model to crop: identifying the genetic basis of the staygreen mutation in the Lolium/Festuca forage and amenity grasses. New Phytologist. 2006, 172: 592-597.

    Google Scholar 

  20. Aubry S: Biochemical and molecular analysis of chlorophyll degradation in higher plants. PhD Thesis. 2008, University of Zürich, Institute of Plant Biology

    Google Scholar 

  21. Strother PK: A speculative review of factors controlling the evolution of phytoplankton during Paleozoic time. Revue de Micropaleontologie. 2008, 51: 9-21.

    Google Scholar 

  22. Bengtson S: Origins and early evolution of predation. Paleontological Society Papers. 2002, 8: 289-317.

    Google Scholar 

  23. Cavalier-Smith T: Chloroplast evolution: secondary symbiogenesis and multiple losses. Current Biology. 2002, 12: R62-R64.

    CAS  Google Scholar 

  24. Huang J, Gogarten JP: Did an ancient chlamydial endosymbiosis facilitate the establishment of primary plastids?. Genome Biology. 2007, 8: R99-

    PubMed Central  Google Scholar 

  25. Huang J, Gogarten JP: Concerted gene recruitment in early plant evolution. Genome Biology. 2008, 9: R109-

    PubMed Central  Google Scholar 

  26. Raven JA, Andrews M, Quigg A: The evolution of oligotrophy: implications for the breeding of crop plants for low input agricultural systems. Annals of Applied Biology. 2005, 146: 261-280.

    CAS  Google Scholar 

  27. Dangl JL, Dietrich RA, Thomas H: Senescence and programmed cell death. Biochemistry and Molecular Biology of Plants. Edited by: Buchanan B, Gruissem W, Jones R. 2000, Rockville: ASPP, 1044-1100.

    Google Scholar 

  28. Hörtensteiner S, Chinner J, Matile P, Thomas H, Donnison IS: Chlorophyll breakdown in Chlorella protothecoides: characterization of degreening and cloning of degreening-related genes. Plant Molecular Biology. 2000, 42: 439-450.

    Google Scholar 

  29. Aoki SM, Matsuka M, Hase E: De- and re-generation of chloroplasts in the cells of Chlorella protothecoides. V. Degeneration of chloroplasts induced by different carbon sources, and effects of some antimetabolites upon the process induced by glucose. Plant and Cell Physiology. 1965, 6: 487-497.

    CAS  Google Scholar 

  30. Oshio Y, Hase E: Studies on red pigments excreted by cells of Chlorella protothecoides during the process of bleaching induced by glucose or acetate. I. Chemical properties of the red pigments. Plant and Cell Physiology. 1969, 10: 41-49.

    CAS  Google Scholar 

  31. Oshio Y, Hase E: Studies on red pigments excreted by cells of Chlorella protothecoides during the process of bleaching induced by glucose or acetate. II. Mode of formation of the red pigments. Plant and Cell Physiology. 1969, 10: 51-59.

    CAS  Google Scholar 

  32. Engel N, Jenny TA, Mooser V, Gossauer A: Chlorophyll catabolism in Chlorella protothecoides. Isolation and structural elucidation of a red bilin derivative. FEBS Letters. 1991, 293: 131-133.

    CAS  Google Scholar 

  33. Hörtensteiner S: Chlorophyll degradation during senescence. Annual Review of Plant Biology. 2006, 57: 55-77.

    Google Scholar 

  34. Doi M, Shima S, Egashira T, Nakamura K, Okayama S: New bile pigment excreted by a Chlamydomonas reinhardtii mutant: A possible breakdown catabolite of chlorophyll a. Journal of Plant Physiology. 1997, 150: 504-508.

    CAS  Google Scholar 

  35. Wua C, Akimotoa H, Ohmiya Y: Tracer studies on dinoflagellate luciferin with [15N]-glycine and [15N]-L-glutamic acid in the dinoflagellate Pyrocystis lunula. Tetrahedron Letters. 2003, 44: 1263-1266.

    Google Scholar 

  36. Nakamura H, Musicki B, Kishi Y: Structure of the light emitter in krill (Euphausia pacifica) bioluminescence. Journal of the American Chemical Society. 1988, 110: 2683-2685.

    CAS  Google Scholar 

  37. Zhang C, Barthelson RA, Lambert GM, Galbraith DW: Global characterization of cell-specific gene expression through fluorescence-activated sorting of nuclei. Plant Physiology. 2008, 147: 30-40.

    PubMed Central  CAS  Google Scholar 

  38. Bleecker AB: The evolutionary basis of leaf senescence: Method to the madness?. Current Opinion in Plant Biology. 1998, 1: 73-78.

    CAS  Google Scholar 

  39. Kaplan DR: The science of plant morphology: definition, history, and role in modern biology. American Journal of Botany. 2001, 88: 1711-1741.

    CAS  Google Scholar 

  40. Raven JA: Long-distance transport in non-vascular plants. Plant, Cell and Environment. 2003, 26: 73-85.

    Google Scholar 

  41. Uauy C, Distelfeld A, Fahima T, Blechl A, Dubcovsky J: A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science. 2006, 314: 1298-1301.

    CAS  Google Scholar 

  42. Guo Y, Gan S: AtNAP, a NAC family transcription factor, has an important role in leaf senescence. Plant Journal. 2006, 46: 601-612.

    CAS  Google Scholar 

  43. McCourt RM, Delwiche CF, Karol KG: Charophyte algae and land plant origins. Trends in Ecology and Evolution. 2004, 19: 661-666.

    Google Scholar 

  44. Hinderhofer K, Zentgraf U: Identification of a transcription factor specifically expressed at the onset of leaf senescence. Planta. 2001, 213: 469-73.

    CAS  Google Scholar 

  45. Niyogi KK: Photoprotection revisited: genetic and molecular approaches. Annual Review of Plant Physiology and Plant Molecular Biology. 1999, 50: 333-59.

    CAS  Google Scholar 

  46. Schaefer HM, Rolshausen G: Plants on red alert: do insects pay attention?. Bio Essays. 2005, 28: 65-71.

    Google Scholar 

  47. Grotewold E: The genetics and biochemistry of floral pigments. Annual Review of Plant Biology. 2006, 57: 761-780.

    CAS  Google Scholar 

  48. Boucher Y, Doolittle WF: The role of lateral gene transfer in the evolution of isoprenoid biosynthesis pathways. Molecular Microbiology. 2000, 37: 703-16.

    CAS  Google Scholar 

  49. Lange BM, Rujan T, Martin W, Croteau R: Isoprenoid biosynthesis: the evolution of two ancient and distinct pathways across genomes. Proceedings of the National Academy of Sciences, USA. 2000, 97: 13172-77.

    CAS  Google Scholar 

  50. Steinmüller D, Tevini M: Composition and function of plastoglobuli. I. Isolation and purification from chloroplasts and chromoplasts. Planta. 1985, 163: 201-207.

    Google Scholar 

  51. Shihira-Ishikawa I, Hase E: Nutritional control of cell pigmentation in Chlorella protothecoides with special reference to the degeneration of chloroplast induced by glucose. Plant and Cell Physiology. 1964, 5: 227-240.

    Google Scholar 

  52. Wältermann M, Steinbüchel A: Neutral lipid bodies in prokaryotes: recent insights into structure, formation, and relationship to eukaryotic lipid depots. Journal of Bacteriology. 2005, 187: 3607-3619.

    PubMed Central  Google Scholar 

  53. Tevini M, Steinmüller D: Composition and function of plastoglobuli. II. Lipid composition of leaves and plastoglobuli during beech leaf senescence. Planta. 1985, 163: 91-96.

    CAS  Google Scholar 

  54. Camara B, Hugueney P, Bouvier F, Kuntz M, Moneger R: Biochemistry and molecular biology of chromoplast development. International Review of Cytology. 1995, 163: 175-247.

    CAS  Google Scholar 

  55. Brehelin C, Kessler F, van Wijk KJ: Plastoglobules: versatile lipoprotein particles in plastids. Trends in Plant Science. 2007, 12: 260-266.

    CAS  Google Scholar 

  56. Giuliano G, Diretto G: Of chromoplasts and chaperones. Trends in Plant Science. 2007, 12: 529-531.

    CAS  Google Scholar 

  57. Snowden KC, Simkin AJ, Janssen BJ, Templeton KR, Loucas HM, Simons JL, Karunairetnam S, Gleave AP, Clark DG, Klee HJ: The Decreased Apical Dominance1/Petunia hybrida CAROTENOID CLEAVAGE DIOXYGENASE8 gene affects branch production and plays a role in leaf senescence, root growth, and flower development. Plant Cell. 2005, 17: 746-759.

    PubMed Central  CAS  Google Scholar 

  58. Tanaka Y, Sasaki N, Ohmiya A: Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant Journal. 2008, 54: 733-749.

    CAS  Google Scholar 

  59. Ougham H, Thomas H, Archetti M: The adaptive value of leaf colour. New Phytologist. 2008, 179: 9-13.

    Google Scholar 

  60. Alfenito MR, Souer E, Goodman CD, Buell R, Mol J, Koes R, Walbot V: Functional complementation of anthocyanin sequestration in the vacuole by widely divergent glutathione S-transferases. Plant Cell. 1998, 10: 1135-1149.

    PubMed Central  CAS  Google Scholar 

  61. Furtek D, Schiefelbein JW, Johnston F, Nelson OE: Sequence comparisons of 3 wild-type bronze-1 alleles from Zea mays. Plant Molecular Biology. 1988, 11: 473-481.

    CAS  Google Scholar 

  62. Paz-Ares J, Ghosal D, Saedler H: Molecular analysis of the C1-I allele from Zea mays: a dominant mutant of the regulatory C1 locus. EMBO Journal. 1990, 9: 315-321.

    PubMed Central  CAS  Google Scholar 

  63. Lu Y-P, Li Z-S, Drozdowicz Y-M, Hörtensteiner S, Martinoia E, Rea PA: AtMRP2, an Arabidopsis ATP binding cassette transporter able to transport glutathione S-conjugates and chlorophyll catabolites: functional comparisons with AtMRP1. Plant Cell. 1998, 10: 267-282.

    PubMed Central  CAS  Google Scholar 

  64. Carter C, Pan S, Zouhar J, Avila EL, Girke T, Raikhel NV: The vegetative vacuole proteome of Arabidopsis thaliana reveals predicted and unexpected proteins. Plant Cell. 2004, 16: 3285-3303.

    PubMed Central  CAS  Google Scholar 

  65. Donnison IS, Gay AP, Thomas H, Edwards KJ, Edwards D, James CL, Thomas AM, Ougham HJ: Modification of nitrogen remobilisation, grain fill and leaf senescence in maize (Zea mays L.) by transposon insertional mutagenesis in a protease gene. New Phytologist. 2007, 173: 481-494.

    CAS  Google Scholar 

  66. Koonin EV: Orthologs, paralogs, and evolutionary genomics. Annual Review of Genetics. 2005, 39: 309-38.

    CAS  Google Scholar 

  67. Degnan JH, Rosenberg NA: Discordance of species trees with their most likely gene trees. PLoS Genetics. 2006, 2 (5): e68-

    PubMed Central  Google Scholar 

  68. Wüthrich KL, Bovet L, Hunziker PE, Donnison IS, Hörtensteiner S: Molecular cloning, functional expression and characterisation of RCC reductase involved in chlorophyll catabolism. Plant Journal. 2000, 21: 189-198.

    Google Scholar 

  69. Park S-Y, Yu J-W, Park J-S, Li J, Yoo S-C, Lee N-Y, Lee S-K, Jeong S-W, Seo HS, Koh H-J, Jeon J-S, Park Y-I, Paek N-C: The senescence-induced staygreen protein regulates chlorophyll degradation. Plant Cell. 2007, 19: 1649-1664.

    PubMed Central  CAS  Google Scholar 

  70. Sato Y, Morita R, Nishimura M, Yamaguchi H, Kusaba M: Mendel's green cotyledon gene encodes a positive regulator of the chlorophyll-degrading pathway. Proceedings of the National Academy of Sciences, USA. 2007, 104: 14169-14174.

    CAS  Google Scholar 

  71. Collins MD, Lawson PA, Willems A, Cordoba JJ, Fernandez-Garayzabal J, Garcia P, Cai J, Hippe H, Farrow JAE: The phylogeny of the genus Clostridium: proposal of five new genera and eleven new species combinations. International Journal of Systematic Bacteriology. 1994, 44: 812-826.

    CAS  Google Scholar 

  72. Chen Y, Korkeala H, Aarnikunnas J, Lindstrom M: Sequencing the botulinum neurotoxin gene and related genes in Clostridium botulinum type E strains reveals orfx3 and a novel type E neurotoxin subtype. Journal of Bacteriology. 2007, 189: 8643-8650.

    PubMed Central  CAS  Google Scholar 

  73. Rea PA: Plant ATP-Binding Cassette transporters. Annual Review of Plant Biology. 2007, 58: 347-375.

    CAS  Google Scholar 

  74. Zhang Y, Wang L: The WRKY transcription factor superfamily: its origin in eukaryotes and expansion in plants. BMC Evolutionary Biology. 2005, 5: 1-

    PubMed Central  CAS  Google Scholar 

  75. Thomson WW, Whatley JM: Development of nongreen plastids. Annual Review of Plant Physiology. 1980, 31: 375-394.

    Google Scholar 

  76. Parthier B: Gerontoplasts: the yellow end in the ontogenesis of chloroplasts. Endocytobiosis and Cell Research. 1988, 5: 163-190.

    Google Scholar 

  77. Zavaleta-Mancera HA, Thomas BJ, Thomas H, Scott IM: Regreening of senescent Nicotiana leaves. II. Redifferentiation of plastids. Journal of Experimental Botany. 1999, 50: 1683-1689.

    CAS  Google Scholar 

  78. Staehelin LA, Newcomb EH: Membrane structure and membranous organelles. Biochemistry and Molecular Biology of Plants. Edited by: Buchanan B, Gruissem W, Jones R. 2000, Rockville: ASPP, 2-50.

    Google Scholar 

  79. Siddique MA, Grossmann J, Gruissem W, Baginsky S: Proteome analysis of bell pepper (Capsicum annuum L.) chromoplasts. Plant and Cell Physiology. 2006, 47: 1663-1673.

    CAS  Google Scholar 

  80. Lu S, Van Eck J, Zhou X, Lopez AB, O'Halloran DM, Cosman KM, Conlin BJ, Paolillo DJ, Garvin DF, Vrebalov J, Kochian LV, Kupper H, Earle ED, Cao J, Li L: The cauliflower Or gene encodes a DNAJ cysteine-rich domain-containing protein that mediates high levels of beta-carotene accumulation. Plant Cell. 2006, 18: 3594-3605.

    PubMed Central  CAS  Google Scholar 

  81. Diaz C, Saliba-Colombani V, Loudet O, Belluomo P, Moreau L, Daniel-Vedele F, Morot-Gaudry J-F, Masclaux-Daubresse C: Leaf yellowing and anthocyanin accumulation are two genetically independent strategies in response to nitrogen limitation in Arabidopsis thaliana. Plant and Cell Physiology. 2006, 47: 74-83.

    CAS  Google Scholar 

  82. Smart CM, Hosken SE, Thomas H, Greaves JA, Blair BG, Schuch W: The timing of maize leaf senescence and characterization of senescence-related cDNAs. Physiologia Plantarum. 1995, 93: 673-682.

    CAS  Google Scholar 

  83. Thomas H, Donnison I: Back from the brink: plant senescence and its reversibility. Programmed Cell Death in Animals and Plants. Edited by: Bryant J, Hughes SG, Garland JM. 2000, Oxford: Bios, 149-162.

    Google Scholar 

  84. Doolittle WF, Bapteste E: Pattern pluralism and the Tree of Life hypothesis. Proceedings of the National Academy of Sciences, USA. 2007, 104: 2043-2049.

    CAS  Google Scholar 

  85. Rakhimberdieva MG, Stadnichuk IN, Elanskaya IV, Karapetyan NV: Carotenoid-induced quenching of the phycobilisome fluorescence in photosystem II-deficient mutant of Synechocystis sp. FEBS Letters. 2004, 574: 85-88.

    CAS  Google Scholar 

  86. Jansson S: A protein family saga: from photoprotection to light-harvesting (and back?). Advances in Photosynthesis and Respiration. 2006, 21: 145-153.

    Google Scholar 

  87. Minai L, Wostrikoff K, Wollman F-A, Choquet Y: Chloroplast biogenesis of Photosystem II cores involves a series of assembly-controlled steps that regulate translation. Plant Cell. 2006, 18: 159-175.

    PubMed Central  CAS  Google Scholar 

  88. Göhre V, Ossenbühl F, Crèvecoeur M, Eichacker LA, Rochaix J-D: One of two Alb3 proteins is essential for the assembly of the photosystems and for cell survival in Chlamydomonas. Plant Cell. 2006, 18: 1454-1466.

    PubMed Central  Google Scholar 

  89. Thomas H: Chlorophyll: a symptom and a regulator of plastid development. New Phytologist. 1997, 136: 163-181.

    CAS  Google Scholar 

  90. Sakamoto W: Protein degradation machineries in plastids. Annual Review of Plant Biology. 2006, 57: 599-621.

    CAS  Google Scholar 

  91. Ma J, Peng L, Guo J, Lu Q, Lu C, Zhang L: LPA2 is required for efficient assembly of Photosystem II in Arabidopsis thaliana. Plant Cell. 2007, 19: 1980-1993.

    PubMed Central  CAS  Google Scholar 

  92. Fu A, He Z, Cho HS, Lima A, Buchanan BB, Luan S: A chloroplast cyclophilin functions in the assembly and maintenance of photosystem II in Arabidopsis thaliana. Proceedings of the National Academy of Sciences, USA. 2007, 104: 15947-15952.

    CAS  Google Scholar 

  93. Martinez M, Diaz I: The origin and evolution of plant cystatins and their target cysteine proteinases indicate a complex functional relationship. BMC Evolutionary Biology. 2008, 8: 198-

    PubMed Central  Google Scholar 

  94. Cavalier-Smith T: Rooting the tree of life by transition analyses. Biology Direct. 2006, 1: 19-

    PubMed Central  Google Scholar 

  95. Melis A: Photosystem-II damage and repair cycle in chloroplasts: what modulates the rate of photodamage in vivo?. Trends in Plant Science. 1999, 4: 130-135.

    Google Scholar 

  96. Ougham H, Armstead I, Howarth C, Galyuon I, Donnison I, Thomas H: The genetic control of senescence revealed by mapping quantitative trait loci. Annual Plant Reviews. 2007, 26: 171-201.

    CAS  Google Scholar 

  97. Padham AK, Hopkins MT, Wang T-W, McNamara LM, Lo M, Richardson LGL, Smith MD, Taylor CA, Thompson JE: Characterization of a plastid triacylglycerol lipase from Arabidopsis. Plant Physiology. 2007, 143: 1372-1384.

    PubMed Central  CAS  Google Scholar 

  98. Hedges SB, Blair JE, Venturi M, Shoe JL: A molecular timescale of eukaryote evolution and the rise of complex multicellular life. BMC Evolutionary Biology. 2004, 4: 2-

    PubMed Central  Google Scholar 

  99. Keech O, Pesquet E, Ahad A, Askne A, Nordvall D, Vodnala SM, Tuominen H, Hurry V, Dizengremel P, Gardeström P: The different fates of mitochondria and chloroplasts during dark-induced senescence in Arabidopsis leaves. Plant, Cell and Environment. 2007, 30: 1523-1534.

    CAS  Google Scholar 

  100. Plaxton WC, Podestá FE: The functional organization and control of plant respiration. Critical Reviews in Plant Sciences. 2006, 25: 159-198.

    CAS  Google Scholar 

  101. Mae T: Leaf senescence and nitrogen metabolism. Plant Cell Death Processes. Edited by: Noodén LD. 2004, NY: Academic Press, 157-168.

    Google Scholar 

  102. Feller U, Anders I, Mae T: Rubiscolytics: fate of Rubisco after its enzymatic function in a cell is terminated. Journal of Experimental Botany. 2007, 59: 1615-1624.

    Google Scholar 

  103. Guiamét JJ, Pichersky E, Noodén LD: Mass exodus from senescing soybean chloroplasts. Plant and Cell Physiology. 1999, 40: 986-992.

    Google Scholar 

  104. Otegui MS, Noh Y-S, Martínez DE, Vila Petroff MG, Staehelin LA, Amasino RM, Guiamét JG: Senescence-associated vacuoles with intense proteolytic activity develop in leaves of Arabidopsis and soybean. Plant Journal. 2005, 41: 831-844.

    CAS  Google Scholar 

  105. Bassham DC: Plant autophagy–more than a starvation response. Current Opinion in Plant Biology. 2007, 10: 587-593.

    CAS  Google Scholar 

  106. Grbic V: SAG2 and SAG12 protein expression in senescing Arabidopsis plants. Physiologia Plantarum. 2003, 119: 263-269.

    CAS  Google Scholar 

  107. Trobacher CP, Senatore A, Greenwood JS: Masterminds or minions? Cysteine proteinases in plant programmed cell death. Canadian Journal of Botany. 2006, 84: 651-667.

    CAS  Google Scholar 

  108. Derelle E, Ferraz C, Rombauts S, Rouze P, Worden AZ, Robbens S, Partensky F, Degroeve S, Echeynie S, Cooke R, Saeys Y, Wuyts J, Jabbari K, Bowler C, Panaud O, Piegu B, Ball S, Ral JP, Bouget FY, Piganeau G, De Baets B, Picard A, Delseny M, Demaille J, Peer Van de Y, Moreau H: Genome analysis of the smallest free-living eukaryote Ostreococcus tauri unveils many unique features. Proceedings of the National Academy of Sciences, USA. 2006, 103: 11647-11652.

    CAS  Google Scholar 

  109. Bidle KD, Bender SJ: Iron starvation and culture age activate metacaspases and programmed cell death in the marine diatom Thalassiosira pseudonana. Eukaryotic Cell. 2008, 7: 223-236.

    PubMed Central  CAS  Google Scholar 

  110. Thomas H, Ougham H, Wagstaff C, Stead AJ: Defining senescence and death. Journal of Experimental Botany. 2003, 54: 1127-1132.

    CAS  Google Scholar 

  111. Fukuda H: Xylogenesis: initiation, progression, and cell death. Annual Review of Plant Physiology and Plant Molecular Biology. 1996, 47: 299-325.

    CAS  Google Scholar 

  112. Samuels AL, Kaneda M, Rensing KH: The cell biology of wood formation: from cambial divisions to mature secondary xylem. Canadian Journal of Botany. 2006, 84: 631-639.

    CAS  Google Scholar 

  113. Matile P, Hörtensteiner S, Thomas H: Chlorophyll degradation. Annual Review of Plant Physiology and Plant Molecular Biology. 1999, 50: 67-95.

    CAS  Google Scholar 

  114. Delwiche CF, Karol KG, Cimino MT, Sytsma KJ: Phylogeny of the genus Coleochaete, Coleochaetales, Charophyta and related taxa inferred by analysis of the chloroplast gene rbcL. Journal of Phycology. 2002, 38: 394-403.

    CAS  Google Scholar 

  115. Delwiche CF, Graham LE, Thomson N: Lignin-like compounds and sporopollenin in Coleochaete, an algal model for land plant ancestry. Science. 1989, 245: 399-401.

    CAS  Google Scholar 

  116. Proseus TE, Boyer JS: Periplasm turgor pressure controls wall deposition and assembly in growing Chara corallina cells. Annals of Botany. 2006, 98: 93-105.

    PubMed Central  CAS  Google Scholar 

  117. Beilby MJ: Action potential in Charophytes. International Review of Cytology. 2007, 257: 43-83.

    CAS  Google Scholar 

  118. Limbach C, Staehelin LA, Sievers A, Braun M: Electron tomographic characterization of a vacuolar reticulum and of six vesicle types that occupy different cytoplasmic domains in the apex of tip-growing Chara rhizoids. Planta. 2008, 227: 1101-1114.

    CAS  Google Scholar 

  119. Doi E, Ohtsuru C, Matoba T: Lysosomal nature of plant vacuoles. II. Acid hydrolases in the central vacuole of internodal cells of Charophyta. Plant and Cell Physiology. 1975, 16: 581-588.

    CAS  Google Scholar 

  120. Moriyasu Y, Tazawa M: Distribution of several proteases inside and outside the central vacuole of Chara australis. Cell Structure and Function. 1986, 11: 81-86.

    CAS  Google Scholar 

  121. Moriyasu Y: Examination of the contribution of vacuolar proteases to intracellular protein degradation in Chara corallina. Plant Physiology. 1995, 109: 1309-1315.

    PubMed Central  CAS  Google Scholar 

  122. Taiz L: The plant vacuole. Journal of Experimental Biology. 1992, 172: 113-122.

    CAS  Google Scholar 

  123. Corner EJH: The Life of Plants. 1964, London: Weidenfeld and Nicolson

    Google Scholar 

  124. Edwards D, Li C-S, Raven JA: Tracheids in an early vascular plant: a tale of two branches. Botanical Journal of the Linnaean Society. 2006, 150: 115-130.

    Google Scholar 

  125. Edwards D: Cells and tissues in the vegetative sporophytes of early land plants. New Phytologist. 1993, 125: 225-247.

    Google Scholar 

  126. Beerling DJ, Berner RA: Feedbacks and the coevolution of plants and atmospheric CO2. Proceedings of the National Academy of Sciences, USA. 2005, 102: 1302-1305.

    CAS  Google Scholar 

  127. Harper JL: The Population Biology of Plants. 1977, London: Academic Press

    Google Scholar 

  128. Raven JA: Evolution of plant life forms. On the economy of plant form and function. Edited by: Givnish TJ. 1986, New York: Cambridge University Press, 421-492.

    Google Scholar 

  129. Pigg KB, Trivett ML: Evolution of the glossopterid gymnosperms from Permian Gondwana. Journal of Plant Research. 1994, 107: 461-477.

    Google Scholar 

  130. Niklas KJ: The chemistry of fossil plants. Bio Science. 1981, 31: 820-825.

    CAS  Google Scholar 

  131. Thomas H, Sadras VO: The capture and gratuitous disposal of resources by plants. Functional Ecology. 2001, 15: 3-12.

    Google Scholar 

  132. Pružinská A, Anders I, Aubry S, Schenk N, Tapernoux-Lüthi E, Müller T, Kräutler B, Hörtensteiner S: In vivo participation of red chlorophyll catabolite reductase in chlorophyll breakdown. Plant Cell. 2007, 19: 369-387.

    PubMed Central  Google Scholar 

  133. Yao N, Eisfelder BJ, Marvin J, Greenberg JT: The mitochondrion, an organelle commonly involved in programmed cell death in Arabidopsis thaliana. Plant Journal. 2004, 40: 596-610.

    CAS  Google Scholar 

  134. Campbell EL, Summers ML, Christman H, Martin ME, Meeks JC: Global gene expression patterns of Nostoc punctiforme in steady-state dinitrogen-grown heterocyst-containing cultures and at single time points during the differentiation of akinetes and hormogonia. Journal of Bacteriology. 2007, 189: 5247-5256.

    PubMed Central  CAS  Google Scholar 

  135. Peirce GJ: Anthoceros and its Nostoc colonies. Botanical Gazette. 1906, 42: 55-59.

    Google Scholar 

  136. Johansson C, Bergman B: Early events during the establishment of the Gunnera/Nostoc symbiosis. Planta. 1992, 188: 403-413.

    CAS  Google Scholar 

  137. Peters GA, Meeks JC: The Azolla-Anabaena symbiosis: basic biology. Annual Review of Plant Physiology and Plant Molecular Biology. 1989, 40: 193-210.

    Google Scholar 

  138. Kreuz K, Tommasini R, Martinoia E: Old enzymes for a new job. Herbicide detoxification in plants. Plant Physiology. 1996, 111: 349-353.

    PubMed Central  CAS  Google Scholar 

  139. Hinder B, Schellenberg M, Rodoni S, Ginsburg S, Vogt E, Martinoia E, Matile P, Hörtensteiner S: How plants dispose of chlorophyll catabolites. Directly energized uptake of tetrapyrrolic breakdown products into isolated vacuoles. Journal of Biological Chemistry. 1996, 271: 27233-27236.

    CAS  Google Scholar 

  140. Beerling DJ: Leaf evolution: gases, genes and geochemistry. Annals of Botany. 2005, 96: 345-352.

    PubMed Central  CAS  Google Scholar 

  141. Classen-Bockhoff R: Plant morphology: the historic concepts of Wilhelm Troll, Walter Zimmermann and Agnes Arber. Annals of Botany. 2001, 88: 1153-1172.

    Google Scholar 

  142. Raven JA, Edwards D: Roots: evolutionary origins and biogeochemical significance. Journal of Experimental Botany. 2001, 52: 382-401.

    Google Scholar 

  143. Feller U, Fischer A: Nitrogen metabolism in senescing leaves. Critical Reviews in Plant Science. 1994, 13: 241-273.

    CAS  Google Scholar 

  144. Hörtensteiner S: The loss of green color during chlorophyll degradation – a prerequisite to prevent cell death?. Planta. 2004, 219: 191-194.

    Google Scholar 

  145. Koeberl MJ, Maravolo NC: The influence of spermine and 6-benzylaminopurine on putative caspase activity during senescence in Marchantia polymorpha thalli. International Journal of Plant Sciences. 2006, 167: 1115-1122.

    CAS  Google Scholar 

  146. Eckstein RL, Karlsson PS: Recycling of nitrogen among segments of Hylocomium splendens as compared with Polytrichum commune: implications for clonal integration in an ectohydric bryophyte. Oikos. 1999, 86: 87-96.

    Google Scholar 

  147. Heintz D, Erxleben A, High AA, Wurtz V, Reski R, Van Dorsselaer A, Sarnighausen E: Rapid alteration of the phosphoproteome in the moss Physcomitrella patens after cytokinin treatment. Journal of Proteome Research. 2006, 5: 2283-93.

    CAS  Google Scholar 

  148. Zimmermann P, Zentgraf U: The correlation between oxidative stress and leaf senescence during plant development. Cellular and Molecular Biology Letters. 2005, 10: 515-534.

    CAS  Google Scholar 

  149. Miao Y, Laun T, Zimmermann P, Zentgraf U: Targets of the WRKY53 transcription factor and its role during leaf senescence in Arabidopsis. Plant Molecular Biology. 2004, 55: 853-867.

    CAS  Google Scholar 

  150. Miao Y, Zentgraf U: The antagonist function of Arabidopsis WRKY53 and ESR/ESP in leaf senescence is modulated by the jasmonic and salicylic acid equilibrium. Plant Cell. 2007, 19: 819-830.

    PubMed Central  CAS  Google Scholar 

  151. Federoff N: Redox regulatory mechanisms in cellular stress responses. Annals of Botany. 2006, 98: 289-300.

    Google Scholar 

  152. Feild TS, Arens NC: Form, function and environments of the early angiosperms: merging extant phylogeny and ecophysiology with fossils. New Phytologist. 2005, 166: 383-408.

    Google Scholar 

  153. Rausher MD: Evolutionary transitions in floral color. International Journal of Plant Science. 2008, 169: 7-21.

    CAS  Google Scholar 

  154. Ackerly DD, Monson RK: Waking the sleeping giant: the evolutionary foundations of plant function. International Journal of Plant Science. 2003, 164 (3 supplement): S1-S6.

    Google Scholar 

  155. Jahren AH, Arens NC, Sarmiento G, Guerrero J, Amundson R: Terrestrial record of methane hydrate dissociation in the Early Cretaceous. Geology. 2001, 29: 159-162.

    CAS  Google Scholar 

  156. Murata N, Takahashi S, Nishiyama Y, Allakhverdiev SI: Photoinhibition of photosystem II under environmental stress. Biochimica et Biophysica Acta – Bioenergetics. 2007, 1767: 414-421.

    CAS  Google Scholar 

  157. Demmig-Adams B, Adams WW: Antioxidants in photosynthesis and human nutrition. Science. 2002, 298: 2149-2153.

    CAS  Google Scholar 

  158. Davison PA, Hunter CN, Horton P: Overexpression of β-carotene hydroxylase enhances stress tolerance in Arabidopsis. Nature. 2002, 418: 203-206.

    CAS  Google Scholar 

  159. Yang Y, Sulpice R, Himmelbach A, Meinhard M, Christmann A, Grill E: Fibrillin expression is regulated by abscisic acid response regulators and is involved in abscisic acid-mediated photoprotection. Proceedings of the National Academy of Sciences, USA. 2006, 103: 6061-6066.

    CAS  Google Scholar 

  160. Willson MF, Whelan CJ: The evolution of fruit color in fleshy-fruited plants. American Naturalist. 1990, 136: 790-809.

    Google Scholar 

  161. Crepet WL, Friis EM, Nixon KC: Fossil evidence for the evolution of biotic pollination. Philosophical Transactions of the Royal Society, London B. 1991, 333: 187-195.

    Google Scholar 

  162. Grimaldi D: The co-radiations of pollinating insects and angiosperms in the Cretaceous. Annals of the Missouri Botanical Garden. 1999, 86: 373-406.

    Google Scholar 

  163. Crepet WL: Progress in understanding angiosperm history, success, and relationships: Darwin's abominably "perplexing phenomenon". Proceedings of the National Academy of Sciences, USA. 2000, 97: 12939-12941.

    CAS  Google Scholar 

  164. Dong X, Braun EL, Grotewold E: Functional conservation of plant secondary metabolic enzymes revealed by complementation of Arabidopsis flavonoid mutants with maize genes. Plant Physiology. 2001, 127: 46-57.

    PubMed Central  CAS  Google Scholar 

  165. Stafford HA: Flavonoid evolution: an enzymic approach. Plant Physiology. 2001, 96: 680-685.

    Google Scholar 

  166. Archetti M, Döring TF, Hagen SB, Hughes NM, Leather SR, Lee DW, Lev-Yadun S, Manetas Y, Ougham HJ, Schaberg PG, Thomas H: Unravelling the evolution of autumn colours – an interdisciplinary approach. Trends in Ecology and Evolution. 2009, 24: 166-173.

    Google Scholar 

  167. Archetti M: The origin of autumn colours by coevolution. Journal of Theoretical Biology. 2000, 205: 625-630.

    CAS  Google Scholar 

  168. Lee DW, Gould KS: Anthocyanins in leaves and other vegetative organs: An introduction. Advances in Botanical Research. 2002, 37: 1-16.

    CAS  Google Scholar 

  169. Ruddiman WF: The anthropogenic greenhouse era began thousands of years ago. Climatic Change. 2003, 61: 261-293.

    CAS  Google Scholar 

  170. Crutzen PJ, Stoermer EF: The "Anthropocene". Global Change Newsletter. 2000, 41: 17-18.

    Google Scholar 

  171. Evans LT: Crop Evolution, Adaptation and Yield. 1993, Cambridge: Cambridge University Press

    Google Scholar 

  172. Thomas H: Canopy survival. Crop Photosynthesis: Spatial and Temporal Determinants. Edited by: Baker N, Thomas H. 1992, Amsterdam: Elsevier, 11-41.

    Google Scholar 

  173. Ooka H, Satoh K, Doi K, Nagata T, Otomo Y, Murakami K, Matsubara K, Osato N, Kawai J, Carninci P, Hayashizaki Y, Suzuki K, Kojima K, Takahara Y, Yamamoto K, Kikuchi S: Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana. DNA Research. 2003, 10: 239-247.

    CAS  Google Scholar 

  174. Hirose N, Makita N, Kojima N, Kamada-Nobusada T, Sakakibara H: Overexpression of a Type-A response regulator alters rice morphology and cytokinin metabolism. Plant and Cell Physiology. 2007, 48: 523-539.

    CAS  Google Scholar 

  175. Thomas H, Smart CM: Crops that stay green. Annals of Applied Biology. 1993, 123: 193-129.

    Google Scholar 

  176. Thomas H, Howarth CJ: Five ways to stay green. Journal of Experimental Botany. 2000, 51: 329-337.

    CAS  Google Scholar 

  177. Giovannoni JJ: Fruit ripening mutants yield insights into ripening control. Current Opinion in Plant Biology. 2007, 10: 283-289.

    CAS  Google Scholar 

  178. Thompson DW: On Growth and Form. 1992, Mineola NY: Dover

    Google Scholar 

  179. van Doorn WG, Woltering EJ: Senescence and programmed cell death: substance or semantics?. Journal of Experimental Botany. 2004, 55: 2147-2153.

    CAS  Google Scholar 

  180. Delph LF, Lively CM: The evolution of floral color change: pollinator attraction versus physiological constraints in Fuchsia excorticata. Evolution. 1989, 43: 1252-1262.

    Google Scholar 

  181. Pyke KA, Page AM: Plastid ontogeny during petal development in Arabidopsis. Plant Physiology. 1998, 116: 797-803.

    PubMed Central  CAS  Google Scholar 

  182. Stead AD, Duckett JG: Plastid ontogeny in the corolla cells of Digitalis purpurea L. cv. Foxy. Annals of Botany. 1980, 46: 549-555.

    Google Scholar 

  183. Schelbert S, Aubry S, Burla B, Agne B, Kessler F, Krupinska K, Hörtensteiner S: Pheophytin pheophorbide hydrolase (pheophytinase) is involved in chlorophyll breakdown during leaf senescence in Arabidopsis. Plant Cell. 2009, 21: 767-785.

    PubMed Central  CAS  Google Scholar 

  184. Breeze E, Harrison E, Page T, Warner N, Shen C, Zhang C, Buchanan-Wollaston V: Transcriptional regulation of plant senescence: from functional genomics to systems biology. Plant Biology. 2008, 10 (supplement 1): 99-109.

    CAS  Google Scholar 

  185. Felsenstein J: Phylogenies from molecular sequences: inference and reliability. Annual Review of Genetics. 1988, 22: 521-565.

    CAS  Google Scholar 

  186. Felsenstein J: PHYLIP – Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.

    Google Scholar 

  187. Yang Z: A space-time process model for the evolution of DNA sequences. Genetics. 1995, 139: 993-1005.

    PubMed Central  CAS  Google Scholar 

  188. Felsenstein J, Churchill GA: A Hidden Markov Model approach to variation among sites in rate of evolution. Molecular Biology of Evolution. 1996, 13: 93-104.

    CAS  Google Scholar 

  189. Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504.

    CAS  Google Scholar 

  190. Schmidt HA, von Haeseler A: Maximum-Likelihood Analysis Using TREE-PUZZLE. Current Protocols in Bioinformatics. Edited by: Baxevanis AD, Davison DB, Page RDM, Stormo G, Stein L. 2003, New York: Wiley and Sons, Unit 6.6

    Google Scholar 

  191. Grennan AK: Genevestigator. Facilitating web-based gene-expression analysis. Plant Physiology. 2006, 141: 1164-1166.

    PubMed Central  CAS  Google Scholar 

  192. Gentleman RC, Carey VJ, Douglas MB, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber H, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang YHY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biology. 2004, 5: R80-

    PubMed Central  Google Scholar 

Download references

Acknowledgements

Howard Thomas is grateful to the Leverhulme Trust for the award of an Emeritus Fellowship. Lin Huang and Helen Ougham acknowledge financial support from a UK Biotechnology and Biological Sciences Research Council Strategic Programme Grant. We thank Joe Ironside for helpful discussion.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Howard Thomas.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HT conceived of the study, participated in sequence alignments and analyses of expression profiles and substantially drafted the manuscript. LH carried out the phylogeny and bootstrap analyses. MY was responsible for alignment and interpretation of Sgr-like sequences in Clostridium. HO devised the bioinformatics strategy, contributed to alignment and expression analyses and carried out the investigations of CCD8, NAC and cyanobacterial RCCR-like sequences. All authors read and approved the final manuscript.

Electronic supplementary material

12862_2008_1077_MOESM1_ESM.doc

Additional file 1: Outfiles from Consensus tree program, version 3.67. Bootstrapping support values for protein consensus trees generated through the use of SEQBOOT, NEIGHBOR and CONSENSE and PUZZLEBOOT programs. (DOC 73 KB)

12862_2008_1077_MOESM2_ESM.doc

Additional file 2: Outfiles from PROML. Amino acid sequence Maximum Likelihood method, version 3.67. Tree branch lengths with error estimates, generated using PROML with Gamma distribution. (DOC 122 KB)

12862_2008_1077_MOESM3_ESM.doc

Additional file 3: Synopsis of significant events in the evolution of plant senescence in relation to the geological timeline. Innovations in the geological time-line relevant to the evolution of a senescence program integrated into plant development and adaptation. (DOC 46 KB)

12862_2008_1077_MOESM4_ESM.doc

Additional file 4: Log Likelihood (Ln) values corresponding to Hidden Markov Model, Gamma Distribution Model and Constant Rate Variation options in PHYLIP for phylogenetic analysis of senescence-related protein sequences. Results of analysing all amino acid sequence data with the PROML program run with three different settings for the R (rate variation among sites) parameter. (DOC 38 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Authors’ original file for figure 8

Authors’ original file for figure 9

Authors’ original file for figure 10

Authors’ original file for figure 11

Authors’ original file for figure 12

Authors’ original file for figure 13

Authors’ original file for figure 14

Authors’ original file for figure 15

Authors’ original file for figure 16

Authors’ original file for figure 17

Authors’ original file for figure 18

Authors’ original file for figure 19

Authors’ original file for figure 20

Authors’ original file for figure 21

Authors’ original file for figure 22

Rights and permissions

This article is published under license to 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.

Reprints and permissions

About this article

Cite this article

Thomas, H., Huang, L., Young, M. et al. Evolution of plant senescence. BMC Evol Biol 9, 163 (2009). https://doi.org/10.1186/1471-2148-9-163

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-9-163

Keywords