Developmental stages and episode-specific regulatory genes in andromonoecious melon flower development.

BACKGROUND AND AIMS: Given the lack of specific studies on floral development in melon (Cucumis melo L.), we carried out an extensive study involving morphological and transcriptomic analyses to characterize floral development in this species. METHODS: Using an andromonoecious line, we analysed the development of floral buds in male and hermaphrodite flowers with both light microscopy and scanning electron microscopy. Based on flower lengths, we established a correlation between the developmental stages and four main episodes of floral development and conducted an extensive RNA sequencing analysis of these episodes. KEY RESULTS: We identified 12 stages of floral development, from the appearance of the floral meristems to anthesis. The main structural differences between male and hermaphrodite flowers appeared between stages 6 and 7; later stages of development leading to the formation of organs and structures in both types of flowers were also described. We analysed the gene expression patterns of the four episodes in flower development to find the genes that were specific to each given episode. Among others, we identified genes that defined the passage from one episode to the next according to the ABCDE model of floral development. CONCLUSIONS: This work combines a detailed morphological analysis and a comprehensive transcriptomic study to enable characterization of the structural and molecular mechanisms that determine the floral development of an andromonoecious genotype in melon. Taken together, our results provide a first insight into gene regulation networks in melon floral development that are crucial for flowering and pollen formation, highlighting potential targets for genetic manipulation to improve crop yield of melon in the future.

The tomato calcium-permeable channel 4.1 (SlOSCA4.1) is a susceptibility factor for pepino mosaic virus.

The hyperosmolality-gated calcium permeable channel 4.1 (OSCA4.1) belongs to an evolutionarily conserved small family of mechano-sensitive channels. OSCA members may represent key players in plant resistance to drought and to pathogen infection but are scarcely studied. After screening for resistance to pepino mosaic virus (PepMV) a collection of 1000 mutagenized tomato families, we identified a mutant showing no symptoms and reduced virus accumulation. Resistance was mapped to chromosome 2 between positions 46 309 531 to 47 044 163, where a missense mutation caused the putative truncation of the OSCA4.1 protein. A CRISPR/Cas9 slosca4.1 mutant was resistant to PepMV, but not to tobacco mosaic virus or potato virus X. Inoculation of mutant and wild type tomato protoplasts showed that resistance was expressed in single cells, suggesting a role for SlOSCA4.1 in early viral function(s); congruently, SlOSCA4.1 re-localized to structures reminiscent of viral replication complexes. We propose that SlOSCA4.1 contributes to the correct regulation of the Ca2+ homeostasis necessary for optimal PepMV infection. PepMV is a pandemic virus that causes significant losses in tomato crops worldwide. In spite of its importance, no tomato-resistant varieties have been deployed yet; the mutant identified here has great potential to breed tomato varieties resistant to PepMV.

Tomato SlGSTU38 interacts with the PepMV coat protein and promotes viral infection.

Pepino mosaic virus (PepMV) is pandemic in tomato crops, causing important economic losses world-wide. No PepMV-resistant varieties have been developed yet. Identification of host factors interacting with PepMV proteins is a promising source of genetic targets to develop PepMV-resistant varieties. The interaction between the PepMV coat protein (CP) and the tomato glutathione S-transferase (GST) SlGSTU38 was identified in a yeast two-hybrid (Y2H) screening and validated by directed Y2H and co-immunoprecipitation assays. SlGSTU38-knocked-out Micro-Tom plants (gstu38) generated by the CRISPR/Cas9 technology together with live-cell imaging were used to understand the role of SlGSTU38 during infection. The transcriptomes of healthy and PepMV-infected wild-type (WT) and gstu38 plants were profiled by RNA-seq analysis. SlGSTU38 functions as a PepMV-specific susceptibility factor in a cell-autonomous manner and relocalizes to the virus replication complexes during infection. Besides, knocking out SlGSTU38 triggers reactive oxygen species accumulation in leaves and the deregulation of stress-responsive genes. SlGSTU38 may play a dual role: On the one hand, SlGSTU38 may exert a proviral function depending on its specific interaction with the PepMV CP; and on the other hand, SlGSTU38 may delay PepMV-infection sensing by participating in the redox intracellular homeostasis in a nonspecific manner.

Cohombrillo-associated virus: a novel virus infecting Ecballium elaterium plants.

We determined the complete genome sequence of a new virus infecting Ecballium elaterium ('cohombrillo amargo') plants, a weed species common on the borders of cultivated fields in the Mediterranean region. The genome of this virus is composed of two molecules of monocistronic positive-sense RNA, 6,934 and 3,501 nucleotides in length, excluding their poly(A) tails. The highest amino acid sequence similarity (50 % identity) in the Pro-Pol core region encoded by RNA 1 was observed in the corresponding protein of strawberry latent ringspot virus. Based on pairwise comparisons and phylogenetic analysis, this virus, tentatively named "cohombrillo-associated virus" (CoAV), appears to be a member of a new species in the genus Stralarivirus (family Secoviridae), for which the name "Stralarivirus elaterii" is proposed. This new virus has different putative cleavage patterns from members of other species belonging to this genus.

Genetic Differentiation and Migration Fluxes of Viruses from Melon Crops and Crop Edge Weeds.

Weeds surrounding crops may act as alternative hosts, playing important epidemiological roles as virus reservoirs and impacting virus evolution. We used high-throughput sequencing to identify viruses in Spanish melon crops and plants belonging to three pluriannual weed species, Ecballium elaterium, Malva sylvestris, and Solanum nigrum, sampled at the edges of the crops. Melon and E. elaterium, both belonging to the family Cucurbitaceae, shared three virus species, whereas there was no virus species overlap between melon and the other two weeds. The diversity of cucurbit aphid-borne yellows virus (CABYV) and tomato leaf curl New Delhi virus (ToLCNDV), both in melon and E. elaterium, was further studied by amplicon sequencing. Phylogenetic and population genetics analyses showed that the CABYV population was structured by the host, identifying three sites in the CABYV RNA-dependent RNA polymerase under positive selection, perhaps reflecting host adaptation. The ToLCNDV population was much less diverse than the CABYV one, likely as a consequence of the relatively recent introduction of ToLCNDV in Spain. In spite of its low diversity, we identified geographical but no host differentiation for ToLCNDV. Potential virus migration fluxes between E. elaterium and melon plants were also analyzed. For CABYV, no evidence of migration between the populations of the two hosts was found, whereas important fluxes were identified between geographically distant subpopulations for each host. For ToLCNDV, in contrast, evidence of migration from melon to E. elaterium was found, but not the other way around. IMPORTANCE It has been reported that about half of the emerging diseases affecting plants are caused by viruses. Alternative hosts often play critical roles in virus emergence as virus reservoirs, bridging host species that are otherwise unconnected and/or favoring virus diversification. In spite of this, the viromes of potential alternative hosts remain largely unexplored. In the case of crops, pluriannual weeds at the crop edges may play these roles. Here, we took advantage of the power of high-throughput sequencing to characterize the viromes of three weed species frequently found at the edges of melon crops. We identified three viruses shared by melon and the cucurbit weed, with two of them being epidemiologically relevant for melon crops. Further genetic analyses showed that these two viruses had contrasting patterns of diversification and migration, providing an interesting example on the role that weeds may play in the ecology and evolution of viruses affecting crops.

Editing melon eIF4E associates with virus resistance and male sterility.

The cap-binding protein eIF4E, through its interaction with eIF4G, constitutes the core of the eIF4F complex, which plays a key role in the circularization of mRNAs and their subsequent cap-dependent translation. In addition to its fundamental role in mRNA translation initiation, other functions have been described or suggested for eIF4E, including acting as a proviral factor and participating in sexual development. We used CRISPR/Cas9 genome editing to generate melon eif4e knockout mutant lines. Editing worked efficiently in melon, as we obtained transformed plants with a single-nucleotide deletion in homozygosis in the first eIF4E exon already in a T0 generation. Edited and non-transgenic plants of a segregating F2 generation were inoculated with Moroccan watermelon mosaic virus (MWMV); homozygous mutant plants showed virus resistance, while heterozygous and non-mutant plants were infected, in agreement with our previous results with plants silenced in eIF4E. Interestingly, all homozygous edited plants of the T0 and F2 generations showed a male sterility phenotype, while crossing with wild-type plants restored fertility, displaying a perfect correlation between the segregation of the male sterility phenotype and the segregation of the eif4e mutation. Morphological comparative analysis of melon male flowers along consecutive developmental stages showed postmeiotic abnormal development for both microsporocytes and tapetum, with clear differences in the timing of tapetum degradation in the mutant versus wild-type. An RNA-Seq analysis identified critical genes in pollen development that were down-regulated in flowers of eif4e/eif4e plants, and suggested that eIF4E-specific mRNA translation initiation is a limiting factor for male gametes formation in melon.

Transcriptome analyses unveiled differential regulation of AGO and DCL genes by pepino mosaic virus strains.

Pepino mosaic virus (PepMV) is a single-stranded (ss), positive-sense (+) RNA potexvirus that affects tomato crops worldwide. We have described an in planta antagonistic interaction between PepMV isolates of two strains in which the EU isolate represses the accumulation of the CH2 isolate during mixed infections. Reports describing transcriptomic responses to mixed infections are scant. We carried out transcriptomic analyses of tomato plants singly and mixed-infected with two PepMV isolates of both strains. Comparison of the transcriptomes of singly infected plants showed that deeper transcriptomic alterations occurred at early infection times, and also that each of the viral strains modulated the host transcriptome differentially. Mixed infections caused transcriptomic alterations similar to those for the sum of single infections at early infection times, but clearly differing at later times postinfection. We next tested the hypothesis that PepMV-EU, in either single or mixed infections, deregulates host gene expression differentially so that virus accumulation of both strains gets repressed. That seemed to be the case for the genes AGO1a, DCL2d, AGO2a, and DCL2b, which are involved in the antiviral silencing pathway and were upregulated by PepMV-EU but not by PepMV-CH2 at early times postinfection. The pattern of AGO2a expression was validated by reverse transcription-quantitative PCR in tomato and Nicotiana benthamiana plants. Using an N. benthamiana ago2 mutant line, we showed that AGO2 indeed plays an important role in the antiviral defence against PepMV, but it is not the primary determinant of the outcome of the antagonistic interaction between the two PepMV strains.

Complete genome sequence of malva-associated soymovirus 1: a novel virus infecting common mallow.

In this work, a novel viral genomic sequence with a gene organization typical of members of the genus Soymovirus was identified using high-throughput sequencing data from common mallow. This species is a vigorous wild weed native to the Mediterranean region, commonly found in borders and edges of cultivated fields, making it a suitable reservoir for plant pests and pathogens. Indeed, plant viruses belonging to different genera have been previously found infecting common malva. This new viral genome consists of a single molecule of circular double-stranded DNA of 8391 base pairs and contains eight open reading frames encoding polymerase, movement, coat, translational transactivator protein typical of caulimoviruses, and four hypothetical proteins. Phylogenetic and pairwise distance analyses showed its close relationship with soybean chlorotic mottle virus. Interestingly, a small intergenic region was detected between ORFs Ib and II. Based on the demarcation criteria of the genus Soymovirus, the new virus, provisionally named malva-associated soymovirus 1, could be a member of a new species Soymovirus masolus. To our knowledge, this is the first report of a soymovirus infecting common mallow.

Use of High-Throughput Sequencing and Two RNA Input Methods to Identify Viruses Infecting Tomato Crops.

We used high-throughput sequencing to identify viruses on tomato samples showing virus-like symptoms. Samples were collected from crops in the Iberian Peninsula. Either total RNA or double-stranded RNA (dsRNA) were used as starting material to build the cDNA libraries. In total, seven virus species were identified, with pepino mosaic virus being the most abundant one. The dsRNA input provided better coverage and read depth but missed one virus species compared with the total RNA input. By performing in silico analyses, we determined a minimum sequencing depth per sample of 0.2 and 1.5 million reads for dsRNA and rRNA-depleted total RNA inputs, respectively, to detect even the less abundant viruses. Primers and TaqMan probes targeting conserved regions in the viral genomes were designed and/or used for virus detection; all viruses were detected by qRT-PCR/RT-PCR in individual samples, with all except one sample showing mixed infections. Three virus species (Olive latent virus 1, Lettuce ring necrosis virus and Tomato fruit blotch virus) are herein reported for the first time in tomato crops in Spain.

Virus-Infected Melon Plants Emit Volatiles that Induce Gene Deregulation in Neighboring Healthy Plants.

It is well described that viral infections stimulate the emission of plant volatiles able to recruit viral vectors thereby promoting virus spread. In contrast, much less is known on the effects that emitted volatiles may have on the metabolism of healthy neighboring plants, which are potential targets for new infections through vector transmission. Watermelon mosaic virus (WMV) (genus Potyvirus, family Potyviridae) is an aphid-transmitted virus endemic in cucurbit crops worldwide. We have compared gene expression profiles of WMV-infected melon plants with those of healthy or healthy-but-cohabited-with-infected plants. Pathogenesis-related (PR) and small heat shock protein encoding genes were deregulated in cohabited plants, and PR deregulation depended on the distance to the infected plant. The signaling was short distance in the experimental conditions used, and cohabiting had a moderate effect on the plant susceptibility to WMV. Static headspace experiments showed that benzaldehyde and γ-butyrolactone were significantly over-emitted by WMV-infected plants. Altogether, our data suggest that perception of a volatile signal encoded by WMV-infected tissues triggers a response to prepare healthy tissues or/and healthy neighboring plants for the incoming infections.

Turnip mosaic virus in oilseed rape activates networks of sRNA-mediated interactions between viral and host genomes.

Virus-induced plant diseases in cultivated plants cause important damages in yield. Although the mechanisms of virus infection are intensely studied at the cell biology level, only little is known about the molecular dialog between the invading virus and the host genome. Here we describe a combinatorial genome-wide approach to identify networks of sRNAs-guided post-transcriptional regulation within local Turnip mosaic virus (TuMV) infection sites in Brassica napus leaves. We show that the induction of host-encoded, virus-activated small interfering RNAs (vasiRNAs) observed in virus-infected tissues is accompanied by site-specific cleavage events on both viral and host RNAs that recalls the activity of small RNA-induced silencing complexes (RISC). Cleavage events also involve virus-derived siRNA (vsiRNA)-directed cleavage of target host transcripts as well as cleavage of viral RNA by both host vasiRNAs and vsiRNAs. Furthermore, certain coding genes act as virus-activated regulatory hubs to produce vasiRNAs for the targeting of other host genes. The observations draw an advanced model of plant-virus interactions and provide insights into the complex regulatory networking at the plant-virus interface within cells undergoing early stages of infection.

ICTV Virus Taxonomy Profile: Botourmiaviridae.

The family Botourmiaviridae includes viruses infecting plants and filamentous fungi containing a positive-sense, ssRNA genome that can be mono- or multi-segmented. Genera in the family include: Ourmiavirus (plant viruses), and Botoulivirus, Magoulivirus and Scleroulivirus (fungal viruses). This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the taxonomy of the family Botourmiaviridae, which is available at ictv.global/report/botourmiaviridae.

Computational Workflow for Small RNA Profiling in Virus-Infected Plants.

In this chapter we describe a series of computational pipelines for the in silico analysis of small RNAs (sRNA) produced in response to viral infections in plants. Our workflow is primarily focused on the analysis of sRNA populations derived from known or previously undescribed viruses infecting host plants. Furthermore, we provide an additional pipeline to examine host-specific endogenous sRNAs activated or specifically expressed during viral infections in plants. We present some key points for a successful and cost-efficient processing of next generation sequencing sRNA libraries, from purification of high quality RNA to guidance for library preparation and sequencing strategies. We report a series of free available tools and programs as well as in-house Perl scripts to perform customized sRNA-seq data mining. Previous bioinformatic background is not required, but experience with basic Unix commands is desirable.

The immune repressor BIR1 contributes to antiviral defense and undergoes transcriptional and post-transcriptional regulation during viral infections.

BIR1 is a receptor-like kinase that functions as a negative regulator of basal immunity and cell death in Arabidopsis. Using Arabidopsis thaliana and Tobacco rattle virus (TRV), we investigate the antiviral role of BIR1, the molecular mechanisms of BIR1 gene expression regulation during viral infections, and the effects of BIR1 overexpression on plant immunity and development. We found that SA acts as a signal molecule for BIR1 activation during infection. Inactivating mutations of BIR1 in the bir1-1 mutant cause strong antiviral resistance independently of constitutive cell death or SA defense priming. BIR1 overexpression leads to severe developmental defects, cell death and premature death, which correlate with the constitutive activation of plant immune responses. Our findings suggest that BIR1 acts as a negative regulator of antiviral defense in plants, and indicate that RNA silencing contributes, alone or in conjunction with other regulatory mechanisms, to define a threshold expression for proper BIR1 function beyond which an autoimmune response may occur. This work provides novel mechanistic insights into the regulation of BIR1 homeostasis that may be common for other plant immune components.

Characterization of Begomoviruses Sampled during Severe Epidemics in Tomato Cultivars Carrying the Ty-1 Gene.

Tomato yellow leaf curl virus (TYLCV, genus Begomovirus, family Geminiviridae) is a major species that causes a tomato disease for which resistant tomato hybrids (mainly carriers of the Ty-1/Ty-3 gene) are being used widely. We have characterized begomoviruses severely affecting resistant tomato crops in Southeast Spain. Circular DNA was prepared from samples by rolling circle amplification, and sequenced by massive sequencing (2015) or cloning and Sanger sequencing (2016). Thus, 23 complete sequences were determined, all belonging to the TYLCV Israel strain (TYLCV-IL). Massive sequencing also revealed the absence of other geminiviral and beta-satellite sequences. A phylogenetic analysis showed that the Spanish isolates belonged to two groups, one related to early TYLCV-IL isolates in the area (Group 1), and another (Group 2) closely related to El Jadida (Morocco) isolates, suggesting a recent introduction. The most parsimonious evolutionary scenario suggested that the TYLCV isolates of Group 2 are back recombinant isolates derived from TYLCV-IS76, a recombinant virus currently predominating in Moroccan epidemics. Thus, an infectious Group 2 clone (TYLCV-Mu15) was constructed and used in in planta competition assays against TYLCV-IS76. TYLCV-Mu15 predominated in single infections, whereas TYLCV-IS76 did so in mixed infections, providing credibility to a scenario of co-occurrence of both types of isolates.

Small RNA-Seq to Characterize Viruses Responsible of Lettuce Big Vein Disease in Spain.

The emerging lettuce big-vein disease (LBVD) is causing losses in lettuce production ranging from 30 to 70% worldwide. Several studies have associated this disease with Mirafiori lettuce big-vein virus (MiLBVV) alone or in mixed infection with lettuce big-vein associated virus (LBVaV). We used Illumina small RNA sequencing (sRNA-seq) to identify viruses present in symptomatic lettuce plants from commercial fields in Southern Spain. Data analysis using the VirusDetect tool showed the consistent presence of MiLBVV and LBVaV in diseased plants. Populations of MiLBVV and LBVaV viral small RNAs (sRNAs) were characterized, showing features essentially similar to those of other viruses, with the peculiarity of an uneven asymmetric distribution of MiLBVV virus-derived small RNAs (vsRNAs) for the different polarities of genomic RNA4 vs. RNAs1 to 3. Sanger sequencing of coat protein genes was used to study MiLBVV and LBVaV phylogenetic relationships and population genetics. The Spanish MiLBVV population was composed of isolates from three well-differentiated lineages and reflected almost all of the diversity reported for the MiLBVV species, whereas the LBVaV population showed very little genetic differentiation at the regional scale but lineage differentiation at a global geographical scale. Universal primers were used to detect and quantify the accumulation of MiLBVV and LBVaV in field samples; both symptomatic and asymptomatic plants from affected fields carried equal viral loads, with LBVaV accumulating at higher levels than MiLBVV.

Potato Virus Y HCPro Suppression of Antiviral Silencing in Nicotiana benthamiana Plants Correlates with Its Ability To Bind In Vivo to 21- and 22-Nucleotide Small RNAs of Viral Sequence.

We have investigated short and small RNAs (sRNAs) that were bound to a biologically active hexahistidine-tagged Potato virus Y (PVY) HCPro suppressor of silencing, expressed from a heterologous virus vector in Nicotiana benthamiana plants, and purified under nondenaturing conditions. We found that RNAs in purified preparations were differentially enriched in 21-nucleotide (nt) and, to a much lesser extent, 22-nt sRNAs of viral sequences (viral sRNAs [vsRNAs]) compared to those found in a control plant protein background bound to nickel resin in the absence of HCPro or in a purified HCPro alanine substitution mutant (HCPro mutB) control that lacked suppressor-of-silencing activity. In both controls, sRNAs were composed almost entirely of molecules of plant sequence, indicating that the resin-bound protein background had no affinity for vsRNAs and also that HCPro mutB failed to bind to vsRNAs. Therefore, PVY HCPro suppressor activity correlated with its ability to bind to 21- and 22-nt vsRNAs. HCPro constituted at least 54% of the total protein content in purified preparations, and we were able to calculate its contribution to the 21- and the 22-nt pools of sRNAs present in the purified samples and its binding strength relative to the background. We also found that in the 21-nt vsRNAs of the HCPro preparation, 5'-terminal adenines were overrepresented relative to the controls, but this was not observed in vsRNAs of other sizes or of plant sequences.IMPORTANCE It was previously shown that HCPro can bind to long RNAs and small RNAs (sRNAs) in vitro and, in the case of Turnip mosaic virus HCPro, also in vivo in arabidopsis AGO2-deficient plants. Our data show that PVY HCPro binds in vivo to sRNAs during infection in wild-type Nicotiana benthamiana plants when expressed from a heterologous virus vector. Using a suppression-of-silencing-deficient HCPro mutant that can accumulate in this host when expressed from a virus vector, we also show that sRNA binding correlates with silencing suppression activity. We demonstrate that HCPro binds at least to sRNAs with viral sequences of 21 nucleotides (nt) and, to a much lesser extent, of 22 nt, which were are also differentially enriched in 5'-end adenines relative to the purified controls. Together, our results support the physical binding of HCPro to vsRNAs of 21 and 22 nt as a means to interfere with antiviral silencing.

Deep sequencing of mycovirus-derived small RNAs from Botrytis species.

RNA silencing is an ancient regulatory mechanism operating in all eukaryotic cells. In fungi, it was first discovered in Neurospora crassa, although its potential as a defence mechanism against mycoviruses was first reported in Cryphonectria parasitica and, later, in several fungal species. There is little evidence of the antiviral potential of RNA silencing in the phytopathogenic species of the fungal genus Botrytis. Moreover, little is known about the RNA silencing components in these fungi, although the analysis of public genome databases identified two Dicer-like genes in B. cinerea, as in most of the ascomycetes sequenced to date. In this work, we used deep sequencing to study the virus-derived small RNA (vsiRNA) populations from different mycoviruses infecting field isolates of Botrytis spp. The mycoviruses under study belong to different genera and species, and have different types of genome [double-stranded RNA (dsRNA), (+)single-stranded RNA (ssRNA) and (-)ssRNA]. In general, vsiRNAs derived from mycoviruses are mostly of 21, 20 and 22 nucleotides in length, possess sense or antisense orientation, either in a similar ratio or with a predominance of sense polarity depending on the virus species, have predominantly U at their 5' end, and are unevenly distributed along the viral genome, showing conspicuous hotspots of vsiRNA accumulation. These characteristics reveal striking similarities with vsiRNAs produced by plant viruses, suggesting similar pathways of viral targeting in plants and fungi. We have shown that the fungal RNA silencing machinery acts against the mycoviruses used in this work in a similar manner independent of their viral or fungal origin.

Characterization of Botrytis cinerea negative-stranded RNA virus 1, a new mycovirus related to plant viruses, and a reconstruction of host pattern evolution in negative-sense ssRNA viruses.

The molecular characterization of a novel negative single-stranded RNA virus infecting the plant pathogenic fungus Botrytis cinerea is reported here. Comparison of the sequence of Botrytis cinerea negative-stranded RNA virus 1 (BcNSRV-1) showed a strong identity with RNA dependent RNA polymerases (RdRps) of plant pathogenic emaraviruses and tospoviruses. We have also found all the molecular signatures present in the RdRp of the genus Emaravirus and in other genera of family Bunyaviridae: the conserved TPD triplet and RY dinucleotide, the three basic residues in premotif A and the conserved motifs A, B, C, D, and E. Our results showed that BcNSRV-1 is phylogenetically close to members of the genus Emaravirus and of the family Bunyaviridae, and an ancestral state reconstruction using the conserved RdRp motifs of type members of each family of (-)ssRNA viruses indicated that BcNSRV-1 could possibly derive from an invertebrate and vertebrate-infecting virus.