Cuticular wax biosynthesis is up-regulated by the MYB94 transcription factor in Arabidopsis

Transcription factors and plant hormones mediate wax metabolism in response to drought stress

Plants have, throughout evolution, developed a hydrophobic cuticle to protect them from various stresses in the terrestrial environment. The cuticle layer is mainly composed of cutin and cuticular wax, a mixture of very-long-chain fatty acids and their derivatives. With the progress of transcriptome sequencing and other research methods, the key enzymes, transporters and regulatory factors in wax synthesis and metabolism have been gradually identified, especially the study on the regulation of wax metabolism by transcription factors and others in response to plant stress has become a hot topic. Drought is a major abiotic stress that limits plant growth and crop productivity. Plant epidermal wax prevents non-stomatal water loss and improves water use efficiency to adapt to arid environments. In this study, the ways of wax synthesis, transport, metabolism and regulation at different levels are reviewed. At the same time, the regulation of wax by different transcription factors and plant hormones in response to drought is elaborated, and key research questions and important directions for future solutions are proposed to enhance the potential application of epidermal wax in agriculture and the environment.

Wounding Triggers Wax Biosynthesis in Arabidopsis Leaves in an Abscisic Acid-Dependent and Jasmonoyl-Isoleucine-Dependent Manner

Wounding caused by insects or abiotic factors such as wind and hail can cause severe stress for plants. Intrigued by the observation that wounding induces expression of genes involved in surface wax synthesis in a jasmonoyl-isoleucine (JA-Ile)-independent manner, the role of wax biosynthesis and respective genes upon wounding was investigated. Wax, a lipid-based barrier, protects plants both from environmental threats and from an uncontrolled loss of water. Its biosynthesis is described to be regulated by abscisic acid (ABA), whereas the main wound signal is the hormone JA-Ile. We show in this study that genes coding for enzymes of surface wax synthesis are induced upon wounding in Arabidopsis thaliana leaves in a JA-Ile-independent but an ABA-dependent manner. Furthermore, the ABA-dependent transcription factor MYB96 is a key regulator of wax biosynthesis upon wounding. On the metabolite level, wound-induced wax accumulation is strongly reduced in JA-Ile-deficient plants, but this induction is only slightly decreased in ABA-reduced plants. To further analyze the ABA-dependent wound response, we conducted wounding experiments in high humidity. They show that high humidity prevents the wound-induced wax accumulation in A. thaliana leaves. Together the data presented in this study show that wound-induced wax accumulation is JA-Ile-dependent on the metabolite level, but the expression of genes coding for enzymes of wax synthesis is regulated by ABA.

Wheat MIXTA-like Transcriptional Activators Positively Regulate Cuticular Wax Accumulation

MIXTA-like transcription factors AtMYB16 and AtMYB106 play important roles in the regulation of cuticular wax accumulation in dicot model plant Arabidopsis thaliana, but there are very few studies on the MIXTA-like transcription factors in monocot plants. Herein, wheat MIXTA-like transcription factors TaMIXTA1 and TaMIXTA2 were characterized as positive regulators of cuticular wax accumulation. The virus-induced gene silencing experiments showed that knock-down of wheat TaMIXTA1 and TaMIXTA2 expressions resulted in the decreased accumulation of leaf cuticular wax, increased leaf water loss rate, and potentiated chlorophyll leaching. Furthermore, three wheat orthologous genes of ECERIFERUM 5 (TaCER5-1A, 1B, and 1D) and their function in cuticular wax deposition were reported. The silencing of TaCER5 by BSMV-VIGS led to reduced loads of leaf cuticular wax and enhanced rates of leaf water loss and chlorophyll leaching, indicating the essential role of the TaCER5 gene in the deposition of wheat cuticular wax. In addition, we demonstrated that TaMIXTA1 and TaMIXTA2 function as transcriptional activators and could directly stimulate the transcription of wax biosynthesis gene TaKCS1 and wax deposition gene TaCER5. The above results strongly support that wheat MIXTA-Like transcriptional activators TaMIXTA1 and TaMIXTA2 positively regulate cuticular wax accumulation via activating TaKCS1 and TaCER5 gene transcription.

AtMYB41 acts as a dual-function transcription factor that regulates the formation of lipids in an organ- and development-dependent manner

The plant cuticle controls non-stomatal water loss and can serve as a barrier against biotic agents, whereas the heteropolymer suberin and its associated waxes are deposited constitutively at specific cell wall locations. While several transcription factors controlling cuticle formation have been identified, those involved in the transcriptional regulation of suberin biosynthesis remain poorly characterized. The major goal of this study was to further analyse the function of the R2R3-Myeloblastosis (MYB) transcription factor AtMYB41 in formation of the cuticle, suberin, and suberin-associated waxes throughout plant development. For functional analysis, the organ-specific expression pattern of AtMYB41 was analysed and Atmyb41ge alleles were generated using the CRISPR/Cas9 system. These were investigated for root growth and water permeability upon stress. In addition, the fatty acid, wax, cutin, and suberin monomer composition of different organs was evaluated by gas chromatography. The characterization of Atmyb41ge mutants revealed that AtMYB41 negatively regulates the production of cuticular lipids and fatty acid biosynthesis in leaves and seeds, respectively. Remarkably, biochemical analyses indicate that AtMYB41 also positively regulates the formation of cuticular waxes in stems of Arabidopsis thaliana. Overall, these results suggest that the AtMYB41 acts as a negative regulator of cuticle and fatty acid biosynthesis in leaves and seeds, respectively, but also as a positive regulator of wax production in A. thaliana stems.

Deciphering variations, identification of marker-trait associations and candidate genes for seed oil content under terminal heat stress in Indian mustard (Brassica juncea L. Czern & Coss) germplasm stock

This research paper investigates the variability in seed oil content (SOC) in Indian mustard (Brassica juncea L.) under terminal heat stress (THS) conditions. A genetic stock of 488 genotypes of B. juncea was evaluated over two years and grouped into five classes based on the reduction in oil content under THS compared to normal sown crop. Based on heat susceptibility index (HSI), a diverse panel of 96 genotypes was selected and evaluated under THS. Twenty-two heat-tolerant donor genotypes were identified, including introgression lines derived from B. tournefortii, B. carinata and Erucastrum cardaminoides. This study is the first to report on marker-trait associations for SOC in B. juncea under THS using a GWAS approach. Furthermore, candidate genes associated with abiotic stress tolerance and lipid metabolism were identified near the significant SNPs, emphasizing their role in SOC regulation under stress. Notable candidate genes include BjuA003240 (encoding for alcohol-forming fatty acyl-CoA reductase), BjuA003242 (involving in lipid biosynthesis), BjuA003244 (associated with mitochondrial functions and stress tolerance), and BjuA003245 (related to MYB transcription factors regulating lipid biosynthesis). This study provides valuable insights into the genetic basis of SOC variation under THS in B. juncea, highlighting potential breeding targets for improved heat stress resilience in Indian mustard cultivation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03985-w.

Characterization of the Fatty Acyl-CoA Reductase (FAR) Gene Family and Its Response to Abiotic Stress in Rice (Oryza sativa L.)

Fatty acyl-CoA reductase (FAR) is an important NADPH-dependent enzyme that can produce primary alcohol from fatty acyl-CoA or fatty acyl-carrier proteins as substrates. It plays a pivotal role in plant growth, development, and stress resistance. Herein, we performed genome-wide identification and expression analysis of FAR members in rice using bioinformatics methods. A total of eight OsFAR genes were identified, and the OsFARs were comprehensively analyzed in terms of phylogenetic relationships, duplication events, protein motifs, etc. The cis-elements of the OsFARs were predicted to respond to growth and development, light, hormones, and abiotic stresses. Gene ontology annotation analysis revealed that OsFAR proteins participate in biological processes as fatty acyl-CoA reductase during lipid metabolism. Numerous microRNA target sites were present in OsFARs mRNAs. The expression analysis showed that OsFARs were expressed at different levels during different developmental periods and in various tissues. Furthermore, the expression levels of OsFARs were altered under abiotic stresses, suggesting that FARs may be involved in abiotic stress tolerance in rice. The findings presented here serve as a solid basis for further exploring the functions of OsFARs.

A genome-wide association study identifies genes associated with cuticular wax metabolism in maize

The plant cuticle is essential in plant defense against biotic and abiotic stresses. To systematically elucidate the genetic architecture of maize (Zea mays L.) cuticular wax metabolism, 2 cuticular wax-related traits, the chlorophyll extraction rate (CER) and water loss rate (WLR) of 389 maize inbred lines, were investigated and a genome-wide association study (GWAS) was performed using 1.25 million single nucleotide polymorphisms (SNPs). In total, 57 nonredundant quantitative trait loci (QTL) explaining 5.57% to 15.07% of the phenotypic variation for each QTL were identified. These QTLs contained 183 genes, among which 21 strong candidates were identified based on functional annotations and previous publications. Remarkably, 3 candidate genes that express differentially during cuticle development encode β-ketoacyl-CoA synthase (KCS). While ZmKCS19 was known to be involved in cuticle wax metabolism, ZmKCS12 and ZmKCS3 functions were not reported. The association between ZmKCS12 and WLR was confirmed by resequencing 106 inbred lines, and the variation of WLR was significant between different haplotypes of ZmKCS12. In this study, the loss-of-function mutant of ZmKCS12 exhibited wrinkled leaf morphology, altered wax crystal morphology, and decreased C32 wax monomer levels, causing an increased WLR and sensitivity to drought. These results confirm that ZmKCS12 plays a vital role in maize C32 wax monomer synthesis and is critical for drought tolerance. In sum, through GWAS of 2 cuticular wax-associated traits, this study reveals comprehensively the genetic architecture in maize cuticular wax metabolism and provides a valuable reference for the genetic improvement of stress tolerance in maize.

MdDEWAX decreases plant drought resistance by regulating wax biosynthesis

Apple epidermal wax protects plants from environmental stresses, determines fruit gloss and improves postharvest storage quality. However, the molecular mechanisms underlying the biosynthesis and regulation of apple epidermal waxes are not fully understood. In this study, we isolated a MdDEWAX gene from apple, which localized in the nucleus, expressed mainly in apple fruit, and induced by drought. We transformed the MdDEWAX gene into Arabidopsis, and found that heterologous expression of MdDEWAX reduced the accumulation of cuticular waxes in leaves and stems, increased epidermal permeability, the rate of water loss, and the rate of chlorophyll extraction of leaves and stems, altered the sensitivity to ABA, and reduced drought tolerance. Meanwhile, overexpression or silencing of the gene in the epidermis of apple fruits decreased or increased wax content, respectively. This study provides candidate genes for breeding apple cultivars and rootstocks with better drought tolerance.

Genome-Wide Analysis of MYB Genes in Primulina eburnea (Hance) and Identification of Members in Response to Drought Stress

Due to periodic water deficiency in karst environments, Primulina eburnea experiences sporadic drought stress in its habitat. Despite being one of the largest gene families and functionally diverse in terms of plant growth and development, MYB transcription factors in P. eburnea have not been studied. Here, a total of 230 MYB genes were identified in P. eburnea, including 67 1R-MYB, 155 R2R3-MYB, six 3R-MYB, and two 4R-MYB genes. The R2R3-type PebMYB genes could be classified into 16 subgroups, while the remaining PebMYB genes (1R-MYB, 3R-MYB, and 4R-MYB genes) were divided into 10 subgroups. Notably, the results of the phylogenetic analysis were further supported by the motif and gene structure analysis, which showed that individuals in the same subgroup had comparable motif and structure organization. Additionally, gene duplication and synteny analyses were performed to better understand the evolution of PebMYB genes, and 291 pairs of segmental duplicated genes were found. Moreover, RNA-seq analysis revealed that the PebMYB genes could be divided into five groups based on their expression characteristics. Furthermore, 11 PebMYB genes that may be involved in drought stress response were identified through comparative analysis with Arabidopsis thaliana. Notably, seven of these genes (PebMYB3, PebMYB13, PebMYB17, PebMYB51, PebMYB142, PebMYB69, and PebMYB95) exhibited significant differences in expression between the control and drought stress treatments, suggesting that they may play important roles in drought stress response. These findings clarified the characteristics of the MYB gene family in P. eburnea, augmenting our comprehension of their potential roles in drought stress adaptation.

The R2R3-MYB transcription factor EVER controls the emission of petunia floral volatiles by regulating epicuticular wax biosynthesis in the petal epidermis

The epidermal cells of petunia (Petunia × hybrida) flowers are the main site of volatile emission. However, the mechanisms underlying the release of volatiles into the environment are still being explored. Here, using cell-layer-specific transcriptomic analysis, reverse genetics by virus-induced gene silencing and clustered regularly interspaced short palindromic repeat (CRISPR), and metabolomics, we identified EPIDERMIS VOLATILE EMISSION REGULATOR (EVER)-a petal adaxial epidermis-specific MYB activator that affects the emission of volatiles. To generate ever knockout lines, we developed a viral-based CRISPR/Cas9 system for efficient gene editing in plants. These knockout lines, together with transient-suppression assays, revealed EVER's involvement in the repression of low-vapor-pressure volatiles. Internal pools and annotated scent-related genes involved in volatile production and emission were not affected by EVER. RNA-Seq analyses of petals of ever knockout lines and EVER-overexpressing flowers revealed enrichment in wax-related biosynthesis genes. Liquid chromatography/gas chromatography-MS analyses of petal epicuticular waxes revealed substantial reductions in wax loads in ever petals, particularly of monomers of fatty acids and wax esters. These results implicate EVER in the emission of volatiles by fine-tuning the composition of petal epicuticular waxes. We reveal a petunia MYB regulator that interlinks epicuticular wax composition and volatile emission, thus unraveling a regulatory layer in the scent-emission machinery in petunia flowers.

Genome-wide identification and expression analysis of the KCS gene family in soybean (Glycine max) reveal their potential roles in response to abiotic stress

Very long chain fatty acids (VLCFAs) are fatty acids with chain lengths of 20 or more carbon atoms, which are the building blocks of various lipids that regulate developmental processes and plant stress responses. 3-ketoacyl-CoA synthase encoded by the KCS gene is the key rate-limiting enzyme in VLCFA biosynthesis, but the KCS gene family in soybean (Glycine max) has not been adequately studied thus far. In this study, 31 KCS genes (namely GmKCS1 - GmKCS31) were identified in the soybean genome, which are unevenly distributed on 14 chromosomes. These GmKCS genes could be phylogenetically classified into seven groups. A total of 27 paralogous GmKCS gene pairs were identified with their Ka/Ks ratios indicating that they had undergone purifying selection during soybean genome expansion. Cis-acting element analysis revealed that GmKCS promoters contained multiple hormone- and stress-responsive elements, indicating that GmKCS gene expression levels may be regulated by various developmental and environmental stimuli. Expression profiles derived from RNA-seq data and qRT-PCR experiments indicated that GmKCS genes were diversely expressed in different organs/tissues, and many GmKCS genes were found to be differentially expressed in the leaves under cold, heat, salt, and drought stresses, suggesting their critical role in soybean resistance to abiotic stress. These results provide fundamental information about the soybean KCS genes and will aid in their further functional elucidation and exploitation.

The Accumulation of Volatile Compounds and the Change in the Morphology of the Leaf Wax Cover Accompanied the "Anti-Aging" Effect in Anethum graveolens L. Plants Sprayed with 6-Benzylaminopurine

Essential oils (EOs) are of commercial importance for medicine, food, cosmetics, the perfume industry, and agriculture. In plants, EOs, like the wax cover, serve as protection against abiotic stresses, such as high temperatures and water deficiency. The use of spraying with exogenous hormones of aromatic plants affects the accumulation and composition of volatile compounds, as well as tolerance to abiotic stress. As a result of cytokinin treatment with 6-BAP (6-benzylaminopurine) (200 mg L-l) of Anetum graveolens L. "Uzory" and "Rusich" varieties, several responses to its action were revealed: a change in the division of leaf blades, inhibition of flowering, an increase in the content of EO and its main components α-phellandrene and p-cymene in leaves, and limonene in umbels and fruits. It was revealed that the increased accumulation of EO in dill leaves was longer with sufficient moisture. In contrast, under conditions of heat and water deficiency, the effect of 6-BAP treatment on accumulations of the EO in leaves was short-lived and did not appear on umbels and fruits. The study of the cytokinin effect on a fine structure of a wax cover on the adaxial side of leaves by scanning electron microscopy revealed a change in its elements (from amorphous layers with scales to thin tubules), which probably increased the sensitivity of leaves to water deficiency and, consequently, led to a decrease in the biosynthetic activity of leaf tissue. Thus, 6-BAP had an impact on the adaptive properties of dill plants, prolonging the "youth" of vegetative organs and the ability to EO biosynthesis under conditions of sufficient moisture.

Time course of changes in the transcriptome during russet induction in apple fruit

BACKGROUND

Russeting is a major problem in many fruit crops. Russeting is caused by environmental factors such as wounding or moisture exposure of the fruit surface. Despite extensive research, the molecular sequence that triggers russet initiation remains unclear. Here, we present high-resolution transcriptomic data by controlled russet induction at very early stages of fruit development. During Phase I, a patch of the fruit surface is exposed to surface moisture. For Phase II, moisture exposure is terminated, and the formerly exposed surface remains dry. We targeted differentially expressed transcripts as soon as 24 h after russet induction.

RESULTS

During moisture exposure (Phase I) of 'Pinova' apple, transcripts associated with the cell cycle, cell wall, and cuticle synthesis (SHN3) decrease, while those related to abiotic stress increase. NAC35 and MYB17 were the earliest induced genes during Phase I. They are therefore linked to the initial processes of cuticle microcracking. After moisture removal (Phase II), the expression of genes related to meristematic activity increased (WOX4 within 24 h, MYB84 within 48 h). Genes related to lignin synthesis (MYB52) and suberin synthesis (MYB93, WRKY56) were upregulated within 3 d after moisture removal. WOX4 and AP2B3 are the earliest differentially expressed genes induced in Phase II. They are therefore linked to early events in periderm formation. The expression profiles were consistent between two different seasons and mirrored differences in russet susceptibility in a comparison of cultivars. Furthermore, expression profiles during Phase II of moisture induction were largely identical to those following wounding.

CONCLUSIONS

The combination of a unique controlled russet induction technique with high-resolution transcriptomic data allowed for the very first time to analyse the formation of cuticular microcracks and periderm in apple fruit immediately after the onset of triggering factors. This data provides valuable insights into the spatial-temporal dynamics of russeting, including the synthesis of cuticles, dedifferentiation of cells, and impregnation of cell walls with suberin and lignin.

The Over-Expression of Two R2R3-MYB Genes, PdMYB2R089 and PdMYB2R151, Increases the Drought-Resistant Capacity of Transgenic Arabidopsis

The R2R3-MYB genes in plants play an essential role in the drought-responsive signaling pathway. Plenty of R2R3-MYB S21 and S22 subgroup genes in Arabidopsis have been implicated in dehydration conditions, yet few have been covered in terms of the role of the S21 and S22 subgroup genes in poplar under drought. PdMYB2R089 and PdMYB2R151 genes, respectively belonging to the S21 and S22 subgroups of NL895 (Populus deltoides × P. euramericana cv. 'Nanlin895'), were selected based on the previous expression analysis of poplar R2R3-MYB genes that are responsive to dehydration. The regulatory functions of two target genes in plant responses to drought stress were studied and speculated through the genetic transformation of Arabidopsis thaliana. PdMYB2R089 and PdMYB2R151 could promote the closure of stomata in leaves, lessen the production of malondialdehyde (MDA), enhance the activity of the peroxidase (POD) enzyme, and shorten the life cycle of transgenic plants, in part owing to their similar conserved domains. Moreover, PdMYB2R089 could strengthen root length and lateral root growth. These results suggest that PdMYB2R089 and PdMYB2R151 genes might have the potential to improve drought adaptability in plants. In addition, PdMYB2R151 could significantly improve the seed germination rate of transgenic Arabidopsis, but PdMYB2R089 could not. This finding provides a clue for the subsequent functional dissection of S21 and S22 subgroup genes in poplar that is responsive to drought.

Roles of the MYB94/FUSED LEAVES1 (ZmFDL1) and GLOSSY2 (ZmGL2) genes in cuticle biosynthesis and potential impacts on Fusarium verticillioides growth on maize silks

Maize silks, the stigmatic portions of the female flowers, have an important role in reproductive development. Silks also provide entry points for pathogens into host tissues since fungal hyphae move along the surface of the silks to reach the site of infection, i.e., the developing kernel. The outer extracellular surface of the silk is covered by a protective hydrophobic cuticle, comprised of a complex array of long-chain hydrocarbons and small amounts of very long chain fatty acids and fatty alcohols. This work illustrates that two previously characterized cuticle-related genes separately exert roles on maize silk cuticle deposition and function. ZmMYB94/FUSED LEAVES 1 (ZmFDL1) MYB transcription factor is a key regulator of cuticle deposition in maize seedlings. The ZmGLOSSY2 (ZmGL2) gene, a putative member of the BAHD superfamily of acyltransferases with close sequence similarity to the Arabidopsis AtCER2 gene, is involved in the elongation of the fatty acid chains that serve as precursors of the waxes on young leaves. In silks, lack of ZmFDL1 action generates a decrease in the accumulation of a wide number of compounds, including alkanes and alkenes of 20 carbons or greater and affects the expression of cuticle-related genes. These results suggest that ZmFDL1 retains a regulatory role in silks, which might be exerted across the entire wax biosynthesis pathway. Separately, a comparison between gl2-ref and wild-type silks reveals differences in the abundance of specific cuticular wax constituents, particularly those of longer unsaturated carbon chain lengths. The inferred role of ZmGL2 is to control the chain lengths of unsaturated hydrocarbons. The treatment of maize silks with Fusarium verticillioides conidia suspension results in altered transcript levels of ZmFDL1 and ZmGL2 genes. In addition, an increase in fungal growth was observed on gl2-ref mutant silks 72 hours after Fusarium infection. These findings suggest that the silk cuticle plays an active role in the response to F. verticillioides infection.

Transcription Factor TaMYB30 Activates Wheat Wax Biosynthesis

The waxy cuticle covers a plant's aerial surface and contributes to environmental adaptation in land plants. Although past decades have seen great advances in understanding wax biosynthesis in model plants, the mechanisms underlying wax biosynthesis in crop plants such as bread wheat remain to be elucidated. In this study, wheat MYB transcription factor TaMYB30 was identified as a transcriptional activator positively regulating wheat wax biosynthesis. The knockdown of TaMYB30 expression using virus-induced gene silencing led to attenuated wax accumulation, increased water loss rates, and enhanced chlorophyll leaching. Furthermore, TaKCS1 and TaECR were isolated as essential components of wax biosynthetic machinery in bread wheat. In addition, silencing TaKCS1 and TaECR resulted in compromised wax biosynthesis and potentiated cuticle permeability. Importantly, we showed that TaMYB30 could directly bind to the promoter regions of TaKCS1 and TaECR genes by recognizing the MBS and Motif 1 cis-elements, and activate their expressions. These results collectively demonstrated that TaMYB30 positively regulates wheat wax biosynthesis presumably via the transcriptional activation of TaKCS1 and TaECR.

Wound-induced triacylglycerol biosynthesis is jasmonoy-l-isoleucin and abscisic acid independent

Triacylglycerol (TAG) plays a significant role during plant stress - it maintains lipid homeostasis. Upon wounding plants accumulate TAG, likely as a storage form of fatty acids (FAs) that originate from damaged membranes. This study asked if this process depends on the two phytohormones jasmonoyl-isoleucine (JA-Ile) and abscisic acid (ABA), which are involved in wound signalling. To analyse regulation of wound-induced TAG accumulation, we used mutants deficient in JA-Ile, with reduced ABA and the myb96 mutant, which is deficient in an ABA-dependent transcription factor. The expression of genes involved in TAG biosynthesis, and TAG content after wounding were analysed via LC-MS and GC-FID, plastidial lipid content in all mentioned mutant lines was also determined. The localization of newly synthesized TAG was investigated using lipid droplet staining. TAG accumulation upon wounding was confirmed as well as the fact that the newly synthesized TAG are mostly composed of polyunsaturated fatty acids. Nevertheless, all tested mutant lines were able to accumulate TAG similar to the WT. We observed differences in reduction of plastidial lipids - in WT plants this was higher than in mutant lines. Newly synthesized TAGs were stored in lipid droplets at and around the wounded area. Our results show that TAG accumulation upon wounding is not dependent on JA-Ile or ABA. The newly synthesized TAG species are composed of unsaturated fatty acids of membrane origin, and most likely serves as a transient energy store.

Comprehensive Genome-Wide Analyses of Poplar R2R3-MYB Transcription Factors and Tissue-Specific Expression Patterns under Drought Stress

R2R3-type MYB transcription factors are implicated in drought stress, which is a primary factor limiting the growth and development of woody plants. The identification of R2R3-MYB genes in the Populus trichocarpa genome has been previously reported. Nevertheless, the diversity and complexity of the conserved domain of the MYB gene caused inconsistencies in these identification results. There is still a lack of drought-responsive expression patterns and functional studies of R2R3-MYB transcription factors in Populus species. In this study, we identified a total of 210 R2R3-MYB genes in the P. trichocarpa genome, of which 207 genes were unevenly distributed across all 19 chromosomes. These poplar R2R3-MYB genes were phylogenetically divided into 23 subgroups. Collinear analysis demonstrated that the poplar R2R3-MYB genes underwent rapid expansion and that whole-genome duplication events were a dominant factor in the process of rapid gene expansion. Subcellular localization assays indicated that poplar R2R3-MYB TFs mainly played a transcriptional regulatory role in the nucleus. Ten R2R3-MYB genes were cloned from P. deltoides × P. euramericana cv. Nanlin895, and their expression patterns were tissue-specific. A majority of the genes showed similar drought-responsive expression patterns in two out of three tissues. This study provides a valid cue for further functional characterization of drought-responsive R2R3-MYB genes in poplar and provides support for the development of new poplar genotypes with elevated drought tolerance.

Genome-wide analysis of R2R3-MYB transcription factors reveals their differential responses to drought stress and ABA treatment in desert poplar (Populus euphratica)

The R2R3-MYB transcription factors are widely involved in the regulation of plant growth, biotic and abiotic stress responses. Meanwhile, seed germination, which is stimulated by internal and external environments, is a critical stage in the plant life cycle. However, the identification, characterization, and expression profiling of the Populus euphratica R2R3-MYB family in drought response during seed germination have been unknown. Our study attempted to identify and characterize the R2R3-MYB genes in P. euphratica (PeR2R3-MYBs) and explore how R2R3-MYBs trigger the drought and abscisic acid (ABA) response mechanism in its seedlings. Based on the analysis of comparative genomics, 174 PeR2R3-MYBs were identified and expanded driven by whole genome duplication or segment duplication events. The analysis of Ka/Ks ratios showed that, in contrast to most PeR2R3-MYBs, the other PeR2R3-MYBs were subjected to positive selection in P. euphratica. Further, the expression data of PeR2R3-MYBs under drought stress and ABA treatment, together with available functional data for Arabidopsis thaliana MYB genes, supported the hypothesis that PeR2R3-MYBs involved in response to drought are dependent or independent on ABA signaling pathway during seed germination, especially PeR2R3-MYBs with MYB binding sites (MBS) cis-element and/or tandem duplication. This study is the first report on the genome-wide analysis of PeR2R3-MYBs, as well as the other two Salicaceae species. The duplication events and differential expressions of PeR2R3-MYBs play important roles in enhancing the adaptation to drought desert environment. Our results provide a reference for prospective functional studies of R2R3-MYBs of poplars and lay the foundation for new breeding strategies to improve the drought tolerance of P. euphratica.

Genome-Wide Identification and Transcriptional Analysis of the MYB Gene Family in Pearl Millet (Pennisetum glaucum)

The MYB gene family widely exists in the plant kingdom and participates in the regulation of plant development and stress response. Pearl millet (Pennisetum glaucum (L.) R. Br.), as one of the most important cereals, is not only considered a good source of protein and nutrients but also has excellent tolerances to various abiotic stresses (e.g., salinity, water deficit, etc.). Although the genome sequence of pearl millet was recently published, bioinformatics and expression pattern analysis of the MYB gene family are limited. Here, we identified 208 PgMYB genes in the pearl millet genome and employed 193 high-confidence candidates for downstream analysis. Phylogenetic and structural analysis classified these PgMYBs into four subgroups. Eighteen pairs of segmental duplications of the PgMYB gene were found using synteny analysis. Collinear analysis revealed pearl millet had the closest evolutionary relationship with foxtail millet. Nucleotide substitution analysis (Ka/Ks) revealed PgMYB genes were under purifying positive selection pressure. Reverse transcription-quantitative PCR analysis of eleven R2R3-type PgMYB genes revealed they were preferentially expressed in shoots and seeds and actively responded to various environment stimuli. Current results provide insightful information regarding the molecular features of the MYB family in pearl millet to support further functional characterizations.

The Plant Fatty Acyl Reductases

Fatty acyl reductase (FAR) is a crucial enzyme that catalyzes the NADPH-dependent reduction of fatty acyl-CoA or acyl-ACP substrates to primary fatty alcohols, which in turn acts as intermediate metabolites or metabolic end products to participate in the formation of plant extracellular lipid protective barriers (e.g., cuticular wax, sporopollenin, suberin, and taproot wax). FARs are widely present across plant evolution processes and play conserved roles during lipid synthesis. In this review, we provide a comprehensive view of FAR family enzymes, including phylogenetic analysis, conserved structural domains, substrate specificity, subcellular localization, tissue-specific expression patterns, their varied functions in lipid biosynthesis, and the regulation mechanism of FAR activity. Finally, we pose several questions to be addressed, such as the roles of FARs in tryphine, the interactions between transcription factors (TFs) and FARs in various environments, and the identification of post-transcriptional, translational, and post-translational regulators.

Domestication and selection footprints in Persian walnuts (Juglans regia)

Walnut (Juglans) species are economically important hardwood trees cultivated worldwide for both edible nuts and high-quality wood. Broad-scale assessments of species diversity, evolutionary history, and domestication are needed to improve walnut breeding. In this study, we sequenced 309 walnut accessions from around the world, including 55 Juglans relatives, 98 wild Persian walnuts (J. regia), 70 J. regia landraces, and 86 J. regia cultivars. The phylogenetic tree indicated that J. regia samples (section Dioscaryon) were monophyletic within Juglans. The core areas of genetic diversity of J. regia germplasm were southwestern China and southern Asia near the Qinghai-Tibet Plateau and the Himalayas, and the uplift of the Himalayas was speculated to be the main factor leading to the current population dynamics of Persian walnut. The pattern of genomic variation in terms of nucleotide diversity, linkage disequilibrium, single nucleotide polymorphisms, and insertions/deletions revealed the domestication and selection footprints in Persian walnut. Selective sweep analysis, GWAS, and expression analysis further identified two transcription factors, JrbHLH and JrMYB6, that influence the thickness of the nut diaphragm as loci under selection during domestication. Our results elucidate the domestication and selection footprints in Persian walnuts and provide a valuable resource for the genomics-assisted breeding of this important crop.

Transcriptome Analysis Reveals a Comprehensive Virus Resistance Response Mechanism in Pecan Infected by a Novel Badnavirus Pecan Virus

Pecan leaf-variegated plant, which was infected with a novel badnavirus named pecan mosaic virus (PMV) detected by small RNA deep sequencing, is a vital model plant for studying the molecular mechanism of retaining green or chlorosis of virus-infected leaves. In this report, PMV infection in pecan leaves induced PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI). PMV infection suppressed the expressions of key genes of fatty acid, oleic acid (C18:1), and very-long-chain fatty acids (VLCFA) biosynthesis, indicating that fatty acids-derived signaling was one of the important defense pathways in response to PMV infection in pecan. PMV infection in pecans enhanced the expressions of pathogenesis-related protein 1 (PR1). However, the transcripts of phenylalanine ammonia-lyase (PAL) and isochorismate synthase (ICS) were downregulated, indicating that salicylic acid (SA) biosynthesis was blocked in pecan infected with PMV. Meanwhile, disruption of auxin signaling affected the activation of the jasmonic acid (JA) pathway. Thus, C18:1 and JA signals are involved in response to PMV infection in pecan. In PMV-infected yellow leaves, damaged chloroplast structure and activation of mitogen-activated protein kinase 3 (MPK3) inhibited photosynthesis. Cytokinin and SA biosynthesis was blocked, leading to plants losing immune responses and systemic acquired resistance (SAR). The repression of photosynthesis and the induction of sink metabolism in the infected tissue led to dramatic changes in carbohydrate partitioning. On the contrary, the green leaves of PMV infection in pecan plants had whole cell tissue structure and chloroplast clustering, establishing a strong antiviral immunity system. Cytokinin biosynthesis and signaling transductions were remarkably strengthened, activating plant immune responses. Meanwhile, cytokinin accumulation in green leaves induced partial SA biosynthesis and gained comparatively higher SAR compared to that of yellow leaves. Disturbance of the ribosome biogenesis might enhance the resistance to PMV infection in pecan and lead to leaves staying green.

MYB44-ENAP1/2 restricts HDT4 to regulate drought tolerance in Arabidopsis

Histone acetylation has been shown to involve in stress responses. However, the detailed molecular mechanisms that how histone deacetylases and transcription factors function in drought stress response remain to be understood. In this research, we show that ENAP1 and ENAP2 are positive regulators of drought tolerance in plants, and the enap1enap2 double mutant is more sensitive to drought stress. Both ENAP1 and ENAP2 interact with MYB44, a transcription factor that interacts with histone deacetylase HDT4. Genetics data show that myb44 null mutation enhances the sensitivity of enap1enap2 to drought stress. Whereas, HDT4 negatively regulates plant drought response, the hdt4 mutant represses enap1enap2myb44 drought sensitive phenotype. In the normal condition, ENAP1/2 and MYB44 counteract the HDT4 function for the regulation of H3K27ac. Upon drought stress, the accumulation of MYB44 and reduction of HDT4 leads to the enrichment of H3K27ac and the activation of target gene expression. Overall, this research provides a novel molecular mechanism by which ENAP1, ENAP2 and MYB44 form a complex to restrict the function of HDT4 in the normal condition; under drought condition, accumulated MYB44 and reduced HDT4 lead to the elevation of H3K27ac and the expression of drought responsive genes, as a result, plants are drought tolerant.

ECERIFERUM1-6A is required for the synthesis of cuticular wax alkanes and promotes drought tolerance in wheat

Cuticular waxes cover the aerial surfaces of land plants and protect them from various environmental stresses. Alkanes are major wax components and contribute to plant drought tolerance, but the biosynthesis and regulation of alkanes remain largely unknown in wheat (Triticum aestivum L.). Here, we identified and functionally characterized a key alkane biosynthesis gene ECERIFERUM1-6A (TaCER1-6A) from wheat. The clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9)-mediated knockout mutation in TaCER1-6A greatly reduced the contents of C27, C29, C31, and C33 alkanes in wheat leaves, while TaCER1-6A overexpression significantly increased the contents of these alkanes in wheat leaves, suggesting that TaCER1-6A is specifically involved in the biosynthesis of C27, C29, C31, and C33 alkanes on wheat leaf surfaces. TaCER1-6A knockout lines exhibited increased cuticle permeability and reduced drought tolerance, whereas TaCER1-6A overexpression lines displayed reduced cuticle permeability and enhanced drought tolerance. TaCER1-6A was highly expressed in flag leaf blades and seedling leaf blades and could respond to abiotic stresses and abscisic acid. TaCER1-6A was located in the endoplasmic reticulum, which is the subcellular compartment responsible for wax biosynthesis. A total of three haplotypes (HapI/II/III) of TaCER1-6A were identified in 43 wheat accessions, and HapI was the dominant haplotype (95%) in these wheat varieties. Additionally, we identified two R2R3-MYB transcription factors TaMYB96-2D and TaMYB96-5D that bound directly to the conserved motif CAACCA in promoters of the cuticular wax biosynthesis genes TaCER1-6A, TaCER1-1A, and fatty acyl-CoA reductase4. Collectively, these results suggest that TaCER1-6A is required for C27, C29, C31, and C33 alkanes biosynthesis and improves drought tolerance in wheat.

Genome-wide identification and characterization of the KCS gene family in sorghum (Sorghum bicolor (L.) Moench)

The aboveground parts of plants are covered with cuticle, a hydrophobic layer composed of cutin polyester and cuticular wax that can protect plants from various environmental stresses. β-Ketoacyl-CoA synthase (KCS) is the key rate-limiting enzyme in plant wax synthesis. Although the properties of KCS family genes have been investigated in many plant species, the understanding of this gene family in sorghum is still limited. Here, a total of 25 SbKCS genes were identified in the sorghum genome, which were named from SbKCS1 to SbKCS25. Evolutionary analysis among different species divided the KCS family into five subfamilies and the SbKCSs were more closely related to maize, implying a closer evolutionary relationship between sorghum and maize. All SbKCS genes were located on chromosomes 1, 2, 3, 4, 5, 6, 9 and 10, respectively, while Chr 1 and Chr 10 contained more KCS genes than other chromosomes. The prediction results of subcellular localization showed that SbKCSs were mainly expressed in the plasma membrane and mitochondria. Gene structure analysis revealed that there was 0-1 intron in the sorghum KCS family and SbKCSs within the same subgroup were similar. Multiple cis-acting elements related to abiotic stress, light and hormone response were enriched in the promoters of SbKCS genes, which indicated the functional diversity among these genes. The three-dimensional structure analysis showed that a compact spherical space structure was formed by various secondary bonds to maintain the stability of SbKCS proteins, which was necessary for their biological activity. qRT-PCR results revealed that nine randomly selected SbKCS genes expressed differently under drought and salt treatments, among which SbKCS8 showed the greatest fold of expression difference at 12 h after drought and salt stresses, which suggested that the SbKCS genes played a potential role in abiotic stress responses. Taken together, these results provided an insight into investigating the functions of KCS family in sorghum and in response to abiotic stress.

Cloning, characterization and functional analysis of NtMYB306a gene reveals its role in wax alkane biosynthesis of tobacco trichomes and stress tolerance

Trichomes are specialized hair-like organs found on epidermal cells of many terrestrial plants, which protect plant from excessive transpiration and numerous abiotic and biotic stresses. However, the genetic basis and underlying mechanisms are largely unknown in Nicotiana tabacum (common tobacco), an established model system for genetic engineering and plant breeding. In present study, we identified, cloned and characterized an unknown function transcription factor NtMYB306a from tobacco cultivar K326 trichomes. Results obtained from sequence phylogenetic tree analysis showed that NtMYB306a-encoded protein belonged to S1 subgroup of the plants' R2R3-MYB transcription factors (TFs). Observation of the green fluorescent signals from NtMYB306a-GFP fusion protein construct exhibited that NtMYB306a was localized in nucleus. In yeast transactivation assays, the transformed yeast containing pGBKT7-NtMYB306a construct was able to grow on SD/-Trp-Ade+X-α-gal selection media, signifying that NtMYB306a exhibits transcriptional activation activity. Results from qRT-PCR, in-situ hybridization and GUS staining of transgenic tobacco plants revealed that NtMYB306a is primarily expressed in tobacco trichomes, especially tall glandular trichomes (TGTs) and short glandular trichomes (SGTs). RNA sequencing (RNA-seq) and qRT-PCR analysis of the NtMYB306a-overexpressing transgenic tobacco line revealed that NtMYB306a activated the expression of a set of key target genes which were associated with wax alkane biosynthesis. Gas Chromatography-Mass Spectrometry (GC-MS) exhibited that the total alkane contents and the contents of n-C28, n-C29, n-C31, and ai-C31 alkanes in leaf exudates of NtMYB306a-OE lines (OE-3, OE-13, and OE-20) were significantly greater when compared to WT. Besides, the promoter region of NtMYB306a contained numerous stress-responsive cis-acting elements, and their differential expression towards salicylic acid and cold stress treatments reflected their roles in signal transduction and cold-stress tolerance. Together, these results suggest that NtMYB306a is necessarily a positive regulator of alkane metabolism in tobacco trichomes that does not affect the number and morphology of tobacco trichomes, and that it can be used as a candidate gene for improving stress resistance and the quality of tobacco.

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

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

Mutation of OsMYB60 reduces rice resilience to drought stress by attenuating cuticular wax biosynthesis

The cuticular wax layer on leaf surfaces limits non-stomatal water loss to the atmosphere and protects against pathogen invasion. Although many genes associated with wax biosynthesis and wax transport in plants have been identified, their regulatory mechanisms remain largely unknown. Here, we show that the MYB transcription factor OsMYB60 positively regulates cuticular wax biosynthesis and this helps rice (Oryza sativa) plants tolerate drought stress. Compared with the wild type (japonica cultivar 'Dongjin'), osmyb60 null mutants (osmyb60-1 and osmyb60-2) exhibited increased drought sensitivity, with more chlorophyll leaching and higher rates of water loss. Quantitative reverse-transcription PCR showed that the loss of function of OsMYB60 led to downregulation of wax biosynthesis genes, leading to reduced amounts of total wax components on leaf surfaces under normal conditions. Yeast one-hybrid, luciferase transient transcriptional activity, and chromatin immunoprecipitation assays revealed that OsMYB60 directly binds to the promoter of OsCER1 (a key gene involved in very-long-chain alkane biosynthesis) and upregulates its expression. Taken together, these results demonstrate that OsMYB60 enhances rice resilience to drought stress by promoting cuticular wax biosynthesis on leaf surfaces.

Genome-wide analysis of the WSD family in sunflower and functional identification of HaWSD9 involvement in wax ester biosynthesis and osmotic stress

The wax esters are important cuticular wax composition that cover the outer surface of plant organs and play a critical role in protection and energy metabolism. Wax ester synthesis in plant is catalyzed by a bifunctional wax ester synthase/acyl-CoA: diacylglycerol acyltransferase (WSD). Sunflower (Helianthus annuus L.) is an important oil crop in the world; however, little is known about WSD in sunflower. In this study, we identified and performed a functional analysis of twelve HaWSD genes from sunflower genome. Tissue-specific expression revealed that 12 HaWSD genes were differentially expressed in various organs and tissues of sunflower, except seeds. HaWSD genes were highly induced by salinity, drought, cold, and abscisic acid (ABA) in sunflower. To ascertain their function, HaWSD9, with highly expressed levels in stems and leaves, was cloned and expressed in a yeast mutant defective in triacylglycerol (TAG) biosynthesis. HaWSD9 complemented the phenotype by producing wax ester but not TAG in vivo, indicating that it functions as a wax ester synthase. Subcellular localization analysis indicated that HaWSD9 was located in the endoplasmic reticulum (ER). Heterologous introduction of HaWSD9 into Arabidopsis wsd1 mutant exhibited increased epicuticular wax crystals and cuticular wax contents on the stems. As compared with the wsd1 mutant, HaWSD9 overexpressing transgenic Arabidopsis showed less cuticle permeability, chlorophyll leaching and water loss rate. Further analysis showed that the HaWSD9 transgenics enhanced tolerance to ABA, mannitol, drought and salinity, and maintained higher leaf relative water content (RWC) than the wsd1 mutant under drought stress, suggesting that HaWSD9 play an important physiological role in stress response as well as wax synthase. These results contribute to understanding the function of HaWSD genes in wax ester synthesis and stress tolerance in sunflower.

Impact of water deficiency on leaf cuticle lipids and gene expression networks in cotton (Gossypium hirsutum L.)

BACKGROUND

Water deficit (WD) has serious effect on the productivity of crops. Formation of cuticular layer with increased content of wax and cutin on leaf surfaces is closely related to drought tolerance. Identification of drought tolerance associated wax components and cutin monomers and the genes responsible for their biosynthesis is essential for understanding the physiological and genetic mechanisms underlying drought tolerance and improving crop drought resistance.

RESULT

In this study, we conducted comparative phenotypic and transcriptomic analyses of two Gossypium hirsutum varieties that are tolerant (XL22) or sensitive (XL17) to drought stress. XL17 consumed more water than XL22, particularly under the WD conditions. WD significantly induced accumulation of most major wax components (C29 and C31 alkanes) and cutin monomers (palmitic acid and stearic acid) in leaves of both XL22 and XL17, although accumulation of the major cutin monomers, i.e., polyunsaturated linolenic acid (C18:3n-3) and linoleic acid (C18:2n-6), were significantly repressed by WD in both XL22 and XL17. According to the results of transcriptome analysis, although many genes and their related pathways were commonly induced or repressed by WD in both XL22 and XL17, WD-induced differentially expressed genes specific to XL22 or XL17 were also evident. Among the genes that were commonly induced by WD were the GhCER1 genes involved in biosynthesis of alkanes, consistent with the observation of enhanced accumulation of alkanes in cotton leaves under the WD conditions. Interestingly, under the WD conditions, several GhCYP86 genes, which encode enzymes catalyzing the omega-hydroxylation of fatty acids and were identified to be the hub genes of one of the co-expression gene modules, showed a different expression pattern between XL22 and XL17 that was in agreement with the WD-induced changes of the content of hydroxyacids or fatty alcohols in these two varieties.

CONCLUSION

The results contribute to our comprehending the physiological and genetic mechanisms underlying drought tolerance and provide possible solutions for the difference of drought resistance of different cotton varieties.

Transcriptional reprogramming during floral fate acquisition

Coordinating growth and patterning is essential for eukaryote morphogenesis. In plants, auxin is a key regulator of morphogenesis implicated throughout development. Despite this central role, our understanding of how auxin coordinates cell fate and growth changes is still limited. Here, we addressed this question using a combination of genomic screens to delve into the transcriptional network induced by auxin at the earliest stage of flower development, prior to morphological changes. We identify a shoot-specific network suggesting that auxin initiates growth through an antagonistic regulation of growth-promoting and growth-repressive hormones, quasi-synchronously to floral fate specification. We further identify two DNA-binding One Zinc Finger (DOF) transcription factors acting in an auxin-dependent network that could interface growth and cell fate from the early stages of flower development onward.

•

Pharmacological approach to probe transcriptional responses in shoot meristems

•

Analysis of a shoot-specific network regulated by auxin during flower initiation

•

Two DOF transcription factors are induced in flower primordia

•

The DOF genes potentially link growth to organ identity acquisition

AtMYB31 is a wax regulator associated with reproductive development in Arabidopsis

AtMYB31, a R2R3-MYB transcription factor that modulates wax biosynthesis in reproductive tissues, is involved in seed development in Arabidopsis. R2R3-MYB transcription factors play important roles in plant development; yet, the exact role of each of them remains to be resolved. Here we report that the Arabidopsis AtMYB31 is required for wax biosynthesis in epidermis of reproductive tissues, and is involved in seed development. AtMYB31 was ubiquitously expressed in both vegetative and reproductive tissues with higher expression levels in siliques and seeds, while AtMYB31 was localized to the nucleus and cytoplasm. Loss of function of AtMYB31 reduced wax accumulation in the epidermis of silique and flower tissues, disrupted seed coat epidermal wall development and mucilage production, altered seed proanthocyanidin and polyester content. AtMYB31 could direct activate expressions of several wax biosynthetic target genes. Altogether, AtMYB31, a R2R3-MYB transcription factor, regulates seed development in Arabidopsis.

Genome-wide identification and expression analysis of the R2R3-MYB gene family in tobacco (Nicotiana tabacum L.)

BACKGROUND

The R2R3-MYB transcription factor is one of the largest gene families in plants and involved in the regulation of plant development, hormone signal transduction, biotic and abiotic stresses. Tobacco is one of the most important model plants. Therefore, it will be of great significance to investigate the R2R3-MYB gene family and their expression patterns under abiotic stress and senescence in tobacco.

RESULTS

A total of 174 R2R3-MYB genes were identified from tobacco (Nicotiana tabacum L.) genome and were divided into 24 subgroups based on phylogenetic analysis. Gene structure (exon/intron) and protein motifs were especially conserved among the NtR2R3-MYB genes, especially members within the same subgroup. The NtR2R3-MYB genes were distributed on 24 tobacco chromosomes. Analysis of gene duplication events obtained 3 pairs of tandem duplication genes and 62 pairs of segmental duplication genes, suggesting that segmental duplications is the major pattern for R2R3-MYB gene family expansion in tobacco. Cis-regulatory elements of the NtR2R3-MYB promoters were involved in cellular development, phytohormones, environmental stress and photoresponsive. Expression profile analysis showed that NtR2R3-MYB genes were widely expressed in different maturity tobacco leaves, and however, the expression patterns of different members appeared to be diverse. The qRT-PCR analysis of 15 NtR2R3-MYBs confirmed their differential expression under different abiotic stresses (cold, salt and drought), and notably, NtMYB46 was significantly up-regulated under three treatments.

CONCLUSIONS

In summary, a genome-wide identification, evolutionary and expression analysis of R2R3-MYB gene family in tobacco were conducted. Our results provided a solid foundation for further biological functional study of NtR2R3-MYB genes in tobacco.

The R2R3 MYB Transcription Factor MYB71 Regulates Abscisic Acid Response in Arabidopsis

Abscisic acid (ABA) regulates plant responses to abiotic stresses via regulating the expression of downstream genes, yet the functions of many ABA responsive genes remain unknown. We report here the characterization of MYB71, a R2R3 MYB transcription factor in regulating ABA responses in Arabidopsis. RT-PCR results show that the expression level of MYB71 was increased in response to ABA treatment. Arabidopsis protoplasts transfection results show that MYB71 was specifically localized in nucleus and it activated the Gal4:GUS reporter gene when recruited to the Gal4 promoter by a fused DNA binding domain GD. Roles of MYB71 in regulating plant response to ABA were analyzed by generating Arabidopsis transgenic plants overexpression MYB71 and gene edited mutants of MYB71. The results show that ABA sensitivity was increased in the transgenic plants overexpression MYB71, but decreased in the MYB71 mutants. By using a DEX inducible system, we further identified genes are likely regulated by MYB71, and found that they are enriched in biological process to environmental stimuli including abiotic stresses, suggesting that MYB71 may regulate plant response to abiotic stresses. Taken together, our results suggest that MYB71 is an ABA responsive gene, and MYB71 functions as a transcription activator and it positively regulates ABA response in Arabidopsis.

Regulatory mechanisms underlying cuticular wax biosynthesis

Plants are sessile organisms that have developed hydrophobic cuticles that cover their aerial epidermal cells to protect them from terrestrial stresses. The cuticle layer is mainly composed of cutin, a polyester of hydroxy and epoxy fatty acids, and cuticular wax, a mixture of very-long-chain fatty acids (>20 carbon atoms) and their derivatives, aldehydes, alkanes, ketones, alcohols, and wax esters. During the last 30 years, forward and reverse genetic, transcriptomic, and biochemical approaches have enabled the identification of key enzymes, transporters, and regulators involved in the biosynthesis of cutin and cuticular waxes. In particular, cuticular wax biosynthesis is significantly influenced in an organ-specific manner or by environmental conditions, and is controlled using a variety of regulators. Recent studies on the regulatory mechanisms underlying cuticular wax biosynthesis have enabled us to understand how plants finely control carbon metabolic pathways to balance between optimal growth and development and defense against abiotic and biotic stresses. In this review, we summarize the regulatory mechanisms underlying cuticular wax biosynthesis at the transcriptional, post-transcriptional, post-translational, and epigenetic levels.

Synthesis of C20-38 Fatty Acids in Plant Tissues

Very-long-chain fatty acids (VLCFA) are involved in a number of important plant physiological functions. Disorders in the expression of genes involved in the synthesis of VLCFA lead to a number of phenotypic consequences, ranging from growth retardation to the death of embryos. The elongation of VLCFA in the endoplasmic reticulum (ER) is carried out by multiple elongase complexes with different substrate specificities and adapted to the synthesis of a number of products required for a number of metabolic pathways. The information about the enzymes involved in the synthesis of VLCFA with more than 26 atoms of Carbon is rather poor. Recently, genes encoding enzymes involved in the synthesis of both regular-length fatty acids and VLCFA have been discovered and investigated. Polyunsaturated VLCFA in plants are formed mainly by 20:1 elongation into new monounsaturated acids, which are then imported into chloroplasts, where they are further desaturated. The formation of saturated VLCFA and their further transformation into a number of aliphatic compounds included in cuticular waxes and suberin require the coordinated activity of a large number of different enzymes.

ArabidopsisKCS5 and KCS6 Play Redundant Roles in Wax Synthesis

3-ketoacyl-CoA synthases (KCSs), as components of a fatty acid elongase (FAE) complex, play key roles in determining the chain length of very-long-chain fatty acids (VLCFAs). KCS6, taking a predominate role during the elongation from C26 to C28, is well known to play an important role in wax synthesis. KCS5 is one paralog of KCS6 and its role in wax synthesis remains unknown. Wax phenotype analysis showed that in kcs5 mutants, the total amounts of wax components derived from carbon 32 (C32) and C34 were apparently decreased in leaves, and those of C26 to C32 derivatives were obviously decreased in flowers. Heterologous yeast expression analysis showed that KCS5 alone displayed specificity towards C24 to C28 acids, and its coordination with CER2 and CER26 catalyzed the elongation of acids exceeding C28, especially displaying higher activity towards C28 acids than KCS6. BiLC experiments identified that KCS5 physically interacts with CER2 and CER26. Wax phenotype analysis of different organs in kcs5 and kcs6 single or double mutants showed that KCS6 mutation causes greater effects on the wax synthesis than KCS5 mutation in the tested organs, and simultaneous repression of both protein activities caused additive effects, suggesting that during the wax biosynthesis process, KCS5 and KCS6 play redundant roles, among which KCS6 plays a major role. In addition, simultaneous mutations of two genes nearly block drought-induced wax production, indicating that the reactions catalyzed by KCS5 and KCS6 play a critical role in the wax biosynthesis in response to drought.

Exploring Camelina sativa lipid metabolism regulation by combining gene co-expression and DNA affinity purification analyses

Camelina (Camelina sativa) is an annual oilseed plant that is gaining momentum as a biofuel cover crop. Understanding gene regulatory networks is essential to deciphering plant metabolic pathways, including lipid metabolism. Here, we take advantage of a growing collection of gene expression datasets to predict transcription factors (TFs) associated with the control of Camelina lipid metabolism. We identified approximately 350 TFs highly co-expressed with lipid-related genes (LRGs). These TFs are highly represented in the MYB, AP2/ERF, bZIP, and bHLH families, including a significant number of homologs of well-known Arabidopsis lipid and seed developmental regulators. After prioritizing the top 22 TFs for further validation, we identified DNA-binding sites and predicted target genes for 16 out of the 22 TFs tested using DNA affinity purification followed by sequencing (DAP-seq). Enrichment analyses of targets supported the co-expression prediction for most TF candidates, and the comparison to Arabidopsis revealed some common themes, but also aspects unique to Camelina. Within the top potential lipid regulators, we identified CsaMYB1, CsaABI3AVP1-2, CsaHB1, CsaNAC2, CsaMYB3, and CsaNAC1 as likely involved in the control of seed fatty acid elongation and CsaABI3AVP1-2 and CsabZIP1 as potential regulators of the synthesis and degradation of triacylglycerols (TAGs), respectively. Altogether, the integration of co-expression data and DNA-binding assays permitted us to generate a high-confidence and short list of Camelina TFs involved in the control of lipid metabolism during seed development.

Increased cuticular wax deposition does not change residual foliar transpiration

The effect of contrasting environmental growth conditions (in vitro tissue culture, ex vitro acclimatisation, climate chamber, greenhouse and outdoor) on leaf development, cuticular wax composition, and foliar transpiration of detached leaves of the Populus × canescens clone 84 K were investigated. Our results show that total amounts of cuticular wax increased more than 10-fold when cultivated in different growth conditions, whereas qualitative wax composition did not change. With exception of plants directly taken from tissue culture showing rapid dehydration, rates of water loss (residual foliar transpiration) of intact but detached leaves were constant and independent from growth conditions and thus independent from increasing wax amounts. Since cuticular transpiration measured with isolated astomatous P. × canescens cuticles was identical to residual foliar transpiration rates of detached leaves, our results confirm that cuticular transpiration of P. × canescens leaves can be predicted with high accuracy from residual transpiration of detached leaves after stomatal closure. Our results convincingly show that more than 10-fold increased wax amounts in P. × canescens cuticles do not lead to decreased rates of residual (cuticular) transpiration.

Molecular insights into AabZIP1-mediated regulation on artemisinin biosynthesis and drought tolerance in Artemisia annua

Artemisia annua is the main natural source of artemisinin production. In A. annua, extended drought stress severely reduces its biomass and artemisinin production while short-term water-withholding or abscisic acid (ABA) treatment can increase artemisinin biosynthesis. ABA-responsive transcription factor AabZIP1 and JA signaling AaMYC2 have been shown in separate studies to promote artemisinin production by targeting several artemisinin biosynthesis genes. Here, we found AabZIP1 promote the expression of multiple artemisinin biosynthesis genes including AaDBR2 and AaALDH1, which AabZIP1 does not directly activate. Subsequently, it was found that AabZIP1 up-regulates AaMYC2 expression through direct binding to its promoter, and that AaMYC2 binds to the promoter of AaALDH1 to activate its transcription. In addition, AabZIP1 directly transactivates wax biosynthesis genes AaCER1 and AaCYP86A1. The biosynthesis of artemisinin and cuticular wax and the tolerance of drought stress were significantly increased by AabZIP1 overexpression, whereas they were significantly decreased in RNAi-AabZIP1 plants. Collectively, we have uncovered the AabZIP1-AaMYC2 transcriptional module as a point of cross-talk between ABA and JA signaling in artemisinin biosynthesis, which may have general implications. We have also identified AabZIP1 as a promising candidate gene for the development of A. annua plants with high artemisinin content and drought tolerance in metabolic engineering breeding.

Tomato SlCER1-1 catalyzes the synthesis of wax alkanes which increases the drought tolerance and fruit storability

Very-long-chain (VLC) alkanes are the main wax compounds of tomato fruit and leaf. ECERIFERUM1 (CER1) and ECERIFERUM3 (CER3) are the two key genes involved in VLC alkane biosynthesis in Arabidopsis thaliana. However, the CER1 and CER3 homologous genes in tomato have not been investigated and their exact biological function remains unknown. We analyzed the wax profiles in tomato leaves and fruits at different growth stages, and characterized the CER1 and CER3 homologous genes. VLC alkanes were the predominant wax compounds both in the leaf and fruit at all developmental stages. We identified five CER1 homologs and two CER3 homologs in tomato, which were designated as SlCER1-1 to SlCER1-5 and SlCER3-1 and SlCER3-2 respectively. The genes exhibited tissue- and organ-dependent expression patterns and were induced by abiotic stresses. SlCER1-1 was localized to the endoplasmic reticulum (ER), which is also the main site of wax biosynthesis. Silencing the SlCER1-1 gene in tomato significantly reduced the amounts of n-Alkanes and branched alkanes, whereas its overexpression in Arabidopsis had the opposite effect. Under drought stress, both n-Alkanes and branched alkanes increased significantly in wild-type but not the SlCER1-1 RNAi tomato plants. Furthermore, SlCER1-1 silencing also increased the cuticular permeabilities of the leaves and fruits. In conclusion, SlCER1-1 is involved in wax alkane biosynthesis in tomato and plays an important role in the drought tolerance and fruit storability.

Function and transcriptional regulation of CsKCS20 in the elongation of very-long-chain fatty acids and wax biosynthesis in Citrus sinensis flavedo

Cuticular wax on plant aerial surfaces plays a vital role in the defense against various stresses, and the genes related to wax metabolism have been well documented in several model plants. However, there is very limited research on the key enzymes and transcription factors (TFs) associated with carbon chain distribution and wax biosynthesis in citrus fruit. In this study, an analysis of wax metabolites indicated that even carbon-chain (C24-C28) metabolites are the dominant wax components in citrus fruit, and a 3-ketoacyl-CoA synthase (KCS) family gene (CsKCS20) plays an important role in the carbon chain distribution during wax biosynthesis in a wax-deficient mutant (MT). Expression of CsKCS20 in yeast indicated that CsKCS20 can catalyze the biosynthesis of C22 and C24 very-long-chain fatty acids (VLCFAs). In addition, transcriptome and sequence analysis indicated that the differential expression of CsKCS20 between the wild-type (WT) and MT fruit can be partly attributed to the regulation of CsMYB96, which was further confirmed by yeast one-hybrid (Y1H) assays, electrophoretic mobility shift assays (EMSAs) and dual luciferase assays. The functions of CsMYB96 and CsKCS20 in wax biosynthesis were further validated by heterologous expression in Arabidopsis. In summary, this study elucidates the important roles of CsKCS20 and CsMYB96 in regulating VLCFA elongation and cuticular wax biosynthesis, which provides new directions for the improvement of citrus fruit wax quality in genetic breeding programs.

Ectopic Overexpression of CsECR From Navel Orange Increases Cuticular Wax Accumulation in Tomato and Enhances Its Tolerance to Drought Stress

Drought stress often occurred in citrus to limit its growth, distribution, and fruit quality. Cuticular waxes play an important role in regulating plant tolerance to drought stress. Plant enoyl-CoA reductase (ECR) is involved in the biosynthesis of cuticular waxes and catalyzes the last step of very long-chain fatty acids (VLCFAs) elongation. In this study, a putative ECR gene, named CsECR, was cloned from "Newhall" navel orange. CsECR protein has high identities with other plant ECR proteins and contained a conserved NADP/NAD-binding motif and three conserved functional sites. The highest expression of CsECR was observed in leaves, followed by stems, flavedos, ovaries, juice sacs, stigmas, stamens, albedos, and petals. Besides, the expression of CsECR was significantly induced by PEG6000 and ABA treatments. Ectopic overexpression of CsECR increased the contents of total waxes and aliphatic wax fractions (n-fatty acids, unsaturated fatty acids, n-alkanes, alkenes, iso-, and anteiso-alkanes) in the leaves and fruits of the transgenic tomato. Furthermore, ectopic overexpression of CsECR reduced the cuticle permeability in the leaves and fruits of the transgenic tomato and increased its tolerance to drought stress. Taken together, our results revealed that CsECR plays an important role in plant response to drought stresses by regulating cuticular wax biosynthesis.

Molecular and Physiological Perspectives of Abscisic Acid Mediated Drought Adjustment Strategies

Drought is the most prevalent unfavorable condition that impairs plant growth and development by altering morphological, physiological, and biochemical functions, thereby impeding plant biomass production. To survive the adverse effects, water limiting condition triggers a sophisticated adjustment mechanism orchestrated mainly by hormones that directly protect plants via the stimulation of several signaling cascades. Predominantly, water deficit signals cause the increase in the level of endogenous ABA, which elicits signaling pathways involving transcription factors that enhance resistance mechanisms to combat drought-stimulated damage in plants. These responses mainly include stomatal closure, seed dormancy, cuticular wax deposition, leaf senescence, and alteration of the shoot and root growth. Unraveling how plants adjust to drought could provide valuable information, and a comprehensive understanding of the resistance mechanisms will help researchers design ways to improve crop performance under water limiting conditions. This review deals with the past and recent updates of ABA-mediated molecular mechanisms that plants can implement to cope with the challenges of drought stress.

The ARRE RING-Type E3 Ubiquitin Ligase Negatively Regulates Cuticular Wax Biosynthesis in Arabidopsis thaliana by Controlling ECERIFERUM1 and ECERIFERUM3 Protein Levels

The outer epidermal cell walls of plant shoots are covered with a cuticle, a continuous lipid structure that provides protection from desiccation, UV light, pathogens, and insects. The cuticle is mostly composed of cutin and cuticular wax. Cuticular wax synthesis is synchronized with surface area expansion during plant development and is associated with plant responses to biotic and abiotic stresses. Cuticular wax deposition is tightly regulated by well-established transcriptional and post-transcriptional regulatory mechanisms, as well as post-translationally via the ubiquitin-26S proteasome system (UPS). The UPS is highly conserved in eukaryotes and involves the covalent attachment of polyubiquitin chains to the target protein by an E3 ligase, followed by the degradation of the modified protein by the 26S proteasome. A large number of E3 ligases are encoded in the Arabidopsis genome, but only a few have been implicated in the regulation of cuticular wax deposition. In this study, we have conducted an E3 ligase reverse genetic screen and identified a novel RING-type E3 ubiquitin ligase, AtARRE, which negatively regulates wax biosynthesis in Arabidopsis. Arabidopsis plants overexpressing AtARRE exhibit glossy stems and siliques, reduced fertility and fusion between aerial organs. Wax load and wax compositional analyses of AtARRE overexpressors showed that the alkane-forming branch of the wax biosynthetic pathway is affected. Co-expression of AtARRE and candidate target proteins involved in alkane formation in both Nicotiana benthamiana and stable Arabidopsis transgenic lines demonstrated that AtARRE controls the levels of wax biosynthetic enzymes ECERIFERUM1 (CER1) and ECERIFERUM3 (CER3). CER1 has also been confirmed to be a ubiquitination substrate of the AtARRE E3 ligase by an in vivo ubiquitination assay using a reconstituted Escherichia coli system. The AtARRE gene is expressed throughout the plant, with the highest expression detected in fully expanded rosette leaves and oldest stem internodes. AtARRE gene expression can also be induced by exposure to pathogens. These findings reveal that wax biosynthesis in mature plant tissues and in response to pathogen infection is controlled post-translationally.

Peridermal fruit skin formation in Actinidia sp. (kiwifruit) is associated with genetic loci controlling russeting and cuticle formation

BACKGROUND

The skin (exocarp) of fleshy fruit is hugely diverse across species. Most fruit types have a live epidermal skin covered by a layer of cuticle made up of cutin while a few create an outermost layer of dead cells (peridermal layer).

RESULTS

In this study we undertook crosses between epidermal and peridermal skinned kiwifruit, and showed that epidermal skin is a semi-dominant trait. Furthermore, backcrossing these epidermal skinned hybrids to a peridermal skinned fruit created a diverse range of phenotypes ranging from epidermal skinned fruit, through fruit with varying degrees of patches of periderm (russeting), to fruit with a complete periderm. Quantitative trait locus (QTL) analysis of this population suggested that periderm formation was associated with four loci. These QTLs were aligned either to ones associated with russet formation on chromosome 19 and 24, or cuticle integrity and coverage located on chromosomes 3, 11 and 24.

CONCLUSION

From the segregation of skin type and QTL analysis, it appears that skin development in kiwifruit is controlled by two competing factors, cuticle strength and propensity to russet. A strong cuticle will inhibit russeting while a strong propensity to russet can create a continuous dead skinned periderm.

Evolution and Characterization of Acetyl Coenzyme A: Diacylglycerol Acyltransferase Genes in Cotton Identify the Roles of GhDGAT3D in Oil Biosynthesis and Fatty Acid Composition

Cottonseed oil is rich in unsaturated fatty acids (UFAs) and serves as an edible oil in human nutrition. Reports suggest that acyl-coenzyme A: diacylglycerol acyltransferases (DGAT) and wax ester synthase/DGAT (WSD1) genes encode a key group of enzymes that catalyze the final step for triacylglycerol biosynthesis and enable an important rate-limiting process. However, their roles in oil biosynthesis and the fatty acid profile of cotton seed are poorly understood. Therefore, the aim of this study was to identify and characterize DGAT and WSD1 genes in cotton plants and examine their roles in oil biosynthesis, the fatty acid profile of cotton seeds, and abiotic stress responses. In this study, 36 GhDGAT and GhWSD1 genes were identified in upland cotton (G. hirsutum) and found to be clustered into four groups: GhDGAT1, GhDGAT2, GhDGAT3, and GhWSD1. Gene structure and domain analyses showed that the GhDGAT and GhWSD1 genes in each group are highly conserved. Gene synteny analysis indicated that segmental and tandem duplication events occurred frequently during cotton evolution. Expression analysis revealed that GhDGAT and GhWSD1 genes function widely in cotton development and stress responses; moreover, several environmental stress and hormone response-related cis-elements were detected in the GhDGAT and GhWSD1 promoter regions. The predicted target transcription factors and miRNAs imply an extensive role of GhDGAT and GhWSD1 genes in stress responses. Increases in GhDGAT3 gene expression with increases in cottonseed oil accumulation were observed. Transformation study results showed that there was an increase in C18:1 content and a decrease in C18:2 and C18:3 contents in seeds of Arabidopsis transgenic plants overexpressing GhDGAT3D compared with that of control plants. Overall, these findings contributed to the understanding of the functions of GhDGAT and GhWSD1 genes in upland cotton, providing basic information for further research.

Multiple Functions of MYB Transcription Factors in Abiotic Stress Responses

Plants face a more volatile environment than other organisms because of their immobility, and they have developed highly efficient mechanisms to adapt to stress conditions. Transcription factors, as an important part of the adaptation process, are activated by different signals and are responsible for the expression of stress-responsive genes. MYB transcription factors, as one of the most widespread transcription factor families in plants, participate in plant development and responses to stresses by combining with MYB cis-elements in promoters of target genes. MYB transcription factors have been extensively studied and have proven to be critical in the biosynthesis of secondary metabolites in plants, including anthocyanins, flavonols, and lignin. Multiple studies have now shown that MYB proteins play diverse roles in the responses to abiotic stresses, such as drought, salt, and cold stresses. However, the regulatory mechanism of MYB proteins in abiotic stresses is still not well understood. In this review, we will focus mainly on the function of Arabidopsis MYB transcription factors in abiotic stresses, especially how MYB proteins participate in these stress responses. We also pay attention to how the MYB proteins are regulated in these processes at both the transcript and protein levels.