Skip to main content

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

Abstract

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.

Peer Review reports

Background

Hop (Humulus lupulus L.) is a diploid, dioecious, perennial climbing plant belonging to the Cannabaceae family. Female plants are cultivated for the commercial production of inflorescences (cones), which are mainly used in the brewing industry but also commonly used in the production of pharmaceuticals and cosmetics [1]. In cones, flavour-active secondary metabolites are biosynthesized and accumulate in glandular trichomes and lupulin glands. There are three key classes of natural products: alpha and beta bitter acids (humulone and lupulone derivatives, respectively), prenylated flavonoids (primarily xanthohumol and desmethylxanthohumol) and essential oils composed mainly of myrcene, α-humulene and β-caryophyllene [2].

The contents of alpha and beta bitter acids in hop cones (5 to 30% of the dry weight) are major chemical characteristics and economical traits of different cultivars during hop production. Their biosynthesis requires precursor sources that directly originate from primary sucrose metabolism and amino acid precursors such as leucine, isoleucine, valine and phenylalanine [3]. Phenylalanine is biosynthesized by the shikimate pathway and others by the branched-chain amino acid (BCAA) pathway in chloroplasts. The final step in the BCAA pathway is catalysed by the enzyme branched-chain amino transferase 2 (HlBCAT2). For bitter acid biosynthesis, BCAAs are then degraded in the mitochondria by branched-chain amino transferase 1 (HlBCAT1) and branched-chain keto-acid dehydrogenase (HlBCKDH) to isovaleric, 2-methylbutyric and isobutyric acids, respectively [4]. Coenzyme A is ligated to these precursors by two branched short-chain fatty acid CoA ligases (HlCCL2 and HlCCL4) in the cytosol [5]. Valerophenone synthase (HlVPS) synthesizes phloroisovalerophenone (PIVP), phloro-2-methylbutyrophenone (PMBP) and phloroisobutyrophenone (PIBP) from precursors and three malonyl-CoAs from phenylalanine [6].

Isoprenoids are a large and highly diverse group synthesized from common prenyl diphosphate precursors. Plants synthesize isopentenyl diphosphate (IPP) via mevalonate and methylerythritol phosphate (MEP) pathways [7]. The MEP pathway is more efficient in hop lupulin glands according to transcriptome [4] and proteome [3] studies. Terpenic compounds are the main components of hop essential oils.

Bitter acids and prenylflavonoids are additional specific compounds of prenylation in lupulin glands. Both are prenylated by a functional heterocomplex of two aromatic prenyltransferases, HlPT1L and HlPT2, in chloroplasts [8,9,10]. HlPT1L catalyses the first prenylation step, and HlPT2 catalyses the next two prenylation steps. The first step produces alpha bitter acid precursors, such as deoxyhumulone, deoxyadhumulone and deoxycohumulone, and the second step produces beta bitter acids, such as lupulone, adlupulone and colupulone. Alpha bitter acids (humulone, adhumulone and cohumulone) are finally produced by humulone synthase (HlHS1 and HlHS2) via the oxygenation of deoxy-precursors [3]. The content of individual analogues of alpha and beta bitter acids is determined by the ratios of precursors and is genetically imprinted in the genomes of individual genotypes; for example, cohumulone (from 12 to 55% of alpha acids) and colupulone (from 30 to 80% of beta acids) contents are very stable and heritable traits [2].

Since the expression of biosynthetic genes is the highest in lupulin glands, there is a tissue-specific regulatory network of transcription factors mainly from the bHLH, MYB, WDR and WRKY families [1]. Promoters of bitter acid and prenylflavonoid biosynthesis pathway genes are regulated by heterotrimeric ternary MBW (HlMyb3/HlbHLH2/HlWDR1 or HlMYB2/HlbHLH2/HlWDR1) complexes [11], binary (HlbHLH2/HlWDR1 or HlWRKY1/HlWDR1) transcription activation complexes [11, 12] or individual transcription factors (HlWRKY1 or HlMyb8) [13, 14]. This regulation allows hop plants to react to developmental changes, environmental stresses and weather conditions, which influence bitter acid contents.

The bitter acid content in hop cones varies across hop cultivars due to their genetic background. Hop cultivars are divided into aromatic hops, with low alpha bitter acid contents in dry cones ranging from 0.5 to 8%, and bitter hops, with alpha bitter acid contents in dry cones over 9 to 23% [15]. It was found that the content of alpha bitter acids very strongly and significantly correlated with the number and size of lupulin glands in hop cones [2, 16, 17]. Therefore, bitter acid contents can be increased either by increasing their production inside secretory cells or by increasing the density of glands on individual cones. The development of multicellular glandular trichomes includes the enlargement of single epidermal cells, followed by several polarized and localized cell divisions and remodelling to generate branched or unbranched structures perpendicular to the epidermal surface. The lupulin glands develop on the abaxial surface of the leaf, bract or bracteole primordia and continue to form until expansion stops [1]. There are already molecular data for genes that play a specific role in glandular trichome development, especially transcription factors, cell cycle regulators, and receptors involved in phytohormone-induced signalling cascades [18].

The model plant A. thaliana is the most utilized for the developmental study of trichomes, which are single-celled and non-glandular [19]. The current model includes transcription factors that function as positive or negative regulators of trichome formation, upstream regulators of these two groups of regulators, and downstream components [20]. The positive regulators belong to three protein classes and include a WD40-repeat protein TRANSPARENT TESTA GLABRA1 (TTG1), three R2R3 MYB-related transcription factors GLABRA1 (GL1, MYB23, MYB5) and four basic helix-loop-helix (bHLH)-like transcription factors GLABRA3 (GL3), ENHANCER OF GLABRA3 (EGL3), TRANSPARENT TESTA (TT8), and MYC-1. They act partially redundantly and form a multimeric activator complex, also known as the MYB-bHLH-WD40 (MBW) complex, which binds the promoter of GLABRA2 (GL2) and several homeodomain (HD-ZIP class IV) transcription factors [19]. GL2 encodes a homeodomain protein required for subsequent phases of trichome morphogenesis, such as endoreduplication, branching, and maturation of the cell wall. The MBW activator complex is usually negatively regulated by single-repeat R3 MYBs, such as CAPRICE (CPC), TRIPTYCHON (TRY), ENHANCER OF TRY and CPC1 (ETC1), ETC2, TRICHOMELESS1 (TCL1) and TCL2/CPL4 [20]. Additionally, several C2H2 zinc finger protein transcription factor genes, including GLABROUS INFLORESCENCE STEMS (GIS, GIS2) and ZINC FINGER PROTEIN 8 and 5 (ZFP8, ZFP5), controlling GL1 and GL3 of the core MBW complex, have also been identified [19, 20].

In contrast, knowledge of the development of glandular trichomes is very limited, and some data have shows that glandular trichome pathways are not as conserved as the nonglandular trichome pathways known in Arabidopsis. What we know about glandular trichome development is mostly from work on Artemisia, tomato and tobacco [19].

In tomato, several genes have been reported that are required for the proper development and function of different types of glandular trichomes. The initiation and development of type I trichomes was controlled by a woolly (Wo) gene encoding a class IV homeodomain-leucine zipper protein homologue of Arabidopsis GL2 and a B-type cyclin gene, SlCycB2 (possibly regulated by Wo). A hairy phenotype was caused by the overproduction of mutant alleles of Wo in type I trichomes. However, the suppression of Wo or SlCycB2 expression by RNA interference decreased trichome density in tomato [21, 22]. Addititonally, three transcription factors, Woolly, Hair, and Mixta-like, were found to be necessary for the development of type I glandular trichomes, while Woolly, Hair, and Myc1 were found to be required for type VI glandular trichome development [19, 23]. SlMYC1 is also involved in the transcriptional regulation of several terpene synthases [24]. Trichome formation is also inhibited by a tomato TRY orthologue. Therefore, some of the MYB-bHLH-WD40 complex elements might be conserved in tomato together with HD-ZIP and C2H2 factors [19].

A recent RNA-sequencing analysis showed that orthologues of Arabidopsis genes involved in trichome formation via the MYB-bHLH-WD40 regulatory mechanism were expressed in tobacco trichomes. Nevertheless, the overexpression of a Wo transgene had no influence on their transcription and trichome induction. In contrast, Wo and SlCycB2 homologues involved in trichome formation in asterids were significantly upregulated in Wo transgenic tobacco plants [18, 25]. The expression of a C2H2-zinc-finger transcription factor, NbGIS, or the tomato Wo gene both led to more glandular trichomes being produced [26].

In Artemisia annua, AaMYB1 overexpression induces a higher density of trichomes [20]; AaHD1, a homeodomain-leucine zipper TF, is necessary for JA-mediated glandular trichome initiation [18]; and the AaHD8/AaMIXTA1 complex regulates trichome initiation and cuticle development [27]. Increased expression of AmMIXTA from Antirrhinum majus also leads to longer glandular trichomes in tobacco leaves [20]. Additionally, some other positive regulators of glandular trichome formation have been found, such as GoPGF, a bHLH TF from cotton, and three HD-Zip TFs, CsGL1, CsGL2 and CsGL3, in cucumber [1, 18]. Collectively, results from these species indicate that similar types of transcription factors control glandular trichome initiation and development [19].

In this study, we evaluated the expression levels of 37 genes involved in bitter acid biosynthesis and the morphological and developmental differentiation of lupulin glands in aromatic and bitter hops to identify key regulatory factors for bitter acid content differences.

Results

Genes involved in bitter acid biosynthesis

Glandular trichomes and lupulin glands are generally dedicated to the synthesis of specialized hop secondary metabolites. Bitter acids are major chemical compounds with economic impacts on hop production. Therefore, we evaluated expression differences of known bitter acid biosynthetic pathway enzyme genes between aromatic and bitter hop cultivars. We found that the gene expression of the HlBCAT1, HlVPS, HlPT1L and HlPT2 genes (Fig. 1A-D) was upregulated in lupulin glands (from 170 to 1000 times compared to leaves) and downregulated in leaves and bracts. There were no significant differences in these genes between aromatic and bitter hop cultivars, except for HlBCAT1 and HlPT1L. Two monooxygenase 2 genes, HlHS1 and HlHS2, were 5.8 times more highly expressed in the lupulin glands of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 1E and F). This trend was also evident in flower and young cone samples.

Fig. 1
figure 1

The relative expression of bitter acid biosynthesis genes A: HlBCAT1, B: HlVPS, C: HlPT1L, D: HlPT2, E: HlHS1 and F: HlHS2 (see Table 1) relative to reference genes in different tissues of aromatic and bitter hop cultivars. Tissues: A – apex, B – leaf, C – bract without lupulin, D – flower, E – young cone, and F – lupulin gland. Significant t-test group differences at the following probability levels: * - P < 0,1, ** - P < 0,05, and *** - P < 0,01. straight line - median, dashed line – average, and box – 95% percentile ± standard deviation

Transcription factors involved in bitter acid biosynthesis pathways

Relative expression differences of biosynthesis pathway genes can be caused by the regulatory networks of involved transcription factors. All of the studied transcription factors, which had previously been suggested to regulate bitter acid biosynthesis, were highly expressed in lupulin glands (Fig. 2A-F). Upregulation was not as evident for the HlMYB3, HlMYB7 and HlMYB8 transcription factors as for HlMYB78, HlbHLH2 and HlWRKY1. HlMYB78 was essentially downregulated in leaf tissue (Fig. 2D) and HlbHLH2 in apex tissue (Fig. 2E). There were no significant differences between aromatic and bitter hop cultivars, except for HlMYB3 in leaves (Fig. 2A) and HlbHLH2 in lupulin glands (Fig. 2E).

Fig. 2
figure 2

The relative expression of bitter acid biosynthesis transcription factors A: HlMYB3, B: HlMYB7, C: HlMYB8, D: HlMYB78, E: HlbHLH2 and F: HlWRKY1 relative to reference genes (see Table 1) in different hop tissues. Tissues: A – apex, B – leaf, C – bract without lupulin, D – flower, E – young cone, and F – lupulin gland. Significant t-test group differences at the following probability levels: * - P < 0,1, ** - P < 0,05, and *** - P < 0,01. straight line - median, dashed line – average, and box – 95% percentile ± standard deviation

Transcription factors involved in trichome initiation and development

The alpha bitter acid content is dependent on lupulin gland density; therefore, the regulation of glandular trichome initiation and development is important for breeding and hop cultivar improvement. Knowledge of the trichome developmental regulation network is limited in hop; therefore, we tried to find key regulatory factors for relative differential expression within aromatic and bitter hop cultivars.

We found that the negative regulator of trichome development HlETC1 was 3 times more highly expressed in flowers of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 3A). However, we did not find any differences for similar HlCPC negative regulators, and their expression profiles were different (Fig. 3B). Two other MYB-related transcription factors, HlMYB61 and HlMYB5, were both 2.2 times more highly expressed in young cones of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 3C and D). Differential expression was also found in leaves for HlMYB61 (Fig. 3C) and apexes for HlMYB5 (Fig. 3D).

Fig. 3
figure 3

The relative expression of reference genes for trichome-regulating transcription factors A: HlETC1, B: HlCPC, C: HlMYB61, D: HlMYB5, E: HlGL2, F: HlCYCB2–4, G: HlHDG2, H: HlPDF2, I: HlHB51, J: HlZHD6, K: HlHB14, L: HlCYCNB1, M: HlHDG11, N: HlZFP4, O: HlZFP5, P: HlZFP6 Q: HlZFP8, R: HlGIS3, S: HlYABBY1, T: HlMYB106, U: HlWRKY44, V: HlGL3, W: HlRAV1, X: HlABIL2 and Y: HlABIL3 (see Table 1) in different tissues of aromatic and bitter hop cultivars. Tissues: A – apex, B – leaf, C – bract without lupulin, D – flower, E – young cone, and F – lupulin gland. Significant t-test group differences at the following probability levels: * - P < 0,1, ** - P < 0,05, and *** - P < 0,01. straight line - median, dashed line – average, and box – 95% percentile ± standard deviation

For homeobox-leucine zipper proteins, HlGLABRA2 was 4.2 times more highly expressed in young cones and 2.8 times more highly expressed in leaves of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 3E). The relative expression of the HlCYCB2–4 homologue was 3 times higher in young cones of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 3F). However, the Wo gene homologues HlHDG2 and HlPDF2 were not differentially expressed in any tissue (Fig. 3G and H). HlHB51 was 5 times more highly expressed in young cones of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 3I). HlZHD6, HlHB14 and HlCYCNB1 had similar expression profiles across hop tissues (Fig. 3J, K and L). HlHDG11 was upregulated in lupulin glands (Fig. 3M).

Of the five studied C2H2 zinc finger proteins (Fig. 3N-R), HlZFP8 was 4.6 times more highly expressed in young cones of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 3P). HlZFP4 (Fig. 3N) and HlZFP5 (Fig. 3O) were also nonsignificantly upregulated in young cones.

Regarding other regulatory factors, the axial regulator HlYABBY 1 was 3.7 times more highly expressed in young cones of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 3S). The MIXTA homologue HlMYB106 was upregulated in young cones compared to other tissues with no differences between aromatic and bitter hop cultivars (Fig. 3T). The relative expression of the HlWRKY44 factor was 3.8 times higher in the lupulin glands of aromatic hop cultivars than in those of bitter hop cultivars (Fig. 3U). There were no significant differences in the bHLH factor HlGL3 (Fig. 3V) and AP2/ERF and B3 domain-containing transcription factor HlRAV1 (Fig. 3W) between aromatic and bitter hop cultivars. For Abl interactor-like proteins, HlABIL2 was 3.3 and 3.2 times more highly expressed in flowers and lupulin glands of aromatic hop cultivars than in those of bitter hop cultivars (Fig. 3X). However, HlABIL3 was 3.4 and 2.1 times more highly expressed in the apexes and young cones of bitter hop cultivars than in those of aromatic hop cultivars (Fig. 3Y).

Discussion

Genes involved in bitter acid biosynthesis

Today, all bitter acid biosynthetic pathway enzyme genes [3,4,5,6,7,8,9,10] that are predominantly and highly expressed in lupulin glands are known. Expression differences between aromatic and bitter hop cultivars have not yet been studied. We confirmed that bitter acid biosynthetic genes are upregulated in lupulin glands and downregulated in leaves and bracts. Significant differences between aromatic and bitter hop cultivars were found for the HlBCAT1 and HlPT1L genes. HlBCAT1 is localized to the mitochondria and is involved in amino acid biosynthesis [4]. Upregulation of this gene in the apex, bract and lupulin gland of aromatic cultivars is connected with the amino acid composition. HlPT1L is also involved in prenylflavonoid biosynthesis [8, 9], and upregulation of this gene in flowers of bitter cultivars is connected with higher levels of these compounds. Recently, two monooxygenase 2 genes [3] (HlHS1 and HlHS2) required for the last step of alpha bitter acid biosynthesis were found to be upregulated in the lupulin glands, flowers and young cones of bitter hop cultivars. From these findings, it can be supposed that HlHS genes are important for the alpha bitter acid content in hop cones together with lupulin gland density [2, 16, 17].

Transcription factors involved in bitter acid biosynthesis pathways

Bitter acid biosynthesis pathway genes are regulated by transcription factor networks [11,12,13,14]. We studied the relative expression of six known transcription factors and confirmed their upregulation in lupulin glands. Two factors, HlMYB3 and HlbHLH2, were differentially expressed in aromatic and bitter hop cultivars in the leaves or lupulin glands, respectively. Both genes are involved in prenylflavonoid biosynthesis [11,12,13,14], and differences in this pathway are connected to these compounds. Differences in HS gene expression have not been found to be related to any of the studied transcription factors. It is likely that either other transcription factors or changes in cis-acting regulatory DNA elements on promoter sequences can produce this upregulation. A more detailed study will be necessary in the future.

Transcription factors involved in trichome initiation and development

Because lupulin gland density is correlated with alpha bitter acid content [2, 16, 17], we studied 25 candidate regulatory factors and genes involved in glandular trichome initiation and development.

One of the essential regulators responsible for glandular trichome density formation [19, 20], the negative regulator of trichome development HlETC1, was upregulated in flowers of bitter hop cultivars. Its expression was downregulated in lupulin glands, similar to the results found in a previous study [1]. A similar negative regulator of HlCPC was differentially expressed in hop tissues. The amino acid (AA) sequence of HlETC1 has 55.7% homology to Arabidopsis thaliana AtCPC, 54.2% homology to AtETC3, 48% homology to AtETC1 and 43.2% homology to AtETC2. Although the AA homology between HlETC1 and HlCPC was only 43.6%, the homology of CPC to A. thaliana genes was similar to that of ETC1 (61.8% to AtCPC, 61.5% to AtETC2, 56.9% to AtETC1 and 44.7% to AtETC3). Therefore, we cannot determine the exact role of this HlETC1 gene in controlling trichome formation.

Other MYB-related transcription factors, such as HlMYB61 and HlMYB5, were also upregulated in young cones of bitter hop cultivars. They are act as GLABRA1 transcription factors in the regulation of mucilage synthesis, seed coat development and trichome morphogenesis in A. thaliana [28] or as AaMYB1 transcription factors in the regulation of terpene synthesis and trichome development in Artemisia annua [29]. Their pleiotropic function was identified based on their differential expression in leaf and apex tissues. These genes may be part of generally conserved MYB-bHLH-WD40 complex regulation, which is involved in lupulin gland development. This result was confirmed by the upregulation of HlGLABRA2 in young cones and leaves of bitter hop cultivars. The GL2 gene possibly regulates the expression of SlCycB2 during the initiation and development of type I trichomes in tomato [21, 22], and its homologue HlCYCB2–4 was also upregulated in young cones of bitter hop cultivars.

The previously identified lupulin-upregulated homeobox-leucine zipper protein HlHB51 was also upregulated in young cones of bitter hop cultivars [1]. ATHB-51 interacts with the meristem regulator LEAFY and participates in trichome formation, leaf morphogenesis and floral meristem determinacy in Arabidopsis [30], and it is an essential regulator of multicellular trichome development in Cucumis sativum [31]. Similarly expressed homeobox-leucine zipper proteins HlZHD6, HlHB14 and HlCYCNB1 probably participate in meristem determinacy and development. HlHDG11 was similarly upregulated in lupulin glands in a previous study [1]. This gene is involved in epidermal cell differentiation, and HDG11 mutants show excess branching of the trichome in A. thaliana [32].

In a previous study, the axial regulator HlYABBY1 was significantly downregulated in lupulin glands and young cones [1], as in our results. We can suppose from its upregulation in young cones of bitter hop cultivars that this gene may be involved in lupulin gland development at abaxial sites of bracts and bracteoles [33].

The MIXTA transcription factor is also considered a positive regulator of glandular trichome formation [19, 20]. It was obvious from our results for its homologue HlMYB106 that this gene is a general transcription factor for glandular trichome development.

HlWRKY44, TRANSPARENT TESTA GLABRA2 (TTG2) factor, regulates trichome-specific TTG1, GL2, GL3 and TRY factors [34]. Its downregulation in the lupulin glands of bitter hop cultivars probably results in downregulation of the MYB-related transcription factor TRY, a negative regulator of trichome development. However, WRKY44 is also involved in the biosynthesis of proanthocyanins in fruits [35] or mucilage and tannins in seed coats [36].

The homologue of the bHLH factor HlGL3 was downregulated in lupulin glands without differences between aromatic and bitter hop cultivars. Therefore, we cannot confirm that this gene is a part of the MYB-bHLH-WD40 complex, and thorough combinatorial analysis will be necessary [11].

Expression differences in hop cone tissues were previously found [1] for AP2/ERF and the B3 domain-containing transcription factor HlRAV1, which were confirmed despite a lack of influence on glandular trichome density or alpha bitter acid content. Other genes that were overexpressed in lupulin glands include the Abl interactor-like proteins HlABIL2 or HlABIL3 [1], which are involved in the regulation of actin and microtubule organization as part of a SCAR/WAVE complex that activates the ARP2/3 complex [37]. Based on the expression differences in apexes, flowers, young cones or lupulin glands, we cannot exactly determine what aspect of trichome, cone or seed development is influenced, even if these genes play a role in the development of plant cell shape and can cause a distortion in the trichome phenotype [38].

Conclusions

Gene expression analyses enabled us to identify 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. Humulone synthase genes (two analogues) are important for alpha bitter acid content in the glandular trichomes of hop cones. Differential gene expression analyses showed that the MYB transcription factors HlETC1, HlMYB61 and HlMYB5, homeobox-leucine zipper protein HlGLABRA 2, C2H2 zinc finger protein HlZFP8 and axial regulator HlYABBY1 may be key regulatory factors for lupulin gland density and morphological and developmental differentiation of glandular trichomes inside hop cones. Analyses also showed that glandular trichome initiation and development are controlled by a large regulatory network and should be studied in detail in the future.

Methods

Plant materials

The hop plants used were grown under standard agronomic conditions in experimental fields on the Steknik farm of the Hop Research Institute in Zatec (Saaz), CR. Based on general long-term knowledge, we selected the “aromatic” hop cultivars Saaz, Fuggle, Hallertauer and Kazbek (alpha bitter acid contents varying from 2.5 to 8.0% and beta bitter acid contents varying from 2.0 to 6.0% in cones) and the “bitter” hop cultivars Vital, Herkules, Columbus, and Magnum (alpha bitter acid contents varying from 11.0 to 18.0% and beta bitter acid contents varying from 4.0 to 10.0% in cones) as contrasting hop genotypes for group analyses [2]. In 2020, three pooled samples of eight plants for each hop cultivar (Saaz, Fuggle, Hallertauer, Vital, Herkules, and Columbus) were collected from different tissues during development: apexes in April (only Saaz and Vital), flowers in July, leaves and young cones in August, and mature cones in August and September. Samples collected in August 2016 were used to increase the tissue variability. Three pooled samples containing one kilogram fresh weight of mature cones were used for bract and lupulin gland analyses in the hop cultivars Saaz, Fuggle, Hallertauer, Kazbek, Vital, Columbus and Magnum. Three pooled samples of four plants of the hop cultivars Saaz, Fuggle Columbus and Magnum collected in May 2013 were used to increase the variability of leaf and apex tissues. All samples were immediately frozen in liquid nitrogen and stored in a deep freezer (− 80 °C). Lupulin glands were separated from mature cones by agitation in liquid nitrogen followed by filtration through a 1 mm metal sieve to remove cone debris [1]. The glands were recovered from liquid nitrogen, and the rest of the cones were analysed as bracts without lupulin. We did not measure either the bitter acid contents or numbers of lupulin glands in hop cone samples.

RNA isolation and gene expression analyses

RNA was isolated using PureLink™ Plant RNA Reagent (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s protocol and purified by DNaseI treatment on a column (RNeasy Plant Mini Kit, Qiagen, Hilden, FRG) [11]. RNA samples were reverse transcribed by oligo (dT)18 primer and a First Strand cDNA Synthesis Kit (Roche Diagnostics, Mannheim, Germany) at 50 °C for 60 min. NGS genome information [39] and the HopBase database [38, 40] were used to search for bitter acid biosynthesis [3] and trichome-specific transcription factor [1] gene sequences (Table 1). Genes were selected based on previous works [1, 3,4,5,6,7,8,9,10,11,12,13,14, 18,19,20]. Advanced BLAST 2.0 (http://www.ncbi.nlm.nih.gov/blast/blast.cgi) and the MegAlign module (LASERGENE system v. 7.1, DNAStar, Madison, WI, USA) were used to evaluate sequence homology. RealTimeDesign software (LGC Biosearch Technologies, Petaluma, CA, USA) was used to design real-time PCR primers (Table S1). A total of 2 μl of 50 × diluted cDNA was used in a 20 μl PCR with iTaq universal SYBR green Supermix (Bio–Rad Laboratories, Hercules, CA, USA) in a CFX Connect real-time PCR cycler (Bio–Rad Laboratories, Hercules, CA, USA). Five reference genes (Table 1), which were found to be constitutively expressed in different tissues [11] by NormFinder in GenEx v.6.0.1.612 (MultiD Analyzes AB, Gothenburg, Sweden), were used for normalization of the samples. The relative expression of five reference genes was calculated by the “delta-delta method” (RE = 2-∆CT) for each sample. SigmaPlot for Windows v.10.0.0.54 (Systat Software Inc., San Jose, CA) was used for statistical group (average, standard deviation) and unpaired two-sample t-test analyses of relative expression.

Table 1 List of analysed bitter acid biosynthesis gene, trichome-specific transcription factor and reference gene sequences

Availability of data and materials

The datasets used and/or analysed during the current study are available from the author on reasonable request.

Abbreviations

A. thaliana :

Arabidopsis thaliana

BCAA:

Branched-chain amino acids

bHLH:

Basic-helix-loop-helix

BLAST:

Basic local alignment search tool

C2H2:

Cys2His2

MBW:

MYB-bHLH-WD40

MEP:

Methylerythritol phosphate

MYB:

Myeloblastosis

MYC:

Myelocytomatosis

NGS:

Next generation sequencing

qRT–PCR:

Quantitative real-time polymerase chain reaction

RNA:

Ribonucleic acid

WD:

Tryptophan-aspartic acid dipeptide

References

  1. Mishra AK, Kocábek T, Nath VS, Awasthi PA, Shreshta A, Killi UK, et al. Dissection of dynamic transcriptome landscape of leaf, bract, and lupulin gland in hop (Humulus lupulus L.). Int J Mol Sci. 2020;21:233.

    Article  CAS  Google Scholar 

  2. Patzak J, Krofta K, Henychová A, Nesvadba V. Number and size of lupulin glands, glandular trichomes of hop (Humulus lupulus L.), play a key role in contents of bitter acids and polyphenols in hop cone. Int J Food Sci Technol. 2015;50:1864–72.

    Article  CAS  Google Scholar 

  3. Champagne A, Boutry M. A comprehensive proteome map of glandular trichomes of hop (Humulus lupulus L.) female cones: Identification of biosynthetic pathways of the major terpenoid-related compounds and possible transport proteins. Proteomics. 2017;17:1600411.

    Article  CAS  Google Scholar 

  4. Clark SM, Vaitheeswaran V, Ambrose SJ, Purves RW, Page JE. Transcriptome analysis of bitter acid biosynthesis and precursor pathways in hop (Humulus lupulus). BMC Plant Biol. 2013;13:12.

    Article  CAS  Google Scholar 

  5. Xu H, Zhang F, Liu B, Huhman DV, Sumner LW, Dixon RA, et al. Characterization of the formation of branched short-chain fatty acid:CoAs for bitter acid biosynthesis in hop glandular Trichomes. Mol Plant. 2013;6:1301–17.

    Article  CAS  Google Scholar 

  6. Okada Y, Ito K. Cloning and analysis of valerophenone synthase gene expressed specifically in lupulin gland of hop (Humulus lupulus L.). Biosci Biotech Biochem. 2001;65:150–5.

    Article  CAS  Google Scholar 

  7. Okada K. The biosynthesis of isoprenoids and the mechanisms regulating it in plants. Biosci Biotechnol Biochem. 2011;75:1219–25.

    Article  CAS  Google Scholar 

  8. Tsurumaru Y, Sasaki K, Miyawaki T, Momma T, Umemoto N, Yazaki K. An aromatic prenyltransferase-like gene HlPT-1 preferentially expressed in lupulin glands of hop. Plant Biotechnol. 2010;27:199–204.

    Article  CAS  Google Scholar 

  9. Tsurumaru Y, Sasaki K, Miyawaki T, Uto Y, Momma T, Umemoto N, et al. HlPT-1, a membrane-bound prenyltransferase responsible for the biosynthesis of bitter acids in hops. Biochem Biophys Res Comm. 2012;417:393–8.

    Article  CAS  Google Scholar 

  10. Li H, Ban Z, Qin H, Ma L, King AJ, Wang G. A heteromeric membrane-bound prenyltransferase complex from hop catalyzes three sequential aromatic prenylations in the bitter acid pathway. Plant Physiol. 2015;167:650–9.

    Article  CAS  Google Scholar 

  11. Matoušek J, Kocábek T, Patzak J, Füssy Z, Procházková J, Heyerick A. Combinatorial analysis of lupulin gland transcription factors from R2R3Myb, bHLH and WDR families indicates a complex regulation of chs_H1 genes essential for prenylflavonoid biosynthesis in hop (Humulus lupulus L.). BMC Plant Biol. 2012;12:27.

    Article  CAS  Google Scholar 

  12. Mishra AK, Duraisamy GS, Khare M, Kocábek T, Jakse J, Bříza J, et al. Genome-wide transcriptome profiling of transgenic hop (Humulus lupulus L.) constitutively overexpressing HlWRKY1 and HlWDR1 transcription factors. BMC Genom. 2018;19:739.

    Article  CAS  Google Scholar 

  13. Matoušek J, Kocábek T, Patzak J, Bříza J, Siglová K, Mishra AK, et al. The “putative” role of transcription factors from HlWRKY family in the regulation of the final steps of prenylflavonid and bitter acids biosynthesis in hop (Humulus lupulus L.). Plant Mol Biol. 2016;92:263–77.

    Article  CAS  Google Scholar 

  14. Kocábek T, Mishra AK, Matoušek J, Patzak J, Lomnická A, Khare M, et al. The R2R3 transcription factor HlMYB8 and its role in flavonoid biosynthesis in hop (Humulus lupulus L.). Plant Sci. 2018;269:32–46.

    Article  CAS  Google Scholar 

  15. Barth-Haas Group. Barth-Haas Hops Companion: A brewer's guide to hop varieties and hop products. In: Kostelecky T, editor. . 3rd ed. Yakima: John I. Haas Incorporated; 2016. p. 183. ISBN: 978-5323-1361-5.

    Google Scholar 

  16. Sabo J, Kisgeci J, Ikic I. Content of active components in dependence on the number of lupulin glands in the hop cones. Rostlinna Vyroba. 2001;47:201–4.

    CAS  Google Scholar 

  17. Srecec S, Zechner-Krpan V, Marag S, Spoljaric I, Kvaternjak I, Mrsic G. Morphogenesis, volume and number of hop (Humulus lupulus L.) glandular trichomes, and their influence on alpha-acid accumulation in fresh bracts of hop cones. Acta Botanica Croatica. 2011;70:1–8.

    Article  Google Scholar 

  18. Huchelmann A, Boutry M, Hachez C. Plant glandular Trichomes: natural cell factories of high biotechnological interest. Plant Physiol. 2017;175:6–22.

    Article  CAS  Google Scholar 

  19. Schuurink R, Tissier A. Glandular trichomes: micro-organs with model status? New Phytol. 2020;225:2251–66.

    Article  Google Scholar 

  20. Hauser MT. Molecular basis of natural variation and environmental control of trichome patterning. Front Plant Sci. 2014;5:320.

    Article  Google Scholar 

  21. Yang C, Li H, Zhang J, Wang T, Ye Z. Fine-mapping of the woolly gene controlling multicellular trichome formation and embryonic development in tomato. Theor Appl Genet. 2011;123:625–33.

    Article  CAS  Google Scholar 

  22. Gao S, Gao Y, Xiong C, Yu G, Chang J, Yang Q, et al. The tomato B-type cyclin gene, SlCycB2, plays key roles in reproductive organ development, trichome initiation, terpenoids biosynthesis and Prodenia litura defense. Plant Sci. 2017;262:103–14.

    Article  CAS  Google Scholar 

  23. Chang J, Yu T, Yang Q, Li C, Xiong C, Gao S, et al. Hair, encoding a single C2H2 zinc-finger protein, regulates multicellular trichome formation in tomato. Plant J. 2018;96:90–102.

    Article  CAS  Google Scholar 

  24. Xu J, van Herwijnen ZO, Dräger DB, Sui C, Haring MA, Schuurink RC. SlMYC1 regulates type VI glandular Trichome formation and terpene biosynthesis in tomato glandular cells. Plant Cell. 2018;30:2988–3005.

    Article  CAS  Google Scholar 

  25. Yang C, Gao Y, Gao S, Yu G, Xiong C, Chang J, et al. Transcriptome profile analysis of cell proliferation molecular processes during multicellular trichome formation induced by tomato Wov gene in tobacco. BMC Genomics. 2015;16:e868.

    Article  CAS  Google Scholar 

  26. Liu Y, Liu D, Khan AR, Liu B, Wu M, Huang L, et al. NbGIS regulates glandular trichome initiation through GA signaling in tobacco. Plant Mol Biol. 2018;98:153–67.

    Article  CAS  Google Scholar 

  27. Yan T, Li L, Xie L, Chen M, Shen Q, Pan Q, et al. A novel HD-ZIP IV/MIXTA complex promotes glandular trichome initiation and cuticle development in Artemisia annua. New Phytol. 2018;218:567–78.

    Article  CAS  Google Scholar 

  28. Li SF, Milliken ON, Pham H, Seyit R, Napoli R, Preston J, et al. The Arabidopsis MYB5 transcription factor regulates mucilage synthesis, seed coat development, and trichome morphogenesis. Plant Cell. 2009;21:72–89.

    Article  CAS  Google Scholar 

  29. Matías-Hernández L, Jiang W, Yang K, Tang K, Brodelius PE, Pelaz S. AaMYB1 and its orthologue AtMYB61 affect terpene metabolism and trichome development in Artemisia annua and Arabidopsis thaliana. Plant J. 2017;90:520–34.

    Article  CAS  Google Scholar 

  30. Saddic LA, Huvermann B, Bezhani S, Su Y, Winter CM, Kwon CS, et al. The LEAFY target LMI1 is a meristem identity regulator and acts together with LEAFY to regulate expressionof CAULIFLOWER. Development. 2006;133:1673–82.

    Article  CAS  Google Scholar 

  31. Zhao J-L, Pan J-S, Guan Y, Nie J-T, Yang J-J, Qu M-L, et al. Transcriptome analysis in Cucumis sativus identifies genes involved in multicellular trichome development. Genomics. 2015;105:296–303.

    Article  CAS  Google Scholar 

  32. Nakamura M, Katsumata H, Abe M, Yabe N, Komeda Y, Yamamoto KT, et al. Characterization of the class IV homeodomain-leucine zipper gene family in Arabidopsis. Plant Physiol. 2006;141:1363–75.

    Article  CAS  Google Scholar 

  33. Filyushin MA, Slugina MA, Shchennikova AV, Kochieva EZ. Identification and expression analysis of the YABBY1 Gene in wild tomato species. Russ J Genet. 2018;54:536–47.

    Article  CAS  Google Scholar 

  34. Pesch M, Dartan B, Birkenbihl R, Somssich IE, Hülskamp M. Arabidopsis TTG2 regulates TRY expression through enhancement of activator complex-triggered activation. Plant Cell. 2014;26:4067–83.

    Article  CAS  Google Scholar 

  35. Peng Y, Thrimawithana AH, Cooney JM, Jensen JD, Espley RV, Allan AC. The proanthocyanin-related transcription factors MYBC1 and WRKY44 regulate branch points in the kiwifruit anthocyanin pathway. Sci Rep. 2020;10:14161.

    Article  CAS  Google Scholar 

  36. Li SF, Allen PJ, Napoli RS, Browne RG, Pham H, Parish RW. MYB–bHLH–TTG1 regulates Arabidopsis seed coat biosynthesis pathways directly and indirectly via multiple tiers of transcription factors. Plant Cell Physiol. 2020;61:1005–18.

    Article  CAS  Google Scholar 

  37. Jörgens CI, Grünewald N, Hülskamp M, Uhrig JF. A role for ABIL3 in plant cell morphogenesis. Plant J. 2010;62:925–35.

    PubMed  Google Scholar 

  38. Padgitt-Cobb LK, Kingan SB, Wells J, Elser J, Kronmiller B, Moore D, et al. A draft phased assembly of the diploid Cascade hop (Humulus lupulus) genome. Plant Genome. 2021:e20072. https://doi.org/10.1002/tpg2.20072.

  39. Natsume S, Takagi H, Shiraishi A, Murata J, Toyonaga H, Patzak J, et al. The draft genome of hop (Humulus lupulus), an essence for brewing. Plant Cell Physiol. 2015;56:428–41.

    Article  CAS  Google Scholar 

  40. Hill ST, Sudarsanam R, Henning J, Hendrix D. HopBase: a unified resource for Humulus genomics. Database. 2017, 2017:bax009.

Download references

Acknowledgements

We thank the Breeding Department staff of the Hop Research Institute in Žatec for the maintenance of plant materials on the experimental farm in Steknik.

Funding

This work was supported by the Grant Agency of the Czech Republic in project 19-19629S.

Author information

Authors and Affiliations

Authors

Contributions

JP contributed to the project development, gene sequence searching, and primer design, performed the statistical analysis and interpreted the results. AH performed the gene expression analyses and interpreted the results. JM contributed to the project development, candidate gene selection and interpretation of the results. All authors have read and agreed to the published version of the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Josef Patzak.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Patzak, J., Henychová, A. & Matoušek, J. Developmental regulation of lupulin gland-associated genes in aromatic and bitter hops (Humulus lupulus L.). BMC Plant Biol 21, 534 (2021). https://doi.org/10.1186/s12870-021-03292-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-021-03292-z

Keywords