MYB58 and MYB63 are transcriptional activators of the lignin biosynthetic pathway during secondary cell wall formation in Arabidopsis

Allele-specific expression and chromatin accessibility contribute to heterosis in tea plants (Camellia sinensis)

Heterosis is extensively used to improve crop productivity, yet its allelic and chromatin regulation remains unclear. Based on our resolved genomes of the maternal TGY and paternal HD, we analyzed the contribution of allele-specific expression (ASE) and chromatin accessibility of JGY and HGY, the artificial hybrids of oolong tea with the largest cultivated area in China. The ASE genes (ASEGs) of tea hybrids with maternal-biased were mainly related to the energy and terpenoid metabolism pathways, whereas the ASEGs with paternal-biased tend to be enriched in glutathione metabolism, and these parental bias of hybrids may coordinate and lead to the acquisition of heterosis in more biological pathways. ATAC-seq results showed that hybrids have significantly higher accessible chromatin regions (ACRs) compared with their parents, which may confer broader and stronger transcriptional activity of genes in hybrids. The number of ACRs with significantly increased accessibility in hybrids was much greater than decreased, and the associated alleles were also affected by differential ACRs across different parents, suggesting enhanced positive chromatin regulation and potential genetic effects in hybrids. Core ASEGs of terpene and purine alkaloid metabolism pathways with significant positive heterosis have greater chromatin accessibility in hybrids, and were potentially regulated by several members of the MYB, DOF and TRB families. The binding motif of CsMYB85 in the promoter ACR of the rate-limiting enzyme CsDXS was verified by DAP-seq. These results suggest that higher numbers and more accessible ACRs in hybrids contribute to the regulation of ASEGs, thereby affecting the formation of heterotic metabolites.

CC-type glutaredoxin, MeGRXC3, associates with catalases and negatively regulates drought tolerance in cassava (Manihot esculenta Crantz)

Glutaredoxins (GRXs) are essential for reactive oxygen species (ROS) homeostasis in responses of plants to environment changes. We previously identified several drought-responsive CC-type GRXs in cassava, an important tropical crop. However, how CC-type GRX regulates ROS homeostasis of cassava under drought stress remained largely unknown. Here, we report that a drought-responsive CC-type GRX, namely MeGRXC3, was associated with activity of catalase in the leaves of 100 cultivars (or unique unnamed genotypes) of cassava under drought stress. MeGRXC3 negatively regulated drought tolerance by modulating drought- and abscisic acid-induced stomatal closure in transgenic cassava. It antagonistically regulated hydrogen peroxide (H₂ O₂ ) accumulation in epidermal cells and guard cells. Moreover, MeGRXC3 interacted with two catalases of cassava, MeCAT1 and MeCAT2, and regulated their activity in vivo. Additionally, MeGRXC3 interacts with a cassava TGA transcription factor, MeTGA2, in the nucleus, and regulates the expression of MeCAT7 through a MeTGA2-MeMYB63 pathway. Overall, we demonstrated the roles of MeGRXC3 in regulating activity of catalase at both transcriptional and post-translational levels, therefore involving in ROS homeostasis and stomatal movement in responses of cassava to drought stress. Our study provides the first insights into how MeGRXC3 may be used in molecular breeding of cassava crops.

MYB transcription factors-master regulators of phenylpropanoid biosynthesis and diverse developmental and stress responses

Phenylpropanoids, the largest class of natural products including flavonoids, anthocyanins, monolignols and tannins perform multiple functions ranging from photosynthesis, nutrient uptake, regulating growth, cell division, maintenance of redox homeostasis and biotic and abiotic stress responses. Being sedentary life forms, plants possess several regulatory modules that increase their performance in varying environments by facilitating activation of several signaling cascades upon perception of developmental and stress signals. Of the various regulatory modules, those involving MYB transcription factors are one of the extensive groups involved in regulating the phenylpropanoid metabolic enzymes in addition to other genes. R2R3 MYB transcription factors are a class of plant-specific transcription factors that regulate the expression of structural genes involved in anthocyanin, flavonoid and monolignol biosynthesis which are indispensable to several developmental pathways and stress responses. The aim of this review is to present the regulation of the phenylpropanoid pathway by MYB transcription factors via Phospholipase D/phosphatidic acid signaling, downstream activation of the structural genes, leading to developmental and/or stress responses. Specific MYB transcription factors inducing or repressing specific structural genes of anthocyanin, flavonoid and lignin biosynthetic pathways are discussed. Further the roles of MYB in activating biotic and abiotic stress responses are delineated. While several articles have reported the role of MYB's in stress responses, they are restricted to two or three specific MYB factors. This review is a consolidation of the diverse roles of different MYB transcription factors involved both in induction and repression of anthocyanin, flavonoid, and lignin biosynthesis.

Mechanism of Stone (Hardened Endocarp) Formation in Fruits: An Attempt toward Pitless Fruits, and Its Advantages and Disadvantages

Stone (hardened endocarp) has a very important role in the continuity of plant life. Nature has gifted plants with various seed protection and dispersal strategies. Stone-fruit-bearing species have evolved a unique adaptation in which the seed is encased in an extremely hard wood-like shell called the stone. The lignification of the fruit endocarp layer produces the stone, a feature that separates drupes from other plants. Stone cells emerge from parenchyma cells after programmed cell death and the deposition of cellulose and lignin in the secondary cell wall. Generally, the deposition of lignin in primary cell walls is followed by secondary thickening of cell walls to form stone cells. This review article describes the molecular mechanisms and factors that influence the production of stone in the fruit. This is the first review article that describes the molecular mechanisms regulating stone (harden endocarp) formation in fruits. This article will help breeders understand the molecular and genetic basis for the stone formation in fruit, and this could lead to new and innovative directions to breed stoneless fruit cultivars in the future.

MINI-EX: Integrative inference of single-cell gene regulatory networks in plants

Multicellular organisms, such as plants, are characterized by highly specialized and tightly regulated cell populations, establishing specific morphological structures and executing distinct functions. Gene regulatory networks (GRNs) describe condition-specific interactions of transcription factors (TFs) regulating the expression of target genes, underpinning these specific functions. As efficient and validated methods to identify cell-type-specific GRNs from single-cell data in plants are lacking, limiting our understanding of the organization of specific cell types in both model species and crops, we developed MINI-EX (Motif-Informed Network Inference based on single-cell EXpression data), an integrative approach to infer cell-type-specific networks in plants. MINI-EX uses single-cell transcriptomic data to define expression-based networks and integrates TF motif information to filter the inferred regulons, resulting in networks with increased accuracy. Next, regulons are assigned to different cell types, leveraging cell-specific expression, and candidate regulators are prioritized using network centrality measures, functional annotations, and expression specificity. This embedded prioritization strategy offers a unique and efficient means to unravel signaling cascades in specific cell types controlling a biological process of interest. We demonstrate the stability of MINI-EX toward input data sets with low number of cells and its robustness toward missing data, and show that it infers state-of-the-art networks with a better performance compared with other related single-cell network tools. MINI-EX successfully identifies key regulators controlling root development in Arabidopsis and rice, leaf development in Arabidopsis, and ear development in maize, enhancing our understanding of cell-type-specific regulation and unraveling the roles of different regulators controlling the development of specific cell types in plants.

Overexpression of PtoMYB115 improves lignocellulose recalcitrance to enhance biomass digestibility and bioethanol yield by specifically regulating lignin biosynthesis in transgenic poplar

Background

Woody plants provide the most abundant biomass resource that is convertible for biofuels. Since lignin is a crucial recalcitrant factor against lignocellulose hydrolysis, genetic engineering of lignin biosynthesis is considered as a promising solution. Many MYB transcription factors have been identified to involve in the regulation of cell wall formation or phenylpropanoid pathway. In a previous study, we identified that PtoMYB115 contributes to the regulation of proanthocyanidin pathway, however, little is known about its role in lignocellulose biosynthesis and biomass saccharification in poplar.

Results

Here, we detected the changes of cell wall features and examined biomass enzymatic saccharification for bioethanol production under various chemical pretreatments in PtoMYB115 transgenic plants. We reported that PtoMYB115 might specifically regulate lignin biosynthesis to affect xylem development. Overexpression of PtoMYB115 altered lignin biosynthetic gene expression, resulting in reduced lignin deposition, raised S/G and beta-O-4 linkage, resulting in a significant reduction in cellulase adsorption with lignin and an increment in cellulose accessibility. These alterations consequently improved lignocellulose recalcitrance for significantly enhanced biomass saccharification and bioethanol yield in the PtoMYB115-OE transgenic lines. In contrast, the knockout of PtoMYB115 by CRISPR/Cas9 showed reduced woody utilization under various chemical pretreatments.

Conclusions

This study shows that PtoMYB115 plays an important role in specifically regulating lignin biosynthesis and improving lignocellulose features. The enhanced biomass saccharification and bioethanol yield in the PtoMYB115-OE lines suggests that PtoMYB115 is a candidate gene for genetic modification to facilitate the utilization of biomass.

Identification of quantitative trait loci for related traits of stalk lodging resistance using genome-wide association studies in maize (Zea mays L.)

BACKGROUND

Stalk lodging is one of the main factors affecting maize (Zea mays L.) yield and limiting mechanized harvesting. Developing maize varieties with high stalk lodging resistance requires exploring the genetic basis of lodging resistance-associated agronomic traits. Stalk strength is an important indicator to evaluate maize lodging and can be evaluated by measuring stalk rind penetrometer resistance (RPR) and stalk buckling strength (SBS). Along with morphological traits of the stalk for the third internodes length (TIL), fourth internode length (FIL), third internode diameter (TID), and the fourth internode diameter (FID) traits are associated with stalk lodging resistance.

RESULTS

In this study, a natural population containing 248 diverse maize inbred lines genotyped with 83,057 single nucleotide polymorphism (SNP) markers was used for genome-wide association study (GWAS) for six stalk lodging resistance-related traits. The heritability of all traits ranged from 0.59 to 0.72 in the association mapping panel. A total of 85 significant SNPs were identified for the association mapping panel using best linear unbiased prediction (BLUP) values of all traits. Additionally, five candidate genes were associated with stalk strength traits, which were either directly or indirectly associated with cell wall components.

CONCLUSIONS

These findings contribute to our understanding of the genetic basis of maize stalk lodging and provide valuable theoretical guidance for lodging resistance in maize breeding in the future.

Biosynthetic regulatory network of flavonoid metabolites in stems and leaves of Salvia miltiorrhiza

Flavonoid secondary metabolites can treat and prevent many diseases, but systematic studies on regulation of the biosynthesis of such metabolites in aboveground parts of Salvia miltiorrhiza are lacking. In this study, metabonomic and transcriptomic analyses of different S. miltiorrhiza phenotypes were conducted to explore pathways of synthesis, catalysis, accumulation, and transport of the main flavonoid secondary metabolites regulating pigment accumulation. Tissue localization and quantitative analysis of flavonoid secondary metabolites were conducted by laser scanning confocal microscopy (LSCM). A total 3090 differentially expressed genes were obtained from 114,431 full-length unigenes in purple and green phenotypes, and 108 functional genes were involved in flavonoid biosynthesis. Five key phenylpropane structural genes (PAL, 4CL, ANS, 3AT, HCT) were highly differentially expressed, and four transcription factor genes (MYB, WRKY, bHLH, bZiP) were identified. In addition, six GST genes, nine ABC transporters, 22 MATE genes, and three SNARE genes were detected with key roles in flavonoid transport. According to LSCM, flavonoids were mainly distributed in epidermis, cortex, and collenchyma. Thus, comprehensive and systematic analyses were used to determine biosynthesis, accumulation, and transport of flavonoids in stems and leaves of different S. miltiorrhiza phenotypes. The findings will provide a reference for flavonoid production and cultivar selection.

Deciphering the roles of tobacco MYB transcription factors in environmental stress tolerance

The MYB members play important roles in development, metabolism, and stress tolerance in plants. In the current study, a total of 246 tobacco R2R3-MYB transcription factors were identified and systemically analyzed from the latest genome annotation. The newly identified tobacco members were divided into 33 subgroups together with the Arabidopsis members. Furthermore, 44 NtMYB gene pairs were identified to arise from duplication events, which might lead to the expansion of tobacco MYB genes. The expression patterns were revealed by transcriptomic analysis. Notably, the results from phylogenetic analysis, synthetic analysis, and expression analysis were integrated to predict the potential functions of these members. Particularly, NtMYB102 was found to act as the homolog of AtMYB70 and significantly induced by drought and salt treatments. The further assays revealed that NtMYB102 had transcriptional activities, and the overexpression of the encoding gene enhanced the drought and salt stress tolerance in transgenic tobacco. The results of this study may be relevant for future functional analyses of the MYB genes in tobacco.

Transcriptional Reactivation of Lignin Biosynthesis for the Heterologous Production of Etoposide Aglycone in Nicotiana benthamiana

Nicotiana benthamiana is a valuable plant chassis for heterologous production of medicinal plant natural products. This host is well suited for the processing of organelle-localized plant enzymes, and the conservation of the primary metabolism across the plant kingdom often provides required plant-specific precursor molecules that feed a given pathway. Despite this commonality in metabolism, limited precursor supply and/or competing host pathways can interfere with yields of heterologous products. Here, we use transient transcriptional reprogramming of endogenous N. benthamiana metabolism to drastically improve flux through the etoposide pathway derived from the medicinal plant Podophyllum spp. Specifically, coexpression of a single lignin-associated transcription factor, MYB85, with pathway genes results in unprecedented levels of heterologous product accumulation in N. benthamiana leaves: 1 mg/g dry weight (DW) of the etoposide aglycone, 35 mg/g DW (-)-deoxypodophyllotoxin, and 3.5 mg/g DW (-)-epipodophyllotoxin─up to two orders of magnitude above previously reported biosynthetic yields for the etoposide aglycone and eight times higher than what is observed for (-)-deoxypodophyllotoxin in the native medicinal plant. Unexpectedly, transient activation of lignin metabolism by transcription factor overexpression also reduces the production of undesired side products that likely result from competing N. benthamiana metabolism. Our work demonstrates that synthetic activation of lignin biosynthesis in leaf tissue is an effective strategy for optimizing the production of medicinal compounds derived from phenylpropanoid precursors in the plant chassis N. benthamiana. Furthermore, our results highlight the engineering value of MYB85, an early switch in lignin biosynthesis, for on-demand modulation of monolignol flux and support the role of MYB46 as a master regulator of lignin polymer deposition.

Transcriptional and metabolic changes associated with internode development and reduced cinnamyl alcohol dehydrogenase activity in sorghum

The molecular mechanisms associated with secondary cell wall (SCW) deposition in sorghum remain largely uncharacterized. Here, we employed untargeted metabolomics and large-scale transcriptomics to correlate changes in SCW deposition with variation in global gene expression profiles and metabolite abundance along an elongating internode of sorghum, with a major focus on lignin and phenolic metabolism. To gain deeper insight into the metabolic and transcriptional changes associated with pathway perturbations, a bmr6 mutant [with reduced cinnamyl alcohol dehydrogenase (CAD) activity] was analyzed. In the wild type, internode development was accompanied by an increase in the content of oligolignols, p-hydroxybenzaldehyde, hydroxycinnamate esters, and flavonoid glucosides, including tricin derivatives. We further identified modules of genes whose expression pattern correlated with SCW deposition and the accumulation of these target metabolites. Reduced CAD activity resulted in the accumulation of hexosylated forms of hydroxycinnamates (and their derivatives), hydroxycinnamaldehydes, and benzenoids. The expression of genes belonging to one specific module in our co-expression analysis correlated with the differential accumulation of these compounds and contributed to explaining this metabolic phenotype. Metabolomics and transcriptomics data further suggested that CAD perturbation activates distinct detoxification routes in sorghum internodes. Our systems biology approach provides a landscape of the metabolic and transcriptional changes associated with internode development and with reduced CAD activity in sorghum.

Comparative analysis of R2R3-MYB transcription factors in the flower of Iris laevigata identifies a novel gene regulating tobacco cold tolerance

Breeding for flower cold resistance is a priority for flower breeding research in northern China. The identification of cold resistance genes will not only provide genetic resources for cold resistance breeding, but also form a basis for the study of plant cold resistance mechanisms. Based on the flower transcriptome of Iris laevigata, 20 R2R3-MYBs were identified and comprehensive analysis, including conservative domain, phylogenetic analyses and functional distribution, were performed for R2R3-MYBs. Expression patterns of the abiotic stress genes under cold stress were detected, the upregulated gene was genetically transformed into tobacco, and the related physiological indicators of the transgenic tobacco were measured. A novel cold resistance gene, IlMYB306, was obtained. qRT-PCR indicated that IlMYB306 was dramatically induced by cold stress and was significantly upregulated in roots. The free proline content, MDA, SOD and POD activity of the transgenic tobacco improved after cold stress, and the chlorophyll content decreased slowly. In addition, overexpression of IlMYB306 improved cold resistance of the seeds. SEM results showed leaves of transgenic tobacco had obvious folds, more grooves and bulges on the lower leaf surface. Overall, we report a novel cold resistance R2R3-MYB gene, IlMYB306, in the flower of I. laevigata, which could improve tobacco cold stress tolerance by thickening the waxy layer, increasing antioxidant activity and the content of proline.

The Regulation of Xylem Development by Transcription Factors and Their Upstream MicroRNAs

Xylem, as a unique organizational structure of vascular plants, bears water transport and supports functions necessary for plant survival. Notably, secondary xylem in the stem (i.e., wood) also has important economic and ecological value. In view of this, the regulation of xylem development has been widely concerned. In recent years, studies on model plants Arabidopsis and poplar have shown that transcription factors play important regulatory roles in various processes of xylem development, including the directional differentiation of procambium and cambium into xylem, xylem arrangement patterns, secondary cell wall formation and programmed cell death. This review focuses on the regulatory roles of widely and thoroughly studied HD-ZIP, MYB and NAC transcription factor gene families in xylem development, and it also pays attention to the regulation of their upstream microRNAs. In addition, the existing questions in the research and future research directions are prospected.

Genome-Wide Identification and Analysis of the R2R3-MYB Gene Family in Theobroma cacao

The MYB gene family is involved in the regulation of plant growth, development and stress responses. In this paper, to identify Theobroma cacao R2R3-MYB (TcMYB) genes involved in environmental stress and phytohormones, we conducted a genome-wide analysis of the R2R3-MYB gene family in Theobroma cacao (cacao). A total of 116 TcMYB genes were identified, and they were divided into 23 subgroups according to the phylogenetic analysis. Meanwhile, the conserved motifs, gene structures and cis-acting elements of promoters were analyzed. Moreover, these TcMYB genes were distributed on 10 chromosomes. We conducted a synteny analysis to understand the evolution of the cacao R2R3-MYB gene family. A total of 37 gene pairs of TcMYB genes were identified through tandem or segmental duplication events. Additionally, we also predicted the subcellular localization and physicochemical properties. All the studies showed that TcMYB genes have multiple functions, including responding to environmental stresses. The results provide an understanding of R2R3-MYB in Theobroma cacao and lay the foundation for a further functional analysis of TcMYB genes in the growth of cacao.

Analyzing lignin biosynthesis pathways in rattan using improved co-expression networks of NACs and MYBs

BACKGROUND

The rattan is a valuable plant resource with multiple applications in tropical forests. Calamus simplicifolius and Daemonorops jenkinsiana are the two most representative rattan species, supplying over 95% of the raw materials for the rattan industry. Hence, the wood properties of both rattans have always attracted researchers' attention.

RESULTS

We re-annotated the genomes, obtained 81 RNA-Seq datasets, and developed an improved pipeline to increase the reliability of co-expression networks of both rattans. Based on the data and pipeline, co-expression relationships were detected in 11 NACs, 49 MYBs, and 86 lignin biosynthesis genes in C. simplicifolius and four NACs, 59 MYBs, and 76 lignin biosynthesis genes in D. jenkinsiana, respectively. Among these co-expression pairs, several genes had a close relationship to the development of wood properties. Additionally, we detected the enzyme gene on the lignin biosynthesis pathway was regulated by either NAC or MYB, while LACCASES was regulated by both NAC and MYB. For D. jenkinsiana, the lignin biosynthesis regulatory network was characterized by positive regulation, and MYB possible negatively regulate non-expressed lignin biosynthesis genes in stem tissues. For C. simplicifolius, NAC may positively regulate highly expressed genes and negatively regulate non-expressed lignin biosynthesis genes in stem tissues. Furthermore, we established core regulatory networks of NAC and MYB for both rattans.

CONCLUSIONS

This work improved the accuracy of rattan gene annotation by integrating an efficient co-expression network analysis pipeline, enhancing gene coverage and accuracy of the constructed network, and facilitating an understanding of co-expression relationships among NAC, MYB, and lignin biosynthesis genes in rattan and other plants.

Isolation and Characterization of an LBD Transcription Factor CsLBD39 from Tea Plant (Camellia sinensis) and Its Roles in Modulating Nitrate Content by Regulating Nitrate-Metabolism-Related Genes

Nitrate nitrogen is an important nitrogen source for tea plants’ growth and development. LBD transcription factors play important roles in response to the presence of nitrate in plants. The functional study of LBD transcription factors in tea plants remains limited. In this study, the LBD family gene CsLBD39 was isolated and characterized from tea plants. Sequence analysis indicated that CsLBD39 contained a highly conserved CX₂CX₆CX₃CX domain. The phylogenetic tree assay showed that CsLBD39 belonged to class II subfamily of the LBD family. CsLBD39 was highly expressed in flowers and root; we determined that its expression could be induced by nitrate treatment. The CsLBD39 protein was located in the nucleus and has transcriptional activation activity in yeast. Compared with the wild type, overexpression of CsLBD39 gene in Arabidopsis resulted in smaller rosettes, shorter main roots, reduced lateral roots and lower plant weights. The nitrate content and the expression levels of genes related to nitrate transport and regulation were decreased in transgenic Arabidopsis hosting CsLBD39 gene. Compared with the wild type, CsLBD39 overexpression in transgenic Arabidopsis had smaller cell structure of leaves, shorter diameter of stem cross section, and slender and compact cell of stem longitudinal section. Under KNO₃ treatment, the contents of nitrate, anthocyanins, and chlorophyll in leaves, and the content of nitrate in roots of Arabidopsis overexpressing CsLBD39 were reduced, the expression levels of nitrate transport and regulation related genes were decreased. The results revealed that CsLBD39 may be involved in nitrate signal transduction in tea plants as a negative regulator and laid the groundwork for future studies into the mechanism of nitrate response.

The MicroRNA397a-LACCASE17 module regulates lignin biosynthesis in Medicago ruthenica (L.)

Mechanical strength is essential for the upright growth habit, which is one of the most important characteristics of terrestrial plants. Lignin, a phenylpropanoid-derived polymer mainly present in secondary cell walls plays critical role in providing mechanical support. Here, we report that the prostrate-stem cultivar of the legume forage Medicago ruthenica cultivar 'Mengnong No. 1' shows compromised mechanical strength compared with the erect-stem cultivar 'Zhilixing'. The erect-stem cultivar, 'Zhilixing' has significantly higher lignin content, leading to higher mechanical strength than the prostrate-stem cultivar. The low abundance of miRNA397a in the Zhiixing cultivar causes reduced cleavage of MrLAC17 transcript, which results in enhanced expression level of MrLAC17 compared to that in the prostrate-stem cultivar Mengnong No. 1. Complementation of the Arabidopsis lac4 lac17 double mutants with MrLAC17 restored the lignin content to wild-type levels, confirming that MrLAC17 perform an exchangeable role with Arabidopsis laccases. LAC17-mediated lignin polymerization is therefore increased in the 'Zhilixing', causing the erect stem phenotype. Our data reveal the importance of the miR397a in the lignin biosynthesis and suggest a strategy for molecular breeding targeting plant architecture in legume forage.

MicroRNA828 negatively regulates lignin biosynthesis in stem of Populus tomentosa through MYB targets

Lignin biosynthesis in the sclerenchyma cells is strictly controlled by a complex network of genetic and environmental signals. In the last decades, the transcriptional regulation of lignin synthesis in woody species has been established. However, the role of microRNA-mediated post-transcriptional modulation in secondary cell wall biosynthesis remains poorly understood. Here, we identified a microRNA, miR828, involved in the regulation specific to lignin biosynthesis during stem development in Populus tomentosa Carr. miR828 is preferentially expressed in the secondary vascular tissues during stem development. Two MYB genes (MYB171 and MYB011) were validated as direct targets of miR828 by degradome analysis and green fluorescent protein signal detection. Overexpression of miR828 in poplar downregulated genes for lignin biosynthesis, resulting in reduced lignin content in cell walls. Conversely, suppression of miR828 in plants by the short tandem target mimics elevated the expression of lignin biosynthetic genes and increased lignin deposition. We further revealed that poplar MYB171, as the most abundant miR828 target in the stem, is a positive regulator for lignin biosynthesis. Transient expression assays showed that both MYB171 and MYB011 activated PAL1 and CCR2 transcription, whereas the introduction of miR828 significantly suppressed their expression that was induced by MYB171 or MYB011. Collectively, our results demonstrate that the miR828-MYBs module precisely regulates lignin biosynthesis during the stem development in P. tomentosa through transcriptional and post-transcriptional manners.

Tissue-specific transcriptomics reveals a central role of CcNST1 in regulating the fruit lignification pattern in Camellia chekiangoleosa, a woody oil-crop

Fruit lignification is of significant economic importance because it affects the quality of fruit and the production of seed oil. The specified lignification pattern in Camellia chekiangoleosa fruits plays critical roles in its seed oil yield, but little is known about how this lignification process is regulated. Here, we report on a comprehensive tissue-specific transcriptomics analysis conducted for C. chekiangoleosa fruit. By mining the differentially expressed genes, we found that lignin biosynthesis and transcriptional regulation pathways were significantly enriched in the lignified tissues. The homolog of NST-like transcription factor, CcNST1, was highly expressed in lignified seed coat and endocarp tissues; transgenic analyses of CcNST1 in Arabidopsis and hybrid poplar revealed the enhanced lignification levels of various tissues. Gene expression analysis of the transgenic lines uncovered potential downstream genes involved in the regulation of lignin biosynthesis. This work provides a valuable gene expression resource and identified the pivotal role of CcNST1 in regulating the lignin biosynthesis underlying fruit lignification.

High temperature induces male sterility via MYB66-MYB4-Casein kinase I signaling in cotton

High temperature (HT) causes male sterility and decreases crop yields. Our previous works have demonstrated that sugar and auxin signaling pathways, Gossypium hirsutum Casein kinase I (GhCKI), and DNA methylation are all involved in HT-induced male sterility in cotton. However, the signaling mechanisms leading to distinct GhCKI expression patterns induced by HT between HT-tolerant and HT-sensitive cotton anthers remain largely unknown. Here, we identified a GhCKI promoter (ProGhCKI) region that functions in response to HT in anthers and found the transcription factor GhMYB4 binds to this region to act as an upstream positive regulator of GhCKI. In the tapetum of early-stage cotton anthers, upregulated expression of GhMYB4 under HT and overexpressed GhMYB4 under normal temperature both led to severe male sterility phenotypes, coupled with enhanced expression of GhCKI. We also found that GhMYB4 interacts with GhMYB66 to form a heterodimer to enhance its binding to ProGhCKI. However, GhMYB66 showed an expression pattern similar to GhMYB4 under HT but did not directly bind to ProGhCKI. Furthermore, HT reduced siRNA-mediated CHH DNA methylations in the GhMYB4 promoter, which enhanced the expression of GhMYB4 in tetrad stage anthers and promoted the formation of the GhMYB4/GhMYB66 heterodimer, which in turn elevated the transcription of GhCKI in the tapetum, leading to male sterility. Overall, we shed light on the GhMYB66-GhMYB4-GhCKI regulatory pathway in response to HT in cotton anthers.

Structural and Functional Analyses of Hub MicroRNAs in An Integrated Gene Regulatory Network of Arabidopsis

MicroRNAs (miRNAs) are trans-acting small regulatory RNAs that work coordinately with transcription factors (TFs) to shape the repertoire of cellular mRNAs available for translation. Despite our growing knowledge of individual plant miRNAs, their global roles in gene regulatory networks remain mostly unassessed. Based on interactions obtained from public databases and curated from the literature, we reconstructed an integrated miRNA network in Arabidopsis that includes 66 core TFs, 318 miRNAs, and 1712 downstream genes. We found that miRNAs occupy distinct niches and enrich miRNA-containing feed-forward loops (FFLs), particularly those with miRNAs as intermediate nodes. Further analyses revealed that miRNA-containing FFLs coordinate TFs located in different hierarchical layers and that intertwined miRNA-containing FFLs are associated with party and date miRNA hubs. Using the date hub MIR858A as an example, we performed detailed molecular and genetic analyses of three interconnected miRNA-containing FFLs. These analyses revealed individual functions of the selected miRNA-containing FFLs and elucidated how the date hub miRNA fulfills multiple regulatory roles. Collectively, our findings highlight the prevalence and importance of miRNA-containing FFLs, and provide new insights into the design principles and control logics of miRNA regulatory networks governing gene expression programs in plants.

Biomass formation and sugar release efficiency of Populus modified by altered expression of a NAC transcription factor

Woody biomass is an important feedstock for biofuel production. Manipulation of wood properties that enable efficient conversion of biomass to biofuel reduces cost of biofuel production. Wood cell wall composition is regulated at several levels that involve expression of transcription factors such as wood-/secondary cell wall-associated NAC domains (WND or SND). In Arabidopsis thaliana, SND1 regulates cell wall composition through activation of its down-stream targets such as MYBs. The functional aspects of SND1 homologs in the woody Populus have been studied through transgenic manipulation. In this study, we investigated the role of PdWND1B, Populus SND1 sequence ortholog, in wood formation using transgenic manipulation through over-expression or silencing under the control of a vascular-specific 4-coumarate-CoA ligase (4CL) promoter. As compared with control plants, PdWND1B-RNAi plants were shorter in height, with significantly reduced stem diameter and dry biomass, whereas there were no significant differences in growth and productivity of PdWND1B over-expression plants. Conversely, PdWND1B over-expression lines showed a significant reduction in cellulose and increase in lignin content, whereas there was no significant impact on lignin content of downregulated lines. Stem carbohydrate composition analysis revealed a decrease in glucose, mannose, arabinose, and galactose, but an increase in xylose in the over-expression lines. Transcriptome analysis revealed upregulation of several downstream transcription factors and secondary cell wall related structural genes in the PdWND1B over-expression lines, partly explaining the observed phenotypic changes in cell wall chemistry. Relative to the control, glucose release efficiency and ethanol production from stem biomass was significantly reduced in over-expression lines. Our results show that PdWND1B is an important factor determining biomass productivity, cell wall chemistry and its conversion to biofuels in Populus.

New Insights into Bacillus-Primed Plant Responses to a Necrotrophic Pathogen Derived from the Tomato-Botrytis Pathosystem

Induced systemic resistance (ISR) is one of the most studied mechanisms of plant−microbe interaction and is considered a very promising alternative for integrated pest management programs. In our study, we explored the plant defense response induced by Bacillus velezensis BBC047 in relation to its application before or after Botrytis cinerea infection of tomato plants. The inoculation of BBC047 did not considerably alter the gene expression of the tomato tissues, whereas infection with B. cinerea in BBC047-primed plants induced expression of LRR and NBS-LRR receptors, which are highly related to the ISR response. As expected, B. cinerea infection generated molecular patterns typical of a defense response to pathogen infection as the overexpression of pathogenesis-related proteins (PRs) in leaflets distant to the point of infection. The curative treatment (P + F + B) allowed us to gain insights into plant response to an inverted priming. In this treatment, B. cinerea caused the m tissue damage, extending nearly entirely across the entire infected leaves. Additionally, genes generally associated with early SAR response (<16 h) were overexpressed, and apparently, the beneficial strain was not perceived as such. Therefore, we infer that the plant defense to the curative treatment represents a higher degree of biological stress triggered by the incorporation of strain BBC047 as second arriving microorganism. We highlight the importance the phytosanitary status of plants prior to inoculation of beneficial microorganism for the biocontrol of pathogens.

Genome-wide analysis of the CAD gene family reveals two bona fide CAD genes in oil palm

Cinnamyl alcohol dehydrogenase (CAD) is the key enzyme for lignin biosynthesis in plants. In this study, genome-wide analysis was performed to identify CAD genes in oil palm (Elaeis guineensis). Phylogenetic analysis was then conducted to select the bona fide EgCADs. The bona fide EgCAD genes and their respective 5' flanking regions were cloned and analysed. Their expression profiles were evaluated in various organs using RT-PCR. Seven EgCAD genes (EgCAD1-7) were identified and divided into four phylogenetic groups. EgCAD1 and EgCAD2 display high sequence similarities with other bona fide CADs and possess all the signature motifs of the bona fide CAD. They also display similar 3D protein structures. Gene expression analysis showed that EgCAD1 was expressed most abundantly in the root tissues, while EgCAD2 was expressed constitutively in all the tissues studied. EgCAD1 possesses only one transcription start site, while EgCAD2 has five. Interestingly, a TC microsatellite was found in the 5' flanking region of EgCAD2. The 5' flanking regions of EgCAD1 and EgCAD2 contain lignin-associated regulatory elements i.e. AC-elements, and other defence-related motifs, including W-box, GT-1 motif and CGTCA-motif. Altogether, these results imply that EgCAD1 and EgCAD2 are bona fide CAD involved in lignin biosynthesis during the normal development of oil palm and in response to stresses. Our findings shed some light on the roles of the bona fide CAD genes in oil palm and pave the way for manipulating lignin content in oil palm through a genetic approach.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03208-0.

PtrWOX13A Promotes Wood Formation and Bioactive Gibberellins Biosynthesis in Populus trichocarpa

WUSCHEL-related homeobox (WOX) genes are plant-specific transcription factors (TFs) involved in multiple processes of plant development. However, there have hitherto no studies on the WOX TFs involved in secondary cell wall (SCW) formation been reported. In this study, we identified a Populus trichocarpa WOX gene, PtrWOX13A, which was predominantly expressed in SCW, and then characterized its functions through generating PtrWOX13A overexpression poplar transgenic lines; these lines exhibited not only significantly enhanced growth potential, but also remarkably increased SCW thicknesses, fiber lengths, and lignin and hemicellulose contents. However, no obvious change in cellulose content was observed. We revealed that PtrWOX13A directly activated its target genes through binding to two cis-elements, ATTGATTG and TTAATSS, in their promoter regions. The fact that PtrWOX13A responded to the exogenous GAs implies that it is responsive to GA homeostasis caused by GA inactivation and activation genes (e.g., PtrGA20ox4, PtrGA2ox1, and PtrGA3ox1), which were regulated by PtrWOX13A directly or indirectly. Since the master switch gene of SCW formation, PtrWND6A, and lignin biosynthesis regulator, MYB28, significantly increased in PtrWOX13A transgenic lines, we proposed that PtrWOX13A, as a higher hierarchy TF, participated in SCW formation through controlling the genes that are components of the known hierarchical transcription regulation network of poplar SCW formation, and simultaneously triggering a gibberellin-mediated signaling cascade. The discovery of PtrWOX13A predominantly expressed in SCW and its regulatory functions in the poplar wood formation has important implications for improving the wood quality of trees via genetic engineering.

Time-course transcriptome analysis identifies rewiring patterns of transcriptional regulatory networks in rice under Rhizoctonia solani infection

Sheath Blight (SB) disease in rice is caused by the infection from the fungal pathogen Rhizoctonia solani (R. solani). SB is one of the most severe rice diseases that can cause up to 50% yield losses in rice. Naturally occurring rice varieties resistant to SB have not been reported yet. We have performed a Time-Series RNA-Seq analysis on a widely cultivated rice variety BPT-5204 for identifying transcriptome level response signatures during R. solani infection at 1st, 2nd and 5th day post infection (dpi). In total, 428, 3225 and 1225 genes were differentially expressed in the treated rice plants on 1, 2 and 5 dpi, respectively. GO and KEGG enrichment analysis identified significant processes and pathways differentially altered in the rice plants during the fungal infection. Machine learning and network based integrative approach was used to construct rice Transcriptional Regulatory Networks (TRNs) for the three time points. TRN analysis identified SUB1B, MYB30 and CCA1 as important regulatory hub transcription factors in rice during R. solani infection. Jasmonic acid, salicylic acid, ethylene biogenesis and signaling were induced on infection. SAR was up regulated, while photosynthesis and carbon fixation processes were significantly down regulated. Involvement of MAPK, CYPs, peroxidase, PAL, chitinase genes were also observed in response to the fungal infection. The integrative analysis identified seven putative SB resistance genes differentially regulated in rice during R. solani infection.

A natural mutation of the NST1 gene arrests secondary cell wall biosynthesis in the seed coat of a hull-less pumpkin accession

Hull-less pumpkins (Cucurbita pepo L.) are naturally occurring novel variants known as oilseed or naked-seeded pumpkins, and are characterized by the absence of a normal lignified seed coat. Due to a specialized seed coat structure, these variants serve as a good model for studying seed coat formation and simplify the processing of pumpkin seeds. However, causal genes for this hull-less trait still remain unknown. Here, by bulked segregant analysis and fine mapping, we found that mutation of a single gene, NAC SECONDARY WALL THICKENING PROMOTING FACTOR 1 (NST1), accounts for the hull-less trait. A 14-bp sequence insertion in the CpNST1 gene causes premature termination of CpNST1 translation, leading to lack of secondary cell wall (SCW) biosynthesis in hull-less seed coats. In situ hybridization analysis provided further evidence for the role of CpNST1 in pumpkin seed coat SCW biosynthesis. Interestingly, through secondary cell wall compositional analysis, we found that the main SCW components differed among cell layers in the seed coat. RNA-seq analysis indicated an upstream role of CpNST1 in the SCW biosynthesis network. Collectively, our findings provide mechanistic insight into seed coat SCW biosynthesis, and a target gene for breeders to introduce this hull-less trait for commercial exploitation.

Comprehensive analysis of the MYB transcription factor gene family in Morus alba

BACKGROUND

The V-myb myeloblastosis viral oncogene homolog (MYB) family of proteins is large, containing functionally diverse transcription factors. However, MYBs in Morus are still poorly annotated and a comprehensive functional analysis of these transcription factors is lacking.

RESULTS

In the present study, a genome-wide identification of MYBs in Morus alba was performed. In total 166 MaMYBs were identified, including 103 R2R3-MYBs and four 3R-MaMYBs. Comprehensive analyses, including the phylogenetic analysis with putative functional annotation, motif and structure analysis, gene structure organization, promoter analysis, chromosomal localization, and syntenic relationships of R2R3-MaMYBs and 3R-MaMYBs, provided primary characterization for these MaMYBs. R2R3-MaMYBs covered the subgroups reported for R2R3-MYBs in Arabidopsis and Populus, and had two Morus-specific subgroups, indicating the high retention of MYBs in Morus. Motif analysis revealed high conservative residues at the start and end of each helix and residues consisting of the third helix in R2 and R3 repeats. Thirteen intron/exon patterns (a-m) were summarized, and the intron/exon pattern of two introns with phase numbers of 0 and 2 was the prevalent pattern for R2R3-MaMYBs. Various cis-elements in promoter regions were identified, and were mainly related to light response, development, phytohormone response, and abiotic and biotic stress response and secondary metabolite production. Expression patterns of R2R3-MaMYBs in different organs showed that MaMYBs involved in secondary cell wall components and stress responsiveness were preferentially expressed in roots or stems. R2R3-MaMYBs involved in flavonoid biosynthesis and anthocyanin accumulation were identified and characterized based on functional annotation and correlation of their expression levels with anthocyanin contents.

CONCLUSION

Based on a comprehensive analysis, this work provided functional annotation for R2R3-MYBs and an informative reference for further functional dissection of MYBs in Morus.

Dissection of the genetic architecture of peduncle vascular bundle-related traits in maize by a genome-wide association study

The peduncle vascular system of maize is critical for the transport of photosynthetic products, nutrients, and water from the roots and leaves to the ear. Accordingly, it positively affects the grain yield. However, the genetic basis of peduncle vascular bundle (PVB)-related traits in maize remains unknown. Thus, 15 PVB-related traits of 386 maize inbred lines were investigated at three locations (Yongcheng, 17YC; Kaifeng, 20KF; and Yuanyang, 20YY). The repeatability for the 15 traits ranged from 35.53% to 92.13%. A genome-wide association study was performed and 69 non-redundant quantitative trait loci (QTL) were detected, including 9, 41, and 27 QTL identified at 17YC, 20KF, and 20YY, respectively. These QTL jointly explained 4.72% (SLL) to 37.30% (NSVB) of the phenotypic variation. Eight QTL were associated with the same trait at two locations. Furthermore, four pleiotropic QTL were identified. Moreover, one QTL (qPVB44), associated with NSVB_20KF, was co-localized with a previously reported locus related to kernel width, implying qPVB44 may affect the kernel width by modulating the number of small vascular bundles. Examinations of the 69 QTL identified 348 candidate genes that were classified in five groups. Additionally, 26 known VB-related homologous genes (e.g. VLN2, KNOX1, and UGT72B3) were detected in 20 of the 69 QTL. A comparison of the NSVB between a Zmvln2 EMS mutant and its wild type elucidated the function of the candidate gene ZmVLN2. These results are important for clarifying the genetic basis of PVB-related traits and may be useful for breeding new high-yielding maize cultivars.

Comprehensive analysis of 1R- and 2R-MYBs reveals novel genic and protein features, complex organisation, selective expansion and insights into evolutionary tendencies

Myeloblastosis (MYB) family, the largest plant transcription factor family, has been subcategorised based on the number and type of repeats in the MYB domain. In spite of several reports, evolution of MYB genes and repeats remains enigmatic. Brassicaceae members are endowed with complex genomes, including dysploidy because of its unique history with multiple rounds of polyploidisation, genomic fractionations and rearrangements. The present study is an attempt to gain insights into the complexities of MYB family diversity, understand impacts of genome evolution on gene families and develop an evolutionary framework to understand the origin of various subcategories of MYB gene family. We identified and analysed 1129 MYBs that included 1R-, 2R-, 3R- and atypical-MYBs across sixteen species representing protists, fungi, animals and plants and exclude MYB identified from Brassicaceae except Arabidopsis thaliana; in addition, a total of 1137 2R-MYB genes from six Brassicaceae species were also analysed. Comparative analysis revealed predominance of 1R-MYBs in protists, fungi, animals and lower plants. Phylogenetic reconstruction and analysis of selection pressure suggested ancestral nature of R1-type repeat containing 1R-MYBs that might have undergone intragenic duplication to form multi-repeat MYBs. Distinct differences in gene structure between 1R-MYB and 2R-MYBs were observed regarding intron number, the ratio of gene length to coding DNA sequence (CDS) length and the length of exons encoding the MYB domain. Conserved as well as novel and lineage-specific intron phases were identified. Analyses of physicochemical properties revealed drastic differences indicating functional diversification in MYBs. Phylogenetic reconstruction of 1R- and 2R-MYB genes revealed a shared structure-function relationship in clades which was supported when transcriptome data was analysed in silico. Comparative genomics to study distribution pattern and mapping of 2R-MYBs revealed congruency and greater degree of synteny and collinearity among closely related species. Micro-synteny analysis of genomic segments revealed high conservation of genes that are immediately flanking the surrounding tandemly organised 2R-MYBs along with instances of local duplication, reorganisations and genome fractionation. In summary, polyploidy, dysploidy, reshuffling and genome fractionation were found to cause loss or gain of 2R-MYB genes. The findings need to be supported with functional validation to understand gene structure-function relationship along the evolutionary lineage and adaptive strategies based on comparative functional genomics in plants.

Full-Length Transcriptome Sequencing and Comparative Transcriptomic Analyses Provide Comprehensive Insight Into Molecular Mechanisms of Cellulose and Lignin Biosynthesis in Cunninghamia lanceolata

Cunninghamia lanceolata is an essential timber species that provide 20%-30% raw materials for China's timber industry. Although a few transcriptomes have been published in C. lanceolata, full-length mRNA transcripts and regulatory mechanisms behind the cellulose and lignin biosynthesis have not been thoroughly investigated. Here, PacBio Iso-seq and RNA-seq analyses were adapted to identify the full-length and differentially expressed transcripts along a developmental gradient from apex to base of C. lanceolata shoots. A total of 48,846 high-quality full-length transcripts were obtained, of which 88.0% are completed transcriptome based on benchmarking universal single-copy orthologs (BUSCO) assessment. Along stem developmental gradient, 18,714 differentially expressed genes (DEGs) were detected. Further, 28 and 125 DEGs were identified as enzyme-coding genes of cellulose and lignin biosynthesis, respectively. Moreover, 57 transcription factors (TFs), including MYB and NAC, were identified to be involved in the regulatory network of cellulose and lignin biosynthesis through weighted gene co-expression network analysis (WGCNA). These TFs are composed of a comparable regulatory network of secondary cell wall formation in angiosperms, revealing a similar mechanism may exist in gymnosperms. Further, through qRT-PCR, we also investigated eight specific TFs involved in compression wood formation. Our findings provide a comprehensive and valuable source for molecular genetics breeding of C. lanceolata and will be beneficial for molecular-assisted selection.

Transcriptome Analysis Reveals Candidate Lignin-Related Genes and Transcription Factors during Fruit Development in Pomelo (Citrus maxima)

Juice sac granulation (a physiological disorder) leads to large postharvest losses of pomelo (Citrus maxima). Previous studies have shown that juice sac granulation is closely related to lignin accumulation, while the molecular mechanisms underlying this disorder remain elusive in pomelo. Our results showed that the lignin content in NC (near the core) and FC (far away from the core) juice sacs overall increased from 157 DPA (days post anthesis) to 212 DPA and reached a maximum at 212 DPA. Additionally, the lignin content of NC juice sacs was higher than that of FC juice sacs. In this study, we used transcriptome-based weighted gene co-expression network analysis (WGCNA) to address how lignin formation in NC and FC juice sacs is generated during the development of pomelo. After data assembly and bioinformatic analysis, we found a most correlated module (black module) to the lignin content, then we used the 11 DEGs in this module as hub genes for lignin biosynthesis. Among these DEGs, PAL (phenylalanine ammonia lyase), HCT (hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase), 4CL2 (4-coumarate: CoA ligase), C4H (cinnamate 4-hydroxylase), C3'H (p-coumarate 3-hydroxylase), and CCoAOMT1 (caffeoyl CoA 3-Omethyltransferase) were the most distinct DEGs in granulated juice sacs. Co-expression analysis revealed that the expression patterns of several transcription factors such as MYB, NAC, OFP6, and bHLH130 are highly correlated with lignin formation. In addition, the expression patterns of the DEGs related to lignin biosynthesis and transcription factors were validated by qRT-PCR, and the results were highly concordant with the RNA-seq results. These results would be beneficial for further studies on the molecular mechanism of lignin accumulation in pomelo juice sacs and would help with citrus breeding.

Genetic loci associated with freezing tolerance in a European rapeseed (Brassica napus L.) diversity panel identified by genome-wide association mapping

Winter biotypes of rapeseed (Brassica napus L.) require a vernalization treatment to enter the reproductive phase and generally produce greater yields than spring rapeseed. To find genetic loci associated with freezing tolerance in rapeseed, we first performed genotyping-by-sequencing (GBS) on a diversity panel consisting of 222 rapeseed accessions originating primarily from Europe, which identified 69,554 high-quality single-nucleotide polymorphisms (SNPs). Model-based cluster analysis suggested that there were eight subgroups. The diversity panel was then phenotyped for freezing survival (visual damage and Fv/Fo and Fv/Fm) after 2 months of cold acclimation (5°C) and a freezing treatment (-15°C for 4 h). The genotypic and phenotypic data for each accession in the rapeseed diversity panel was then used to conduct a genome-wide association study (GWAS). GWAS results showed that 14 significant markers were mapped to seven chromosomes for the phenotypes scored. Twenty-four candidate genes located within the mapped loci were identified as previously associated with lipid, photosynthesis, flowering, ubiquitination, and cytochrome P450 in rapeseed or other plant species.

Diverse roles of MYB transcription factors in regulating secondary metabolite biosynthesis, shoot development, and stress responses in tea plants (Camellia sinensis)

Tea (Camellia sinensis) is concocted from tea plant shoot tips that produce catechins, caffeine, theanine, and terpenoids, which collectively determine the rich flavors and health benefits of the infusion. However, little is known about the integrated regulation of shoot tip development and characteristic secondary metabolite biosynthesis in tea plants. Here, we demonstrate that MYB transcription factors (TFs) play key and yet diverse roles in regulating leaf and stem development, secondary metabolite biosynthesis, and environmental stress responses in tea plants. By integrating transcriptomic and metabolic profiling data in different tissues at a series of developmental stages or under various stress conditions, alongside biochemical and genetic analyses, we predicted the MYB TFs involved in regulating shoot development (CsMYB2, 98, 107, and 221), epidermal cell initiation (CsMYB184, 41, 139, and 219), stomatal initiation (CsMYB113 and 153), and the biosynthesis of flavonoids (including catechins, anthocyanins, and flavonols; CsMYB8 and 99), caffeine (CsMYB85 and 86), theanine (CsMYB9 and 49), carotenoids (CsMYB110), mono-/sesquiterpenoid volatiles (CsMYB68, 147, 148, and 193), lignin (CsMYB164 and 192), and indolic compounds (CsMYB139, 162, and 198), as well as the MYB TFs that are likely involved in hormone signaling-mediated environmental stress and defense responses. We characterized the functions of some key MYBs in regulating flavonoid and carotenoid biosynthesis for tea quality and flavor. This study provides a cross-family analysis of MYBs in tea alongside new insights into the coordinated regulation of tea plant shoot development and secondary metabolism, paving the way towards understanding of tea quality trait formation and genetic improvement of quality tea plants.

MYB43 as a novel substrate for CRL4PRL1 E3 ligases negatively regulates cadmium tolerance through transcriptional inhibition of HMAs in Arabidopsis

Controlled stability of proteins is a highly efficient mechanism to direct diverse processes in plants. A key regulatory system for protein stability is given by the CULLIN-RING E3 ligases (CRLs). In this work, MYB43 is identified as a novel target of a CUL4-DDB1-PRL1 (PLEIOTROPIC REGULATORY LOCUS 1)-RING E3 ligase (CRL4^PRL1 E3 ligase). Its stability depends on the presence of PRL1, a WD40-containing protein functioning as a substrate receptor of the CRL4 E3 ligases. Genetic studies have indicated that MYB43 is a negative regulator of cadmium (Cd) tolerance in Arabidopsis by transcriptional inhibition of important Cd transporters (HMA2, HMA3 and HMA4), while PRL1 and CUL4 positively regulate Cd tolerance. Expression of CUL4 and PRL1 was enhanced in response to Cd stress, and PRL1 can interact with and target MYB43 for degradation depending on assembly of CRL4^PRL1 E3 ligase, and consequently increase the expression of HMA2, HMA3 and HMA4 through attenuating the transcriptional inhibition. HMA2 and HMA4 are shown to transport cadmium ion (Cd²⁺ ) from the roots of plants to the shoots through the xylem, ultimately increasing the plants' tolerance to Cd stress.

Effects of Salt Stress on Transcriptional and Physiological Responses in Barley Leaves with Contrasting Salt Tolerance

Salt stress negatively impacts crop production worldwide. Genetic diversity among barley (Hordeum vulgare) landraces adapted to adverse conditions should provide a valuable reservoir of tolerance genes for breeding programs. To identify molecular and biochemical differences between barley genotypes, transcriptomic and antioxidant enzyme profiles along with several morpho-physiological features were compared between salt-tolerant (Boulifa) and salt-sensitive (Testour) genotypes subjected to salt stress. Decreases in biomass, photosynthetic parameters, and relative water content were low in Boulifa compared to Testour. Boulifa had better antioxidant protection against salt stress than Testour, with greater antioxidant enzymes activities including catalase, superoxide dismutase, and guaiacol peroxidase. Transcriptome assembly for both genotypes revealed greater accumulation of differentially expressed transcripts in Testour compared to Boulifa, emphasizing the elevated transcriptional response in Testour following salt exposure. Various salt-responsive genes, including the antioxidant catalase 3, the osmoprotectant betaine aldehyde dehydrogenase 2, and the transcription factors MYB20 and MYB41, were induced only in Boulifa. By contrast, several genes associated with photosystems I and II, and light receptor chlorophylls A and B, were more repressed in Testour. Co-expression network analysis identified specific gene modules correlating with differences in genotypes and morpho-physiological traits. Overall, salinity-induced differential transcript accumulation underlies the differential morpho-physiological response in both genotypes and could be important for breeding salt tolerance in barley.

Overexpression of PnMYB2 from Panax notoginseng induces cellulose and lignin biosynthesis during cell wall formation

Panax notoginseng PnMYB2 is a transcriptional activator of primary and secondary cell wall formation by promoting the PCW-specific gene CesA3 and key lignin biosynthetic gene CCoAOMT1, respectively. R2R3-MYB transcription factors play important roles in regulation secondary cell wall (SCW) formation. However, there are few reports on the functions of MYB transcription factors which involved in both primary cell wall (PCW) and SCW formation. Here, we isolated an R2R3-MYB transcription factor, PnMYB2, from Panax notoginseng roots which are widely used in Chinese traditional medicines and contain abundant cellulose and lignin. The expression pattern of PnMYB2 was similar to the accumulation pattern of cellulose and lignin contents in different organs. PnMYB2 localized in the nucleus and may function as a transcriptional activator. Overexpression of PnMYB2 in Arabidopsis thaliana enhanced cellulose and lignin biosynthesis, and remarkably increased thickness of PCW and SCW in the stem of transgenic plants compared with wild-type plants. The expression levels of genes associated with PCW-specific cellulose synthase (CesA) genes and key SCW-specific lignin biosynthetic genes were significantly increased in PnMYB2-overexpressing plants compared to the wild type plants. Furthermore, yeast one-hybrid, dual-luciferase reporter assays and electrophoretic mobility shift assays (EMSA) results verified that PnMYB2 could bind and activate the promoters of AtCesA3 and PnCesA3, which are the PCW-specific cellulose biosynthetic genes, and AtCCoAOMT1 and PnCCoAOMT1, which are the key lignin biosynthetic genes. These results demonstrated the central role of PnMYB2 in PCW-specific cellulose formation and SCW-specific lignin biosynthesis.

Insights into the Molecular Regulation of Lignin Content in Triploid Poplar Leaves

After polyploidization, plants usually undergo some morphological and physiological changes, including the lignin content of polyploids usually becoming lower than that of diploids. However, the regulatory mechanism of the variation of lignin content in polyploid plants remains unclear. Therefore, in this research, we used full-sib poplar triploids and diploids to explore the molecular regulatory basis of lignin content in poplar triploid leaves through the determination of lignin content, the observation of xylem cells, and transcriptome sequencing. The results showed that the lignin content of triploid leaves was significantly lower than that of diploid leaves. The xylem cells of triploid leaves were significantly larger than those of diploids. Transcriptome sequencing data show that most lignin biosynthesis genes were significantly downregulated, and genes related to cell growth were mostly upregulated in triploid leaves compared with diploid leaves. In addition, co-expression network analysis showed that several transcription factors might be involved in the regulation of lignin biosynthesis. Consequently, the altered expression of genes related to lignin might lead to the reduced lignin content in triploids. These results provide a theoretical basis for further exploring the molecular mechanism of the variation of polyploid lignin content and the utilization of polyploid lignocellulosic resources.

Genome-Wide Association Study Reveals That PvGUX1_1 Is Associated with Pod Stringlessness in Snap Bean (Phaseolus vulgaris L.)

Simple Summary

Using 138 snap bean accessions as plant materials, we investigated their suture strings across two years. With the goal of identifying the gene(s) responsible for the formation of suture strings, we conducted a genome-wide association study. A strong association signal was found in a 266.19 kb region on Chr02. Within the region, 23 candidate genes were identified. Importantly, the sequence and gene expression of PvGUX1_1 differed significantly between sutured pods and non-sutured pods. In addition, PvGUX1_1 was also a domesticated locus that diverged from PvGUX1_2 during an early stage. The results obtained in this study can provide important information for the improvement of pod quality in snap beans.

Abstract

Suture strings are a particularly important pod trait that determine the quality and texture of snap beans (Phaseolus vulgaris L.). The St locus on chromosome 2 has been described as a major locus associated with suture strings. However, the gene and genetic basis underlying this locus remain unknown. Here, we investigated the suture strings of 138 snap bean accessions across two years. A total of 3.66 million single-nucleotide polymorphisms (SNPs) were obtained by deep resequencing. Based on these SNPs, we identified a strong association signal on Chr02 and a promising candidate gene, PvGUX1_1. Further analysis revealed that the 2 bp deletion in the exon of PvGUX1_1 was significantly associated with stringlessness. Comparative mapping indicated that PvGUX1_1 was a domesticated locus and diverged from PvGUX1_2 during an early stage. Our study provides important insights into the genetic mechanism of suture string formation and useful information for snap bean improvement.

Multi-Dimensional Molecular Regulation of Trichome Development in Arabidopsis and Cotton

Plant trichomes are specialized epidermal cells that are widely distributed on plant aerial tissues. The initiation and progression of trichomes are controlled in a coordinated sequence of multiple molecular events. During the past decade, major breakthroughs in the molecular understanding of trichome development were achieved through the characterization of various trichomes defective mutants and trichome-associated genes, which revealed a highly complex molecular regulatory network underlying plant trichome development. This review focuses on the recent millstone in plant trichomes research obtained using genetic and molecular studies, as well as 'omics' analyses in model plant Arabidopsis and fiber crop cotton. In particular, we discuss the latest understanding and insights into the underlying molecular mechanisms of trichomes formation at multiple dimensions, including at the chromatin, transcriptional, post-transcriptional, and post-translational levels. We summarize that the integration of multi-dimensional trichome-associated genes will enable us to systematically understand the molecular regulation network that landscapes the development of the plant trichomes. These advances will enable us to address the unresolved questions regarding the molecular crosstalk that coordinate concurrent and ordered the changes in cotton fiber initiation and progression, together with their possible implications for genetic improvement of cotton fiber.

Transcriptome sequencing and differential expression analysis of natural and BTH-treated wound healing in potato tubers (Solanum tuberosum L.)

Background

Wound healing is a representative phenomenon of potato tubers subjected to mechanical injuries. Our previous results found that benzo-(1,2,3)-thiadiazole-7-carbothioic acid S-methyl ester (BTH) promoted the wound healing of potato tubers. However, the molecular mechanism related to inducible wound healing remains unknown.

Results

Transcriptomic evaluation of healing tissues from potato tubers at three stages, namely, 0 d (nonhealing), 5 d (wounded tubers healed for 5 d) and 5 d (BTH-treated tubers healed for 5 d) using RNA-Seq and differentially expressed genes (DEGs) analysis showed that more than 515 million high-quality reads were generated and a total of 7665 DEGs were enriched, and 16 of these DEGs were selected by qRT-PCR analysis to further confirm the RNA sequencing data. Gene ontology (GO) enrichment analysis indicated that the most highly DEGs were involved in metabolic and cellular processes, and KEGG enrichment analysis indicated that a large number of DEGs were associated with plant hormones, starch and sugar metabolism, fatty acid metabolism, phenylpropanoid biosynthesis and terpenoid skeleton biosynthesis. Furthermore, a few candidate transcription factors, including MYB, NAC and WRKY, and genes related to Ca²⁺-mediated signal transduction were also found to be differentially expressed during wound healing. Most of these enriched DEGs were upregulated after BTH treatment.

Conclusion

This comparative expression profile provided useful resources for studies of the molecular mechanism via these promising candidates involved in natural or elicitor-induced wound healing in potato tubers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08480-1.

The R2R3-MYB transcription factor FaMYB63 participates in regulation of eugenol production in strawberry

The biosynthetic pathway of volatile phenylpropanoids, including 4-allyl-2-methoxyphenol (eugenol), has been investigated in petunia (Petunia hybrida). However, the regulatory network for eugenol accumulation in strawberry (Fragaria × ananassa Duch.) fruit remains unclear. Here, an R2R3-type MYB transcription factor (TF; FaMYB63) was isolated from strawberry by yeast one-hybrid (Y1H) screening using the promoter of the FaEGS1 (eugenol synthase 1 [EGS 1]) gene, which encodes the enzyme responsible for the last step in eugenol biosynthesis. FaMYB63 is phylogenetically distinct from other R2R3-MYB TFs, including FaEOBІІ (EMISSION OF BENZENOID II [EOBII]), which also participates in regulating eugenol biosynthesis in strawberry receptacles. Reverse transcription quantitative PCR (RT-qPCR) assays showed that the expression of FaMYB63 was tissue-specific and consistent with eugenol content through strawberry fruit development, was repressed by abscisic acid, and was activated by auxins (indole-3-acetic acid). Overexpression and RNA interference-mediated silencing of FaMYB63 resulted in marked changes in the transcript levels of the biosynthetic genes FaEGS1, FaEGS2, and FaCAD1 (cinnamyl alcohol dehydrogenase 1 [CAD1]) and, thereby, the accumulation of eugenol. Electrophoretic mobility shift, Y1H, GUS activity, and dual-luciferase activity assays demonstrated that the transcript levels of FaEOBІІ and FaMYB10 were regulated by FaMYB63, but not the other way around. Together, these results demonstrate that FaMYB63 directly activates FaEGS1, FaEGS2, FaCAD1, FaEOBІІ, and FaMYB10 to induce eugenol biosynthesis during strawberry fruit development. These findings deepen the understanding of the regulatory network that influences eugenol metabolism in an edible fruit crop.

An MKP-MAPK protein phosphorylation cascade controls vascular immunity in plants

Global crop production is greatly reduced by vascular diseases. These diseases include bacterial blight of rice and crucifer black rot caused by Xanthomonas oryzae pv. oryzae (Xoo) and Xanthomonas campestris pv. campestris (Xcc). The molecular mechanisms that activate vascular defense against such pathogens remains underexplored. Here, we show that an Arabidopsis MAPK phosphatase 1 (MKP1) mutant has increased host susceptibility to the adapted pathogen Xcc and is compromised in nonhost resistance to the rice pathogen Xoo. MKP1 regulates MAPK-mediated phosphorylation of the transcription factor MYB4 that negatively regulates vascular lignification through inhibiting lignin biosynthesis. Induction of lignin biosynthesis is, therefore, an important part of vascular-specific immunity. The role of MKP-MAPK-MYB signaling in lignin biosynthesis and vascular resistance to Xoo is conserved in rice, indicating that these factors form a tissue-specific defense regulatory network. Our study likely reveals a major vascular immune mechanism that underlies tissue-specific disease resistance against bacterial pathogens in plants.

A Genome-Wide Association Study Dissects the Genetic Architecture of the Metaxylem Vessel Number in Maize Brace Roots

Metaxylem vessels in maize brace roots are key tissue, and their number (MVN) affects plant water and inorganic salt transportation and lodging resistance. Dissecting the genetic basis of MVN in maize brace roots can help guide the genetic improvement of maize drought resistance and lodging resistance during late developmental stages. In this study, we used 508 inbred lines with tropical, subtropical, and temperate backgrounds to analyze the genetic architecture of MVN in maize brace roots. The phenotypic variation in MVN in brace roots was evaluated in three environments, which revealed broad natural variation and relative low levels of heritability (h ² = 0.42). Stiff-stalk lines with a temperate background tended to have higher MVNs than plants in other genetic backgrounds. MVN was significantly positively correlated with plant height, tassel maximum axis length, ear length, and kernel number per row, which indicates that MVN may affect plant morphological development and yield. In addition, MVN was extremely significantly negatively correlated with brace root radius, but significantly positively correlated with brace root angle (BRA), diameter, and number, thus suggesting that the morphological function of some brace root traits may be essentially determined by MVN. Association analysis of MVN in brace roots combined 1,253,814 single nucleotide polymorphisms (SNPs) using FarmCPU revealed a total of nine SNPs significantly associated with MVN at P < 7.96 × 10^-7. Five candidate genes for MVN that may participate in secondary wall formation (GRMZM2G168365, GRMZM2G470499, and GRMZM2G028982) and regulate flowering time (GRMZM2G381691 and GRMZM2G449165). These results provide useful information for understanding the genetic basis of MVN in brace root development. Further functional studies of identified candidate genes should help elucidate the molecular pathways that regulate MVN in maize brace roots.

ODORANT1 targets multiple metabolic networks in petunia flowers

Scent bouquets produced by the flowers of Petunia spp. (petunia) are composed of a complex mixture of floral volatile benzenoid and phenylpropanoid compounds (FVBPs), which are specialized metabolites derived from phenylalanine (Phe) through an interconnected network of enzymes. The biosynthesis and emission of high levels of these volatiles requires coordinated transcriptional activation of both primary and specialized metabolic networks. The petunia R2R3-MYB transcription factor ODORANT 1 (ODO1) was identified as a master regulator of FVBP production and emission; however, our knowledge of the direct regulatory targets of ODO1 has remained limited. Using chromatin immunoprecipitation followed by sequencing (ChIP-seq) in petunia flowers, we identify genome-wide ODO1-bound genes that are enriched not only in genes involved in the biosynthesis of the Phe precursor, as previously reported, but also genes associated with the specialized metabolic pathways involved in generating phenylpropanoid intermediates for FVBPs. ODO1-bound genes are also involved in methionine and S-adenosylmethionine metabolism, which could modulate methyl group supplies for certain FVBPs. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) and RNA-seq analysis in an ODO1 RNAi knockdown line revealed that ODO1-bound targets are expressed at lower levels when ODO1 is suppressed. A cis-regulatory motif, CACCAACCCC, was identified as a potential binding site for ODO1 in the promoters of genes that are both bound and activated by ODO1, which was validated by in planta promoter reporter assays with wild-type and mutated promoters. Overall, our work presents a mechanistic model for ODO1 controlling an extensive gene regulatory network that contributes to FVBP production to give rise to floral scent.

PtrHAT22, as a higher hierarchy regulator, coordinately regulates secondary cell wall component biosynthesis in Populus trichocarpa

Homeodomain-leucine zipper (HD-Zip) II transcription factors (TFs) have been reported to play vital roles in diverse biological processes of plants. However, it remains unclear whether HD-Zip II TFs regulate secondary cell wall (SCW) in woody plants. In this study, we performed the functional characterization of a Populus trichocarpa HD-Zip II TF, PtrHAT22, which encodes a nuclear localized transcription repressor predominantly expressing in secondary developing tissues. Overexpression of PtrHAT22 showed arrested growths, including reduced heights and diameters above the ground, small leaves, and decreased biomass. Meanwhile, the contents of lignin, cellulose, and thickness of SCW significantly decreased, whilst the content of hemicellulose obviously increased in PtrHAT22 transgenic poplar. The expressions of some wood-associated TFs and structural genes significantly changed accordingly with the alternations of SCW characteristics in PtrHAT22 transgenic poplar. Furthermore, PtrHAT22 directly repressed the promoter activities of PtrMYB20, PtrMYB28, and PtrCOMT2, and bind two cis-acting elements that were specifically enriched in their promoter regions. Taken together, our results suggested that PtrHAT22, as a higher hierarchy TF like PtrWNDs, exerted coordination regulation of poplar SCW component biosynthesis through directly and indirectly regulating structural genes and different hierarchy TFs of SCW formation network.

Transcriptomic Analysis Reveal the Molecular Mechanisms of Seed Coat Development in Cucurbita pepo L

Cucurbita pepo is one of the earliest cultivated crops. It is native to Central and South America and is now widely cultivated all over the world for its rich nutrition, short growth period, and high yield, which make it suitable for intercropping. Hull-less C. pepo L. (HLCP) is a rare variant in nature that is easier to consume. Its seed has a seed kernel but lacks a seed coat. The molecular mechanism underlying the lack of seed coat development in the HLCP variety is not clear yet. The BGISEQ-500 sequencing platform was used to sequence 18 cDNA libraries of seed coats from hulled C. pepo (CP) and HLCP at three developmental stages (8, 18, and 28 days) post-pollination. We found that lignin accumulation in the seed coat of the HLCP variety was much lower than that of the CP variety. A total of 2,099 DEGs were identified in the CP variety, which were enriched mainly in the phenylpropanoid biosynthesis pathway, amino sugar, and nucleotide sugar metabolism pathways. A total of 1,831 DEGs were identified in the HLCP variety and found to be enriched mainly in the phenylpropanoid biosynthesis and metabolism pathways of starch and sucrose. Among the DEGs, hub proteins (FusA), protein kinases (IRAK4), and several transcription factors related to seed coat development (MYB, bHLH, NAC, AP2/EREBP, WRKY) were upregulated in the CP variety. The relative expression levels of 12 randomly selected DEGs were determined using quantitative real-time PCR analysis and found to be consistent with those obtained using RNA-Seq, with a correlation coefficient of 0.9474. We found that IRAK4 protein kinases, AP2/EREBP, MYB, bHLH, and NAC transcription factors may play important roles in seed coat development, leading to the formation of HLCP.

Comparative transcriptomics analysis of contrasting varieties of Eucalyptus camaldulensis reveals wind resistance genes

BACKGROUND

Wind, an important abiotic stress factor, affects forests in coastal areas, causes tree damage and timber loss.

METHODS

Two genotypes of Eucalyptus camaldulensis-strong wind-resistant CA5 and weak wind-resistant C037 were used for RNA-seq analysis to screen for candidate wind-resistance genes and transcription factors (TFs) by comparing the transcriptome analysis of the two varieties in response to wind stress.

RESULTS

It showed that 7061 differentially expressed unigenes could be annotated including 4,110 up-regulated unigenes and 2,951 down-regulated unigenes. Gene Ontology (GO) analysis revealed that six cellulose pathways were involved in response to wind stress. The unigenes in phenylpropanoid biosynthesis, phenylalanine metabolism, and flavonoid biosynthesis pathways were found to be differentially expressed based on Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Moreover, 37 differentially expressed genes were functionally annotated to be involved in the secondary metabolism of phenylalanine (ko00940). Seventy-eight TFs related to the regulating cellulose and lignin synthesis were expressed differently from the various treatments. The expressions of C3H, POX, MYB, NAC, Gene008307, and Gene011799 were significantly upregulated in CA5. Overall, the main response of Eucalyptus to wind stress was associated with cell wall biosynthesis; key genes of cellulose and lignin biosynthesis pathways and related TFs were involved in the tree response to wind stress.

Genome-Wide Identification of MYB Transcription Factors and Screening of Members Involved in Stress Response in Actinidia

MYB transcription factors (TFs) play an active role in plant responses to abiotic stresses, but they have not been systematically studied in kiwifruit (Actinidia chinensis). In this study, 181 AcMYB TFs were identified from the kiwifruit genome, unevenly distributed on 29 chromosomes. The high proportion (97.53%) of segmental duplication events (Ka/Ks values less than 1) indicated that AcMYB TFs underwent strong purification selection during evolution. According to the conservative structure, 91 AcR2R3-MYB TFs could be divided into 34 subgroups. A combination of transcriptomic data under drought and high temperature from four AcMYB TFs (AcMYB2, AcMYB60, AcMYB61 and AcMYB102) was screened out in response to stress and involvement in the phenylpropanoid pathway. They were highly correlated with the expression of genes related to lignin biosynthesis. qRT-PCR analysis showed that they were highly correlated with the expression of genes related to lignin biosynthesis in different tissues or under stress, which was consistent with the results of lignin fluorescence detection. The above results laid a foundation for further clarifying the role of MYB in stress.