The MYB46/MYB83-mediated transcriptional regulatory programme is a gatekeeper of secondary wall biosynthesis

Genomic insights into the origin, domestication and genetic basis of agronomic traits of castor bean

BACKGROUND

Castor bean (Ricinus communis L.) is an important oil crop, which belongs to the Euphorbiaceae family. The seed oil of castor bean is currently the only commercial source of ricinoleic acid that can be used for producing about 2000 industrial products. However, it remains largely unknown regarding the origin, domestication, and the genetic basis of key traits of castor bean.

RESULTS

Here we perform a de novo chromosome-level genome assembly of the wild progenitor of castor bean. By resequencing and analyzing 505 worldwide accessions, we reveal that the accessions from East Africa are the extant wild progenitors of castor bean, and the domestication occurs ~ 3200 years ago. We demonstrate that significant genetic differentiation between wild populations in Kenya and Ethiopia is associated with past climate fluctuation in the Turkana depression ~ 7000 years ago. This dramatic change in climate may have caused the genetic bottleneck in wild castor bean populations. By a genome-wide association study, combined with quantitative trait locus analysis, we identify important candidate genes associated with plant architecture and seed size.

CONCLUSIONS

This study provides novel insights of domestication and genome evolution of castor bean, which facilitates genomics-based breeding of this important oilseed crop and potentially other tree-like crops in future.

MYB Transcription Factors and Its Regulation in Secondary Cell Wall Formation and Lignin Biosynthesis during Xylem Development

The secondary wall is the main part of wood and is composed of cellulose, xylan, lignin, and small amounts of structural proteins and enzymes. Lignin molecules can interact directly or indirectly with cellulose, xylan and other polysaccharide molecules in the cell wall, increasing the mechanical strength and hydrophobicity of plant cells and tissues and facilitating the long-distance transportation of water in plants. MYBs (v-myb avian myeloblastosis viral oncogene homolog) belong to one of the largest superfamilies of transcription factors, the members of which regulate secondary cell-wall formation by promoting/inhibiting the biosynthesis of lignin, cellulose, and xylan. Among them, MYB46 and MYB83, which comprise the second layer of the main switch of secondary cell-wall biosynthesis, coordinate upstream and downstream secondary wall synthesis-related transcription factors. In addition, MYB transcription factors other than MYB46/83, as well as noncoding RNAs, hormones, and other factors, interact with one another to regulate the biosynthesis of the secondary wall. Here, we discuss the biosynthesis of secondary wall, classification and functions of MYB transcription factors and their regulation of lignin polymerization and secondary cell-wall formation during wood formation.

NAC and MYB Families and Lignin Biosynthesis-Related Members Identification and Expression Analysis in Melilotus albus

Melilotus albus is an annual or biennial legume species that adapts to extreme environments via its high stress tolerance. NAC and MYB transcription factors (TFs) are involved in the regulation of lignin biosynthesis, which has not been studied in M. albus. A total of 101 MaNAC and 299 MaMYB members were identified based on M. albus genome. Chromosome distribution and synteny analysis indicated that some genes underwent tandem duplication. Ka/Ks analysis suggested that MaNACs and MaMYBs underwent strong purifying selection. Stress-, hormone- and development-related cis-elements and MYB-binding sites were identified in the promoter regions of MaNACs and MaMYBs. Five MaNACs, two MaMYBs and ten lignin biosynthesis genes were identified as presenting coexpression relationships according to weighted gene coexpression network analysis (WGCNA). Eleven and thirteen candidate MaNAC and MaMYB genes related to lignin biosynthesis were identified, respectively, and a network comprising these genes was constructed which further confirmed the MaNAC and MaMYB relationship. These candidate genes had conserved gene structures and motifs and were highly expressed in the stems and roots, and qRT-PCR further verified the expression patterns. Overall, our results provide a reference for determining the precise role of NAC and MYB genes in M. albus and may facilitate efforts to breed low-lignin-content forage cultivars in the future.

Comparative gene expression analysis reveals that multiple mechanisms regulate the weeping trait in Prunus mume

Prunus mume (also known as Mei) is an important ornamental plant that is popular with Asians. The weeping trait in P. mume has attracted the attention of researchers for its high ornamental value. However, the formation of the weeping trait of woody plants is a complex process and the molecular basis of weeping stem development is unclear. Here, the morphological and histochemical characteristics and transcriptome profiles of upright and weeping stems from P. mume were studied. Significant alterations in the histochemical characteristics of upright and weeping stems were observed, and the absence of phloem fibres and less xylem in weeping stems might be responsible for their inability to resist gravity and to grow downward. Transcriptome analysis showed that differentially expressed genes (DEGs) were enriched in phenylpropanoid biosynthesis and phytohormone signal transduction pathways. To investigate the differential responses to hormones, upright and weeping stems were treated with IAA (auxin) and GA₃ (gibberellin A3), respectively, and the results revealed that weeping stems had a weaker IAA response ability and reduced upward bending angles than upright stems. On the contrary, weeping stems had increased upward bending angles than upright stems with GA₃ treatment. Compared to upright stems, interestingly, DEGs associated with diterpenoid biosynthesis and phenylpropanoid biosynthesis were significantly enriched after being treated with IAA, and expression levels of genes associated with phenylpropanoid biosynthesis, ABC transporters, glycosylphosphatidylinositol (GPI)—anchor biosynthesis were altered after being treated with GA₃ in weeping stems. Those results reveal that multiple molecular mechanisms regulate the formation of weeping trait in P. mume, which lays a theoretical foundation for the cultivation of new varieties.

Unravelling the Complex Interplay of Transcription Factors Orchestrating Seed Oil Content in Brassica napus L

Transcription factors (TFs) and their complex interplay are essential for directing specific genetic programs, such as responses to environmental stresses, tissue development, or cell differentiation by regulating gene expression. Knowledge regarding TF-TF cooperations could be promising in gaining insight into the developmental switches between the cultivars of Brassica napus L., namely Zhongshuang11 (ZS11), a double-low accession with high-oil- content, and Zhongyou821 (ZY821), a double-high accession with low-oil-content. In this regard, we analysed a time series RNA-seq data set of seed tissue from both of the cultivars by mainly focusing on the monotonically expressed genes (MEGs). The consideration of the MEGs enables the capturing of multi-stage progression processes that are orchestrated by the cooperative TFs and, thus, facilitates the understanding of the molecular mechanisms determining seed oil content. Our findings show that TF families, such as NAC, MYB, DOF, GATA, and HD-ZIP are highly involved in the seed developmental process. Particularly, their preferential partner choices as well as changes in their gene expression profiles seem to be strongly associated with the differentiation of the oil content between the two cultivars. These findings are essential in enhancing our understanding of the genetic programs in both cultivars and developing novel hypotheses for further experimental studies.

Major latex protein-like encoding genes contribute to Rhizoctonia solani defense responses in sugar beet

Sugar beets are attacked by several pathogens that cause root damages. Rhizoctonia (Greek for “root killer”) is one of them. Rhizoctonia root rot has become an increasing problem for sugar beet production and to decrease yield losses agronomical measures are adopted. Here, two partially resistant and two susceptible sugar beet genotypes were used for transcriptome analysis to discover new defense genes to this fungal disease, information to be implemented in molecular resistance breeding. Among 217 transcripts with increased expression at 2 days post-infection (dpi), three resistance-like genes were found. These genes were not significantly elevated at 5 dpi, a time point when increased expression of three Bet v I/Major latex protein (MLP) homologous genes BvMLP1, BvMLP2 and BvML3 was observed in the partially resistant genotypes. Quantitative RT-PCR analysis on diseased sugar beet seedlings validated the activity of BvMLP1 and BvMLP3 observed in the transcriptome during challenge by R. solani. The three BvMLP genes were cloned and overexpressed in Arabidopsis thaliana to further dissect their individual contribution. Transgenic plants were also compared to T-DNA mutants of orthologous MLP genes. Plants overexpressing BvMLP1 and BvMLP3 showed significantly less infection whereas additive effects were seen on Atmlp1/Atmlp3 double mutants. The data suggest that BvMLP1 and BvMLP3 may contribute to the reduction of the Rhizoctonia root rot disease in sugar beet. Impact on the defense reaction from other differential expressed genes observed in the study is discussed.

Information on the move: vascular tissue development in space and time during postembryonic root growth

Cascades of temporal and spatial regulation of gene expression play crucial roles in the vascular development in plant roots. Once vascular cell fates are determined, the timing of their differentiation is tightly controlled in a cell-autonomous manner. In contrast, extensive cell-to-cell communication contributes to the positioning and specifying of vascular cell types in the root meristem. Diverse factors moving short distances in a radial direction were found to be key contributors to these processes. Furthermore, signals from differentiated phloem were found to influence the phloem precursor and determine how the corresponding asymmetric cell division proceeded. These findings highlight the potential importance of underexplored types of intercellular communication in relation to vascular tissue development during postembryonic root growth.

The class II KNOX transcription factors KNAT3 and KNAT7 synergistically regulate monolignol biosynthesis in Arabidopsis

The function of the transcription factor KNOTTED ARABIDOPSIS THALIANA7 (KNAT7) is still unclear since it appears to be either a negative or a positive regulator for secondary cell wall deposition with its loss-of-function mutant displaying thicker interfascicular and xylary fiber cell walls but thinner vessel cell walls in inflorescence stems. To explore the exact function of KNAT7, class II KNOTTED1-LIKE HOMEOBOX (KNOX II) genes in Arabidopsis including KNAT3, KNAT4, and KNAT5 were studied together. By chimeric repressor technology, we found that both KNAT3 and KNAT7 repressors exhibited a similar dwarf phenotype. Both KNAT3 and KNAT7 genes were expressed in the inflorescence stems and the knat3 knat7 double mutant exhibited a dwarf phenotype similar to the repressor lines. A stem cross-section of knat3 knat7 displayed an enhanced irregular xylem phenotype as compared with the single mutants, and its cell wall thickness in xylem vessels and interfascicular fibers was significantly reduced. Analysis of cell wall chemical composition revealed that syringyl lignin was significantly decreased while guaiacyl lignin was increased in the knat3 knat7 double mutant. Coincidently, the knat3 knat7 transcriptome showed that most lignin pathway genes were activated, whereas the syringyl lignin-related gene Ferulate 5-Hydroxylase (F5H) was down-regulated. Protein interaction analysis revealed that KNAT3 and KNAT7 can form a heterodimer, and KNAT3, but not KNAT7, can interact with the key secondary cell wall formation transcription factors NST1/2, which suggests that the KNAT3-NST1/2 heterodimer complex regulates F5H to promote syringyl lignin synthesis. These results indicate that KNAT3 and KNAT7 synergistically work together to promote secondary cell wall biosynthesis.

Grapevine VlbZIP30 improves drought resistance by directly activating VvNAC17 and promoting lignin biosynthesis through the regulation of three peroxidase genes

Drought stress severely affects grapevine quality and yield, and recent reports have revealed that lignin plays an important role in protection from drought stress. Since little is known about lignin-mediated drought resistance in grapevine, we investigated its significance. Herein, we show that VlbZIP30 mediates drought resistance by activating the expression of lignin biosynthetic genes and increasing lignin deposition. Transgenic grapevine plants overexpressing VlbZIP30 exhibited lignin deposition (mainly G and S monomers) in the stem secondary xylem under control conditions, which resulted from the upregulated expression of VvPRX4 and VvPRX72. Overexpression of VlbZIP30 improves drought tolerance, characterized by a reduction in the water loss rate, maintenance of an effective photosynthesis rate, and increased lignin content (mainly G monomer) in leaves under drought conditions. Electrophoretic mobility shift assay, luciferase reporter assays, and chromatin immunoprecipitation-qPCR assays indicated that VlbZIP30 directly binds to the G-box cis-element in the promoters of lignin biosynthetic (VvPRX N1) and drought-responsive (VvNAC17) genes to regulate their expression. In summary, we report a novel VlbZIP30-mediated mechanism linking lignification and drought tolerance in grapevine. The results of this study may be of value for the development of molecular breeding strategies to produce drought-resistant fruit crops.

Characterization of the UDP-glycosyltransferase UGT72 Family in Poplar and Identification of Genes Involved in the Glycosylation of Monolignols

Monolignols are the building blocks for lignin polymerization in the apoplastic domain. Monolignol biosynthesis, transport, storage, glycosylation, and deglycosylation are the main biological processes partaking in their homeostasis. In Arabidopsis thaliana, members of the uridine diphosphate-dependent glucosyltransferases UGT72E and UGT72B subfamilies have been demonstrated to glycosylate monolignols. Here, the poplar UGT72 family, which is clustered into four groups, was characterized: Group 1 UGT72AZ1 and UGT72AZ2, homologs of Arabidopsis UGT72E1-3, as well as group 4 UGT72B37 and UGT72B39, homologs of Arabidopsis UGT72B1-3, glycosylate monolignols. In addition, promoter-GUS analyses indicated that poplar UGT72 members are expressed within vascular tissues. At the subcellular level, poplar UGT72s belonging to group 1 and group 4 were found to be associated with the nucleus and the endoplasmic reticulum. However, UGT72A2, belonging to group 2, was localized in bodies associated with chloroplasts, as well as possibly in chloroplasts. These results show a partial conservation of substrate recognition between Arabidopsis and poplar homologs, as well as divergent functions between different groups of the UGT72 family, for which the substrates remain unknown.

Identification and Functional Characterization of PtoMYB055 Involved in the Regulation of the Lignin Biosynthesis Pathway in Populus tomentosa

Wood, which is mainly composed of lignified secondary cell wall, is the most abundant biomass in woody plants. Previous studies have revealed that R2R3-type MYB transcription factors are important regulators of the formation of the secondary cell wall in vascular plants. In this study, we isolated the R2R3-type MYB transcription factor gene PtoMYB055, which is mainly expressed in xylem and phloem tissue, from Populus tomentosa and demonstrate that PtoMYB055 is a key regulator of lignin biosynthesis. PtoMYB055 as a transcriptional activator is localized to the nucleus. Overexpression of PtoMYB055 upregulates expression of lignin biosynthetic genes in transgenic poplar plants, resulting in ectopic deposition of lignin in phloem tissue and an increase in thickness of the secondary cell wall. In sum, PtoMYB055 is a transcriptional activator that is involved in regulating lignin biosynthesis during the formation of the secondary cell wall in poplar.

Double knockout of OsWRKY36 and OsWRKY102 boosts lignification with altering culm morphology of rice

Breeding to enrich lignin, a major component of lignocelluloses, in plants contributes to enhanced applications of lignocellulosic biomass into solid biofuels and valuable aromatic chemicals. To collect information on enhancing lignin deposition in grass species, important lignocellulose feedstocks, we generated rice (Oryza sativa) transgenic lines deficient in OsWRKY36 and OsWRKY102, which encode putative transcriptional repressors for secondary cell wall formation. We used CRISPR/Cas9-mediated targeted mutagenesis and closely characterized their altered cell walls using chemical and nuclear magnetic resonance (NMR) methods. Both OsWRKY36 and OsWRKY102 mutations significantly increased lignin content by up to 28 % and 32 %, respectively. Additionally, OsWRKY36/OsWRKY102-double-mutant lines displayed lignin enrichment of cell walls (by up to 41 %) with substantially altered culm morphology over the single-mutant lines as well as the wild-type controls. Our chemical and NMR analyses showed that relative abundances of guaiacyl and p-coumarate units were slightly higher and lower, respectively, in the WRKY mutant lignins compared with those in the wild-type lignins. Our results provide evidence that both OsWRKY36 and OsWRKY102 are associated with repression of rice lignification.

Involvement of VNS NAC-domain transcription factors in tracheid formation in Pinus taeda

Vascular plants have two types of water-conducting cells, xylem vessel cells (in angiosperms) and tracheid cells (in ferns and gymnosperms). These cells are commonly characterized by secondary cell wall (SCW) formation and programmed cell death (PCD), which increase the efficiency of water conduction. The differentiation of xylem vessel cells is regulated by a set of NAC (NAM, ATAF1/2 and CUC2) transcription factors, called the VASCULAR-RELATED NAC-DOMAIN (VND) family, in Arabidopsis thaliana Linne. The VNDs regulate the transcriptional induction of genes required for SCW formation and PCD. However, information on the transcriptional regulation of tracheid cell differentiation is still limited. Here, we performed functional analysis of loblolly pine (Pinus taeda Linne) VND homologs (PtaVNS, for VND, NST/SND, SMB-related protein). We identified five PtaVNS genes in the loblolly pine genome, and four of these PtaVNS genes were highly expressed in tissues with tracheid cells, such as shoot apices and developing xylem. Transient overexpression of PtaVNS genes induced xylem vessel cell-like patterning of SCW deposition in tobacco (Nicotiana benthamiana Domin) leaves, and up-regulated the promoter activities of loblolly pine genes homologous to SCW-related MYB transcription factor genes and cellulose synthase genes, as well as to cysteine protease genes for PCD. Collectively, our data indicated that PtaVNS proteins possess transcriptional activity to induce the molecular programs required for tracheid formation, i.e., SCW formation and PCD. Moreover, these findings suggest that the VNS-MYB-based transcriptional network regulating water-conducting cell differentiation in angiosperm and moss plants is conserved in gymnosperms.

Construction of a hierarchical gene regulatory network centered around a transcription factor

We have modified a multitude of transcription factors (TFs) in numerous plant species and some animal species, and obtained transgenic lines that exhibit phenotypic alterations. Whenever we observe phenotypic changes in a TF's transgenic lines, we are always eager to identify its target genes, collaborative regulators and even upstream high hierarchical regulators. This issue can be addressed by establishing a multilayered hierarchical gene regulatory network (ML-hGRN) centered around a given TF. In this article, a practical approach for constructing an ML-hGRN centered on a TF using a combined approach of top-down and bottom-up network construction methods is described. Strategies for constructing ML-hGRNs are vitally important, as these networks provide key information to advance our understanding of how biological processes are regulated.

MYB20, MYB42, MYB43, and MYB85 Regulate Phenylalanine and Lignin Biosynthesis during Secondary Cell Wall Formation

Lignin is a phenylpropanoid-derived polymer that functions as a major component of cell walls in plant vascular tissues. Biosynthesis of the aromatic amino acid Phe provides precursors for many secondary metabolites, including lignins and flavonoids. Here, we discovered that MYB transcription factors MYB20, MYB42, MYB43, and MYB85 are transcriptional regulators that directly activate lignin biosynthesis genes and Phe biosynthesis genes during secondary wall formation in Arabidopsis (Arabidopsis thaliana). Disruption of MYB20, MYB42, MYB43, and MYB85 resulted in growth development defects and substantial reductions in lignin biosynthesis. In addition, our data showed that these MYB proteins directly activated transcriptional repressors that specifically inhibit flavonoid biosynthesis, which competes with lignin biosynthesis for Phe precursors. Together, our results provide important insights into the molecular framework for the lignin biosynthesis pathway.

SiMYB56 Confers Drought Stress Tolerance in Transgenic Rice by Regulating Lignin Biosynthesis and ABA Signaling Pathway

Foxtail millet (Setaria italica) originated in China and is generally cultivated in arid and barren soil. Through long-term harsh environmental selection, foxtail millet has acquired significant drought resistance. However, the molecular mechanism of foxtail millet drought resistance is still unknown. Here, we identified a drought-induced R2R3-MYB transcription factor SiMYB56 in foxtail millet. Overexpression of SiMYB56 significantly enhances tolerance to drought stress in transgenic rice plants at both the vegetative and the reproductive stage and has no adverse effect on its normal growth. Compared with wild-type controls, SiMYB56-overexpressing rice plants had lower MDA content and higher lignin content under drought conditions. Quantitative real-time PCR and Transcriptional activity assays demonstrated that SiMYB56 could activate expression of lignin biosynthesis genes under drought conditions. Also, we found that overexpression of SiMYB56 can led to ABA accumulation in the seeds transgenic rice plants. Further experiments showed that Overexpression of SiMYB56 can upregulate the expression of ABA synthesis and response related genes under drought conditions. In conclusion, SiMYB56 may enhance the drought resistance of transgenic rice plants by regulating lignin biosynthesis and ABA signaling pathway, making SiMYB56 a candidate gene for drought resistance improvement in gramineous crops.

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.

Transcriptional and Post-transcriptional Regulation of Lignin Biosynthesis Pathway Genes in Populus

Lignin is a heterogeneous polymer of aromatic subunits derived from phenylalanine. It is polymerized in intimate proximity to the polysaccharide components in plant cell walls and provides additional rigidity and compressive strength for plants. Understanding the regulatory mechanisms of lignin biosynthesis is important for genetic modification of the plant cell wall for agricultural and industrial applications. Over the past 10 years the transcriptional regulatory model of lignin biosynthesis has been established in plants. However, the role of post-transcriptional regulation is still largely unknown. Increasing evidence suggests that lignin biosynthesis pathway genes are also regulated by alternative splicing, microRNA, and long non-coding RNA. In this review, we briefly summarize recent progress on the transcriptional regulation, then we focus on reviewing progress on the post-transcriptional regulation of lignin biosynthesis pathway genes in the woody model plant Populus.

The Progression of Xylem Vessel Cell Differentiation is Dependent on the Activity Level of VND7 in Arabidopsis thaliana

Xylem vessels are important for water conduction in vascular plants. The VASCULAR-RELATED NAC-DOMAIN (VND) family proteins, master regulators of xylem vessel cell differentiation in Arabidopsis thaliana, can upregulate a set of genes required for xylem vessel cell differentiation, including those involved in secondary cell wall (SCW) formation and programmed cell death (PCD); however, it is not fully understood how VND activity levels influence these processes. Here, we examined the Arabidopsis VND7-VP16-GR line, in which VND7 activity is post-translationally activated by treatments with different concentrations of dexamethasone (DEX), a synthetic glucocorticoid. Our observations showed that 1 nM DEX induced weak SCW deposition, but not PCD, whereas 10 or 100 nM DEX induced both SCW deposition and PCD. The decreased chlorophyll contents and SCW deposition were apparent after 24 h of 100 nM DEX treatment, but became evident only after 48 h of 10 nM DEX treatment. Moreover, the lower DEX concentrations delayed the upregulation of VND7 downstream genes, and decreased their induction levels. They collectively suggest that the regulation of VND activity is important not only to initiate xylem vessel cell differentiation, but also regulate the quality of the xylem vessels through VND-activity-dependent upregulation of the PCD- and SCW-related genes.

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 Profiling Unravels a Vital Role of Pectin and Pectinase in Anther Dehiscence in Chrysanthemum

Chrysanthemum (Chrysanthemum morifolium (Ramat.) Kitamura) plants have great ornamental value, but their flowers can also be a source of pollen contamination. Previously, morphological and cytological studies have shown that anthers of some chrysanthemum cultivars such as 'Qx-115' fail to dehisce, although the underlying mechanism is largely unknown. In this study, we investigated the molecular basis of anther indehiscence in chrysanthemum via transcriptome analysis of a dehiscent cultivar ('Qx-097') and an indehiscent cultivar ('Qx-115'). We also measured related physiological indicators during and preceding the period of anther dehiscence. Our results showed a difference in pectinase accumulation and activity between the two cultivars during dehiscence. Detection of de-esterified pectin and highly esterified pectin in anthers during the period preceding anther dehiscence using LM19 and LM20 monoclonal antibodies showed that both forms of pectin were absent in the stomium region of 'Qx-097' anthers but were abundant in that of 'Qx-115' anthers. Analysis of transcriptome data revealed a significant difference in the expression levels of two transcription factor-encoding genes, CmLOB27 and CmERF72, between 'Qx-097' and 'Qx-115' during anther development. Transient overexpression of CmLOB27 and CmERF72 separately in tobacco leaves promoted pectinase biosynthesis. We conclude that CmLOB27 and CmERF72 are involved in the synthesis of pectinase, which promotes the degradation of pectin. Our results lay a foundation for further investigation of the role of CmLOB27 and CmERF72 transcription factors in the process of anther dehiscence in chrysanthemum.

Transcriptomic analysis of early fruit development in Chinese white pear (Pyrus bretschneideri Rehd.) and functional identification of PbCCR1 in lignin biosynthesis

BACKGROUND

The content of stone cells and lignin is one of the key factors affecting the quality of pear fruit. In a previous study, we determined the developmental regularity of stone cells and lignin in 'Dangshan Su' pear fruit 15-145 days after pollination (DAP). However, the development of fruit stone cells and lignin before 15 DAP has not been heavily researched.

RESULTS

In this study, we found that primordial stone cells began to appear at 7 DAP and that the fruit had formed a large number of stone cells at 15 DAP. Subsequently, transcriptome sequencing was performed on fruits at 0, 7, and 15 DAP and identified 3834 (0 vs. 7 DAP), 4049 (7 vs. 15 DAP) and 5763 (0 vs. 15 DAP) DEGs. During the 7-15 DAP period, a large number of key enzyme genes essential for lignin biosynthesis are gradually up-regulated, and their expression pattern is consistent with the accumulation of lignin in this period. Further analysis found that the biosynthesis of S-type lignin in 'Dangshan Su' pear does not depend on the catalytic activity of PbSAD but is primarily generated by the catalytic activity of caffeoyl-CoA through CCoAOMT, CCR, F5H, and CAD. We cloned PbCCR1, 2 and analysed their functions in Chinese white pear lignin biosynthesis. PbCCR1 and 2 have a degree of functional redundancy; both demonstrate the ability to participate in lignin biosynthesis. However, PbCCR1 may be the major gene for lignin biosynthesis, while PbCCR2 has little effect on lignin biosynthesis.

CONCLUSIONS

Our results revealed that 'Dangshan Su' pear began to form a large number of stone cells and produce lignin after 7 DAP and mainly accumulated materials from 0 to 7 DAP. PbCCR1 is mainly involved in the biosynthesis of lignin in 'Dangshan Su' pear and plays a positive role in lignin biosynthesis.

Mining MYB transcription factors from the genomes of orchids (Phalaenopsis and Dendrobium) and characterization of an orchid R2R3-MYB gene involved in water-soluble polysaccharide biosynthesis

Members of the MYB superfamily act as regulators in a wide range of biological processes in plants. Despite this, the MYB superfamily from the Orchidaceae has not been identified, and MYB genes related to bioactive water-soluble polysaccharide (WSP) biosynthesis are relatively unknown. In this study, we identified 159 and 165 MYB genes from two orchids, Phalaenopsis equestris and Dendrobium officinale, respectively. The MYB proteins were classified into four MYB classes in both orchids: MYB-related (MYBR), R2R3-MYB, 3R-MYB and atypical MYB proteins. The MYBR proteins in both orchids were classified into five subfamilies and 12 genes were strongly up-regulated in response to cold stress in D. officinale. The R2R3-MYB proteins were both divided into 31 clades in P. equestris and D. officinale. Among these clades, nine contained MYB TFs related to secondary cell wall biosynthesis or testa mucilage biosynthesis in Arabidopsis thaliana. In D. officinale, 10 candidate genes showed an expression pattern corresponding to changes in the WSP content. Overexpression of one of these candidate genes (DoMYB75) in A. thaliana increased seed WSP content by about 14%. This study provides information about MYB genes in two orchids that will further help to understand the transcriptional regulation of WSP biosynthesis in these orchids as well as other plant species.

Mining MYB transcription factors from the genomes of orchids (Phalaenopsis and Dendrobium) and characterization of an orchid R2R3-MYB gene involved in water-soluble polysaccharide biosynthesis

Members of the MYB superfamily act as regulators in a wide range of biological processes in plants. Despite this, the MYB superfamily from the Orchidaceae has not been identified, and MYB genes related to bioactive water-soluble polysaccharide (WSP) biosynthesis are relatively unknown. In this study, we identified 159 and 165 MYB genes from two orchids, Phalaenopsis equestris and Dendrobium officinale, respectively. The MYB proteins were classified into four MYB classes in both orchids: MYB-related (MYBR), R2R3-MYB, 3R-MYB and atypical MYB proteins. The MYBR proteins in both orchids were classified into five subfamilies and 12 genes were strongly up-regulated in response to cold stress in D. officinale. The R2R3-MYB proteins were both divided into 31 clades in P. equestris and D. officinale. Among these clades, nine contained MYB TFs related to secondary cell wall biosynthesis or testa mucilage biosynthesis in Arabidopsis thaliana. In D. officinale, 10 candidate genes showed an expression pattern corresponding to changes in the WSP content. Overexpression of one of these candidate genes (DoMYB75) in A. thaliana increased seed WSP content by about 14%. This study provides information about MYB genes in two orchids that will further help to understand the transcriptional regulation of WSP biosynthesis in these orchids as well as other plant species.

Identification and functional evaluation of accessible chromatin associated with wood formation in Eucalyptus grandis

Accessible chromatin changes dynamically during development and harbours functional regulatory regions which are poorly understood in the context of wood development. We explored the importance of accessible chromatin in Eucalyptus grandis in immature xylem generally, and MYB transcription factor-mediated transcriptional programmes specifically. We identified biologically reproducible DNase I Hypersensitive Sites (DHSs) and assessed their functional significance in immature xylem through their associations with gene expression, epigenomic data and DNA sequence conservation. We identified in vitro DNA binding sites for six secondary cell wall-associated Eucalyptus MYB (EgrMYB) transcription factors using DAP-seq, reconstructed protein-DNA networks of predicted targets based on binding sites within or outside DHSs and assessed biological enrichment of these networks with published datasets. 25 319 identified immature xylem DHSs were associated with increased transcription and significantly enriched for various epigenetic signatures (H3K4me3, H3K27me3, RNA pol II), conserved noncoding sequences and depleted single nucleotide variants. Predicted networks built from EgrMYB binding sites located in accessible chromatin were significantly enriched for systems biology datasets relevant to wood formation, whereas those occurring in inaccessible chromatin were not. Our study demonstrates that DHSs in E. grandis immature xylem, most of which are intergenic, are of functional significance to gene regulation in this tissue.

Function, dynamics and evolution of network motif modules in integrated gene regulatory networks of worm and plant

Gene regulatory networks (GRNs) consist of different molecular interactions that closely work together to establish proper gene expression in time and space. Especially in higher eukaryotes, many questions remain on how these interactions collectively coordinate gene regulation. We study high quality GRNs consisting of undirected protein–protein, genetic and homologous interactions, and directed protein–DNA, regulatory and miRNA–mRNA interactions in the worm Caenorhabditis elegans and the plant Arabidopsis thaliana. Our data-integration framework integrates interactions in composite network motifs, clusters these in biologically relevant, higher-order topological network motif modules, overlays these with gene expression profiles and discovers novel connections between modules and regulators. Similar modules exist in the integrated GRNs of worm and plant. We show how experimental or computational methodologies underlying a certain data type impact network topology. Through phylogenetic decomposition, we found that proteins of worm and plant tend to functionally interact with proteins of a similar age, while at the regulatory level TFs favor same age, but also older target genes. Despite some influence of the duplication mode difference, we also observe at the motif and module level for both species a preference for age homogeneity for undirected and age heterogeneity for directed interactions. This leads to a model where novel genes are added together to the GRNs in a specific biological functional context, regulated by one or more TFs that also target older genes in the GRNs. Overall, we detected topological, functional and evolutionary properties of GRNs that are potentially universal in all species.

Wood forming tissue-specific bicistronic expression of PdGA20ox1 and PtrMYB221 improves both the quality and quantity of woody biomass production in a hybrid poplar

With the exponential growth of the human population and industrial developments, research on renewable energy resources is required to alleviate environmental and economic impacts caused by the consumption of fossil fuels. In this study, we present a synthetic biological application of a wood forming tissue-specific bicistronic gene expression system to improve both the quantity and quality of woody biomass to minimize undesirable growth penalties. Our transgenic poplars, designed to express both PdGA20ox1 (a GA20-oxidase from Pinus densiflora producing bioactive gibberellin, GA) and PtrMYB221 (a MYB transcription factor negatively regulating lignin biosynthesis) under the developing xylem (DX) tissue-specific promoter (i.e., DX15::PdGA20ox1-2A-PtrMYB221 poplar), resulted in a 2-fold increase in biomass quantity compared to wild-type (WT), without undesirable growth defects. A similar phenotype was observed in transgenic Arabidopsis plants harboring the same gene constructs. These phenotypic consequences were further verified in the field experiments. Importantly, our transgenic poplars exhibited an improved quality of biomass with reduced lignin content (~16.0 wt%) but increased holocellulose content (~6.6 wt%). Furthermore, the saccharification efficiency of our transgenic poplar increased significantly by up to 8%. Our results demonstrate that the controlled production of both GA and a secondary wall modifying regulator in the same spatio-temporal manner can be utilized as an efficient biotechnological tool for producing the desired multi-purpose woody biomass.

OsMYB108 loss-of-function enriches p-coumaroylated and tricin lignin units in rice cell walls

Breeding approaches to enrich lignins in biomass could be beneficial to improving the biorefinery process because lignins increase biomass heating value and represent a potent source of valuable aromatic chemicals. However, despite the fact that grasses are promising lignocellulose feedstocks, limited information is yet available for molecular-breeding approaches to upregulate lignin biosynthesis in grass species. In this study, we generated lignin-enriched transgenic rice (Oryza sativa), a model grass species, via targeted mutagenesis of the transcriptional repressor OsMYB108 using CRISPR/Cas9-mediated genome editing. The OsMYB108-knockout rice mutants displayed increased expressions of lignin biosynthetic genes and enhanced lignin deposition in culm cell walls. Chemical and two-dimensional nuclear magnetic resonance (NMR) analyses revealed that the mutant cell walls were preferentially enriched in γ-p-coumaroylated and tricin lignin units, both of which are typical and unique components in grass lignins. NMR analysis also showed that the relative abundances of major lignin linkage types were altered in the OsMYB108 mutants.

Genome-wide identification and expression analysis of R2R3, 3R- and 4R-MYB transcription factors during lignin biosynthesis in flax (Linum usitatissimum)

MYB transcription factors (TFs) have vital roles in regulating lignin or fiber development. Flax (Linum usitatissimum) is known as one of the plants with high fiber production capacity. However, no studies have been conducted to identify and characterize MYB TFs in the flax genome. Results showed that flax genome harbours 167 R2R3, seven 3R, and one 4R-type MYB TFs. 22 MYB genes (%13) were estimated to be tandem duplicated dated around 13.3-86.98 Mya. 130 flax MYB members have apparent orthologous with Arabidopsis, in which 17 R2R3 MYBs are associated with lignin biosynthesis. MYB062, MYB072, MYB096, MYB141, and MYB146 genes were up-regulated in tissues having higher lignin production capacity. In opposite, MYB012 and MYB113 genes were down-regulated which points out the involvement of those genes in the lignin biosynthesis mechanism. This comprehensive study can provide a basis for understanding the role of MYBs in fiber or lignin production in flax.

Regulation of Gene Expression of Methionine Sulfoxide Reductases and Their New Putative Roles in Plants

Oxidation of methionine to methionine sulfoxide is a type of posttranslational modification reversed by methionine sulfoxide reductases (Msrs), which present an exceptionally high number of gene copies in plants. The side-form general antioxidant function-specific role of each Msr isoform has not been fully studied. Thirty homologous genes of Msr type A (MsrA) and type B (MsrB) that originate from the genomes of Arabidopsis thaliana, Populus trichocarpa, and Oryza sativa were analyzed in silico. From 109 to 201 transcription factors and responsive elements were predicted for each gene. Among the species, 220 and 190 common transcription factors and responsive elements were detected for the MsrA and MsrB isoforms, respectively. In a comparison of 14 MsrA and 16 MsrB genes, 424 transcription factors and responsive elements were reported in both types of genes, with almost ten times fewer unique elements. The transcription factors mainly comprised plant growth and development regulators, transcription factors important in stress responses with significant overrepresentation of the myeloblastosis viral oncogene homolog (MYB) and no apical meristem, Arabidopsis transcription activation factor and cup-shaped cotyledon (NAC) families and responsive elements sensitive to ethylene, jasmonate, sugar, and prolamine. Gene Ontology term-based functional classification revealed that cellular, metabolic, and developmental process terms and the response to stimulus term dominated in the biological process category. Available experimental transcriptomic and proteomic data, in combination with a set of predictions, gave coherent results validating this research. Thus, new manners Msr gene expression regulation, as well as new putative roles of Msrs, are proposed.

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

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

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.

Arabidopsis Group IIId ERF proteins positively regulate primary cell wall-type CESA genes

The cell wall determines morphology and the environmental responses of plant cells. The primary cell wall (PCW) is produced during cell division and expansion, determining the cell shape and volume. After cell expansion, specific types of plant cells produce a lignified wall, known as a secondary cell wall (SCW). We functionally analyzed Group IIId Arabidopsis AP2/EREBP genes, namely ERF34, ERF35, ERF38, and ERF39, which are homologs of a rice ERF gene previously proposed to be related to SCW biosynthesis. Expression analysis revealed that these four genes are expressed in regions related to cell division and/or cell differentiation in seedlings (i.e., shoot apical meristems, the primordia of leaves and lateral roots, trichomes, and central cylinder of primary roots) and flowers (i.e., vascular tissues of floral organs and replums and/or valve margins of pistils). Overexpression of ERF genes significantly upregulated PCW-type, but not SCW-type, CESA genes encoding cellulose synthase catalytic subunits in Arabidopsis seedlings. Transient co-expression reporter analysis indicated that ERF35, ERF38, and ERF39 possess transcriptional activator activity, and that ERF34, ERF35, ERF38, and ERF39 upregulated the promoter activity of CESA1, a PCW-type CESA gene, through the DRECRTCOREAT elements, the core cis-acting elements known to be recognized by AP2/ERF proteins. Together, our findings show that Group IIId ERF genes are positive transcriptional regulators of PCW-type CESA genes in Arabidopsis and are possibly involved in modulating cellulose biosynthesis in response to developmental requirements and environmental stimuli.

SECONDARY WALL ASSOCIATED MYB1 is a positive regulator of secondary cell wall thickening in Brachypodium distachyon and is not found in the Brassicaceae

Grass biomass is comprised chiefly of secondary walls that surround fiber and xylem cells. A regulatory network of interacting transcription factors in part regulates cell wall thickening. We identified Brachypodium distachyon SECONDARY WALL ASSOCIATED MYB1 (SWAM1) as a potential regulator of secondary cell wall biosynthesis based on gene expression, phylogeny, and transgenic plant phenotypes. SWAM1 interacts with cellulose and lignin gene promoters with preferential binding to AC-rich sequence motifs commonly found in the promoters of cell wall-related genes. SWAM1 overexpression (SWAM-OE) lines had greater above-ground biomass with only a slight change in flowering time while SWAM1 dominant repressor (SWAM1-DR) plants were severely dwarfed with a striking reduction in lignin of sclerenchyma fibers and stem epidermal cell length. Cellulose, hemicellulose, and lignin genes were significantly down-regulated in SWAM1-DR plants and up-regulated in SWAM1-OE plants. There was no reduction in bioconversion yield in SWAM1-OE lines; however, it was significantly increased for SWAM1-DR samples. Phylogenetic and syntenic analyses strongly suggest that the SWAM1 clade was present in the last common ancestor between eudicots and grasses, but is not in the Brassicaceae. Collectively, these data suggest that SWAM1 is a transcriptional activator of secondary cell wall thickening and biomass accumulation in B. distachyon.

MdMYB88 and MdMYB124 Enhance Drought Tolerance by Modulating Root Vessels and Cell Walls in Apple

Water deficit is one of the main limiting factors in apple (Malus × domestica Borkh.) cultivation. Root architecture plays an important role in the drought tolerance of plants; however, research efforts to improve drought tolerance of apple trees have focused on aboveground targets. Due to the difficulties associated with visualization and data analysis, there is currently a poor understanding of the genetic players and molecular mechanisms involved in the root architecture of apple trees under drought conditions. We previously observed that MdMYB88 and its paralog MdMYB124 regulate apple tree root morphology. In this study, we found that MdMYB88 and MdMYB124 play important roles in maintaining root hydraulic conductivity under long-term drought conditions and therefore contribute toward adaptive drought tolerance. Further investigation revealed that MdMYB88 and MdMYB124 regulate root xylem development by directly binding MdVND6 and MdMYB46 promoters and thus influence expression of their target genes under drought conditions. In addition, MdMYB88 and MdMYB124 were shown to regulate the deposition of cellulose and lignin root cell walls in response to drought. Taken together, our results provide novel insights into the importance of MdMYB88 and MdMYB124 in root architecture, root xylem development, and secondary cell wall deposition in response to drought in apple trees.

The quest for transcriptional hubs of lignin biosynthesis: beyond the NAC-MYB-gene regulatory network model

Lignin is an important secondary metabolite in plants. The biosynthesis of lignin is initiated by the transcriptional upregulation of genes encoding enzymes involved in monolignol biosynthesis and lignin polymerization. Based on studies of xylem differentiation over the last two decades, the NAC-MYB-based gene regulatory network (NAC-MYB-GRN) model is widely considered to underpin developmental lignin biosynthesis. We are now standing on the threshold of a new direction in transcriptional regulation research; the search for novel molecular hubs connecting developmental/environmental signals in lignin biosynthesis. Emerging genome-wide 'omics' technologies are a promising approach for understanding such hubs. Elucidating these molecular hubs may enable us to control lignification in harmony with plant development and environmental adaptation.

Ectopic Expression of PtoMYB74 in Poplar and Arabidopsis Promotes Secondary Cell Wall Formation

In vascular plants, R2R3-MYB transcription factors are important regulators of secondary cell wall formation. Although 192 annotated R2R3 MYB genes were identified in the poplar genome, only a few members were characterized to participate in the regulation of the secondary cell wall biosynthesis. In this paper, we identify an R2R3-MYB transcription factor, PtoMYB74, which is predicted to be an ortholog of Arabidopsis AtMYB61, a transcription activator that regulates the secondary cell wall formation, lignin biosynthesis, stomatal aperture, and the mucilage of seed coat. PtoMYB74 is mainly expressed in the stems, especially in the xylem tissues and organs. PtoMYB74 as a transcriptional activator is localized to the nucleus. The overexpression of PtoMYB74 increased the secondary cell wall thickness of vessels in transgenic poplar and changed the secondary cell wall compositions. The expression levels of lignin and cellulose biosynthetic genes were elevated in the transgenic poplar overexpressing PtoMYB74 compared to the wild type, while there was no change in the xylan biosynthetic genes. Transcriptional activation assays demonstrated that PtoMYB74 could activate the promoters of structural genes in the lignin and cellulose biosynthetic pathways. Taken together, our data show that PtoMYB74 positively regulates the secondary cell wall biosynthesis in poplar.

Sweet Sorghum Originated through Selection of Dry, a Plant-Specific NAC Transcription Factor Gene

Sorghum (Sorghum bicolor) is the fifth most popular crop worldwide and a C₄ model plant. Domesticated sorghum comes in many forms, including sweet cultivars with juicy stems and grain sorghum with dry, pithy stems at maturity. The Dry locus, which controls the pithy/juicy stem trait, was discovered over a century ago. Here, we found that Dry gene encodes a plant-specific NAC transcription factor. Dry was either deleted or acquired loss-of-function mutations in sweet sorghum, resulting in cell collapse and altered secondary cell wall composition in the stem. Twenty-three Dry ancestral haplotypes, all with dry, pithy stems, were found among wild sorghum and wild sorghum relatives. Two of the haplotypes were detected in domesticated landraces, with four additional dry haplotypes with juicy stems detected in improved lines. These results imply that selection for Dry gene mutations was a major step leading to the origin of sweet sorghum. The Dry gene is conserved in major cereals; fine-tuning its regulatory network could provide a molecular tool to control crop stem texture.

Plant cell wall sugars: sweeteners for a bio-based economy

Global warming and the consequent climate change is one of the major environmental challenges we are facing today. The driving force behind the rise in temperature is our fossil-based economy, which releases massive amounts of the greenhouse gas carbon dioxide into the atmosphere. In order to reduce greenhouse gas emission, we need to scale down our dependency on fossil resources, implying that we need other sources for energy and chemicals to feed our economy. Here, plants have an important role to play; by means of photosynthesis, plants capture solar energy to split water and fix carbon derived from atmospheric carbon dioxide. A significant fraction of the fixed carbon ends up as polysaccharides in the plant cell wall. Fermentable sugars derived from cell wall polysaccharides form an ideal carbon source for the production of bio-platform molecules. However, a major limiting factor in the use of plant biomass as feedstock for the bio-based economy is the complexity of the plant cell wall and its recalcitrance towards deconstruction. To facilitate the release of fermentable sugars during downstream biomass processing, the composition and structure of the cell wall can be engineered. Different strategies to reduce cell wall recalcitrance will be described in this review. The ultimate goal is to obtain a tailor-made biomass, derived from plants with a cell wall optimized for particular industrial or agricultural applications, without affecting plant growth and development.

Temporal analysis of Arabidopsis genes activated by Eucalyptus grandis NAC transcription factors associated with xylem fibre and vessel development

Secondary cell wall (SCW) deposition in Arabidopsis is regulated among others by NAC transcription factors, where SND1 chiefly initiates xylem fibre differentiation while VND6 controls metaxylem vessel SCW development, especially programmed cell death and wall patterning. The translational relevance of Arabidopsis SCW regulation theory and the utility of characterized transcription factors as modular synthetic biology tools for improving commercial fibre crops is unclear. We investigated inter-lineage gene activation dynamics for potential fibre and vessel differentiation regulators from the widely grown hardwood Eucalyptus grandis (Myrtales). EgrNAC26, a VND6 homolog, and EgrNAC61, an SND1 homolog, were transiently expressed in Arabidopsis mesophyll protoplasts in parallel to determine early and late (i.e. 7 and 14 hours post-transfection) gene targets. Surprisingly, across the time series EgrNAC26 activated only a subset of SCW-related transcription factors and biosynthetic genes activated by EgrNAC61, specializing instead in targeting vessel-specific wall pit and programmed cell death markers. Promoters of EgrNAC26 and EgrNAC61 both induced reporter gene expression in vessels of young Arabidopsis plants, with EgrNAC61 also conferring xylem- and cork cambium-preferential expression in Populus. Our results demonstrate partial conservation, with notable exceptions, of SND1 and VND6 homologs in Eucalyptus and a first report of cork cambium expression for EgrNAC61.

A sorghum NAC gene is associated with variation in biomass properties and yield potential

Sorghum bicolor is a C4 grass widely cultivated for grain, forage, sugar, and biomass. The sorghum Dry Stalk (D) locus controls a qualitative difference between juicy green (dd) and dry white (D-) stalks and midribs, and co-localizes with a quantitative trait locus for sugar yield. Here, we apply fine-mapping and genome-wide association study (GWAS) to identify a candidate gene underlying D, and use nearly isogenic lines (NILs) to characterize the transcriptional, compositional, and agronomic effects of variation at the D locus. The D locus was fine-mapped to a 36 kb interval containing four genes. One of these genes is a NAC transcription factor that contains a stop codon in the NAC domain in the recessive (dd) parent. Allelic variation at D affects grain yield, sugar yield, and biomass composition in NILs. Green midrib (dd) NILs show reductions in lignin in stalk tissue and produce higher sugar and grain yields under well-watered field conditions. Increased yield potential in dd NILs is associated with increased stalk mass and moisture, higher biomass digestibility, and an extended period of grain filling. Transcriptome profiling of midrib tissue at the 4-6 leaf stages, when NILs first become phenotypically distinct, reveals that dd NILs have increased expression of a miniature zinc finger (MIF) gene. MIF genes dimerize with and suppress zinc finger homeodomain (ZF-HD) transcription factors, and a ZF-HD gene is associated with midrib color variation in a GWAS analysis across 1,694 diverse sorghum inbreds. A premature stop codon in a NAC gene is the most likely candidate polymorphism underlying the sorghum D locus. More detailed understanding of the sorghum D locus could help improve agronomic potential in cereals.

V. AM, Solanke AU, Singh NK, Sharma TR, Gaikwad K. Genome-wide discovery of tissue-specific miRNAs in clusterbean (Cyamopsis tetragonoloba) indicates their association with galactomannan biosynthesis

Owing to the presence of 80% soluble dietary fibre, high protein content and high value gum, clusterbean (Cyamopsis tetragonoloba) has recently emerged as an economically important legume. The developing clusterbean seeds accumulate 90% galactomannans in the endosperm and, therefore, can be used as a model crop to understand galactomannan biosynthesis and its regulation. miRNAs are tiny master regulators of their corresponding target genes, resulting in variations in the amounts of their metabolic end products. To understand the role of these regulators in galactomannan biosynthesis regulation, small RNA libraries were prepared and sequenced from five tissues of clusterbean genotype RGC-936, and miRanalyzer and DSAP programs were used to identify conserved miRNAs and novel small RNAs. A total of 187 known and 171 novel miRNAs were found to be differentially expressed, of which 10 miRNAs were validated. A complicated network topology and 35% sharing of the target mRNAs between known and novel miRNAs suggest random evolution of novel miRNAs. The gene ontology (GO) annotation of potential target genes revealed the genes coding for signalling and carbohydrate metabolism (50.10%), kinases and other enzymes (20.75%), transcription factors (10.20%), transporters (8.35%) and other targets (10.6%). Two novel unigenes were annotated as ManS (mannosyltransferase/mannan synthase) and UGE (UDP- D-glucose 4-epimerase) and validated as targets for three novel miRNAs, that is Ct-miR3130, Ct-miR3135 and Ct-miR3157. Our findings reveal that these novel miRNAs could play an important role in the regulation of the galactomannan pathway in C. tetragonoloba and possibly other galactomannan-producing species.

Ethylene-Related Gene Expression Networks in Wood Formation

Thickening of tree stems is the result of secondary growth, accomplished by the meristematic activity of the vascular cambium. Secondary growth of the stem entails developmental cascades resulting in the formation of secondary phloem outwards and secondary xylem (i.e., wood) inwards of the stem. Signaling and transcriptional reprogramming by the phytohormone ethylene modifies cambial growth and cell differentiation, but the molecular link between ethylene and secondary growth remains unknown. We addressed this shortcoming by analyzing expression profiles and co-expression networks of ethylene pathway genes using the AspWood transcriptome database which covers all stages of secondary growth in aspen (Populus tremula) stems. ACC synthase expression suggests that the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) is synthesized during xylem expansion and xylem cell maturation. Ethylene-mediated transcriptional reprogramming occurs during all stages of secondary growth, as deduced from AspWood expression profiles of ethylene-responsive genes. A network centrality analysis of the AspWood dataset identified EIN3D and 11 ERFs as hubs. No overlap was found between the co-expressed genes of the EIN3 and ERF hubs, suggesting target diversification and hence independent roles for these transcription factor families during normal wood formation. The EIN3D hub was part of a large co-expression gene module, which contained 16 transcription factors, among them several new candidates that have not been earlier connected to wood formation and a VND-INTERACTING 2 (VNI2) homolog. We experimentally demonstrated Populus EIN3D function in ethylene signaling in Arabidopsis thaliana. The ERF hubs ERF118 and ERF119 were connected on the basis of their expression pattern and gene co-expression module composition to xylem cell expansion and secondary cell wall formation, respectively. We hereby establish data resources for ethylene-responsive genes and potential targets for EIN3D and ERF transcription factors in Populus stem tissues, which can help to understand the range of ethylene targeted biological processes during secondary growth.

QTL delineation for five fiber quality traits based on an intra-specific Gossypium hirsutum L. recombinant inbred line population

Gossypium hirsutum L. is the most important fiber crop worldwide and contributes to more than 95% of global cotton production. Marker-assisted selection (MAS) is an effective approach for improving fiber quality, and quantitative trait loci (QTL) mapping of fiber quality traits is important for cotton breeding. In this study, a permanent intra-specific recombinant inbred line (RIL) population containing 137 families was used for fiber quality testing. Based on a previously reported high-density genetic map with an average marker distance of 0.63 cM, 186 additive QTLs were obtained for five fiber quality traits over five consecutive years, including 39 for fiber length (FL), 36 for fiber strength (FS), 50 for fiber uniformity (FU), 33 for micronaire (MC) and 28 for fiber elongation (FE). Three stable QTLs, qMC-A4-1, qMC-D2-3 and qFS-D9-1, were detected in four datasets, and another eight stable QTLs, qMC-A4-2, qMC-D11-2, qFU-A9-1, qFU-A10-4, qFS-D11-1, qFL-D9-2, qFL-D11-1 and qFE-A3-2, were detected in three datasets. The annotated genes in these 11 stable QTLs were collected, and these genes included many transcription factors with functions during fiber development. 33 QTL coincidence regions were found, and these involved nearly half of the total QTLs. Four chromosome regions containing at least 6 QTLs were promising for fine mapping. In addition, 41 pairs of epistatic QTLs (e-QTLs) were screened, including 6 for FL, 30 for FS, 2 for FU and 3 for MC. The identification of stable QTLs adds valuable information for further QTL fine mapping and gene positional cloning for fiber quality genetic detection and provides useful markers for further molecular breeding in enhancing fiber quality.

Functional Genomics in the Study of Metabolic Pathways in Medicago truncatula: An Overview

In addition to its value as a model system for studies on symbiotic nitrogen fixation, Medicago truncatula has recently become an organism of choice for dissection of complex pathways of secondary metabolism. This work has been driven by two main reasons, both with practical implications. First Medicago species possess a wide range of flavonoid and terpenoid natural products, many of which, for example, the isoflavonoids and triterpene saponins, have important biological activities impacting both plant and animal (including human) health. Second, M. truncatula serves as an excellent model for alfalfa, the world's major forage legume, and forage quality is determined in large part by the concentrations of products of secondary metabolism, particularly lignin and condensed tannins. We here review recent progress in understanding the pathways leading to flavonoids, lignin, and triterpene saponins through utilization of genetic resources in M. truncatula.

Recent Advances in the Transcriptional Regulation of Secondary Cell Wall Biosynthesis in the Woody Plants

Plant cell walls provide structural support for growth and serve as a barrier for pathogen attack. Plant cell walls are also a source of renewable biomass for conversion to biofuels and bioproducts. Understanding plant cell wall biosynthesis and its regulation is of critical importance for the genetic modification of plant feedstocks for cost-effective biofuels and bioproducts conversion and production. Great progress has been made in identifying enzymes involved in plant cell wall biosynthesis, and in Arabidopsis it is generally recognized that the regulation of genes encoding these enzymes is under a transcriptional regulatory network with coherent feedforward and feedback loops. However, less is known about the transcriptional regulation of plant secondary cell wall (SCW) biosynthesis in woody species despite of its high relevance to biofuels and bioproducts conversion and production. In this article, we synthesize recent progress on the transcriptional regulation of SCW biosynthesis in Arabidopsis and contrast to what is known in woody species. Furthermore, we evaluate progress in related emerging regulatory machineries targeting transcription factors in this complex regulatory network of SCW biosynthesis.

Análisis in silico y expresión génica del factor de transcripción TgNAC01 implicado en xilogénesis y estrés abiótico en Tectona grandis

El xilema secundario es el componente más abundante de la biomasa vegetal. Por tanto, conocer los genes que regulan su formación ayudaría a diseñar estrategias para el mejoramiento genético de la madera. Así, el objetivo de este trabajo fue realizar el análisis computacional de la estructura primaria y secundaria del factor de transcripción (FT) TgNAC01 de Tectona grandis, además de evaluar su historia evolutiva, dominios conservados y expresión génica en tejidos lignificados de árboles de 12 y 60 años. Para ello, se realizó una evaluación del potencial de interacción ion-electrón (PIIE), mediante el método del espectro de la información (MEI) utilizando la librería SFAPS de R-Project, seguido del modelamiento estructural utilizando el software MODELLER y visualizado mediante PyMol. Además, el análisis de alineamiento de secuencia múltiple y filogenia fue mediante el software Bioedit y MrBayes respectivamente. También se evaluó los niveles de síntesis del FT TgNAC01 mediante qRT-PCR. Como resultados, se evidencio que el FT mantiene una estructura β-hoja antiparalela retorcida, que se compacta contra una α-hélice en la región N-terminal, teniendo así tres dominios α hélice y siete dominios β plegada. Asimismo, mediante el MEI se demostró que tiene alrededor de cinco funciones biológicas y mutaciones sobre los aminoácidos con mayor PIIE, lo que conlleva a evoluciones sobre las redes de regulación genética. Finalmente, el FT TgNAC01 juega un papel fundamental en la organización y desarrollo de las partes que componen la albura, como las células radiales de la zona cambial, los vasos, fibras y los anillos de crecimiento.

Tissue and cell-specific transcriptomes in cotton reveal the subtleties of gene regulation underlying the diversity of plant secondary cell walls

BACKGROUND

Knowledge of plant secondary cell wall (SCW) regulation and deposition is mainly based on the Arabidopsis model of a 'typical' lignocellulosic SCW. However, SCWs in other plants can vary from this. The SCW of mature cotton seed fibres is highly cellulosic and lacks lignification whereas xylem SCWs are lignocellulosic. We used cotton as a model to study different SCWs and the expression of the genes involved in their formation via RNA deep sequencing and chemical analysis of stem and seed fibre.

RESULTS

Transcriptome comparisons from cotton xylem and pith as well as from a developmental series of seed fibres revealed tissue-specific and developmentally regulated expression of several NAC transcription factors some of which are likely to be important as top tier regulators of SCW formation in xylem and/or seed fibre. A so far undescribed hierarchy was identified between the top tier NAC transcription factors SND1-like and NST1/2 in cotton. Key SCW MYB transcription factors, homologs of Arabidopsis MYB46/83, were practically absent in cotton stem xylem. Lack of expression of other lignin-specific MYBs in seed fibre relative to xylem could account for the lack of lignin deposition in seed fibre. Expression of a MYB103 homolog correlated with temporal expression of SCW CesAs and cellulose synthesis in seed fibres. FLAs were highly expressed and may be important structural components of seed fibre SCWs. Finally, we made the unexpected observation that cell walls in the pith of cotton stems contained lignin and had a higher S:G ratio than in xylem, despite that tissue's lacking many of the gene transcripts normally associated with lignin biosynthesis.

CONCLUSIONS

Our study in cotton confirmed some features of the currently accepted gene regulatory cascade for 'typical' plant SCWs, but also revealed substantial differences, especially with key downstream NACs and MYBs. The lignocellulosic SCW of cotton xylem appears to be achieved differently from that in Arabidopsis. Pith cell walls in cotton stems are compositionally very different from that reported for other plant species, including Arabidopsis. The current definition of a 'typical' primary or secondary cell wall might not be applicable to all cell types in all plant species.

Inferring Gene Regulatory Networks in the Arabidopsis Root Using a Dynamic Bayesian Network Approach

Gene regulatory network (GRN) models have been shown to predict and represent interactions among sets of genes. Here, we first show the basic steps to implement a simple but computationally efficient algorithm to infer GRNs based on dynamic Bayesian networks (DBNs), and we then explain how to approximate DBN-based GRN models with continuous models. In addition, we show a MATLAB implementation of the key steps of this method, which we use to infer an Arabidopsis root GRN.

RNA-seq Profiling Reveals Defense Responses in a Tolerant Potato Cultivar to Stem Infection by Pectobacterium carotovorum ssp. brasiliense

Pectobacterium carotovorum subsp. brasiliense is a member of the soft rot Enterobacteriaceae (SRE) family that causes tuber soft rot and blackleg diseases of stems in potato plants. Currently, there are no effective chemical strategies for the control of members of the SRE. Thus, an understanding of the inducible defense responses in stems of potato plants is important, particularly during colonization of the vascular system. Here, time-course RNA-sequencing analysis was used to compare expressed genes between a susceptible potato cultivar (Solanum tuberosum cv Valor) and a tolerant cultivar (S. tuberosum cv BP1) at 0, 6, 12, 24, and 72 h post-inoculation with P. c. brasiliense. In total, we identified 6139 and 8214 differentially expressed genes (DEGs) in the tolerant and susceptible cultivars, compared to mock-inoculated controls, respectively. Key DEGs distinguishing between tolerance and susceptibility were associated with negative regulation of cell death and plant-type cell wall organization/biogenesis biological processes in the tolerant and susceptible cultivars, respectively. Among these were DEGs involved in signaling (mainly MAPK cascade and ethylene pathway), defense-related transcription regulation including WRKY transcription factors, and downstream secondary cell biosynthesis. Together, our results suggest that S. tuberosum cv BP1 likely employs quantitative defense response against P. c. brasiliense. Overall, our study provides the first transcriptome-wide insight into the molecular basis of tolerance and/or resistance of potato stems to SRE infection.