A systems-oriented analysis of the grapevine R2R3-MYB transcription factor family uncovers new insights into the regulation of stilbene accumulation

NAC61 regulates late- and post-ripening osmotic, oxidative, and biotic stress responses in grapevine

During late- and post-ripening stages, grape berry undergoes profound biochemical and physiological changes whose molecular control is poorly understood. Here, we report the role of NAC61, a grapevine NAC transcription factor, in regulating different processes involved in berry ripening progression. NAC61 is highly expressed during post-harvest berry dehydration and its expression pattern is closely related to sugar concentration. The ectopic expression of NAC61 in Nicotiana benthamiana leaves resulted in low stomatal conductance, high leaf temperature, tissue collapse and a higher relative water content. Transcriptome analysis of grapevine leaves transiently overexpressing NAC61 and DNA affinity purification and sequencing analyses allowed us to narrow down a list of NAC61-regulated genes. Direct regulation of the stilbene synthase regulator MYB14, the osmotic stress-related gene DHN1b, the Botrytis cinerea susceptibility gene WRKY52, and NAC61 itself was validated. We also demonstrate that NAC61 interacts with NAC60, a proposed master regulator of grapevine organ maturation, in the activation of MYB14 and NAC61 expression. Overall, our findings establish NAC61 as a key player in a regulatory network that governs stilbenoid metabolism and osmotic, oxidative, and biotic stress responses that are the hallmark of late- and post-ripening grape stages.

Exploring large-scale gene coexpression networks in peach (Prunus persica L.): a new tool for predicting gene function

Peach is a model for Prunus genetics and genomics, however, identifying and validating genes associated to peach breeding traits is a complex task. A gene coexpression network (GCN) capable of capturing stable gene-gene relationships would help researchers overcome the intrinsic limitations of peach genetics and genomics approaches and outline future research opportunities. In this study, we created four GCNs from 604 Illumina RNA-Seq libraries. We evaluated the performance of every GCN in predicting functional annotations using an algorithm based on the 'guilty-by-association' principle. The GCN with the best performance was COO300, encompassing 21 956 genes. To validate its performance predicting gene function, we performed two case studies. In case study 1, we used two genes involved in fruit flesh softening: the endopolygalacturonases PpPG21 and PpPG22. Genes coexpressing with both genes were extracted and referred to as melting flesh (MF) network. Finally, we performed an enrichment analysis of MF network and compared the results with the current knowledge regarding peach fruit softening. The MF network mostly included genes involved in cell wall expansion and remodeling, and with expressions triggered by ripening-related phytohormones, such as ethylene, auxin, and methyl jasmonate. In case study 2, we explored potential targets of the anthocyanin regulator PpMYB10.1 by comparing its gene-centered coexpression network with that of its grapevine orthologues, identifying a common regulatory network. These results validated COO300 as a powerful tool for peach and Prunus research. This network, renamed as PeachGCN v1.0, and the scripts required to perform a function prediction analysis are available at https://github.com/felipecobos/PeachGCN.

Overexpression of a 'Beta' MYB Factor Gene, VhMYB15, Increases Salinity and Drought Tolerance in Arabidopsis thaliana

'Beta' is a hybrid of Vitis riparia L. and V. labrusca and has a strong ability to adapt to adverse growth environments and is mainly cultivated and used as a resistant rootstock. At present, the most extensively studied MYB TFs are R2R3-type, which have been found to be involved in plant growth, development, and stress response processes. In the present research, VhMYB15, a key transcription factor for abiotic stress tolerance, was screened by bioinformatics in 'Beta' rootstock, and its function under salinity and drought stresses was investigated. VhMYB15 was highly expressed in roots and mature leave under salinity and drought stresses. Observing the phenotype and calculating the survival rate of plants, it was found that VhMYB15-overexpressing plants exhibited relatively less yellowing and wilting of leaves and a higher survival rate under salinity and drought stresses. Consistent with the above results, through the determination of stress-related physiological indicators and the expression analysis of stress-related genes (AtSOS2, AtSOS3, AtSOS1, AtNHX1, AtSnRK2.6, AtNCED3, AtP5CS1, and AtCAT1), it was found that transgenic Arabidopsis showed better stress tolerance and stronger adaptability under salinity and drought stresses. Based on the above data, it was preliminarily indicated that VhMYB15 may be a key factor in salinity and drought regulation networks, enhancing the adaptability of 'Beta' to adverse environments.

MYB24 orchestrates terpene and flavonol metabolism as light responses to anthocyanin depletion in variegated grape berries

Variegation is a rare type of mosaicism not fully studied in plants, especially fruits. We examined red and white sections of grape (Vitis vinifera cv. 'Béquignol') variegated berries and found that accumulation of products from branches of the phenylpropanoid and isoprenoid pathways showed an opposite tendency. Light-responsive flavonol and monoterpene levels increased in anthocyanin-depleted areas in correlation with increasing MYB24 expression. Cistrome analysis suggested that MYB24 binds to the promoters of 22 terpene synthase (TPS) genes, as well as 32 photosynthesis/light-related genes, including carotenoid pathway members, the flavonol regulator HY5 HOMOLOGUE (HYH), and other radiation response genes. Indeed, TPS35, TPS09, the carotenoid isomerase gene CRTISO2, and HYH were activated in the presence of MYB24 and MYC2. We suggest that MYB24 modulates ultraviolet and high-intensity visible light stress responses that include terpene and flavonol synthesis and potentially affects carotenoids. The MYB24 regulatory network is developmentally triggered after the onset of berry ripening, while the absence of anthocyanin sunscreens accelerates its activation, likely in a dose-dependent manner due to increased radiation exposure. Anthocyanins and flavonols in variegated berry skins act as effective sunscreens but for different wavelength ranges. The expression patterns of stress marker genes in red and white sections of 'Béquignol' berries strongly suggest that MYB24 promotes light stress amelioration but only partly succeeds during late ripening.

The flavour of grape colour: anthocyanin content tunes aroma precursor composition by altering the berry microenvironment

Anthocyaninless (white) instead of black/red (coloured) fruits develop in grapevine cultivars without functional VviMYBA1 and VviMYBA2 genes, and this conditions the colour of wines that can be produced. To evaluate whether this genetic variation has additional consequences on fruit ripening and composition, we performed comparisons of microenvironment, transcriptomics, and metabolomics of developing grapes between near-isogenic white- and black-berried somatic variants of Garnacha and Tempranillo cultivars. Berry temperature was as much as 3.5 ºC lower in white- compared to black-berried Tempranillo. An RNA-seq study combined with targeted and untargeted metabolomics revealed that ripening fruits of white-berried variants were characterized by the up-regulation of photosynthesis-related and other light-responsive genes and by their higher accumulation of specific terpene aroma precursors, fatty acid-derived aldehyde volatiles, and phenylpropanoid precursor amino acids. MYBA1-MYBA2 function proved essential for flavonol trihydroxylation in black-berried somatic variants, which were also characterized by enhanced expression of pathogen defence genes in the berry skin and increased accumulation of C6-derived alcohol and ester volatiles and γ-aminobutyric acid. Collectively, our results indicate that anthocyanin depletion has side-effects on grape composition by altering the internal microenvironment of the berry and the partitioning of the phenylpropanoid pathway. Our findings show how fruit colour can condition other fruit features, such as flavour potential and stress homeostasis.

Transcriptomic Analysis Reveals Novel Regulators of the Scots Pine Stilbene Pathway

Stilbenes accumulate in Scots pine heartwood where they have important roles in protecting wood from decaying fungi. They are also part of active defense responses, and their production is induced by different (a)biotic stressors. The specific transcriptional regulators as well as the enzyme responsible for activating the stilbene precursor cinnamate in the pathway are still unknown. UV-C radiation was the first discovered artificial stress activator of the pathway. Here, we describe a large-scale transcriptomic analysis of pine needles in response to UV-C and treatment with translational inhibitors, both activating the transcription of stilbene pathway genes. We used the data to identify putative candidates for the missing CoA ligase and for pathway regulators. We further showed that the pathway is transcriptionally activated by phosphatase inhibitor, ethylene and jasmonate treatments, as in grapevine, and that the stilbene synthase promoter retains its inducibility in some of the tested conditions in Arabidopsis, a species that normally does not synthesize stilbenes. Shared features between gymnosperm and angiosperm regulation and partially retained inducibility in Arabidopsis suggest that pathway regulation occurs not only via ancient stress-response pathway(s) but also via species-specific regulators. Understanding which genes control the biosynthesis of stilbenes in Scots pine aids breeding of more resistant trees.

Genomic-wide identification and expression analysis of R2R3-MYB transcription factors related to flavonol biosynthesis in Morinda officinalis

BACKGROUND

The R2R3-MYB transcription factors are a crucial and extensive gene family in plants, which participate in diverse processes, including development, metabolism, defense, differentiation, and stress response. In the Lingnan region of China, Morinda officinalis is extensively grown and is renowned for its use as both a medicinal herb and food source. However, there are relatively few reports on the R2R3-MYB transcription factor family in M.officinalis.

RESULTS

In this study, we identified 97 R2R3-MYB genes in the genome of Morinda officinalis and classified them into 32 subgroups based on phylogenetic comparison with Arabidopsis thaliana. The lack of recent whole-genome duplication events in M.officinalis may be the reason for the relatively few members of the R2R3-MYB family. We also further analyzed the physical and chemical characteristics, conserved motifs, gene structure, and chromosomal location. Gene duplication events found 21 fragment duplication pairs and five tandem duplication event R2R3-MYB genes in M.officinalis may also affect gene family expansion. Based on phylogenetic analysis, cis-element analysis, co-expression analysis and RT-qPCR, we concluded that MoMYB33 might modulate flavonol levels by regulating the expression of 4-coumarate-CoA ligase Mo4CL2, chalcone isomerase MoCHI3, and flavonol synthase MoFLS4/11/12. MoMYB33 and AtMYB111 showed the highest similarity of 79% and may be involved in flavonol synthase networks by the STRING database. Moreover, we also identified MoMYB genes that respond to methyl Jasmonate (MeJA) and abscisic acid (ABA) stress by RT-qPCR.

CONCLUSIONS

This study offers a thorough comprehension of R2R3-MYB in M.officinalis, which lays the foundation for the regulation of flavonol synthesis and the response of MoMYB genes to phytohormones in M.officinalis.

Function of the R2R3-MYB Transcription Factors in Dalbergia odorifera and Their Relationship with Heartwood Formation

R2R3-MYB transcription factors (TFs) form one of the most important TF families involved in regulating various physiological functions in plants. The heartwood of Dalbergia odorifera is a kind of high-grade mahogany and valuable herbal medicine with wide application. However, the role of R2R3-MYB genes in the growth and development of D. odorifera, especially their relevance to heartwood formation, has not been revealed. A total of 126 R2R3-MYBs were screened from the D. odorifera genome and named DodMYB1-126 based on their location on 10 chromosomes. The collinearity results showed that purification selection was the main driving force for the evolution of the R2R3-MYB TFs family, and whole genome/fragment replication event was the main form for expanding the R2R3-MYB family, generating a divergence of gene structure and function. Comparative phylogenetic analysis classified the R2R3-MYB TFs into 33 subfamilies. S3-7,10,12-13,21 and N4-7 were extensively involved in the metabolic process; S9,13,16-19,24-25 and N1-3,8 were associated with the growth and development of D. odorifera. Based on the differential transcriptional expression levels of R2R3-MYBs in different tissues, DodMYB32, DodMYB55, and DodMYB89 were tentatively screened for involvement in the regulatory process of heartwood. Further studies have shown that the DodMYB89, localized in the nucleus, has transcriptional activation activity and is involved in regulating the biosynthesis of the secondary metabolites of heartwood by activating the promoters of the structural genes DodI2'H and DodCOMT. This study aimed to comprehensively analyze the functions of the R2R3-MYB TFs and screen for candidate genes that might be involved in heartwood formation of D. odorifera.

Overexpression of a Grape MYB Transcription Factor Gene VhMYB2 Increases Salinity and Drought Tolerance in Arabidopsis thaliana

In viticulture, the highly resistant rootstock 'Beta' is widely used in Chinese grape production to avoid the effects of soil salinization and drought on grape growth. However, the mechanism of high resistance to abiotic stress in the 'Beta' rootstock is not clear. In this study, we demonstrated that VhMYB2 as a transcription factor made a significant contribution to salinity and drought stress, which was isolated from the 'Beta' rootstock. The coding sequence of the VhMYB2 gene was 858 bp, encoding 285 amino acids. The subcellular localization of VhMYB2 was located in the nucleus of tobacco epidermal cells. Moreover, RT-qPCR found that VhMYB2 was predominantly expressed in the mature leaf and root of the grape. Under salinity and drought stress, overexpressing VhMYB2 showed a higher resistant phenotype and survival rates in A. thaliana while the transgenic lines had a survival advantage by measuring the contents of proline, chlorophyll, and MDA, and activities of POD, SOD, and CAT, and expression levels of related stress response genes. The results reveal that VhMYB2 may be an important transcription factor regulating 'Beta' resistance in response to abiotic stress.

Identification of Vitis vinifera MYB transcription factors and their response against grapevine berry inner necrosis virus

BACKGROUND

The myeloblastosis (MYB) superfamily is the largest transcription factor family in plants that play diverse roles during stress responses. However, the biotic stress-responsive MYB transcription factors of the grapevine have not been systematically studied. In China, grapevine berries are often infected with the grapevine berry inner necrosis virus (GINV), which eventually reduces the nutritional quality and commodity value.

RESULTS

The present study identified and characterized 265 VvMYB or VvMYB-related genes of the "Crimson seedless" grapevine. Based on DNA-binding domain analysis, these VvMYB proteins were classified into four subfamilies, including MYB-related, 2R-MYB, 3R-MYB, and 4R-MYB. Phylogenetic analysis divided the MYB transcription factors into 26 subgroups. Overexpression of VvMYB58 suppressed GINV abundance in the grapevine. Further qPCR indicated that among 41 randomly selected VvMYB genes, 12 were induced during GINV infection, while 28 were downregulated. These findings suggest that VvMYB genes actively regulate defense response in the grapevine.

CONCLUSION

A deeper understanding of the MYB TFs engaged in GINV defense response will help devise better management strategies. The present study also provides a foundation for further research on the functions of the MYB transcription factors.

Dissecting the effect of soil on plant phenology and berry transcriptional plasticity in two Italian grapevine varieties (Vitis vinifera L.)

Grapevine embodies a fascinating species as regards phenotypic plasticity and genotype-per-environment interactions. The terroir, namely the set of agri-environmental factors to which a variety is subjected, can influence the phenotype at the physiological, molecular, and biochemical level, representing an important phenomenon connected to the typicality of productions. We investigated the determinants of plasticity by conducting a field-experiment where all terroir variables, except soil, were kept as constant as possible. We isolated the effect of soils collected from different areas, on phenology, physiology, and transcriptional responses of skin and flesh of a red and a white variety of great economic value: Corvina and Glera. Molecular results, together with physio-phenological parameters, suggest a specific effect of soil on grapevine plastic response, highlighting a higher transcriptional plasticity of Glera in respect to Corvina and a marked response of skin compared to flesh. Using a novel statistical approach, we identified clusters of plastic genes subjected to the specific influence of soil. These findings could represent an issue of applicative value, posing the basis for targeted agricultural practices to enhance the desired characteristics for any soil/cultivar combination, to improve vineyards management for a better resource usage and to valorize vineyards uniqueness maximizing the terroir-effect.

Identification and Phylogenetic Analysis of the R2R3-MYB Subfamily in Brassica napus

The R2R3-MYB sub-family proteins are composed of most members of MYB (v-Myb avian myeloblastosis viral oncogene homolog) protein, a plant-specific transcription factor (TF) that is classified into four classes depending on the number of MYB repeats. R2R3-MYB TFs are involved in physiological and biochemical processes. However, the functions of the Brassica napus R2R3-MYB genes are still mainly unknown. In this study, 35 Brassica napus MYB (BnaMYB) genes were screened in the genome of Brassica napus, and details about their physical and chemical characteristics, evolutionary relationships, chromosome locations, gene structures, three-dimensional protein structures, cis-acting promoter elements, and gene duplications were uncovered. The BnaMYB genes have undergone segmental duplications and positive selection pressure, according to evolutionary studies. The same subfamilies have similar intron-exon patterns and motifs, according to the genes' structure and conserved motifs. Additionally, through cis-element analysis, many drought-responsive and other stress-responsive cis-elements have been found in the promoter regions of the BnaMYB genes. The expression of the BnaMYB gene displays a variety of tissue-specific patterns. Ten lignin-related genes were chosen for drought treatment. Our research screened four genes that showed significant upregulation under drought stress, and thus may be important drought-responsive genes. The findings lay a new foundation for understanding the complex mechanisms of BnaMYB in multiple developmental stages and pathways related to drought stress in rapeseed.

MYB30 and MYB14 form a repressor-activator module with WRKY8 that controls stilbene biosynthesis in grapevine

When exposed to pathogen infection or ultraviolet (UV) radiation, grapevine (Vitis vinifera) plants rapidly accumulate the stilbenoid resveratrol (Res) with concomitant increase of stilbene synthase (STS), the key enzyme in stilbene biosynthesis. Although a few transcription factors have been shown to regulate STSs, the molecular mechanism governing the regulation of STSs is not well elucidated. Our previous work showed that a VvMYB14-VvWRKY8 regulatory loop fine-tunes stilbene biosynthesis in grapevine through protein-protein interaction; overexpression of VvWRKY8 down-regulates VvMYB14 and VvSTS15/21; and application of exogenous Res up-regulates WRKY8 expression. Here, we identified an R2R3-MYB repressor, VvMYB30, which competes with the activator VvMYB14 for binding to the common binding sites in the VvSTS15/21 promoter. Similar to VvMYB14, VvMYB30 physically interacts with VvWRKY8 through their N-termini, forming a complex that does not bind DNA. Exposure to UV-B/C stress induces VvMYB14, VvWRKY8, and VvSTS15/21, but represses VvMYB30 in grapevine leaves. In addition, MYB30 expression is up-regulated by VvWRKY8-overexpression or exogenous Res. These findings suggest that the VvMYB14-VvWRKY8-VvMYB30 regulatory circuit allows grapevine to respond to UV stress by producing Res and prevents over-accumulation of Res to balance metabolic costs. Our work highlights the stress-mediated induction and feedback inhibition of stilbene biosynthesis through a complex regulatory network involving multiple positive and negative transcriptional regulators.

Comparative Genomics and Functional Studies of Putative m6A Methyltransferase (METTL) Genes in Cotton

N6-methyladenosine (m6A) RNA modification plays important regulatory roles in plant development and adapting to the environment, which requires methyltransferases to achieve the methylation process. However, there has been no research regarding m6A RNA methyltransferases in cotton. Here, a systematic analysis of the m6A methyltransferase (METTL) gene family was performed on twelve cotton species, resulting in six METTLs identified in five allotetraploid cottons, respectively, and three to four METTLs in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that METTL genes from cottons, Arabidopsis thaliana, and Homo sapiens could be classified into three clades (METTL3, METTL14, and METTL-like clades). Cis-element analysis predicated the possible functions of METTL genes in G. hirsutum. RNA-seq data revealed that GhMETTL14 (GH_A07G0817/GH_D07G0819) and GhMETTL3 (GH_A12G2586/GH_D12G2605) had high expressions in root, stem, leaf, torus, petal, stamen, pistil, and calycle tissues. GhMETTL14 also had the highest expression in 20 and 25 dpa fiber cells, implying a potential role at the cell wall thickening stage. Suppressing GhMETTL3 and GhMETTL14 by VIGS caused growth arrest and even death in G. hirsutum, along with decreased m6A abundance from the leaf tissues of VIGS plants. Overexpression of GhMETTL3 and GhMETTL14 produced distinct differentially expressed genes (DEGs) in A. thaliana, indicating their possible divergent functions after gene duplication. Overall, GhMETTLs play indispensable but divergent roles during the growth of cotton plants, which provides the basis for the systematic investigation of m6A in subsequent studies to improve the agronomic traits in cotton.

Signal transduction in non-climacteric fruit ripening

Fleshy fruit ripening involves changes in numerous cellular processes and metabolic pathways, resulting from the coordinated actions of diverse classes of structural and regulatory proteins. These include enzymes, transporters and complex signal transduction systems. Many aspects of the signaling machinery that orchestrates the ripening of climacteric fruits, such as tomato (Solanum lycopersicum), have been elucidated, but less is known about analogous processes in non-climacteric fruits. The latter include strawberry (Fragaria x ananassa) and grape (Vitis vinifera), both of which are used as non-climacteric fruit experimental model systems, although they originate from different organs: the grape berry is a true fruit derived from the ovary, while strawberry is an accessory fruit that is derived from the floral receptacle. In this article, we summarize insights into the signal transduction events involved in strawberry and grape berry ripening. We highlight the mechanisms underlying non-climacteric fruit ripening, the multiple primary signals and their integrated action, individual signaling components, pathways and their crosstalk, as well as the associated transcription factors and their signaling output.

Effect of VaMyb40 and VaMyb60 Overexpression on Stilbene Biosynthesis in Cell Cultures of Grapevine Vitis amurensis Rupr

Stilbenes are plant defense compounds known to rapidly accumulate in grapevine and some other plant species in response to microbial infection and several abiotic stresses. Stilbenes have attracted considerable attention due to valuable biological effects with multi-spectrum therapeutic application. However, there is a lack of information on natural signaling pathways and transcription factors regulating stilbene biosynthesis. It has been previously shown that MYB R2R3 transcription factor genes VaMyb40 and VaMyb60 were up-regulated in cell cultures of wild-growing grapevine Vitis amurensis Rupr. in response to UV irradiation. In this study, the effects of VaMyb40 or VaMyb60 overexpression in cell cultures of V. amurensis on their capability to produce stilbenes were investigated. Overexpression of the VaMyb60 gene led to a considerable increase in the content of stilbenes in three independently transformed transgenic lines in 5.9–13.9 times, while overexpression of the VaMyb40 gene also increased the content of stilbenes, although to a lesser extent (in 3.4–4.0 times) in comparison with stilbene levels in the empty vector-transformed calli. Stilbene content and stilbene production in the VaMyb60-transgenic calli reached 18.8 mg/g of dry weight (DW) and 150.8 mg/L, respectively. Using HPLC analysis, we detected eight individual stilbenes: t-resveratrol diglucoside, t-piceid, t-resveratrol, ε-viniferin, δ-viniferin, cis-resveratrol, cis-piceid, t-piceatannol. T-resveratrol prevailed over other stilbenoid compounds (53.1–89.5% of all stilbenes) in the VaMyb-overexpressing cell cultures. Moreover, the VaMyb40- and VaMyb60-transformed calli were capable of producing anthocyanins up to 0.035 mg/g DW, while the control calli did not produce anthocyanins. These findings show that the VaMyb40 and VaMyb60 genes positively regulate the stilbene biosynthesis as strong positive transcription regulators and can be used in biotechnological applications for stilbene production or high-quality viticulture and winemaking.

Comparative transcriptomics of tropical woody plants supports fast and furious strategy along the leaf economics spectrum in lianas

Lianas, climbing woody plants, influence the structure and function of tropical forests. Climbing traits have evolved multiple times, including ancestral groups such as gymnosperms and pteridophytes, but the genetic basis of the liana strategy is largely unknown. Here, we use a comparative transcriptomic approach for 47 tropical plant species, including ten lianas of diverse taxonomic origins, to identify genes that are consistently expressed or downregulated only in lianas. Our comparative analysis of full-length transcripts enabled the identification of a core interactomic network common to lianas. Sets of transcripts identified from our analysis reveal features related to functional traits pertinent to leaf economics spectrum in lianas, include upregulation of genes controlling epidermal cuticular properties, cell wall remodeling, carbon concentrating mechanism, cell cycle progression, DNA repair and a large suit of downregulated transcription factors and enzymes involved in ABA-mediated stress response as well as lignin and suberin synthesis. All together, these genes are known to be significant in shaping plant morphologies through responses such as gravitropism, phyllotaxy and shade avoidance.

Summary: The full-length fraction of liana transcriptomes mapped on a protein–protein interactome revealed the nature of their convergence through distinct sets of expressed and downregulated genes not observed in free-standing plants.

Direct regulation of shikimate, early phenylpropanoid, and stilbenoid pathways by Subgroup 2 R2R3-MYBs in grapevine

The stilbenoid pathway is responsible for the production of resveratrol in grapevine (Vitis vinifera L.). A few transcription factors (TFs) have been identified as regulators of this pathway but the extent of this control has not been deeply studied. Here we show how DNA affinity purification sequencing (DAP-Seq) allows for the genome-wide TF-binding site interrogation in grape. We obtained 5190 and 4443 binding events assigned to 4041 and 3626 genes for MYB14 and MYB15, respectively (approximately 40% of peaks located within −10 kb of transcription start sites). DAP-Seq of MYB14/MYB15 was combined with aggregate gene co-expression networks (GCNs) built from more than 1400 transcriptomic datasets from leaves, fruits, and flowers to narrow down bound genes to a set of high confidence targets. The analysis of MYB14, MYB15, and MYB13, a third uncharacterized member of Subgroup 2 (S2), showed that in addition to the few previously known stilbene synthase (STS) targets, these regulators bind to 30 of 47 STS family genes. Moreover, all three MYBs bind to several PAL, C4H, and 4CL genes, in addition to shikimate pathway genes, the WRKY03 stilbenoid co-regulator and resveratrol-modifying gene candidates among which ROMT2-3 were validated enzymatically. A high proportion of DAP-Seq bound genes were induced in the activated transcriptomes of transient MYB15-overexpressing grapevine leaves, validating our methodological approach for delimiting TF targets. Overall, Subgroup 2 R2R3-MYBs appear to play a key role in binding and directly regulating several primary and secondary metabolic steps leading to an increased flux towards stilbenoid production. The integration of DAP-Seq and reciprocal GCNs offers a rapid framework for gene function characterization using genome-wide approaches in the context of non-model plant species and stands up as a valid first approach for identifying gene regulatory networks of specialized metabolism.

Comprehensive identification and expression analysis of CRY gene family in Gossypium

BACKGROUND

The cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance.

RESULTS

In the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments.

CONCLUSION

We report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.

Anthocyanin and Flavonol Glycoside Metabolic Pathways Underpin Floral Color Mimicry and Contrast in a Sexually Deceptive Orchid

Sexually deceptive plants secure pollination by luring specific male insects as pollinators using a combination of olfactory, visual, and morphological mimicry. Flower color is a key component to this attraction, but its chemical and genetic basis remains poorly understood. Chiloglottis trapeziformis is a sexually deceptive orchid which has predominantly dull green-red flowers except for the central black callus projecting from the labellum lamina. The callus mimics the female of the pollinator and the stark color contrast between the black callus and dull green or red lamina is thought to enhance the visibility of the mimic. The goal of this study was to investigate the chemical composition and genetic regulation of temporal and spatial color patterns leading to visual mimicry, by integrating targeted metabolite profiling and transcriptomic analysis. Even at the very young bud stage, high levels of anthocyanins were detected in the dark callus, with peak accumulation by the mature bud stage. In contrast, anthocyanin levels in the lamina peaked as the buds opened and became reddish-green. Coordinated upregulation of multiple genes, including dihydroflavonol reductase and leucoanthocyanidin dioxygenase, and the downregulation of flavonol synthase genes (FLS) in the callus at the very young bud stage underpins the initial high anthocyanin levels. Conversely, within the lamina, upregulated FLS genes promote flavonol glycoside over anthocyanin production, with the downstream upregulation of flavonoid O-methyltransferase genes further contributing to the accumulation of methylated flavonol glycosides, whose levels peaked in the mature bud stage. Finally, the peak anthocyanin content of the reddish-green lamina of the open flower is underpinned by small increases in gene expression levels and/or differential upregulation in the lamina in select anthocyanin genes while FLS patterns showed little change. Differential expression of candidate genes involved in specific transport, vacuolar acidification, and photosynthetic pathways may also assist in maintaining the distinct callus and contrasting lamina color from the earliest bud stage through to the mature flower. Our findings highlight that flower color in this sexually deceptive orchid is achieved by complex tissue-specific coordinated regulation of genes and biochemical pathways across multiple developmental stages.

Genome-Wide Identification and Expression Analysis of the R2R3-MYB Transcription Factor Family Revealed Their Potential Roles in the Flowering Process in Longan (Dimocarpus longan)

Longan (Dimocarpus longan Lour.) is a productive fruit crop with high nutritional and medical value in tropical and subtropical regions. The MYB gene family is one of the most widespread plant transcription factor (TF) families participating in the flowering regulation. However, little is known about the MYB TFs involved in the flowering process in longan and its regulatory network. In this study, a total of 119 DlR2R3-MYB genes were identified in the longan genome and were phylogenetically grouped into 28 subgroups. The groupings were supported by highly conserved gene structures and motif composition of DlR2R3-MYB genes in each subgroup. Collinearity analysis demonstrated that segmental replications played a more crucial role in the expansion of the DlR2R3-MYB gene family compared to tandem duplications, and all tandem/segmental duplication gene pairs have evolved under purifying selection. Interspecies synteny analysis among longan and five representative species implied the occurrence of gene duplication events was one of the reasons contributing to functional differentiation among species. RNA-seq data from various tissues showed DlR2R3-MYB genes displayed tissue-preferential expression patterns. The pathway of flower development was enriched with six DlR2R3-MYB genes. Cis-acting element prediction revealed the putative functions of DlR2R3-MYB genes were related to the plant development, phytohormones, and environmental stresses. Notably, the orthologous counterparts between Arabidopsis and longan R2R3-MYB members tended to play conserved roles in the flowering regulation and stress responses. Transcriptome profiling on off-season flower induction (FI) by KClO₃ indicated two up-regulated and four down-regulated DlR2R3-MYB genes involved in the response to KClO₃ treatment compared with control groups. Additionally, qRT-PCR confirmed certain genes exhibited high expression in flowers/flower buds. Subcellular localization experiments revealed that three predicted flowering-associated MYB proteins were localized in the nucleus. Future functional studies on these potential candidate genes involved in the flowering development could further the understanding of the flowering regulation mechanism.

Vitis OneGenE: A Causality-Based Approach to Generate Gene Networks in Vitis vinifera Sheds Light on the Laccase and Dirigent Gene Families

The abundance of transcriptomic data and the development of causal inference methods have paved the way for gene network analyses in grapevine. Vitis OneGenE is a transcriptomic data mining tool that finds direct correlations between genes, thus producing association networks. As a proof of concept, the stilbene synthase gene regulatory network obtained with OneGenE has been compared with published co-expression analysis and experimental data, including cistrome data for MYB stilbenoid regulators. As a case study, the two secondary metabolism pathways of stilbenoids and lignin synthesis were explored. Several isoforms of laccase, peroxidase, and dirigent protein genes, putatively involved in the final oxidative oligomerization steps, were identified as specifically belonging to either one of these pathways. Manual curation of the predicted sequences exploiting the last available genome assembly, and the integration of phylogenetic and OneGenE analyses, identified a group of laccases exclusively present in grapevine and related to stilbenoids. Here we show how network analysis by OneGenE can accelerate knowledge discovery by suggesting new candidates for functional characterization and application in breeding programs.

The R2R3 Transcription Factor CsMYB59 Regulates Polyphenol Oxidase Gene CsPPO1 in Tea Plants (Camellia sinensis)

Polyphenol oxidase (PPO) plays a role in stress response, secondary metabolism, and other physiological processes during plant growth and development, and is also a critical enzyme in black tea production. However, the regulatory mechanisms of PPO genes and their activity in tea plants are still unclear. In this study, we measured PPO activity in two different tea cultivars, Taoyuandaye (TYDY) and Bixiangzao (BXZ), which are commonly used to produce black tea and green tea, respectively. The expression pattern of CsPPO1 was assessed and validated via transcriptomics and quantitative polymerase chain reaction in both tea varieties. In addition, we isolated and identified an R2R3-MYB transcription factor CsMYB59 that may regulate CsPPO1 expression. CsMYB59 was found to be a nuclear protein, and its expression in tea leaves was positively correlated with CsPPO1 expression and PPO activity. Transcriptional activity analysis showed that CsMYB59 was a transcriptional activator, and the dual-luciferase assay indicated that CsMYB59 could activate the expression of CsPPO1 in tobacco leaves. In summary, our study demonstrates that CsMYB59 represents a transcriptional activator in tea plants and may mediate the regulation of PPO activity by activating CsPPO1 expression. These findings provide novel insights into the regulatory mechanism of PPO gene in Camellia sinensis, which might help to breed tea cultivars with high PPO activity.

Metabolite profiling and transcriptome analyses reveal novel regulatory mechanisms of melatonin biosynthesis in hickory

Studies have shown that melatonin regulates the expression of various elements in the biosynthesis and catabolism of plant hormones. In contrast, the effects of these different plant hormones on the biosynthesis and metabolism of melatonin and their underlying molecular mechanisms are still unclear. In this study, the melatonin biosynthesis pathway was proposed from constructed metabolomic and transcriptomic libraries from hickory (Carya cathayensis Sarg.) nuts. The candidate pathway genes were further identified by phylogenetic analysis, amino-acid sequence alignment, and subcellular localization. Notably, most of the transcription factor-related genes coexpressed with melatonin pathway genes were hormone-responsive genes. Furthermore, dual-luciferase and yeast one-hybrid assays revealed that CcEIN3 (response to ethylene) and CcAZF2 (response to abscisic acid) could activate melatonin biosynthesis pathway genes, a tryptophan decarboxylase coding gene (CcTDC1) and an N-acetylserotonin methyltransferase coding gene (CcASMT1), by directly binding to their promoters, respectively. Our results provide a molecular basis for the characterization of novel melatonin biosynthesis regulatory mechanisms and demonstrate for the first time that abscisic acid and ethylene can regulate melatonin biosynthesis.

The grapevine (Vitis vinifera L.) floral transcriptome in Pinot noir variety: identification of tissue-related gene networks and whorl-specific markers in pre- and post-anthesis phases

The comprehension of molecular processes underlying the development and progression of flowering in plants is a hot topic, not only because that often the products of interest for human and animal nutrition are linked to the development of fruits or seeds, but also because the processes of gametes formation occurring in sexual organs are at the basis of recombination and genetic variability which constitutes the matter on which evolution acts, whether understood as natural or human driven. In the present study, we used an NGS approach to produce a grapevine flower transcriptome snapshot in different whorls and tissues including calyx, calyptra, filament, anther, stigma, ovary, and embryo in both pre- and post-anthesis phases. Our investigation aimed at identifying hub genes that unequivocally distinguish the different tissues providing insights into the molecular mechanisms that are at the basis of floral whorls and tissue development. To this end we have used different analytical approaches, some now consolidated in transcriptomic studies on plants, such as pairwise comparison and weighted-gene coexpression network analysis, others used mainly in studies on animals or human's genomics, such as the tau (τ) analysis aimed at isolating highly and absolutely tissue-specific genes. The intersection of data obtained by these analyses allowed us to gradually narrow the field, providing evidence about the molecular mechanisms occurring in those whorls directly involved in reproductive processes, such as anther and stigma, and giving insights into the role of other whorls not directly related to reproduction, such as calyptra and calyx. We believe this work could represent an important genomic resource for functional analyses of grapevine floral organ growth and fruit development shading light on molecular networks underlying grapevine reproductive organ determination.

Genome-wide analysis of MYB transcription factors of Vaccinium corymbosum and their positive responses to drought stress

Background

Blueberry (Vaccinium corymbosum L.) is an important species with a high content of flavonoids in fruits. As a perennial shrub, blueberry is characterized by shallow-rooted property and susceptible to drought stress. MYB transcription factor was reported to be widely involved in plant response to abiotic stresses, however, the role of MYB family in blueberry responding to drought stress remains elusive.

Results

In this study, we conducted a comprehensive analysis of VcMYBs in blueberry based on the genome data under drought stress, including phylogenetic relationship, identification of differentially expressed genes (DEGs), expression profiling, conserved motifs, expression correlation and protein-protein interaction prediction, etc. The results showed that 229 non-redundant MYB sequences were identified in the blueberry genome, and divided into 23 subgroups. A total of 102 MYB DEGs with a significant response to drought stress were identified, of which 72 in leaves and 69 in roots, and 8 differential expression genes with a > 20-fold change in the level of expression. 17 DEGs had a higher expression correlation with other MYB members. The interaction partners of the key VcMYB proteins were predicted by STRING analysis and in combination with physiological and morphological observation. 10 key VcMYB genes such as VcMYB8, VcMYB102 and VcMYB228 were predicted to be probably involved in reactive oxygen species (ROS) pathway, and 7 key VcMYB genes (VcMYB41, VcMYB88 and VcMYB100, etc..) probably participated in leaf regulation under drought treatment.

Conclusions

Our studies provide a new understanding of the regulation mechanism of VcMYB family in blueberry response to drought stress, and lay fundamental support for future studies on blueberry grown in regions with limited water supply for this crop.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07850-5.

Transcriptomic analysis of temporal shifts in berry development between two grapevine cultivars of the Pinot family reveals potential genes controlling ripening time

BACKGROUND

Grapevine cultivars of the Pinot family represent clonally propagated mutants with major phenotypic and physiological differences, such as different colour or shifted ripening time, as well as changes in important viticultural traits. Specifically, the cultivars 'Pinot Noir' (PN) and 'Pinot Noir Precoce' (PNP, early ripening) flower at the same time, but vary in the beginning of berry ripening (veraison) and, consequently, harvest time. In addition to genotype, seasonal climatic conditions (i.e. high temperatures) also affect ripening times. To reveal possible regulatory genes that affect the timing of veraison onset, we investigated differences in gene expression profiles between PN and PNP throughout berry development with a closely meshed time series and over two separate years.

RESULTS

The difference in the duration of berry formation between PN and PNP was quantified to be approximately two weeks under the growth conditions applied, using plant material with a proven PN and PNP clonal relationship. Clusters of co-expressed genes and differentially expressed genes (DEGs) were detected which reflect the shift in the timing of veraison onset. Functional annotation of these DEGs fit to observed phenotypic and physiological changes during berry development. In total, we observed 3,342 DEGs in 2014 and 2,745 DEGs in 2017 between PN and PNP, with 1,923 DEGs across both years. Among these, 388 DEGs were identified as veraison-specific and 12 were considered as berry ripening time regulatory candidates. The expression profiles revealed two candidate genes for ripening time control which we designated VviRTIC1 and VviRTIC2 (VIT_210s0071g01145 and VIT_200s0366g00020, respectively). These genes likely contribute the phenotypic differences observed between PN and PNP.

CONCLUSIONS

Many of the 1,923 DEGs show highly similar expression profiles in both cultivars if the patterns are aligned according to developmental stage. In our work, putative genes differentially expressed between PNP and PN which could control ripening time as well as veraison-specific genes were identified. We point out connections of these genes to molecular events during berry development and discuss potential candidate genes which may control ripening time. Two of these candidates were observed to be differentially expressed in the early berry development phase. Several down-regulated genes during berry ripening are annotated as auxin response factors / ARFs. Conceivably, general changes in auxin signaling may cause the earlier ripening phenotype of PNP.

VqMYB154 promotes polygene expression and enhances resistance to pathogens in Chinese wild grapevine

Resveratrol plays a crucial phytoalexin role in the grapevine and is beneficial to human health. However, the molecular mechanism of resveratrol accumulation in the enhancement of disease resistance is unclear. Here, we report that the transcription factor VqMYB154 from Vitis quinquangularis accession Danfeng-2 is strongly expressed under artificial inoculation with Uncinula necator and regulates resveratrol accumulation. Unlike its homolog, VqMYB154 has a pathogen-induced promoter and responds to stimulation by U. necator, Pseudomonas syringae, and other treatments. Yeast one-hybrid and GUS activity assays confirmed that VqMYB154 can activate the stilbene synthase genes VqSTS9, VqSTS32, and VqSTS42 by directly binding to their promoters. Overexpression of VqMYB154 in grape leaves resulted in activation of the stilbene pathway, upregulation of STS genes, and accumulation of stilbenoids. In addition, heterologous overexpression of VqMYB154 in Arabidopsis activated resistance-related genes and resulted in greater programmed cell death and accumulation of reactive oxygen species, which led to resistance against P. syringae. These results suggest that the transcription factor VqMYB154 from V. quinquangularis accession Danfeng-2 participates in the regulatory mechanism that improves the biosynthesis and accumulation of stilbenes and enhances resistance to disease in grapevine.

The MdHY5-MdWRKY41-MdMYB transcription factor cascade regulates the anthocyanin and proanthocyanidin biosynthesis in red-fleshed apple

Red-fleshed apple fruits are popular because of their high flavonoid content. Although MdMYB10 and its homologs have been identified as crucial regulators of the fruit coloring process, other transcription factors (TFs) contributing to the differences in flesh coloration have not been fully characterized. In this study, we investigated the regulatory effects of MdWRKY41 on anthocyanin and proanthocyanidin (PA) synthesis in red-fleshed apples. The overexpression of MdWRKY41 in red-fleshed apple calli inhibited anthocyanin and PA accumulation by downregulating the expression of a MYB TF gene (MdMYB12) and specific structural genes (MdLAR, MdUFGT, and MdANR). Furthermore, MdWRKY41 was shown to interact with MdMYB16 to form a complex that can further suppress MdANR and MdUFGT expression. Interestingly, MdWRKY41 was targeted by the photoresponse factor MdHY5 and inhibited its transcription. Overall, our findings provide insights into a novel MdHY5-MdWRKY41-MdMYB regulatory module influencing anthocyanin and PA synthesis in red-fleshed apple fruits.

Genetic variation and association analyses identify genes linked to fruit set-related traits in grapevine

Grapevine is one of the most valuable fruit crops in the world. Adverse environmental conditions reduce fruit quality and crop yield, so understanding the genetic and molecular mechanisms determining crop yield components is essential to optimize grape production. The analysis of a diverse collection of grapevine cultivars allowed us to evaluate the relationship between fruit set-related components of yield, including the incidence of reproductive disorders such as coulure and millerandage. The collection displayed a great phenotypic variation that we surveyed in a genetics association study using 15,309 single nucleotide polymorphisms (SNPs) detected in the sequence of 289 candidate genes scattered across the 19 grapevine linkage groups. After correcting statistical models for population structure and linkage disequilibrium effects, 164 SNPs from 34 of these genes were found to associate with fruit set-related traits, supporting a complex polygenic determinism. Many of them were found in the sequence of different putative MADS-box transcription factors, a gene family related with plant reproductive development control. In addition, we observed an additive effect of some of the associated SNPs on the phenotype, suggesting that advantageous alleles from different loci could be pyramided to generate superior cultivars with optimized fruit production.

R2R3-MYB transcription factor family in tea plant (Camellia sinensis): Genome-wide characterization, phylogeny, chromosome location, structure and expression patterns

MYB transcription factors play essential roles in many biological processes and environmental stimuli. However, the functions of the MYB transcription factor family in tea plants have not been elucidated. Here, a total of 122 CsR2R3-MYB genes were identified from the chromosome level genome of tea plant (Camellia sinensis). The CsR2R3-MYB genes were phylogenetically classified into 25 groups. Results from the structure analysis of the gene, conserved motifs, and chromosomal distribution supported the relative conservation of the R2R3-MYB genes family in the tea plant. Synteny analysis indicated that 122, 34, and 112 CsR2R3-MYB genes were orthologous to Arabidopsis thaliana, Oryza sativa and C. sinensis var. 'huangdan' (HD), respectively. Tissue-specific expression showed that all CsR2R3-MYB genes had different expression patterns in the tea plant tissues, indicating that these genes may perform diverse functions. The expression patterns of representative R2R3-MYB genes and the regulatory network of the main anthocyanin components were analyzed, which suggested that CsMYB17 may played a key role in the regulation of cya-3-O-gal, del-3-O-gal, cya-3-O-glu and pel-3-O-glu. Results from the qRT-PCR validation of selected genes suggested that CsR2R3-MYB genes were induced in response to drought, cold, GA, and ABA treatments. Overall, this study provides comprehensive and systematic information for research on the function of R2R3-MYB genes in tea plants.

Genomic analysis uncovers functional variation in the C-terminus of anthocyanin-activating MYB transcription factors

MYB transcription factors regulate diverse aspects of plant development and secondary metabolism, often by partnering in transcriptional regulatory complexes. Here, we harness genomic resources to identify novel MYBs, thereby producing an updated eudicot MYB phylogeny with revised relationships among subgroups as well as new information on sequence variation in the disordered C-terminus of anthocyanin-activating MYBs. BLAST® and hidden Markov model scans of gene annotations identified a total of 714 MYB transcription factors across the genomes of four crops that span the eudicots: apple, grape, kiwifruit and tomato. Codon model-based phylogenetic inference identified novel members of previously defined subgroups, and the function of specific anthocyanin-activating subgroup 6 members was assayed transiently in tobacco leaves. Sequence conservation within subgroup 6 highlighted one previously described and two novel short linear motifs in the disordered C-terminal region. The novel motifs have a mix of hydrophobic and acidic residues and are predicted to be relatively ordered compared with flanking protein sequences. Comparison of motifs with the Eukaryotic Linear Motif database suggests roles in protein-protein interaction. Engineering of motifs and their flanking regions from strong anthocyanin activators into weak activators, and vice versa, affected function. We conclude that, although the MYB C-terminal sequence diverges greatly even within MYB clades, variation within the C-terminus at and near relatively ordered regions offers opportunities for exploring MYB function and developing superior alleles for plant breeding.

Characterization of cotton ARF factors and the role of GhARF2b in fiber development

BACKGROUND

Cotton fiber is a model system for studying plant cell development. At present, the functions of many transcription factors in cotton fiber development have been elucidated, however, the roles of auxin response factor (ARF) genes in cotton fiber development need be further explored.

RESULTS

Here, we identify auxin response factor (ARF) genes in three cotton species: the tetraploid upland cotton G. hirsutum, which has 73 ARF genes, and its putative extent parental diploids G. arboreum and G. raimondii, which have 36 and 35 ARFs, respectively. Ka and Ks analyses revealed that in G. hirsutum ARF genes have undergone asymmetric evolution in the two subgenomes. The cotton ARFs can be classified into four phylogenetic clades and are actively expressed in young tissues. We demonstrate that GhARF2b, a homolog of the Arabidopsis AtARF2, was preferentially expressed in developing ovules and fibers. Overexpression of GhARF2b by a fiber specific promoter inhibited fiber cell elongation but promoted initiation and, conversely, its downregulation by RNAi resulted in fewer but longer fiber. We show that GhARF2b directly interacts with GhHOX3 and represses the transcriptional activity of GhHOX3 on target genes.

CONCLUSION

Our results uncover an important role of the ARF factor in modulating cotton fiber development at the early stage.

Cytokinin but not gibberellin application had major impact on the phenylpropanoid pathway in grape

Cytokinin and gibberellic acid (GA) are growth regulators used to increase berry size in seedless grapes and it is of interest to understand their effects on the phenylpropanoid pathway and on ripening processes. GA₃ and synthetic cytokinin forchlorfenuron (N-(2-chloro-4-pyridyl)-N'-phenylurea, CPPU) and their combination were applied to 6 mm diameter fruitlets of 'Sable Seedless', and berries were sampled 51 and 70 days (d) following application. All treatments increased berry size and delayed sugar accumulation and acid degradation with a stronger effect of CPPU. CPPU, but not GA, reduced berry color and the levels of anthocyanins. While CPPU reduced the levels of anthocyanins by more than 50%, the combined treatment of GA+CPPU reduced the levels by about 25% at 51 d. CPPU treatment had minor effects on flavonols content but increased the levels of monomeric flavan-3-ols by more than two-fold. Phloroglucinol analysis using HPLC showed that proanthocyanidin content was significantly increased by CPPU, whereas mean degree of polymerization was reduced from 26 to 19. Volatile analysis by GC-MS showed changes in composition with CPPU or GA treatment with potential impact on flavor. RNA-seq analysis showed that GA had a minor overall effect on the transcriptome whereas CPPU had pronounced effects on gene expression at both 51 and 70 d. Comparing the control and CPPU at similar Brix of ca. 19.7°, a reduced expression of stilbene synthases (STSs) including their regulators MYB14 and MYB15, and other phenylpropanoid-related genes was observed in CPPU-treated grapes. Overall, our study shows that CPPU had a major influence on the phenylpropanoid pathway and affected multiple ripening-related processes.

Xylella fastidiosa causes transcriptional shifts that precede tylose formation and starch depletion in xylem

Pierce's disease (PD) in grapevine (Vitis vinifera) is caused by the bacterial pathogen Xylella fastidiosa. X. fastidiosa is limited to the xylem tissue and following infection induces extensive plant-derived xylem blockages, primarily in the form of tyloses. Tylose-mediated vessel occlusions are a hallmark of PD, particularly in susceptible V. vinifera. We temporally monitored tylose development over the course of the disease to link symptom severity to the level of tylose occlusion and the presence/absence of the bacterial pathogen at fine-scale resolution. The majority of vessels containing tyloses were devoid of bacterial cells, indicating that direct, localized perception of X. fastidiosa was not a primary cause of tylose formation. In addition, we used X-ray computed microtomography and machine-learning to determine that X. fastidiosa induces significant starch depletion in xylem ray parenchyma cells. This suggests that a signalling mechanism emanating from the vessels colonized by bacteria enables a systemic response to X. fastidiosa infection. To understand the transcriptional changes underlying these phenotypes, we integrated global transcriptomics into the phenotypes we tracked over the disease spectrum. Differential gene expression analysis revealed that considerable transcriptomic reprogramming occurred during early PD before symptom appearance. Specifically, we determined that many genes associated with tylose formation (ethylene signalling and cell wall biogenesis) and drought stress were up-regulated during both Phase I and Phase II of PD. On the contrary, several genes related to photosynthesis and carbon fixation were down-regulated during both phases. These responses correlate with significant starch depletion observed in ray cells and tylose synthesis in vessels.

Genome-Wide Comparative Analysis of R2R3 MYB Gene Family in Populus and Salix and Identification of Male Flower Bud Development-Related Genes

The MYB transcription factor (TF) family is one of the largest plant transcription factor gene family playing vital roles in plant growth and development, including defense, cell differentiation, secondary metabolism, and responses to biotic and abiotic stresses. As a model tree species of woody plants, in recent years, the identification and functional prediction of certain MYB family members in the poplar genome have been reported. However, to date, the characterization of the gene family in the genome of the poplar's sister species willow has not been done, nor are the differences and similarities between the poplar and willow genomes understood. In this study, we conducted the first genome-wide investigation of the R2R3 MYB subfamily in the willow, identifying 216 R2R3 MYB gene members, and combined with the poplar R2R3 MYB genes, performed the first comparative analysis of R2R3 MYB genes between the poplar and willow. We identified 81 and 86 pairs of R2R3 MYB paralogs in the poplar and willow, respectively. There were 17 pairs of tandem repeat genes in the willow, indicating active duplication of willow R2R3 MYB genes. A further 166 pairs of poplar and willow orthologs were identified by collinear and synonymous analysis. The findings support the duplication of R2R3 MYB genes in the ancestral species, with most of the R2R3 MYB genes being retained during the evolutionary process. The phylogenetic trees of the R2R3 MYB genes of 10 different species were drawn. The functions of the poplar and willow R2R3 MYB genes were predicted using reported functional groupings and clustering by OrthoFinder. Identified 5 subgroups in general expanded in woody species, three subgroups were predicted to be related to lignin synthesis, and we further speculate that the other two subgroups also play a role in wood formation. We analyzed the expression patterns of the GAMYB gene of subgroup 18 (S18) related to pollen development in the male flower buds of poplar and willow at different developmental stages by qRT-PCR. The results showed that the GAMYB gene was specifically expressed in the male flower bud from pollen formation to maturity, and that the expression first increased and then decreased. Both the specificity of tissue expression specificity and conservation indicated that GAMYB played an important role in pollen development in both poplar and willow and was an ideal candidate gene for the analysis of male flower development-related functions of the two species.

Grape Berry Secondary Metabolites and Their Modulation by Abiotic Factors in a Climate Change Scenario-A Review

Temperature, water, solar radiation, and atmospheric CO₂ concentration are the main abiotic factors that are changing in the course of global warming. These abiotic factors govern the synthesis and degradation of primary (sugars, amino acids, organic acids, etc.) and secondary (phenolic and volatile flavor compounds and their precursors) metabolites directly, via the regulation of their biosynthetic pathways, or indirectly, via their effects on vine physiology and phenology. Several hundred secondary metabolites have been identified in the grape berry. Their biosynthesis and degradation have been characterized and have been shown to occur during different developmental stages of the berry. The understanding of how the different abiotic factors modulate secondary metabolism and thus berry quality is of crucial importance for breeders and growers to develop plant material and viticultural practices to maintain high-quality fruit and wine production in the context of global warming. Here, we review the main secondary metabolites of the grape berry, their biosynthesis, and how their accumulation and degradation is influenced by abiotic factors. The first part of the review provides an update on structure, biosynthesis, and degradation of phenolic compounds (flavonoids and non-flavonoids) and major aroma compounds (terpenes, thiols, methoxypyrazines, and C13 norisoprenoids). The second part gives an update on the influence of abiotic factors, such as water availability, temperature, radiation, and CO₂ concentration, on berry secondary metabolism. At the end of the paper, we raise some critical questions regarding intracluster berry heterogeneity and dilution effects and how the sampling strategy can impact the outcome of studies on the grapevine berry response to abiotic factors.

Phytostilbenes as agrochemicals: biosynthesis, bioactivity, metabolic engineering and biotechnology

Covering: 1976 to 2020. Although constituting a limited chemical family, phytostilbenes represent an emblematic group of molecules among natural compounds. Ever since their discovery as antifungal compounds in plants and their ascribed role in human health and disease, phytostilbenes have never ceased to arouse interest for researchers, leading to a huge development of the literature in this field. Owing to this, the number of references to this class of compounds has reached the tens of thousands. The objective of this article is thus to offer an overview of the different aspects of these compounds through a large bibliography analysis of more than 500 articles. All the aspects regarding phytostilbenes will be covered including their chemistry and biochemistry, regulation of their biosynthesis, biological activities in plants, molecular engineering of stilbene pathways in plants and microbes as well as their biotechnological production by plant cell systems.

RNA-Seq Time Series of Vitis vinifera Bud Development Reveals Correlation of Expression Patterns with the Local Temperature Profile

Plants display sophisticated mechanisms to tolerate challenging environmental conditions and need to manage their ontogenesis in parallel. Here, we set out to generate an RNA-Seq time series dataset throughout grapevine (Vitis vinifera) early bud development. The expression of the developmental regulator VviAP1 served as an indicator of the progression of development. We investigated the impact of changing temperatures on gene expression levels during the time series and detected a correlation between increased temperatures and a high expression level of genes encoding heat-shock proteins. The dataset also allowed the exemplary investigation of expression patterns of genes from three transcription factor (TF) gene families, namely MADS-box, WRKY, and R2R3-MYB genes. Inspection of the expression profiles from all three TF gene families indicated that a switch in the developmental program takes place in July which coincides with increased expression of the bud dormancy marker gene VviDRM1.

Drought activates MYB41 orthologs and induces suberization of grapevine fine roots

The permeability of roots to water and nutrients is controlled through a variety of mechanisms and one of the most conspicuous is the presence of the Casparian strips and suberin lamellae. Roots actively regulate the creation of these structures developmentally, along the length of the root, and in response to the environment, including drought. In the current study, we characterized the suberin composition along the length of grapevine fine roots during development and in response to water deficit, and in the same root systems we quantified changes in expression of suberin biosynthesis- and deposition-related gene families (via RNAseq) allowing the identification of drought-responsive suberin-related genes. Grapevine suberin composition did not differ between primary and lateral roots, and was similar to that of other species. Under water deficit there was a global upregulation of suberin biosynthesis which resulted in an increase of suberin specific monomers, but without changes in their relative abundances, and this upregulation took place across all the developmental stages of fine roots. These changes corresponded to the upregulation of numerous suberin biosynthesis- and export-related genes which included orthologs of the previously characterized AtMYB41 transcriptional factor. Functional validation of two grapevine MYB41 orthologs, VriMYB41 and VriMYB41-like, confirmed their ability to globally upregulate suberin biosynthesis, export, and deposition. This study provides a detailed characterization of the developmental and water deficit induced suberization of grapevine fine roots and identifies important orthologs responsible for suberin biosynthesis, export, and its regulation in grape.

Genome-wide characterization of the GRF family and their roles in response to salt stress in Gossypium

BACKGROUND

Cotton (Gossypium spp.) is the most important world-wide fiber crop but salt stress limits cotton production in coastal and other areas. Growth regulation factors (GRFs) play regulatory roles in response to salt stress, but their roles have not been studied in cotton under salt stress.

RESULTS

We identified 19 GRF genes in G. raimondii, 18 in G. arboreum, 34 in G. hirsutum and 45 in G. barbadense, respectively. These GRF genes were phylogenetically analyzed leading to the recognition of seven GRF clades. GRF genes from diploid cottons (G. raimondii and G. arboreum) were largely retained in allopolyploid cotton, with subsequent gene expansion in G. barbadense relative to G. hirsutum. Most G. hirsutum GRF (GhGRF) genes are preferentially expressed in young and growing tissues. To explore their possible role in salt stress, we used qRT-PCR to study expression responses to NaCl treatment, showing that five GhGRF genes were down-regulated in leaves. RNA-seq experiments showed that seven GhGRF genes exhibited decreased expression in leaves under NaCl treatment, three of which (GhGRF3, GhGRF4, and GhGRF16) were identified by both RNA-seq and qRT-PCR. We also identified six and three GRF genes that exhibit decreased expression under salt stress in G. arboreum and G. barbadense, respectively. Consistent with its lack of leaf withering or yellowing under the salt treatment conditions, G. arboreum had better salt tolerance than G. hirsutum and G. barbadense. Our results suggest that GRF genes are involved in salt stress responses in Gossypium.

CONCLUSION

In summary, we identified candidate GRF genes that were involved in salt stress responses in cotton.

A Novel R2R3-MYB Transcription Factor PqMYB4 Inhibited Anthocyanin Biosynthesis in Paeonia qiui

Paeonia qiui is a wild tree peony native to China. Its leaves show a clear purple-red color from the germination to the flowering stage, and it has high leaf-viewing value. A MYB transcription factor gene, designated as PqMYB4, was isolated from leaves of P. qiui based on transcriptome datas. The full-length cDNA of PqMYB4 was 693 bp, encoding 230 amino acids. Sequence alignment and phylogenetic analysis revealed that PqMYB4 was a R2R3-MYB transcription factor clustered with AtMYB4 in Arabidopsis thaliana. Moreover, it contained a C1 motif, an EAR repression motif and a TLLLFR motif in the C-terminal domains, which were unique in transcription repression MYB. Subcellular location analysis showed that PqMYB4 was located in the cell nucleus. PqMYB4 was highly expressed in the late stage of leaf development, and was negatively correlated with the anthocyanin content. The petiole of wild-type Arabidopsis seedlings was deeper in color than the transgenic lines of PqMYB4 and showed a little purple-red color. The seed coat color of Arabidopsis seeds that overexpressed PqMYB4 gene was significantly lighter than that of wild-type seeds. In transgenic Arabidopsis, the expression level of AtCHS, AtCHI, AtDFR and AtANS were down-regulated significantly. These results showed that PqMYB4 was involved in the negative regulation of anthocyanin biosynthesis in tree peony leaves, which can control the anthocyanin pathway genes. Together, these findings provide a valuable resource with which to further study the regulatory mechanism of anthocyanin biosynthesis in the leaf of P. qiui. They also benefit the molecular breeding of tree peony cultivars with colored leaf.

Network aggregation improves gene function prediction of grapevine gene co-expression networks

Aggregation across multiple networks highlights robust co-expression interactions and improves the functional connectivity of grapevine gene co-expression networks. In recent years, the rapid accumulation of transcriptome datasets from diverse experimental conditions has enabled the widespread use of gene co-expression network (GCN) analysis in plants. In grapevine, GCN analysis has shown great promise for gene function prediction, however, measurable progress is currently lacking. Using accumulated microarray datasets from the grapevine whole-genome array (33 experiments, 1359 samples), we explored how meta-analysis through aggregation influences the functional connectivity (performance) of derived networks using guilt-by-association neighbor voting. Two annotation schemes, i.e. MapMan BIN and Pfam, at two sparsity thresholds, i.e. top 100 (stringent) and 300 (relaxed) ranked genes were evaluated. We observed that aggregating across multiple networks improves performance dramatically, with the aggregate outperforming the majority of functional terms across individual networks. Network sparsity and size (i.e. the number of samples and aggregates) were key factors influencing performance while the choice of annotation scheme had little. Systematic comparison with various state-of-the-art microarray and RNA-seq networks was also performed, however, none outperformed the aggregate microarray network despite having good predictive performance. Repeating these series of tests using a functional enrichment-based performance metric also showed remarkably consistent findings with guilt-by-association neighbor voting. To demonstrate its functionality, we explore the function and transcriptional regulation of grapevine EXPANSIN genes. We envisage that network aggregation will offer new and unique opportunities for gene function prediction in future grapevine functional genomics studies. To this end, we make the aggregate networks and associated metadata publicly available at VTC-Agg (https://sites.google.com/view/vtc-agg).

The WRKY53 transcription factor enhances stilbene synthesis and disease resistance by interacting with MYB14 and MYB15 in Chinese wild grape

Resveratrol is notable not only for its functions in disease resistance in plants but also for its health benefits when it forms part of the human diet. Identification of new transcription factors helps to reveal the regulatory mechanisms of stilbene synthesis. Here, the WRKY53 transcription factor was isolated from the Chinese wild grape, Vitis quinquangularis. Vqwrky53 was expressed in a variety of tissues and responded to powdery mildew infection and to exogenous hormone application. VqWRKY53 was located in the nucleus and had transcriptional activation activity in yeast. A yeast two-hybrid assay and a bimolecular fluorescence complementation assay confirmed that VqWRKY53 interacted physically with VqMYB14 and VqMYB15, which have previously been reported to regulate stilbene synthesis. When Vqwrky53 was overexpressed in grape leaves, the expression of VqSTS32 and VqSTS41 and the content of stilbenes were increased. A yeast one-hybrid assay demonstrated that VqWRKY53 could bind directly to the promoters of STS genes. Overexpression of Vqwrky53 activated β-glucuronidase expression, driven by STS promoters, and co-expressing Vqwrky53 with VqMYB14 and VqMYB15 showed stronger regulatory functions. Heterologous overexpression of Vqwrky53 in Arabidopsis accelerated leaf senescence and disease resistance to PstDC3000.

Combined Analysis of the Metabolome and Transcriptome Identified Candidate Genes Involved in Phenolic Acid Biosynthesis in the Leaves of Cyclocarya paliurus

To assess changes of metabolite content and regulation mechanism of the phenolic acid biosynthesis pathway at different developmental stages of leaves, this study performed a combined metabolome and transcriptome analysis of Cyclocarya paliurus leaves at different developmental stages. Metabolite and transcript profiling were conducted by ultra-performance liquid chromatography quadrupole time-of-flight tandem mass spectrometer and high-throughput RNA sequencing, respectively. Transcriptome identification showed that 58 genes were involved in the biosynthesis of phenolic acid. Among them, 10 differentially expressed genes were detected between every two developmental stages. Identification and quantification of metabolites indicated that 14 metabolites were located in the phenolic acid biosynthetic pathway. Among them, eight differentially accumulated metabolites were detected between every two developmental stages. Association analysis between metabolome and transcriptome showed that six differentially expressed structural genes were significantly positively correlated with metabolite accumulation and showed similar expression trends. A total of 128 transcription factors were identified that may be involved in the regulation of phenolic acid biosynthesis; these include 12 MYBs and 10 basic helix–loop–helix (bHLH) transcription factors. A regulatory network of the phenolic acid biosynthesis was established to visualize differentially expressed candidate genes that are involved in the accumulation of metabolites with significant differences. The results of this study contribute to the further understanding of phenolic acid biosynthesis during the development of leaves of C. paliurus.

The Effect of Transcription Factor MYB14 on Defense Mechanisms in Vitis quinquangularis-Pingyi

In the current study, we identified a transcription factor, MYB14, from Chinese wild grape, Vitis quinquangularis-Pingyi (V. quinquangularis-PY), which could enhance the main stilbene contents and expression of stilbene biosynthesis genes (StSy/RS) by overexpression of VqMYB14. The promoter of VqMYB14 (pVqMYB14) was shown to be induced as part of both basal immunity (also called pathogen-associated molecular pattern (PAMP)-triggered immunity, PTI) and effector-triggered immunity (ETI), triggered by the elicitors flg22 and harpin, respectively. This was demonstrated by expression of pVqMYB14 in Nicotiana benthamiana and Vitis. We identified sequence differences, notably an 11 bp segment in pVqMYB14 that is important for the PTI/ETI, and particularly for the harpin-induced ETI response. In addition, we showed that activation of the MYB14 promoter correlates with differences in the expression of MYB14 and stilbene pattern induced by flg22 and harpin. An experimental model of upstream signaling in V. quinquangularis-PY is presented, where early defense responses triggered by flg22 and harpin partially overlap, but where the timing and levels differ. This translates into a qualitative difference with respect to patterns of stilbene accumulation.

Genetic analysis of a white-to-red berry skin color reversion and its transcriptomic and metabolic consequences in grapevine (Vitis vinifera cv. 'Moscatel Galego')

BACKGROUND

Somatic mutations occurring within meristems of vegetative propagation material have had a major role in increasing the genetic diversity of the domesticated grapevine (Vitis vinifera subsp. vinifera). The most well studied somatic variation in this species is the one affecting fruit pigmentation, leading to a plethora of different berry skin colors. Color depletion and reversion are often observed in the field. In this study we analyzed the origin of a novel white-to-red skin color reversion and studied its possible metabolic and transcriptomic consequences on cv. 'Muscat à Petits Grains Blancs' (synonym cv. 'Moscatel Galego Branco'), a member of the large family of Muscats.

RESULTS

The mild red-skinned variant (cv. 'Muscat à Petits Grains Rouge', synonym cv. 'Moscatel Galego Roxo'), characterized by a preferential accumulation of di-hydroxylated anthocyanins, showed in heterozygosis a partially-excised Gret1 retrotransposon in the promoter region of the MYBA1 anthocyanin regulator, while MYBA2 was still in homozygosis for its non-functional allele. Through metabolic (anthocyanin, resveratrol and piceid quantifications) and transcriptomic (RNA-Seq) analyses, we show that within a near-isogenic background, the transcriptomic consequences of color reversion are largely associated to diminished light/UV-B responses probably as a consequence of the augment of metabolic sunscreens (i.e. anthocyanins).

CONCLUSIONS

We propose that the reduced activity of the flavonoid tri-hydroxylated sub-branch and decreased anthocyanin synthesis and modification (e.g. methylation and acylation) are the potential causes for the mild red-skinned coloration in the pigmented revertant. The observed positive relation between anthocyanins and stilbenes could be attributable to an increased influx of phenylpropanoid intermediaries due to the replenished activity of MYBA1, an effect yet to be demonstrated in other somatic variants.

Molecular Analysis of UV-C Induced Resveratrol Accumulation in Polygonum cuspidatum Leaves

Resveratrol is one of the most studied plant secondary metabolites owing to its numerous health benefits. It is accumulated in some plants following biotic and abiotic stress pressures, including UV-C irradiation. Polygonum cuspidatum represents the major natural source of concentrated resveratrol but the underlying mechanisms as well as the effects of UV-C irradiation on resveratrol content have not yet been documented. Herein, we found that UV-C irradiation significantly increased by 2.6-fold and 1.6-fold the resveratrol content in irradiated leaf samples followed by a dark incubation for 6 h and 12 h, respectively, compared to the untreated samples. De novo transcriptome sequencing and assembly resulted into 165,013 unigenes with 98 unigenes mapped to the resveratrol biosynthetic pathway. Differential expression analysis showed that P.cuspidatum strongly induced the genes directly involved in the resveratrol synthesis, including phenylalanine ammonia-lyase, cinnamic acid 4-hydroxylase, 4-coumarate-CoA ligase and stilbene synthase (STS) genes, while strongly decreased the chalcone synthase (CHS) genes after exposure to UV-C. Since CHS and STS share the same substrate, P. cuspidatum tends to preferentially divert the substrate to the resveratrol synthesis pathway under UV-C treatment. We identified several members of the MYB, bHLH and ERF families as potential regulators of the resveratrol biosynthesis genes.

Coordinated Regulation of Grape Berry Flesh Color by Transcriptional Activators and Repressors

Yan73 is a teinturier Vitis vinifera variety with red berry flesh, but the molecular mechanisms underlying its flesh coloration remain unclear. We analyzed the flavonoid metabolic and transcriptome profiles of Yan73 berry red and white flesh using HPLC-ESI-MS/MS and RNA-sequencing technologies. Anthocyanins are the main flavonoids responsible for Yan73 berry flesh color, and the coloration is coordinately regulated by the VvMYBA1 transcriptional activator and VvMYBC2-L1 transcriptional repressor. Furthermore, yeast one- and two-hybrid, dual luciferase, and bimolecular fluorescence complementation assays suggested that VvMYBA1 positively regulates Yan73 berry flesh color via interactions with VvWDR1 and the activation of the VvCHI3, VvOMT, and VvGST4 promoters, whereas VvMYBC2-L1 negatively regulates Yan73 berry flesh color, possibly by competing with the R2R3-MYB transcriptional activators for bHLH partners or by repressing VvOMT and VvGST4 expression. Our findings provide new insights into the molecular mechanisms regulating grape flesh color.