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

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

Novel Identification of the Collection of Pathogenic Fungal Species Verticillium with the Development of Species-Specific SSR Markers.

The genus Verticillium is a group of ascomycete fungi that includes several pathogenic plant species. In 2011, a new taxonomic classification, proposed by Inderbitzin and coworkers (2011), re-defined the genus as Verticillium sensu stricto. The objective of our study was the re-classification of the fungal species held in the culture collection in the Slovenian Institute of Hop Research and Brewing in accordance with the newly established taxonomy. With the PCR marker system proposed by Inderbitzin and coworkers in 2011, we re-classified 88 Verticillium isolates out of the 105 samples that are held in the institute's bank, which were obtained from different geographic locations in Europe, North America, and Japan, and from different host plants, including alfalfa, cotton, hop, olive, potato, and tomato. However, the PCR marker for the V. dahliae identification proved to be less specific, and it resulted in the positive amplification of Gibellulopsis nigrescens, V. isaacii, and V. longisporum. To enable the accurate distinction of the fungi, the SSR and LAMP markers were added to the analyses. The 12 newly identified SSR markers, which were used in simplex PCR reactions or in combination, enabled the accurate identification of all included Verticillium isolates and could potentially be used as biomarkers for rapid and easy species identification.

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.

Small RNA Sequencing and Multiplex RT-PCR for Diagnostics of Grapevine Viruses and Virus-like Organisms.

Metagenomic approaches used for virus diagnostics allow for rapid and accurate detection of all viral pathogens in the plants. In order to investigate the occurrence of viruses and virus-like organisms infecting grapevine from the Ampelographic collection Kromberk in Slovenia, we used Ion Torrent small RNA sequencing (sRNA-seq) and the VirusDetect pipeline to analyze the sRNA-seq data. The used method revealed the presence of: Grapevine leafroll-associated virus 1 (GLRaV-1), Grapevine leafroll-associated virus 2 (GLRaV-2), Grapevine leafroll-associated virus 3 (GLRaV-3), Grapevine rupestris stem pitting-associated virus (GRSPaV), Grapevine fanleaf virus (GFLV) and its satellite RNA (satGFLV), Grapevine fleck virus (GFkV), Grapevine rupestris vein feathering virus (GRVFV), Grapevine Pinot gris virus (GPGV), Grapevine satellite virus (GV-Sat), Hop stunt viroid (HSVd), and Grapevine yellow speckle viroid 1 (GYSVd-1). Multiplex reverse transcription-polymerase chain reaction (mRT-PCR) was developed for validation of sRNA-seq predicted infections, including various combinations of viruses or viroids and satellite RNA. mRT-PCR could further be used for rapid and cost-effective routine molecular diagnosis, including widespread, emerging, and seemingly rare viruses, as well as viroids which testing is usually overlooked.

Elimination of Eight Viruses and Two Viroids from Preclonal Candidates of Six Grapevine Varieties (Vitis vinifera L.) through In Vivo Thermotherapy and In Vitro Meristem Tip Micrografting.

Viruses and virus-like organisms are a major problem in viticulture worldwide. They cannot be controlled by standard plant protection measures, and once infected, plants remain infected throughout their life; therefore, the propagation of healthy vegetative material is crucial. In vivo thermotherapy at 36-38 °C for at least six weeks, followed by meristem tip micrografting (0.1-0.2 mm) onto in vitro-growing seedling rootstocks of Vialla (Vitis labrusca × Vitis riparia), was successfully used to eliminate eight viruses (grapevine rupestris stem pitting-associated virus (GRSPaV), grapevine Pinot gris virus (GPGV), grapevine fanleaf virus (GFLV), grapevine leafroll-associated virus 3 (GLRaV-3), grapevine fleck virus (GFkV), grapevine rupestris vein feathering virus (GRVFV), grapevine Syrah virus-1 (GSyV-1), and raspberry bushy dwarf virus (RBDV)), as well as two viroids (hop stunt viroid (HSVd) and grapevine yellow speckle viroid 1 (GYSVd-1)) from preclonal candidates of six grapevine varieties (Vitis vinifera L.). A half-strength MS medium including vitamins supplemented with 30 g/L of sucrose and solidified with 8 g/L of agar, without plant growth regulators, was used for the growth and root development of micrografts and the subsequently micropropagated plants; no callus formation, hyperhydricity, or necrosis of shoot tips was observed. Although the overall regeneration was low (higher in white than in red varieties), a 100% elimination was achieved for all eight viruses, whereas the elimination level for viroids was lower, reaching only 39.2% of HSVd-free and 42.6% GYSVd-1-free vines. To the best of our knowledge, this is the first report of GPGV, GRVFV, GSyV-1, HSVd, and GYSVd-1 elimination through combining in vivo thermotherapy and in vitro meristem tip micrografting, and the first report of RBDV elimination from grapevines. The virus-free vines were successfully acclimatized in rockwool plugs and then transferred to soil.

Virome Status of Preclonal Candidates of Grapevine Varieties (Vitis vinifera L.) From the Slovenian Wine-Growing Region Primorska as Determined by High-Throughput Sequencing.

Diseases caused by viruses and virus-like organisms are one of the major problems in viticulture and grapevine marketing worldwide. Therefore, rapid and accurate diagnosis and identification is crucial. In this study, we used HTS of virus- and viroid-derived small RNAs to determine the virome status of Slovenian preclonal candidates of autochthonous and local grapevine varieties (Vitis vinifera L.). The method applied to the studied vines revealed the presence of nine viruses and two viroids. All viral entities were validated and more than 160 Sanger sequences were generated and deposited in NCBI. In addition, a complete description into the co-infections in each plant studied was obtained. No vine was found to be virus- and viroid-free, and no vine was found to be infected with only one virus or viroid, while the highest number of viral entities in a plant was eight.

Discovery of microRNA-like Small RNAs in Pathogenic Plant Fungus Verticillium nonalfalfae Using High-Throughput Sequencing and qPCR and RLM-RACE Validation.

Verticillium nonalfalfae (V. nonalfalfae) is one of the most problematic hop (Humulus lupulus L.) pathogens, as the highly virulent fungal pathotypes cause severe annual yield losses due to infections of entire hop fields. In recent years, the RNA interference (RNAi) mechanism has become one of the main areas of focus in plant-fungal pathogen interaction studies and has been implicated as one of the major contributors to fungal pathogenicity. MicroRNA-like RNAs (milRNAs) have been identified in several important plant pathogenic fungi; however, to date, no milRNA has been reported in the V. nonalfalfae species. In the present study, using a high-throughput sequencing approach and extensive bioinformatics analysis, a total of 156 milRNA precursors were identified in the annotated V. nonalfalfae genome, and 27 of these milRNA precursors were selected as true milRNA candidates, with appropriate microRNA hairpin secondary structures. The stem-loop RT-qPCR assay was used for milRNA validation; a total of nine V. nonalfalfae milRNAs were detected, and their expression was confirmed. The milRNA expression patterns, determined by the absolute quantification approach, imply that milRNAs play an important role in the pathogenicity of highly virulent V. nonalfalfae pathotypes. Computational analysis predicted milRNA targets in the V. nonalfalfae genome and in the host hop transcriptome, and the activity of milRNA-mediated RNAi target cleavage was subsequently confirmed for two selected endogenous fungal target gene models using the 5' RLM-RACE approach.

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.

Plant genera Cannabis and Humulus share the same pair of well-differentiated sex chromosomes.

We recently described, in Cannabis sativa, the oldest sex chromosome system documented so far in plants (12-28 Myr old). Based on the estimated age, we predicted that it should be shared by its sister genus Humulus, which is known also to possess XY chromosomes. Here, we used transcriptome sequencing of an F1 family of H. lupulus to identify and study the sex chromosomes in this species using the probabilistic method SEX-DETector. We identified 265 sex-linked genes in H. lupulus, which preferentially mapped to the C. sativa X chromosome. Using phylogenies of sex-linked genes, we showed that a region of the sex chromosomes had already stopped recombining in an ancestor of both species. Furthermore, as in C. sativa, Y-linked gene expression reduction is correlated to the position on the X chromosome, and highly Y degenerated genes showed dosage compensation. We report, for the first time in Angiosperms, a sex chromosome system that is shared by two different genera. Thus, recombination suppression started at least 21-25 Myr ago, and then (either gradually or step-wise) spread to a large part of the sex chromosomes (c. 70%), leading to a degenerated Y chromosome.

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.

Genetic Structure and Core Collection of Olive Germplasm from Albania Revealed by Microsatellite Markers.

Olive is considered one of the oldest and the most important cultivated fruit trees in Albania. In the present study, the genetic diversity and structure of Albanian olive germplasm is represented by a set of 194 olive genotypes collected in-situ in their natural ecosystems and in the ex-situ collection. The study was conducted using 26 microsatellite markers (14 genomic SSR and 12 Expressed Sequence Tag microsatellites). The identity analysis revealed 183 unique genotypes. Genetic distance-based and model-based Bayesian analyses were used to investigate the genetic diversity, relatedness, and the partitioning of the genetic variability among the Albanian olive germplasm. The genetic distance-based analysis grouped olives into 12 clusters, with an average similarity of 50.9%. Albanian native olives clustered in one main group separated from introduced foreign cultivars, which was also supported by Principal Coordinate Analysis (PCoA) and model-based methods. A core collection of 57 genotypes representing all allelic richness found in Albanian germplasm was developed for the first time. Herein, we report the first extended genetic characterization and structure of olive germplasm in Albania. The findings suggest that Albanian olive germplasm is a unique gene pool and provides an interesting genetic basis for breeding programs.

Population genetic analysis in old Montenegrin vineyards reveals ancient ways currently active to generate diversity in Vitis vinifera.

Global viticulture has evolved following market trends, causing loss of cultivar diversity and traditional practices. In Montenegro, modern viticulture co-exists with a traditional viticulture that still maintains ancient practices and exploits local cultivars. As a result, this region provides a unique opportunity to explore processes increasing genetic diversity. To evaluate the diversity of Montenegrin grapevines and the processes involved in their diversification, we collected and analyzed 419 samples in situ across the country (cultivated plants from old orchards and vines growing in the wild), and 57 local varieties preserved in a grapevine collection. We obtained 144 different genetic profiles, more than 100 corresponding to cultivated grapevines, representing a surprising diversity for one of the smallest European countries. Part of this high diversity reflects historical records indicating multiple and intense introduction events from diverse viticultural regions at different times. Another important gene pool includes many autochthonous varieties, some on the edge of extinction, linked in a complex parentage network where two varieties (Razaklija and Kratosija) played a leading role on the generation of indigenous varieties. Finally, analyses of genetic structure unveiled several putative proto-varieties, likely representing the first steps involved in the generation of new cultivars or even secondary domestication events.

The Potential of HTS Approaches for Accurate Genotyping in Grapevine (Vitis vinifera L.).

The main challenge associated with genotyping based on conventional length polymorphisms is the cross-laboratory standardization of allele sizes. This step requires the inclusion of standards and manual sizing to avoid false results. Capillary electrophoresis (CE) approaches limit the information to the length polymorphism and do not allow the determination of a complete marker sequence. As an alternative, high-throughput sequencing (HTS) offers complete information regarding marker sequences and their flanking regions. In this work, we investigated the suitability of a semi-quantitative sequencing approach for microsatellite genotyping using Illumina paired-end technology. Twelve microsatellite loci that are well established for grapevine CE typing were analysed on 96 grapevine samples from six different countries. We redesigned primers to the length of the amplicon for short sequencing (~100 bp). The primer pair was flanked with a 10 bp overhang for the introduction of barcodes on both sides of the amplicon to enable high multiplexing. The highest data peaks were determined as simple sequence repeat (SSR) alleles and compared with the CE dataset based on 12 reference samples. The comparison showed that HTS SSR genotyping can successfully replace the CE system in further experiments. We believe that, with next-generation sequencing, genotyping can be improved in terms of its speed, accuracy, and price.

Changes in the Phenolic Compounds of Hop (Humulus lupulus L.) Induced by Infection with Verticillium nonalfalfae, the Causal Agent of Hop Verticillium Wilt.

Phenolic compounds are involved in plant responses to various biotic and abiotic stress factors, with many studies suggesting their role in defense mechanisms against fungal pathogens. Soilborne vascular pathogen Verticillium nonalfalfae causes severe wilting and consequent dieback in a wide range of economically important crops, including hops (Humulus lupulus L.). In this study, we investigated the differential accumulation of phenolics in the susceptible "Celeia" and resistant "Wye Target" hop cultivars during the pathogenesis of Verticillium wilt. Quantitative polymerase chain reaction showed that colonization in the roots of both cultivars was intensive, but decreased continuously throughout the experiment in the resistant cultivar, while the relative fungal amount continuously increased in the stems of the susceptible cultivar. In response to colonization in the roots of the resistant cultivar, a significant increase in total flavanols was detected at three days postinoculation (dpi), suggesting a possible role in preventing fungus spread into the stems. The accumulation of phenolic compounds was less pronounced in the stems of the resistant cultivar since, compared to the latter, significant increases in flavonols at 3 and 15 dpi and hydroxycinnamic acids at 6 dpi were observed in the stems of the susceptible cultivar.

Proposal of a New Hybrid Breeding Method Based on Genotyping, Inter-Pollination, Phenotyping and Paternity Testing of Selected Elite F1 Hybrids.

Testing inbred lines for their combining ability is, due to high numbers of line to line testing needed for determination of hybrid performance, the most limiting factor in the F1 hybrid breeding procedure. We propose a novel method of F1 hybrid breeding that enables evaluation of large number of line to line crosses for their hybrid performance. Inbred lines (preferably doubled haploid - DH) are produced from heterozygous populations, genotyped and maintained. A group of lines is inter-pollinated randomly and their progeny examined. To identify elite F1 hybrids, these individual plants are selected by their superior phenotypic characteristics. Finally using paternity testing only of selected hybrids, the origin of paternal lines is revealed. To predict the number of F1 offspring needed in relation to the number of inbred lines being inter-pollinated, a mathematical formula was developed. For instance, using this formula for the inter-pollination of 60 distinct lines, the probability of obtaining all descendants of paternal-parent lines in a maternal-parent row represented at least once is achieved with 420 F1 plants in a row (p = 0.95). In a practical experiment with white cabbage, DH lines were produced using microspore culture; plants were grown to maturity and genotyped at eight polymorphic SSR loci. Two groups of lines (36 and 33 lines per group) were inter-pollinated by two methods, either using cage pollination with bumblebees or using open pollination in isolated field. A total of 9,858 F1 plants were planted and based on their phenotypic characteristics 213 were selected as elite phenotypes. 99 of them were genetically diverse and 5 of them were selected as super elite. Selected plants were analysed by the same SSR markers and the paternal origin of selected F1 plants was determined. Out of 213 selected elite plants 48 were reciprocals thus exhibiting power of selection based on single plant. We demonstrate that this new approach to hybrid development is efficient in white cabbage and we propose breeders to test it in various vegetable and crop species. Moreover, some other aspects of the proposed technique need to be tested and verified both for practical and economic criteria.

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.

RNA interference core components identified and characterised in Verticillium nonalfalfae, a vascular wilt pathogenic plant fungi of hops.

The conserved RNA interference mechanism (RNAi) in the fungal kingdom has become a focus of intense scientific investigation. The three catalytic core components, Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RdRP), and their associated small interfering RNA molecules (siRNAs) have been identified and characterised in several fungal species. Recent studies have proposed that RNAi is a major contributor to the virulence of fungal pathogens as a result of so-called trans-kingdom RNA silencing. In the present study, we report on the existence of three core RNAi proteins in the pathogenic plant fungus Verticillium nonalfalfae, which is a soilborne plant pathogen that causes severe wilting disease in hops (Humulus lupulus L.). Two DCL proteins, two AGO proteins, and two RdRP proteins were identified, and their conserved RNAi domains were characterised. Our phylogeny results confirm the existing taxonomic relationships in the Ascomycete fungal phylum and show that the fungi of the Hypocreomycetidae subclass of the Sordariomycetes class have high amino acid sequence similarity. The expression analysis revealed a potential role of RNAi in the pathogenicity of the fungi, since all the RNAi genes were highly upregulated in the highly virulent isolate T2 and were also differentially expressed in the V. nonalfalfae-susceptible Celeia and V. nonalfalfae-resistant Wye Target cultivars.

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

Viroids, the smallest known pathogens, unable to encode any proteins, can cause severe diseases in their host plants. One of the proposed mechanisms of their pathogenicity includes silencing the host's genes via viroid-derived small RNAs, which are products of the host's immune response to the viroid's double stranded RNA. Humulus lupulus (hop) plants are hosts to several viroids; two of them, HLVd and CBCVd, are interesting models for studying host-viroid interactions, due to the symptomless infection of the former and severe stunting disease caused by the latter. To study these interactions, we constructed a deep hop NGS transcriptome based on 35 Gb paired-end sequencing data assembled into over 74 Mb of contigs. These transcripts were used for in-silico prediction of target transcripts of vd-sRNA of the two aforementioned viroids, using two different software tools. Prediction models revealed that 1062 and 1387 hop transcripts share nucleotide similarities with HLVd- and CBCVd-derived small RNAs, respectively, so they could be silenced in an RNA interference process. Furthermore, we selected 17 transcripts from 4 groups of targets involved in the metabolism of plant hormones, small RNA biogenesis, transcripts with high complementarity with viroid-derived small RNAs and transcripts targeted by CBCVd-derived small RNAs with high cellular concentrations. Their expression was monitored by reverse transcription quantitative PCR performed using leaf, flower and cone samples. Additionally, the expression of 5 pathogenesis related genes was monitored. Expression analysis confirmed high expression levels of four pathogenesis related genes in leaves of HLVd and CBCVd infected hop plants. Expression fluctuations were observed for the majority of targets, with possible evidence of downregulation of GATA transcription factor by CBCVd- and of linoleate 13S-lipoxygenase by HLVd-derived small RNAs. These results provide a deep transcriptome of hop and the first insights into complex viroid-hop plant interactions.

Comparative transcriptional analysis of hop responses to infection with Verticillium nonalfalfae.

KEY MESSAGE: Dynamic transcriptome profiling revealed excessive, yet ineffective, immune response to V. nonalfalfae infection in susceptible hop, global gene downregulation in shoots of resistant hop and only a few infection-associated genes in roots. Hop (Humulus lupulus L.) production is hampered by Verticillium wilt, a disease predominantly caused by the soil-borne fungus Verticillium nonalfalfae. Only a few hop cultivars exhibit resistance towards it and mechanisms of this resistance have not been discovered. In this study, we compared global transcriptional responses in roots and shoots of resistant and susceptible hop plants infected by a lethal strain of V. nonalfalfae. Time-series differential gene expression profiles between infected and mock inoculated plants were determined and subjected to network-based analysis of functional enrichment. In the resistant hop cultivar, a remarkably low number of genes were differentially expressed in roots in response to V. nonalfalfae infection, while the majority of differentially expressed genes were down-regulated in shoots. The most significantly affected genes were related to cutin biosynthesis, cell wall biogenesis, lateral root development and terpenoid biosynthesis. On the other hand, susceptible hop exhibited a strong defence response in shoots and roots, including increased expression of genes associated with plant responses, such as innate immunity, wounding, jasmonic acid pathway and chitinase activity. Strong induction of defence-associated genes in susceptible hop and a low number of infection-responsive genes in the roots of resistant hop are consistent with previous findings, confirming the pattern of excessive response of the susceptible cultivar, which ultimately fails to protect the plant from V. nonalfalfae. This research offers a multifaceted overview of transcriptional responses of susceptible and resistant hop cultivars to V. nonalfalfae infection and represents a valuable resource in the study of this plant-pathogen interaction.

Identification of Novel Virulence-Associated Proteins Secreted to Xylem by Verticillium nonalfalfae During Colonization of Hop Plants.

Plant pathogens employ various secreted proteins to suppress host immunity for their successful host colonization. Identification and characterization of pathogen-secreted proteins can contribute to an understanding of the pathogenicity mechanism and help in disease control. We used proteomics to search for proteins secreted to xylem by the vascular pathogen Verticillium nonalfalfae during colonization of hop plants. Three highly abundant fungal proteins were identified: two enzymes, α-N-arabinofuranosidase (VnaAbf4.216) and peroxidase (VnaPRX1.1277), and one small secreted hypothetical protein (VnaSSP4.2). These are the first secreted proteins so far identified in xylem sap following infection with Verticillium spp. VnaPRX1.1277, classified as a heme-containing peroxidase from Class II, similar to other Verticillium spp. lignin-degrading peroxidases, and VnaSSP4.2, a 14-kDa cysteine-containing protein with unknown function and with a close homolog in related V. alfalfae strains, were further examined. The in planta expression of VnaPRX1.1277 and VnaSSP4.2 genes increased with the progression of colonization, implicating their role in fungal virulence. Indeed, V. nonalfalfae deletion mutants of both genes exhibited attenuated virulence on hop plants, which returned to the level of the wild-type pathogenicity in the knockout complementation lines, supporting VnaPRX1.1277 and VnaSSP4.2 as virulence factors required to promote V. nonalfalfae colonization of hop plants.