The Draft Genome of Hop (Humulus lupulus), an Essence for Brewing

Unveiling the phenology and associated floral regulatory pathways of Humulus lupulus L. in subtropical conditions

MAIN CONCLUSION

The hop phenological cycle was described in subtropical condition of Brazil showing that flowering can happen at any time of year and this was related to developmental molecular pathways. Hops are traditionally produced in temperate regions, as it was believed that vernalization was necessary for flowering. Nevertheless, recent studies have revealed the potential for hops to flower in tropical and subtropical climates. In this work, we observed that hops in the subtropical climate of Minas Gerais, Brazil grow and flower multiple times throughout the year, independently of the season, contrasting with what happens in temperate regions. This could be due to the photoperiod consistently being inductive, with daylight hours below the described threshold (16.5 h critical). We observed that when the plants reached 7-9 nodes, the leaves began to transition from heart-shaped to trilobed-shaped, which could be indicative of the juvenile to adult transition. This could be related to the fact that the 5th node (in plants with 10 nodes) had the highest expression of miR156, while two miR172s increased in the 20th node (in plants with 25 nodes). Hop flowers appeared later, in the 25th or 28th nodes, and the expression of HlFT3 and HlFT5 was upregulated in plants between 15 and 20 nodes, while the expression of HlTFL3 was upregulated in plants with 20 nodes. These results indicate the role of axillary meristem age in regulating this process and suggest that the florigenic signal should be maintained until the hop plants bloom. In addition, it is possible that the expression of TFL is not sufficient to inhibit flowering in these conditions and promote branching. These findings suggest that the reproductive transition in hop under inductive photoperiodic conditions could occur in plants between 15 and 20 nodes. Our study sheds light on the intricate molecular mechanisms underlying hop floral development, paving the way for potential advancements in hop production on a global scale.

Diversity and Pathobiology of an Ilarvirus Unexpectedly Detected in Diverse Plants and Global Sequencing Data

High-throughput sequencing (HTS) and sequence mining tools revolutionized virus detection and discovery in recent years, and implementing them with classical plant virology techniques results in a powerful approach to characterize viruses. An example of a virus discovered through HTS is Solanum nigrum ilarvirus 1 (SnIV1) (Bromoviridae), which was recently reported in various solanaceous plants from France, Slovenia, Greece, and South Africa. It was likewise detected in grapevines (Vitaceae) and several Fabaceae and Rosaceae plant species. Such a diverse set of source organisms is atypical for ilarviruses, thus warranting further investigation. In this study, modern and classical virological tools were combined to accelerate the characterization of SnIV1. Through HTS-based virome surveys, mining of sequence read archive datasets, and a literature search, SnIV1 was further identified from diverse plant and non-plant sources globally. SnIV1 isolates showed relatively low variability compared with other phylogenetically related ilarviruses. Phylogenetic analyses showed a distinct basal clade of isolates from Europe, whereas the rest formed clades of mixed geographic origin. Furthermore, systemic infection of SnIV1 in Solanum villosum and its mechanical and graft transmissibility to solanaceous species were demonstrated. Near-identical SnIV1 genomes from the inoculum (S. villosum) and inoculated Nicotiana benthamiana were sequenced, thus partially fulfilling Koch's postulates. SnIV1 was shown to be seed-transmitted and potentially pollen-borne, has spherical virions, and possibly induces histopathological changes in infected N. benthamiana leaf tissues. Overall, this study provides information to better understand the diversity, global presence, and pathobiology of SnIV1; however, its possible emergence as a destructive pathogen remains uncertain. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

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.

An improved assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and refines time of divergence estimates for the Cannabaceae family

We present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly were fragmented. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur. The chromosome-level assembly of Cascade was developed by scaffolding the previously reported Cascade assembly generated with PacBio long-read sequencing and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assessed molecular evolution in gene sequences, gene family expansion and contraction, and time of divergence from Cannabis sativa and other closely related plant species using Bayesian inference. We identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the estimate for the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families. The Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.

In Vitro Propagation of Humulus lupulus through the Induction of Axillary Bud Development

Humulus lupulus (hop) is a necessary material for beer brewing. Improved breeding cultivars of hops with enhanced tolerance to environmental stresses, such as drought and heat stress, accompanying climate change have been developed. However, a propagation system, which is needed for the proliferation of new cultivars, is not currently available for hops. In this study, we found that treatment of stem explants with 0.01–0.05 ppm gibberellic acid (GA3) induced the development of axillary buds in the hop cultivar Kirin-2, resulting in the proliferation of shoot branching. Additionally, 0.01 ppm benzyl adenine (BA) enhanced the development of axillary buds formed in response to 0.05 ppm GA3 in various hop cultivars, particularly Nugget. The development of axillary buds was strongly repressed by the application of 0.05 ppm BA at a concentration equal to the 0.05 ppm GA3 concentration, which showed the possibility that a high concentration of cytokinin preferentially prevents the effect of GA3 on the development of axillary buds in hops. These results indicated that combined treatment of stem explants with GA3 and cytokinin at appropriate concentrations is effective for the propagation of proliferated hop cultivars through shoot branching.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Key Enzymes Involved in the Synthesis of Hops Phytochemical Compounds: From Structure, Functions to Applications

Humulus lupulus L. is an essential source of aroma compounds, hop bitter acids, and xanthohumol derivatives mainly exploited as flavourings in beer brewing and with demonstrated potential for the treatment of certain diseases. To acquire a comprehensive understanding of the biosynthesis of these compounds, the primary enzymes involved in the three major pathways of hops' phytochemical composition are herein critically summarized. Hops' phytochemical components impart bitterness, aroma, and antioxidant activity to beers. The biosynthesis pathways have been extensively studied and enzymes play essential roles in the processes. Here, we introduced the enzymes involved in the biosynthesis of hop bitter acids, monoterpenes and xanthohumol derivatives, including the branched-chain aminotransferase (BCAT), branched-chain keto-acid dehydrogenase (BCKDH), carboxyl CoA ligase (CCL), valerophenone synthase (VPS), prenyltransferase (PT), 1-deoxyxylulose-5-phosphate synthase (DXS), 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (HDR), Geranyl diphosphate synthase (GPPS), monoterpene synthase enzymes (MTS), cinnamate 4-hydroxylase (C4H), chalcone synthase (CHS_H1), chalcone isomerase (CHI)-like proteins (CHIL), and O-methyltransferase (OMT1). Furthermore, research advancements of each enzyme in terms of reaction conditions, substrate recognition, enzyme structures, and use in engineered microbes are described in depth. Hence, an extensive review of the key enzymes involved in the phytochemical compounds of hops will provide fundamentals for their applications in beer production.

Large-scale whole-genome resequencing unravels the domestication history of Cannabis sativa

Cannabis sativa has long been an important source of fiber extracted from hemp and both medicinal and recreational drugs based on cannabinoid compounds. Here, we investigated its poorly known domestication history using whole-genome resequencing of 110 accessions from worldwide origins. We show that C. sativa was first domesticated in early Neolithic times in East Asia and that all current hemp and drug cultivars diverged from an ancestral gene pool currently represented by feral plants and landraces in China. We identified candidate genes associated with traits differentiating hemp and drug cultivars, including branching pattern and cellulose/lignin biosynthesis. We also found evidence for loss of function of genes involved in the synthesis of the two major biochemically competing cannabinoids during selection for increased fiber production or psychoactive properties. Our results provide a unique global view of the domestication of C. sativa and offer valuable genomic resources for ongoing functional and molecular breeding research.

The hop downy mildew pathogen Pseudoperonospora humuli

Pseudoperonospora humuli is an obligate biotrophic oomycete that causes downy mildew, one of the most devastating diseases of cultivated hop, Humulus lupulus. Downy mildew occurs in all production areas of the crop in the Northern Hemisphere and Argentina. The pathogen overwinters in hop crowns and roots, and causes considerable crop loss. Downy mildew is managed by sanitation practices, planting of resistant cultivars, and fungicide applications. However, the scarcity of sources of host resistance and fungicide resistance in pathogen populations complicates disease management. This review summarizes the current knowledge on the symptoms of the disease, life cycle, virulence factors, and management of hop downy mildew, including various forecasting systems available in the world. Additionally, recent developments in genomics and effector discovery, and the future prospects of using such resources in successful disease management are also discussed.

TAXONOMY

Class: Oomycota; Order: Peronosporales; Family: Peronosporaceae; Genus: Pseudoperonospora; Species: Pseudoperonospora humuli.

DISEASE SYMPTOMS

The disease is characterized by systemically infected chlorotic shoots called "spikes". Leaf symptoms and signs include angular chlorotic lesions and profuse sporulation on the abaxial side of the leaf. Under severe disease pressure, dark brown discolouration or lesions are observed on cones. Infected crowns have brown to black streaks when cut open. Cultivars highly susceptible to crown rot may die at this phase of the disease cycle without producing shoots. However, foliar symptoms may not be present on plants with systemically infected root systems.

INFECTION PROCESS

Pathogen mycelium overwinters in buds and crowns, and emerges on infected shoots in spring. Profuse sporulation occurs on infected tissues and sporangia are released and dispersed by air currents. Under favourable conditions, sporangia germinate and produce biflagellate zoospores that infect healthy tissue, thus perpetuating the infection cycle. Though oospores are produced in infected tissues, their role in the infection cycle is not defined.

CONTROL

Downy mildew on hop is managed by a combination of sanitation practices and timely fungicide applications. Forecasting systems are used to time fungicide applications for successful management of the disease. USEFUL WEBSITES: https://content.ces.ncsu.edu/hop-downy-mildew (North Carolina State University disease factsheet), https://www.canr.msu.edu/resources/michigan-hop-management-guide (Michigan Hop Management Guide), http://uspest.org/risk/models (Oregon State University Integrated Plant Protection Center degree-day model for hop downy mildew), https://www.usahops.org/cabinet/data/Field-Guide.pdf (Field Guide for Integrated Pest Management in Hops).

Tracing foreign sequences in plant transcriptomes and genomes using OCT4, a POU domain protein

Contaminations in sequencing data, especially in reference genomes, lead to inevitable errors in downstream analyses. Similarly, presence of contaminants in transcriptomes, misrepresents the molecular basis of various interactions. In this study, we report the presence of a large number of plant transcriptomes contaminated with RNAs encoding POU domain proteins; a family of proteins that has not been reported in plants and fungi. Besides, our findings illustrated that there are four POU domain protein-coding sequences in the reference genome of Rhodamnia argentea. It turned out that the existing foreign fragments are related to arthropods that are considered as plant pests. We also identified two contaminated draft genomes, Humulus lupulus and Cannabis sativa that contained complete rDNA sequences originating from Tetranychus species. As a result, careful screening of sequencing data before releasing them in public databases or checking existing genomes for possible contaminations is recommended.

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

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

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

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

A draft phased assembly of the diploid Cascade hop (Humulus lupulus) genome

Hop (Humulus lupulus L. var Lupulus) is a diploid, dioecious plant with a history of cultivation spanning more than one thousand years. Hop cones are valued for their use in brewing and contain compounds of therapeutic interest including xanthohumol. Efforts to determine how biochemical pathways responsible for desirable traits are regulated have been challenged by the large (2.8 Gb), repetitive, and heterozygous genome of hop. We present a draft haplotype-phased assembly of the Cascade cultivar genome. Our draft assembly and annotation of the Cascade genome is the most extensive representation of the hop genome to date. PacBio long-read sequences from hop were assembled with FALCON and partially phased with FALCON-Unzip. Comparative analysis of haplotype sequences provides insight into selective pressures that have driven evolution in hop. We discovered genes with greater sequence divergence enriched for stress-response, growth, and flowering functions in the draft phased assembly. With improved resolution of long terminal retrotransposons (LTRs) due to long-read sequencing, we found that hop is over 70% repetitive. We identified a homolog of cannabidiolic acid synthase (CBDAS) that is expressed in multiple tissues. The approaches we developed to analyze the draft phased assembly serve to deepen our understanding of the genomic landscape of hop and may have broader applicability to the study of other large, complex genomes.

The Effector Repertoire of the Hop Downy Mildew Pathogen Pseudoperonospora humuli

Pseudoperonospora humuli is an obligate biotrophic oomycete that causes downy mildew (DM), one of the most destructive diseases of cultivated hop that can lead to 100% crop loss in susceptible cultivars. We used the published genome of P. humuli to predict the secretome and effectorome and analyze the transcriptome variation among diverse isolates and during infection of hop leaves. Mining the predicted coding genes of the sequenced isolate OR502AA of P. humuli revealed a secretome of 1,250 genes. We identified 296 RXLR and RXLR-like effector-encoding genes in the secretome. Among the predicted RXLRs, there were several WY-motif-containing effectors that lacked canonical RXLR domains. Transcriptome analysis of sporangia from 12 different isolates collected from various hop cultivars revealed 754 secreted proteins and 201 RXLR effectors that showed transcript evidence across all isolates with reads per kilobase million (RPKM) values > 0. RNA-seq analysis of OR502AA-infected hop leaf samples at different time points after infection revealed highly expressed effectors that may play a relevant role in pathogenicity. Quantitative RT-PCR analysis confirmed the differential expression of selected effectors. We identified a set of P. humuli core effectors that showed transcript evidence in all tested isolates and elevated expression during infection. These effectors are ideal candidates for functional analysis and effector-assisted breeding to develop DM resistant hop cultivars.

Identification of tandem repeat families from long-read sequences of Humulus lupulus

Hop (Humulus lupulus L.) is known for its use as a bittering agent in beer and has a rich history of cultivation, beginning in Europe and now spanning the globe. There are five wild varieties worldwide, which may have been introgressed with cultivated varieties. As a dioecious species, its obligate outcrossing, non-Mendelian inheritance, and genomic structural variability have confounded directed breeding efforts. Consequently, understanding the hop genome represents a considerable challenge, requiring additional resources. In order to facilitate investigations into the transmission genetics of hop, we report here a tandem repeat discovery pipeline developed using k-mer filtering and dot plot analysis of PacBio long-read sequences from the hop cultivar Apollo. From this we identified 17 new and distinct tandem repeat sequence families, which represent candidates for FISH probe development. For two of these candidates, HuluTR120 and HuluTR225, we produced oligonucleotide FISH probes from conserved regions of and demonstrated their utility by staining meiotic chromosomes from wild hop, var. neomexicanus to address, for example, questions about hop transmission genetics. Collectively, these tandem repeat sequence families represent new resources suitable for development of additional cytogenomic tools for hop research.

Population Diversity and Structure of Podosphaera macularis in the Pacific Northwestern United States and Other Populations

Powdery mildew, caused by Podosphaera macularis, is one of the most important diseases of hop. The disease was first reported in the Pacific Northwestern United States, the primary hop-growing region in this country, in the mid-1990s. More recently, the disease has reemerged in newly planted hopyards of the eastern United States, as hop production has expanded to meet demands of local craft brewers. The spread of strains virulent on previously resistant cultivars, the paucity of available fungicides, and the potential introduction of the MAT1-2 mating type to the western United States, all threaten sustainability of hop production. We sequenced the transcriptome of 104 isolates of P. macularis collected throughout the western United States, eastern United States, and Europe to quantify genetic diversity of pathogen populations and elucidate the possible origins of pathogen populations in the western United States. Discriminant analysis of principal components grouped isolates within three to five geographic populations, dependent on stringency of grouping criteria. Isolates from the western United States were phenotyped and categorized into one of three pathogenic races based on disease symptoms generated on differential cultivars. Western U.S. populations were clonal, irrespective of pathogenic race, and grouped with isolates originating from Europe. Isolates originating from wild hop plants in the eastern United States were genetically differentiated from all other populations, whereas isolates from cultivated hop plants in the eastern United States mostly grouped with isolates originating from the west, consistent with origins from nursery sources. Mating types of isolates originating from cultivated western and eastern U.S. hop plants were entirely MAT1-1. In contrast, a 1:1 ratio of MAT1-1 and MAT1-2 was observed with isolates sampled from wild plants or Europe. Within the western United States a set of highly differentiated loci were identified in P. macularis isolates associated with virulence to the powdery mildew R-gene R6. The weight of genetic and phenotypic evidence suggests a European origin of the P. macularis populations in the western United States, followed by spread of the pathogen from the western United States to re-emergent production regions in the eastern United States. Furthermore, R6 compatibility appears to have been selected from an extant isolate within the western United States. Greater emphasis on sanitation measures during propagation and quarantine policies should be considered to limit further spread of novel genotypes of the pathogen, both between and within production areas.

The Genomics of Cannabis and Its Close Relatives

Cannabis sativa L. is an important yet controversial plant with a long history of recreational, medicinal, industrial, and agricultural use, and together with its sister genus Humulus, it represents a group of plants with a myriad of academic, agricultural, pharmaceutical, industrial, and social interests. We have performed a meta-analysis of pooled published genomics data, andwe present a comprehensive literature review on the evolutionary history of Cannabis and Humulus, including medicinal and industrial applications. We demonstrate that current Cannabis genome assemblies are incomplete, with ∼10% missing, 10-25% unmapped, and 45S and 5S ribosomal DNA clusters as well as centromeres/satellite sequences not represented. These assemblies are also ordered at a low resolution, and their consensus quality clouds the accurate annotation of complete, partial, and pseudogenized gene copies. Considering the importance of genomics in the development of any crop, this analysis underlines the need for a coordinated effort to quantify the genetic and biochemical diversity of this species.

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.

Phytochemical Characterization of Wild Hops (Humulus lupulus ssp. lupuloides) Germplasm Resources From the Maritimes Region of Canada

A survey was conducted in the Maritimes region of eastern Canada to measure the phytochemical diversity of prenylchalcone, soft resins (alpha & beta acids), and flavonol constituents from 30 unique wild-growing populations of hops (Humulus lupulus L.). Based on cone chemometrics, the majority of accessions (63.3%) are native Humulus lupulus ssp. lupoloides, with cones containing both xanthogalenol and 4'-O-methyl xanthohumol as chemotaxonomic indicator molecules. Interestingly, the leaves of all verified Humulus lupulus ssp. lupulus accessions accumulated high proportions (>0.20 total flavonols) of two acylated flavonol derivatives (kaempferol-3-O-(6''-O-malonyl)-β-D-glucopyranoside; quercetin-3-O-(6''-O-malonyl)-β-D-glucopyranoside), both previously unreported from hops leaves. The native lupuloides accessions examined possess only trace amounts of this compound in their leaves (<0.10 total flavonols), suggesting its potential utility as a novel, leaf-derived chemotaxonomic marker for subspecies identification purposes. A leaf-derived taxonomic marker is useful for identifying wild-growing accessions, as leaves are present throughout the entire growing season, whereas cones are only produced late in summer. Additionally, the collection of cones from 10-meter tall wild plants in overgrown riparian habitats is often difficult. The total levels of alpha acids, beta acids, and prenylchalcones in wild-collected Maritimes lupuloides cones are markedly higher than those previously reported for lupuloides individuals in the westernmost extent of its native range and show potentially valuable traits for future cultivar development, while some may be worthy of immediate commercial release. The accessions will be maintained as a core germplasm resource for future cultivar development.

Draft genome of Tanacetum cinerariifolium, the natural source of mosquito coil

Pyrethrum (Tanacetum cinerariifolium), which is a perennial Asteraceae plant with white daisy-like flowers, is the original source of mosquito coils and is known for the biosynthesis of the pyrethrin class of natural insecticides. However, the molecular basis of the production of pyrethrins by T. cinerariifolium has yet to be fully elucidated. Here, we present the 7.1-Gb draft genome of T. cinerariifolium, consisting of 2,016,451 scaffolds and 60,080 genes predicted with high confidence. Notably, analyses of transposable elements (TEs) indicated that TEs occupy 33.84% of the genome sequence. Furthermore, TEs of the sire and oryco clades were found to be enriched in the T. cinerariifolium-specific evolutionary lineage, occupying a total of 13% of the genome sequence, a proportion approximately 8-fold higher than that in other plants. InterProScan analysis demonstrated that biodefense-related toxic proteins (e.g., ribosome inactivating proteins), signal transduction-related proteins (e.g., histidine kinases), and metabolic enzymes (e.g., lipoxygenases, acyl-CoA dehydrogenases/oxygenases, and P450s) are also highly enriched in the T. cinerariifolium genome. Molecular phylogenetic analysis detected a variety of enzymes with genus-specific multiplication, including both common enzymes and others that appear to be specific to pyrethrin biosynthesis. Together, these data identify possible novel components of the pyrethrin biosynthesis pathway and provide new insights into the unique genomic features of T. cinerariifolium.

The spatio-temporal expression of some genes involved in the biosynthetic pathways of terpenes/phenylpropanoids in yarrow (Achillea millefolium)

Yarrow (Achillea millefolium) is a medicinal plant from the Asteracea which biosynthesize different secondary metabolites especially terpenes and phenylpropanoids. To improve our understanding of the regulatory mechanisms behind the biosynthesis of these compounds we analyzed the expression of some genes associated with the biosynthesis of terpenes and phenylpropanoids in different tissues and in response to trans-cinnamic acid (tCA) as an inhibitor of PAL activity. Isolation and expression analysis of DXR, GPPS, PAL and CHS genes together with linalool synthase (LIS) as monoterpene synthase was conducted in different developmental stages of leaves, flowers and in response to trans-cinnamic acid (tCA). Differential expression of these genes observed in different tissues. tCA up-regulated the biosynthetic genes of monterpenes and down-regulated the biosynthetic genes of phenylpropanoids. Gene expression analysis in intact leaves and leaves without glandular trichomes showed that DXR, LIS, PAL and CHS are highly expressed in glandular trichomes while GPPS expressed ubiquitously. Analysis of essential oils composition showed that sesquiterpenes and monoterpenes are main compounds; in which from 57 identified compounds the highest were germacreneD (% 11.5), guaiol (%10.38), spatulenol (%8.73) and caryophyllene oxide (%7.48).

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.

Genome Sequencing and Transcriptome Analysis of the Hop Downy Mildew Pathogen Pseudoperonospora humuli Reveal Species-Specific Genes for Molecular Detection

Pseudoperonospora humuli is an obligate oomycete pathogen of hop (Humulus lupulus) that causes downy mildew, an important disease in most production regions in the Northern Hemisphere. The pathogen can cause a systemic infection in hop, overwinter in the root system, and infect propagation material. Substantial yield loss may occur owing to P. humuli infection of strobiles (seed cones), shoots, and cone-bearing branches. Fungicide application and cultural practices are the primary methods to manage hop downy mildew. However, effective, sustainable, and cost-effective management of downy mildew can be improved by developing early detection systems to inform on disease risk and timely fungicide application. However, no species-specific diagnostic assays or genomic resources are available for P. humuli. The genome of the P. humuli OR502AA isolate was partially sequenced using Illumina technology and assembled with ABySS. The assembly had a minimum scaffold length of 500 bp and an N50 (median scaffold length of the assembled genome) of 19.2 kbp. A total number of 18,656 genes were identified using MAKER standard gene predictions. Additionally, transcriptome assemblies were generated using RNA-seq and Trinity for seven additional P. humuli isolates. Bioinformatics analyses of next generation sequencing reads of P. humuli and P. cubensis (a closely related sister species) identified 242 candidate species-specific P. humuli genes that could be used as diagnostic molecular markers. These candidate genes were validated using polymerase chain reaction against a diverse collection of isolates from P. humuli, P. cubensis, and other oomycetes. Overall, four diagnostic markers were found to be uniquely present in P. humuli. These candidate markers identified through comparative genomics can be used for pathogen diagnostics in propagation material, such as rhizomes and vegetative cuttings, or adapted for biosurveillance of airborne sporangia, an important source of inoculum in hop downy mildew epidemics.

De novo assembly of the Platycladus orientalis (L.) Franco transcriptome provides insight into the development and pollination mechanism of female cone based on RNA-Seq data

For seed-bearing plants, the basis of seed and fruit formation is pollination. The normal progression of pollination is through advances in continuous signal exchange and material transfer, which occur mainly in female reproductive organs; thus, the molecular mechanism of development in female reproductive organs is vital for understanding the principle of pollination. However, molecular biology studies on the development of female cones related to pollination are rare and unclear in gymnosperms, especially in Cupressaceae. In this study, Platycladus orientalis, a monotypic genus within Cupressaceae, was chosen to examine female cone transcriptomes at pre-pollination and pollination stages by Illumina paired-end sequencing technology to de novo sequence six libraries with 3 biological replicates. These libraries were used to construct a P. orientalis transcriptome database containing 71,669 unigenes (4,963 upregulated unigenes and 11,747 downregulated unigenes at the pollination stage) for subsequent analysis. Based on the annotations and expression levels, the functions of differentially expressed unigenes and enriched pathways between the developmental processes of female cones were analysed to detail the preliminary development and pollination mechanism of the female cone. Targeted investigations were specifically performed to determine the elementary mechanism of secretion and functioning of the pollination drop, a vital ovule secretion at the pollination stage. Ultimately, the expression of 15 unigenes selected between two stages were further assessed and confirmed using qRT-PCR, which demonstrated reliable data and significant differences in the expression profiles of key genes. As one of the largest available transcriptomic resources of this species, the database is constructed to prospectively adapt to the physiological and genomic data of woody plants. This work provided the first transcriptome profile of P. orientalis female cones at different developmental stages, and will promote the illumination of the pollination mechanism of P. orientalis, and will serve as the basis for in-depth genomic study in the Cupressaceae family. This initiative will arouse the interest and attention of scholars and pave the way for future studies.

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.

The molecular biology of fruity and floral aromas in beer and other alcoholic beverages

Aroma compounds provide attractiveness and variety to alcoholic beverages. We discuss the molecular biology of a major subset of beer aroma volatiles, fruity and floral compounds, originating from raw materials (malt and hops), or formed by yeast during fermentation. We introduce aroma perception, describe the most aroma-active, fruity and floral compounds in fruits and their presence and origin in beer. They are classified into categories based on their functional groups and biosynthesis pathways: (1) higher alcohols and esters, (2) polyfunctional thiols, (3) lactones and furanones, and (4) terpenoids. Yeast and hops are the main sources of fruity and flowery aroma compounds in beer. For yeast, the focus is on higher alcohols and esters, and particularly the complex regulation of the alcohol acetyl transferase ATF1 gene. We discuss the release of polyfunctional thiols and monoterpenoids from cysteine- and glutathione-S-conjugated compounds and glucosides, respectively, the primary biological functions of the yeast enzymes involved, their mode of action and mechanisms of regulation that control aroma compound production. Furthermore, we discuss biochemistry and genetics of terpenoid production and formation of non-volatile precursors in Humulus lupulus (hops). Insight in these pathways provides a toolbox for creating innovative products with a diversity of pleasant aromas.

Single Nucleotide Polymorphisms and Biochemical Markers As Complementary Tools To Characterize Hops ( Humulus lupulus L.) in Brewing Practice

In brewing practice, the use of the appropriate hop variety is essential to produce consistent and high-quality beers. Yet, hop batches of the same variety cultivated in different geographical regions can display significant biochemical differences, resulting in specific taste- and aroma-related characteristics in beer. In this study, we illustrate the complementarity of genetic and biochemical fingerprinting methods to fully characterize hop batches. Using genotyping-by-sequencing (GBS), a set of 1 830 polymorphic single nucleotide polymorphism (SNP) markers generated 48 unique genetic fingerprints for a collection of 56 commercial hop varieties. Three groups of varieties, consisting of somaclonal variants, could not be further differentiated using this set of markers. Biochemical marker information offered added value to characterize hop samples from a given variety grown at different geographical locations. We demonstrate the power of combining genetic and biochemical fingerprints for quality control of hop batches in the brewing industry.

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

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.

Comparative genomics of the nonlegume Parasponia reveals insights into evolution of nitrogen-fixing rhizobium symbioses

Nodules harboring nitrogen-fixing rhizobia are a well-known trait of legumes, but nodules also occur in other plant lineages, with rhizobia or the actinomycete Frankia as microsymbiont. It is generally assumed that nodulation evolved independently multiple times. However, molecular-genetic support for this hypothesis is lacking, as the genetic changes underlying nodule evolution remain elusive. We conducted genetic and comparative genomics studies by using Parasponia species (Cannabaceae), the only nonlegumes that can establish nitrogen-fixing nodules with rhizobium. Intergeneric crosses between Parasponia andersonii and its nonnodulating relative Trema tomentosa demonstrated that nodule organogenesis, but not intracellular infection, is a dominant genetic trait. Comparative transcriptomics of P. andersonii and the legume Medicago truncatula revealed utilization of at least 290 orthologous symbiosis genes in nodules. Among these are key genes that, in legumes, are essential for nodulation, including NODULE INCEPTION (NIN) and RHIZOBIUM-DIRECTED POLAR GROWTH (RPG). Comparative analysis of genomes from three Parasponia species and related nonnodulating plant species show evidence of parallel loss in nonnodulating species of putative orthologs of NIN, RPG, and NOD FACTOR PERCEPTION Parallel loss of these symbiosis genes indicates that these nonnodulating lineages lost the potential to nodulate. Taken together, our results challenge the view that nodulation evolved in parallel and raises the possibility that nodulation originated ∼100 Mya in a common ancestor of all nodulating plant species, but was subsequently lost in many descendant lineages. This will have profound implications for translational approaches aimed at engineering nitrogen-fixing nodules in crop plants.

Evidence for contemporary plant mitoviruses

Mitoviruses have small RNA(+) genomes, replicate in mitochondria, and have been shown to infect only fungi to date. For this report, sequences that appear to represent nearly complete plant mitovirus genomes were recovered from publicly available transcriptome data. Twenty of the refined sequences, 2684-2898 nt long and derived from 10 different species of land plants, appear to encompass the complete coding regions of contemporary plant mitoviruses, which furthermore constitute a monophyletic cluster within genus Mitovirus. Complete coding sequences of several of these viruses were recovered from multiple transcriptome (but not genome) studies of the same plant species and also from multiple plant tissues. Crop plants among implicated hosts include beet and hemp. Other new results suggest that such genuine plant mitoviruses were immediate ancestors to endogenized mitovirus elements now widespread in land plant genomes. Whether these mitoviruses are wholly cryptic with regard to plant health remains to be investigated.

Diversity and evolution of the repetitive genomic content in Cannabis sativa

Background

The repetitive content of the genome, once considered to be “junk DNA”, is in fact an essential component of genomic architecture and evolution. In this study, we used the genomes of three varieties of Cannabis sativa, three varieties of Humulus lupulus and one genotype of Morus notabilis to explore their repetitive content using a graph-based clustering method, designed to explore and compare repeat content in genomes that have not been fully assembled.

Results

The repetitive content in the C. sativa genome is mainly composed of the retrotransposons LTR/Copia and LTR/Gypsy (14% and 14.8%, respectively), ribosomal DNA (2%), and low-complexity sequences (29%). We observed a recent copy number expansion in some transposable element families. Simple repeats and low complexity regions of the genome show higher intra and inter species variation.

Conclusions

As with other sequenced genomes, the repetitive content of C. sativa’s genome exhibits a wide range of evolutionary patterns. Some repeat types have patterns of diversity consistent with expansions followed by losses in copy number, while others may have expanded more slowly and reached a steady state. Still, other repetitive sequences, particularly ribosomal DNA (rDNA), show signs of concerted evolution playing a major role in homogenizing sequence variation.

Electronic supplementary material

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

Non-Mendelian Single-Nucleotide Polymorphism Inheritance and Atypical Meiotic Configurations are Prevalent in Hop

Hop ( L.) breeding programs seek to exploit genetic resources for bitter flavor, aroma, and disease resistance. However, these efforts have been thwarted by segregation distortion including female-biased sex ratios. To better understand the transmission genetics of hop, we genotyped 4512 worldwide accessions of hop, including cultivars, landraces, and over 100 wild accessions using a genotyping-by-sequencing (GBS) approach. From the resulting ∼1.2 million single-nucleotide polymorphisms (SNPs), prequalified GBS markers were validated by inferences in population structures and phylogeny. Analysis of pseudo-testcross (Pt) mapping data from F families revealed mixed patterns of Mendelian and non-Mendelian segregation. Three-dimensional (3D) cytogenetic analysis of late meiotic prophase nuclei from two wild and two cultivated hop revealed conspicuous and prevalent occurrences of multiple, atypical, nondisomic chromosome complexes including autosomes. We used genome-wide association studies (GWAS) and fixation index (F) analysis to demonstrate selection mapping of genetic loci for key traits including sex, bitter acids, and drought tolerance. Among the possible mechanisms underlying the observed segregation distortion from the genomic data analysis, the cytogenetic analysis points to meiotic chromosome behavior as one of the contributing factors. The findings shed light on long-standing questions on the unusual transmission genetics and phenotypic variation in hop, with major implications for breeding, cultivation, and the natural history of .

Challenges towards Revitalizing Hemp: A Multifaceted Crop

Hemp has been an important crop throughout human history for food, fiber, and medicine. Despite significant progress made by the international research community, the basic biology of hemp plants remains insufficiently understood. Clear objectives are needed to guide future research. As a semi-domesticated plant, hemp has many desirable traits that require improvement, including eliminating seed shattering, enhancing the quantity and quality of stem fiber, and increasing the accumulation of phytocannabinoids. Methods to manipulate the sex of hemp plants will also be important for optimizing yields of seed, fiber, and cannabinoids. Currently, research into trait improvement is hindered by the lack of molecular techniques adapted to hemp. Here we review how addressing these limitations will help advance our knowledge of plant biology and enable us to fully domesticate and maximize the agronomic potential of this promising crop.

Heap: a highly sensitive and accurate SNP detection tool for low-coverage high-throughput sequencing data

Recent availability of large-scale genomic resources enables us to conduct so called genome-wide association studies (GWAS) and genomic prediction (GP) studies, particularly with next-generation sequencing (NGS) data. The effectiveness of GWAS and GP depends on not only their mathematical models, but the quality and quantity of variants employed in the analysis. In NGS single nucleotide polymorphism (SNP) calling, conventional tools ideally require more reads for higher SNP sensitivity and accuracy. In this study, we aimed to develop a tool, Heap, that enables robustly sensitive and accurate calling of SNPs, particularly with a low coverage NGS data, which must be aligned to the reference genome sequences in advance. To reduce false positive SNPs, Heap determines genotypes and calls SNPs at each site except for sites at the both ends of reads or containing a minor allele supported by only one read. Performance comparison with existing tools showed that Heap achieved the highest F-scores with low coverage (7X) restriction-site associated DNA sequencing reads of sorghum and rice individuals. This will facilitate cost-effective GWAS and GP studies in this NGS era. Code and documentation of Heap are freely available from https://github.com/meiji-bioinf/heap (29 March 2017, date last accessed) and our web site (http://bioinf.mind.meiji.ac.jp/lab/en/tools.html (29 March 2017, date last accessed)).

Hop stunt viroid: Effect on Host (Humulus lupulus) Transcriptome and Its Interactions With Hop Powdery Mildew (Podospheara macularis)

Viroids are the smallest known plant pathogens that exploit host systems for their replication and cause diseases in many hosts. In this study, the host response of hop plants to Hop stunt viroid (HSVd) infection was studied through transcriptome analysis. RNA sequence analysis of hop leaves infected with HSVd revealed dynamic changes in hop gene expression. Defense-related genes and genes involved in lipid and terpenoid metabolism are the major categories that showed differential expression due to HSVd infection. Additionally, the effect of HSVd on development of hop powdery mildew (Podospheara macularis) (HPM) was studied. Transcriptome analysis followed by quantitative reverse transcription-polymerase chain reaction analysis showed that transcript levels of pathogenesis-related (PR) genes such as PR protein 1, chitinase, and thaumatin-like protein genes are induced in leaves infected with HPM alone. The response in these genes to HPM is significantly down-regulated in leaves with HSVd-HPM mixed infection. These results confirm that HSVd alters host metabolism, physiology, and plant defense responses. Nevertheless, in detached leaf assays, HPM consistently expanded faster on HSVd-negative leaves relative to HSVd-positive leaves. Although HSVd infection suppresses elements associated with the host immunity response, infection by HSVd is antagonistic to HPM infection of hops.

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.

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.