Development of hop transcriptome to support research into host-viroid interactions

Optimization and Validation of Singleplex and Multiplex RT-qPCR for Detection of Citrus bark cracking viroid (CBCVd), Hop latent viroid (HLVd), and Hop stunt viroid (HSVd) in Hops (Humulus lupulus)

Direct crop losses due to plant diseases and the measures used to control them have significant agricultural and economic impacts. The shift from diverse small-scale to large-scale genetically uniform monoculture production, along with agricultural intensification and climate change, has led to several known epidemics in man-made agroecosystems that have been rendered more vulnerable to pathogens. One such example is hop growing, which is threatened by highly aggressive hop viroids. Since 2007, almost one-third (about 500 ha) of Slovenian hop gardens have been affected by severe hop stunt disease caused by Citrus bark cracking viroid (CBCVd), which continues to spread despite strict prevention measures. We have developed and validated a multiplex RT-qPCR (mRT-qPCR) for the sensitive detection of CBCVd, Hop latent viroid (HLVd), and Hop stunt viroid (HSVd). Singleplex RT-qPCR assays were designed individually and subsequently combined in a one-step mRT-qPCR assay. Hop-specific mRNA170 and mRNA1192 internal controls were also developed to detect possible PCR inhibition. Analytical specificity was tested on 35 samples from different hosts, geographic regions, and combinations of viroids. Method validation showed that mRT-qPCR had lower sensitivity than singleplex RT-qPCR, while specificity, selectivity, repeatability, and reproducibility remained unchanged. The newly developed assays were found to be robust, reliable, and suitable for large-scale screening of hop viroids.

Cytosine Methylation in Genomic DNA and Characterization of DNA Methylases and Demethylases and Their Expression Profiles in Viroid-Infected Hop Plants (Humulus lupulus Var. 'Celeia')

Abiotic and biotic stresses can lead to changes in host DNA methylation, which in plants is also mediated by an RNA-directed DNA methylation mechanism. Infections with viroids have been shown to affect DNA methylation dynamics in different plant hosts. The aim of our research was to determine the content of 5-methylcytosine (5-mC) in genomic DNA at the whole genome level of hop plants (Humulus lupulus Var. 'Celeia') infected with different viroids and their combinations and to analyse the expression of the selected genes to improve our understanding of DNA methylation dynamics in plant-viroid systems. The adapted HPLC-UV method used proved to be suitable for this purpose, and thus we were able to estimate for the first time that the cytosine methylation level in viroid-free hop plants was 26.7%. Interestingly, the observed 5-mC level was the lowest in hop plants infected simultaneously with CBCVd, HLVd and HSVd (23.7%), whereas the highest level was observed in plants infected with HLVd (31.4%). In addition, we identified three DNA methylases and one DNA demethylase gene in the hop's draft genome. The RT-qPCR revealed upregulation of all newly identified genes in hop plants infected with all three viroids, while no altered expression was observed in any of the other hop plants tested, except for CBCVd-infected hop plants, in which one DNA methylase was also upregulated.

Transcriptomics Advancement in the Complex Response of Plants to Viroid Infection

Viroids are the smallest plant pathogens, consisting of a single-stranded circular RNA of less than 500 ribonucleotides in length. Despite their noncoding nature, viroids elicit disease symptoms in many economically important plant hosts, and are, thus, a class of pathogens of great interest. How these viroids establish disease within host plants, however, is not yet fully understood. Recent transcriptomic studies have revealed that viroid infection influences the expression of genes in several pathways and processes in plants, including defence responses, phytohormone signalling, cell wall modification, photosynthesis, secondary metabolism, transport, gene expression and protein modification. There is much debate about whether affected pathways signify a plant response to viroid infection, or are associated with the appearance of disease symptoms in these interactions. In this review, we consolidate the findings of viroid–host transcriptome studies to provide an overview of trends observed in the data. When considered together, changes in the gene expression of different hosts upon viroid infection reveal commonalities and differences in diverse interactions. Here, we discuss whether trends in host gene expression can be correlated to plant defence or disease development during viroid infection, and highlight avenues for future research in this field.

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.

Identification and Characterization of Verticillium nonalfalfae-Responsive MicroRNAs in the Roots of Resistant and Susceptible Hop Cultivars

MicroRNAs are 21- to 24-nucleotide-long, non-coding RNA molecules that regulate gene expression at the post-transcriptional level. They can modulate various biological processes, including plant response and resistance to fungal pathogens. Hops are grown for use in the brewing industry and, recently, also for the pharmaceutical industry. Severe Verticillium wilt caused by the phytopathogenic fungus Verticillium nonalfalfae, is the main factor in yield loss in many crops, including hops (Humulus lupulus L.). In our study, we identified 56 known and 43 novel miRNAs and their expression patterns in the roots of susceptible and resistant hop cultivars after inoculation with V. nonalfalfae. In response to inoculation with V. nonalfalfae, we found five known and two novel miRNAs that are differentially expressed in the susceptible cultivar and six known miRNAs in the resistant cultivar. Differentially expressed miRNAs target 49 transcripts involved in protein localization and pigment synthesis in the susceptible cultivar, whereas they are involved in transcription factor regulation and hormone signalling in the resistant cultivar. The results of our study suggest that the susceptible and resistant hop cultivars respond differently to V. nonalfalfae inoculation at the miRNA level and that miRNAs may contribute to the successful defence of the resistant cultivar.

Core RNA Interference Genes Involved in miRNA and Ta-siRNA Biogenesis in Hops and Their Expression Analysis after Challenging with Verticillium nonalfalfae

RNA interference is an evolutionary conserved mechanism by which organisms regulate the expression of genes in a sequence-specific manner to modulate defense responses against various abiotic or biotic stresses. Hops are grown for their use in brewing and, in recent years, for the pharmaceutical industry. Hop production is threatened by many phytopathogens, of which Verticillium, the causal agent of Verticillium wilt, is a major contributor to yield losses. In the present study, we performed identification, characterization, phylogenetic, and expression analyses of three Argonaute, two Dicer-like, and two RNA-dependent RNA polymerase genes in the susceptible hop cultivar Celeia and the resistant cultivar Wye Target after infection with Verticillium nonalfalfae. Phylogeny results showed clustering of hop RNAi proteins with their orthologues from the closely related species Cannabis sativa, Morus notabilis and Ziziphus jujuba which form a common cluster with species of the Rosaceae family. Expression analysis revealed downregulation of argonaute 2 in both cultivars on the third day post-inoculation, which may result in reduced AGO2-siRNA-mediated posttranscriptional gene silencing. Both cultivars may also repress ta-siRNA biogenesis at different dpi, as we observed downregulation of argonaute 7 in the susceptible cultivar on day 1 and downregulation of RDR6 in the resistant cultivar on day 3 after inoculation.

Pilot RNA-seq data from 24 species of vascular plants at Harvard Forest

PREMISE

Large-scale projects such as the National Ecological Observatory Network (NEON) collect ecological data on entire biomes to track climate change. NEON provides an opportunity to launch community transcriptomic projects that ask integrative questions in ecology and evolution. We conducted a pilot study to investigate the challenges of collecting RNA-seq data from diverse plant communities.

METHODS

We generated >650 Gbp of RNA-seq for 24 vascular plant species representing 12 genera and nine families at the Harvard Forest NEON site. Each species was sampled twice in 2016 (July and August). We assessed transcriptome quality and content with TransRate, BUSCO, and Gene Ontology annotations.

RESULTS

Only modest differences in assembly quality were observed across multiple k-mers. On average, transcriptomes contained hits to >70% of loci in the BUSCO database. We found no significant difference in the number of assembled and annotated transcripts between diploid and polyploid transcriptomes.

DISCUSSION

We provide new RNA-seq data sets for 24 species of vascular plants in Harvard Forest. Challenges associated with this type of study included recovery of high-quality RNA from diverse species and access to NEON sites for genomic sampling. Overcoming these challenges offers opportunities for large-scale studies at the intersection of ecology and genomics.

Impact of Nucleic Acid Sequencing on Viroid Biology

The early 1970s marked two breakthroughs in the field of biology: (i) The development of nucleotide sequencing technology; and, (ii) the discovery of the viroids. The first DNA sequences were obtained by two-dimensional chromatography which was later replaced by sequencing using electrophoresis technique. The subsequent development of fluorescence-based sequencing method which made DNA sequencing not only easier, but many orders of magnitude faster. The knowledge of DNA sequences has become an indispensable tool for both basic and applied research. It has shed light biology of viroids, the highly structured, circular, single-stranded non-coding RNA molecules that infect numerous economically important plants. Our understanding of viroid molecular biology and biochemistry has been intimately associated with the evolution of nucleic acid sequencing technologies. With the development of the next-generation sequence method, viroid research exponentially progressed, notably in the areas of the molecular mechanisms of viroids and viroid diseases, viroid pathogenesis, viroid quasi-species, viroid adaptability, and viroid–host interactions, to name a few examples. In this review, the progress in the understanding of viroid biology in conjunction with the improvements in nucleotide sequencing technology is summarized. The future of viroid research with respect to the use of third-generation sequencing technology is also briefly envisaged.

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.

Root Transcriptomic Analysis Reveals Global Changes Induced by Systemic Infection of Solanum lycopersicum with Mild and Severe Variants of Potato Spindle Tuber Viroid

Potato spindle tuber viroid (PSTVd) causes systemic infection in plant hosts. There are many studies on viroid-host plant interactions, but they have predominantly focused on the aboveground part of the plant. Here, we investigated transcriptomic profile changes in tomato roots systemically infected with mild or severe PSTVd variants using a combined microarray/RNA-seq approach. Analysis indicated differential expression of genes related to various Gene Ontology categories depending on the stage of infection and PSTVd variant. A majority of cell-wall-related genes were down-regulated at early infection stages, but at the late stage, the number of up-regulated genes increased significantly. Along with observed alterations of many lignin-related genes, performed lignin quantification indicated their disrupted level in PSTVd-infected roots. Altered expression of genes related to biosynthesis and signaling of auxin and cytokinin, which are crucial for lateral root development, was also identified. Comparison of both PSTVd infections showed that transcriptional changes induced by the severe variant were stronger than those caused by the mild variant, especially at the late infection stage. Taken together, we showed that similarly to aboveground plant parts, PSTVd infection in the underground tissues activates the plant immune response.

New Male Specific Markers for Hop and Application in Breeding Program

Male specific DNA sequences were selected from a Diversity Arrays Technology (DArT) mapping study to evaluate their suitability for determination of the sex phenotype among young seedlings in a hop (Humulus lupulus L.) breeding program. Ten male specific DArT markers showed complete linkage with male sex phenotype in three crossing families. Following optimization, four were successfully converted into PCR markers and a multiplex PCR approach for their use was developed. Among 197 plants (97 from the world collection; 100 from three segregating families), 94-100% positive correlation with sex phenotypic data was achieved for the single PCR amplification, whereas the multiplex approach showed 100% correlation. To develop a fast and low-cost method, crude sample multiplex PCR was evaluated in 253 progenies from 14 segregating populations without losing accuracy. The study describes, for the first time, the routine application of molecular markers linked to male sex in an intensive Slovenian hop breeding program. The methods described could be employed for screening of sex at the seedling stage in other hop programs worldwide, thereby saving resources for desirable female plants.

Temporal and spatial assessment of defence responses in resistant and susceptible hop cultivars during infection with Verticillium nonalfalfae

Hop (Humulus lupulus L.) is an important industrial plant providing ingredients for brewing and pharmaceutical industry worldwide. Its intensive production is challenged by numerous diseases. One of the most lethal and difficult to control is verticillium wilt, a vascular disease caused by the fungal pathogen Verticillium nonalfalfae. The disease can be successfully controlled by the host resistance. Despite various studies that already researched resistance mechanisms of hops, only limited number of resistance genes and markers that could be utilized for efficient resistance breeding has been identified. In this study we aimed to follow fungus colonization pattern and the differential expression of selected genes during pre-symptomatic period of susceptible (Celeia) and resistant (Wye Target) hop cultivars. Results of gene expressions and fungal colonisation of compatible and incompatible interactions with V. nonalfalfae suggest that the hop plant is challenged already at the very early fungal colonisation stages. In total, nine out of 17 gene targets investigated in our study resulted in differential expression between inoculated and control plants of susceptible and resistant cultivars. The difference was the most evident in stems at an early stage of colonisation (6 dpi), showing relatively stronger changes in targeted gene expression to infection in the resistant cultivar than in the susceptible one. Analysed gene targets are involved in the overall defence response processes of nucleic acid binding, signalling, protein ubiquitination, cell oxidative burst, hydroxylation, peroxidation, alternative splicing, and metabolite biosynthesis. The up-regulation of some genes (e.g. glycine-rich RNA-binding family protein, protein phosphatase, cysteine-rich receptor-like protein kinase, zinc finger CCCH domain-containing protein 40, cinnamic acid 4-hydroxylase, class III peroxidase, putative MAPK2, peroxiredoxin-2F) upon infection in incompatible interactions might reflect defence activation, restriction of disease spreading throughout the plant and successful response of resistant genotype.

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.

Global Transcriptomic Analysis Reveals Insights into the Response of 'Etrog' Citron (Citrus medica L.) to Citrus Exocortis Viroid Infection

Citrus exocortis viroid (CEVd) is the causal agent of citrus exocortis disease. We employed CEVd-infected ‘Etrog’ citron as a system to study the feedback regulation mechanism using transcriptome analysis in this study. Three months after CEVd infection, the transcriptome of fresh leaves was analyzed, and 1530 differentially expressed genes were detected. The replication of CEVd in citron induced upregulation of genes encoding key proteins that were involved in the RNA silencing pathway such as Dicer-like 2, RNA-dependent RNA polymerase 1, argonaute 2, argonaute 7, and silencing defective 3, as well as those genes encoding proteins that are related to basic defense responses. Many genes involved in secondary metabolite biosynthesis and chitinase activity were upregulated, whereas other genes related to cell wall and phytohormone signal transduction were downregulated. Moreover, genes encoding disease resistance proteins, pathogenicity-related proteins, and heat shock cognate 70 kDa proteins were also upregulated in response to CEVd infection. These results suggest that basic defense and RNA silencing mechanisms are activated by CEVd infection, and this information improves our understanding of the pathogenesis of viroids in woody plants.

Revisiting the Role of Transcription Factors in Coordinating the Defense Response Against Citrus Bark Cracking Viroid Infection in Commercial Hop (Humulus Lupulus L.)

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.

How sequence variants of a plastid-replicating viroid with one single nucleotide change initiate disease in its natural host

Understanding how viruses and subviral agents initiate disease is central to plant pathology. Whether RNA silencing mediates the primary lesion triggered by viroids (small non-protein-coding RNAs), or just intermediate-late steps of a signaling cascade, remains unsolved. While most variants of the plastid-replicating peach latent mosaic viroid (PLMVd) are asymptomatic, some incite peach mosaics or albinism (peach calico, PC). We have previously shown that two 21-nt small RNAs (PLMVd-sRNAs) containing a 12-13-nt PC-associated insertion guide cleavage, via RNA silencing, of the mRNA encoding a heat-shock protein involved in chloroplast biogenesis. To gain evidence supporting that such event is the initial lesion, and more specifically, that different chloroses have different primary causes, here we focused on a PLMVd-induced peach yellow mosaic (PYM) expressed in leaf sectors interspersed with others green. First, sequencing PLMVd-cDNAs from both sectors and bioassays mapped the PYM determinant at one nucleotide, a notion further sustained by the phenotype incited by other natural and artificial PLMVd variants. And second, sRNA deep-sequencing and RNA ligase-mediated RACE identified one PLMVd-sRNA with the PYM-associated change that guides cleavage, as predicted by RNA silencing, of the mRNA encoding a thylakoid translocase subunit required for chloroplast development. RT-qPCR showed lower accumulation of this mRNA in PYM-expressing tissues. Remarkably, PLMVd-sRNAs triggering PYM and PC have 5'-terminal Us, involving Argonaute 1 in what likely are the initial alterations eliciting distinct chloroses.

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.

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

Background

The hop plant (Humulus lupulus L.) is a valuable source of several secondary metabolites, such as flavonoids, bitter acids, and essential oils. These compounds are widely implicated in the beer brewing industry and are having potential biomedical applications. Several independent breeding programs around the world have been initiated to develop new cultivars with enriched lupulin and secondary metabolite contents but met with limited success due to several constraints. In the present work, a pioneering attempt has been made to overexpress master regulator binary transcription factor complex formed by HlWRKY1 and HlWDR1 using a plant expression vector to enhance the level of prenylflavonoid and bitter acid content in the hop. Subsequently, we performed transcriptional profiling using high-throughput RNA-Seq technology in leaves of resultant transformants and wild-type hop to gain in-depth information about the genome-wide functional changes induced by HlWRKY1 and HlWDR1 overexpression.

Results

The transgenic WW-lines exhibited an elevated expression of structural and regulatory genes involved in prenylflavonoid and bitter acid biosynthesis pathways. In addition, the comparative transcriptome analysis revealed a total of 522 transcripts involved in 30 pathways, including lipids and amino acids biosynthesis, primary carbon metabolism, phytohormone signaling and stress responses were differentially expressed in WW-transformants. It was apparent from the whole transcriptome sequencing that modulation of primary carbon metabolism and other pathways by HlWRKY1 and HlWDR1 overexpression resulted in enhanced substrate flux towards secondary metabolites pathway. The detailed analyses suggested that none of the pathways or genes, which have a detrimental effect on physiology, growth and development processes, were induced on a genome-wide scale in WW-transgenic lines.

Conclusions

Taken together, our results suggest that HlWRKY1 and HlWDR1 simultaneous overexpression positively regulates the prenylflavonoid and bitter acid biosynthesis pathways in the hop and thus these transgenes are presented as prospective candidates for achieving enhanced secondary metabolite content in the hop.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5125-8) contains supplementary material, which is available to authorized users.

Time-Course Microarray Analysis Reveals Differences between Transcriptional Changes in Tomato Leaves Triggered by Mild and Severe Variants of Potato Spindle Tuber Viroid

Viroids are small non-capsidated non-coding RNA replicons that utilize host factors for efficient propagation and spread through the entire plant. They can incite specific disease symptoms in susceptible plants. To better understand viroid-plant interactions, we employed microarray analysis to observe the changes of gene expression in “Rutgers” tomato leaves in response to the mild (M) and severe (S23) variants of potato spindle tuber viroid (PSTVd). The changes were analyzed over a time course of viroid infection development: (i) the pre-symptomatic stage; (ii) early symptoms; (iii) full spectrum of symptoms and (iv) the so-called ‘recovery’ stage, when stem regrowth was observed in severely affected plants. Gene expression profiles differed depending on stage of infection and variant. In S23-infected plants, the expression of over 3000 genes was affected, while M-infected plants showed 3-fold fewer differentially expressed genes, only 20% of which were specific to the M variant. The differentially expressed genes included many genes related to stress; defense; hormone metabolism and signaling; photosynthesis and chloroplasts; cell wall; RNA regulation, processing and binding; protein metabolism and modification and others. The expression levels of several genes were confirmed by nCounter analysis.

Global Transcriptomic Changes Induced by Infection of Cucumber (Cucumis sativus L.) with Mild and Severe Variants of Hop Stunt Viroid

Fifteen years after transfer to hops, hop stunt viroid-grapevine (HSVd-g) was replaced by HSVd-hop (HSVd-h), a sequence variant that contains changes at five different positions. HSVd-g54 is a laboratory mutant derived from HSVd-g that differs from its progenitor by a single G to A substitution at position 54. While infection by HSVd-h induces only mild stunting in cucumber (Cucumis sativus L.), HSVd-g54 induces much more severe symptoms in this indicator host. Comparison of transcriptome profiles of cucumber infected with HSVd-h or HSVd-g54 with those of mock-inoculated controls obtained by whole transcriptome shotgun sequencing revealed that many genes related to photosynthesis were down-regulated following infection. In contrast, genes encoding RNA-dependent RNA polymerase 1 (CsRDR1), especially CsRDR1c1 and CsRDR1c2, as well as those related to basal defense responses were up-regulated. Expression of genes associated with phytohormone signaling pathways were also altered, indicating that viroid infection initiates a complex array of changes in the host transcriptome. HSVd-g54 induced an earlier and stronger response than HSVd-h, and further examination of these differences will contribute to a better understanding of the mechanisms that determine viroid pathogenicity.