High-throughput single-cell RNA sequencing (scRNA-seq) has advantages over traditional RNA-seq to explore spatiotemporal information on gene dynamic expressions in heterogenous tissues. We performed Drop-seq, a method for the dropwise sequestration of single cells for sequencing, on protoplasts from the differentiating xylem of Populus alba × Populus glandulosa. The scRNA-seq profiled 9,798 cells, which were grouped into 12 clusters. Through characterization of differentially expressed genes in each cluster and RNA in situ hybridizations, we identified vessel cells, fiber cells, ray parenchyma cells and xylem precursor cells. Diffusion pseudotime analyses revealed the differentiating trajectory of vessels, fiber cells and ray parenchyma cells and indicated a different differentiation process between vessels and fiber cells, and a similar differentiation process between fiber cells and ray parenchyma cells. We identified marker genes for each cell type (cluster) and key candidate regulators during developmental stages of xylem cell differentiation. Our study generates a high-resolution expression atlas of wood formation at the single cell level and provides valuable information on wood formation.

Deciphering protein‐protein interactions (PPIs) is fundamental for understanding signal transduction pathways in plants. The split firefly luciferase (Fluc) complementation (SLC) assay has been widely used for analyzing PPIs. However, concern has risen about the bulky halves of Fluc interfering with the functions of their fusion partners. Nano luciferase (Nluc) is the smallest substitute for Fluc with improved stability and luminescence. Here, we developed a dual‐use system enabling the detection of PPIs through the Nluc‐based SLC and co‐immunoprecipitation assays. This was realized by coexpression of two proteins under investigation in fusion with the HA‐ or FLAG‐tagged Nluc halves, respectively. We validated the robustness of this system by reproducing multiple previously documented PPIs in protoplasts or Agrobacterium‐transformed plants. We next applied this system to evaluate the homodimerization of Arabidopsis CERK1, a coreceptor of fungal elicitor chitin, and its heterodimerization with other homologs in the absence or presence of chitin. Moreover, split fragments of Nluc were fused to two cytosolic ends of Arabidopsis calcium channels CNGC2 and CNGC4 to help sense the allosteric change induced by the bacterial elicitor flg22. Collectively, these results demonstrate the usefulness of the Nluc‐based SLC assay for probing constitutive or inducible PPIs and protein allostery in plant cells.

The tremendous functional, spatial, and temporal diversity of the plant proteome is regulated by multiple factors that continuously modify protein abundance, modifications, interactions, localization, and activity to meet the dynamic needs of plants. Dissecting the proteome complexity and its underlying genetic variation is attracting increasing research attention. Mass spectrometry (MS)-based proteomics has become a powerful approach in the global study of protein functions and their relationships on a systems level. Here, we review recent breakthroughs and strategies adopted to unravel the diversity of the proteome, with a specific focus on the methods used to analyze posttranslational modifications (PTMs), protein localization, and the organization of proteins into functional modules. We also consider PTM crosstalk and multiple PTMs temporally regulating the life cycle of proteins. Finally, we discuss recent quantitative studies using MS to measure protein turnover rates and examine future directions in the study of the plant proteome.

In eukaryotes, MEDIATOR is a conserved multi‐subunit complex that links transcription factors and RNA polymerase II and that thereby facilitates transcriptional initiation. Although the composition of MEDIATOR has been well studied in yeast and mammals, relatively little is known about the composition of MEDIATOR in plants. By affinity purification followed by mass spectrometry, we identified 28 conserved MEDIATOR subunits in Arabidopsis thaliana, including putative MEDIATOR subunits that were not previously validated. Our results indicated that MED34, MED35, MED36, and MED37 are not Arabidopsis MEDIATOR subunits, as previously proposed. Our results also revealed that two homologous CBP/p300 histone acetyltransferases, HAC1 and HAC5 (HAC1/5) are in fact plant‐specific MEDIATOR subunits. The MEDIATOR subunits MED8 and MED25 (MED8/25) are partially responsible for the association of MEDIATOR with HAC1/5, MED8/25 and HAC1/5 co‐regulate gene expression and thereby affect flowering time and floral development. Our in vitro observations indicated that MED8 and HAC1 form liquid‐like droplets by phase separation, and our in vivo observations indicated that these droplets co‐localize in the nuclear bodies at a subset of nuclei. The formation of liquid‐like droplets is required for MED8 to interact with RNA polymerase II. In summary, we have identified all of the components of Arabidopsis MEDIATOR and revealed the mechanism underlying the link of histone acetylation and transcriptional regulation.

Circular RNAs (circRNAs) are a recently discovered type of non-coding RNA derived from pre-mRNAs. R-loops consist of a DNA:RNA hybrid and the associated single-stranded DNA. In Arabidopsis thaliana, circRNA:DNA R-loops regulate alternative splicing (AS) of SEPALLATA3 (SEP3). However, the occurrence and functions of circRNAs and R-loops in Populus trichocarpa are largely unexplored. Here, we performed circRNA-enriched sequencing in the stem-differentiating xylem (SDX) of P. trichocarpa and identified 2,742 distinct circRNAs, including circ-CESA4, circ-IRX7, and circ-GUX1, which are generated from genes involved in cellulose, and hemicellulose biosynthesis, respectively. To investigate the roles of circRNAs in modulating alternative splicing (AS), we detected 7,836 AS events using PacBio Iso-Seq and identified 634 circRNAs that overlapped with 699 AS events. Furthermore, using DNA:RNA hybrid immunoprecipitation followed by sequencing (DRIP-seq), we identified 8,932 R-loop peaks that overlapped with 181 circRNAs and 672 AS events. Notably, several SDX-related circRNAs overlapped with R-loop peaks, pointing to their possible roles in modulating AS in SDX. Indeed, overexpressing circ-IRX7 increased the levels of R-loop structures and decreased the frequency of intron retention in linear IRX7 transcripts. This study provides a valuable R-loop atlas resource and uncovers the interplay between circRNAs and AS in SDX of P. trichocarpa.

Influenza epidemics frequently and unpredictably break out all over the world, and seriously affect the breeding industry and human activity. Inactivated and live attenuated viruses have been used as protective vaccines but exhibit high risks for biosafety. Subunit vaccines enjoy high biosafety and specificity but have a few weak points compared to inactivated virus or live attenuated virus vaccines, especially in low immunogenicity. In this study, we developed a new subunit vaccine platform for a potent, adjuvant-free, and multivalent vaccination. The ectodomains of hemagglutinins (HAs) of influenza viruses were expressed in plants as trimers (tHAs) to mimic their native forms. tHAs in plant extracts were directly used without purification for binding to inactivated Lactococcus (iLact) to produce iLact-tHAs, an antigen-carrying bacteria-like particle (BLP). tHAs BLP showed strong immune responses in mice and chickens without adjuvants. Moreover, simultaneous injection of two different antigens by two different formulas, tHA^{H5N6 + H9N2} BLP or a combination of tHA^H5N6 BLP and tHA^H9N2 BLP, led to strong immune responses to both antigens. Based on these results, we propose combinations of plant-based antigen production and BLP-based delivery as a highly potent and cost-effective platform for multivalent vaccination for subunit vaccines.

Lonicera japonica Thunb., a traditional Chinese herb, has been used for treating human diseases for thousands of years. Recently, the genome of L. japonica has been decoded, providing valuable information for research into gene function. However, no comprehensive database for gene functional analysis and mining is available for L. japonica. We therefore constructed LjaFGD (www.gzybioinformatics.cn/LjaFGD and bioinformatics.cau.edu.cn/LjaFGD), a database for analyzing and comparing gene function in L. japonica. We constructed a gene co-expression network based on 77 RNA-seq samples, and then annotated genes of L. japonica by alignment against protein sequences from public databases. We also introduced several tools for gene functional analysis, including Blast, motif analysis, gene set enrichment analysis, heatmap analysis, and JBrowse. Our co-expression network revealed that MYB and WRKY transcription factor family genes were co-expressed with genes encoding key enzymes in the biosynthesis of chlorogenic acid and luteolin in L. japonica. We used flavonol synthase 1 (LjFLS1) as an example to show the reliability and applicability of our database. LjaFGD and its various associated tools will provide researchers with an accessible platform for retrieving functional information on L. japonica genes to further biological discovery.