A method for targeting a specified segment of DNA to a bacterial microorganelle.

Encapsulation of a selected DNA molecule in a cell has important implications for bionanotechnology. Non-viral proteins that can be used as nucleic acid containers include proteinaceous subcellular bacterial microcompartments (MCPs) that self-assemble into a selectively permeable protein shell containing an enzymatic core. Here, we adapted a propanediol utilization (Pdu) MCP into a synthetic protein cage to package a specified DNA segment in vivo, thereby enabling subsequent affinity purification. To this end, we engineered the LacI transcription repressor to be routed, together with target DNA, into the lumen of a Strep-tagged Pdu shell. Sequencing of extracted DNA from the affinity-isolated MCPs shows that our strategy results in packaging of a DNA segment carrying multiple LacI binding sites, but not the flanking regions. Furthermore, we used LacI to drive the encapsulation of a DNA segment containing operators for LacI and for a second transcription factor.

First report of apple hammerhead viroid infecting apple trees in Montenegro.

Apple hammerhead viroid (AHVd, Pelamoviroid, Avsunviroidae) is one of the five viroids infecting apples. It has been identified on all continents except Australia since its viroid nature was confirmed (DiSerio et al. 2018; CABI and EPPO 2022). AHVd has been found in apple trees showing leaf mosaic, ringspot and dieback (Hamdi et al., 2021). Apple (Malus domestica Borkh.) and its wild relatives are traditionally grown in Montenegro. With an annual production of 7767 tons on 216 ha, it is the second most important fruit tree (after plum) in the country (Anonymous 2022). In a 2020-2022 survey, 29 apple trees exhibiting virus-like symptoms (e.g. mosaic, necrosis) were sampled throughout Montenegro, including 16 locations in eight municipalities (Podgorica, Danilovgrad, Niksic, Mojkovac, Bijelo Polje, Berane, Pljevlja and Savnik). Small RNAs were isolated using the mirVana miRNA Isolation Kit (Ambion, Life Technologies) and pooled into three bulk samples. Each bulk contained 9 to 10 samples. Libraries of sRNAs were constructed using the Ion Total RNA-Seq Kit v2 and barcoded using the Xpress RNA-Seq Barcode 1-16 Kit (Ion Torrent) according to the manufacturer's instructions. Small RNA library sequencing was performed on Illumina platform (Novogene Europe) yielding 9.9, 9.8 and 18.6 million reads in the three libraries. The CLC Genomics Workbench software was used to demultiplex the reads into pools using the 'Demultiplex Reads' tool. The online program VirusDetect (Zheng et al. 2017) was used for virus/viroid detection and identification. Besides viruses known to infect apple (apple stem grooving virus, apple stem pitting virus, apple mosaic virus), contigs mapping to AHVd were identified in all three bulks enabling full AHVd genomes reconstruction. To verify AHVd presence, all 29 apple samples were tested by reverse transcription-polymerase chain reaction (RT-PCR) using the AHVd PG13f/PG12r primers (Messmer et al. 2017). AHVd amplicons were obtained in three samples (30/21, 32/21 and 38/21) from bulk 1 and two samples (47/21 and 55/21) from bulk 2, while all samples from bulk 3 tested negative potentially due to the low titer of the pathogen or nucleotide mismatches at the 3' end of the primers. The three amplicons from bulk 1 were Sanger sequenced and partial AHVd genomes over 200 nts were obtained from two of them (30/21 and 32/21) (GenBank acc. nos. OQ863319 and OR020603). Furthermore, three full consensus AHVd genomes were assembled in Geneious Prime by mapping Sanger sequences onto contigs from Virus Detect and named 30/21, 32/21 and 38/21 (acc. nos. PP133245, -46, and -47, respectively). All three genomes exhibited conserved hammerhead motifs (Messmer et al. 2017). In BLASTn analysis, the isolate 30/21 from Montenegro shared the highest nt identity (98.8%) with the isolate SA-36 (ON564299) from Czechia, while 32/21 and 38/21 showed the highest identities (95.4% and 92.3%) with isolates SD17_2-3 (MK188691) from Canada and JF2 (ON564298) from Czechia, respectively. To the best of our knowledge, this is the first report of AHVd infecting Malus domestica in Montenegro. The AHVd-positive samples 30/21 and 32/21 originated from at least two-decade-old apple trees from Niksic, whilst 38/21 came from a 40-year-old tree from Mojkovac district, suggesting that this viroid has long been present in different parts of the country. The AHVd discovery in Montenegro should be considered in any phytosanitary regulations and pome fruit certification program in the country.

Dig-up Primers: A Pipeline for Identification of Polymorphic Microsatellites Loci within Assemblies of Related Species.

Simple sequence repeats (SSRs) have become one of the most popular molecular markers and are used in numerous fields, including conservation genetics, population genetic studies, and genetic mapping. Advances in next-generation sequencing technology and the growing amount of genomic data are driving the development of bioinformatics tools for SSR marker design. These tools work with different combinations of input data, which can be raw reads or assemblies, and with one or more input datasets. We present here a new strategy and implementation of a simple standalone pipeline that utilizes more than one assembly for the in silico design of PCR primers for microsatellite loci in more than one species. Primers are tested in silico to determine if they are polymorphic, eliminating the need to test time-consuming cross-species amplification in the laboratory. The end result is a set of markers that are in silico polymorphic in all analyzed species and have great potential for the identification of interspecies hybrids. The efficiency of the tool is demonstrated using two examples at different taxonomic levels and with different numbers of input assemblies to generate promising, high-quality SSR markers.

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

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

Towards the Investigation of the Adaptive Divergence in a Species of Exceptional Ecological Plasticity: Chromosome-Scale Genome Assembly of Chouardia litardierei (Hyacinthaceae).

One of the central goals of evolutionary biology is to understand the genomic basis of adaptive divergence. Different aspects of evolutionary processes should be studied through genome-wide approaches, therefore maximizing the investigated genomic space. However, in-depth genome-scale analyses often are restricted to a model or economically important species and their closely related wild congeners with available reference genomes. Here, we present the high-quality chromosome-level genome assembly of Chouardia litardierei, a plant species with exceptional ecological plasticity. By combining PacBio and Hi-C sequencing technologies, we generated a 3.7 Gbp genome with a scaffold N50 size of 210 Mbp. Over 80% of the genome comprised repetitive elements, among which the LTR retrotransposons prevailed. Approximately 86% of the 27,257 predicted genes were functionally annotated using public databases. For the comparative analysis of different ecotypes' genomes, the whole-genome sequencing of two individuals, each from a distinct ecotype, was performed. The detected above-average SNP density within coding regions suggests increased adaptive divergence-related mutation rates, therefore confirming the assumed divergence processes within the group. The constructed genome presents an invaluable resource for future research activities oriented toward the investigation of the genetics underlying the adaptive divergence that is likely unfolding among the studied species' ecotypes.

Identification and characterization of long non-coding RNA and their response against citrus bark cracking viroid infection in Humulus lupulus.

Long non-coding RNAs (lncRNAs) are a highly heterogeneous class of non-protein-encoding transcripts that play an essential regulatory role in diverse biological processes, including stress responses. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) poses a major threat to the production of Humulus lupulus (hop) plants. In this study, we systematically investigate the characteristics of the lncRNAs in hop and their role in CBCVd-infection using RNA-sequencing data. Following a stringent filtration criterion, a total of 3598 putative lncRNAs were identified with a high degree of certainty, of which 19% (684) of the lncRNAs were significantly differentially expressed (DE) in CBCVd-infected hop, which were predicted to be mainly involved in plant-pathogen interactions, kinase cascades, secondary metabolism and phytohormone signal transduction. Besides, several lncRNAs and CBCVd-responsive lncRNAs were identified as the precursor of microRNAs and predicted as endogenous target mimics (eTMs) for hop microRNAs involved in CBCVd-infection.

Chloroplast Genome Annotation Tools: Prolegomena to the Identification of Inverted Repeats.

The development of next-generation sequencing technology and the increasing amount of sequencing data have brought the bioinformatic tools used in genome assembly into focus. The final step of the process is genome annotation, which works on assembled genome sequences to identify the location of genome features. In the case of organelle genomes, specialized annotation tools are used to identify organelle genes and structural features. Numerous annotation tools target chloroplast sequences. Most chloroplast DNA genomes have a quadripartite structure caused by two copies of a large inverted repeat. We investigated the strategies of six annotation tools (Chloë, Chloroplot, GeSeq, ORG.Annotate, PGA, Plann) for identifying inverted repeats and analyzed their success using publicly available complete chloroplast sequences of taxa belonging to the asterid and rosid clades. The annotation tools use two different approaches to identify inverted repeats, using existing general search tools or implementing stand-alone solutions. The chloroplast sequences studied show that there are different types of imperfections in the assembled data and that each tool performs better on some sequences than the others.

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

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

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

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

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

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

Integrated Proteo-Transcriptomic Analyses Reveal Insights into Regulation of Pollen Development Stages and Dynamics of Cellular Response to Apple Fruit Crinkle Viroid (AFCVd)-Infection in Nicotiana tabacum.

Tobacco (Nicotiana tabacum) pollen is a well-suited model for studying many fundamental biological processes owing to its well-defined and distinct development stages. It is also one of the major agents involved in the transmission of infectious viroids, which is the primary mechanism of viroid pathogenicity in plants. However, some viroids are non-transmissible and may be possibly degraded or eliminated during the gradual process of pollen development maturation. The molecular details behind the response of developing pollen against the apple fruit crinkle viroid (AFCVd) infection and viroid eradication is largely unknown. In this study, we performed an integrative analysis of the transcriptome and proteome profiles to disentangle the molecular cascade of events governing the three pollen development stages: early bicellular pollen (stage 3, S3), late bicellular pollen (stage 5, S5), and 6 h-pollen tube (PT6). The integrated analysis delivered the molecular portraits of the developing pollen against AFCVd infection, including mechanistic insights into the viroid eradication during the last steps of pollen development. The isobaric tags for label-free relative quantification (iTRAQ) with digital gene expression (DGE) experiments led us to reliably identify subsets of 5321, 5286, and 6923 proteins and 64,033, 60,597, and 46,640 expressed genes in S3, S5, and PT6, respectively. In these subsets, 2234, 2108 proteins and 9207 and 14,065 mRNAs were differentially expressed in pairwise comparisons of three stages S5 vs. S3 and PT6 vs. S5 of control pollen in tobacco. Correlation analysis between the abundance of differentially expressed mRNAs (DEGs) and differentially expressed proteins (DEPs) in pairwise comparisons of three stages of pollen revealed numerous discordant changes in mRNA/protein pairs. Only a modest correlation was observed, indicative of divergent transcription, and its regulation and importance of post-transcriptional events in the determination of the fate of early and late pollen development in tobacco. The functional and enrichment analysis of correlated DEGs/DEPs revealed the activation in pathways involved in carbohydrate metabolism, amino acid metabolism, lipid metabolism, and cofactor as well as vitamin metabolism, which points to the importance of these metabolic pathways in pollen development. Furthermore, the detailed picture of AFCVd-infected correlated DEGs/DEPs was obtained in pairwise comparisons of three stages of infected pollen. The AFCVd infection caused the modulation of several genes involved in protein degradation, nuclear transport, phytohormone signaling, defense response, and phosphorylation. Intriguingly, we also identified several factors including, DNA-dependent RNA-polymerase, ribosomal protein, Argonaute (AGO) proteins, nucleotide binding proteins, and RNA exonucleases, which may plausibly involve in viroid stabilization and eradication during the last steps of pollen development. The present study provides essential insights into the transcriptional and translational dynamics of tobacco pollen, which further strengthens our understanding of plant-viroid interactions and support for future mechanistic studies directed at delineating the functional role of candidate factors involved in viroid elimination.

Mapping the Gene Expression Spectrum of Mediator Subunits in Response to Viroid Infection in Plants.

The mediator (MED) represents a large, conserved, multi-subunit protein complex that regulates gene expression through interactions with RNA polymerase II and enhancer-bound transcription factors. Expanding research accomplishments suggest the predominant role of plant MED subunits in the regulation of various physiological and developmental processes, including the biotic stress response against bacterial and fungal pathogens. However, the involvement of MED subunits in virus/viroid pathogenesis remains elusive. In this study, we investigated for the first time the gene expression modulation of selected MED subunits in response to five viroid species (Apple fruit crinkle viroid (AFCVd), Citrus bark cracking viroid (CBCVd), Hop latent viroid (HLVd), Hop stunt viroid (HSVd), and Potato spindle tuber viroid (PSTVd)) in two model plant species (Nicotiana tabacum and N. benthamiana) and a commercially important hop (Humulus lupulus) cultivar. Our results showed a differential expression pattern of MED subunits in response to a viroid infection. The individual plant MED subunits displayed a differential and tailored expression pattern in response to different viroid species, suggesting that the MED expression is viroid- and plant species-dependent. The explicit evidence obtained from our results warrants further investigation into the association of the MED subunit with symptom development. Together, we provide a comprehensive portrait of MED subunit expression in response to viroid infection and a plausible involvement of MED subunits in fine-tuning transcriptional reprogramming in response to viroid infection, suggesting them as a potential candidate for rewiring the defense response network in plants against pathogens.

Global transcriptome profiling and functional analysis reveal that tissue-specific constitutive overexpression of cytochrome P450s confers tolerance to imidacloprid in palm weevils in date palm fields.

BACKGROUND: Cytochrome P450-dependent monooxygenases (P450s), constituting one of the largest and oldest gene superfamilies found in many organisms from bacteria to humans, play a vital role in the detoxification and inactivation of endogenous toxic compounds. The use of various insecticides has increased over the last two decades, and insects have developed resistance to most of these compounds through the detoxifying function of P450s. In this study, we focused on the red palm weevil (RPW), Rhynchophorus ferrugineus, the most devastating pest of palm trees worldwide, and demonstrated through functional analysis that upregulation of P450 gene expression has evolved as an adaptation to insecticide stress arising from exposure to the neonicotinoid-class systematic insecticide imidacloprid. RESULTS: Based on the RPW global transcriptome analysis, we identified 101 putative P450 genes, including 77 likely encoding protein coding genes with ubiquitous expression. A phylogenetic analysis revealed extensive functional and species-specific diversification of RPW P450s, indicating that multiple CYPs actively participated in the detoxification process. We identified highly conserved paralogs of insect P450s that likely play a role in the development of resistance to imidacloprid: Drosophila Cyp6g1 (CYP6345J1) and Bemisia tabaci CYP4C64 (CYP4LE1). We performed a toxicity bioassay and evaluated the induction of P450s, followed by the identification of overexpressed P450s, including CYP9Z82, CYP6fra5, CYP6NR1, CYP6345J1 and CYP4BD4, which confer cross-resistance to imidacloprid. In addition, under imidacloprid insecticide stress in a date palm field, we observed increased expression of various P450 genes, with CYP9Z82, CYP4BD4, CYP6NR1 and CYP6345J1 being the most upregulated detoxification genes in RPWs. Expression profiling and cluster analysis revealed P450 genes with multiple patterns of induction and differential expression. Furthermore, we used RNA interference to knock down the overexpressed P450s, after which a toxicity bioassay and quantitative expression analysis revealed likely candidates involved in metabolic resistance against imidacloprid in RPW. Ingestion of double-stranded RNA (dsRNA) successfully knocked down the expression of CYP9Z82, CYP6NR1 and CYP345J1 and demonstrated that silencing of CYP345J1 and CYP6NR1 significantly decreased the survival rate of adult RPWs treated with imidacloprid, indicating that overexpression of these two P450s may play an important role in developing tolerance to imidacloprid in a date palm field. CONCLUSION: Our study provides useful background information on imidacloprid-specific induction and overexpression of P450s, which may enable the development of diagnostic tools/markers for monitoring the spread of insecticide resistant RPWs. The observed trend of increasing tolerance to imidacloprid in the date palm field therefore indicated that strategies for resistance management are urgently needed.

Evaluation of Disease Severity and Global Transcriptome Response Induced by Citrus bark cracking viroid, Hop latent viroid, and Their Co-Infection in Hop (Humulus lupulus L.).

Viroids are small non-capsidated, single-stranded, covalently-closed circular noncoding RNA replicons of 239-401 nucleotides that exploit host factors for their replication, and some cause disease in several economically important crop plants, while others appear to be benign. The proposed mechanisms of viroid pathogenesis include direct interaction of the genomic viroid RNA with host factors and post-transcriptional or transcriptional gene silencing via viroid-derived small RNAs (vd-sRNAs) generated by the host defensive machinery. Humulus lupulus (hop) plants are hosts to several viroids among which Hop latent viroid (HLVd) and Citrus bark cracking viroid (CBCVd) are attractive model systems for the study of viroid-host interactions due to the symptomless infection of the former and severe symptoms induced by the latter in this indicator host. To better understand their interactions with hop plant, a comparative transcriptomic analysis based on RNA sequencing (RNA-seq) was performed to reveal the transcriptional alterations induced as a result of single HLVd and CBCVd infection in hop. Additionally, the effect of HLVd on the aggressiveness of CBCVd that underlies severe stunting in hop in a mixed infection was studied by transcriptomic analysis. Our analysis revealed that CBCVd infection resulted in dynamic changes in the activity of genes as compared to single HLVd infection and their mixed infection. The differentially expressed genes that are involved in defense, phytohormone signaling, photosynthesis and chloroplasts, RNA regulation, processing and binding; protein metabolism and modification; and other mechanisms were more modulated in the CBCVd infection of hop. Nevertheless, Gene Ontology (GO) classification and pathway enrichment analysis showed that the expression of genes involved in the proteolysis mechanism is more active in a mixed infection as compared to a single one, suggesting co-infecting viroids may result in interference with host factors more prominently. Collectively, our results provide a deep transcriptome of hop and insight into complex single HLVd, CBCVd, and their coinfection in hop-plant interactions.

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.

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.

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

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

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

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

Identification and characterization of microRNAs in Humulus lupulus using high-throughput sequencing and their response to Citrus bark cracking viroid (CBCVd) infection.

BACKGROUND: Hop (Humulus lupulus L.) plants are grown primarily for the brewing industry and have been used as a traditional medicinal herb for a long time. Severe hop stunt disease caused by the recently discovered Citrus bark cracking viroid (CBCVd) is one of the most devastating diseases among other viroid infections in hop. MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in gene expression regulation. To identify miRNAs in hop and their response to CBCVd-infection, two small RNA (sRNA) libraries were prepared from healthy and CBCVd-infected hop plants and were investigated by high throughput sequencing. RESULTS: A total of 67 conserved and 49 novel miRNAs were identified. Among them, 36 conserved and 37 novel miRNAs were found to be differentially recovered in response to CBCVd-infection. A total of 311 potential targets was predicted for conserved and novel miRNAs based on a sequence homology search using hop transcriptome data. The majority of predicted targets significantly belonged to transcriptional factors that may regulate hop leaf, root and cone growth and development. In addition, the identified miRNAs might also play an important roles in other cellular and metabolic processes, such as signal transduction, stress response and other physiological processes, including prenylflavonoid biosynthesis pathways. Quantitative real time PCR analysis of selected targets revealed their negative correlation with their corresponding CBCVd-responsive miRNAs. CONCLUSIONS: Based on the results, we concluded that CBCVd-responsive miRNAs modulate several hormone pathways and transcriptional factors that play important roles in the regulation of metabolism, growth and development. These results provide a framework for further analysis of regulatory roles of sRNAs in plant defense mechanism including other hop infecting viroids in particular.

Identification of the genes involved in odorant reception and detection in the palm weevil Rhynchophorus ferrugineus, an important quarantine pest, by antennal transcriptome analysis.

BACKGROUND: The Red Palm Weevil (RPW) Rhynchophorus ferrugineus (Oliver) is one of the most damaging invasive insect species in the world. This weevil is highly specialized to thrive in adverse desert climates, and it causes major economic losses due to its effects on palm trees around the world. RPWs locate palm trees by means of plant volatile cues and use an aggregation pheromone to coordinate a mass-attack. Here we report on the high throughput sequencing of the RPW antennal transcriptome and present a description of the highly expressed chemosensory gene families. RESULTS: Deep sequencing and assembly of the RPW antennal transcriptome yielded 35,667 transcripts with an average length of 857 bp and identified a large number of highly expressed transcripts of odorant binding proteins (OBPs), chemosensory proteins (CSPs), odorant receptors/co-receptors (ORs/Orcos), sensory neuron membrane proteins (SNMPs), gustatory receptors (GRs) and ionotropic receptors (IRs). In total, 38 OBPs, 12 CSPs, 76 ORs, 1 Orco, 6 SNMPs, 15 GRs and 10 IRs were annotated in the R. ferrugineus antennal transcriptome. A comparative transcriptome analysis with the bark beetle showed that 25% of the blast hits were unique to R. ferrugineus, indicating a higher, more complete transcript coverage for R. ferrugineus. We categorized the RPW ORs into seven subfamilies of coleopteran ORs and predicted two new subfamilies of ORs. The OR protein sequences were compared with those of the flour beetle, the cerambycid beetle and the bark beetle, and we identified coleopteran-specific, highly conserved ORs as well as unique ORs that are putatively involved in RPW aggregation pheromone detection. We identified 26 Minus-C OBPs and 8 Plus-C OBPs and grouped R. ferrugineus OBPs into different OBP-subfamilies according to phylogeny, which indicated significant species-specific expansion and divergence in R. ferrugineus. We also identified a diverse family of CSP proteins, as well as a coleopteran-specific CSP lineage that diverged from Diptera and Lepidoptera. We identified several extremely diverged IR orthologues as well as highly conserved insect IR co-receptor orthologous transcripts in R. ferrugineus. Notably, GR orthologous transcripts for CO2-sensing and sweet tastants were identified in R. ferrugineus, and we found a great diversity of GRs within the coleopteran family. With respect to SNMP-1 and SNMP-2 orthologous transcripts, one SNMP-1 orthologue was found to be strikingly highly expressed in the R. ferrugineus antennal transcriptome. CONCLUSION: Our study presents the first comprehensive catalogue of olfactory gene families involved in pheromone and general odorant detection in R. ferrugineus, which are potential novel targets for pest control strategies.