Merit of integrating in situ transcriptomics and anatomical information for cell annotation and lineage construction in single-cell analyses of Populus.

Cell type annotation and lineage construction are two of the most critical tasks conducted in the analyses of single-cell RNA sequencing (scRNA-seq). Four recent scRNA-seq studies of differentiating xylem propose four models on differentiating xylem development in Populus. The differences are mostly caused by the use of different strategies for cell type annotation and subsequent lineage interpretation. Here, we emphasize the necessity of using in situ transcriptomes and anatomical information to construct the most plausible xylem development model.

Multipathway regulation induced by 4-(phenylsulfonyl)morpholine derivatives against triple-negative breast cancer.

Phenotypic drug discovery (PDD) is an effective drug discovery approach by observation of therapeutic effects on disease phenotypes, especially in complex disease systems. Triple-negative breast cancer (TNBC) is composed of several complex disease features, including high tumor heterogeneity, high invasive and metastatic potential, and a lack of effective therapeutic targets. Therefore, identifying effective and novel agents through PDD is a current trend in TNBC drug development. In this study, 23 novel small molecules were synthesized using 4-(phenylsulfonyl)morpholine as a pharmacophore. Among these derivatives, GL24 (4m) exhibited the lowest half-maximal inhibitory concentration value (0.90 µM) in MDA-MB-231 cells. To investigate the tumor-suppressive mechanisms of GL24, transcriptomic analyses were used to detect the perturbation for gene expression upon GL24 treatment. Followed by gene ontology (GO) analysis, gene set enrichment analysis (GSEA), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, multiple ER stress-dependent tumor suppressive signals were identified, such as unfolded protein response (UPR), p53 pathway, G2/M checkpoint, and E2F targets. Most of the identified pathways triggered by GL24 eventually led to cell-cycle arrest and then to apoptosis. In summary, we developed a novel 4-(phenylsulfonyl)morpholine derivative GL24 with a strong potential for inhibiting TNBC cell growth through ER stress-dependent tumor suppressive signals.

Woody plant cell walls: Fundamentals and utilization.

Cell walls in plants, particularly forest trees, are the major carbon sink of the terrestrial ecosystem. Chemical and biosynthetic features of plant cell walls were revealed early on, focusing mostly on herbaceous model species. Recent developments in genomics, transcriptomics, epigenomics, transgenesis, and associated analytical techniques are enabling novel insights into formation of woody cell walls. Here, we review multilevel regulation of cell wall biosynthesis in forest tree species. We highlight current approaches to engineering cell walls as potential feedstock for materials and energy and survey reported field tests of such engineered transgenic trees. We outline opportunities and challenges in future research to better understand cell type biogenesis for more efficient wood cell wall modification and utilization for biomaterials or for enhanced carbon capture and storage.

Wood cellulose microfibrils have a 24-chain core-shell nanostructure in seed plants.

Wood cellulose microfibril (CMF) is the most abundant organic substance on Earth but its nanostructure remains poorly understood. There are controversies regarding the glucan chain number (N) of CMFs during initial synthesis and whether they become fused afterward. Here, we combined small-angle X-ray scattering, solid-state nuclear magnetic resonance and X-ray diffraction analyses to resolve CMF nanostructures in native wood. We developed small-angle X-ray scattering measurement methods for the cross-section aspect ratio and area of the crystalline-ordered CMF core, which has a higher scattering length density than the semidisordered shell zone. The 1:1 aspect ratio suggested that CMFs remain mostly segregated, not fused. The area measurement reflected the chain number in the core zone (Ncore). To measure the ratio of ordered cellulose over total cellulose (Roc) by solid-state nuclear magnetic resonance, we developed a method termed global iterative fitting of T1ρ-edited decay (GIFTED), in addition to the conventional proton spin relaxation editing method. Using the formula N = Ncore/Roc, most wood CMFs were found to contain 24 glucan chains, conserved between gymnosperm and angiosperm trees. The average CMF has a crystalline-ordered core of ~2.2 nm diameter and a semidisordered shell of ~0.5 nm thickness. In naturally and artificially aged wood, we observed only CMF aggregation (contact without crystalline continuity) but not fusion (forming a conjoined crystalline unit). This further argued against the existence of partially fused CMFs in new wood, overturning the recently proposed 18-chain fusion hypothesis. Our findings are important for advancing wood structural knowledge and more efficient use of wood resources in sustainable bio-economies.

Single-cell transcriptomics unveils xylem cell development and evolution.

BACKGROUND: Xylem, the most abundant tissue on Earth, 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 comprise vessel elements for water transportation and libriform fibers for mechanical support, while both functions are performed by tracheids in other vascular plants such as gymnosperms. Little is known about the developmental programs and evolutionary relationships of these xylem cell types. RESULTS: Through both single-cell and laser capture microdissection transcriptomic profiling, we determine the developmental lineages of ray and fusiform cells in stem-differentiating xylem across four divergent woody angiosperms. Based on cross-species analyses of single-cell clusters and overlapping trajectories, we reveal highly conserved ray, yet variable fusiform, lineages across angiosperms. Core eudicots Populus trichocarpa and Eucalyptus grandis share nearly identical fusiform lineages, whereas the more basal angiosperm Liriodendron chinense has a fusiform lineage distinct from that in core eudicots. The tracheids in the basal eudicot Trochodendron aralioides, an evolutionarily reversed trait, exhibit strong transcriptomic similarity to vessel elements rather than libriform fibers. CONCLUSIONS: This evo-devo 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.

Cell-type-specific PtrWOX4a and PtrVCS2 form a regulatory nexus with a histone modification system for stem cambium development in Populus trichocarpa.

Stem vascular cambium cells in forest trees produce wood for materials and energy. WOX4 affects the proliferation of such cells in Populus. Here we show that PtrWOX4a is the most highly expressed stem vascular-cambium-specific (VCS) gene in P. trichocarpa, and its expression is controlled by the product of the second most highly expressed VCS gene, PtrVCS2, encoding a zinc finger protein. PtrVCS2 binds to the PtrWOX4a promoter as part of a PtrWOX13a-PtrVCS2-PtrGCN5-1-PtrADA2b-3 protein tetramer. PtrVCS2 prevented the interaction between PtrGCN5-1 and PtrADA2b-3, resulting in H3K9, H3K14 and H3K27 hypoacetylation at the PtrWOX4a promoter, which led to fewer cambium cell layers. These effects on cambium cell proliferation were consistent across more than 20 sets of transgenic lines overexpressing individual genes, gene-edited mutants and RNA interference lines in P. trichocarpa. We propose that the tetramer-PtrWOX4a system may coordinate genetic and epigenetic regulation to maintain normal vascular cambium development for wood formation.

Development of a petal protoplast transfection system for Sinningia speciosa.

Premise: Transient gene expression systems are powerful tools for studying gene interactions in plant species without available or stable genetic transformation protocols. We optimized a petal protoplast transformation protocol for Sinningia speciosa, a model plant, to study the development of floral symmetry. Methods and Results: A high yield of petal protoplasts was obtained using a 6-h enzyme digestion in a solution of 1.5% cellulase and 0.4% macerozyme. Modest transfection efficiency (average 41.4%) was achieved. The viability of the transfected protoplasts remained at more than 90%. A fusion of green fluorescent protein and CYCLOIDEA (SsCYC), the Teosinte branched 1/Cincinnata/Proliferating cell factor transcription factor responsible for floral symmetry, was subcellularly localized inside the nuclei of the protoplasts. Transiently overexpressing SsCYC indicates the success of this system, which resulted in the predicted increased (but nonsignificant) expression of its known target RADIALIS (SsRAD1), consistent with gene network expectations. Conclusions: The transient transfection system presented herein can be effectively used to study gene-regulatory interactions in Gesneriaceae species.

A novel AMPK activator shows therapeutic potential in hepatocellular carcinoma by suppressing HIF1α-mediated aerobic glycolysis.

Hepatocellular carcinoma (HCC) is characterized by rapid growth, early vascular invasion, and high metastasis. Currently available US Food and Drug Administration (FDA)-approved drugs show low therapeutic efficacy, limiting HCC treatment to chemotherapy. We designed and synthesized a novel small molecule, SCT-1015, that allosterically activated adenosine monophosphate-activated protein kinase (AMPK) to suppress the aerobic glycolysis in HCC. SCT-1015 was shown to bind the AMPK α and ß-subunit interface, thereby exposing the kinase α domain to the upstream kinases, resulting in the increased AMPK activity. SCT-1015 dramatically reduced HCC cell growth in vitro and tumor growth in vivo. We further found that AMPK formed protein complexes with hypoxia-inducible factor 1-alpha (HIF1α) and that SCT-1015-activated AMPK promoted hydroxylation of HIF1α (402P and 564P), resulting in HIF1α degradation by the ubiquitin-proteasome system. With declined HIF1α abundance, many glycolysis-related enzymes were downregulated, suppressing aerobic glycolysis, and promoting oxidative phosphorylation. These results indicated that SCT-1015 channeled HCC cells into an unfavorable metabolic status. Overall, we reported SCT-1015 as a direct activator of AMPK signaling that held therapeutic potential in HCC.

Enzyme Complexes of Ptr4CL and PtrHCT Modulate Co-enzyme A Ligation of Hydroxycinnamic Acids for Monolignol Biosynthesis in Populus trichocarpa.

Co-enzyme A (CoA) ligation of hydroxycinnamic acids by 4-coumaric acid:CoA ligase (4CL) is a critical step in the biosynthesis of monolignols. Perturbation of 4CL activity significantly impacts the lignin content of diverse plant species. In Populus trichocarpa, two well-studied xylem-specific Ptr4CLs (Ptr4CL3 and Ptr4CL5) catalyze the CoA ligation of 4-coumaric acid to 4-coumaroyl-CoA and caffeic acid to caffeoyl-CoA. Subsequently, two 4-hydroxycinnamoyl-CoA:shikimic acid hydroxycinnamoyl transferases (PtrHCT1 and PtrHCT6) mediate the conversion of 4-coumaroyl-CoA to caffeoyl-CoA. Here, we show that the CoA ligation of 4-coumaric and caffeic acids is modulated by Ptr4CL/PtrHCT protein complexes. Downregulation of PtrHCTs reduced Ptr4CL activities in the stem-differentiating xylem (SDX) of transgenic P. trichocarpa. The Ptr4CL/PtrHCT interactions were then validated in vivo using biomolecular fluorescence complementation (BiFC) and protein pull-down assays in P. trichocarpa SDX extracts. Enzyme activity assays using recombinant proteins of Ptr4CL and PtrHCT showed elevated CoA ligation activity for Ptr4CL when supplemented with PtrHCT. Numerical analyses based on an evolutionary computation of the CoA ligation activity estimated the stoichiometry of the protein complex to consist of one Ptr4CL and two PtrHCTs, which was experimentally confirmed by chemical cross-linking using SDX plant protein extracts and recombinant proteins. Based on these results, we propose that Ptr4CL/PtrHCT complexes modulate the metabolic flux of CoA ligation for monolignol biosynthesis during wood formation in P. trichocarpa.

Large-scale data analysis for robotic yeast one-hybrid platforms and multi-disciplinary studies using GateMultiplex.

BACKGROUND: Yeast one-hybrid (Y1H) is a common technique for identifying DNA-protein interactions, and robotic platforms have been developed for high-throughput analyses to unravel the gene regulatory networks in many organisms. Use of these high-throughput techniques has led to the generation of increasingly large datasets, and several software packages have been developed to analyze such data. We previously established the currently most efficient Y1H system, meiosis-directed Y1H; however, the available software tools were not designed for processing the additional parameters suggested by meiosis-directed Y1H to avoid false positives and required programming skills for operation. RESULTS: We developed a new tool named GateMultiplex with high computing performance using C++. GateMultiplex incorporated a graphical user interface (GUI), which allows the operation without any programming skills. Flexible parameter options were designed for multiple experimental purposes to enable the application of GateMultiplex even beyond Y1H platforms. We further demonstrated the data analysis from other three fields using GateMultiplex, the identification of lead compounds in preclinical cancer drug discovery, the crop line selection in precision agriculture, and the ocean pollution detection from deep-sea fishery. CONCLUSIONS: The user-friendly GUI, fast C++ computing speed, flexible parameter setting, and applicability of GateMultiplex facilitate the feasibility of large-scale data analysis in life science fields.

Transcriptional reprogramming of xylem cell wall biosynthesis in tension wood.

Tension wood (TW) is a specialized xylem tissue developed under mechanical/tension stress in angiosperm trees. TW development involves transregulation of secondary cell wall genes, which leads to altered wood properties for stress adaptation. We induced TW in the stems of black cottonwood (Populus trichocarpa, Nisqually-1) and identified two significantly repressed transcription factor (TF) genes: class B3 heat-shock TF (HSFB3-1) and MYB092. Transcriptomic analysis and chromatin immunoprecipitation (ChIP) were used to identify direct TF-DNA interactions in P. trichocarpa xylem protoplasts overexpressing the TFs. This analysis established a transcriptional regulatory network in which PtrHSFB3-1 and PtrMYB092 directly activate 8 and 11 monolignol genes, respectively. The TF-DNA interactions were verified for their specificity and transactivator roles in 35 independent CRISPR-based biallelic mutants and overexpression transgenic lines of PtrHSFB3-1 and PtrMYB092 in P. trichocarpa. The gene-edited trees (mimicking the repressed PtrHSFB3-1 and PtrMYB092 under tension stress) have stem wood composition resembling that of TW during normal growth and under tension stress (i.e., low lignin and high cellulose), whereas the overexpressors showed an opposite effect (high lignin and low cellulose). Individual overexpression of the TFs impeded lignin reduction under tension stress and restored high levels of lignin biosynthesis in the TW. This study offers biological insights to further uncover how metabolism, growth, and stress adaptation are coordinately regulated in trees.

MYB Transcription Factor161 Mediates Feedback Regulation of Secondary wall-associated NAC-Domain1 Family Genes for Wood Formation.

Wood formation is a complex process that involves cell differentiation, cell expansion, secondary wall deposition, and programmed cell death. We constructed a four-layer wood formation transcriptional regulatory network (TRN) in Populus trichocarpa (black cottonwood) that has four Secondary wall-associated NAC-Domain1 (PtrSND1) transcription factor (TF) family members as the top-layer regulators. We characterized the function of a MYB (PtrMYB161) TF in this PtrSND1-TRN, using transgenic P trichocarpa cells and whole plants. PtrMYB161 is a third-layer regulator that directly transactivates five wood formation genes. Overexpression of PtrMYB161 in P. trichocarpa (OE-PtrMYB161) led to reduced wood, altered cell type proportions, and inhibited growth. Integrative analysis of wood cell-based chromatin-binding assays with OE-PtrMYB161 transcriptomics revealed a feedback regulation system in the PtrSND1-TRN, where PtrMYB161 represses all four top-layer regulators and one second-layer regulator, PtrMYB021, possibly affecting many downstream TFs in, and likely beyond, the TRN, to generate the observed phenotypic changes. Our data also suggested that the PtrMYB161's repressor function operates through interaction of the base PtrMYB161 target-binding system with gene-silencing cofactors. PtrMYB161 protein does not contain any known negative regulatory domains. CRISPR-based mutants of PtrMYB161 in P. trichocarpa exhibited phenotypes similar to the wild type, suggesting that PtrMYB161's activator functions are redundant among many TFs. Our work demonstrated that PtrMYB161 binds to multiple sets of target genes, a feature that allows it to function as an activator as well as a repressor. The balance of the two functions may be important to the establishment of regulatory homeostasis for normal growth and development.

A novel synthetic-genetic-array-based yeast one-hybrid system for high discovery rate and short processing time.

Eukaryotic gene expression is often tightly regulated by interactions between transcription factors (TFs) and their DNA cis targets. Yeast one-hybrid (Y1H) is one of the most extensively used methods to discover these interactions. We developed a high-throughput meiosis-directed yeast one-hybrid system using the Magic Markers of the synthetic genetic array analysis. The system has a transcription factor-DNA interaction discovery rate twice as high as the conventional diploid-mating approach and a processing time nearly one-tenth of the haploid-transformation method. The system also offers the highest accuracy in identifying TF-DNA interactions that can be authenticated in vivo by chromatin immunoprecipitation. With these unique features, this meiosis-directed Y1H system is particularly suited for constructing novel and comprehensive genome-scale gene regulatory networks for various organisms.

Hierarchical Transcription Factor and Chromatin Binding Network for Wood Formation in Black Cottonwood (Populus trichocarpa).

Wood remains the world's most abundant and renewable resource for timber and pulp and is an alternative to fossil fuels. Understanding the molecular regulation of wood formation can advance the engineering of wood for more efficient material and energy productions. We integrated a black cottonwood (Populus trichocarpa) wood-forming cell system with quantitative transcriptomics and chromatin binding assays to construct a transcriptional regulatory network (TRN) directed by a key transcription factor (TF), PtrSND1-B1 (secondary wall-associated NAC-domain protein). The network consists of four layers of TF-target gene interactions with quantitative regulatory effects, describing the specificity of how the regulation is transduced through these interactions to activate cell wall genes (effector genes) for wood formation. PtrSND1-B1 directs 57 TF-DNA interactions through 17 TFs transregulating 27 effector genes. Of the 57 interactions, 55 are novel. We tested 42 of these 57 interactions in 30 genotypes of transgenic P. trichocarpa and verified that ∼90% of the tested interactions function in vivo. The TRN reveals common transregulatory targets for distinct TFs, leading to the discovery of nine TF protein complexes (dimers and trimers) implicated in regulating the biosynthesis of specific types of lignin. Our work suggests that wood formation may involve regulatory homeostasis determined by combinations of TF-DNA and TF-TF (protein-protein) regulations.

The AREB1 Transcription Factor Influences Histone Acetylation to Regulate Drought Responses and Tolerance in Populus trichocarpa.

Plants develop tolerance to drought by activating genes with altered levels of epigenetic modifications. Specific transcription factors are involved in this activation, but the molecular connections within the regulatory system are unclear. Here, we analyzed genome-wide acetylated lysine residue 9 of histone H3 (H3K9ac) enrichment and examined its association with transcriptomes in Populus trichocarpa under drought stress. We revealed that abscisic acid-Responsive Element (ABRE) motifs in promoters of the drought-responsive genes PtrNAC006, PtrNAC007, and PtrNAC120 are involved in H3K9ac enhancement and activation of these genes. Overexpressing these PtrNAC genes in P trichocarpa resulted in strong drought-tolerance phenotypes. We showed that the ABRE binding protein PtrAREB1-2 binds to ABRE motifs associated with these PtrNAC genes and recruits the histone acetyltransferase unit ADA2b-GCN5, forming AREB1-ADA2b-GCN5 ternary protein complexes. Moreover, this recruitment enables GCN5-mediated histone acetylation to enhance H3K9ac and enrich RNA polymerase II specifically at these PtrNAC genes for the development of drought tolerance. CRISPR editing or RNA interference-mediated downregulation of any of the ternary members results in highly drought-sensitive P trichocarpa Thus, the combinatorial function of the ternary proteins establishes a coordinated histone acetylation and transcription factor-mediated gene activation for drought response and tolerance in Populus species.

CAD1 and CCR2 protein complex formation in monolignol biosynthesis in Populus trichocarpa.

Lignin is the major phenolic polymer in plant secondary cell walls and is polymerized from monomeric subunits, the monolignols. Eleven enzyme families are implicated in monolignol biosynthesis. Here, we studied the functions of members of the cinnamyl alcohol dehydrogenase (CAD) and cinnamoyl-CoA reductase (CCR) families in wood formation in Populus trichocarpa, including the regulatory effects of their transcripts and protein activities on monolignol biosynthesis. Enzyme activity assays from stem-differentiating xylem (SDX) proteins showed that RNAi suppression of PtrCAD1 in P. trichocarpa transgenics caused a reduction in SDX CCR activity. RNAi suppression of PtrCCR2, the only CCR member highly expressed in SDX, caused a reciprocal reduction in SDX protein CAD activities. The enzyme assays of mixed and coexpressed recombinant proteins supported physical interactions between PtrCAD1 and PtrCCR2. Biomolecular fluorescence complementation and pull-down/co-immunoprecipitation experiments supported a hypothesis of PtrCAD1/PtrCCR2 heterodimer formation. These results provide evidence for the formation of PtrCAD1/PtrCCR2 protein complexes in monolignol biosynthesis in planta.

Improving wood properties for wood utilization through multi-omics integration in lignin biosynthesis.

A multi-omics quantitative integrative analysis of lignin biosynthesis can advance the strategic engineering of wood for timber, pulp, and biofuels. Lignin is polymerized from three monomers (monolignols) produced by a grid-like pathway. The pathway in wood formation of Populus trichocarpa has at least 21 genes, encoding enzymes that mediate 37 reactions on 24 metabolites, leading to lignin and affecting wood properties. We perturb these 21 pathway genes and integrate transcriptomic, proteomic, fluxomic and phenomic data from 221 lines selected from ~2000 transgenics (6-month-old). The integrative analysis estimates how changing expression of pathway gene or gene combination affects protein abundance, metabolic-flux, metabolite concentrations, and 25 wood traits, including lignin, tree-growth, density, strength, and saccharification. The analysis then predicts improvements in any of these 25 traits individually or in combinations, through engineering expression of specific monolignol genes. The analysis may lead to greater understanding of other pathways for improved growth and adaptation.

Reciprocal cross-regulation of VND and SND multigene TF families for wood formation in Populus trichocarpa.

Secondary cell wall (SCW) biosynthesis is the biological process that generates wood, an important renewable feedstock for materials and energy. NAC domain transcription factors, particularly Vascular-Related NAC-Domain (VND) and Secondary Wall-Associated NAC Domain (SND) proteins, are known to regulate SCW differentiation. The regulation of VND and SND is important to maintain homeostasis for plants to avoid abnormal growth and development. We previously identified a splice variant, PtrSND1-A2IR , derived from PtrSND1-A2 as a dominant-negative regulator, which suppresses the transactivation of all PtrSND1 family members. PtrSND1-A2IR also suppresses the self-activation of the PtrSND1 family members except for its cognate transcription factor, PtrSND1-A2, suggesting the existence of an unknown factor needed to regulate PtrSND1-A2 Here, a splice variant, PtrVND6-C1IR , derived from PtrVND6-C1 was discovered that suppresses the protein functions of all PtrVND6 family members. PtrVND6-C1IR also suppresses the expression of all PtrSND1 members, including PtrSND1-A2, demonstrating that PtrVND6-C1IR is the previously unidentified regulator of PtrSND1-A2 We also found that PtrVND6-C1IR cannot suppress the expression of its cognate transcription factor, PtrVND6-C1PtrVND6-C1 is suppressed by PtrSND1-A2IR Both PtrVND6-C1IR and PtrSND1-A2IR cannot suppress their cognate transcription factors but can suppress all members of the other family. The results indicate that the splice variants from the PtrVND6 and PtrSND1 family may exert reciprocal cross-regulation for complete transcriptional regulation of these two families in wood formation. This reciprocal cross-regulation between families suggests a general mechanism among NAC domain proteins and likely other transcription factors, where intron-retained splice variants provide an additional level of regulation.