Genome-wide transcriptome profiling of transgenic hop (Humulus lupulus L.) constitutively overexpressing HlWRKY1 and HlWDR1 transcription factors

CRISPR/Cas9-mediated mutagenesis of the mediator complex subunits MED5a and MED5b genes impaired secondary metabolite accumulation in hop (Humulus lupulus)

Hop (Humulus lupulus L.) is an important commercial crop known for the biosynthesis of valuable specialized secondary metabolites in glandular trichomes (lupulin glands), which are used for the brewing industry. To achieve burgeoning market demands is the essentiality of comprehensive understanding of the mechanisms of biosynthesis of secondary metabolites in hop. Over the past year, several studies using structural biology and functional genomics approaches have shown that Mediator (MED) serves as an integrative hub for RNAP II-mediated transcriptional regulation of various physiological and cellular processes, including involvement of MED5a and MED5b in hyperaccumulation of phenylpropanoid in A. thaliana. In the present work, an unprecedented attempt was made to generate Hlmed5a/med5b double loci mutant lines in hop using a CRISPR/Cas9-based genome editing system. The Hlmed5a/med5b double loci mutant lines showed reduced expression of structural genes of the flavonoid, humulone, and terpenoid biosynthetic pathways, which was more pronounced in the lupulin gland compared to leaf tissue and was consistent with their reduced accumulation. Phenotypic and anatomical observations revealed that Hlmed5a/med5b double loci mutant line exhibited robust growth, earlier flowering, earlier cone maturity, reduced cone size, variations in floral structure patterns, and distorted lupulin glands without any remarkable changes in leaf morphology, intensity of leaf color, and chlorophyll content. Comparative transcriptome analysis of leaf and lupulin gland tissues indicates that the expression of enzymatic genes related to secondary metabolite biosynthesis, phytohormone biosynthesis, floral organs, flowering time, and trichome development, including other genes related to starch and sucrose metabolism and defense mechanisms, were differentially modulated in the Hlmed5a/med5b lines. The combined results from functional and transcriptomic analyses illuminates the pivotal function of HlMED5a and HlMED5b in homeostasis of secondary meatbolites accumulation in hop.

Advances in the study of the function and mechanism of the action of flavonoids in plants under environmental stresses

MAIN CONCLUSION

This review summarizes the anti-stress effects of flavonoids in plants and highlights its role in the regulation of polar auxin transport and free radical scavenging mechanism. As secondary metabolites widely present in plants, flavonoids play a vital function in plant growth, but also in resistance to stresses. This review introduces the classification, structure and synthetic pathways of flavonoids. The effects of flavonoids in plant stress resistance were enumerated, and the mechanism of flavonoids in plant stress resistance was discussed in detail. It is clarified that plants under stress accumulate flavonoids by regulating the expression of flavonoid synthase genes. It was also determined that the synthesized flavonoids are transported in plants through three pathways: membrane transport proteins, vesicles, and bound to glutathione S-transferase (GST). At the same time, the paper explores that flavonoids regulate polar auxin transport (PAT) by acting on the auxin export carrier PIN-FORMED (PIN) in the form of ATP-binding cassette subfamily B/P-glycoprotein (ABCB/PGP) transporter, which can help plants to respond in a more dominant form to stress. We have demonstrated that the number and location of hydroxyl groups in the structure of flavonoids can determine their free radical scavenging ability and also elucidated the mechanism by which flavonoids exert free radical removal in cells. We also identified flavonoids as signaling molecules to promote rhizobial nodulation and colonization of arbuscular mycorrhizal fungi (AMF) to enhance plant-microbial symbiosis in defense to stresses. Given all this knowledge, we can foresee that the in-depth study of flavonoids will be an essential way to reveal plant tolerance and enhance plant stress resistance.

The multifaceted roles of R2R3 transcription factor HlMYB7 in the regulation of flavonoid and bitter acids biosynthesis, development and biotic stress tolerance in hop (Humulus lupulus L.)

Hop (Humulus lupulus) biosynthesizes the highly economically valuable secondary metabolites, which include flavonoids, bitter acids, polyphenols and essential oils. These compounds have important pharmacological properties and are widely implicated in the brewing industry owing to bittering flavor, floral aroma and preservative activity. Our previous studies documented that ternary MYB-bHLH-WD40 (MBW) and binary WRKY1-WD40 (WW) protein complexes transcriptionally regulate the accumulation of bitter acid (BA) and prenylflavonoids (PF). In the present study, we investigated the regulatory functions of the R2R3-MYB repressor HlMYB7 transcription factor, which contains a conserved N-terminal domain along with the repressive motif EAR, in regulating the PF- and BA-biosynthetic pathway and their accumulation in hop. Constitutive expression of HlMYB7 resulted in transcriptional repression of structural genes involved in the terminal steps of biosynthesis of PF and BA, as well as stunted growth, delayed flowering, and reduced tolerance to viroid infection in hop. Furthermore, yeast two-hybrid and transient reporter assays revealed that HlMYB7 targets both PF and BA pathway genes and suppresses MBW and WW protein complexes. Heterologous expression of HlMYB7 leads to down-regulation of structural genes of flavonoid pathway in Arabidopsis thaliana, including a decrease in anthocyanin content in Nicotiana tabacum. The combined results from functional and transcriptomic analyses highlight the important role of HlMYB7 in fine-tuning and balancing the accumulation of secondary metabolites at the transcriptional level, thus offer a plausible target for metabolic engineering in hop.

Transcriptome analysis of mulberry (Morus alba L.) leaves to identify differentially expressed genes associated with post-harvest shelf-life elongation

Present study deals with molecular expression patterns responsible for post-harvest shelf-life extension of mulberry leaves. Quantitative profiling showed retention of primary metabolite and accumulation of stress markers in NS7 and CO7 respectively. The leaf mRNA profiles was sequenced using the Illumina platform to identify DEGs. A total of 3413 DEGs were identified between the treatments. Annotation with Arabidopsis database has identified 1022 DEGs unigenes. STRING generated protein-protein interaction, identified 1013 DEGs nodes with p < 1.0e-16. KEGG classifier has identified genes and their participating biological processes. MCODE and BiNGO detected sub-networking and ontological enrichment, respectively at p ≤ 0.05. Genes associated with chloroplast architecture, photosynthesis, detoxifying ROS and RCS, and innate-immune response were significantly up-regulated, responsible for extending shelf-life in NS7. Loss of storage sucrose, enhanced activity of senescence-related hormones, accumulation of xenobiotics, and development of osmotic stress inside tissue system was the probable reason for tissue deterioration in CO7. qPCR validation of DEGs was in good agreement with RNA sequencing results, indicating the reliability of the sequencing platform. Present outcome provides a molecular insight regarding involvement of genes in self-life extension, which might help the sericulture industry to overcome their pre-existing problems related to landless farmers and larval feeding during monsoon.

Developmental regulation of lupulin gland-associated genes in aromatic and bitter hops (Humulus lupulus L.)

BACKGROUND

Hop (Humulus lupulus L.) bitter acids are valuable metabolites for the brewing industry. They are biosynthesized and accumulate in glandular trichomes of the female inflorescence (hop cone). The content of alpha bitter acids, such as humulones, in hop cones can differentiate aromatic from bitter hop cultivars. These contents are subject to genetic and environmental control but significantly correlate with the number and size of glandular trichomes (lupulin glands).

RESULTS

We evaluated the expression levels of 37 genes involved in bitter acid biosynthesis and morphological and developmental differentiation of glandular trichomes to identify key regulatory factors involved in bitter acid content differences. For bitter acid biosynthesis genes, upregulation of humulone synthase genes, which are important for the biosynthesis of alpha bitter acids in lupulin glands, could explain the higher accumulation of alpha bitter acids in bitter hops. Several transcription factors, including HlETC1, HlMYB61 and HlMYB5 from the MYB family, as well as HlGLABRA2, HlCYCB2-4, HlZFP8 and HlYABBY1, were also more highly expressed in the bitter hop cultivars; therefore, these factors may be important for the higher density of lupulin glands also seen in the bitter hop cultivars.

CONCLUSIONS

Gene expression analyses enabled us to investigate the differences between aromatic and bitter hops. This study confirmed that the bitter acid content in glandular trichomes (lupulin glands) is dependent on the last step of alpha bitter acid biosynthesis and glandular trichome density.

The Influence of Hop Latent Viroid (HLVd) Infection on Gene Expression and Secondary Metabolite Contents in Hop (Humulus lupulus L.) Glandular Trichomes

Viroids are small infectious pathogens, composed of a short single-stranded circular RNA. Hop (Humulus lupulus L.) plants are hosts to four viroids from the family Pospiviroidae. Hop latent viroid (HLVd) is spread worldwide in all hop-growing regions without any visible symptoms on infected hop plants. In this study, we evaluated the influence of HLVd infection on the content and the composition of secondary metabolites in maturated hop cones, together with gene expression analyses of involved biosynthesis and regulation genes for Saaz, Sládek, Premiant and Agnus cultivars. We confirmed that the contents of alpha bitter acids were significantly reduced in the range from 8.8% to 34% by viroid infection. New, we found that viroid infection significantly reduced the contents of xanthohumol in the range from 3.9% to 23.5%. In essential oils of Saaz cultivar, the contents of monoterpenes, terpene epoxides and terpene alcohols were increased, but the contents of sesquiterpenes and terpene ketones were decreased. Secondary metabolites changes were supported by gene expression analyses, except essential oils. Last-step biosynthesis enzyme genes, namely humulone synthase 1 (HS1) and 2 (HS2) for alpha bitter acids and O-methytransferase 1 (OMT1) for xanthohumol, were down-regulated by viroid infection. We found that the expression of ribosomal protein L5 (RPL5) RPL5 and the splicing of transcription factor IIIA-7ZF were affected by viroid infection and a disbalance in proteosynthesis can influence transcriptions of biosynthesis and regulatory genes involved in of secondary metabolites biosynthesis. We suppose that RPL5/TFIIIA-7ZF regulatory cascade can be involved in HLVd replication as for other viroids of the family Pospiviroidae.

Gene expression for secondary metabolite biosynthesis in hop (Humulus lupulus L.) leaf lupulin glands exposed to heat and low-water stress

Hops are valued for their secondary metabolites, including bitter acids, flavonoids, oils, and polyphenols, that impart flavor in beer. Previous studies have shown that hop yield and bitter acid content decline with increased temperatures and low-water stress. We looked at physiological traits and differential gene expression in leaf, stem, and root tissue from hop (Humulus lupulus) cv. USDA Cascade in plants exposed to high temperature stress, low-water stress, and a compound treatment of both high temperature and low-water stress for six weeks. The stress conditions imposed in these experiments caused substantial changes to the transcriptome, with significant reductions in the expression of numerous genes involved in secondary metabolite biosynthesis. Of the genes involved in bitter acid production, the critical gene valerophenone synthase (VPS) experienced significant reductions in expression levels across stress treatments, suggesting stress-induced lability in this gene and/or its regulatory elements may be at least partially responsible for previously reported declines in bitter acid content. We also identified a number of transcripts with homology to genes shown to affect abiotic stress tolerance in other plants that may be useful as markers for breeding improved abiotic stress tolerance in hop. Lastly, we provide the first transcriptome from hop root tissue.

Establishment of CRISPR/Cas9 mediated targeted mutagenesis in hop (Humulus lupulus)

The CRISPR/Cas9-based targeted genome editing has emerged as a versatile technique, widely employed in plant genome engineering, both to decipher gene function and as an alternative to classical breeding technique for traits improvement in plants. However, to date, no such platform has been developed for hop (Humulus lupulus L.), which is an economically important crop producing valuable secondary metabolites utilized in the brewing and pharmaceutical industries. Here, we present the first report on the successful establishment of efficient CRISPR/Cas9-based genome editing using the visible endogenous marker gene phytoene desaturase (PDS) involved in carotenoid biosynthesis to demonstrate successful genome editing in hop. Agrobacterium tumefaciens-mediated transformation of in vitro generated internodal explants was used for the stable integration of constructs expressing plant codon-optimized Cas9 and a pair of co-expressed guide RNAs to target the distinct genomic sites of the PDS gene of hop. Analysis of RNA-guided genome-editing events, including mutant lines screening and homozygosity assessment using the T7 endonuclease assay showed that 33.3% of transformed plants were successfully edited at the target site, displaying albino and mosaic regenerants. Intriguingly, the detected mutations were ranges of deletions (16 bp to 39 bp) which led to disruption of the exon-intron boundary, few base substitutions, and a 1 bp insertion at 3 bp upstream of the PAM region of the target site. The decrease in chlorophyll a/b, and carotenoid content in the mutant lines further confirmed the functional disruption of the HlPDS gene. Taken together, our results demonstrate that the CRISPR/Cas9 system can precisely edit the targeted genome sequences, which may revolutionize our way to overcome some of the obstacles that have plagued the traits improvement in hop.

Dissection of Dynamic Transcriptome Landscape of Leaf, Bract, and Lupulin Gland in Hop (Humulus lupulus L.)

The hop plant (Humulus lupulus L.) produces several valuable secondary metabolites, such as prenylflavonoid, bitter acids, and essential oils. These compounds are biosynthesized in glandular trichomes (lupulin glands) endowed with pharmacological properties and widely implicated in the beer brewing industry. The present study is an attempt to generate exhaustive information of transcriptome dynamics and gene regulatory mechanisms involved in biosynthesis and regulation of these compounds, developmental changes including trichome development at three development stages, namely leaf, bract, and mature lupulin glands. Using high-throughput RNA-Seq technology, a total of 61.13, 50.01, and 20.18 Mb clean reads in the leaf, bract, and lupulin gland libraries, respectively, were obtained and assembled into 43,550 unigenes. The putative functions were assigned to 30,996 transcripts (71.17%) based on basic local alignment search tool similarity searches against public sequence databases, including GO, KEGG, NR, and COG families, which indicated that genes are principally involved in fundamental cellular and molecular functions, and biosynthesis of secondary metabolites. The expression levels of all unigenes were analyzed in leaf, bract, and lupulin glands tissues of hop. The expression profile of transcript encoding enzymes of BCAA metabolism, MEP, and shikimate pathway was most up-regulated in lupulin glands compared with leaves and bracts. Similarly, the expression levels of the transcription factors and structural genes that directly encode enzymes involved in xanthohumol, bitter acids, and terpenoids biosynthesis pathway were found to be significantly enhanced in lupulin glands, suggesting that production of these metabolites increases after the leaf development. In addition, numerous genes involved in primary metabolism, lipid metabolism, photosynthesis, generation of precursor metabolites/energy, protein modification, transporter activity, and cell wall component biogenesis were differentially regulated in three developmental stages, suggesting their involvement in the dynamics of the lupulin gland development. The identification of differentially regulated trichome-related genes provided a new foundation for molecular research on trichome development and differentiation in hop. In conclusion, the reported results provide directions for future functional genomics studies for genetic engineering or molecular breeding for augmentation of secondary metabolite content in hop.

Evaluation of Disease Severity and Global Transcriptome Response Induced by Citrus bark cracking viroid, Hop latent viroid, and Their Co-Infection in Hop (Humulus lupulus L.)

Viroids are small non-capsidated, single-stranded, covalently-closed circular noncoding RNA replicons of 239-401 nucleotides that exploit host factors for their replication, and some cause disease in several economically important crop plants, while others appear to be benign. The proposed mechanisms of viroid pathogenesis include direct interaction of the genomic viroid RNA with host factors and post-transcriptional or transcriptional gene silencing via viroid-derived small RNAs (vd-sRNAs) generated by the host defensive machinery. Humulus lupulus (hop) plants are hosts to several viroids among which Hop latent viroid (HLVd) and Citrus bark cracking viroid (CBCVd) are attractive model systems for the study of viroid-host interactions due to the symptomless infection of the former and severe symptoms induced by the latter in this indicator host. To better understand their interactions with hop plant, a comparative transcriptomic analysis based on RNA sequencing (RNA-seq) was performed to reveal the transcriptional alterations induced as a result of single HLVd and CBCVd infection in hop. Additionally, the effect of HLVd on the aggressiveness of CBCVd that underlies severe stunting in hop in a mixed infection was studied by transcriptomic analysis. Our analysis revealed that CBCVd infection resulted in dynamic changes in the activity of genes as compared to single HLVd infection and their mixed infection. The differentially expressed genes that are involved in defense, phytohormone signaling, photosynthesis and chloroplasts, RNA regulation, processing and binding; protein metabolism and modification; and other mechanisms were more modulated in the CBCVd infection of hop. Nevertheless, Gene Ontology (GO) classification and pathway enrichment analysis showed that the expression of genes involved in the proteolysis mechanism is more active in a mixed infection as compared to a single one, suggesting co-infecting viroids may result in interference with host factors more prominently. Collectively, our results provide a deep transcriptome of hop and insight into complex single HLVd, CBCVd, and their coinfection in hop-plant interactions.

Genome-Wide Transcriptomic Analysis Reveals Insights into the Response to Citrus bark cracking viroid (CBCVd) in Hop (Humulus lupulus L.)

Viroids are smallest known pathogen that consist of non-capsidated, single-stranded non-coding RNA replicons and they exploits host factors for their replication and propagation. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) is a serious threat, which spreads rapidly within hop gardens. In this study, we employed comprehensive transcriptome analyses to dissect host-viroid interactions and identify gene expression changes that are associated with disease development in hop. Our analysis revealed that CBCVd-infection resulted in the massive modulation of activity of over 2000 genes. Expression of genes associated with plant immune responses (protein kinase and mitogen-activated protein kinase), hypersensitive responses, phytohormone signaling pathways, photosynthesis, pigment metabolism, protein metabolism, sugar metabolism, and modification, and others were altered, which could be attributed to systemic symptom development upon CBCVd-infection in hop. In addition, genes encoding RNA-dependent RNA polymerase, pathogenesis-related protein, chitinase, as well as those related to basal defense responses were up-regulated. The expression levels of several genes identified from RNA sequencing analysis were confirmed by qRT-PCR. Our systematic comprehensive CBCVd-responsive transcriptome analysis provides a better understanding and insights into complex viroid-hop plant interaction. This information will assist further in the development of future measures for the prevention of CBCVd spread in hop fields.