Long-read transcriptome sequencing provides insight into lignan biosynthesis during fruit development in Schisandra chinensis

Background Schisandra chinensis, an ancient member of the most basal angiosperm lineage which is known as the ANITA, is a fruit-bearing vine with the pharmacological effects of a multidrug system, such as antioxidant, anti-inflammatory, cardioprotective, neuroprotective, anti-osteoporosis effects. Its major bioactive compound is represented by lignans such as schisandrin. Molecular characterization of lignan biosynthesis in S. chinensis is of great importance for improving the production of this class of active compound. However, the biosynthetic mechanism of schisandrin remains largely unknown. Results To understand the potential key catalytic steps and their regulation of schisandrin biosynthesis, we generated genome-wide transcriptome data from three different tissues of S. chinensis cultivar Cheongsoon, including leaf, root, and fruit, via long- and short-read sequencing technologies. A total of 132,856 assembled transcripts were generated with an average length of 1.9 kb and high assembly completeness. Overall, our data presented effective, accurate gene annotation in the prediction of functional pathways. In particular, the annotation revealed the abundance of transcripts related to phenylpropanoid biosynthesis. Remarkably, transcriptome profiling during fruit development of S. chinensis cultivar Cheongsoon revealed that the phenylpropanoid biosynthetic pathway, specific to coniferyl alcohol biosynthesis, showed a tendency to be upregulated at the postfruit development stage. Further the analysis also revealed that the pathway forms a transcriptional network with fruit ripening-related genes, especially the ABA signaling-related pathway. Finally, candidate unigenes homologous to isoeugenol synthase 1 (IGS1) and dirigent-like protein (DIR), which are subsequently activated by phenylpropanoid biosynthesis and thus catalyze key upstream steps in schisandrin biosynthesis, were identified. Their expression was increased at the postfruit development stage, suggesting that they may be involved in the regulation of schisandrin biosynthesis in S. chinensis. Conclusions Our results provide new insights into the production and accumulation of schisandrin in S. chinensis berries and will be utilized as a valuable transcriptomic resource for improving the schisandrin content.

Evolution of 14-3-3 Proteins in Angiosperm Plants: Recurring Gene Duplication and Loss

14-3-3 proteins are key regulatory factors in plants and are involved in a broad range of physiological processes. We addressed the evolutionary history of 14-3-3s from 46 angiosperm species, including basal angiosperm Amborella and major lineage of monocotyledons and eudicotyledons. Orthologs of Arabidopsis isoforms were detected. There were several rounds of duplication events in the evolutionary history of the 14-3-3 protein family in plants. At least four subfamilies (iota, epsilon, kappa, and psi) formed as a result of ancient duplication in a common ancestor of angiosperm plants. Recent duplication events followed by gene loss in plant lineage, among others Brassicaceae, Fabaceae, and Poaceae, further shaped the high diversity of 14-3-3 isoforms in plants. Coexpression data showed that 14-3-3 proteins formed different functional groups in different species. In some species, evolutionarily related groups of 14-3-3 proteins had coexpressed together under certain physiological conditions, whereas in other species, closely related isoforms expressed in the opposite manner. A possible explanation is that gene duplication and loss is accompanied by functional plasticity of 14-3-3 proteins.

Core circadian clock and light signalling genes brought into genetic linkage across the green lineage

The circadian clock ensures that biological processes are phased to the correct time of day. In plants the circadian clock is conserved at both the level of transcriptional networks as well as core genes. In the model plant Arabidopsis thaliana, the core circadian singleMYB (sMYB) genes CCA1 and RVE4 are in genetic linkage with the PSEUDO-RESPONSE REGULATOR (PRR) genes PRR9 and PRR7 respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis it was possible to trace this genetic linkage back to the basal angiosperm Amborella and identify an additional evolutionarily conserved genetic linkage between PIF3 and PHYA. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7 and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of higher plant families representing distinct whole genome duplication and fractionation events. Soybean maintains all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae) such as maize where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The evolutionary conservation of the genetic linkage as well as its loss in the grasses provides new insight in the plant circadian clock, which has been a critical target of breeding and domestication.

Single-cell transcriptomics unveils xylem cell development and evolution

As the most abundant tissue on Earth, xylem is responsible for lateral growth in plants. Typical xylem has a radial system composed of ray parenchyma cells and an axial system of fusiform cells. In most angiosperms, fusiform cells are a combination of vessel elements for water transportation and libriform fibers for mechanical support, while both functions are performed together by tracheids in other vascular plants. However, little is known about the developmental programs and evolutionary relationships of these xylem cell types. Through both single-cell and laser-capture microdissection transcriptomic profiling, here we demonstrate the developmental lineages of ray and fusiform cells in stem-differentiating xylem across four divergent woody angiosperms. Cross-species analyses of single-cell trajectories reveal highly conserved ray, yet variable fusiform, lineages across angiosperms. Core eudicots Populus trichocarpa and Eucalyptus grandis share nearly identical fusiform lineages. The tracheids in the basal eudicot Trochodendron aralioides, an evolutionarily reversed character, exhibit strong transcriptomic similarity to vessel elements but not libriform fibers, suggesting that water transportation, instead of mechanical support, is the major feature. We also found that the more basal angiosperm Liriodendron chinense has a fusiform lineage distinct from that in core eudicots. This evo-developmental framework provides a comprehensive understanding of the formation of xylem cell lineages across multiple plant species spanning over a hundred million years of evolutionary history.

Organelle Genomes and Transcriptomes of Nymphaea Reveal the Interplay between Intron Splicing and RNA Editing

Posttranscriptional modifications, including intron splicing and RNA editing, are common processes during regulation of gene expression in plant organelle genomes. However, the intermediate products of intron-splicing, and the interplay between intron-splicing and RNA-editing were not well studied. Most organelle transcriptome analyses were based on the Illumina short reads which were unable to capture the full spectrum of transcript intermediates within an organelle. To fully investigate the intermediates during intron splicing and the underlying relationships with RNA editing, we used PacBio DNA-seq and Iso-seq, together with Illumina short reads genome and transcriptome sequencing data to assemble the chloroplast and mitochondrial genomes of Nymphaea ‘Joey Tomocik’ and analyze their posttranscriptional features. With the direct evidence from Iso-seq, multiple intermediates partially or fully intron-spliced were observed, and we also found that both cis- and trans-splicing introns were spliced randomly. Moreover, by using rRNA-depleted and non-Oligo(dT)-enrichment strand-specific RNA-seq data and combining direct SNP-calling and transcript-mapping methods, we identified 98 and 865 RNA-editing sites in the plastome and mitogenome of N. ‘Joey Tomocik’, respectively. The target codon preference, the tendency of increasing protein hydrophobicity, and the bias distribution of editing sites are similar in both organelles, suggesting their common evolutionary origin and shared editing machinery. The distribution of RNA editing sites also implies that the RNA editing sites in the intron and exon regions may splice synchronously, except those exonic sites adjacent to intron which could only be edited after being intron-spliced. Our study provides solid evidence for the multiple intermediates co-existing during intron-splicing and their interplay with RNA editing in organelle genomes of a basal angiosperm.

Functional and Predictive Structural Characterization of WRINKLED2, A Unique Oil Biosynthesis Regulator in Avocado

WRINKLED1 (WRI1), a member of the APETALA2 (AP2) class of transcription factors regulates fatty acid biosynthesis and triacylglycerol (TAG) accumulation in plants. Among the four known Arabidopsis WRI1 paralogs, only WRI2 was unable to complement and restore fatty acid content in wri1-1 mutant seeds. Avocado (Persea americana) mesocarp, which accumulates 60-70% dry weight oil content, showed high expression levels for orthologs of WRI2, along with WRI1 and WRI3, during fruit development. While the role of WRI1 as a master regulator of oil biosynthesis is well-established, the function of WRI1 paralogs is poorly understood. Comprehensive and comparative in silico analyses of WRI1 paralogs from avocado (a basal angiosperm) with higher angiosperms Arabidopsis (dicot), maize (monocot) revealed distinct features. Predictive structural analyses of the WRI orthologs from these three species revealed the presence of AP2 domains and other highly conserved features, such as intrinsically disordered regions associated with predicted PEST motifs and phosphorylation sites. Additionally, avocado WRI proteins also contained distinct features that were absent in the nonfunctional Arabidopsis ortholog AtWRI2. Through transient expression assays, we demonstrated that both avocado WRI1 and WRI2 are functional and drive TAG accumulation in Nicotiana benthamiana leaves. We predict that the unique features and activities of ancestral PaWRI2 were likely lost in orthologous genes such as AtWRI2 during evolution and speciation, leading to at least partial loss of function in some higher eudicots. This study provides us with new targets to enhance oil biosynthesis in plants.

Diversified amino acid-mediated allosteric regulation of phosphoglycerate dehydrogenase for serine biosynthesis in land plants

The phosphorylated pathway of serine biosynthesis is initiated with 3-phosphoglycerate dehydrogenase (PGDH). The liverwort Marchantia polymorpha possesses an amino acid-sensitive MpPGDH which is inhibited by L-serine and activated by five proteinogenic amino acids, while the eudicot Arabidopsis thaliana has amino acid-sensitive AtPGDH1 and AtPGDH3 as well as amino acid-insensitive AtPGDH2. In this study, we analyzed PGDH isozymes of the representative land plants: the monocot Oryza sativa (OsPGDH1‒3), basal angiosperm Amborella trichopoda (AmtriPGDH1‒2), and moss Physcomitrium (Physcomitrella) patens (PpPGDH1‒4). We demonstrated that OsPGDH1, AmtriPGDH1, PpPGDH1, and PpPGDH3 were amino acid-sensitive, whereas OsPGDH2, OsPGDH3, AmtriPGDH2, PpPGDH2 and PpPGDH4 were either sensitive to only some of the six effector amino acids or insensitive to all effectors. This indicates that PGDH sensitivity to effectors has been diversified among isozymes and that the land plant species examined, except for M. polymorpha, possess different isozyme types in terms of regulation. Phylogenetic analysis suggested that the different sensitivities convergently evolved in the bryophyte and angiosperm lineages. Site-directed mutagenesis of AtPGDH1 revealed that Asp538 and Asn556 residues in the ACT domain are involved in allosteric regulation by the effectors. These findings provide insight into the evolution of PGDH isozymes, highlighting the functional diversification of allosteric regulation in land plants.

A Reappraisal of the Evolutionary and Developmental Pathway of Apomixis and Its Genetic Control in Angiosperms

Apomixis sensu stricto (agamospermy) is asexual reproduction by seed. In angiosperms it represents an easy byway of life cycle renewal through gamete-like cells that give rise to maternal embryos without ploidy reduction (meiosis) and ploidy restitution (syngamy). The origin of apomixis still represents an unsolved problem, as it may be either evolved from sex or the other way around. This review deals with a reappraisal of the origin of apomixis in order to deepen knowledge on such asexual mode of reproduction which seems mainly lacking in the most basal angiosperm orders (i.e., Amborellales, Nymphaeales and Austrobaileyales, also known as ANA-grade), while it clearly occurs in different forms and variants in many unrelated families of monocots and eudicots. Overall findings strengthen the hypothesis that apomixis as a whole may have evolved multiple times in angiosperm evolution following different developmental pathways deviating to different extents from sexuality. Recent developments on the genetic control of apomixis in model species are also presented and adequately discussed in order to shed additional light on the antagonist theories of gain- and loss-of-function over sexuality.