NAC-MYB-based transcriptional regulation of secondary cell wall biosynthesis in land plants

Tracing post-domestication historical events and screening pre-breeding germplasm from large gene pools in wheat in the absence of phenotype data

Wheat, particularly common wheat (Triticum aestivum L.), is a major crop accounting for 25% of the world cereal production and thriving in diverse ecogeographic regions. Its adaptation to diverse environments arises from its three distinct genomes adapted to different environments and post-domestication anthropogenic interventions. In search of key genomic regions revealing historic events and breeding significance to common wheat, we performed genome scan and genome-environment association (GEA) analyses using high-marker density genotype datasets. Whole-genome scans revealed highly differentiated regions on chromosomes 2A, 3B, and 4A. In-depth analyses corroborated our previous prediction of the 4A differentiated region signifying the separation between Spelt/Macha and other wheat types. Individual chromosome scans captured key introgressions, including one from T. timopheevii and one from Thinopyrum ponticum on 2B and 3D, respectively, as well as known genes such as Vrn-A1 on 5A. GEA highlighted loci linked to latitude-induced environmental variations, influencing traits such as photoperiodism and responses to abiotic stress. Variation at the Vrn-A1 locus on 5A assigned accessions to two haplotypes (6% and 94%). Further analysis on Vrn-A1 coding gene revealed four subgroups of the major haplotype, while the minor haplotype remained undifferentiated. Analyses at differentiated loci mostly dichotomized the population, illustrating the possibility of isolating pre-breeding materials with desirable traits from large gene pools in the absence of phenotype data. Given the current availability of broad genetic data, the genome-scan-GEA hybrid can be an efficient and cost-effective approach for pinpointing environmentally resilient pre-breeding germplasm from vast gene pools, including gene banks regardless of trait characterization.

Active protein ubiquitination regulates xylem vessel functionality

Xylem vessels function in the long-distance conduction of water in land plants. The NAC transcription factor VASCULAR-RELATED NAC-DOMAIN7 (VND7) is a master regulator of xylem vessel cell differentiation in Arabidopsis (Arabidopsis thaliana). We previously isolated suppressor of ectopic xylem vessel cell differentiation induced by VND7 (seiv) mutants. Here, we report that the responsible genes for seiv3, seiv4, seiv6, and seiv9 are protein ubiquitination-related genes encoding PLANT U-BOX46 (PUB46), an uncharacterized F-BOX protein (FBX), PUB36, and UBIQUITIN-SPECIFIC PROTEASE1 (UBP1), respectively. We also found decreased expression of genes downstream of VND7 and abnormal xylem transport activity in the seiv mutants. Upon VND7 induction, ubiquitination levels from 492 and 180 protein groups were upregulated and downregulated, respectively. VND7 induction resulted in the ubiquitination of proteins for cell wall biosynthesis and protein transport, whereas such active protein ubiquitination did not occur in the seiv mutants. We detected the ubiquitination of three lysine residues in VND7: K94, K105, and K260. Substituting K94 with arginine significantly decreased the transactivation activity of VND7, suggesting that the ubiquitination of K94 is crucial for regulating VND7 activity. Our findings highlight the crucial roles of target protein ubiquitination in regulating xylem vessel activity.

Mining Candidate Genes for Leaf Angle in Brassica napus L. by Combining QTL Mapping and RNA Sequencing Analysis

Leaf angle (LA) is an important trait of plant architecture, and individuals with narrow LA can better capture canopy light under high-density planting, which is beneficial for increasing the overall yield per unit area. To study the genetic basis and molecular regulation mechanism of leaf angle in rapeseed, we carried out a series of experiments. Quantitative trait loci (QTL) mapping was performed using the RIL population, and seven QTLs were identified. Transcriptome analysis showed that the cell wall formation/biogenesis processes and biosynthesis/metabolism of cell wall components were the most enrichment classes. Most differentially expressed genes (DEGs) involved in the synthesis of lignin, xylan, and cellulose showed down-regulated expression in narrow leaf material. Microscopic analysis suggested that the cell size affected by the cell wall in the junction area of the stem and petiole was the main factor in leaf petiole angle (LPA) differences. Combining QTL mapping and RNA sequencing, five promising candidate genes BnaA01G0125600ZS, BnaA01G0135700ZS, BnaA01G0154600ZS, BnaA10G0154200ZS, and BnaC03G0294200ZS were identified in rapeseed, and most of them were involved in cell wall biogenesis and the synthesis/metabolism of cell wall components. The results of QTL, transcriptome analysis, and cytological analysis were highly consistent, collectively revealing that genes related to cell wall function played a crucial role in regulating the LA trait in rapeseed. The study provides further insights into LA traits, and the discovery of new QTLs and candidate genes is highly beneficial for genetic improvement.

A GhLac1-centered transcriptional regulatory cascade mediates cotton resistance to Verticillium dahliae through the lignin biosynthesis pathway

The lignin biosynthesis pathway plays a crucial role in the defense response against V. dahliae in cotton, and it is essential to identify the key regulators in this pathway for disease-resistant breeding. In a previous study, the cotton laccase gene GhLac1 was identified as mediating plant broad-spectrum biotic stress tolerance by manipulating phenylpropanoid metabolism. However, the upstream master regulators and regulatory mechanism of lignin are still largely unknown. This study aims to identify the upstream regulators of GhLac1 and explore the molecular mechanism underlying cotton's disease resistance response to V. dahliae. Through the study, three WRKY, three MYB, and one APETALA2/ETHYLENE RESPONSIVE FACTOR (ERF) TFs were identified as differentially responding to V. dahliae infection in cotton. Among these TFs, GhWRKY30, GhWRKY41, GhMYB42, and GhTINY2 were found to directly bind to the GhLac1 promoter and activate its expression. Transient overexpression of these four TFs in cotton led to increased expression of GhLac1 and other the laccase family members, while knockdown of these TFs resulted in reduced lignin accumulation and increased susceptibility to V. dahliae. Additionally, GhWRKY30 and GhWRKY41 were observed to interact with themselves and with each other, synergistically transactivating the GhLac1 promoter. This study reveals a GhLac1-centered transcriptional regulatory cascade of lignin synthesis that contributes to cotton's defense response by modulating lignin metabolism.

QTL mapping for pod quality and yield traits in snap bean (Phaseolus vulgaris L.)

Pod quality and yield traits in snap bean (Phaseolus vulgaris L.) influence consumer preferences, crop adoption by farmers, and the ability of the product to be commercially competitive locally and globally. The objective of the study was to identify the quantitative trait loci (QTL) for pod quality and yield traits in a snap × dry bean recombinant inbred line (RIL) population. A total of 184 F6 RILs derived from a cross between Vanilla (snap bean) and MCM5001 (dry bean) were grown in three field sites in Kenya and one greenhouse environment in Davis, CA, USA. They were genotyped at 5,951 single nucleotide polymorphisms (SNPs), and composite interval mapping was conducted to identify QTL for 16 pod quality and yield traits, including pod wall fiber, pod string, pod size, and harvest metrics. A combined total of 44 QTL were identified in field and greenhouse trials. The QTL for pod quality were identified on chromosomes Pv01, Pv02, Pv03, Pv04, Pv06, and Pv07, and for pod yield were identified on Pv08. Co-localization of QTL was observed for pod quality and yield traits. Some identified QTL overlapped with previously mapped QTL for pod quality and yield traits, with several others identified as novel. The identified QTL can be used in future marker-assisted selection in snap bean.

Stem lodging Resistance-1 controls stem strength by positively regulating the biosynthesis of cell wall components in Capsicum annuum L

Lodging presents a significant challenge in cultivating high-yield crops with extensive above-ground biomass, yet the molecular mechanisms underlying this phenomenon in the Solanaceae family remain largely unexplored. In this study, we identified a gene, CaSLR1 (Capsicum annuum Stem Lodging Resistance 1), which encodes a MYELOBLASTOSIS (MYB) family transcription factor, from a lodging-affected C. annuum EMS mutant. The suppression of CaSLR1 expression in pepper led to notable stem lodging, reduced thickness of the secondary cell wall, and decreased stem strength. A similar phenotype was observed in tomato with the knockdown of SlMYB61, the orthologous gene to CaSLR1. Further investigations demonstrated that CaNAC6, a gene involved in secondary cell wall (SCW) formation, is co-expressed with CaSLR1 and acts as a positive regulator of its expression, as confirmed through yeast one-hybrid, dual-luciferase reporter assays, and electrophoretic mobility shift assays. These findings elucidate the CaNAC6-CaSLR1 module that contributes to lodging resistance, emphasizing the critical role of CaSLR1 in the lodging resistance regulatory network.

PagMYB128 regulates secondary cell wall formation by direct activation of cell wall biosynthetic genes during wood formation in poplar

The biosynthesis of cellulose, lignin, and hemicelluloses in plant secondary cell walls (SCWs) is regulated by a hierarchical transcriptional regulatory network. This network features orthologous transcription factors shared between poplar and Arabidopsis, highlighting a foundational similarity in their genetic regulation. However, knowledge on the discrepant behavior of the transcriptional-level molecular regulatory mechanisms between poplar and Arabidopsis remains limited. In this study, we investigated the function of PagMYB128 during wood formation and found it had broader impacts on SCW formation compared to its Arabidopsis ortholog, AtMYB103. Transgenic poplar trees overexpressing PagMYB128 exhibited significantly enhanced xylem development, with fiber cells and vessels displaying thicker walls, and an increase in the levels of cellulose, lignin, and hemicelluloses in the wood. In contrast, plants with dominant repression of PagMYB128 demonstrated the opposite phenotypes. RNA sequencing and reverse transcription - quantitative polymerase chain reaction showed that PagMYB128 could activate SCW biosynthetic gene expression, and chromatin immunoprecipitation along with yeast one-hybrid, and effector-reporter assays showed this regulation was direct. Further analysis revealed that PagSND1 (SECONDARY WALL-ASSOCIATED NAC-DOMAIN PROTEIN1) directly regulates PagMYB128 but not cell wall metabolic genes, highlighting the pivotal role of PagMYB128 in the SND1-driven regulatory network for wood development, thereby creating a feedforward loop in SCW biosynthesis.

Mechanistic analysis of the hardening process of the thorns on stems of Bougainvillea glabra "Elizabeth Angus"

Introduction: Bougainvillea glabra "Elizabeth Angus" is a thorny woody vine or shrub. However, the hard thorns are considered a deficiency in its ornamental value. Methods: To find the genes and pathways related to the hardening process of the thorns on the stems of B. glabra, the eukaryotic unreferenced transcriptome sequencing analysis was conducted to explore the 3 stages of the thorn-hardening process. Total RNA was extracted from thorns and stems, and transcriptome libraries were constructed and sequenced using unreferenced Illumina sequencing. Results: Gene function annotation was performed using various databases, resulting in 8937 co-annotated genes. The density distribution of Fragments Per Kilobase of transcript per Million mapped reads (FPKM) depicted the overall gene expression patterns. The study found that stage 2 as the period of highest gene expression activity during the thorns hardening process in B. glabra. Differential expression analysis revealed that during thorn-hardening, 1045 genes up-regulated and 391 genes down-regulated significantly in thorns at stage 2 compared to stage 1 (early stage of thorns formation). Meanwhile, 938 genes up-regulated and 784 genes down-regulated significantly in stems. At stage 3, as thorns became harder, 63 genes exhibited notable expression increase and 98 genes' expression decreased obviously within thorns, and 46 genes up-regulated and 29 genes down-regulated in stems, compared to stage 2. Phenylpropanoid biosynthesis was the key step in the hardening process of the thorns of B. glabra. The formation and hardening of thorns on the stem of B. glabra was a process in which lignin gradually accumulated in the thorns, and several genes were involved in this process. They include PAL (EC:4.3.1.24), CYP73A (EC:1.14.14.91), 4CL (EC:6.2.1.12), CCR (EC:1.2.1.44), CAD (EC:1.1.1.195) and POX (EC:1.11.1.7). Discussion: This transcriptome analysis offers insights into the molecular mechanisms underlying thorns development in this plant species.

Xylem cell size regulation is a key adaptive response to water deficit in Eucalyptus grandis

Future climatic scenarios forecast increasingly frequent droughts that will pose substantial consequences on tree mortality. In light of this, drought-tolerant eucalypts have been propagated; however, the severity of these conditions will invoke adaptive responses, impacting the commercially valuable wood properties. To determine what mechanisms govern the wood anatomical adaptive response, highly controlled drought experiments were conducted in Eucalyptus grandis W. Hill ex Maiden, with the tree physiology and transcriptome closely monitored. In response to water deficit, E. grandis displays an isohydric stomatal response to conserve water and enable stem growth to continue, albeit at a reduced rate. Maintaining gaseous exchange is likely a critical short-term response that drives the formation of hydraulically safer xylem. For instance, the development of significantly smaller fibers and vessels was found to increase cellular density, thereby promoting drought tolerance through improved functional redundancy, as well as implosion and cavitation resistance. The transcriptome was explored to identify the molecular mechanisms responsible for controlling xylem cell size during prolonged water deficit. Downregulation of genes associated with cell wall remodeling and the biosynthesis of cellulose, hemicellulose and pectin appeared to coincide with a reduction in cellular enlargement during drought. Furthermore, transcript levels of NAC and MYB transcription factors, vital for cell wall component biosynthesis, were reduced, while those linked to lignification increased. The upregulation of EgCAD and various peroxidases under water deficit did not correlate with an increased lignin composition. However, with the elevated cellular density, a higher lignin content per xylem cross-sectional area was observed, potentially enhancing hydraulic safety. These results support the requirement for higher density, drought-adapted wood as a long-term adaptive response in E. grandis, which is largely influenced by the isohydric stomatal response coupled with cellular expansion-related molecular processes.

PbARF19-mediated auxin signaling regulates lignification in pear fruit stone cells

The stone cells in pear fruits cause rough flesh and low juice, seriously affecting the taste. Lignin has been demonstrated as the main component of stone cells. Auxin, one of the most important plant hormone, regulates most physiological processes in plants including lignification. However, the concentration effect and regulators of auxin on pear fruits stone cell formation remains unclear. Here, endogenous indole-3-acetic acid (IAA) and stone cells were found to be co-localized in lignified cells by immunofluorescence localization analysis. The exogenous treatment of different concentrations of IAA demonstrated that the application of 200 µM IAA significantly reduced stone cell content, while concentrations greater than 500 µM significantly increased stone cell content. Besides, 31 auxin response factors (ARFs) were identified in pear genome. Putative ARFs were predicted as critical regulators involved in the lignification of pear flesh cells by phylogenetic relationship and expression analysis. Furthermore, the negative regulation of PbARF19 on stone cell formation in pear fruit was demonstrated by overexpression in pear fruitlets and Arabidopsis. These results illustrated that the PbARF19-mediated auxin signal plays a critical role in the lignification of pear stone cell by regulating lignin biosynthetic genes. This study provides theoretical and practical guidance for improving fruit quality in pear production.

Cloning and functional mechanism of the dwarf gene gba affecting stem elongation and cellulose biosynthesis in jute (Corchorus olitorius)

Plant height (PH) is an important factor affecting bast fiber yield in jute. Here, we report the mechanism of dwarfism in the 'Guangbaai' (gba) of jute. The mutant gba had shorter internode length and cell length compared to the standard cultivar 'TaiZi 4' (TZ4). Exogenous GA3 treatment indicated that gba is a GA-insensitive dwarf mutant. Quantitative trait locus (QTL) analysis of three PH-related traits via a high-density genetic linkage map according to re-seq showed that a total of 25 QTLs were identified, including 13 QTLs for PH, with phenotypic variation explained ranging from 2.42 to 74.16%. Notably, the functional mechanism of the candidate gene CoGID1a, the gibberellic acid receptor, of the major locus qPHIL5 was evaluated by transgenic analysis and virus-induced gene silencing. A dwarf phenotype-related single nucleotide mutation in CoGID1a was identified in gba, which was also unique to the dwarf phenotype of gba among 57 cultivars. Cogid1a was unable to interact with the growth-repressor DELLA even in the presence of highly accumulated gibberellins in gba. Differentially expressed genes between transcriptomes of gba and TZ4 after GA3 treatment indicated up-regulation of genes involved in gibberellin and cellulose synthesis in gba. Interestingly, it was found that up-regulation of CoMYB46, a key transcription factor in the secondary cell wall, by the highly accumulated gibberellins in gba promoted the expression of cellulose synthase genes CoCesA4 and CoCesA7. These findings provide valuable insights into fiber development affected by endogenous gibberellin accumulation in plants.

Methylation of microRNA genes and its effect on secondary xylem development of stem in poplar

MicroRNAs (miRNAs) and DNA methylation are both vital regulators of gene expression. DNA methylation can affect the transcription of miRNAs, just like coding genes, through methylating the CpG islands in the gene regions of miRNAs. Although previous studies have shown that DNA methylation and miRNAs can each be involved in the process of wood formation, the relationship between the two has been relatively little studied in plant wood formation. Studies have shown that the second internode (IN2) (from top to bottom) of 3-month-old poplar trees can represent the primary stage of poplar stem development and IN8 can represent the secondary stage. There were also significant differences in DNA methylation patterns and miRNA expression patterns obtained from PS and SS. In this study, we first interactively analyzed methylation and miRNA sequencing data to identify 43 differentially expressed miRNAs regulated by differential methylation from the primary stage and secondary stage, which were found to be involved in multiple biological processes related to wood formation by enrichment analysis. In addition, six miRNA/target gene modules were finally identified as potentially involved in secondary xylem development of poplar stems through degradome sequencing and functional analysis. In conclusion, this study provides important reference information on the mechanism of interaction between different regulatory pathways of wood formation.

CcNAC6 Acts as a Positive Regulator of Secondary Cell Wall Synthesis in Sudan Grass (Sorghum sudanense S.)

The degree of forage lignification is a key factor affecting its digestibility by ruminants such as cattle and sheep. Sudan grass (Sorghum sudanense S.) is a high-quality sorghum forage, and its lignocellulose is mostly stored in the secondary cell wall. However, the secondary cell wall synthesis mechanism of Sudan grass has not yet been studied in depth. To further study the secondary cell wall synthesis mechanism of Sudan grass using established transcriptome data, this study found that CcNAC6, a homologous gene of Arabidopsis AtSND2, is related to the secondary cell wall synthesis of Sudan grass. Accordingly, we constructed a CcNAC6-overexpressing line of Arabidopsis to investigate the function of the CcNAC6 gene in secondary cell wall synthesis. The results showed that the overexpression of the CcNAC6 gene could significantly increase the lignin content of Arabidopsis. Based on subcellular localization analysis, CcNAC6 is found in the nucleus. In addition, yeast two-hybridization screening showed that CcCP1, associated with secondary cell wall synthesis, can interact with CcNAC6. Therefore, the above results indicate that CcNAC6 has a positive regulatory effect on the secondary cell wall synthesis of Sudan grass, and it is speculated that CcNAC6 may be the main regulator of the secondary cell wall synthesis of Sudan grass through its interaction with another regulatory protein, CcCP1. This study provides a theoretical basis and new genetic resources for the creation of new Sudan grass germplasm with a low lignin content.

Morphological Innovation Drives Sperm Release in Bryophytes

Plant movements for survival are nontrivial. Antheridia in the moss Physcomitrium patens (P. patens) use motion to eject sperm in the presence of water. However, the biological and mechanical mechanisms that actuate the process are unknown. Here, the burst of the antheridium of P. patens, triggered by water, results from elastic instability and is determined by an asymmetric change in cell geometry. The tension generated in jacket cell walls of antheridium arises from turgor pressure, and is further promoted when the inner walls of apex burst in hydration, causing water and cellular contents of apex quickly influx into sperm chamber. The outer walls of the jacket cells are strengthened by NAC transcription factor VNS4 and serve as key morphomechanical innovations to store hydrostatic energy in a confined space in P. patens. However, the antheridium in liverwort Marchantia polymorpha (M. polymorpha) adopts a different strategy for sperm release; like jacket cell outer walls of P. patens, the cells surrounding the antheridium of M. polymorpha appear to play a similar role in the storage of energy. Collectively, the work shows that plants have evolved different ingenious devices for sperm discharge and that morphological innovations can differ.

Impacts of in-vitro zebularine treatment on genome-wide DNA methylation and transcriptomic profiles in Salix purpurea L

Renewable energy resources such as biomass are crucial for a sustainable global society. Trees are a major source of lignocellulosic biomass, which can vary in response to different environmental factors owing to epigenetic regulation, such as DNA C-methylation. To investigate the effects of DNA methylation on plant development and wood formation, and its impacts on gene expression, with a focus on secondary cell wall (SCW)-associated genes, Salix purpurea plantlets were cloned from buds derived from a single hybrid tree for both treatment and control conditions. For the treatment condition, buds were exposed to 50 μM zebularine in vitro and a combined strategy of whole-genome bisulfite sequencing (WGBS) and RNA-seq was employed to examine the methylome and transcriptome profiles of different tissues collected at various time points under both conditions. Transcriptomic and methylome data revealed that most of the promoter and gene body demethylation had no marked effects on the expression profiles of genes. Nevertheless, gene expression tended to decrease with the increased methylation levels of genes with highly methylated promoters. Results indicated that demethylation is less evident in centromeric regions and sex chromosomes. Promoters of secondary cell wall-associated genes, such as 4-coumarate-CoA ligase-like and Rac-like GTP-binding protein RHO, were differentially methylated in the secondary xylem samples collected from two-month potted treated plants compared to control samples. Our results provide novel insights into DNA methylation and gene expression landscapes and a basis for investigating the epigenetic regulation of wood formation in S. purpurea as a model plant for bioenergy species.

Multilevel analysis between Physcomitrium patens and Mortierellaceae endophytes explores potential long-standing interaction among land plants and fungi

The model moss species Physcomitrium patens has long been used for studying divergence of land plants spanning from bryophytes to angiosperms. In addition to its phylogenetic relationships, the limited number of differential tissues, and comparable morphology to the earliest embryophytes provide a system to represent basic plant architecture. Based on plant-fungal interactions today, it is hypothesized these kingdoms have a long-standing relationship, predating plant terrestrialization. Mortierellaceae have origins diverging from other land fungi paralleling bryophyte divergence, are related to arbuscular mycorrhizal fungi but are free-living, observed to interact with plants, and can be found in moss microbiomes globally. Due to their parallel origins, we assess here how two Mortierellaceae species, Linnemannia elongata and Benniella erionia, interact with P. patens in coculture. We also assess how Mollicute-related or Burkholderia-related endobacterial symbionts (MRE or BRE) of these fungi impact plant response. Coculture interactions are investigated through high-throughput phenomics, microscopy, RNA-sequencing, differential expression profiling, gene ontology enrichment, and comparisons among 99 other P. patens transcriptomic studies. Here we present new high-throughput approaches for measuring P. patens growth, identify novel expression of over 800 genes that are not expressed on traditional agar media, identify subtle interactions between P. patens and Mortierellaceae, and observe changes to plant-fungal interactions dependent on whether MRE or BRE are present. Our study provides insights into how plants and fungal partners may have interacted based on their communications observed today as well as identifying L. elongata and B. erionia as modern fungal endophytes with P. patens.

Allelic variations of WAK106-E2Fa-DPb1-UGT74E2 module regulate fibre properties in Populus tomentosa

Wood formation, intricately linked to the carbohydrate metabolism pathway, underpins the capacity of trees to produce renewable resources and offer vital ecosystem services. Despite their importance, the genetic regulatory mechanisms governing wood fibre properties in woody plants remain enigmatic. In this study, we identified a pivotal module comprising 158 high-priority core genes implicated in wood formation, drawing upon tissue-specific gene expression profiles from 22 Populus samples. Initially, we conducted a module-based association study in a natural population of 435 Populus tomentosa, pinpointing PtoDPb1 as the key gene contributing to wood formation through the carbohydrate metabolic pathway. Overexpressing PtoDPb1 led to a 52.91% surge in cellulose content, a reduction of 14.34% in fibre length, and an increment of 38.21% in fibre width in transgenic poplar. Moreover, by integrating co-expression patterns, RNA-sequencing analysis, and expression quantitative trait nucleotide (eQTN) mapping, we identified a PtoDPb1-mediated genetic module of PtoWAK106-PtoDPb1-PtoE2Fa-PtoUGT74E2 responsible for fibre properties in Populus. Additionally, we discovered the two PtoDPb1 haplotypes that influenced protein interaction efficiency between PtoE2Fa-PtoDPb1 and PtoDPb1-PtoWAK106, respectively. The transcriptional activation activity of the PtoE2Fa-PtoDPb1 haplotype-1 complex on the promoter of PtoUGT74E2 surpassed that of the PtoE2Fa-PtoDPb1 haplotype-2 complex. Taken together, our findings provide novel insights into the regulatory mechanisms of fibre properties in Populus, orchestrated by PtoDPb1, and offer a practical module for expediting genetic breeding in woody plants via molecular design.

Comparative transcriptome analysis reveals the adaptive mechanisms of halophyte Suaeda dendroides encountering high saline environment

Suaeda dendroides, a succulent euhalophyte of the Chenopodiaceae family, intermittently spread around northern Xinjiang, China, has the ability to grow and develop in saline and alkali environments. The objective of this study was therefore to investigate the underlying molecular mechanisms of S. dendroides response to high salt conditions. 27 sequencing libraries prepared from low salt (200 mM NaCl) and high salt (800 mM NaCl) treated plants at 5 different stages were sequenced using Illumina Hiseq 2000. A total of 133,107 unigenes were obtained, of which 4,758 were DEGs. The number of DEGs in the high salt group (3,189) was more than the low salt treatment group (733) compared with the control. GO and KEGG analysis of the DEGs at different time points of the high salt treatment group showed that the genes related to cell wall biosynthesis and modification, plant hormone signal transduction, ion homeostasis, organic osmolyte accumulation, and reactive oxygen species (ROS) detoxification were significantly expressed, which indicated that these could be the main mechanisms of S. dendroides acclimate to high salt stress. The study provides a new perspective for understanding the molecular mechanisms of halophytes adapting to high salinity. It also provides a basis for future investigations of key salt-responsive genes in S. dendroides.

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

RtNAC055 promotes drought tolerance via a stomatal closure pathway linked to methyl jasmonate/hydrogen peroxide signaling in Reaumuria trigyna

The stomata regulate CO2 uptake and efficient water usage, thereby promoting drought stress tolerance. NAC proteins (NAM, ATAF1/2, and CUC2) participate in plant reactions following drought stress, but the molecular mechanisms underlying NAC-mediated regulation of stomatal movement are unclear. In this study, a novel NAC gene from Reaumuria trigyna, RtNAC055, was found to enhance drought tolerance via a stomatal closure pathway. It was regulated by RtMYC2 and integrated with jasmonic acid signaling and was predominantly expressed in stomata and root. The suppression of RtNAC055 could improve jasmonic acid and H2O2 production and increase the drought tolerance of transgenic R. trigyna callus. Ectopic expression of RtNAC055 in the Arabidopsis atnac055 mutant rescued its drought-sensitive phenotype by decreasing stomatal aperture. Under drought stress, overexpression of RtNAC055 in poplar promoted ROS (H2O2) accumulation in stomata, which accelerated stomatal closure and maintained a high photosynthetic rate. Drought upregulated the expression of PtRbohD/F, PtP5CS2, and PtDREB1.1, as well as antioxidant enzyme activities in heterologous expression poplars. RtNAC055 promoted H2O2 production in guard cells by directly binding to the promoter of RtRbohE, thus regulating stomatal closure. The stress-related genes RtDREB1.1/P5CS1 were directly regulated by RtNAC055. These results indicate that RtNAC055 regulates stomatal closure by maintaining the balance between the antioxidant system and H2O2 level, reducing the transpiration rate and water loss, and improving photosynthetic efficiency and drought resistance.

Shoring up the base: the development and regulation of cortical sclerenchyma in grass nodal roots

Plants depend on the combined action of a shoot-root-soil system to maintain their anchorage to the soil. Mechanical failure of any component of this system results in lodging, a permanent and irreversible inability to maintain vertical orientation. Models of anchorage in grass crops identify the compressive strength of roots near the soil surface as key determinant of resistance to lodging. Indeed, studies of disparate grasses report a ring of thickened, sclerenchyma cells surrounding the root cortex, present only at the base of nodal roots. Here, in the investigation of the development and regulation of this agronomically important trait, we show that development of these cells is uncoupled from the maturation of other secondary cell wall-fortified cells, and that cortical sclerenchyma wall thickening is stimulated by mechanical forces transduced from the shoot to the root. We also show that exogenous application of gibberellic acid stimulates thickening of lignified cell types in the root, including cortical sclerenchyma, but is not sufficient to establish sclerenchyma identity in cortex cells. Leveraging the ability to manipulate cortex development via mechanical stimulus, we show that cortical sclerenchyma development alters root mechanical properties and improves resistance to lodging. We describe transcriptome changes associated with cortical sclerenchyma development under both ambient and mechanically stimulated conditions and identify SECONDARY WALL NAC7 as a putative regulator of mechanically responsive cortex cell wall development at the root base.

PtoMYB031, the R2R3 MYB transcription factor involved in secondary cell wall biosynthesis in poplar

Introduction

The biosynthesis of the secondary cell wall (SCW) is orchestrated by an intricate hierarchical transcriptional regulatory network. This network is initiated by first-layer master switches, SCW-NAC transcription factors, which in turn activate the second-layer master switches MYBs. These switches play a crucial role in regulating xylem specification and differentiation during SCW formation. However, the roles of most MYBs in woody plants are yet to be fully understood.

Methods

In this study, we identified and isolated the R2R3-MYB transcription factor, PtoMYB031, from Populus tomentosa. We explored its expression, mainly in xylem tissues, and its role as a transcriptional repressor in the nucleus. We used overexpression and RNA interference techniques in poplar, along with Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays, to analyze the regulatory effects of PtoMYB031.

Results

Overexpression of PtoMYB031 in poplar significantly reduced lignin, cellulose, and hemicellulose content, and inhibited vascular development in stems, resulting in decreased SCW thickness in xylem tissues. Gene expression analysis showed that structural genes involved in SCW biosynthesis were downregulated in PtoMYB031-OE lines. Conversely, RNA interference of PtoMYB031 increased these compounds. Additionally, PtoMYB031 was found to recruit the repressor PtoZAT11, forming a transcriptional inhibition complex.

Discussion

Our findings provide new insights into how PtoMYB031, through its interaction with PtoZAT11, forms a complex that can suppress the expression of key regulatory genes, PtoWND1A and PtoWND2B, in SCW biosynthesis. This study enhances our understanding of the transcriptional regulation involved in SCW formation in poplar, highlighting the significant role of PtoMYB031.

Transcriptome Profiles Reveals ScDREB10 from Syntrichia caninervis Regulated Phenylpropanoid Biosynthesis and Starch/Sucrose Metabolism to Enhance Plant Stress Tolerance

Desiccation is a kind of extreme form of drought stress and desiccation tolerance (DT) is an ancient trait of plants that allows them to survive tissue water potentials reaching -100 MPa or lower. ScDREB10 is a DREB A-5 transcription factor gene from a DT moss named Syntrichia caninervis, which has strong comprehensive tolerance to osmotic and salt stresses. This study delves further into the molecular mechanism of ScDREB10 stress tolerance based on the transcriptome data of the overexpression of ScDREB10 in Arabidopsis under control, osmotic and salt treatments. The transcriptional analysis of weight gene co-expression network analysis (WGCNA) showed that "phenylpropanoid biosynthesis" and "starch and sucrose metabolism" were key pathways in the network of cyan and yellow modules. Meanwhile, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of differentially expressed genes (DEGs) also showed that "phenylpropanoid biosynthesis" and "starch and sucrose metabolism" pathways demonstrate the highest enrichment in response to osmotic and salt stress, respectively. Quantitative real-time PCR (qRT-PCR) results confirmed that most genes related to phenylpropanoid biosynthesis" and "starch and sucrose metabolism" pathways in overexpressing ScDREB10 Arabidopsis were up-regulated in response to osmotic and salt stresses, respectively. In line with the results, the corresponding lignin, sucrose, and trehalose contents and sucrose phosphate synthase activities were also increased in overexpressing ScDREB10 Arabidopsis under osmotic and salt stress treatments. Additionally, cis-acting promoter element analyses and yeast one-hybrid experiments showed that ScDREB10 was not only able to bind with classical cis-elements, such as DRE and TATCCC (MYBST1), but also bind with unknown element CGTCCA. All of these findings suggest that ScDREB10 may regulate plant stress tolerance by effecting phenylpropanoid biosynthesis, and starch and sucrose metabolism pathways. This research provides insights into the molecular mechanisms underpinning ScDREB10-mediated stress tolerance and contributes to deeply understanding the A-5 DREB regulatory mechanism.

Comparative transcriptome analysis revealed that auxin and cell wall biosynthesis play important roles in the formation of hollow hearts in cucumber

BACKGROUND

Hollow heart is a kind of physiological defect that seriously affects the yield, quality, and economic value of cucumber. However, the formation of hollow hearts may relate to multiple factors in cucumber, and it is necessary to conduct analysis.

RESULTS

In this study, hollow and non-hollow fruits of cucumber K07 were used for comparative transcriptome sequencing and analysis. 253 differentially expressed genes and 139 transcription factors were identified as being associated with the formation of hollow hearts. Hormone (auxin) signaling and cell wall biosynthesis were mainly enriched in GO and KEGG pathways. Expression levels of key genes involved in indole-3-acetic acid biosynthesis in carpel were lower in the hollow fruits than non-hollow fruits, while there was no difference in the flesh. The concentration of indole-3-acetic also showed lower in the carpel than flesh. The biosynthetic pathway and content analysis of the main components of the cell wall found that lignin biosynthesis had obvious regularity with hollow heart, followed by hemicellulose and cellulose. Correlation analysis showed that there may be an interaction between auxin and cell wall biosynthesis, and they collectively participate in the formation of hollow hearts in cucumber. Among the differentially expressed transcription factors, MYB members were the most abundant, followed by NAC, ERF, and bHLH.

CONCLUSIONS

The results and analyses showed that the low content of auxin in the carpel affected the activity of enzymes related to cell wall biosynthesis at the early stage of fruit development, resulting in incomplete development of carpel cells, thus forming a hollow heart in cucumber. Some transcription factors may play regulatory roles in this progress. The results may enrich the theory of the formation of hollow hearts and provide a basis for future research.

The Carbon Flow Shifts from Primary to Secondary Metabolism during Xylem Vessel Cell Differentiation in Arabidopsis thaliana

Xylem vessel cell differentiation is characterized by the deposition of a secondary cell wall (SCW) containing cellulose, hemicellulose and lignin. VASCULAR-RELATED NAC-DOMAIN7 (VND7), a plant-specific NAC (NAM, ATAF1/2, and CUC2) transcription factor, is a master regulator of xylem vessel cell differentiation in Arabidopsis (Arabidopsis thaliana). Previous metabolome analysis using the VND7-inducible system in tobacco BY-2 cells successfully revealed significant quantitative changes in primary metabolites during xylem vessel cell differentiation. However, the flow of primary metabolites is not yet well understood. Here, we performed a metabolomic analysis of VND7-inducible Arabidopsis T87 suspension cells. Capillary electrophoresis-time-of-flight mass spectrometry quantified 57 metabolites, and subsequent data analysis highlighted active changes in the levels of UDP-glucose and phenylalanine, which are building blocks of cellulose and lignin, respectively. In a metabolic flow analysis using stable carbon 13 (13C) isotope, the 13C-labeling ratio specifically increased in 3-phosphoglycerate after 12 h of VND7 induction, followed by an increase in shikimate after 24 h of induction, while the inflow of 13C into lactate from pyruvate was significantly inhibited, indicating an active shift of carbon flow from glycolysis to the shikimate pathway during xylem vessel cell differentiation. In support of this notion, most glycolytic genes involved in the downstream of glyceraldehyde 3-phosphate were downregulated following the induction of xylem vessel cell differentiation, whereas genes for the shikimate pathway and phenylalanine biosynthesis were upregulated. These findings provide evidence for the active shift of carbon flow from primary metabolic pathways to the SCW polymer biosynthetic pathway at specific points during xylem vessel cell differentiation.

NACs, generalist in plant life

Plant-specific NAC proteins constitute a major transcription factor family that is well-known for its roles in plant growth, development, and responses to abiotic and biotic stresses. In recent years, there has been significant progress in understanding the functions of NAC proteins. NAC proteins have a highly conserved DNA-binding domain; however, their functions are diverse. Previous understanding of the structure of NAC transcription factors can be used as the basis for their functional diversity. NAC transcription factors consist of a target-binding domain at the N-terminus and a highly versatile C-terminal domain that interacts with other proteins. A growing body of research on NAC transcription factors helps us comprehend the intricate signalling network and transcriptional reprogramming facilitated by NAC-mediated complexes. However, most studies of NAC proteins have been limited to a single function. Here, we discuss the upstream regulators, regulatory components and targets of NAC in the context of their prospective roles in plant improvement strategies via biotechnology intervention, highlighting the importance of the NAC transcription factor family in plants and the need for further research.

UBC6, a ubiquitin-conjugating enzyme, participates in secondary cell wall thickening in the inflorescence stem of Arabidopsis

Secondary cell wall (SCW) thickening in plant inflorescence stems is a complicated cellular process that is essential for stem strength and biomass. Although Arabidopsis NAC transcription factor (TF) 1 (NST1) regulates the SCW thickening in anther walls, the single T-DNA-insertion mutant (nst1) does not show disrupted SCW thickening in anther endothecium, interfascicular fibers or xylem. To better understand the regulatory mechanism of this process, we generated an ethyl methanesulfonate (EMS)-mutagenized Arabidopsis population with the nst1 background. scd5 (SCW-defective mutant 5) was isolated in a forward genetic screen from the EMS mutant library, which displayed not only less lignin deposition in the interfascicular fiber and xylem than the wild type but also a pendent inflorescence stem. The EMS-induced mutation associated with the scd5 phenotype was found in the 5th exon of At2G46030 that encodes a ubiquitin-conjugating enzyme (UBC6), we thereby renamed the allele nst1 ubc6. Overexpressing UBC6 in nst1 ubc6 rescued the defective SCW, whereas disrupting UBC6 in nst1 by the CRISPR/Cas9 system caused a phenotype similar to that observed in nst1 ubc6. UBC6 was localized to the nucleus and plasma membrane, and possessed E2 ubiquitin-conjugating activity in vitro. MYB7 and MYB32 are considered as transcription repressors in the phenylpropanoid pathway and are involved in NAC TF-related transcriptional regulation in SCW thickening. UBC6 can interact with MYB7 and MYB32 and positively mediate the degradation of MYB7 and MYB32 by the 26S proteasome. Overall, these results indicated the contribution of UBC6 to SCW thickening in Arabidopsis inflorescence stems.

Regulation of secondary growth by poplar BLADE-ON-PETIOLE genes in Arabidopsis

BLADE-ON-PETIOLE (BOP) genes are essential regulators of vegetative and reproductive development in land plants. First characterized in Arabidopsis thaliana (Arabidopsis), members of this clade function as transcriptional co-activators by recruiting TGACG-motif binding (TGA) basic leucine zipper (bZIP) transcription factors. Highly expressed at organ boundaries, these genes are also expressed in vascular tissue and contribute to lignin biosynthesis during secondary growth. How these genes function in trees, which undergo extensive secondary growth to produce wood, remains unclear. Here, we investigate the functional conservation of BOP orthologs in Populus trichocarpa (poplar), a widely-used model for tree development. Within the poplar genome, we identified two BOP-like genes, PtrBPL1 and PtrBPL2, with abundant transcripts in stems. To assess their functions, we used heterologous assays in Arabidopsis plants. The promoters of PtrBPL1 and PtrBPL2, fused with a β-glucuronidase (GUS) reporter gene showed activity at organ boundaries and in secondary xylem and phloem. When introduced into Arabidopsis plants, PtrBPL1 and PtrBPL2 complemented leaf and flower patterning defects in bop1 bop2 mutants. Notably, Arabidopsis plants overexpressing PtrBPL1 and PtrBPL2 showed defects in stem elongation and the lignification of secondary tissues in the hypocotyl and stem. Finally, PtrBPL1 and PtrBPL2 formed complexes with TGA bZIP proteins in yeast. Collectively, our findings suggest that PtrBPL1 and PtrBPL2 are orthologs of Arabidopsis BOP1 and BOP2, potentially contributing to secondary growth regulation in poplar trees. This work provides a foundation for functional studies in trees.

Integrated microRNA and whole-transcriptome sequencing reveals the involvement of small and long non-coding RNAs in the fiber growth of ramie plant

BACKGROUND

MicroRNAs (miRNAs) and long non-coding RNAs (lncRNAs) are the two main types of non-coding RNAs that play crucial roles in plant growth and development. However, their specific roles in the fiber growth of ramie plant (Boehmeria nivea L. Gaud) remain largely unknown.

METHODS

In this study, we performed miRNA and whole-transcriptome sequencing of two stem bark sections exhibiting different fiber growth stages to determine the expression profiles of miRNAs, lncRNAs, and protein-encoding genes.

RESULTS

Among the identified 378 miRNAs and 6,839 lncRNAs, 88 miRNAs and 1,288 lncRNAs exhibited differential expression. Bioinformatics analysis revealed that 29 and 228 differentially expressed protein-encoding genes were targeted by differentially expressed miRNAs and lncRNAs, respectively, constituting eight putative competing endogenous RNA networks. lncR00022274 exhibited downregulated expression in barks with growing fibers. It also had an antisense overlap with the MYB gene, BntWG10016451, whose overexpression drastically increased the xylem fiber number and secondary wall thickness of fibers in the stems of transgenic Arabidopsis, suggesting the potential association of lncR00022274-BntWG10016451 expression with fiber growth.

CONCLUSIONS

These findings provide insights into the roles of ncRNAs in the regulation of fiber growth in ramie, which can be used for the biotechnological improvement of its fiber yield and quality in the future.

Comprehensive Time-Course Transcriptome Reveals the Crucial Biological Pathways Involved in the Seasonal Branch Growth in Siberian Elm (Ulmus pumila)

Timber, the most prevalent organic material on this planet, is the result of a secondary xylem emerging from vascular cambium. Yet, the intricate processes governing its seasonal generation are largely a mystery. To better understand the cyclic growth of vascular tissues in elm, we undertook an extensive study examining the anatomy, physiology, and genetic expressions in Ulmus pumila. We chose three robust 15-year-old elm trees for our study. The cultivars used in this study were collected from the Inner Mongolia Autonomous Region in China and nurtured in the tree farm of Shandong Normal University. Monthly samples of 2-year-old elm branches were taken from the tree from February to September. Marked seasonal shifts in elm branch vascular tissues were observed by phenotypic observation: In February, the cambium of the branch emerged from dormancy, spurring growth. By May, elms began generating secondary xylem, or latewood, recognized by its tiny pores and dense cell structure. From June to August, there was a marked increase in the thickness of the secondary xylem. Transcriptome sequencing provides a potential molecular mechanism for the thickening of elm branches and their response to stress. In February, the tree enhanced its genetic responses to cold and drought stress. The amplified expression of CDKB, CYCB, WOX4, and ARF5 in the months of February and March reinforced their essential role in the development of the vascular cambium in elm. Starting in May, the elm deployed carbohydrates as a carbon resource to synthesize the abundant cellulose and lignin necessary for the formation of the secondary wall. Major genes participating in cellulose (SUC and CESA homologs), xylan (UGD, UXS, IRX9, IRX10, and IRX14), and lignin (PAL, C4H, 4CL, HCT, C3H, COMT, and CAD) biosynthetic pathways for secondary wall formation were up-regulated by May or/and June. In conclusion, our findings provided a foundation for an in-depth exploration of the molecular processes dictating the seasonal growth of elm timber.

Functional investigation of five R2R3-MYB transcription factors associated with wood development in Eucalyptus using DAP-seq-ML

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta.

Comparative Root Transcriptome Profiling and Gene Regulatory Network Analysis between Eastern and Western Carrot (Daucus carota L.) Cultivars Reveals Candidate Genes for Vascular Tissue Patterning

Carrot (Daucus carota L.) is a highly consumed vegetable rich in carotenoids, known for their potent antioxidant, anti-inflammatory, and immune-protecting properties. While genetic and molecular studies have largely focused on wild and Western carrot cultivars (cvs), little is known about the evolutionary interactions between closely related Eastern and Western cvs. In this study, we conducted comparative transcriptome profiling of root tissues from Eastern (UHSBC-23-1) and Western (UHSBC-100) carrot cv. to better understand differentially expressed genes (DEGs) associated with storage root development and vascular cambium (VC) tissue patterning. Through reference-guided TopHat mapping, we achieved an average mapping rate of 73.87% and identified a total of 3544 DEGs (p < 0.05). Functional annotation and gene ontology classification revealed 97 functional categories, including 33 biological processes, 19 cellular components, 45 metabolic processes, and 26 KEGG pathways. Notably, Eastern cv. exhibited enrichment in cell wall, plant-pathogen interaction, and signal transduction terms, while Western cv. showed dominance in photosynthesis, metabolic process, and carbon metabolism terms. Moreover, constructed gene regulatory network (GRN) for both cvs. obtained orthologs with 1222 VC-responsive genes of Arabidopsis thaliana. In Western cv, GRN revealed VC-responsive gene clusters primarily associated with photosynthetic processes and carbon metabolism. In contrast, Eastern cv. exhibited a higher number of stress-responsive genes, and transcription factors (e.g., MYB15, WRKY46, AP2/ERF TF connected via signaling pathways with NAC036) were identified as master regulators of xylem vessel differentiation and secondary cell wall thickening. By elucidating the comparative transcriptome profiles of Eastern and Western cvs. for the first time, our study provides valuable insights into the differentially expressed genes involved in root development and VC tissue patterning. The identification of key regulatory genes and their roles in these processes represents a significant advancement in our understanding of the evolutionary relations and molecular mechanisms underlying secondary growth of carrot and regulation by vascular cambium.

Identification of a novel viral factor inducing tumorous symptoms by disturbing vascular development in planta

Plant viruses induce various disease symptoms that substantially impact agriculture, but the underlying mechanisms of viral disease in plants are poorly understood. Kobu-sho is a disease in gentian that shows gall formation with ectopic development of lignified cells and vascular tissues such as xylem. Here, we show that a gene fragment of gentian Kobu-sho-associated virus, which is designated as Kobu-sho-inducing factor (KOBU), induces gall formation accompanied by ectopic development of lignified cells and xylem-like tissue in Nicotiana benthamiana. Transgenic gentian expressing KOBU exhibited tumorous symptoms, confirming the gall-forming activity of KOBU. Surprisingly, KOBU expression can also induce differentiation of an additional leaf-like tissue on the abaxial side of veins in normal N. benthamiana and gentian leaves. Transcriptome analysis with Arabidopsis thaliana expressing KOBU revealed that KOBU activates signaling pathways that regulate xylem development. KOBU protein forms granules and plate-like structures and co-localizes with mRNA splicing factors within the nucleus. Our findings suggest that KOBU is a novel pleiotropic virulence factor that stimulates vascular and leaf development. IMPORTANCE While various mechanisms determine disease symptoms in plants depending on virus-host combinations, the details of how plant viruses induce symptoms remain largely unknown in most plant species. Kobu-sho is a disease in gentian that shows gall formation with ectopic development of lignified cells and vascular tissues such as xylem. Our findings demonstrate that a gene fragment of gentian Kobu-sho-associated virus (GKaV), which is designated as Kobu-sho-inducing factor, induces the gall formation accompanied by the ectopic development of lignified cells and xylem-like tissue in Nicotiana benthamiana. The molecular mechanism by which gentian Kobu-sho-associated virus induces the Kobu-sho symptoms will provide new insight into not only plant-virus interactions but also the regulatory mechanisms underlying vascular and leaf development.

Transcriptome features of stone cell development in weevil-resistant and susceptible Sitka spruce

Stone cells are a specialized, highly lignified cell type found in both angiosperms and gymnosperms. In conifers, abundance of stone cells in the cortex provides a robust constitutive physical defense against stem feeding insects. Stone cells are a major insect-resistance trait in Sitka spruce (Picea sitchensis), occurring in dense clusters in apical shoots of trees resistant (R) to spruce weevil (Pissodes strobi) but being rare in susceptible (S) trees. To learn more about molecular mechanisms of stone cell formation in conifers, we used laser microdissection and RNA sequencing to develop cell-type-specific transcriptomes of developing stone cells from R and S trees. Using light, immunohistochemical, and fluorescence microscopy, we also visualized the deposition of cellulose, xylan, and lignin associated with stone cell development. A total of 1293 genes were differentially expressed at higher levels in developing stone cells relative to cortical parenchyma. Genes with potential roles in stone cell secondary cell wall formation (SCW) were identified and their expression evaluated over a time course of stone cell formation in R and S trees. The expression of several transcriptional regulators was associated with stone cell formation, including a NAC family transcription factor and several genes annotated as MYB transcription factors with known roles in SCW formation.

NAC domain transcription factors VNI2 and ATAF2 form protein complexes and regulate leaf senescence

The NAM, ATAF1/2, and CUC2 (NAC) domain transcription factor VND-INTERACTING2 (VNI2) negatively regulates xylem vessel formation by interacting with another NAC domain transcription factor, VASCULAR-RELATED NAC-DOMAIN7 (VND7), a master regulator of xylem vessel formation. Here, we screened interacting proteins with VNI2 using yeast two-hybrid assay and isolated two NAC domain transcription factors, Arabidopsis thaliana ACTIVATION FACTOR 2 (ATAF2) and NAC DOMAIN CONTAINING PROTEIN 102 (ANAC102). A transient gene expression assay showed that ATAF2 upregulates the expression of genes involved in leaf senescence, and VNI2 effectively inhibits the transcriptional activation activity of ATAF2. vni2 mutants accelerate leaf senescence, whereas ataf2 mutants delay leaf senescence. In addition, the accelerated leaf senescence phenotype of the vni2 mutant is recovered by simultaneous mutation of ATAF2. Our findings strongly suggest that VNI2 interacts with and inhibits ATAF2, resulting in negatively regulating leaf senescence.

Association Mapping and Expression Analysis of the Genes Involved in the Wood Formation of Poplar

Xylogenesis is a complex and sequential biosynthetic process controlled by polygenes. Deciphering the genetic architecture of this complex quantitative trait could provide valuable information for increasing wood biomass and improving its properties. Here, we performed genomic resequencing of 64 24-year-old trees (64 hybrids of section Aigeiros and their parents) grown in the same field and conducted full-sib family-based association analyses of two growth and six woody traits using GEMMA as a choice of association model selection. We identified 1342 significantly associated single nucleotide polymorphisms (SNPs), 673 located in the region upstream and downstream of 565 protein-encoding genes. The transcriptional regulation network of secondary cell wall (SCW) biosynthesis was further constructed based on the published data of poplar miRNA, transcriptome, and degradome. These provided a certain scientific basis for the in-depth understanding of the mechanism of poplar timber formation and the molecular-assisted breeding in the future.

Genome-Wide Identification, Characterization, and Expression Analysis of the NAC Gene Family in Litchi chinensis

NAC proteins play an essential role in the growth and development of litchi, especially during reproductive development. However, a comprehensive analysis of the litchi NAC gene family is currently absent. Based on information from the litchi genome, we found that the 112 NAC genes of litchi show an uneven distribution on the chromosomes. Phylogenetic and conserved structural domain analyses indicated that different types of variability were exhibited in the family of litchi NACs (LcNACs). Gene covariance analysis showed that the LcNACs showed better similarity in the same genus than with Arabidopsis. We further investigated the differential expression patterns of LcNACs in buds and rudimentary leaves of litchi. qRT-PCR results implied that they were involved in the process. Profiling of LcNAC promoter elements in litchi showed that they were extensively involved in light response, phytohormone regulation, abiotic stress response, and plant growth and development processes. This study provides new insights into the identification, structural characterization, tissue-specific expression analysis, and promoter response elements of LcNACs. It reveals the characteristics of the LcNACs and lays the foundation for the subsequent understanding of its biological functions and molecular regulatory mechanisms.

The OsbHLH002/OsICE1-OSH1 module orchestrates secondary cell wall formation in rice

Transcriptional regulation of secondary cell wall (SCW) formation is strictly controlled by a complex network of transcription factors in vascular plants and has been shown to be mediated by a group of NAC master switches. In this study, we show that in a bHLH transcription factor, OsbHLH002/OsICE1, its loss-of-function mutant displays a lodging phenotype. Further results show that OsbHLH002 and Oryza sativa homeobox1 (OSH1) interact and share a set of common targets. In addition, the DELLA protein SLENDER RICE1, rice ortholog of KNOTTED ARABIDOPSIS THALIANA7, and OsNAC31 interact with OsbHLH002 and OSH1 and regulate their binding capacity on OsMYB61, a key regulatory factor in SCW development. Collectively, our results indicate OsbHLH002 and OSH1 as key regulators in SCW formation and shed light on molecular mechanisms of how active and repressive factors precisely orchestrate SCW synthesis in rice, which may provide a strategy for manipulating plant biomass production.

Transcriptomic and Metabolomic Analyses of the Response of Resistant Peanut Seeds to Aspergillus flavus Infection

Peanut seeds are susceptible to Aspergillus flavus infection, which has a severe impact on the peanut industry and human health. However, the molecular mechanism underlying this defense remains poorly understood. The aim of this study was to analyze the changes in differentially expressed genes (DEGs) and differential metabolites during A. flavus infection between Zhonghua 6 and Yuanza 9102 by transcriptomic and metabolomic analysis. A total of 5768 DEGs were detected in the transcriptomic study. Further functional analysis showed that some DEGs were significantly enriched in pectinase catabolism, hydrogen peroxide decomposition and cell wall tissues of resistant varieties at the early stage of infection, while these genes were differentially enriched in the middle and late stages of infection in the nonresponsive variety Yuanza 9102. Some DEGs, such as those encoding transcription factors, disease course-related proteins, peroxidase (POD), chitinase and phenylalanine ammonialyase (PAL), were highly expressed in the infection stage. Metabolomic analysis yielded 349 differential metabolites. Resveratrol, cinnamic acid, coumaric acid, ferulic acid in phenylalanine metabolism and 13S-HPODE in the linolenic acid metabolism pathway play major and active roles in peanut resistance to A. flavus. Combined analysis of the differential metabolites and DEGs showed that they were mainly enriched in phenylpropane metabolism and the linolenic acid metabolism pathway. Transcriptomic and metabolomic analyses further confirmed that peanuts infected with A. flavus activates various defense mechanisms, and the response to A. flavus is more rapid in resistant materials. These results can be used to further elucidate the molecular mechanism of peanut resistance to A. flavus infection and provide directions for early detection of infection and for breeding peanut varieties resistant to aflatoxin contamination.

Comparative Transcriptome Analysis Reveals the Effect of Lignin on Storage Roots Formation in Two Sweetpotato (Ipomoea batatas (L.) Lam.) Cultivars

Sweet potato (Ipomoea batatas (L.) Lam.) is one of the most important crops with high storage roots yield. The formation and expansion rate of storage root (SR) plays a crucial role in the production of sweet potato. Lignin affects the SR formation; however, the molecular mechanisms of lignin in SR development have been lacking. To reveal the problem, we performed transcriptome sequencing of SR harvested at 32, 46, and 67 days after planting (DAP) to analyze two sweet potato lines, Jishu25 and Jishu29, in which SR expansion of Jishu29 was early and had a higher yield. A total of 52,137 transcripts and 21,148 unigenes were obtained after corrected with Hiseq2500 sequencing. Through the comparative analysis, 9577 unigenes were found to be differently expressed in the different stages in two cultivars. In addition, phenotypic analysis of two cultivars, combined with analysis of GO, KEGG, and WGCNA showed the regulation of lignin synthesis and related transcription factors play a crucial role in the early expansion of SR. The four key genes swbp1, swpa7, IbERF061, and IbERF109 were proved as potential candidates for regulating lignin synthesis and SR expansion in sweet potato. The data from this study provides new insights into the molecular mechanisms underlying the impact of lignin synthesis on the formation and expansion of SR in sweet potatoes and proposes several candidate genes that may affect sweet potato yield.

PagUNE12 encodes a basic helix-loop-helix transcription factor that regulates the development of secondary vascular tissue in poplar

Secondary growth in woody plants generates new cells and tissues via the activity of the vascular cambium and drives the radial expansion of stems and roots. It is regulated by a series of endogenous factors, especially transcription factors. Here, we cloned the basic helix-loop-helix (bHLH) transcription factor gene UNFERTILIZED EMBRYO SAC12 (UNE12) from poplar (Populus alba × Populus glandulosa Uyeki) and used biochemical, molecular, and cytological assays to investigate the biological functions and regulatory mechanism of PagUNE12. PagUNE12 mainly localized in the nucleus and possessed transcriptional activation activity. It was widely expressed in vascular tissues, including primary phloem and xylem and secondary phloem and xylem. Poplar plants overexpressing PagUNE12 showed significantly reduced plant height, shorter internodes, and curled leaves compared with wild-type plants. Optical microscopy and transmission electron microscopy revealed that overexpressing PagUNE12 promoted secondary xylem development, with thicker secondary cell walls than wild-type poplar. Fourier transform infrared spectroscopy, confocal Raman microscopy, and 2D Heteronuclear Single Quantum Correlation analysis indicated that these plants also had increased lignin contents, with a lower relative abundance of syringyl lignin units and a higher relative abundance of guaiacyl lignin units. Therefore, overexpressing PagUNE12 promoted secondary xylem development and increased the lignin contents of secondary xylem in poplar, suggesting that this gene could be used to improve wood quality in the future.

The MADS-box gene EjAGL15 positively regulates lignin deposition in the flesh of loquat fruit during its storage

Introduction

Lignification of fruit flesh is a common physiological disorder that occurs during post-harvest storage, resulting in the deterioration of fruit quality. Lignin deposition in loquat fruit flesh occurs due to chilling injury or senescence, at temperatures around 0°C or 20°C, respectively. Despite extensive research on the molecular mechanisms underlying chilling-induced lignification, the key genes responsible for the lignification process during senescence in loquat fruit remain unknown. MADS-box genes, an evolutionarily conserved transcription factor family, have been suggested to play a role in regulating senescence. However, it is still unclear whether MADS-box genes can regulate the lignin deposition that arises from fruit senescence.

Methods

Both senescence- and chilling-induced flesh lignification were simulated by applying temperature treatments on loquat fruits. The flesh lignin content during the storage was measured. Transcriptomic, quantitative reverse transcription PCR and correlation analysis were employed to identify key MADS-box genes that may be involved in flesh lignification. The Dual-luciferase assay was utilized to identify the potential interactions between MADS-box members and genes in phenylpropanoid pathway.

Results and Discussion

The lignin content of the flesh samples treated at 20°C or 0°C increased during storage, but at different rates. Results from transcriptome analysis, quantitative reverse transcription PCR, and correlation analysis led us to identify a senescence-specific MADS-box gene, EjAGL15, which correlated positively with the variation in lignin content of loquat fruit. Luciferase assay results confirmed that EjAGL15 activated multiple lignin biosynthesis-related genes. Our findings suggest that EjAGL15 functions as a positive regulator of senescence-induced flesh lignification in loquat fruit.

Genomic characterization of the NAC transcription factors, directed at understanding their functions involved in endocarp lignification of iron walnut (Juglans sigillata Dode)

The NAC (NAM, ATAF1/2, and CUC2) transcription factors (TF), one of the largest plant-specific gene families, play important roles in the regulation of plant growth and development, stress response and disease resistance. In particular, several NAC TFs have been identified as master regulators of secondary cell wall (SCW) biosynthesis. Iron walnut (Juglans sigillata Dode), an economically important nut and oilseed tree, has been widely planted in the southwest China. The thick and high lignified shell derived endocarp tissues, however, brings troubles in processing processes of products in industry. It is indispensable to dissect the molecular mechanism of thick endocarp formation for further genetic improvement of iron walnut. In the present study, based on genome reference of iron walnut, 117 NAC genes, in total, were identified and characterized in silico, which involves only computational analysis to provide insight into gene function and regulation. We found that the amino acids encoded by these NAC genes varied from 103 to 1,264 in length, and conserved motif numbers ranged from 2 to 10. The JsiNAC genes were unevenly distributed across the genome of 16 chromosomes, and 96 of these genes were identified as segmental duplication genes. Furthermore, 117 JsiNAC genes were divided into 14 subfamilies (A-N) according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana and common walnut (Juglans regia). Furthermore, tissue-specific expression pattern analysis demonstrated that a majority of NAC genes were constitutively expressed in five different tissues (bud, root, fruit, endocarp, and stem xylem), while a total of 19 genes were specifically expressed in endocarp, and most of them also showed high and specific expression levels in the middle and late stages during iron walnut endocarp development. Our result provided a new insight into the gene structure and function of JsiNACs in iron walnut, and identified key candidate JsiNAC genes involved in endocarp development, probably providing mechanistic insight into shell thickness formation across nut species.

PagERF81 regulates lignin biosynthesis and xylem cell differentiation in poplar

Lignin is a major component of plant cell walls and is essential for plant growth and development. Lignin biosynthesis is controlled by a hierarchical regulatory network involving multiple transcription factors. In this study, we showed that the gene encoding an APETALA 2/ethylene-responsive element binding factor (AP2/ERF) transcription factor, PagERF81, from poplar 84 K (Populus alba × P. glandulosa) is highly expressed in expanding secondary xylem cells. Two independent homozygous Pagerf81 mutant lines created by gene editing, produced significantly more but smaller vessel cells and longer fiber cells with more lignin in cell walls, while PagERF81 overexpression lines had less lignin, compared to non-transgenic controls. Transcriptome and reverse transcription quantitative PCR data revealed that multiple lignin biosynthesis genes including Cinnamoyl CoA reductase 1 (PagCCR1), Cinnamyl alcohol dehydrogenase 6 (PagCAD6), and 4-Coumarate-CoA ligase-like 9 (Pag4CLL9) were up-regulated in Pagerf81 mutants, but down-regulated in PagERF81 overexpression lines. In addition, a transient transactivation assay revealed that PagERF81 repressed the transcription of these three genes. Furthermore, yeast one hybrid and electrophoretic mobility shift assays showed that PagERF81 directly bound to a GCC sequence in the PagCCR1 promoter. No known vessel or fiber cell differentiation related genes were differentially expressed, so the smaller vessel cells and longer fiber cells observed in the Pagerf81 lines might be caused by abnormal lignin deposition in the secondary cell walls. This study provides insight into the regulation of lignin biosynthesis, and a molecular tool to engineer wood with high lignin content, which would contribute to the lignin-related chemical industry and carbon sequestration.

OsPEX1, an extensin-like protein, negatively regulates root growth in a gibberellin-mediated manner in rice

Leucine-rich repeat extensins (LRXs) are required for plant growth and development through affecting cell growth and cell wall formation. LRX gene family can be classified into two categories: predominantly vegetative-expressed LRX and reproductive-expressed PEX. In contrast to the tissue specificity of Arabidopsis PEX genes in reproductive organs, rice OsPEX1 is also highly expressed in roots in addition to reproductive tissue. However, whether and how OsPEX1 affects root growth is unclear. Here, we found that overexpression of OsPEX1 retarded root growth by reducing cell elongation likely caused by an increase of lignin deposition, whereas knockdown of OsPEX1 had an opposite effect on root growth, indicating that OsPEX1 negatively regulated root growth in rice. Further investigation uncovered the existence of a feedback loop between OsPEX1 expression level and GA biosynthesis for proper root growth. This was supported by the facts that exogenous GA₃ application downregulated transcript levels of OsPEX1 and lignin-related genes and rescued the root developmental defects of the OsPEX1 overexpression mutant, whereas OsPEX1 overexpression reduced GA level and the expression of GA biosynthesis genes. Moreover, OsPEX1 and GA showed antagonistic action on the lignin biosynthesis in root. OsPEX1 overexpression upregulated transcript levels of lignin-related genes, whereas exogenous GA₃ application downregulated their expression. Taken together, this study reveals a possible molecular pathway of OsPEX1mediated regulation of root growth through coordinate modulation of lignin deposition via a negative feedback regulation between OsPEX1 expression and GA biosynthesis.

Ectopic Expression of FvVND4c Promotes Secondary Cell Wall Thickening and Flavonoid Accumulation in Fragaria vesca

Secondary cell wall (SCW) thickening has a significant effect on the growth and development of plants, as well as in the resistance to various biotic and abiotic stresses. Lignin accounts for the strength of SCW. It is synthesized through the phenylpropanoid pathway that also leads to flavonoid synthesis. The coupling strategies for lignin and flavonoid syntheses are diverse in plants. How their syntheses are balanced by transcriptional regulation in fleshy fruits is still unclear. The diploid strawberry (Fragaria vesca) is a model for fleshy fruits research due to its small genome and wide scope of genetic transformation. SCW thickening is regulated by a multilevel transcriptional regulatory network wherein vascular-related NAC domains (VNDs) act as key regulators. In this study, we systematically characterized VNDs in Fragaria vesca and explored their functions. The overexpression of FvVND4c in diploid strawberry fruits resulted in SCW thickening and fruit color changes accompanied with the accumulation of lignin and flavonoids. Genes related to these phenotypes were also induced upon FvVND4c overexpression. Among the induced genes, we found FvMYB46 to be a direct downstream regulator of FvVND4c. The overexpression of FvMYB46 resulted in similar phenotypes as FvVND4c, except for the color change. Transcriptomic analyses suggest that both FvVND4c and FvMYB46 act on phenylpropanoid and flavonoid biosynthesis pathways, and induce lignin synthesis for SCW. These results suggest that FvVND4c and FvMYB46 cooperatively regulate SCW thickening and flavonoid accumulation in Fragaria vesca.

Alternative Splicing of NAC Transcription Factor Gene CmNST1 Is Associated with Naked Seed Mutation in Pumpkin, Cucurbita moschata

In pumpkin (Cucurbita moschata), the naked or hull-less seed phenotype has great benefits for breeding this crop for oil or snack use. We previously identified a naked seed mutant in this crop. In this study, we report genetic mapping, identification, and characterization of a candidate gene for this mutation. We showed that the naked seed phenotype is controlled by a single recessive gene (N). The bulked segregant analysis identified a 2.4 Mb region on Chromosome 17 with 15 predicted genes. Multiple lines of evidence suggested that CmoCh17G004790 is the most probable candidate gene for the N locus which encodes a NAC transcription factor WALL THICKENING PROMOTING FACTOR 1 (CmNST1). No nucleotide polymorphism or structural variation was found in the genomic DNA sequences of CmNST1 between the mutant and the wildtype inbred line (hulled seed). However, the cDNA sequence cloned from developing seed coat samples of the naked seed mutant was 112 bp shorter than that from the wildtype which is due to seed coat-specific alternative splicing in the second exon of the mutant CmNST1 transcript. The expression level of CmNST1 in the developing seed coat was higher in the mutant than in the wildtype during early seed coat development which was reversed later. Transcriptomic profiling with RNA-Seq at different stages of seed development in the mutant and wildtype revealed a critical role of CmNST1 as a master regulator for the lignin biosynthesis pathway during seed coat development while other NAC and MYB transcription factors were also involved in forming a regulatory network for the building of secondary cell walls. This work provides a novel mechanism for the well-characterized NST1 transcription factor gene in regulating secondary cell wall development. The cloned gene also provides a useful tool for marker-assisted breeding of hull-less C. moschata varieties.

CcNAC1 by Transcriptome Analysis Is Involved in Sudan Grass Secondary Cell Wall Formation as a Positive Regulator

Sudan grass is a high-quality forage of sorghum. The degree of lignification of Sudan grass is the main factor affecting its digestibility in ruminants such as cattle and sheep. Almost all lignocellulose in Sudan grass is stored in the secondary cell wall, but the mechanism and synthesis of the secondary cell wall in Sudan grass is still unclear. In order to study the mechanism of secondary cell wall synthesis in Sudan grass, we used an in vitro induction system of Sudan grass secondary cell wall. Through transcriptome sequencing, it was found that the NAC transcription factor CcNAC1 gene was related to the synthesis of the Sudan grass secondary cell wall. This study further generated CcNAC1 overexpression lines of Arabidopsis to study CcNAC1 gene function in secondary cell wall synthesis. It was shown that the overexpression of the CcNAC1 gene can significantly increase lignin content in Arabidopsis lines. Through subcellular localization analysis, CcNAC1 genes could be expressed in the nucleus of a plant. In addition, we used yeast two-hybrid screening to find 26 proteins interacting with CcNAC1. GO and KEGG analysis showed that CcNAC1 relates to the metabolic pathways and biosynthesis of secondary metabolites. In summary, the synthesis of secondary cell wall of Sudan grass can be regulated by CcNAC1.

Meliaceae genomes provide insights into wood development and limonoids biosynthesis

Meliaceae is a useful plant family owing to its high-quality timber and its many limonoids that have pharmacological and biological activities. Although some genomes of Meliaceae species have been reported, many questions regarding their unique family features, namely wood quality and natural products, have not been answered. In this study, we provide the whole-genome sequence of Melia azedarach comprising 237.16 Mb with a contig N50 of 8.07 Mb, and an improved genome sequence of Azadirachta indica comprising 223.66 Mb with a contig N50 of 8.91 Mb. Moreover, genome skimming data, transcriptomes and other published genomes were comprehensively analysed to determine the genes and proteins that produce superior wood and valuable limonoids. Phylogenetic analysis of chloroplast genomes, single-copy gene families and single-nucleotide polymorphisms revealed that Meliaceae should be classified into two subfamilies: Cedreloideae and Melioideae. Although the Meliaceae species did not undergo additional whole-genome duplication events, the secondary wall biosynthetic genes of the woody Cedreloideae species, Toona sinensis, expanded significantly compared to those of A. indica and M. azedarach, especially in downstream transcription factors and cellulose/hemicellulose biosynthesis-related genes. Moreover, expanded special oxidosqualene cyclase catalogues can help diversify Sapindales skeletons, and the clustered genes that regulate terpene chain elongation, cyclization and modification would support their roles in limonoid biosynthesis. The expanded clans of terpene synthase, O-methyltransferase and cytochrome P450, which are mainly derived from tandem duplication, are responsible for the different limonoid classes among the species. These results are beneficial for further investigations of wood development and limonoid biosynthesis.

Integrated metabolome, transcriptome analysis, and multi-flux full-length sequencing offer novel insights into the function of lignin biosynthesis as a Sesuvium portulacastrum response to salt stress

Sesuvium portulacastrum is a typical halophyte. However, few studies have investigated its salt-tolerant molecular mechanism. In this study, metabolome, transcriptome, and multi-flux full-length sequencing analysis were conducted to investigate the significantly different metabolites (SDMs) and differentially expressed genes (DEGs) of S. portulacastrum samples under salinity. The complete-length transcriptome of S. portulacastrum was developed, which contained 39,659 non-redundant unigenes. RNA-seq results showed that 52 DEGs involved in lignin biosynthesis may be responsible for S. portulacastrum salt tolerance. Furthermore, 130 SDMs were identified, and the salt response could be attributed to the p-coumaryl alcohol-rich in lignin biosynthesis. The co-expression network that was constructed after comparing the different salt treatment processes showed that the p-Coumaryl alcohol was linked to 30 DEGs. Herein, 8 structures genes, i.e., Sp4CL, SpCAD, SpCCR, SpCOMT, SpF5H, SpCYP73A, SpCCoAOMT, and SpC3'H were identified as significant factors in regulating lignin biosynthesis. Further investigation revealed that 64 putative transcription factors (TFs) may interact with the promoters of the above-mentioned genes. Together, the data revealed a potential regulatory network comprising important genes, putative TFs, and metabolites involved in the lignin biosynthesis of S. portulacastrum roots under salt stress, which could serve as a rich useful genetic resource for breeding excellent salt-tolerant plants.