Identification of DNA Viruses in Ancient DNA from Herbarium Samples

Herbaria encompass millions of plant specimens, mostly collected in the nineteenth and twentieth centuries that can represent a key resource for investigating the history and evolution of phytopathogens. In the last years, the application of high-throughput sequencing technologies for the analysis of ancient nucleic acids has revolutionized the study of ancient pathogens including viruses, allowing the reconstruction of historical genomic viral sequences, improving phylogenetic based molecular dating, and providing essential insight into plant virus ecology. In this chapter, we describe a protocol to reconstruct ancient plant and soil viral sequences starting from highly fragmented ancient DNA extracted from herbarium plants and their associated rhizospheric soil. Following Illumina high-throughput sequencing, sequence data are de novo assembled, and DNA viral sequences are selected, according to their similarity with known viruses.

Modulation of class III peroxidase pathways and phenylpropanoids in Arundo donax under salt and phosphorus stress

Arundo donax L. is an invasive species that has been recently employed for biomass production due to its well-known ability to colonize harsh environment. Based on previous observations, the present study investigated the potential role of phenylpropanoids and class III peroxidases to confer adaptation through biochemical and transcriptomic analysis in A. donax after Na⁺ and P excess supply, both in single stress and in combination, and after growth at low P level. The levels of hydrogen peroxide, flavonoids (i.e., quercetin, apigenin and kaempferol derivatives) and the activity of class III peroxidases, as well as the expression of several genes encoding for their enzymes involved in their biosynthesis, increased when Na⁺ was supplied in combination with P. These results suggest that those biomolecules are involved in the response of A. donax, to the presence of +Na and P in the soil. Moreover, even though at the sampling time no significant accumulation of lignin has been determined, the trend of accumulation of such metabolite and most of all the increase of several transcripts involved in its synthesis was found. This work for the first time indicates the need for further investigation devoted to elucidating whether the strengthening of cell walls via lignin synthesis is one of the mechanisms used by A. donax to adapt to harsh environments.

Women in the European Virus Bioinformatics Center

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

No Evidence for Seed Transmission of Tomato Yellow Leaf Curl Sardinia Virus in Tomato

Seed transmission is an important factor in the epidemiology of plant pathogens. Geminiviruses are serious pests spread in tropical and subtropical regions. They are transmitted by hemipteran insects, but a few cases of transmission through seeds were recently reported. Here, we investigated the tomato seed transmissibility of the begomovirus tomato yellow leaf curl Sardinia virus (TYLCSV), one of the agents inducing the tomato yellow leaf curl disease, heavily affecting tomato crops in the Mediterranean area. None of the 180 seedlings originating from TYLCSV-infected plants showed any phenotypic alteration typical of virus infection. Moreover, whole viral genomic molecules could not be detected in their cotyledons and true leaves, neither by membrane hybridization nor by rolling-circle amplification followed by PCR, indicating that TYLCSV is not a seed-transmissible pathogen for tomato. Examining the localization of TYLCSV DNA in progenitor plants, we detected the virus genome by PCR in all vegetative and reproductive tissues, but viral genomic and replicative forms were found only in leaves, flowers and fruit flesh, not in seeds and embryos. Closer investigations allowed us to discover for the first time that these embryos were superficially contaminated by TYLCSV DNA but whole genomic molecules were not detectable. Therefore, the inability of TYLCSV genomic molecules to colonize tomato embryos during infection justifies the lack of seed transmissibility observed in this host.

A Primer on the Analysis of High-Throughput Sequencing Data for Detection of Plant Viruses

High-throughput sequencing (HTS) technologies have become indispensable tools assisting plant virus diagnostics and research thanks to their ability to detect any plant virus in a sample without prior knowledge. As HTS technologies are heavily relying on bioinformatics analysis of the huge amount of generated sequences, it is of utmost importance that researchers can rely on efficient and reliable bioinformatic tools and can understand the principles, advantages, and disadvantages of the tools used. Here, we present a critical overview of the steps involved in HTS as employed for plant virus detection and virome characterization. We start from sample preparation and nucleic acid extraction as appropriate to the chosen HTS strategy, which is followed by basic data analysis requirements, an extensive overview of the in-depth data processing options, and taxonomic classification of viral sequences detected. By presenting the bioinformatic tools and a detailed overview of the consecutive steps that can be used to implement a well-structured HTS data analysis in an easy and accessible way, this paper is targeted at both beginners and expert scientists engaging in HTS plant virome projects.

The Induction of an Effective dsRNA-Mediated Resistance Against Tomato Spotted Wilt Virus by Exogenous Application of Double-Stranded RNA Largely Depends on the Selection of the Viral RNA Target Region

Tomato spotted wilt virus (TSWV) is a devastating plant pathogen, causing huge crop losses worldwide. Unfortunately, due to its wide host range and emergence of resistance breaking strains, its management is challenging. Up to now, resistance to TSWV infection based on RNA interference (RNAi) has been achieved only in transgenic plants expressing parts of the viral genome or artificial microRNAs targeting it. Exogenous application of double-stranded RNAs (dsRNAs) for inducing virus resistance in plants, namely RNAi-based vaccination, represents an attractive and promising alternative, already shown to be effective against different positive-sense RNA viruses and viroids. In the present study, the protection efficacy of exogenous application of dsRNAs targeting the nucleocapsid (N) or the movement protein (NSm) coding genes of the negative-sense RNA virus TSWV was evaluated in Nicotiana benthamiana as model plant and in tomato as economically important crop. Most of the plants treated with N-targeting dsRNAs, but not with NSm-targeting dsRNAs, remained asymptomatic until 40 (N. benthamiana) and 63 (tomato) dpi, while the remaining ones showed a significant delay in systemic symptoms appearance. The different efficacy of N- and NSm-targeting dsRNAs in protecting plants is discussed in the light of their processing, mobility and biological role. These results indicate that the RNAi-based vaccination is effective also against negative-sense RNA viruses but emphasize that the choice of the target viral sequence in designing RNAi-based vaccines is crucial for its success.

Arbuscular Mycorrhizal Symbiosis Primes Tolerance to Cucumber Mosaic Virus in Tomato

Tomato plants can establish symbiotic interactions with arbuscular mycorrhizal fungi (AMF) able to promote plant nutrition and prime systemic plant defenses against pathogens attack; the mechanism involved is known as mycorrhiza-induced resistance (MIR). However, studies on the effect of AMF on viral infection, still limited and not conclusive, indicate that AMF colonization may have a detrimental effect on plant defenses against viruses, so that the term "mycorrhiza-induced susceptibility" (MIS) has been proposed for these cases. To expand the case studies to a not yet tested viral family, that is, Bromoviridae, we investigated the effect of the colonization by the AMF Funneliformis mosseae on cucumber mosaic virus (CMV) infection in tomato by phenotypic, physiological, biochemical, and transcriptional analyses. Our results showed that the establishment of a functional AM symbiosis is able to limit symptoms development. Physiological and transcriptomic data highlighted that AMF mitigates the drastic downregulation of photosynthesis-related genes and the reduction of photosynthetic CO2 assimilation rate caused by CMV infection. In parallel, an increase of salicylic acid level and a modulation of reactive oxygen species (ROS)-related genes, toward a limitation of ROS accumulation, was specifically observed in CMV-infected mycorrhizal plants. Overall, our data indicate that the AM symbiosis influences the development of CMV infection in tomato plants and exerts a priming effect able to enhance tolerance to viral infection.

Different Genetic Sources Contribute to the Small RNA Population in the Arbuscular Mycorrhizal Fungus Gigaspora margarita

RNA interference (RNAi) is a key regulatory pathway of gene expression in almost all eukaryotes. This mechanism relies on short non-coding RNA molecules (sRNAs) to recognize in a sequence-specific manner DNA or RNA targets leading to transcriptional or post-transcriptional gene silencing. To date, the fundamental role of sRNAs in the regulation of development, stress responses, defense against viruses and mobile elements, and cross-kingdom interactions has been extensively studied in a number of biological systems. However, the knowledge of the “RNAi world” in arbuscular mycorrhizal fungi (AMF) is still limited. AMF are obligate mutualistic endosymbionts of plants, able to provide several benefits to their partners, from improved mineral nutrition to stress tolerance. Here we described the RNAi-related genes of the AMF Gigaspora margarita and characterized, through sRNA sequencing, its complex small RNAome, considering the possible genetic sources and targets of the sRNAs. G. margarita indeed is a mosaic of different genomes since it hosts endobacteria, RNA viruses, and non-integrated DNA fragments corresponding to mitovirus sequences. Our findings show that G. margarita is equipped with a complete set of RNAi-related genes characterized by the expansion of the Argonaute-like (AGO-like) gene family that seems a common trait of AMF. With regards to sRNAs, we detected populations of sRNA reads mapping to nuclear, mitochondrial, and viral genomes that share similar features (25-nt long and 5′-end uracil read enrichments), and that clearly differ from sRNAs of endobacterial origin. Furthermore, the annotation of nuclear loci producing sRNAs suggests the occurrence of different sRNA-generating processes. In silico analyses indicate that the most abundant G. margarita sRNAs, including those of viral origin, could target transcripts in the host plant, through a hypothetical cross-kingdom RNAi.

Impact of high or low levels of phosphorus and high sodium in soils on productivity and stress tolerance of Arundo donax plants

The potential of Arundo donax to grow in degraded soils, characterized by excess of salinity (Na+), and phosphorus deficiency (-P) or excess (+P) also coupled with salinity (+NaP), was investigated by combining in vivo plant phenotyping, quantification of metabolites and ultrastructural imaging of leaves with a transcriptome-wide screening. Photosynthesis and growth were impaired by + Na, -P and + NaP. While + Na caused stomatal closure, enhanced biosynthesis of carotenoids, sucrose and isoprene and impaired anatomy of cell walls, +P negatively affected starch production and isoprene emission, and damaged chloroplasts. Finally, +NaP largely inhibited photosynthesis due to stomatal limitations, increased sugar content, induced/repressed a number of genes 10 time higher with respect to + P and + Na, and caused appearance of numerous and large plastoglobules and starch granules in chloroplasts. Our results show that A. donax is sensitive to unbalances of soil ion content, despite activation of defensive mechanisms that enhance plant resilience, growth and biomass production of A. donax under these conditions.

Nondestructive Raman Spectroscopy as a Tool for Early Detection and Discrimination of the Infection of Tomato Plants by Two Economically Important Viruses

Global population forecasts dictate a rapid adoption of multifaceted approaches to fulfill increasing food requirements, ameliorate food dietary value and security using sustainable and economically feasible agricultural processes. Plant pathogens induce up to 25% losses in vegetable crops and their early detection would contribute to limit their spread and economic impact. As an alternative to time-consuming, destructive, and expensive diagnostic procedures, such as immunological assays and nucleic acid-based techniques, Raman spectroscopy (RS) is a nondestructive rapid technique that generates a chemical fingerprinting of a sample, at low operating costs. Here, we assessed the suitability of RS combined to chemometric analysis to monitor the infection of an important vegetable crop plant, tomato, by two dangerous and peculiarly different viral pathogens, Tomato yellow leaf curl Sardinia virus (TYLCSV) and Tomato spotted wilt virus (TSWV). Experimentally inoculated plants were monitored over 28 days for symptom occurrence and subjected to RS analysis, alongside with measuring the virus amount by quantitative real-time PCR. RS allowed to discriminate mock inoculated (healthy) from virus-infected specimens, reaching an accuracy of >70% after only 14 days after inoculation for TYLCSV and >85% only after 8 days for TSWV, demonstrating its suitability for early detection of virus infection. Importantly, RS also highlighted spectral differences induced by the two viruses, providing specific information on the infecting agent.

Arbuscular Mycorrhizal Symbiosis: Plant Friend or Foe in the Fight Against Viruses?

Plant roots establish interactions with several beneficial soil microorganisms including arbuscular mycorrhizal fungi (AMF). In addition to promoting plant nutrition and growth, AMF colonization can prime systemic plant defense and enhance tolerance to a wide range of environmental stresses and below-ground pathogens. A protective effect of the AMF against above-ground pathogens has also been described in different plant species, but it seems to largely rely on the type of attacker. Viruses are obligate biotrophic pathogens able to infect a large number of plant species, causing massive losses in crop yield worldwide. Despite their economic importance, information on the effect of the AM symbiosis on viral infection is limited and not conclusive. However, several experimental evidences, obtained under controlled conditions, show that AMF colonization may enhance viral infection, affecting susceptibility, symptomatology and viral replication, possibly related to the improved nutritional status and to the delayed induction of pathogenesis-related proteins in the mycorrhizal plants. In this review, we give an overview of the impact of the AMF colonization on plant infection by pathogenic viruses and summarize the current knowledge of the underlying mechanisms. For the cases where AMF colonization increases the susceptibility of plants to viruses, the term "mycorrhiza-induced susceptibility" (MIS) is proposed.

In silico analysis of fungal small RNA accumulation reveals putative plant mRNA targets in the symbiosis between an arbuscular mycorrhizal fungus and its host plant

BACKGROUND

Small RNAs (sRNAs) are short non-coding RNA molecules (20-30 nt) that regulate gene expression at transcriptional or post-transcriptional levels in many eukaryotic organisms, through a mechanism known as RNA interference (RNAi). Recent studies have highlighted that they are also involved in cross-kingdom communication: sRNAs can move across the contact surfaces from "donor" to "receiver" organisms and, once in the host cells of the receiver, they can target specific mRNAs, leading to a modulation of host metabolic pathways and defense responses. Very little is known about RNAi mechanism and sRNAs occurrence in Arbuscular Mycorrhizal Fungi (AMF), an important component of the plant root microbiota that provide several benefits to host plants, such as improved mineral uptake and tolerance to biotic and abiotic stress.

RESULTS

Taking advantage of the available genomic resources for the AMF Rhizophagus irregularis we described its putative RNAi machinery, which is characterized by a single Dicer-like (DCL) gene and an unusual expansion of Argonaute-like (AGO-like) and RNA-dependent RNA polymerase (RdRp) gene families. In silico investigations of previously published transcriptomic data and experimental assays carried out in this work provided evidence of gene expression for most of the identified sequences. Focusing on the symbiosis between R. irregularis and the model plant Medicago truncatula, we characterized the fungal sRNA population, highlighting the occurrence of an active sRNA-generating pathway and the presence of microRNA-like sequences. In silico analyses, supported by host plant degradome data, revealed that several fungal sRNAs have the potential to target M. truncatula transcripts, including some specific mRNA already shown to be modulated in roots upon AMF colonization.

CONCLUSIONS

The identification of RNAi-related genes, together with the characterization of the sRNAs population, suggest that R. irregularis is equipped with a functional sRNA-generating pathway. Moreover, the in silico analysis predicted 237 plant transcripts as putative targets of specific fungal sRNAs suggesting that cross-kingdom post-transcriptional gene silencing may occur during AMF colonization.

A Short Indel-Lacking-Resistance Gene Triggers Silencing of the Photosynthetic Machinery Components Through TYLCSV-Associated Endogenous siRNAs in Tomato

Plant viruses modify gene expression in infected tissues by altering the micro (mi)RNA-mediated regulation of genes. Among conserved miRNA targets there are transcripts coding for transcription factors, RNA silencing core, and disease-resistance proteins. Paralogs in these gene families are widely present in plant genomes and are known to respond differently to miRNA-mediated regulation during plant virus infections. Using genome-wide approaches applied to Solanum lycopersicum infected by a nuclear-replicating virus, we highlighted miRNA-mediated cleavage events that could not be revealed in virus-free systems. Among them we confirmed miR6024 targeting and cleavage of RX-coiled-coil (RX-CC), nucleotide binding site (NBS), leucine-rich (LRR) mRNA. Cleavage of paralogs was associated with short indels close to the target sites, indicating a general functional significance of indels in fine-tuning gene expression in plant-virus interaction. miR6024-mediated cleavage, uniquely in virus-infected tissues, triggers the production of several 21-22 nt secondary siRNAs. These secondary siRNAs, rather than being involved in the cascade regulation of other NBS-LRR paralogs, explained cleavages of several mRNAs annotated as defence-related proteins and components of the photosynthetic machinery. Outputs of these data explain part of the phenotype plasticity in plants, including the appearance of yellowing symptoms in the viral pathosystem.

The interaction of Lolium latent virus major coat protein with ankyrin repeat protein NbANKr redirects it to chloroplasts and modulates virus infection

The Lolium latent virus (LoLV) major coat protein sequence contains a typical chloroplast transit peptide (cTP) domain. In infected Nicotiana benthamiana leaf tissue, LoLV coat proteins can be detected at the chloroplast. In transient expression, several N-terminal deletions of the CP sequence, increasing in length, result in disruption of the domain functionality, markedly affecting intracellular localization. A yeast two-hybrid-based study using LoLV CP as bait identified several potentially interacting Arabidopsis host proteins, most of them with chloroplast-linked pathways. One of them, an ankyrin repeat protein, was studied in detail. The N. benthamiana homologue (NbANKr) targets chloroplasts, is able to co-localize with LoLV CP at chloroplast membranes in transient expression and shows a robust interaction with LoLV CP in vivo by BiFC, which has been confirmed by yeast two-hybrid data. Silencing NbANKr genes in N. benthamiana plants, prior to challenging with LoLV by mechanical inoculation, affects LoLV infection, significantly reducing the level of viral RNA in young leaves, compared to levels in control plants, and suggesting an inhibition of virus movement. Silencing of NbANKr has no obvious effect on plant phenotype, but is able to interfere with LoLV infection, opening the way for a new strategy for virus infection control.

Deep Sequencing Data and Infectivity Assays Indicate that Chickpea Chlorotic Dwarf Virus is the Etiological Agent of the "Hard Fruit Syndrome" of Watermelon

Chickpea chlorotic dwarf virus (CpCDV), a polyphagous mastrevirus, family Geminiviridae, has been recently linked to the onset of the "hard fruit syndrome" of watermelon, first described in Tunisia, that makes fruits unmarketable due to the presence of white hard portions in the flesh, chlorotic mottling on the rind, and an unpleasant taste. To investigate the etiological agent of this disease, total RNA extracted from symptomatic watermelon fruits was subjected to small RNA sequencing through next generation sequencing (NGS) techniques. Data obtained showed the presence of CpCDV and two other viral species. However, following validation through polymerase chain reaction (PCR), CpCDV was the only viral species consistently detected in all samples. Watermelon seedlings were then challenged by an agroinfectious CpCDV clone; several plants proved to be CpCDV-infected, and were able to produce fruits. CpCDV infected and replicated in watermelon fruits and leaves, leading to abnormality in fruits and in seed production, similar to those described in field. These results indicate that CpCDV is the etiological agent of the "hard fruit syndrome" of watermelon.

Novel functional microRNAs from virus-free and infected Vitis vinifera plants under water stress

MicroRNAs (miRNAs) are small non-coding RNAs that regulate the post-transcriptional control of several pathway intermediates, thus playing pivotal roles in plant growth, development and response to biotic and abiotic stresses. In recent years, the grapevine genome release, small(s)-RNAseq and degradome-RNAseq together has allowed the discovery and characterisation of many miRNA species, thus rendering the discovery of additional miRNAs difficult and uncertain. Taking advantage of the miRNA responsiveness to stresses and the availability of virus-free Vitis vinifera plants and those infected only by a latent virus, we have analysed grapevines subjected to drought in greenhouse conditions. The sRNA-seq and other sequence-specific molecular analyses have allowed us to characterise conserved miRNA expression profiles in association with specific eco-physiological parameters. In addition, we here report 12 novel grapevine-specific miRNA candidates and describe their expression profile. We show that latent viral infection can influence the miRNA profiles of V. vinifera in response to drought. Moreover, study of eco-physiological parameters showed that photosynthetic rate, stomatal conductance and hydraulic resistance to water transport were significantly influenced by drought and viral infection. Although no unequivocal cause–effect explanation could be attributed to each miRNA target, their contribution to the drought response is discussed.

Comparison of small RNA profiles in Nicotiana benthamiana and Solanum lycopersicum infected by polygonum ringspot tospovirus reveals host-specific responses to viral infection

Viral small RNAs (vsRNAs) are one of the key elements involved in RNA silencing-based defense against viruses in plants. We analyzed the vsRNA profiles in Nicotiana benthamiana and Solanum lycopersicum infected by polygonum ringspot virus (PolRSV) (Tospovirus, Bunyaviridae). VsRNAs were abundant in both hosts, but a different size profile was observed, with an abundance peak at 21 in N. benthamiana and at 22 nt in tomato. VsRNAs mapping to the PolRSV L genomic segment were under-represented in both hosts, while S and M segments were differentially and highly targeted in N. benthamiana and tomato, respectively. Differences in preferential targeting of single ORFs were observed, with over-representation of NSs ORF-derived reads in N. benthamiana. Intergenic regions (IGRs)-mapping vsRNAs were under-represented, while enrichment of vsRNAs reads mapping to the NSs positive sense strand was observed in both hosts. Comparison with a previous study on tomato spotted wilt virus (TSWV) under the same experimental conditions, showed that the relative accumulation of PolRSV-specific and endogenous sRNAs was similar to the one observed for silencing suppressor-deficient TSWV strains, suggesting possible different properties of PolRSV NSs silencing suppressor compared to that of TSWV.

Small RNA profiles of wild-type and silencing suppressor-deficient tomato spotted wilt virus infected Nicotiana benthamiana

Tospoviruses are plant-infecting viruses belonging to the family Bunyaviridae. We used a collection of wild-type, phylogenetically distinct tomato spotted wilt virus isolates and related silencing-suppressor defective mutants to study the effects on the small RNA (sRNA) accumulation during infection of Nicotiana benthamiana. Our data showed that absence of a functional silencing suppressor determined a marked increase of the total amount of viral sRNAs (vsRNAs), and specifically of the 21 nt class. We observed a common under-representation of vsRNAs mapping to the intergenic region of S and M genomic segments, and preferential mapping of the reads against the viral sense open reading frames, with the exception of the NSs gene. The NSs-mutant strains showed enrichment of NSm-derived vsRNA compared to the expected amount based on gene size. Analysis of 5' terminal nucleotide preference evidenced a significant enrichment in U for the 21 nt- and in A for 24 nt-long endogenous sRNAs in all the samples. Hotspot analysis revealed a common abundant accumulation of reads at the 5' end of the L segment, mostly in the antiviral sense, for the NSs-defective isolates, suggesting that absence of the silencing suppressor can influence preferential targeting of the viral genome.

Real-time PCR protocols for the quantification of the begomovirus tomato yellow leaf curl Sardinia virus in tomato plants and in its insect vector

Tomato yellow leaf curl Sardinia virus (TYLCSV) (Geminiviridae) is an important pathogen, transmitted by the whitefly Bemisia tabaci, that severely affects the tomato production in the Mediterranean basin. Here, we describe real-time PCR protocols suitable for relative and absolute quantification of TYLCSV in tomato plants and in whitefly extracts. Using primers and probe specifically designed for TYLCSV, the protocols for relative quantification allow to compare the amount of TYLCSV present in different plant or whitefly samples, normalized to the amount of DNA present in each sample using endogenous tomato or Bemisia genes as internal references. The absolute quantification protocol allows to calculate the number of genomic units of TYLCSV over the genomic units of the plant host (tomato), with a sensitivity of as few as ten viral genome copies per sample. The described protocols are potentially suitable for several applications, such as plant breeding for resistance, analysis of virus replication, and virus-vector interaction studies.

Bioinformatics approaches for viral metagenomics in plants using short RNAs: model case of study and application to a Cicer arietinum population

Over the past years deep sequencing experiments have opened novel doors to reconstruct viral populations in a high-throughput and cost-effective manner. Currently a substantial number of studies have been performed which employ next generation sequencing techniques to either analyze known viruses by means of a reference-guided approach or to discover novel viruses using a de novo-based strategy. Taking advantage of the well-known Cymbidium ringspot virus we have carried out a comparison of different bioinformatics tools to reconstruct the viral genome based on 21–27 nt short (s)RNA sequencing with the aim to identify the most efficient pipeline. The same approach was applied to a population of plants constituting an ancient variety of Cicer arietinum with red seeds. Among the discovered viruses, we describe the presence of a Tobamovirus referring to the Tomato mottle mosaic virus (NC_022230), which was not yet observed on C. arietinum nor revealed in Europe and a viroid referring to Hop stunt viroid (NC_001351.1) never reported in chickpea. Notably, a reference sequence guided approach appeared the most efficient in such kind of investigation. Instead, the de novo assembly reached a non-appreciable coverage although the most prominent viral species could still be identified. Advantages and limitations of viral metagenomics analysis using sRNAs are discussed.

Drawing siRNAs of viral origin out from plant siRNAs libraries

Viruses are obligate intracellular entities that infect all forms of life. In plants, invading viral nucleic acids trigger RNA silencing machinery and it results in the accumulation of viral short interfering RNAs (v-siRNAs). The study of v-siRNAs population in biological samples has become a major part of many research projects aiming to identify viruses infecting them, including unknown viruses, even at extremely low titer. Currently, siRNA populations are investigated by high-throughput sequencing approaches, which generate very large data sets. The major difficulty in these studies is to properly analyze such huge amount of data. In this regard, easy-to-use bioinformatics tools to groom and decipher siRNA libraries and to draw out v-siRNAs are needed. Here we describe a workflow, which permit users with little experience in bioinformatics to draw out v-siRNAs from raw data sequences obtained by Illumina technology. Such pipeline has been released in the context of Galaxy, an open source Web-based platform for bioinformatics analyses.

The complete genome sequence of polygonum ringspot virus

The complete genome sequence of polygonum ringspot virus (PolRSV), genus Tospovirus, family Bunyaviridae, was determined. This is the first report of the complete genome sequence for a European tospovirus isolate. The large RNA of PolRSV was 8893 nucleotides (nt) in size and contained a single open reading frame of 8628 nucleotides in the viral-complementary sense, coding for a predicted RNA-dependent RNA polymerase of 330.9 kDa. Two untranslated regions of 230 and 32 nucleotides were present at the 5' and 3' termini, respectively, which showed conserved terminal sequences, as commonly observed for tospovirus genomic RNAs. The medium and small (S) RNAs were 4710 and 2485 nucleotides in size, respectively, and showed 99 % homology to the corresponding genomic segment of a previously partially characterized PolRSV isolate, Plg3. Protein sequences for GN/GC, N and NSs were identical in length in the two PolRSV isolates, while an amino acid insertion was observed for the NSm protein of the newly characterized isolate. The noncoding intergenic region of the S RNA was very short (183 nt) and was not predicted to form a hairpin structure, confirming that this unique characteristic within tospoviruses, previously observed for Plg3, is not isolate specific.

The arbuscular mycorrhizal symbiosis attenuates symptom severity and reduces virus concentration in tomato infected by Tomato yellow leaf curl Sardinia virus (TYLCSV)

The arbuscular mycorrhizal (AM) symbiosis is considered a natural instrument to improve plant health and productivity since mycorrhizal plants often show higher tolerance to abiotic and biotic stresses. However, the impact of the AM symbiosis on infection by viral pathogens is still largely uncertain and little explored. In the present study, tomato plants were grown under controlled conditions and inoculated with the AM fungus Funneliformis mosseae. Once the mycorrhizal colonization had developed, plants were inoculated with the Tomato yellow leaf curl Sardinia virus (TYLCSV), a geminivirus causing one of the most serious viral diseases of tomatoes in Mediterranean areas. Biological conditions consisted of control plants (C), TYLCSV-infected plants (V), mycorrhizal plants (M), and TYLCSV-infected mycorrhizal plants (MV). At the time of analysis, the level of mycorrhiza development and the expression profiles of mycorrhiza-responsive selected genes were not significantly modified by virus infection, thus indicating that the AM symbiosis was unaffected by the presence and spread of the virus. Viral symptoms were milder, and both shoot and root concentrations of viral DNA were lower in MV plants than in V plants. Overall F. mosseae colonization appears to exert a beneficial effect on tomato plants in attenuating the disease caused by TYLCSV.

From root to fruit: RNA-Seq analysis shows that arbuscular mycorrhizal symbiosis may affect tomato fruit metabolism

BACKGROUND

Tomato (Solanum lycopersicum) establishes a beneficial symbiosis with arbuscular mycorrhizal (AM) fungi. The formation of the mycorrhizal association in the roots leads to plant-wide modulation of gene expression. To understand the systemic effect of the fungal symbiosis on the tomato fruit, we used RNA-Seq to perform global transcriptome profiling on Moneymaker tomato fruits at the turning ripening stage.

RESULTS

Fruits were collected at 55 days after flowering, from plants colonized with Funneliformis mosseae and from control plants, which were fertilized to avoid responses related to nutrient deficiency. Transcriptome analysis identified 712 genes that are differentially expressed in fruits from mycorrhizal and control plants. Gene Ontology (GO) enrichment analysis of these genes showed 81 overrepresented functional GO classes. Up-regulated GO classes include photosynthesis, stress response, transport, amino acid synthesis and carbohydrate metabolism functions, suggesting a general impact of fungal symbiosis on primary metabolisms and, particularly, on mineral nutrition. Down-regulated GO classes include cell wall, metabolism and ethylene response pathways. Quantitative RT-PCR validated the RNA-Seq results for 12 genes out of 14 when tested at three fruit ripening stages, mature green, breaker and turning. Quantification of fruit nutraceutical and mineral contents produced values consistent with the expression changes observed by RNA-Seq analysis.

CONCLUSIONS

This RNA-Seq profiling produced a novel data set that explores the intersection of mycorrhization and fruit development. We found that the fruits of mycorrhizal plants show two transcriptomic "signatures": genes characteristic of a climacteric fleshy fruit, and genes characteristic of mycorrhizal status, like phosphate and sulphate transporters. Moreover, mycorrhizal plants under low nutrient conditions produce fruits with a nutrient content similar to those from non-mycorrhizal plants under high nutrient conditions, indicating that AM fungi can help replace exogenous fertilizer for fruit crops.

Transcriptomics of the interaction between the monopartite phloem-limited geminivirus tomato yellow leaf curl Sardinia virus and Solanum lycopersicum highlights a role for plant hormones, autophagy and plant immune system fine tuning during infection

Tomato yellow leaf curl Sardinia virus (TYLCSV), a DNA virus belonging to the genus Begomovirus, causes severe losses in tomato crops. It infects only a limited number of cells in the vascular tissues, making difficult to detect changes in host gene expression linked to its presence. Here we present the first microarray study of transcriptional changes induced by the phloem-limited geminivirus TYLCSV infecting tomato, its natural host. The analysis was performed on the midrib of mature leaves, a material naturally enriched in vascular tissues. A total of 2206 genes were up-regulated and 1398 were down-regulated in infected plants, with an overrepresentation of genes involved in hormone metabolism and responses, nucleic acid metabolism, regulation of transcription, ubiquitin-proteasome pathway and autophagy among those up-regulated, and in primary and secondary metabolism, phosphorylation, transcription and methylation-dependent chromatin silencing among those down-regulated. Our analysis showed a series of responses, such as the induction of GA- and ABA-responsive genes, the activation of the autophagic process and the fine tuning of the plant immune system, observed only in TYLCSV-tomato compatible interaction so far. On the other hand, comparisons with transcriptional changes observed in other geminivirus-plant interactions highlighted common host responses consisting in the deregulation of biotic stress responsive genes, key enzymes in the ethylene biosynthesis and methylation cycle, components of the ubiquitin proteasome system and DNA polymerases II. The involvement of conserved miRNAs and of solanaceous- and tomato-specific miRNAs in geminivirus infection, investigated by integrating differential gene expression data with miRNA targeting data, is discussed.

Analysis of small RNAs derived from tomato yellow leaf curl Sardinia virus reveals a cross reaction between the major viral hotspot and the plant host genome

RNA silencing is a defense mechanism exploited by plants against viruses. Upon infection, viral genomes and their transcripts are processed by Dicer-like (DCL) ribonucleases into viral small interfering RNAs (vsRNAs) of 21-24 nucleotides that further guide silencing of viral transcripts. To get an insight into the molecular interaction between tomato and the monopartite phloem-limited begomovirus tomato yellow leaf curl Sardinia virus (TYLCSV), a pathogen inducing a devastating disease of tomato in the Mediterranean region, we characterized by deep sequencing the vsRNA population in virus-infected tomato plants, using a Solexa/Illumina platform. TYLCSV-sRNAs spanned the entire viral genome but were discontinuously distributed throughout it, with a prevalence from the transcribed regions. TYLCSV-sRNAs were mainly 21-22 nucleotides in length and their polarity was asymmetrically distributed along the genome. The most abundant vsRNAs originated from a narrow region overlapping the Rep/C4 genes and from a broader region including the end of the V2 and the beginning of the coat protein genes. Deep sequencing results were validated by different hybridization techniques. Comparisons with the data available on vsRNAs for other begomoviruses highlighted both similarities and differences. Host-derived RNA species cross-reacting with a portion of the viral genome corresponding to the most abundant vsRNAs hotspot were detected. Bioinformatics analyses were carried out to investigate the nature of these host molecules.

Genome-wide identification of viral and host transcripts targeted by viral siRNAs in Vitis vinifera

In plants, RNA silencing is a surveillance mechanism against invading viruses. It involves the production of virus-derived small interfering RNAs (vsiRNAs), which guide the RNA-induced silencing complex (RISC) to inactivate viruses. vsiRNAs may also promote the silencing of host mRNAs in a sequence-specific manner. In this work, vsiRNAs derived from two grapevine-infecting viruses (Grapevine fleck virus and Grapevine rupestris stem pitting-associated virus) were selected from cDNA libraries of short RNAs and were cross-referenced with the remnants of both cleaved host transcripts and viral RNAs from a degradome dataset. We identified dozens of host transcripts targeted by vsiRNAs. Among them, several encode putative proteins involved in ribosome biogenesis and in biotic and abiotic stresses. Moreover, we identified vsiRNAs which explain the cleavage sites in viral genomes. A consistent fraction of vsiRNAs did not apparently account for cleavage, suggesting that only a low percentage of vsiRNAs are involved in the antiviral response.

Recombination profiles between Tomato yellow leaf curl virus and Tomato yellow leaf curl Sardinia virus in laboratory and field conditions: evolutionary and taxonomic implications

Tomato yellow leaf curl Sardinia virus and Tomato yellow leaf curl virus have co-existed in Italian tomato crops since 2002 and have reached equilibrium, with plants hosting molecules of both species plus their recombinants being the most frequent case. Recombination events are studied in field samples, as well as in experimental co-infections, when recombinants were detected as early as 45 days following inoculation. In both conditions, recombination breakpoints were essentially absent in regions corresponding to ORFs V2, CP and C4, whereas density was highest in the 3′-terminal portion of ORF C3, next to the region where the two transcription units co-terminate. The vast majority of breakpoints were mapped at antisense ORFs, supporting speculation that the rolling-circle replication mechanism, and the existence of sense and antisense ORFs on the circular genome, may result in clashes between replication and transcription complexes.

Arbuscular mycorrhizal symbiosis limits foliar transcriptional responses to viral infection and favors long-term virus accumulation

Tomato (Solanum lycopersicum) can establish symbiotic interactions with arbuscular mycorrhizal (AM) fungi, and can be infected by several pathogenic viruses. Here, we investigated the impact of mycorrhization by the fungus Glomus mosseae on the Tomato spotted wilt virus (TSWV) infection of tomato plants by transcriptomic and hormones level analyses. In TSWV-infected mycorrhizal plants, the AM fungus root colonization limited virus-induced changes in gene expression in the aerial parts. The virus-responsive upregulated genes, no longer induced in infected mycorrhizal plants, were mainly involved in defense responses and hormone signaling, while the virus-responsive downregulated genes, no longer repressed in mycorrhizal plants, were involved in primary metabolism. The presence of the AM fungus limits, in a salicylic acid-independent manner, the accumulation of abscissic acid observed in response to viral infection. At the time of the molecular analysis, no differences in virus concentration or symptom severity were detected between mycorrhizal and nonmycorrhizal plants. However, in a longer period, increase in virus titer and delay in the appearance of recovery were observed in mycorrhizal plants, thus indicating that the plant's reaction to TSWV infection is attenuated by mycorrhization.

Deep sequencing analysis of viral short RNAs from an infected Pinot Noir grapevine

Virus-derived short interfering RNAs (vsiRNAs) isolated from grapevine V. vinifera Pinot Noir clone ENTAV 115 were analyzed by high-throughput sequencing using the Illumina Solexa platform. We identified and characterized vsiRNAs derived from grapevine field plants naturally infected with different viruses belonging to the genera Foveavirus, Maculavirus, Marafivirus and Nepovirus. These vsiRNAs were mainly of 21 and 22 nucleotides (nt) in size and were discontinuously distributed throughout Grapevine rupestris stem-pitting associated virus (GRSPaV) and Grapevine fleck virus (GFkV) genomic RNAs. Among the studied viruses, GRSPaV and GFkV vsiRNAs had a 5' terminal nucleotide bias, which differed from that described for experimental viral infections in Arabidopsis thaliana. VsiRNAs were found to originate from both genomic and antigenomic GRSPaV RNA strands, whereas with the grapevine tymoviruses GFkV and Grapevine Red Globe associated virus (GRGV), the large majority derived from the antigenomic viral strand, a feature never observed in other plant-virus interactions.

The complete nucleotide sequence of an isolate of Tomato yellow leaf curl Sardinia virus found in Sicily

Partial sequences of Tomato yellow leaf curl Sardinia virus (TYLCSV) derived from tomato samples collected in Sicily in 1999, 2002 and 2004 indicated the presence of a TYLCSV different from the one previously described as the Sic strain. Here, we report a complete DNA sequence that is classified as belonging to the TYLCSV type strain (Sar strain), confirming the co-existence in Sicily of virus populations of both strains. Moreover, comparisons between this new sequence and those of the two recombinants recently described in Sicily revealed unequivocally (99% identity) that their TYLCSV-derived portion originated from the Sar strain.

ORTom: a multi-species approach based on conserved co-expression to identify putative functional relationships among genes in tomato

Co-expressed genes are often expected to be functionally related and many bioinformatics approaches based on co-expression have been developed to infer their biological role. However, such annotations may be unreliable, whereas the evolutionary conservation of gene co-expression among species may form a basis for more confident predictions. The huge amount of expression data (microarrays, SAGE, ESTs) has already allowed functional studies based on conserved co-expression in animals. Up to now, the implementation of analogous tools for plants has been strongly limited probably by the paucity and heterogeneity of data. Here we present ORTom, a tomato-centred EST data-mining approach based on conserved co-expression in the Solanaceae family. ORTom can be used to predict functional relationships among genes and to prioritize candidate genes for targeted studies. The method consists in ranking ESTs co-expressed with a gene of interest according to the level of expression pattern conservation in phylogenetically-related plants (potato, tobacco and pepper) to obtain lists of putative functionally-related genes. The lists are then analyzed for Gene Ontology keyword enrichment. The web server ORTom has been implemented to make the results publicly-available and searchable. Few biological examples on how the tool can be used are presented.

Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis

In plants, microRNAs (miRNAs) comprise one of three classes of small RNAs regulating gene expression at the post-transcriptional level. Many plant miRNAs are conserved, and play a role in development, abiotic stress responses or pathogen responses. However, some miRNAs have only been found in certain species. Here, we use deep-sequencing, computational and molecular methods to identify, profile, and describe conserved and non-conserved miRNAs in four grapevine (Vitis vinifera) tissues. A total of 24 conserved miRNA families were identified in all four tissues, and 26 known but non-conserved miRNAs were also found. In addition to known miRNAs, we also found 21 new grapevine-specific miRNAs together with their star strands. We have also shown that almost all of them originated from single genes. Furthermore, 21 other plausible miRNA candidates have been described. We have found that many known and new miRNAs showed tissue-specific expression. Finally, 112 target mRNAs of known and 44 target mRNAs of new grapevine-specific miRNAs were identified by genomic-scale high-throughput sequencing of miRNA cleaved mRNAs.

Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis

In plants, microRNAs (miRNAs) comprise one of three classes of small RNAs regulating gene expression at the post-transcriptional level. Many plant miRNAs are conserved, and play a role in development, abiotic stress responses or pathogen responses. However, some miRNAs have only been found in certain species. Here, we use deep-sequencing, computational and molecular methods to identify, profile, and describe conserved and non-conserved miRNAs in four grapevine (Vitis vinifera) tissues. A total of 24 conserved miRNA families were identified in all four tissues, and 26 known but non-conserved miRNAs were also found. In addition to known miRNAs, we also found 21 new grapevine-specific miRNAs together with their star strands. We have also shown that almost all of them originated from single genes. Furthermore, 21 other plausible miRNA candidates have been described. We have found that many known and new miRNAs showed tissue-specific expression. Finally, 112 target mRNAs of known and 44 target mRNAs of new grapevine-specific miRNAs were identified by genomic-scale high-throughput sequencing of miRNA cleaved mRNAs.

Comparative analysis of expression profiles in shoots and roots of tomato systemically infected by Tomato spotted wilt virus reveals organ-specific transcriptional responses

Tomato (Solanum lycopersicon), a model species for the family Solanaceae, is severely affected by Tomato spotted wilt virus (TSWV) worldwide. To elucidate the systemic transcriptional response of plants to TSWV infection, microarray experiments were performed on tomato. Parallel analysis of both shoots and roots revealed organ-specific responses, although the virus was present in similar concentration. In the shoots, genes related to defense and to signal transduction were induced, while there was general repression of genes related to primary and secondary metabolism as well as to amino acid metabolism. In roots, expression of genes involved in primary metabolism and signal transduction appear unaffected by TSWV infection, while those related to the response to biotic stimuli were induced and those associated to the response to abiotic stress were generally repressed or unaltered. Genes related to amino acid metabolism were unaffected, except for those involved in synthesis of secondary compounds, where induction was evident. Differential expression of genes involved in metabolism and response to ethylene and abscisic acid was observed in the two organs. Our results provide new insight into the biology of the economically important interaction between tomato and TSWV.

Global and cell-type gene expression profiles in tomato plants colonized by an arbuscular mycorrhizal fungus

*Arbuscular mycorrhizal symbiosis develops in roots; extensive cellular reorganizations and specific metabolic changes occur, which are mirrored by local and systemic changes in the transcript profiles. *A TOM2 microarray (c. 12 000 probes) has been used to obtain an overview of the transcriptional changes that are triggered in Solanum lycopersicum roots and shoots, as a result of colonization by the arbuscular mycorrhizal fungus Glomus mosseae. The cell-type expression profile of a subset of genes was monitored, using laser microdissection, to identify possible plant determinants of arbuscule development,. *Microarrays revealed 362 up-regulated and 293 down-regulated genes in roots. Significant gene modulation was also observed in shoots: 85 up- and 337 down-regulated genes. The most responsive genes in both organs were ascribed to primary and secondary metabolism, defence and response to stimuli, cell organization and protein modification, and transcriptional regulation. Six genes, preferentially expressed in arbusculated cells, were identified. *A comparative analysis only showed a limited overlap with transcript profiles identified in mycorrhizal roots of Medicago truncatula, probably as a consequence of the largely nonoverlapping probe sets on the microarray tools used. The results suggest that auxin and abscisic acid metabolism are involved in arbuscule formation and/or functioning.

Phospholipase A2 up-regulation during mycorrhiza formation in Tuber borchii

TbSP1 is a secreted and surface-associated phospholipase A(2) previously found to be up-regulated in C- or N-deprived free-living mycelia from the ectomycorrhizal ascomycete Tuber borchii. As nutrient limitation is considered an important environmental factor favouring the transition to symbiotic status, TbSP1 was suggested to be involved in the formation of mycorrhizas. An in vitro symbiosis system between Cistus incanus and T. borchii was set up: TbSP1 mRNA levels in free-living mycelia and in mycorrhizas sampled in different districts of the plant-fungus interaction were examined. In the same samples, TbSP1 protein expression was analysed by immunoelectron microscopy. A substantially enhanced TbSP1 mRNA expression, compared with nutrient-limited but free-living mycelia, was detected in the presence of the plant and reached maximal levels in fully developed mycorrhizas. A similar expression trend was revealed by immunolocalization experiments. We have shown that TbSP1 appears to respond to two partially overlapping yet distinct stimuli: nutrient starvation and mycorrhiza formation.