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 , ) 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 ; [
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 ). 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 ). 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 , ). 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 ) and to components such as See2 (Figures , ) 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 ).
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 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 ). BLAST identifies algal sequences with similarity to AtMRP2, and comparisons of domain structure suggest common functions (Figure ). 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 CO
2 assimilation within chloroplasts are regulated by systems present in
Chlamydomonas and other unicellular chlorophyceans [
45]. Atmospheric CO
2 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 ) 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 ), 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 ). 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 CO
2 concentrations at the time when plants invaded the land persisted for 40–50 my, until the late Devonian, when a steep fall in CO
2 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 CO
2 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 N
2 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 ). 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 ), 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 H
2O
2, 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 ). Other genes listed in Figure and profiled in Figure 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 ).
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 ) 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 ) 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 , Figure ; [
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 ) 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 ). 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 ). 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 ) 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.
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 ). In this regard it is significant that genes for Fibrillin, OrI and OrII, which originated phylogenetically long before chromoplasts had evolved (Figures , , are close to each other when clustered by expression profile (Figure , 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 ) 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 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 ), 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 ), 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 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 ) and the phylogenetic diversity of senescence-related genes (Figure ) 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.