Plants have evolved elegant regulatory programs that can rapidly elicit the expression of defense-related genes in response to pests and pathogens, such as viroids. These defense-related genes include small RNAs, long noncoding RNAs, protein-coding mRNAs, etc. Transcriptomic analyses provide a powerful toolset to uncover those genes that have altered expression profiles in viroid-infected plants. This chapter introduces the basic concept and a brief history of transcriptomic analysis as well as its application in understanding plant-viroid interactions. This chapter intends to focus on introducing the background and methodology of transcriptomic analyses instead of providing a comprehensive literature review. While the chapter focuses on the description of the most popular RNA deep sequencing analyses, readers are referred to other resources in the Further Reading section to gain a better grasp on small RNA deep sequencing, degradome RNA deep sequencing, etc.

More than 50 years ago, viroids were firstly described as the smallest RNA molecules capable to infect certain plants and to autonomously self-replicate in host plants. Viroids are covalently closed circular single-stranded RNAs that are non-coding and depend for most of their infection cycle on host proteins. Today, viroids are subdivided into the two families Avsunviroidae and Pospiviroidae. Members of Avsunviroidae replicate in the chloroplast and have a highly bifurcated structure including hammerhead ribozymes, which cleave oligomeric replication intermediates into monomers and ligate them to mature circles. Members of Pospiviroidae accumulate in the nucleus, have a rod-like structure and depend on host proteins for cleavage and ligation. We will describe our present knowledge on sequence and structural elements of viroids in connection to their replication and trafficking.

Viroids are small, non-coding, pathogenic RNAs with the ability to disturb plant developmental processes. This dysregulation redirects the morphogenesis of plant organs, significantly impairing their functionality. Citrus bark cracking viroid (CBCVd) causes detrimental developmental distortions in infected hops (Humulus lupulus) and causes significant economic losses. CBCVd can infect cells and tissues of the model plant tobacco (Nicotiana tabacum), provided it is delivered via transgenesis. The levels of CBCVd in tobacco were enhanced in plant hybrids expressing CBCVd cDNAs and either the tobacco or hop variant of TFIIIA-7ZF, a viroid-mediated splicing derivative of transcription factor IIIA, which is important for viroid replication by DNA-dependent RNA polymerase II. The TFIIIA-7ZF variants can change the tobacco morphogenesis if expressed in leaves and shoots. In addition to the splitting of shoots, the “pathomorphogenic” network in hybrid plants expressing CBCVd and HlTFIIIA-7ZF induced leaf fusions and malformations. Moreover, CBCVd can dramatically change another morphogenesis into teratomic and petal-like tissues if propagated above some limit in young transgenic tobacco microspores and anthers. By comparative RNA profiling of transgenic tobacco shoots bearing TFIIIA-7ZFs and CBCVd-transformed/infected anthers, we found a differential expression of many genes at p < 0.05. As the main common factor showing the differential up-regulation in shoot and anther tissues, a LITTLE ZIPPER 2-like transcription factor was found. We propose that this factor, which can interact as a competitive inhibitor of the also dysregulated homeobox-leucin zipper family protein (HD-ZIPIII) in apical meristem, is essential for a network responsible for some morphological changes and modifications of plant degradome within shoot meristem regulation and secondary xylem differentiation.

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.

Viroids are small, non-coding, pathogenic RNAs with a significant ability of adaptation to several basic cellular processes in plants. TFIIIA-7ZF, a splicing variant of transcription factor IIIA, is involved in replication of nuclear-replicating viroids by DNA-dependent polymerase II. We overexpressed NbTFIIIA-7ZF from Nicotiana benthamiana in tobacco (Nicotiana tabacum) where it caused morphological and physiological deviations like plant stunting, splitting of leaf petioles, pistils or apexes, irregular branching of shoots, formation of double-blade leaves, deformation of main stems, and modification of glandular trichomes. Plant aging and senescence was dramatically delayed in transgenic lines. Factors potentially involved in viroid degradation and elimination in pollen were transiently depressed in transgenic leaves. This depressed "degradome" in young plants involved NtTudor S-like nuclease, dicers, argonoute 5, and pollen extracellular nuclease I showing expression in tobacco anthers and leaves. Analysis of the "degradome" in tobacco leaves transformed with either of two hop viroids confirmed modifications of the "degradome" and TFIIIA expression. Thus, the regulatory network connected to TFIIIA-7ZF could be involved in plant pathogenesis as well as in viroid adaptation to avoid its degradation. These results support the hypothesis on a significant impact of limited TFIIIA-7ZF on viroid elimination in pollen.

Long non-coding RNAs (lncRNAs) are a highly heterogeneous class of non-protein-encoding transcripts that play an essential regulatory role in diverse biological processes, including stress responses. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) poses a major threat to the production of Humulus lupulus (hop) plants. In this study, we systematically investigate the characteristics of the lncRNAs in hop and their role in CBCVd-infection using RNA-sequencing data. Following a stringent filtration criterion, a total of 3598 putative lncRNAs were identified with a high degree of certainty, of which 19% (684) of the lncRNAs were significantly differentially expressed (DE) in CBCVd-infected hop, which were predicted to be mainly involved in plantpathogen interactions, kinase cascades, secondary metabolism and phytohormone signal transduction. Besides, several lncRNAs and CBCVd-responsive lncRNAs were identified as the precursor of microRNAs and predicted as endogenous target mimics (eTMs) for hop microRNAs involved in CBCVd-infection.

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 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.

Members of the family Pospiviroidae have single-stranded circular RNA genomes that adopt a rod-like or a quasi-rod-like conformation. These genomes contain a central conserved region that is involved in replication in the nucleus through an asymmetric RNA-RNA rolling-circle mechanism. Members of the family Pospiviroidae lack the hammerhead ribozymes that are typical of viroids classified in the family Avsunviroidae. The family Pospiviroidae includes the genera Apscaviroid, Cocadviroid, Coleviroid, Hostuviroid and Pospiviroid, with >25 species. This is a summary of the ICTV Report on the family Pospiviroidae, which is available at ictv.global/report/pospiviroidae.

Viroids are small, non-coding, parasitic RNAs that promote developmental distortions in sensitive plants. We analyzed pollen of Nicotiana benthamiana after infection and/or ectopic transformation with cDNAs of citrus bark cracking viroid (CBCVd), apple fruit crinkle viroid (AFCVd) and potato spindle tuber viroid (PSTVd) variant AS1. These viroids were seed non-transmissible in N. benthamiana. All viroids propagated to high levels in immature anthers similar to leaves, while their levels were drastically reduced by approximately 3.6 × 103, 800 and 59 times in mature pollen of CBCVd, AFCVd and PSTVd infected N. benthamiana, respectively, in comparison to leaves. These results suggest similar elimination processes during male gametophyte development as in the Nicotiana tabacum we presented in our previous study. Mature pollen of N. benthamiana showed no apparent defects in infected plants although all three viroids induced strong pathological symptoms on leaves. While Nicotiana species have naturally bicellular mature pollen, we noted a rare occurrence of mature pollen with three nuclei in CBCVd-infected N. benthamiana. Changes in the expression of ribosomal marker proteins in AFCVd-infected pollen were detected, suggesting some changes in pollen metabolism. N. benthamiana transformed with 35S-driven viroid cDNAs showed strong symptoms including defects in pollen development. A large number of aborted pollen (34% and 62%) and a slight increase of young pollen grains (8% and 15%) were found in mature pollen of AFCVd and CBCVd transformants, respectively, in comparison to control plants (3.9% aborted pollen and 0.3% young pollen). Moreover, pollen grains with malformed nuclei or trinuclear pollen were found in CBCVd-transformed plants. Our results suggest that "forcing" overexpression of seed non-transmissible viroid led to strong pollen pathogenesis. Viroid adaptation to pollen metabolism can be assumed as an important factor for viroid transmissibility through pollen and seeds.

Some viroids—single-stranded, non-coding, circular RNA parasites of plants—are not transmissible through pollen to seeds and to next generation. We analyzed the cause for the elimination of apple fruit crinkle viroid (AFCVd) and citrus bark cracking viroid (CBCVd) from male gametophyte cells of Nicotiana tabacum by RNA deep sequencing and molecular methods using infected and transformed tobacco pollen tissues at different developmental stages. AFCVd was not transferable from pollen to seeds in reciprocal pollinations, due to a complete viroid eradication during the last steps of pollen development and fertilization. In pollen, the viroid replication pathway proceeds with detectable replication intermediates, but is dramatically depressed in comparison to leaves. Specific and unspecific viroid degradation with some preference for (−) chains occurred in pollen, as detected by analysis of viroid-derived small RNAs, by quantification of viroid levels and by detection of viroid degradation products forming “comets” on Northern blots. The decrease of viroid levels during pollen development correlated with mRNA accumulation of several RNA-degrading factors, such as AGO5 nuclease, DICER-like and TUDOR S-like nuclease. In addition, the functional status of pollen, as a tissue with high ribosome content, could play a role during suppression of AFCVd replication involving transcription factors IIIA and ribosomal protein L5.

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.

Tobacco (Nicotiana tabacum) pollen is a well-suited model for studying many fundamental biological processes owing to its well-defined and distinct development stages. It is also one of the major agents involved in the transmission of infectious viroids, which is the primary mechanism of viroid pathogenicity in plants. However, some viroids are non-transmissible and may be possibly degraded or eliminated during the gradual process of pollen development maturation. The molecular details behind the response of developing pollen against the apple fruit crinkle viroid (AFCVd) infection and viroid eradication is largely unknown. In this study, we performed an integrative analysis of the transcriptome and proteome profiles to disentangle the molecular cascade of events governing the three pollen development stages: early bicellular pollen (stage 3, S3), late bicellular pollen (stage 5, S5), and 6 h-pollen tube (PT6). The integrated analysis delivered the molecular portraits of the developing pollen against AFCVd infection, including mechanistic insights into the viroid eradication during the last steps of pollen development. The isobaric tags for label-free relative quantification (iTRAQ) with digital gene expression (DGE) experiments led us to reliably identify subsets of 5321, 5286, and 6923 proteins and 64,033, 60,597, and 46,640 expressed genes in S3, S5, and PT6, respectively. In these subsets, 2234, 2108 proteins and 9207 and 14,065 mRNAs were differentially expressed in pairwise comparisons of three stages S5 vs. S3 and PT6 vs. S5 of control pollen in tobacco. Correlation analysis between the abundance of differentially expressed mRNAs (DEGs) and differentially expressed proteins (DEPs) in pairwise comparisons of three stages of pollen revealed numerous discordant changes in mRNA/protein pairs. Only a modest correlation was observed, indicative of divergent transcription, and its regulation and importance of post-transcriptional events in the determination of the fate of early and late pollen development in tobacco. The functional and enrichment analysis of correlated DEGs/DEPs revealed the activation in pathways involved in carbohydrate metabolism, amino acid metabolism, lipid metabolism, and cofactor as well as vitamin metabolism, which points to the importance of these metabolic pathways in pollen development. Furthermore, the detailed picture of AFCVd-infected correlated DEGs/DEPs was obtained in pairwise comparisons of three stages of infected pollen. The AFCVd infection caused the modulation of several genes involved in protein degradation, nuclear transport, phytohormone signaling, defense response, and phosphorylation. Intriguingly, we also identified several factors including, DNA-dependent RNA-polymerase, ribosomal protein, Argonaute (AGO) proteins, nucleotide binding proteins, and RNA exonucleases, which may plausibly involve in viroid stabilization and eradication during the last steps of pollen development. The present study provides essential insights into the transcriptional and translational dynamics of tobacco pollen, which further strengthens our understanding of plant-viroid interactions and support for future mechanistic studies directed at delineating the functional role of candidate factors involved in viroid elimination.

The mediator (MED) represents a large, conserved, multi-subunit protein complex that regulates gene expression through interactions with RNA polymerase II and enhancer-bound transcription factors. Expanding research accomplishments suggest the predominant role of plant MED subunits in the regulation of various physiological and developmental processes, including the biotic stress response against bacterial and fungal pathogens. However, the involvement of MED subunits in virus/viroid pathogenesis remains elusive. In this study, we investigated for the first time the gene expression modulation of selected MED subunits in response to five viroid species (Apple fruit crinkle viroid (AFCVd), Citrus bark cracking viroid (CBCVd), Hop latent viroid (HLVd), Hop stunt viroid (HSVd), and Potato spindle tuber viroid (PSTVd)) in two model plant species (Nicotiana tabacum and N. benthamiana) and a commercially important hop (Humulus lupulus) cultivar. Our results showed a differential expression pattern of MED subunits in response to a viroid infection. The individual plant MED subunits displayed a differential and tailored expression pattern in response to different viroid species, suggesting that the MED expression is viroid- and plant species-dependent. The explicit evidence obtained from our results warrants further investigation into the association of the MED subunit with symptom development. Together, we provide a comprehensive portrait of MED subunit expression in response to viroid infection and a plausible involvement of MED subunits in fine-tuning transcriptional reprogramming in response to viroid infection, suggesting them as a potential candidate for rewiring the defense response network in plants against pathogens.

Hop ( Humulus lupulus L.) plants are hosts to several viroids; some of them can be highly aggressive, and their infection can manifest as complete plant dieback. Since molecular detection of multiple viroids can be time consuming and cost inefficient, a reliable one-step multiplex RT-PCR (mRT-PCR) was developed to detect simultaneously all four viroids infecting hops: Hop latent viroid (HLVd), Hop stunt viroid (HSVd), Apple fruit crinkle viroid (AFCVd) and Citrus bark cracking viroid (CBCVd). Several primer pairs were tested on different viroid variants from hops, citruses and grapevines, and from among them, specific primer pairs for detection of hop viroids were selected and confirmed in a single-tube assay. To improve mRT-PCR reliability and validate its effectiveness, nad5 and DRH1 genes were included as an internal control. The specificities of single and mRT-PCR assays for all four viroids were comparable. The sensitivity of mRT-PCR was compared with that of dot-blot hybridization and single RT-PCR assay on biolistically infected hop plants. The results show mRT-PCR to be more sensitive than the dot-blot and slightly less sensitive than the single RT-PCR assay. Furthermore, mRT-PCR was validated using field samples and a group of 135 hop plants, which are used in the certification scheme for planting material propagation, and the method proved to be robust, rapid and simple. Additionally, this approach can be applicable to similar methods of systematic surveys of emerging diseases and epidemiological studies.

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

Transcription factors (TFs) play a major role in controlling gene expression by intricately regulating diverse biological processes such as growth and development, the response to external stimuli and the activation of defense responses. The systematic identification and classification of TF genes are essential to gain insight into their evolutionary history, biological roles, and regulatory networks. In this study, we performed a global mining and characterization of hop TFs and their involvement in Citrus bark cracking viroid CBCVd infection by employing a digital gene expression analysis. Our systematic analysis resulted in the identification of a total of 3,818 putative hop TFs that were classified into 99 families based on their conserved domains. A phylogenetic analysis classified the hop TFs into several subgroups based on a phylogenetic comparison with reference TF proteins from Arabidopsis thaliana providing glimpses of their evolutionary history. Members of the same subfamily and subgroup shared conserved motif compositions. The putative functions of the CBCVd-responsive hop TFs were predicted using their orthologous counterparts in A. thaliana. The analysis of the expression profiling of the CBCVd-responsive hop TFs revealed a massive differential modulation, and the expression of the selected TFs was validated using qRT-PCR. Together, the comprehensive integrated analysis in this study provides better insights into the TF regulatory networks associated with CBCVd infections in the hop, and also offers candidate TF genes for improving the resistance in hop against viroids.

The regulation of gene expression involves a multifarious regulatory system. Each gene contains a unique combination of cis-acting regulatory sequence elements in the 5' regulatory region that determines its temporal and spatial expression. Recent release of Humulus lupulus (hop), enabled us to predict the possible cis-acting regulatory elements by by scanning 2 kbp of 5' regulatory regions of genes important for lupulin biosynthesis, using Plant CARE, PLACE and Genomatix Matinspector professional databases. The result reveals the probable cis-acting regulatory elements that possibly are involved in the expression genes important for lupulin biosynthesis. In addition, new generation sequencing of small RNAs in hop leads to significant progress in identification of microRNAs that act as predominat regulators of transcription factors (TFs). This study is the first endeavor on mining of miRNAs and analyzing their involvement in secondary metabolism in hop.

Cannabis sativa, an annual herbaceous plant, produce wide variety of secondary metabolites among which delta-9-tetrahydrocannabinol (THC) is the most important one. The dissection of biosynthetic pathway(s) of this compound and its regulation by transcription factors (TFs) is an important prerequisite for efficient biotechnological manipulation of its secondary metabolome. A polyketide synthase (PKS) of C. sativa catalyzes the first step of cannabinoid biosynthesis, leading to the biosynthesis of olivetolic acid. Cloning and analysis of PKS promoter based on online PLACE, PlantCARE, and Genomatix Matinspector professional databases, indicated that PKS promoter consisted of cis-elements such as TATA-box, CAAT-box, W-box, Myb-box, E-box, and P-box. Plant expression vector PKS::GUS was constructed in such a way that the ATG of the PKS gene was in the frame with the β-glucuronidase (GUS) coding region. Using a combinatorial transient GUS expression system in Nicotiana benthamania leaves, it was shown that heterologous TFs such as HlWRKY1, HlMYB3, HlWDR1 and HlbZIP1 from Humulus lupulus significantly activated PKS promoter. Moreover, Tombusvirus p19 core protein, which is known for silencing suppressor functions, acted in our combinatorial transient expression system as an enhancer of PKS promoter activity along with hop TFs. Our analyses suggested the involvement of the hop derived TFs (HlWRKY1, HlMYB3, HlWDR1 and HlbZIP1A) and p19 in the activation of PKS gene promoter, which could be used for the genetic manipulation of C. sativa to enhance the cannabinoid production.

Table S5. The unigenes differentially expressed between the control and WW-transgenic plants of hop. (XLS 296 kb)