SND2, a NAC transcription factor gene, regulates genes involved in secondary cell wall development in Arabidopsis fibres and increases fibre cell area in Eucalyptus

C2H2 zinc finger protein PagIDD15A regulates secondary wall thickening and lignin biosynthesis in poplar

Wood production is largely determined by the activity of cambial cell proliferation, and the secondary cell wall (SCW) thickening of xylem cells determines the wood property. In this study, we identified an INDETERMINATE DOMAIN (IDD) type C2H2 zinc finger transcription factor PagIDD15A as a regulator of wood formation in Populus alba × Populus glandulosa. Downregulation of PagIDD15A expression by RNA interference (RNAi) inhibited xylem development and xylem cell secondary wall thickening. RNA-seq analysis showed that PagPAL1, PagCCR2 and PagCCoAOMT1 were downregulated in the differentiating xylem of the PagIDD15A-RNAi transgenic plants, showing that PagIDD15A may regulate SCW biosynthesis through inhibiting lignin biosynthesis. The downregulation of PagVND6-B2, PagMYB10 and PagMYC4 and upregulation of PagWRKY12 in the differentiating xylem of RNAi transgenic plants suggest that PagIDD15A may also regulate these transcription factor (TF) genes to affect SCW thickening. RT-qPCR analysis in the phloem-cambium of RNAi transgenic demonstrates that PagIDD15A may regulate the expression of the genes associated with cell proliferation, including, PagSHR (SHORTROOT), PagSCR (SCARECROW), PagCYCD3;1 (CYCLIN D3;1) and PagSMR4 (SIAMESE-RELATED4), to affect the cambial activity. This study provides the knowledge of the IDD-type C2H2 zinc finger protein in regulating wood formation.

Single-cell transcriptome atlases of soybean root and mature nodule reveal new regulatory programs that control the nodulation process

The soybean root system is complex. In addition to being composed of various cell types, the soybean root system includes the primary root, the lateral roots, and the nodule, an organ in which mutualistic symbiosis with N-fixing rhizobia occurs. A mature soybean root nodule is characterized by a central infection zone where atmospheric nitrogen is fixed and assimilated by the symbiont, resulting from the close cooperation between the plant cell and the bacteria. To date, the transcriptome of individual cells isolated from developing soybean nodules has been established, but the transcriptomic signatures of cells from the mature soybean nodule have not yet been characterized. Using single-nucleus RNA-seq and Molecular Cartography technologies, we precisely characterized the transcriptomic signature of soybean root and mature nodule cell types and revealed the co-existence of different sub-populations of B. diazoefficiens-infected cells in the mature soybean nodule, including those actively involved in nitrogen fixation and those engaged in senescence. Mining of the single-cell-resolution nodule transcriptome atlas and the associated gene co-expression network confirmed the role of known nodulation-related genes and identified new genes that control the nodulation process. For instance, we functionally characterized the role of GmFWL3, a plasma membrane microdomain-associated protein that controls rhizobial infection. Our study reveals the unique cellular complexity of the mature soybean nodule and helps redefine the concept of cell types when considering the infection zone of the soybean nodule.

The transcriptional regulation of a putative hemicellulose gene, PtrPARVUS2 in poplar

The plant cell wall serves as a critical interface between the plant and its environment, offering protection against various stresses and contributing to biomass production. Hemicellulose is one of the major components of the cell wall, and understanding the transcriptional regulation of its production is essential to fully understanding cell wall formation. This study explores the regulatory mechanisms underlying one of the genes involved in hemicellulose biosynthesis, PtrPARVUS2. Six transcription factors (TFs) were identified from a xylem-biased library to negatively regulate PtrPARVUS2 expression. These TFs, belonging to diverse TF families, were confirmed to bind to specific cis-elements in the PtrPARVUS2 promoter region, as validated by Yeast One-Hybrid (Y1H) assays, transient expression analysis, and Chromatin Immunoprecipitation sequencing (ChIP-seq) assays. Furthermore, motif analysis identified putative cis-regulatory elements bound by these TFs, shedding light on the transcriptional regulation of SCW biosynthesis genes. Notably, several TFs targeted genes encoding uridine diphosphate glycosyltransferases (UGTs), crucial enzymes involved in hemicellulose glycosylation. Phylogenetic analysis of UGTs regulated by these TFs highlighted their diverse roles in modulating hemicellulose synthesis. Overall, this study identifies a set of TFs that regulate PARVUS2 in poplar, providing insights into the intricate coordination of TFs and PtrPARVUS2 in SCW formation. Understanding these regulatory mechanisms enhances our ability to engineer plant biomass for tailored applications, including biofuel production and bioproduct development.

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.

NACs strike again: NOR-like1 is responsible for cuticle development in tomato fruit

This article comments on:

Liu G-S, Huang H, Grierson D, Gao Y, Ji X, Peng Z-Z, Li H-L, Niu X-L, Jia W, He J-L, Xiang L-T, Gao H-Y, Qu G-Q, Zhu H-L, Zhu B-Z, Luo Y-B, Fu D-Q. 2024. NAC transcription factor SlNOR-like1 plays a dual regulatory role in tomato fruit cuticle formation. Journal of Experimental Botany 75, 1903–1918.

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.

Trees need closure too: Wound-induced secondary vascular tissue regeneration

Trees play a pivotal role in terrestrial ecosystems as well as being an important natural resource. These attributes are primarily associated with the capacity of trees to continuously produce woody tissue from the vascular cambium, a ring of stem cells located just beneath the bark. Long-lived trees are exposed to a myriad of biological and environmental stresses that may result in wounding, leading to a loss of bark and the underlying vascular cambium. This affects both wood formation and the quality of timber arising from the tree. In addition, the exposed wound site is a potential entry point for pathogens that cause disease. In response to wounding, trees have the capacity to regenerate lost or damaged tissues at this site. Investigating gene expression changes associated with different stages of wound healing reveals complex and dynamic changes in the activity of transcription factors, signalling pathways and hormone responses. In this review we summarise these data and discuss how they relate to our current understanding of vascular cambium formation and xylem differentiation during secondary growth. Based on this analysis, a model for wound healing that provides the conceptual foundations for future studies aimed at understanding this intriguing process is proposed.

GhWRKY41 forms a positive feedback regulation loop and increases cotton defence response against Verticillium dahliae by regulating phenylpropanoid metabolism

Despite the established significance of WRKY proteins and phenylpropanoid metabolism in plant immunity, how WRKY proteins modulate aspects of the phenylpropanoid pathway remains undetermined. To understand better the role of WRKY proteins in plant defence, we identified a cotton (Gossypium hirsutum) protein, GhWRKY41, that is, universally and rapidly induced in three disease-resistant cotton cultivars following inoculation with the plant pathogenic fungus, Verticillium dahliae. We show that overexpression of GhWRKY41 in transgenic cotton and Arabidopsis enhances resistance to V. dahliae, while knock-down increases cotton more susceptibility to the fungus. GhWRKY41 physically interacts with itself and directly activates its own transcription. A genome-wide chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq), in combination with RNA sequencing (RNA-seq) analyses, revealed that 43.1% of GhWRKY41-binding genes were up-regulated in cotton upon inoculation with V. dahliae, including several phenylpropanoid metabolism master switches, receptor kinases, and disease resistance-related proteins. We also show that GhWRKY41 homodimer directly activates the expression of GhC4H and Gh4CL, thereby modulating the accumulation of lignin and flavonoids. This finding expands our understanding of WRKY-WRKY protein interactions and provides important insights into the regulation of the phenylpropanoid pathway in plant immune responses by a WRKY protein.

High-resolution anatomical and spatial transcriptome analyses reveal two types of meristematic cell pools within the secondary vascular tissue of poplar stem

The secondary vascular tissue emanating from meristems is central to understanding how vascular plants such as forest trees evolve, grow, and regulate secondary radial growth. However, the overall molecular characterization of meristem origins and developmental trajectories from primary to secondary vascular tissues in woody tree stems is technically challenging. In this study, we combined high-resolution anatomic analysis with a spatial transcriptome (ST) technique to define features of meristematic cells in a developmental gradient from primary to secondary vascular tissues in poplar stems. The tissue-specific gene expression of meristems and derived vascular tissue types were accordingly mapped to specific anatomical domains. Pseudotime analyses were used to track the origins and changes of meristems throughout the development from primary to secondary vascular tissues. Surprisingly, two types of meristematic-like cell pools within secondary vascular tissues were inferred based on high-resolution microscopy combined with ST, and the results were confirmed by in situ hybridization of, transgenic trees, and single-cell sequencing. The rectangle shape procambium-like (PCL) cells develop from procambium meristematic cells and are located within the phloem domain to produce phloem cells, whereas fusiform shape cambium zone (CZ) meristematic cells develop from fusiform metacambium meristematic cells and are located inside the CZ to produce xylem cells. The gene expression atlas and transcriptional networks spanning the primary transition to secondary vascular tissues generated in this work provide new resources for studying the regulation of meristem activities and the evolution of vascular plants. A web server (https://pgx.zju.edu.cn/stRNAPal/) was also established to facilitate the use of ST RNA-seq data.

The NAC transcription factor family in Eucommia ulmoides: Genome-wide identification, characterization, and network analysis in relation to the rubber biosynthetic genes

The NAC transcription factor family is a large plant gene family, participating in plant growth and development, secondary metabolite synthesis, biotic and abiotic stresses responses, and hormone signaling. Eucommia ulmoides is a widely planted economic tree species in China that can produce trans-polyisoprene: Eucommia rubber (Eu-rubber). However, genome-wide identification of the NAC gene family has not been reported in E. ulmoides. In this study, 71 NAC proteins were identified based on genomic database of E. ulmoides. Phylogenetic analysis showed that the EuNAC proteins were distributed in 17 subgroups based on homology with NAC proteins in Arabidopsis, including the E. ulmoides-specific subgroup Eu_NAC. Gene structure analysis suggested that the number of exons varied from 1 to 7, and multitudinous EuNAC genes contained two or three exons. Chromosomal location analysis revealed that the EuNAC genes were unevenly distributed on 16 chromosomes. Three pairs of genes of tandem duplicates genes and 12 segmental duplications were detected, which indicated that segmental duplications may provide the primary driving force of expansion of EuNAC. Prediction of cis-regulatory elements indicated that the EuNAC genes were involved in development, light response, stress response and hormone response. For the gene expression analysis, the expression levels of EuNAC genes in various tissues were quite different. To explore the effect of EuNAC genes on Eu-rubber biosynthesis, a co-expression regulatory network between Eu-rubber biosynthesis genes and EuNAC genes was constructed, which indicated that six EuNAC genes may play an important role in the regulation of Eu-rubber biosynthesis. In addition, this six EuNAC genes expression profiles in E. ulmoides different tissues were consistent with the trend in Eu-rubber content. Quantitative real-time PCR analysis showed that EuNAC genes were responsive to different hormone treatment. These results will provide a useful reference for further studies addressing the functional characteristics of the NAC genes and its potential role in Eu-rubber biosynthesis.

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.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression analysis under Fusarium kyushuense and drought stress conditions in Passiflora edulis

The NAC gene family is one of the largest plant transcription factors (TFs) families and plays important roles in plant growth, development, metabolism, and biotic and abiotic stresses. However, NAC gene family has not been reported in passion fruit (Passiflora edulis). In this study, a total of 105 NAC genes were identified in the passion fruit genome and were unevenly distributed across all nine-passion fruit chromomere, with a maximum of 48 PeNAC genes on chromosome one. The physicochemical features of all 105 PeNAC genes varied including 120 to 3,052 amino acids, 3 to 8 conserved motifs, and 1 to 3 introns. The PeNAC genes were named (PeNAC001-PeNAC105) according to their chromosomal locations and phylogenetically grouped into 15 clades (NAC-a to NAC-o). Most PeNAC proteins were predicted to be localized in the nucleus. The cis-element analysis indicated the possible roles of PeNAC genes in plant growth, development, light, hormones, and stress responsiveness. Moreover, the PeNAC gene duplications including tandem (11 gene pairs) and segmental (12 gene pairs) were identified and subjected to purifying selection. All PeNAC proteins exhibited similar 3D structures, and a protein-protein interaction network analysis with known Arabidopsis proteins was predicted. Furthermore, 17 putative ped-miRNAs were identified to target 25 PeNAC genes. Potential TFs including ERF, BBR-BPC, Dof, and bZIP were identified in promoter region of all 105 PeNAC genes and visualized in a TF regulatory network. GO and KEGG annotation analysis exposed that PeNAC genes were related to different biological, molecular, and cellular terms. The qRT-PCR expression analysis discovered that most of the PeNAC genes including PeNAC001, PeNAC003, PeNAC008, PeNAC028, PeNAC033, PeNAC058, PeNAC063, and PeNAC077 were significantly upregulated under Fusarium kyushuense and drought stress conditions compared to controls. In conclusion, these findings lay the foundation for further functional studies of PeNAC genes to facilitate the genetic improvement of plants to stress resistance.

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.

Genomic Regions Associated With Seed Meal Quality Traits in Brassica napus Germplasm

The defatted Brassica napus (rapeseed) meal can be high-protein feed for livestock as the protein value of rapeseed meal is higher than that of the majority of other vegetable proteins. Extensive work has already been carried out on developing canola rapeseed where the focus was on reducing erucic acid and glucosinolate content, with less consideration to other antinutritional factors such as tannin, phytate, sinapine, crude fiber, etc. The presence of these antinutrients limits the use and marketing of rapeseed meals and a significant amount of it goes unused and ends up as waste. We investigated the genetic architecture of crude protein, methionine, tryptophan, total phenols, β-carotene, glucosinolates (GLSs), phytate, tannins, sinapine, and crude fiber content of defatted seed meal samples by conducting a genome-wide association study (GWAS), using a diversity panel comprising 96 B. napus genotypes. Genotyping by sequencing was used to identify 77,889 SNPs, spread over 19 chromosomes. Genetic diversity and phenotypic variations were generally high for the studied traits. A total of eleven genotypes were identified which showed high-quality protein, high antioxidants, and lower amount of antinutrients. A significant negative correlation between protein and limiting amino acids and a significant positive correlation between GLS and phytic acid were observed. General and mixed linear models were used to estimate the association between the SNP markers and the seed quality traits and quantile-quantile (QQ) plots were generated to allow the best-fit algorithm. Annotation of genomic regions around associated SNPs helped to predict various trait-related candidates such as ASP2 and EMB1027 (amino acid biosynthesis); HEMA2, GLU1, and PGM (tryptophan biosynthesis); MS3, CYSD1, and MTO1 (methionine biosynthesis); LYC (β-carotene biosynthesis); HDR and ISPF (MEP pathway); COS1 (riboflavin synthesis); UGT (phenolics biosynthesis); NAC073 (cellulose and hemicellulose biosynthesis); CYT1 (cellulose biosynthesis); BGLU45 and BGLU46 (lignin biosynthesis); SOT12 and UGT88A1 (flavonoid pathway); and CYP79A2, DIN2, and GSTT2 (GLS metabolism), etc. The functional validation of these candidate genes could confirm key seed meal quality genes for germplasm enhancement programs directed at improving protein quality and reducing the antinutritional components in B. napus.

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.

Identification and Functional Analysis of the CgNAC043 Gene Involved in Lignin Synthesis from Citrusgrandis "San Hong"

Overaccumulation of lignin (a physiological disorder known as granulation) often occurs during fruit ripening and postharvest storage in pomelo (Citrus grandis). It causes an unpleasant fruit texture and taste. Previous studies have shown that lignin metabolism is closely associated with the process of juice sacs granulation. At present, the underlying transcriptional regulatory mechanisms remain unclear. In this study, we identified and isolated a candidate NAC transcription factor, CgNAC043, that is involved in the regulation of lignin biosynthesis in Citrus grandis, which has homologs in Arabidopsis and other plants. We used the fruit juice sacs of 'San hong' as the material, the staining for lignin with HCl-phloroglucinol of fruit juice sacs became dark red from the various developmental stages at 172 to 212 days post anthesis (DPA). The RT-qPCR was used to analyze the gene expression of CgNAC043 and its target gene CgMYB46 in fruit sacs, it was found that the expression trend of CgNAC043 was basically same as CgMYB46, which increased gradually and peaked at 212 DPA. The expression level of CgNAC043 in juice sacs obtained away from the core was the lowest, while those near the core and granulated area were highly expressed. The transcriptional activation activity of CgNAC043 and CgMYB46 was analyzed by a yeast two-hybrid system, with only CgNAC043 showing transcriptional activation activity in Y2H Gold yeast. A transformation vector, p1301- CgNAC043, was transformed into the mesocarp of 'San hong' by Agrobacterium-mediated transformation. Results showed that the expression of transcription factors CgMYB58 and CgMYB46 are all upregulated. Further experiments proved that CgNAC043 not only can directly trans-activate the promoter of CgMYB46 but also trans-activate the promoters for the lignin biosynthesis-related genes CgCCoAOMT and CgC3H by dual luciferase assay. We isolated the CgNAC043 gene in pomelo and found CgNAC043 regulates target genes conferring the regulation of juice sacs granulation.

High-density genetic linkage mapping reveals low stability of QTLs across environments for economic traits in Eucalyptus

Introduction

Eucalyptus urophylla, E. tereticornis and their hybrids are the most important commercial forest tree species in South China where they are grown for pulpwood and solid wood production. Construction of a fine-scale genetic linkage map and detecting quantitative trait loci (QTL) for economically important traits linked to these end-uses will facilitate identification of the main candidate genes and elucidate the regulatory mechanisms.

Method

A high-density consensus map (a total of 2754 SNPs with 1359.18 cM) was constructed using genotyping by sequencing (GBS) on clonal progenies of E. urophylla × tereticornis hybrids. QTL mapping of growth and wood property traits were conducted in three common garden experiments, resulting in a total of 108 QTLs. A total of 1052 candidate genes were screened by the efficient combination of QTL mapping and transcriptome analysis.

Results

Only ten QTLs were found to be stable across two environments, and only one (qSG10Stable mapped on chromosome 10, and associated with lignin syringyl-to-guaiacyl ratio) was stable across all three environments. Compared to other QTLs, qSG10Stable explained a very high level of phenotypic variation (18.4-23.6%), perhaps suggesting that QTLs with strong effects may be more stably inherited across multiple environments. Screened candidate genes were associated with some transcription factor families, such as TALE, which play an important role in the secondary growth of plant cell walls and the regulation of wood formation.

Discussion

While QTLs such as qSG10Stable, found to be stable across three sites, appear to be comparatively uncommon, their identification is likely to be a key to practical QTL-based breeding. Further research involving clonally-replicated populations, deployed across multiple target planting sites, will be required to further elucidate QTL-by-environment interactions.

Deep transcriptomic study reveals the role of cell wall biosynthesis and organization networks in the developing shell of peanut pod

BACKGROUND

Peanut (Arachis hypogaea L.) belongs to an exceptional group of legume plants, wherein the flowers are produced aerially, but the pods develop under the ground. In such a unique environment, the pod's outer shell plays a vital role as a barrier against mechanical damage and soilborne pathogens. Recent studies have reported the uniqueness and importance of gene expression patterns that accompany peanut pods' biogenesis. These studies focused on biogenesis and pod development during the early stages, but the late developmental stages and disease resistance aspects still have gaps. To extend this information, we analyzed the transcriptome generated from four pod developmental stages of two genotypes, Hanoch (Virginia-type) and IGC53 (Peruvian-type), which differs significantly in their pod shell characteristics and pathogen resistance.

RESULTS

The transcriptome study revealed a significant reprogramming of the number and nature of differentially expressed (DE) genes during shell development. Generally, the numbers of DE genes were higher in IGC53 than in Hanoch, and the R5-R6 transition was the most dynamic in terms of transcriptomic changes. Genes related to cell wall biosynthesis, modification and transcription factors (TFs) dominated these changes therefore, we focused on their differential, temporal and spatial expression patterns. Analysis of the cellulose synthase superfamily identified specific Cellulose synthase (CesAs) and Cellulose synthase-like (Csl) genes and their coordinated interplay with other cell wall-related genes during the peanut shell development was demonstrated. TFs were also identified as being involved in the shell development process, and their pattern of expression differed in the two peanut genotypes. The shell component analysis showed that overall crude fiber, cellulose, lignin, hemicelluloses and dry matter increased with shell development, whereas K, N, protein, and ash content decreased. Genotype IGC53 contained a higher level of crude fiber, cellulose, NDF, ADF, K, ash, and dry matter percentage, while Hanoch had higher protein and nitrogen content.

CONCLUSIONS

The comparative transcriptome analysis identified differentially expressed genes, enriched processes, and molecular processes like cell wall biosynthesis/modifications, carbohydrate metabolic process, signaling, transcription factors, transport, stress, and lignin biosynthesis during the peanut shell development between two contrasting genotypes. TFs and other genes like chitinases were also enriched in peanut shells known for pathogen resistance against soilborne major pathogens causing pod wart disease and pod damages. This study will shed new light on the biological processes involved with underground pod development in an important legume crop.

Investigation Into Different Wood Formation Mechanisms Between Angiosperm and Gymnosperm Tree Species at the Transcriptional and Post-transcriptional Level

Enormous distinctions of the stem structure and cell types between gymnosperms and angiosperms tree species are expected to cause quite different wood physical and mechanical attributes, however, the molecular mechanisms underlying the differing wood morphology are still unclear. In this study, we compared the transcriptomes obtained by RNA-Seq between Populus alba × P. glandulosa clone 84K, and Larix kaempferi (Lamb.) Carr trees. Available genome resource served as reference for P. alba × P. glandulosa and the Iso-Seq results of a three-tissues mixture (xylem, phloem, and leaf) were used as the reference for L. kaempferi to compare the xylem-specifically expressed genes and their alternative splicing model. Through screening, we obtained 13,907 xylem-specifically expressed genes (5,954 up-regulated, 7,953 down-regulated) in the xylem of P. alba × P. glandulosa, and 2,596 xylem-specifically expressed genes (1,648 up-regulated, 948 down-regulated) in the xylem of L. kaempferi. From the GO and KEGG analyses, some genes associated with two wood formation-related pathways, namely those for phenylpropanoid biosynthesis, and starch and sucrose metabolism, were successfully screened. Then the distributions and gene expression models between P. alba × P. glandulosa and L. kaempferi in those pathways were compared, which suggested differential wood formation processes between the angiosperm and gymnosperm trees. Furthermore, a Weight Gene Co-expression Network Analysis (WGCNA) for total xylem-specifically expressed genes in two species was conducted, from which wood formation-related modules were selected to build a co-expression network for the two tree species. The genes within this co-expression network showed different co-expression relationships between the angiosperm and gymnosperm woody species. Comparing the alternative splicing events for wood formation-related genes suggests a different post-transcriptional regulation process exists between the angiosperm and gymnosperm trees. Our research thus provides the foundation for the in-depth investigation of different wood formation mechanisms of angiosperm and gymnosperm species.

Comparative anatomical and transcriptomic analyses of the color variation of leaves in Aquilaria sinensis

Color variation in plant tissues is a common phenomenon accompanied with a series of biological changes. In this study, a special-phenotype Aquilaria sinensis (GS) with color variation of leaf was firstly reported, and DNA barcode sequences showed GS samples could not be discriminated clearly with the normal A. sinensis sample (NS), which suggested that the variety was not the cause of the GS formation. To reveal the characteristics of GS compared to NS, the anatomical and transcriptome sequencing studies were carried out. In microscopic observation, the leaves of golden-vein-leaf sample (LGS) and normal-vein-leaf sample (LNS) showed significant differences including the area of the included phloem in midrib and the thickness parameters of palisade and spongy tissues; the stems of golden-vein-leaf sample (SGS) and normal-vein-leaf sample (SNS) were also different in many aspects such as the area of vessels and included phloem. In addition, the structure of chloroplast was more complete in the midrib of LNS than that of LGS, and some particles suspected as virus were found through transmission electron microscope as well. Genes upregulated in LGS in contrast with LNS were mainly enriched in photosynthesis. As for stems, most of the genes upregulated in SGS compared to SNS were involved in translation and metabolism processes. The pathways about photosynthesis and chlorophyll metabolism as well as some important transcription factors may explain the molecular mechanism of the unique phenotypes of leaves and the genes related to suberin biosynthesis may result in the difference of stems. In addition, the genes about defense response especially biotic stress associated with numerous pathogenesis-related (PR) genes upregulated in LGS compared to LNS indicated that the pathogen may be the internal factor. Taken together, our results reveal the macro- and micro-phenotype variations as well as gene expression profiles between GS and NS, which could provide valuable clues for elucidating the mechanism of the color variation of Aquilaria.

The biotechnological importance of the plant-specific NAC transcription factor family in crop improvement

Climate change, malnutrition, and food insecurity are the inevitable challenges being faced by the agriculture sector today. Plants are susceptible to extreme temperatures during the crucial phases of flowering and seed development, and elevated carbon levels also lead to yield losses. Productivity is also affected by floods and droughts. Therefore, increasing plant yield and stress tolerance are the priorities to be met through novel biotechnological interventions. The contributions of NAC genes towards enhancing plant survivability under stress is well known. Here we focus on the potential of NAC genes in the regulation of abiotic stress tolerance, secondary cell wall synthesis, lateral root development, yield potential, seed size and biomass, ROS signaling, leaf senescence, and programmed cell death. Once naturally tolerant candidate NAC genes have been identified, and the nature of their association with growth and fitness against multi-environmental stresses has been determined, they can be exploited for building inherent tolerance in future crops via transgenic technologies. An update on the latest developments is provided in this review, which summarizes the current understanding of the roles of NAC in the establishment of various stress-adaptive mechanisms in model and food crop plants.

Comparative Analyses of Full-Length Transcriptomes Reveal Gnetum luofuense Stem Developmental Dynamics

Genus Gnetum, of which the majority species are pantropical liana, have broad industrial uses including for string, nets, and paper production. Although numerous studies have investigated anatomical structures during stem development, the underlying molecular mechanisms that regulate this developmental trajectory in Gnetum species remain poorly understood. A total of 12 full-length transcriptomes were generated from four stem developmental stages of an arborescent representative of this genus, Gnetum luofuense, using Oxford Nanopore Technologies. The results of this analysis reveal a total of 24,151 alternative splicing (AS) and 134,391 alternative polyadenylation events. A remarkably dynamic pattern of AS events, especially in the case of intron retentions, was found across the four developmental stages while no dynamic pattern was found among transcript numbers with varied poly(A) sites. A total of 728 long non-coding RNAs were also detected; the number of cis-regulated target genes dramatically increased while no changes were found among trans-regulated target genes. In addition, a K-means clustering analysis of all full-length transcripts revealed that primary growth is associated with carbohydrate metabolism and fungi defense, while secondary growth is closely linked with photosynthesis, nitrogen transportation, and leaf ontogenesis. The use of weighted gene co-expression network analysis as well as differentially expressed transcripts reveals that bHLH, GRF, and MYB-related transcription factors are involved in primary growth, while AP2/ERF, MYB, NAC, PLAZ, and bZIP participate in G. luofuense stem secondary growth. The results of this study provide further evidence that Nanopore sequencing technology provides a cost-effective method for generating full-length transcriptome data as well as for investigating seed plant organ development.

Proximal and Distal Parts of Sweetpotato Adventitious Roots Display Differences in Root Architecture, Lignin, and Starch Metabolism and Their Developmental Fates

Sweetpotato is an important food crop globally, serving as a rich source of carbohydrates, vitamins, fiber, and micronutrients. Sweetpotato yield depends on the modification of adventitious roots into storage roots. The underlying mechanism of this developmental switch is not fully understood. Interestingly, storage-root formation is manifested by formation of starch-accumulating parenchyma cells and bulking of the distal part of the root, while the proximal part does not show bulking. This system, where two parts of the same adventitious root display different developmental fates, was used by us in order to better characterize the anatomical, physiological, and molecular mechanisms involved in sweetpotato storage-root formation. We show that, as early as 1 and 2 weeks after planting, the proximal part of the root exhibited enhanced xylem development together with increased/massive lignin deposition, while, at the same time, the distal root part exhibited significantly elevated starch accumulation. In accordance with these developmental differences, the proximal root part exhibited up-regulated transcript levels of sweetpotato orthologs of Arabidopsis vascular-development regulators and key genes of lignin biosynthesis, while the distal part showed up-regulation of genes encoding enzymes of starch biosynthesis. All these recorded differences between proximal and distal root parts were further enhanced at 5 weeks after planting, when storage roots were formed at the distal part. Our results point to down-regulation of fiber formation and lignification, together with up-regulation of starch biosynthesis, as the main events underlying storage-root formation, marking/highlighting several genes as potential regulators, providing a valuable database of genes for further research.

DELLA-NAC Interactions Mediate GA Signaling to Promote Secondary Cell Wall Formation in Cotton Stem

Gibberellins (GAs) promote secondary cell wall (SCW) development in plants, but the underlying molecular mechanism is still to be elucidated. Here, we employed a new system, the first internode of cotton, and the virus-induced gene silencing method to address this problem. We found that knocking down major DELLA genes via VIGS phenocopied GA treatment and significantly enhanced SCW formation in the xylem and phloem of cotton stems. Cotton DELLA proteins were found to interact with a wide range of SCW-related NAC proteins, and virus-induced gene silencing of these NAC genes inhibited SCW development with downregulated biosynthesis and deposition of lignin. The findings indicated a framework for the GA regulation of SCW formation; that is, the interactions between DELLA and NAC proteins mediated GA signaling to regulate SCW formation in cotton stems.

Identification and Characterization of Secondary Wall-Associated NAC Genes and Their Involvement in Hormonal Responses in Tobacco (Nicotiana tabacum)

Secondary wall-associated NAC (SWN) genes are a subgroup of NAC (NAM, ATAF, and CUC) transcription factors (TF) that play a key role in regulating secondary cell wall biosynthesis in plants. However, this gene family has not been systematically characterized, and their potential roles in response to hormones are unknown in Nicotiana tabacum. In this study, a total of 40 SWN genes, of which 12 from Nicotiana tomentosiformis, 13 from Nicotiana sylvestris, and 15 from Nicotiana tabacum, were successfully identified. The 15 SWNs from Nicotiana tabacum were further classified into three groups, namely, vascular-related NAC domain genes (NtVNDs), NAC secondary wall thickening promoting factor genes (NtNSTs), and secondary wall-associated NAC domain genes (NtSNDs). The protein characteristic, gene structure, and chromosomal location of 15 NtSWNs (also named Nt1 to Nt15) were also analyzed. The NtVND and NtNST group genes had five conserved subdomains in their N-terminal regions and a motif (LP[Q/x] L[E/x] S[P/A]) in their diverged C- terminal regions. Some hormones, dark and low-temperature related cis-acting elements, were significantly enriched in the promoters of NtSWN genes. A comprehensive expression profile analysis revealed that Nt4 and Nt12 might play a role in vein development. Others might be important for stem development. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) revealed that in the NtNST group, genes such as Nt7, Nt8, and Nt13 were more sensitive than the genes in NtVND and NtSND groups under abiotic stress conditions. A transactivation assay further suggested that Nt7, Nt8, and Nt13 showed a significant transactivation activity. Overall, SWN genes were finally identified and characterized in diploid and tetraploid tobacco, revealing new insights into their evolution, variation, and homology relationships. Transcriptome, cis-acting element, qRT-PCR, and transactivation assay analysis indicated the roles in hormonal and stress responses, which provided further resources in molecular mechanism and genetic improvement.

NST- and SND-subgroup NAC proteins coordinately act to regulate secondary cell wall formation in cotton

Secondary cell wall (SCW) has a strong impact on plant growth and adaptation to the environments. Previous studies have shown that NAC (NAM, ATAF1/2, and CUC2) transcription factors act as key regulators of SCW biosynthesis. However, the regulatory network triggered by NAC proteins is largely unknown, especially in cotton, a model plant for SCW development studies. Here, we show that several cotton NAC transcription factors are clustered in the same group with Arabidopsis secondary wall NACs (SWNs), including secondary wall-associated NAC domain protein1 (SND1) and NAC secondary wall thickening promoting factor1/2 (NST1/2), so we name these cotton orthologs as SND1s and NST1s. We found that simultaneous silencing of SND1s and NST1s led to severe xylem and phloem developmental defect in cotton stems, however silencing either SND1s or NST1s alone had no visible phenotype. Silencing both SND1s and NST1s but not one subgroup caused decreased expression of a set of SCW-associated genes, while over-expression of cotton SWNs in tobacco leaves resulted in SCW deposition. SWNs could bind the promoter of MYB46 and MYB83, which are highly expressed in SCW-rich tissues of cotton. In total, our data provide evidence that cotton SWNs positively and coordinately regulate SCW formation.

High-order mutants reveal an essential requirement for peroxidases but not laccases in Casparian strip lignification

Lignin is a defining polymer of vascular plants and of great physiological, ecological, and economical importance. Yet, its polymerization in the cell wall is still not understood. Lignin polymerizing enzymes, laccases and peroxidases, exist in vast numbers in plant genomes. By focusing on a specific lignin structure, the ring-like Casparian strips (CSs), we reduced candidate numbers and abolished essentially all laccases with detectable endodermal expression. Yet, not even slight defects in CS formation were detected. By contrast, a quintuple peroxidase mutant displayed a complete absence of CS. Our findings suggest that cells lignify differently depending on whether lignin is localized or ubiquitous and whether cells stay alive during and after lignification, as well as the composition of the cell wall.

Lignin has enabled plants to colonize land, grow tall, transport water within their bodies, and protect themselves against various stresses. Consequently, this polyphenolic polymer, impregnating cellulosic plant cell walls, is the second most abundant polymer on Earth. Yet, despite its great physiological, ecological, and economical importance, our knowledge of lignin biosynthesis in vivo, especially the polymerization steps within the cell wall, remains vague—specifically, the respective roles of the two polymerizing enzymes classes, laccases and peroxidases. One reason for this lies in the very high numbers of laccases and peroxidases encoded by 17 and 73 homologous genes, respectively, in Arabidopsis. Here, we have focused on a specific lignin structure, the ring-like Casparian strips (CSs) within the root endodermis. By reducing candidate numbers using cellular resolution expression and localization data and by boosting stacking of mutants using CRISPR-Cas9, we mutated the majority of laccases in Arabidopsis in a nonuple mutant—essentially abolishing laccases with detectable endodermal expression. Yet, we were unable to detect even slight defects in CS formation. By contrast, we were able to induce a complete absence of CS formation in a quintuple peroxidase mutant. Our findings are in stark contrast to the strong requirement of xylem vessels for laccase action and indicate that lignin in different cell types can be polymerized in very distinct ways. We speculate that cells lignify differently depending on whether lignin is localized or ubiquitous and whether cells stay alive during and after lignification, as well as the composition of the cell wall.

A review of current knowledge about the formation of native peridermal exocarp in fruit

The outer skin layer in any plant is essential in offering a protective barrier against water loss and pathogen attack. Within fleshy fruit, the skin supports internal cell layers and can provide the initial cues in attracting seed-dispersing animals. The skin of a fruit, termed the exocarp, is a key element of consumer preference and a target for many breeding programs. Across fruiting species there is a huge diversity of exocarp types and these range from a simple single living cell layer (epidermis) often covered with a waxy layer, to complex multicellular suberised and dead cell layers (periderm), with various intermediate russet forms in between. Each exocarp can be interspersed with other structures such as hairs or spines. The epidermis has been well characterised and remains pluripotent with the help of the cells immediately under the epidermis. The periderm, in contrast, is the result of secondary meristematic activity, which replaces the epidermal layers, and is not well characterised in fruits. In this review we explore the structure, composition and mechanisms that control the development of a periderm type fruit exocarp. We draw upon literature from non-fleshy fruit species that form periderm tissue, from which a considerable amount of research has been undertaken.

Ectopic expression of a poplar gene NAC13 confers enhanced tolerance to salinity stress in transgenic Nicotiana tabacum

NACs are one of the major transcription factor families in plants which play an important role in plant growth and development, as well as in adverse stress responses. In this study, we cloned a salt-inducible NAC transcription factor gene (NAC13) from a poplar variety 84K, followed by transforming it into both Nicotiana tabacum and Arabidopsis thaliana. Stable expression analysis of 35S::NAC13-GFP fusion protein in Arabidopsis indicated that NAC13 protein was localized to the nucleus. We also obtained five transgenic tobacco lines. Evidence from morphological and physiological characterization and salt treatment analyses indicated that in the transgenic tobacco the salt tolerance was enhanced, suggesting that NAC13 gene may function as a positive regulator in tobacco responses to salt stress.

Transcriptome analysis of salt-responsive and wood-associated NACs in Populus simonii × Populus nigra

BACKGROUND

NAC (NAM, ATAF1-2, and CUC2) family is one of the largest plant-specific transcription factor families known to play significant roles in plant development processes and stress responses.

RESULTS

In the study, a total of 112 NACs were identified to be differentially expressed in the comparisons of leaves and stems, leaves and roots, roots and stems of Populus simonii×P. nigra among 289 members by RNA-Seq. And 148, 144 and 134 NACs were detected to be salt-responsive in the roots, stems and leaves under 150 mM NaCl stress, respectively. Among them, a total of 53 salt-responsive NACs were shared across the three tissues. Under salt stress, 41/37 NACs were identified to be up/down-regulated in the leaves of Populus simonii × P.nigra among 170 non-redundant NACs by RT-qPCR, which was similar with RNA-Seq results. The expression pattern analysis of 6 NACs including four randomly up-regulated genes (NAC86, NAC105, NAC139 and NAC163) and two down-regulated genes (NAC15 and NAC149) indicated a few NACs showed specific temporal and spatial expression patterns in the three tissues of Populus simonii×P.nigra. Based on transcriptome screening and phylogenic analysis of differentially expressed NACs in different tissues under salt stress, 18 potential NACs associated with wood formation and 20 involved in stress responses were identified in Populus simonii×P.nigra.

CONCLUSIONS

The study further gains an understanding of the connection of tissue specificity and gene function in poplar, and lays the foundation of functional analysis of poplar NACs in stress responses.

Over-expression of transcription factor ARK1 gene leads to down-regulation of lignin synthesis related genes in hybrid poplar '717'

Improving wood growth rate and wood quality are worthy goals in forest genetics and breeding research. The ARK1 gene is one member of the ARBORKNOX family in all plants, which play an essential role in the process of plant growth and development, but the mechanism associated with its gene network regulation is poorly investigated. In order to generate over-expression transgenic hybrid poplar, the agrobacterium-mediated transformation was used to obtain transgenic hybrid poplar ‘717’ plants to provide insight into the function of the ARK1 gene in poplar. Moreover, the morphology of transgenic plants was observed, and transcriptome analysis was performed to explore the ARK1 gene function. The results showed that there were significant differences in pitch, stem diameter, petiole length, leaf width, leaf length and seedling height between ARK1 transgenic seedlings and non-transgenic seedlings. The transgenic seedlings usually had multiple branches and slender leaves, with some leaves not being fully developed. The results of transcriptome analysis showed that the differentially expressed genes were involved in the growth of poplars, including proteins, transcription factors and protein kinases. Genes related to the positive regulation in plant hormone signal transduction pathways were up-regulated, and the genes related to lignin synthesis were down-regulated. The RT-qPCR analysis confirmed the expression levels of the genes involved in the plant hormone signal transduction pathways and phenylpropanoid pathway. In conclusion, the ARK1 gene had a positive regulatory effect on plant growth, and the gene’s coding enzymes related to lignin synthesis were down-regulated.

Comparative transcriptome analyses of fruit development among pears, peaches, and strawberries provide new insights into single sigmoid patterns

BACKGROUND

Pear fruit exhibit a single sigmoid pattern during development, while peach and strawberry fruits exhibit a double sigmoid pattern. However, little is known about the differences between these two patterns.

RESULTS

In this study, fruit weights were measured and paraffin sections were made from fruitlet to maturated pear, peach, and strawberry samples. Results revealed that both single and double sigmoid patterns resulted from cell expansion, but not cell division. Comparative transcriptome analyses were conducted among pear, peach, and strawberry fruits at five fruit enlargement stages. Comparing the genes involved in these intervals among peaches and strawberries, 836 genes were found to be associated with all three fruit enlargement stages in pears (Model I). Of these genes, 25 were located within the quantitative trait locus (QTL) regions related to fruit weight and 90 were involved in cell development. Moreover, 649 genes were associated with the middle enlargement stage, but not early or late enlargement in pears (Model II). Additionally, 22 genes were located within the QTL regions related to fruit weight and 63 were involved in cell development. Lastly, dual-luciferase assays revealed that the screened bHLH transcription factors induced the expression of cell expansion-related genes, suggesting that the two models explain the single sigmoid pattern.

CONCLUSIONS

Single sigmoid patterns are coordinately mediated by Models I and II, thus, a potential gene regulation network for the single sigmoid pattern was proposed. These results enhance our understanding of the molecular regulation of fruit size in Rosaceae.

Transcriptional Regulation of Sorghum Stem Composition: Key Players Identified Through Co-expression Gene Network and Comparative Genomics Analyses

Most sorghum biomass accumulates in stem secondary cell walls (SCW). As sorghum stems are used as raw materials for various purposes such as feed, energy and fiber reinforced polymers, identifying the genes responsible for SCW establishment is highly important. Taking advantage of studies performed in model species, most of the structural genes contributing at the molecular level to the SCW biosynthesis in sorghum have been proposed while their regulatory factors have mostly not been determined. Validation of the role of several MYB and NAC transcription factors in SCW regulation in Arabidopsis and a few other species has been provided. In this study, we contributed to the recent efforts made in grasses to uncover the mechanisms underlying SCW establishment. We reported updated phylogenies of NAC and MYB in 9 different species and exploited findings from other species to highlight candidate regulators of SCW in sorghum. We acquired expression data during sorghum internode development and used co-expression analyses to determine groups of co-expressed genes that are likely to be involved in SCW establishment. We were able to identify two groups of co-expressed genes presenting multiple evidences of involvement in SCW building. Gene enrichment analysis of MYB and NAC genes provided evidence that while NAC SECONDARY WALL THICKENING PROMOTING FACTOR NST genes and SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN gene functions appear to be conserved in sorghum, NAC master regulators of SCW in sorghum may not be as tissue compartmentalized as in Arabidopsis. We showed that for every homolog of the key SCW MYB in Arabidopsis, a similar role is expected for sorghum. In addition, we unveiled sorghum MYB and NAC that have not been identified to date as being involved in cell wall regulation. Although specific validation of the MYB and NAC genes uncovered in this study is needed, we provide a network of sorghum genes involved in SCW both at the structural and regulatory levels.

Molecular characterisation and expression analysis of NAC transcription factor genes in wild Medicago falcata under abiotic stresses

The No apical meristem-Arabidopsis transcription activation factor-Cup-shaped cotyledon (NAC) proteins play vital roles in plant development processes and responses to abiotic stress. In this study, 146 unigenes were identified as NAC genes from wild Medicago falcata L. by RNA sequencing. Among these were 30 full-length NACs, which, except for MfNAC63, MfNAC64 and MfNAC91, contained a complete DNA-binding domain and a variable transcriptional activation region. Sequence analyses of MfNACs along with their Arabidopsis thaliana (L.) Heynh. counterparts allowed these proteins to be phylogenetically classified into nine groups. MfNAC35, MfNAC88, MfNAC79, MfNAC26 and MfNAC95 were found to be stress-responsive genes. The eight MfNAC genes that were chosen for further analysis had different expression abilities in the leaves, stems and roots of M. falcata. Additionally, their expression levels were regulated by salinity, drought and cold stress, and ABA. This study will be useful for understanding the roles of MfNACs in wild M. falcata and could provide important information for the selection of candidate genes associated with stress tolerance.

A Large Transposon Insertion in the stiff1 Promoter Increases Stalk Strength in Maize

Stalk lodging, which is generally determined by stalk strength, results in considerable yield loss and has become a primary threat to maize (Zea mays) yield under high-density planting. However, the molecular genetic basis of maize stalk strength remains unclear, and improvement methods remain inefficient. Here, we combined map-based cloning and association mapping and identified the gene stiff1 underlying a major quantitative trait locus for stalk strength in maize. A 27.2-kb transposable element insertion was present in the promoter of the stiff1 gene, which encodes an F-box domain protein. This transposable element insertion repressed the transcription of stiff1, leading to the increased cellulose and lignin contents in the cell wall and consequently greater stalk strength. Furthermore, a precisely edited allele of stiff1 generated through the CRISPR/Cas9 system resulted in plants with a stronger stalk than the unedited control. Nucleotide diversity analysis revealed that the promoter of stiff1 was under strong selection in the maize stiff-stalk group. Our cloning of stiff1 reveals a case in which a transposable element played an important role in maize improvement. The identification of stiff1 and our edited stiff1 allele pave the way for efficient improvement of maize stalk strength.

Identification of functional single nucleotide polymorphism of Populus trichocarpa PtrEPSP-TF and determination of its transcriptional effect

In plants, the phenylpropanoid pathway is responsible for the synthesis of a diverse array of secondary metabolites that include lignin monomers, flavonoids, and coumarins, many of which are essential for plant structure, biomass recalcitrance, stress defense, and nutritional quality. Our previous studies have demonstrated that Populus trichocarpa PtrEPSP-TF, an isoform of 5-enolpyruvylshikimate 3-phosphate (EPSP) synthase, has transcriptional activity and regulates phenylpropanoid biosynthesis in Populus. In this study, we report the identification of single nucleotide polymorphism (SNP) of PtrEPSP-TF that defines its functionality. Populus natural variants carrying this SNP were shown to have reduced lignin content. Here, we demonstrated that the SNP-induced substitution of 142nd amino acid (PtrEPSP-TFD142E) dramatically impairs the DNA-binding and transcriptional activity of PtrEPSP-TF. When introduced to a monocot species rice (Oryza sativa) in which an EPSP synthase isoform with the DNA-binding helix-turn-helix (HTH) motif is absent, the PtrEPSP-TF, but not PtrEPSP-TFD142E, activated genes in the phenylpropanoid pathway. More importantly, heterologous expression of PtrEPSP-TF uncovered five new transcriptional regulators of phenylpropanoid biosynthesis in rice. Collectively, this study identifies the key amino acid required for PtrEPSP-TF functionality and provides a strategy to uncover new transcriptional regulators in phenylpropanoid biosynthesis.

Analysis of NAC Domain Transcription Factor Genes of Tectona grandis L.f. Involved in Secondary Cell Wall Deposition

NAC proteins are one of the largest families of plant-specific transcription factors (TFs). They regulate diverse complex biological processes, including secondary xylem differentiation and wood formation. Recent genomic and transcriptomic studies of Tectona grandis L.f. (teak), one of the most valuable hardwood trees in the world, have allowed identification and analysis of developmental genes. In the present work, T. grandis NAC genes were identified and analyzed regarding to their evolution and expression profile during wood formation. We analyzed the recently published T. grandis genome, and identified 130 NAC proteins that are coded by 107 gene loci. These proteins were classified into 23 clades of the NAC family, together with Populus, Eucalyptus, and Arabidopsis. Data on transcript expression revealed specific temporal and spatial expression patterns for the majority of teak NAC genes. RT-PCR indicated expression of VND genes (Tg11g04450-VND2 and Tg15g08390-VND4) related to secondary cell wall formation in xylem vessels of 16-year-old juvenile trees. Our findings open a way to further understanding of NAC transcription factor genes in T. grandis wood biosynthesis, while they are potentially useful for future studies aiming to improve biomass and wood quality using biotechnological approaches.

Transcriptome landscape of Rafflesia cantleyi floral buds reveals insights into the roles of transcription factors and phytohormones in flower development

Rafflesia possesses unique biological features and known primarily for producing the world’s largest and existing as a single flower. However, to date, little is known about key regulators participating in Rafflesia flower development. In order to further understand the molecular mechanism that regulates Rafflesia cantleyi flower development, RNA-seq data from three developmental stages of floral bud, representing the floral organ primordia initiation, floral organ differentiation, and floral bud outgrowth, were analysed. A total of 89,890 transcripts were assembled of which up to 35% could be annotated based on homology search. Advanced transcriptome analysis using K-mean clustering on the differentially expressed genes (DEGs) was able to identify 12 expression clusters that reflect major trends and key transitional states, which correlate to specific developmental stages. Through this, comparative gene expression analysis of different floral bud stages identified various transcription factors related to flower development. The members of WRKY, NAC, bHLH, and MYB families are the most represented among the DEGs, suggesting their important function in flower development. Furthermore, pathway enrichment analysis also revealed DEGs that are involved in various phytohormone signal transduction events such as auxin and auxin transport, cytokinin and gibberellin biosynthesis. Results of this study imply that transcription factors and phytohormone signalling pathways play major role in Rafflesia floral bud development. This study provides an invaluable resource for molecular studies of the flower development process in Rafflesia and other plant species.

Molecular responses to freshwater limitation in the mangrove tree Avicennia germinans (Acanthaceae)

Environmental variation along the geographical space can shape populations by natural selection. In the context of global warming and changing precipitation regimes, it is crucial to understand the role of environmental heterogeneity in tropical trees adaptation, given their disproportional contribution to water and carbon biogeochemical cycles. Here, we investigated how heterogeneity in freshwater availability along tropical wetlands has influenced molecular variations of the black mangrove (Avicennia germinans). A total of 57 trees were sampled at seven sites differing markedly in precipitation regime and riverine freshwater inputs. Using 2,297 genome-wide single nucleotide polymorphic markers, we found signatures of natural selection by the association between variations in allele frequencies and environmental variables, including the precipitation of the warmest quarter and the annual precipitation. Additionally, we found candidate loci for selection based on statistical deviations from neutral expectations of interpopulation differentiation. Most candidate loci within transcribed sequences were functionally associated with central aspects of drought tolerance or plant response to drought. Moreover, our results suggest the occurrence of the rapid evolution of a population, probably in response to sudden and persistent limitations in plant access to soil water, following a road construction in 1974. Observations supporting rapid evolution included the reduction in tree size and changes in allele frequencies and in transcript expression associated with increased drought tolerance through the accumulation of osmoprotectants and antioxidants, biosynthesis of cuticles, protection against protein degradation, stomatal closure, photorespiration and photosynthesis. We describe a major role of spatial heterogeneity in freshwater availability in the specialization of this typically tropical tree.

Induction of PrMADS10 on the lower side of bent pine tree stems: potential role in modifying plant cell wall properties and wood anatomy

The molecular mechanisms underlying inclination responses in trees are unclear. In this study, we identified a MADS-box transcription factor differentially expressed early after inclination in the stems of Pinus radiata D. Don. PrMADS10 has a CDS of 582 bp and encodes a group II MADS-box transcription factor. We measured highest accumulation of this transcript on the lower side of inclined pine stems. In an effort to identify putative targets, we stably transformed Arabidopsis thaliana with a 35S::PrMADS10 construct. Transcriptome analysis revealed 1,219 genes differentially-expressed, with 690 and 529 genes up- and down-regulated respectively, when comparing the transgenic and wild-type. Differentially-expressed genes belong to different biological processes, but were enriched in cell wall remodeling and phenylpropanoid metabolic functions. Interestingly, lignin content was 30% higher in transgenic as compared to wild-type plants consistent with observed changes in gene expression. Differentially expressed transcription factors and phenylpropanoid genes were analyzed using STRING. Several MYB and NAC transcription factors showed interactions with genes of the phenylpropanoid pathway. Together, these results implicate PrMADS10 as a regulatory factor, triggering the expression of other transcription factors and genes involved in the synthesis of lignin.

Genome-Wide Investigation of the NAC Gene Family and Its Potential Association with the Secondary Cell Wall in Moso Bamboo

NAC (NAM, ATAF, and CUC) transcription factors (TFs) are implicated in the transcriptional regulation of diverse processes and have been characterized in a number of plant species. However, NAC TFs are still not well understood in bamboo, especially their potential association with the secondary cell wall (SCW). Here, 94 PeNACs were identified and characterized in moso bamboo (Phyllostachys edulis). Based on their gene structures and conserved motifs, the PeNACs were divided into 11 groups according to their homologs in Arabidopsis. PeNACs were expressed variously in different tissues of moso bamboo, suggesting their functional diversity. Fifteen PeNACs associated with the SCW were selected for co-expression analysis and validation. It was predicted that 396 genes were co-expressed with the 15 PeNACs, in which 16 and 55 genes were involved in the lignin catabolic process and cellulose biosynthetic process respectively. As the degree of lignification in the growing bamboo shoots increased, all 15 PeNACs were upregulated with a trend of rising first and then decreasing except PeNAC37, which increased continuously. These results indicated that these PeNACs might play important roles in SCW biosynthesis and lignification in bamboo shoots. Seven of 15 PeNACs had been found positively co-expressed with seven PeMYBs, and they had similar expression patterns with those of the PeMYBs in bamboo shoots. The targeted sites of miR164 were found in 16 PeNACs, of which three PeNACs associated with SCW were validated to have an opposite expression trend to that of miR164 in growing bamboo shoots. In addition, three PeNACs were selected and verified to have self-activation activities. These results provide comprehensive information of the NAC gene family in moso bamboo, which will be helpful for further functional studies of PeNACs to reveal the molecular regulatory mechanisms of bamboo wood property.

Integrated Genome-Scale Analysis Identifies Novel Genes and Networks Underlying Senescence in Maize

Premature senescence in annual crops reduces yield, while delayed senescence, termed stay-green, imposes positive and negative impacts on yield and nutrition quality. Despite its importance, scant information is available on the genetic architecture of senescence in maize (Zea mays) and other cereals. We combined a systematic characterization of natural diversity for senescence in maize and coexpression networks derived from transcriptome analysis of normally senescing and stay-green lines. Sixty-four candidate genes were identified by genome-wide association study (GWAS), and 14 of these genes are supported by additional evidence for involvement in senescence-related processes including proteolysis, sugar transport and signaling, and sink activity. Eight of the GWAS candidates, independently supported by a coexpression network underlying stay-green, include a trehalose-6-phosphate synthase, a NAC transcription factor, and two xylan biosynthetic enzymes. Source-sink communication and the activity of cell walls as a secondary sink emerge as key determinants of stay-green. Mutant analysis supports the role of a candidate encoding Cys protease in stay-green in Arabidopsis (Arabidopsis thaliana), and analysis of natural alleles suggests a similar role in maize. This study provides a foundation for enhanced understanding and manipulation of senescence for increasing carbon yield, nutritional quality, and stress tolerance of maize and other cereals.

Genome-Wide Analysis of NAC Gene Family in Betula pendula

NACs (NAM, ATAF1/2, and CUC2) are plant-specific transcription factors that play diverse roles in various plant developmental processes. In this study, we identified the NAC gene family in birch (Betula pendula) and further analyzed the function of BpNACs. Phylogenetic analysis reveals that the 114 BpNACs can be divided into seven subfamilies. We investigated the expression levels of these BpNACs in different tissues of birch including roots, xylem, leaves, and flowers, and the results showed that the BpNACs seem to be expressed higher in xylem and roots than leaves and flowers. In addition to tissue-specific expression analysis, we investigated the expression of BpNACs under low-temperature stress. A total of 21 BpNACs were differentially expressed under low-temperature stress, of which 17 were up-regulated, and four were down-regulated. Using the gene expression data, we reconstructed the gene co-expression network for the 21 low-temperature-responsive BpNACs. In conclusion, our results provide insight into the evolution of NAC genes in the B. pendula genome, and provide a basis for understanding the molecular mechanism for BpNAC-mediated cold responses in birch.

Transcriptional Reprogramming of Pea Leaves at Early Reproductive Stages

Pea (Pisum sativum L.) is an important source of dietary proteins. Nutrient recycling from leaves contributes to the accumulation of seed proteins and is a pivotal determinant of protein yields in this grain legume. The aim of this study was to unveil the transcriptional regulations occurring in pea leaves before the sharp decrease in chlorophyll breakdown. As a prelude to this study, a time-series analysis of ¹⁵N translocation at the whole plant level was performed, which indicated that nitrogen recycling among organs was highly dynamic during this period and varied depending on nitrate availability. Leaves collected on vegetative and reproductive nodes were further analyzed by transcriptomics. The data revealed extensive transcriptome changes in leaves of reproductive nodes during early seed development (from flowering to 14 days after flowering), including an up-regulation of genes encoding transporters, and particularly of sulfate that might sustain sulfur metabolism in leaves of the reproductive part. This developmental period was also characterized by a down-regulation of cell wall-associated genes in leaves of both reproductive and vegetative nodes, reflecting a shift in cell wall structure. Later on, 27 days after flowering, genes potentially switching the metabolism of leaves toward senescence were pinpointed, some of which are related to ribosomal RNA processing, autophagy, or transport systems. Transcription factors differentially regulated in leaves between stages were identified and a gene co-expression network pointed out some of them as potential regulators of the above-mentioned biological processes. The same approach was conducted in Medicago truncatula to identify shared regulations with this wild legume species. Altogether the results give a global view of transcriptional events in leaves of legumes at early reproductive stages and provide a valuable resource of candidate genes that could be targeted by reverse genetics to improve nutrient remobilization and/or delay catabolic processes leading to senescence.

A systems biology view of wood formation in Eucalyptus grandis trees submitted to different potassium and water regimes

Wood production in fast-growing Eucalyptus grandis trees is highly dependent on both potassium (K) fertilization and water availability but the molecular processes underlying wood formation in response to the combined effects of these two limiting factors remain unknown. E. grandis trees were submitted to four combinations of K-fertilization and water supply. Weighted gene co-expression network analysis and MixOmics-based co-regulation networks were used to integrate xylem transcriptome, metabolome and complex wood traits. Functional characterization of a candidate gene was performed in transgenic E. grandis hairy roots. This integrated network-based approach enabled us to identify meaningful biological processes and regulators impacted by K-fertilization and/or water limitation. It revealed that modules of co-regulated genes and metabolites strongly correlated to wood complex traits are in the heart of a complex trade-off between biomass production and stress responses. Nested in these modules, potential new cell-wall regulators were identified, as further confirmed by the functional characterization of EgMYB137. These findings provide new insights into the regulatory mechanisms of wood formation under stressful conditions, pointing out both known and new regulators co-opted by K-fertilization and/or water limitation that may potentially promote adaptive wood traits.

Modified expression of ZmMYB167 in Brachypodium distachyon and Zea mays leads to increased cell wall lignin and phenolic content

One of the challenges to enable targeted modification of lignocellulosic biomass from grasses for improved biofuel and biochemical production lies within our limited understanding of the transcriptional control of secondary cell wall biosynthesis. Here, we investigated the role of the maize MYB transcription factor ZmMYB167 in secondary cell wall biosynthesis and how modified ZmMYB167 expression in two distinct grass model species affects plant biomass and growth phenotypes. Heterologous expression of ZmMYB167 in the C3 model system Brachypodium led to mild dwarf phenotypes, increased lignin (~7% to 13%) and S-lignin monomer (~11% to 16%) content, elevated concentrations of cell wall-bound p-coumaric acid (~15% to 24%) and reduced biomass sugar release (~20%) compared to controls. Overexpression of ZmMYB167 in the C4 model system Zea mays increased lignin (~4% to 13%), p-coumaric acid (~8% to 52%) and ferulic acid (~13% to 38%) content but did not affect plant growth and development nor biomass recalcitrance. Taken together, modifying ZmMYB167 expression represents a target to alter lignin and phenolic content in grasses. The ZmMYB167 expression-induced discrepancies in plant phenotypic and biomass properties between the two grass model systems highlight the challenges and opportunities for MYB transcription factor-based genetic engineering approaches of grass biomass.

Genome-wide analysis of the NAC transcription factor family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

The NAC (NAM, ATAF1/2, and CUC2) transcription factor family represents a group of large plant-specific transcriptional regulators, participating in plant development and response to external stress. However, there is no comprehensive study on the NAC genes of Tartary buckwheat (Fagopyrum tataricum), a large group of extensively cultivated medicinal and edible plants. The recently published Tartary buckwheat genome permits us to explore all the FtNAC genes on a genome-wide basis.

RESULTS

In the present study, 80 NAC (FtNAC) genes of Tartary buckwheat were obtained and named uniformly according to their distribution on chromosomes. Phylogenetic analysis of NAC proteins in both Tartary buckwheat and Arabidopsis showed that the FtNAC proteins are widely distributed in 15 subgroups with one subgroup unclassified. Gene structure analysis found that multitudinous FtNAC genes contained three exons, indicating that the structural diversity in Tartary buckwheat NAC genes is relatively low. Some duplication genes of FtNAC have a conserved structure that was different from others, indicating that these genes may have a variety of functions. By observing gene expression, we found that FtNAC genes showed abundant differences in expression levels in various tissues and at different stages of fruit development.

CONCLUSIONS

In this research, 80 NAC genes were identified in Tartary buckwheat, and their phylogenetic relationships, gene structures, duplication, global expression and potential roles in Tartary buckwheat development were studied. Comprehensive analysis will be useful for a follow-up study of functional characteristics of FtNAC genes and for the development of high-quality Tartary buckwheat varieties.

Genome-wide identification and expression analysis of the MYB transcription factor in moso bamboo (Phyllostachys edulis)

The MYB family, one of the largest transcription factor (TF) families in the plant kingdom, plays vital roles in cell formation, morphogenesis and signal transduction, as well as responses to biotic and abiotic stresses. However, the underlying function of bamboo MYB TFs remains unclear. To gain insight into the status of these proteins, a total of 85 PeMYBs, which were further divided into 11 subgroups, were identified in moso bamboo (Phyllostachys edulis) by using a genome-wide search strategy. Gene structure analysis showed that PeMYBs were significantly different, with exon numbers varying from 4 to 13. Phylogenetic analysis indicated that PeMYBs clustered into 27 clades, of which the function of 18 clades has been predicted. In addition, almost all of the PeMYBs were differently expressed in leaves, panicles, rhizomes and shoots based on RNA-seq data. Furthermore, qRT-PCR analysis showed that 12 PeMYBs related to the biosynthesis and deposition of the secondary cell wall (SCW) were constitutively expressed, and their transcript abundance levels have changed significantly with increasing height of the bamboo shoots, for which the degree of lignification continuously increased. This result indicated that these PeMYBs might play fundamental roles in SCW thickening and bamboo shoot lignification. The present comprehensive and systematic study on the members of the MYB family provided a reference and solid foundation for further functional analysis of MYB TFs in moso bamboo.

Finding New Cell Wall Regulatory Genes in Populus trichocarpa Using Multiple Lines of Evidence

Understanding the regulatory network controlling cell wall biosynthesis is of great interest in Populus trichocarpa, both because of its status as a model woody perennial and its importance for lignocellulosic products. We searched for genes with putatively unknown roles in regulating cell wall biosynthesis using an extended network-based Lines of Evidence (LOE) pipeline to combine multiple omics data sets in P. trichocarpa, including gene coexpression, gene comethylation, population level pairwise SNP correlations, and two distinct SNP-metabolite Genome Wide Association Study (GWAS) layers. By incorporating validation, ranking, and filtering approaches we produced a list of nine high priority gene candidates for involvement in the regulation of cell wall biosynthesis. We subsequently performed a detailed investigation of candidate gene GROWTH-REGULATING FACTOR 9 (PtGRF9). To investigate the role of PtGRF9 in regulating cell wall biosynthesis, we assessed the genome-wide connections of PtGRF9 and a paralog across data layers with functional enrichment analyses, predictive transcription factor binding site analysis, and an independent comparison to eQTN data. Our findings indicate that PtGRF9 likely affects the cell wall by directly repressing genes involved in cell wall biosynthesis, such as PtCCoAOMT and PtMYB.41, and indirectly by regulating homeobox genes. Furthermore, evidence suggests that PtGRF9 paralogs may act as transcriptional co-regulators that direct the global energy usage of the plant. Using our extended pipeline, we show multiple lines of evidence implicating the involvement of these genes in cell wall regulatory functions and demonstrate the value of this method for prioritizing candidate genes for experimental validation.

Rice Genome-Scale Network Integration Reveals Transcriptional Regulators of Grass Cell Wall Synthesis

Grasses have evolved distinct cell wall composition and patterning relative to dicotyledonous plants. However, despite the importance of this plant family, transcriptional regulation of its cell wall biosynthesis is poorly understood. To identify grass cell wall-associated transcription factors, we constructed the Rice Combined mutual Ranked Network (RCRN). The RCRN covers >90% of annotated rice (Oryza sativa) genes, is high quality, and includes most grass-specific cell wall genes, such as mixed-linkage glucan synthases and hydroxycinnamoyl acyltransferases. Comparing the RCRN and an equivalent Arabidopsis network suggests that grass orthologs of most genetically verified eudicot cell wall regulators also control this process in grasses, but some transcription factors vary significantly in network connectivity between these divergent species. Reverse genetics, yeast-one-hybrid, and protoplast-based assays reveal that OsMYB61a activates a grass-specific acyltransferase promoter, which confirms network predictions and supports grass-specific cell wall synthesis genes being incorporated into conserved regulatory circuits. In addition, 10 of 15 tested transcription factors, including six novel Wall-Associated regulators (WAP1, WACH1, WAHL1, WADH1, OsMYB13a, and OsMYB13b), alter abundance of cell wall-related transcripts when transiently expressed. The results highlight the quality of the RCRN for examining rice biology, provide insight into the evolution of cell wall regulation, and identify network nodes and edges that are possible leads for improving cell wall composition.