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

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.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

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.

Era of gapless plant genomes: innovations in sequencing and mapping technologies revolutionize genomics and breeding

Whole-genome sequencing and assembly have revolutionized plant genetics and molecular biology over the last two decades. However, significant shortcomings in first- and second-generation technology resulted in imperfect reference genomes: numerous and large gaps of low quality or undeterminable sequence in areas of highly repetitive DNA along with limited chromosomal phasing restricted the ability of researchers to characterize regulatory noncoding elements and genic regions that underwent recent duplication events. Recently, advances in long-read sequencing have resulted in the first gapless, telomere-to-telomere (T2T) assemblies of plant genomes. This leap forward has the potential to increase the speed and confidence of genomics and molecular experimentation while reducing costs for the research community.

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.

Convergent gene clusters underpin hyperforin biosynthesis in St John's wort

The meroterpenoid hyperforin is responsible for the antidepressant activity of St John's wort extracts, but the genes controlling its biosynthesis are unknown. Using genome mining and biochemical work, we characterize two biosynthetic gene clusters (BGCs) that encode the first three steps in the biosynthesis of hyperforin precursors. The findings of syntenic and phylogenetic analyses reveal the parallel assembly of the two BGCs. The syntenous BGC in Mesua ferrea indicates that the first cluster was assembled before the divergence of the Hypericaceae and Calophyllaceae families. The assembly of the second cluster is the result of a coalescence of genomic fragments after a major duplication event. The differences between the two BGCs - in terms of gene expression, response to methyl jasmonate, substrate specificity and subcellular localization of key enzymes - suggest that the presence of the two clusters could serve to generate separate pools of precursors. The parallel assembly of two BGCs with similar compositions in a single plant species is uncommon, and our work provides insights into how and when these gene clusters form. Our discovery helps to advance our understanding of the evolution of plant specialized metabolism and its genomic organization. Additionally, our results offer a foundation from which hyperforin biosynthesis can be more fully understood, and which can be used in future metabolic engineering applications.

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.

Large structural variations in the haplotype-resolved African cassava genome

Cassava (Manihot esculenta Crantz, 2n = 36) is a global food security crop. It has a highly heterozygous genome, high genetic load, and genotype-dependent asynchronous flowering. It is typically propagated by stem cuttings and any genetic variation between haplotypes, including large structural variations, is preserved by such clonal propagation. Traditional genome assembly approaches generate a collapsed haplotype representation of the genome. In highly heterozygous plants, this results in artifacts and an oversimplification of heterozygous regions. We used a combination of Pacific Biosciences (PacBio), Illumina, and Hi-C to resolve each haplotype of the genome of a farmer-preferred cassava line, TME7 (Oko-iyawo). PacBio reads were assembled using the FALCON suite. Phase switch errors were corrected using FALCON-Phase and Hi-C read data. The ultralong-range information from Hi-C sequencing was also used for scaffolding. Comparison of the two phases revealed >5000 large haplotype-specific structural variants affecting over 8 Mb, including insertions and deletions spanning thousands of base pairs. The potential of these variants to affect allele-specific expression was further explored. RNA-sequencing data from 11 different tissue types were mapped against the scaffolded haploid assembly and gene expression data are incorporated into our existing easy-to-use web-based interface to facilitate use by the broader plant science community. These two assemblies provide an excellent means to study the effects of heterozygosity, haplotype-specific structural variation, gene hemizygosity, and allele-specific gene expression contributing to important agricultural traits and further our understanding of the genetics and domestication of cassava.

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.

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.

Widely assumed phenotypic associations in Cannabis sativa lack a shared genetic basis

The flowering plant Cannabis sativa, cultivated for centuries for multiple purposes, displays extensive variation in phenotypic traits in addition to its wide array of secondary metabolite production. Notably, Cannabis produces two well-known secondary-metabolite cannabinoids: cannabidiolic acid (CBDA) and delta-9-tetrahydrocannabinolic acid (THCA), which are the main products sought by consumers in the medical and recreational market. Cannabis has several suggested subspecies which have been shown to differ in chemistry, branching patterns, leaf morphology and other traits. In this study we obtained measurements related to phytochemistry, reproductive traits, growth architecture, and leaf morphology from 297 hybrid individuals from a cross between two diverse lineages. We explored correlations among these characteristics to inform our understanding of which traits may be causally associated. Many of the traits widely assumed to be strongly correlated did not show any relationship in this hybrid population. The current taxonomy and legal regulation within Cannabis is based on phenotypic and chemical characteristics. However, we find these traits are not associated when lineages are inter-crossed, which is a common breeding practice and forms the basis of most modern marijuana and hemp germplasms. Our results suggest naming conventions based on leaf morphology do not correspond to the chemical properties in plants with hybrid ancestry. Therefore, a new system for identifying variation within Cannabis is warranted that will provide reliable identifiers of the properties important for recreational and, especially, medical use.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.

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.