Systems genetic analysis of lignin biosynthesis in Populus tremula.

The genetic control underlying natural variation in lignin content and composition in trees is not fully understood. We performed a systems genetic analysis to uncover the genetic regulation of lignin biosynthesis in a natural 'SwAsp' population of aspen (Populus tremula) trees. We analyzed gene expression by RNA sequencing (RNA-seq) in differentiating xylem tissues, and lignin content and composition using Pyrolysis-GC-MS in mature wood of 268 trees from 99 genotypes. Abundant variation was observed for lignin content and composition, and genome-wide association study identified proteins in the pentose phosphate pathway and arabinogalactan protein glycosylation among the top-ranked genes that are associated with these traits. Variation in gene expression and the associated genetic polymorphism was revealed through the identification of 312 705 local and 292 003 distant expression quantitative trait loci (eQTL). A co-expression network analysis suggested modularization of lignin biosynthesis and novel functions for the lignin-biosynthetic CINNAMYL ALCOHOL DEHYDROGENASE 2 and CAFFEOYL-CoA O-METHYLTRANSFERASE 3. PHENYLALANINE AMMONIA LYASE 3 was co-expressed with HOMEOBOX PROTEIN 5 (HB5), and the role of HB5 in stimulating lignification was demonstrated in transgenic trees. The systems genetic approach allowed linking natural variation in lignin biosynthesis to trees´ responses to external cues such as mechanical stimulus and nutrient availability.

The effect of nitrogen source and levels on hybrid aspen tree physiology and wood formation.

Nitrogen can be taken up by trees in the form of nitrate, ammonium and amino acids, but the influence of the different forms on tree growth and development is poorly understood in angiosperm species like Populus. We studied the effects of both organic and inorganic forms of nitrogen on growth and wood formation of hybrid aspen trees in experimental conditions that allowed growth under four distinct steady-state nitrogen levels. Increased nitrogen availability had a positive influence on biomass accumulation and the radial dimensions of both xylem vessels and fibers, and a negative influence on wood density. An optimal level of nitrogen availability was identified where increases in biomass accumulation outweighed decreases in wood density. None of these responses depended on the source of nitrogen except for shoot biomass accumulation, which was stimulated more by treatments complemented with nitrate than by ammonium alone or the organic source arginine. The most striking difference between the nitrogen sources was the effect on lignin composition, whereby the abundance of H-type lignin increased only in the presence of nitrate. The differential effect of nitrate is possibly related to the well-known role of nitrate as a signaling compound. RNA-sequencing revealed that while the lignin-biosynthetic genes did not significantly (FDR <0.01) respond to added NO3 - , the expression of several laccases, catalysing lignin polymerization, was dependent on N-availability. These results reveal a unique role of nitrate in wood formation and contribute to the knowledge basis for decision-making in utilizing hybrid aspen as a bioresource.

A phloem-localized Arabidopsis metacaspase (AtMC3) improves drought tolerance.

Increasing drought phenomena pose a serious threat to agricultural productivity. Although plants have multiple ways to respond to the complexity of drought stress, the underlying mechanisms of stress sensing and signaling remain unclear. The role of the vasculature, in particular the phloem, in facilitating inter-organ communication is critical and poorly understood. Combining genetic, proteomic and physiological approaches, we investigated the role of AtMC3, a phloem-specific member of the metacaspase family, in osmotic stress responses in Arabidopsis thaliana. Analyses of the proteome in plants with altered AtMC3 levels revealed differential abundance of proteins related to osmotic stress pointing into a role of the protein in water-stress-related responses. Overexpression of AtMC3 conferred drought tolerance by enhancing the differentiation of specific vascular tissues and maintaining higher levels of vascular-mediated transportation, while plants lacking the protein showed an impaired response to drought and inability to respond effectively to the hormone abscisic acid. Overall, our data highlight the importance of AtMC3 and vascular plasticity in fine-tuning early drought responses at the whole plant level without affecting growth or yield.

ETHYLENE RESPONSE FACTOR 115 integrates jasmonate and cytokinin signaling machineries to repress adventitious rooting in Arabidopsis.

Adventitious root initiation (ARI) is a de novo organogenesis program and a key adaptive trait in plants. Several hormones regulate ARI but the underlying genetic architecture that integrates the hormonal crosstalk governing this process remains largely elusive. In this study, we use genetics, genome editing, transcriptomics, hormone profiling and cell biological approaches to demonstrate a crucial role played by the APETALA2/ETHYLENE RESPONSE FACTOR 115 transcription factor. We demonstrate that ERF115 functions as a repressor of ARI by activating the cytokinin (CK) signaling machinery. We also demonstrate that ERF115 is transcriptionally activated by jasmonate (JA), an oxylipin-derived phytohormone, which represses ARI in NINJA-dependent and independent manners. Our data indicate that NINJA-dependent JA signaling in pericycle cells blocks early events of ARI. Altogether, our results reveal a previously unreported molecular network involving cooperative crosstalk between JA and CK machineries that represses ARI.

PIRIN2 suppresses S-type lignin accumulation in a noncell-autonomous manner in Arabidopsis xylem elements.

PIRIN (PRN) genes encode cupin domain-containing proteins that function as transcriptional co-regulators in humans but that are poorly described in plants. A previous study in xylogenic cell cultures of Zinnia elegans suggested a role for a PRN protein in lignification. This study aimed to identify the function of Arabidopsis (Arabidopsis thaliana) PRN proteins in lignification of xylem tissues. Chemical composition of the secondary cell walls was analysed in Arabidopsis stems and/or hypocotyls by pyrolysis-gas chromatography/mass spectrometry, 2D-nuclear magnetic resonance and phenolic profiling. Secondary cell walls of individual xylem elements were chemotyped by Fourier transform infrared and Raman microspectroscopy. Arabidopsis PRN2 suppressed accumulation of S-type lignin in Arabidopsis stems and hypocotyls. PRN2 promoter activity and PRN2:GFP fusion protein were localised specifically in cells next to the vessel elements, suggesting a role for PRN2 in noncell-autonomous lignification of xylem vessels. Accordingly, PRN2 modulated lignin chemistry in the secondary cell walls of the neighbouring vessel elements. These results indicate that PRN2 suppresses S-type lignin accumulation in the neighbourhood of xylem vessels to bestow G-type enriched lignin composition on the secondary cell walls of the vessel elements. Gene expression analyses suggested that PRN2 function is mediated by regulation of the expression of the lignin-biosynthetic genes.

The chromatin-modifying protein HUB2 is involved in the regulation of lignin composition in xylem vessels.

PIRIN2 (PRN2) was earlier reported to suppress syringyl (S)-type lignin accumulation of xylem vessels of Arabidopsis thaliana. In the present study, we report yeast two-hybrid results supporting the interaction of PRN2 with HISTONE MONOUBIQUITINATION2 (HUB2) in Arabidopsis. HUB2 has been previously implicated in several plant developmental processes, but not in lignification. Interaction between PRN2 and HUB2 was verified by ß-galactosidase enzymatic and co-immunoprecipitation assays. HUB2 promoted the deposition of S-type lignin in the secondary cell walls of both stem and hypocotyl tissues, as analysed by pyrolysis-GC/MS. Chemical fingerprinting of individual xylem vessel cell walls by Raman and Fourier transform infrared microspectroscopy supported the function of HUB2 in lignin deposition. These results, together with a genetic analysis of the hub2 prn2 double mutant, support the antagonistic function of PRN2 and HUB2 in deposition of S-type lignin. Transcriptome analyses indicated the opposite regulation of the S-type lignin biosynthetic gene FERULATE-5-HYDROXYLASE1 by PRN2 and HUB2 as the underlying mechanism. PRN2 and HUB2 promoter activities co-localized in cells neighbouring the xylem vessel elements, suggesting that the S-type lignin-promoting function of HUB2 is antagonized by PRN2 for the benefit of the guaiacyl (G)-type lignin enrichment of the neighbouring xylem vessel elements.

AspWood: High-Spatial-Resolution Transcriptome Profiles Reveal Uncharacterized Modularity of Wood Formation in Populus tremula.

Trees represent the largest terrestrial carbon sink and a renewable source of ligno-cellulose. There is significant scope for yield and quality improvement in these largely undomesticated species, and efforts to engineer elite varieties will benefit from improved understanding of the transcriptional network underlying cambial growth and wood formation. We generated high-spatial-resolution RNA sequencing data spanning the secondary phloem, vascular cambium, and wood-forming tissues of Populus tremula The transcriptome comprised 28,294 expressed, annotated genes, 78 novel protein-coding genes, and 567 putative long intergenic noncoding RNAs. Most paralogs originating from the Salicaceae whole-genome duplication had diverged expression, with the exception of those highly expressed during secondary cell wall deposition. Coexpression network analyses revealed that regulation of the transcriptome underlying cambial growth and wood formation comprises numerous modules forming a continuum of active processes across the tissues. A comparative analysis revealed that a majority of these modules are conserved in Picea abies The high spatial resolution of our data enabled identification of novel roles for characterized genes involved in xylan and cellulose biosynthesis, regulators of xylem vessel and fiber differentiation and lignification. An associated web resource (AspWood, http://aspwood.popgenie.org) provides interactive tools for exploring the expression profiles and coexpression network.

GRIM REAPER peptide binds to receptor kinase PRK5 to trigger cell death in Arabidopsis.

Recognition of extracellular peptides by plasma membrane-localized receptor proteins is commonly used in signal transduction. In plants, very little is known about how extracellular peptides are processed and activated in order to allow recognition by receptors. Here, we show that induction of cell death in planta by a secreted plant protein GRIM REAPER (GRI) is dependent on the activity of the type II metacaspase METACASPASE-9. GRI is cleaved by METACASPASE-9 in vitro resulting in the release of an 11 amino acid peptide. This peptide bound in vivo to the extracellular domain of the plasma membrane-localized, atypical leucine-rich repeat receptor-like kinase POLLEN-SPECIFIC RECEPTOR-LIKE KINASE 5 (PRK5) and was sufficient to induce oxidative stress/ROS-dependent cell death. This shows a signaling pathway in plants from processing and activation of an extracellular protein to recognition by its receptor.

An AP2/ERF transcription factor ERF139 coordinates xylem cell expansion and secondary cell wall deposition.

Differentiation of xylem elements involves cell expansion, secondary cell wall (SCW) deposition and programmed cell death. Transitions between these phases require strict spatiotemporal control. The function of Populus ERF139 (Potri.013G101100) in xylem differentiation was characterized in transgenic overexpression and dominant repressor lines of ERF139 in hybrid aspen (Populus tremula × tremuloides). Xylem properties, SCW chemistry and downstream targets were analyzed in both types of transgenic trees using microscopy techniques, Fourier transform-infrared spectroscopy, pyrolysis-GC/MS, wet chemistry methods and RNA sequencing. Opposite phenotypes were observed in the secondary xylem vessel sizes and SCW chemistry in the two different types of transgenic trees, supporting the function of ERF139 in suppressing the radial expansion of vessel elements and stimulating accumulation of guaiacyl-type lignin and possibly also xylan. Comparative transcriptomics identified genes related to SCW biosynthesis (LAC5, LBD15, MYB86) and salt and drought stress-responsive genes (ANAC002, ABA1) as potential direct targets of ERF139. The phenotypes of the transgenic trees and the stem expression profiles of ERF139 potential target genes support the role of ERF139 as a transcriptional regulator of xylem cell expansion and SCW formation, possibly in response to osmotic changes of the cells.

Transcriptional Roadmap to Seasonal Variation in Wood Formation of Norway Spruce.

Seasonal cues influence several aspects of the secondary growth of tree stems, including cambial activity, wood chemistry, and transition to latewood formation. We investigated seasonal changes in cambial activity, secondary cell wall formation, and tracheid cell death in woody tissues of Norway spruce (Picea abies) throughout one seasonal cycle. RNA sequencing was performed simultaneously in both the xylem and cambium/phloem tissues of the stem. Principal component analysis revealed gradual shifts in the transcriptomes that followed a chronological order throughout the season. A notable remodeling of the transcriptome was observed in the winter, with many genes having maximal expression during the coldest months of the year. A highly coexpressed set of monolignol biosynthesis genes showed high expression during the period of secondary cell wall formation as well as a second peak in midwinter. This midwinter peak in expression did not trigger lignin deposition, as determined by pyrolysis-gas chromatography/mass spectrometry. Coexpression consensus network analyses suggested the involvement of transcription factors belonging to the ASYMMETRIC LEAVES2/LATERAL ORGAN BOUNDARIES and MYELOBLASTOSIS-HOMEOBOX families in the seasonal control of secondary cell wall formation of tracheids. Interestingly, the lifetime of the latewood tracheids stretched beyond the winter dormancy period, correlating with a lack of cell death-related gene expression. Our transcriptomic analyses combined with phylogenetic and microscopic analyses also identified the cellulose and lignin biosynthetic genes and putative regulators for latewood formation and tracheid cell death in Norway spruce, providing a toolbox for further physiological and functional assays of these important phase transitions.

Extracellular peptide Kratos restricts cell death during vascular development and stress in Arabidopsis.

During plant vascular development, xylem tracheary elements (TEs) form water-conducting, empty pipes by genetically regulated cell death. Cell death is prevented from spreading to non-TEs by unidentified intercellular mechanisms, downstream of METACASPASE9 (MC9)-mediated regulation of autophagy in TEs. Here, we identified differentially abundant extracellular peptides in vascular-differentiating wild-type and MC9-down-regulated Arabidopsis cell suspensions. A peptide named Kratos rescued the abnormally high ectopic non-TE death resulting from either MC9 knockout or TE-specific overexpression of the ATG5 autophagy protein during experimentally induced vascular differentiation in Arabidopsis cotyledons. Kratos also reduced cell death following mechanical damage and extracellular ROS production in Arabidopsis leaves. Stress-induced but not vascular non-TE cell death was enhanced by another identified peptide, named Bia. Bia is therefore reminiscent of several known plant cell death-inducing peptides acting as damage-associated molecular patterns. In contrast, Kratos plays a novel extracellular cell survival role in the context of development and during stress response.

The Norway spruce genome sequence and conifer genome evolution.

Conifers have dominated forests for more than 200 million years and are of huge ecological and economic importance. Here we present the draft assembly of the 20-gigabase genome of Norway spruce (Picea abies), the first available for any gymnosperm. The number of well-supported genes (28,354) is similar to the >100 times smaller genome of Arabidopsis thaliana, and there is no evidence of a recent whole-genome duplication in the gymnosperm lineage. Instead, the large genome size seems to result from the slow and steady accumulation of a diverse set of long-terminal repeat transposable elements, possibly owing to the lack of an efficient elimination mechanism. Comparative sequencing of Pinus sylvestris, Abies sibirica, Juniperus communis, Taxus baccata and Gnetum gnemon reveals that the transposable element diversity is shared among extant conifers. Expression of 24-nucleotide small RNAs, previously implicated in transposable element silencing, is tissue-specific and much lower than in other plants. We further identify numerous long (>10,000 base pairs) introns, gene-like fragments, uncharacterized long non-coding RNAs and short RNAs. This opens up new genomic avenues for conifer forestry and breeding.

A multi-omics approach reveals function of Secretory Carrier-Associated Membrane Proteins in wood formation of​ ​​Populus​​ ​trees.

BACKGROUND: Secretory Carrier-Associated Membrane Proteins (SCAMPs) are highly conserved 32-38 kDa proteins that are involved in membrane trafficking. A systems approach was taken to elucidate function of SCAMPs in wood formation of Populus trees. Phenotypic and multi-omics analyses were performed in woody tissues of transgenic Populus trees carrying an RNAi construct for Populus tremula x tremuloides SCAMP3 (PttSCAMP3; Potri.019G104000). RESULTS: The woody tissues of the transgenic trees displayed increased amounts of both polysaccharides and lignin oligomers, indicating increased deposition of both the carbohydrate and lignin components of the secondary cell walls. This coincided with a tendency towards increased wood density as well as significantly increased thickness of the suberized cork in the transgenic lines. Multivariate OnPLS (orthogonal projections to latent structures) modeling of five different omics datasets (the transcriptome, proteome, GC-MS metabolome, LC-MS metabolome and pyrolysis-GC/MS metabolome) collected from the secondary xylem tissues of the stem revealed systemic variation in the different variables in the transgenic lines, including changes that correlated with the changes in the secondary cell wall composition. The OnPLS model also identified a rather large number of proteins that were more abundant in the transgenic lines than in the wild type. Several of these were related to secretion and/or endocytosis as well as both primary and secondary cell wall biosynthesis. CONCLUSIONS: Populus SCAMP proteins were shown to influence accumulation of secondary cell wall components, including polysaccharides and phenolic compounds, in the woody tissues of Populus tree stems. Our multi-omics analyses combined with the OnPLS modelling suggest that this function is mediated by changes in membrane trafficking to fine-tune the abundance of cell wall precursors and/or proteins involved in cell wall biosynthesis and transport. The data provides a multi-level source of information for future studies on the function of the SCAMP proteins in plant stem tissues.

The function of two type II metacaspases in woody tissues of Populus trees.

Metacaspases (MCs) are cysteine proteases that are implicated in programmed cell death of plants. AtMC9 (Arabidopsis thaliana Metacaspase9) is a member of the Arabidopsis MC family that controls the rapid autolysis of the xylem vessel elements, but its downstream targets in xylem remain uncharacterized. PttMC13 and PttMC14 were identified as AtMC9 homologs in hybrid aspen (Populus tremula × tremuloides). A proteomic analysis was conducted in xylem tissues of transgenic hybrid aspen trees which carried either an overexpression or an RNA interference construct for PttMC13 and PttMC14. The proteomic analysis revealed modulation of levels of both previously known targets of metacaspases, such as Tudor staphylococcal nuclease, heat shock proteins and 14-3-3 proteins, as well as novel proteins, such as homologs of the PUTATIVE ASPARTIC PROTEASE3 (PASPA3) and the cysteine protease RD21 by PttMC13 and PttMC14. We identified here the pathways and processes that are modulated by PttMC13 and PttMC14 in xylem tissues. In particular, the results indicate involvement of PttMC13 and/or PttMC14 in downstream proteolytic processes and cell death of xylem elements. This work provides a valuable reference dataset on xylem-specific metacaspase functions for future functional and biochemical analyses.

NorWood: a gene expression resource for evo-devo studies of conifer wood development.

The secondary xylem of conifers is composed mainly of tracheids that differ anatomically and chemically from angiosperm xylem cells. There is currently no high-spatial-resolution data available profiling gene expression during wood formation for any coniferous species, which limits insight into tracheid development. RNA-sequencing data from replicated, high-spatial-resolution section series throughout the cambial and woody tissues of Picea abies were used to generate the NorWood.conGenIE.org web resource, which facilitates exploration of the associated gene expression profiles and co-expression networks. Integration within PlantGenIE.org enabled a comparative regulomics analysis, revealing divergent co-expression networks between P. abies and the two angiosperm species Arabidopsis thaliana and Populus tremula for the secondary cell wall (SCW) master regulator NAC Class IIB transcription factors. The SCW cellulose synthase genes (CesAs) were located in the neighbourhoods of the NAC factors in A. thaliana and P. tremula, but not in P. abies. The NorWood co-expression network enabled identification of potential SCW CesA regulators in P. abies. The NorWood web resource represents a powerful community tool for generating evo-devo insights into the divergence of wood formation between angiosperms and gymnosperms and for advancing understanding of the regulation of wood development in P. abies.

Non-cell-autonomous postmortem lignification of tracheary elements in Zinnia elegans.

Postmortem lignification of xylem tracheary elements (TEs) has been debated for decades. Here, we provide evidence in Zinnia elegans TE cell cultures, using pharmacological inhibitors and in intact Z. elegans plants using Fourier transform infrared microspectroscopy, that TE lignification occurs postmortem (i.e., after TE programmed cell death). In situ RT-PCR verified expression of the lignin monomer biosynthetic cinnamoyl CoA reductase and cinnamyl alcohol dehydrogenase in not only the lignifying TEs but also in the unlignified non-TE cells of Z. elegans TE cell cultures and in living, parenchymatic xylem cells that surround TEs in stems. These cells were also shown to have the capacity to synthesize and transport lignin monomers and reactive oxygen species to the cell walls of dead TEs. Differential gene expression analysis in Z. elegans TE cell cultures and concomitant functional analysis in Arabidopsis thaliana resulted in identification of several genes that were expressed in the non-TE cells and that affected lignin chemistry on the basis of pyrolysis-gas chromatography/mass spectrometry analysis. These data suggest that living, parenchymatic xylem cells contribute to TE lignification in a non-cell-autonomous manner, thus enabling the postmortem lignification of TEs.

PIRIN2 stabilizes cysteine protease XCP2 and increases susceptibility to the vascular pathogen Ralstonia solanacearum in Arabidopsis.

PIRIN (PRN) is a member of the functionally diverse cupin protein superfamily. There are four members of the Arabidopsis thaliana PRN family, but the roles of these proteins are largely unknown. Here we describe a function of the Arabidopsis PIRIN2 (PRN2) that is related to susceptibility to the bacterial plant pathogen Ralstonia solanacearum. Two prn2 mutant alleles displayed decreased disease development and bacterial growth in response to R.  solanacearum infection. We elucidated the underlying molecular mechanism by analyzing PRN2 interactions with the papain-like cysteine proteases (PLCPs) XCP2, RD21A, and RD21B, all of which bound to PRN2 in yeast two-hybrid assays and in Arabidopsis protoplast co-immunoprecipitation assays. We show that XCP2 is stabilized by PRN2 through inhibition of its autolysis on the basis of PLCP activity profiling assays and enzymatic assays with recombinant protein. The stabilization of XCP2 by PRN2 was also confirmed in planta. Like prn2 mutants, an xcp2 single knockout mutant and xcp2 prn2 double knockout mutant displayed decreased susceptibility to R. solanacearum, suggesting that stabilization of XCP2 by PRN2 underlies susceptibility to R. solanacearum in Arabidopsis.

Thermospermine levels are controlled by an auxin-dependent feedback loop mechanism in Populus xylem.

Polyamines are small polycationic amines that are widespread in living organisms. Thermospermine, synthesized by thermospermine synthase ACAULIS5 (ACL5), was recently shown to be an endogenous plant polyamine. Thermospermine is critical for proper vascular development and xylem cell specification, but it is not known how thermospermine homeostasis is controlled in the xylem. We present data in the Populus model system supporting the existence of a negative feedback control of thermospermine levels in stem xylem tissues, the main site of thermospermine biosynthesis. While over-expression of the ACL5 homologue in Populus, POPACAULIS5, resulted in strong up-regulation of ACL5 expression and thermospermine accumulation in leaves, the corresponding levels in the secondary xylem tissues of the stem were similar or lower than those in the wild-type. POPACAULIS5 over-expression had a negative effect on accumulation of indole-3-acetic acid, while exogenous auxin had a positive effect on POPACAULIS5 expression, thus promoting thermospermine accumulation. Further, over-expression of POPACAULIS5 negatively affected expression of the class III homeodomain leucine zipper (HD-Zip III) transcription factor gene PttHB8, a homologue of AtHB8, while up-regulation of PttHB8 positively affected POPACAULIS5 expression. These results indicate that excessive accumulation of thermospermine is prevented by a negative feedback control of POPACAULIS5 transcript levels through suppression of indole-3-acetic acid levels, and that PttHB8 is involved in the control of POPACAULIS5 expression. We propose that this negative feedback loop functions to maintain steady-state levels of thermospermine, which is required for proper xylem development, and that it is dependent on the presence of high concentrations of endogenous indole-3-acetic acid, such as those present in the secondary xylem tissues.

Programmes of cell death and autolysis in tracheary elements: when a suicidal cell arranges its own corpse removal.

Tracheary element (TE) differentiation represents a unique system to study plant developmental programmed cell death (PCD). TE PCD occurs after deposition of the secondary cell walls when an unknown signal induces tonoplast rupture and the arrest of cytoplasmic streaming. TE PCD is tightly followed by autolysis of the protoplast and partial hydrolysis of the primary cell walls. This review integrates TE differentiation, programmed cell death (PCD), and autolysis in a biological and evolutionary context. The collective evidence from the evolutionary and molecular studies suggests that TE differentiation consists primarily of a programme for cell death and autolysis under the direct control of the transcriptional master switches VASCULAR NAC DOMAIN 6 (VND6) and VND7. In this scenario, secondary cell walls represent a later innovation to improve the water transport capacity of TEs which necessitates transcriptional regulators downstream of VND6 and VND7. One of the most fascinating features of TEs is that they need to prepare their own corpse removal by expression and accumulation of hydrolases that are released from the vacuole after TE cell death. Therefore, TE differentiation involves, in addition to PCD, a programmed autolysis which is initiated before cell death and executed post-mortem. It has recently become clear that TE PCD and autolysis are separate processes with separate molecular regulation. Therefore, the importance of distinguishing between the cell death programme per se and autolysis in all plant PCD research and of careful description of the morphological, biochemical, and molecular sequences in each of these processes, is advocated.