Transcriptome analysis of Cucumis sativus infected by Cucurbit chlorotic yellows virus

Insights into the early transcriptomic response against watermelon mosaic virus in melon

BACKGROUND

Watermelon mosaic virus (WMV) is one of the most prevalent viruses affecting melon worldwide. Recessive resistance to WMV in melon has previously been reported in the African accession TGR-1551. Moreover, the genomic regions associated to the resistance have also been described. Nevertheless, the transcriptomic response that might infer the resistance to this potyvirus has not been explored.

RESULTS

We have performed a comparative transcriptomic analysis using mock and WMV-inoculated plants of the susceptible cultivar "Bola de oro" (BO) and a resistant RIL (Recombinant inbred line) derived from the initial cross between "TGR-1551" and BO. In total, 616 genes were identified as differentially expressed and the weighted gene co-expression network analysis (WGCNA) detected 19 gene clusters (GCs), of which 7 were differentially expressed for the genotype x treatment interaction term. SNPs with a predicted high impact on the protein function were detected within the coding regions of most of the detected DEGs. Moreover, 3 and 16 DEGs were detected within the QTL regions previously described in chromosomes 11 and 5, respectively. In addition to these two specific genomic regions, we also observde large transcriptomic changes from genes spread across the genome in the resistant plants in response to the virus infection. This early response against WMV implied genes involved in plant-pathogen interaction, plant hormone signal transduction, the MAPK signaling pathway or ubiquitin mediated proteolysis, in detriment to the photosynthetic and basal metabolites pathways. Moreover, the gene MELO3C021395, which coded a mediator of RNA polymerase II transcription subunit 33A (MED33A), has been proposed as the candidate gene located on chromosome 11 conferring resistance to WMV.

CONCLUSIONS

The comparative transcriptomic analysis presented here showed that, even though the resistance to WMV in TGR-1551 has a recessive nature, it triggers an active defense response at a transcriptomic level, which involves broad-spectrum resistance mechanisms. Thus, this study represents a step forward on our understanding of the mechanisms underlaying WMV resistance in melon. In addition, it sheds light into a broader topic on the mechanisms of recessive resistances.

Comparative transcriptome profiling of potato cultivars infected by late blight pathogen Phytophthora infestans: Diversity of quantitative and qualitative responses

The Estonia potato cultivar Ando has shown elevated field resistance to Phytophthora infestans, even after being widely grown for over 40 years. A comprehensive transcriptional analysis was performed using RNA-seq from plant leaf tissues to gain insight into the mechanisms activated for the defense after infection. Pathogen infection in Ando resulted in about 5927 differentially expressed genes (DEGs) compared to 1161 DEGs in the susceptible cultivar Arielle. The expression levels of genes related to plant disease resistance such as serine/threonine kinase activity, signal transduction, plant-pathogen interaction, endocytosis, autophagy, mitogen-activated protein kinase (MAPK), and others were significantly enriched in the upregulated DEGs in Ando, whereas in the susceptible cultivar, only the pathway related to phenylpropanoid biosynthesis was enriched in the upregulated DEGs. However, in response to infection, photosynthesis was deregulated in Ando. Multi-signaling pathways of the salicylic-jasmonic-ethylene biosynthesis pathway were also activated in response to Phytophthora infestans infection.

RNA-seq analysis reveals an early defense response to tomato leaf curl New Delhi virus in potato cultivar Kufri Bahar

Potatoes in India are very susceptible to apical leaf curl disease, which causes severe symptoms and greater yield losses. Because the majority of potato cultivars are susceptible to the virus, it is crucial to discover sources of resistance and investigate the mechanism of resistance/susceptibility in potato cultivars. In this study, the gene expression profile of two potato cultivars, Kufri Bahar (resistant) and Kufri Pukhraj (susceptible), varying in their level of resistance to ToLCNDV, was analyzed using RNA-Seq. The Ion ProtonTM system was used to sequence eight RiboMinus RNA libraries from inoculated and uninoculated potato plants at 15 and 20 days after inoculation (DAI). The findings indicated that the majority of differentially expressed genes (DEGs) were cultivar-or time-specific. These DEGs included genes for proteins that interact with viruses, genes linked with the cell cycle, genes for proteins involved in defense, transcription and translation initiation factors, and plant hormone signaling pathway genes. Interestingly, defense responses were generated early in Kufri Bahar, at 15 DAI, which may have impeded the replication and spread of ToLCNDV. This research provides a genome-wide transcriptional analysis of two potato cultivars with variable levels of ToLCNDV resistance. At an early stage, we observed suppression of genes that interact with viral proteins, induction of genes associated with restriction of cell division, genes encoding defense proteins, AP2/ERF transcription factors, and altered expression of zinc finger protein genes, HSPs, JA, and SA pathway-related genes. Our findings add to a greater comprehension of the molecular basis of potato resistance to ToLCNDV and may aid in the development of more effective disease management techniques.

Transcriptome mapping related genes encoding PR1 protein involved in necrotic symptoms to soybean mosaic virus infection

Necrosis caused by soybean mosaic virus (SMV) has not been specifically distinguished from susceptible symptoms. The molecular mechanism for the occurrence of necrosis is largely overlooked in soybean genetic research. Field evaluation reveals that SMV disease seriously influences soybean production as indicated by decreasing 22.4% ~ 77.0% and 8.8% ~ 17.0% of yield and quality production, respectively. To expand molecular mechanism behind necrotic reactions, transcriptomic data obtained from the asymptomatic, mosaic, and necrotic pools were assessed. Compared between asymptomatic and mosaic plants, 1689 and 1752 up- and down-regulated differentially expressed genes (DEGs) were specifically found in necrotic plants. Interestingly, the top five enriched pathways with up-regulated DEGs were highly related to the process of the stress response, whereas the top three enriched pathways with down-regulated DEGs were highly related to the process of photosynthesis, demonstrating that defense systems are extensively activated, while the photosynthesis systems were severely destroyed. Further, results of the phylogenetic tree based on gene expression pattern and an amino acid sequence and validation experiments discovered three PR1 genes, Glyma.15G062400, Glyma.15G062500, and Glyma.15G062700, which were especially expressed in necrotic leaves. Meanwhile, exogenous salicylic acid (SA) but not methyl jasmonate (MeJA) could induce the three PR1 gene expressions on healthy leaves. Contrastingly, exogenous SA obviously decreased the expression level of Glyma.15G062400, Glyma.15G062500, and concentration of SMV, but increased Glyma.15G062700 expression in necrotic leaves. These results showed that GmPR1 is associated with the development of SMV-induced necrotic symptoms in soybean. Glyma.15G062400, Glyma.15G062500, and Glyma.15G062700 is up-regulated in necrotic leaves at the transcriptional levels, which will greatly facilitate a better understanding of the mechanism behind necrosis caused by SMV disease.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-022-01351-3.

Transcriptional Analysis of the Differences between ToLCNDV-India and ToLCNDV-ES Leading to Contrary Symptom Development in Cucumber

Tomato leaf curl New Delhi virus-ES (ToLCNDV-ES), a high threat to cucurbits in the Mediterranean Basin, is listed as a different strain from the Asian ToLCNDV isolates. In this study, the infectivity of two clones previously isolated from Italy and Pakistan were compared in cucumbers, which resulted in the opposite symptom appearance. The swapping subgenome was processed; however, the mechanisms related to the disease phenotype remain unclear. To identify the disease-associated genes that could contribute to symptom development under the two ToLCNDV infections, the transcriptomes of ToLCNDV-infected and mock-inoculated cucumber plants were compared 21 days postinoculation. The number of differentially expressed genes in ToLCNDV-India-infected plants was 10 times higher than in ToLCNDV-ES-infected samples. The gene ontology (GO) and pathway enrichment were analyzed using the Cucurbits Genomics Database. The flavonoid pathway-related genes were upregulated in ToLCNDV-ES, but some were downregulated in ToLCNDV-India infection, suggesting their role in resistance to the two ToLCNDV infections. The relative expression levels of the selected candidate genes were validated by qRT-PCR under two ToLCNDV-infected conditions. Our results reveal the different infectivity of the two ToLCNDVs in cucumber and also provide primary information based on RNA-seq for further analysis related to different ToLCNDV infections.

Reverse transcription recombinase polymerase amplification assay for rapid detection of the cucurbit chlorotic yellows virus

The cucurbit chlorotic yellows virus (CCYV) causes severe economic losses in cucurbit plants. Although it has been widely known in various countries for several years, CCYV is rarely recognized due to the lack of rapid and effective detection methods in the early stage of the disease. Recombinase polymerase amplification (RPA) is a new, efficient, and simple technology for nucleic acid detection. In the present study, reverse transcription (RT)-RPA and quantitative RT-RPA were developed and utilized for fast detection of CCYV in field-collected melon samples. The analysis was performed under constant temperature conditions without the necessity for a thermal cycler in just 20 min. Moreover, the detection limit of RT-RPA for CCYV was determined at 10 pg. In the study, 58 field-collected samples were employed to evaluate the performance of the two assays. The positive rates were established at 72.4 % (42/58) and 75.9 % (44/58) by RT-RPA and qRT-RPA, respectively, and were consistent with the RT-PCR results. The successful application of RPA for the detection of CCYV in field-collected melon samples indicated its potential applicability. Thus, the developed RPA assays provide an alternative for fast, efficient, sensitive, and reliable detection of CCYV in diagnostic laboratories, which lack the precise instrumentation, and fields without appropriate equipment.

Severe Stunting Symptoms upon Nepovirus Infection Are Reminiscent of a Chronic Hypersensitive-like Response in a Perennial Woody Fruit Crop

Virus infection of plants can result in various degrees of detrimental impacts and disparate symptom types and severities. Although great strides have been made in our understanding of the virus-host interactions in herbaceous model plants, the mechanisms underlying symptom development are poorly understood in perennial fruit crops. Grapevine fanleaf virus (GFLV) causes variable symptoms in most vineyards worldwide. To better understand GFLV-grapevine interactions in relation to symptom development, field and greenhouse trials were conducted with a grapevine genotype that exhibits distinct symptoms in response to a severe and a mild strain of GFLV. After validation of the infection status of the experimental vines by high-throughput sequencing, the transcriptomic and metabolomic profiles in plants infected with the two viral strains were tested and compared by RNA-Seq and LC-MS, respectively, in the differentiating grapevine genotype. In vines infected with the severe GFLV strain, 1023 genes, among which some are implicated in the regulation of the hypersensitive-type response, were specifically deregulated, and a higher accumulation of resveratrol and phytohormones was observed. Interestingly, some experimental vines restricted the virus to the rootstock and remained symptomless. Our results suggest that GFLV induces a strain- and cultivar-specific defense reaction similar to a hypersensitive reaction. This type of defense leads to a severe stunting phenotype in some grapevines, whereas others are resistant. This work is the first evidence of a hypersensitive-like reaction in grapevine during virus infection.

Transcriptome Profiling of Cucumber (Cucumis sativus L.) Early Response to Pseudomonas syringae pv. lachrymans

Bacterial angular leaf spot disease (ALS) caused by Pseudomonas syringae pv. lachrymans (Psl) is one of the biological factors limiting cucumber open-field production. The goal of this study was to characterize cytological and transcriptomic response of cucumber to this pathogen. Plants of two inbred lines, B10 (susceptible) and Gy14 (resistant), were grown, and leaves were inoculated with highly virulent Psl strain 814/98 under growth chamber conditions. Microscopic and transcriptional evaluations were performed at three time points: before, 1 and 3 days post inoculation (dpi). Investigated lines showed distinct response to Psl. At 1 dpi bacterial colonies were surrounded by necrotized mesophyll cells. At 3 dpi, in the susceptible B10 line bacteria were in contact with degraded cells, whereas cells next to bacteria in the resistant Gy14 line were plasmolyzed, but apparently still alive and functional. Additionally, the level of H₂O₂ production was higher in resistant Gy14 plants than in B10 at both examined time points. In RNA sequencing more than 18,800 transcripts were detected in each sample. As many as 1648 and 2755 differentially expressed genes (DEGs) at 1 dpi as well as 2992 and 3141 DEGs at 3 dpi were identified in B10 and Gy14, respectively. DEGs were characterized in terms of functional categories. Resistant line Gy14 showed massive transcriptomic response to Psl at 1 dpi compared to susceptible line B10, while a similar number of DEGs was detected for both lines at 3 dpi. This suggests that dynamic transcriptomic response to the invading pathogen may be related with host resistance. This manuscript provides the first transcriptomic data on cucumber infected with the pathovar lachrymans and helps to elucidate resistance mechanism against ALS disease.

Comparative RNA-Seq analysis unfolds a complex regulatory network imparting yellow mosaic disease resistance in mungbean [Vigna radiata (L.) R. Wilczek]

Yellow Mosaic Disease (YMD) in mungbean [Vigna radiata (L.) R. Wilczek] is one of the most damaging diseases in Asia. In the northern part of India, the YMD is caused by Mungbean Yellow Mosaic India Virus (MYMIV), while in southern India this is caused by Mungbean Yellow Mosaic Virus (MYMV). The molecular mechanism of YMD resistance in mungbean remains largely unknown. In this study, RNA-seq analysis was conducted between a resistant (PMR-1) and a susceptible (Pusa Vishal) mungbean genotype under infected and control conditions to understand the regulatory network operating between mungbean-YMV. Overall, 76.8 million raw reads could be generated in different treatment combinations, while mapping rate per library to the reference genome varied from 86.78% to 93.35%. The resistance to MYMIV showed a very complicated gene network, which begins with the production of general PAMPs (pathogen-associated molecular patterns), then activation of various signaling cascades like kinases, jasmonic acid (JA) and brassinosteroid (BR), and finally the expression of specific genes (like PR-proteins, virus resistance and R-gene proteins) leading to resistance response. The function of WRKY, NAC and MYB transcription factors in imparting the resistance against MYMIV could be established. The string analysis also revealed the role of proteins involved in kinase, viral movement and phytoene synthase activity in imparting YMD resistance. A set of novel stress-related EST-SSRs are also identified from the RNA-Seq data which may be used to find the linked genes/QTLs with the YMD resistance. Also, 11 defence-related transcripts could be validated through quantitative real-time PCR analysis. The identified gene networks have led to an insight about the defence mechanism operating against MYMIV infection in mungbean which will be of immense use to manage the YMD resistance in mungbean.

RNA-Seq Transcriptome Analysis Provides Candidate Genes for Resistance to Tomato Leaf Curl New Delhi Virus in Melon

Tomato leaf curl New Delhi virus (ToLCNDV) emerged in the Mediterranean Basin in 2012 as the first DNA bipartite begomovirus (Geminiviridae family), causing severe yield and economic losses in cucurbit crops. A major resistance locus was identified in the wild melon accession WM-7 (Cucumis melo kachri group), but the mechanisms involved in the resistant response remained unknown. In this work, we used RNA-sequencing to identify disease-associated genes that are differentially expressed in the course of ToLCNDV infection and could contribute to resistance. Transcriptomes of the resistant WM-7 genotype and the susceptible cultivar Piñonet Piel de Sapo (PS) (C. melo ibericus group) in ToLCNDV and mock inoculated plants were compared at four time points during infection (0, 3, 6, and 12 days post inoculation). Different gene expression patterns were observed over time in the resistant and susceptible genotypes in comparison to their respective controls. Differentially expressed genes (DEGs) in ToLCNDV-infected plants were classified using gene ontology (GO) terms, and genes of the categories transcription, DNA replication, and helicase activity were downregulated in WM-7 but upregulated in PS, suggesting that reduced activity of these functions reduces ToLCNDV replication and intercellular spread and thereby contributes to resistance. DEGs involved in the jasmonic acid signaling pathway, photosynthesis, RNA silencing, transmembrane, and sugar transporters entail adverse consequences for systemic infection in the resistant genotype, and lead to susceptibility in PS. The expression levels of selected candidate genes were validated by qRT-PCR to corroborate their differential expression upon ToLCNDV infection in resistant and susceptible melon. Furthermore, single nucleotide polymorphism (SNPs) with an effect on structural functionality of DEGs linked to the main QTLs for ToLCNDV resistance have been identified. The obtained results pinpoint cellular functions and candidate genes that are differentially expressed in a resistant and susceptible melon line in response to ToLCNDV, an information of great relevance for breeding ToLCNDV-resistant melon cultivars.

Identification of Differentially Expressed Proteins in Sugarcane in Response to Infection by Xanthomonas albilineans Using iTRAQ Quantitative Proteomics

Sugarcane can suffer severe yield losses when affected by leaf scald, a disease caused by Xanthomonas albilineans. This bacterial pathogen colonizes the vascular system of sugarcane, which can result in reduced plant growth and plant death. In order to better understand the molecular mechanisms involved in the resistance of sugarcane to leaf scald, a comparative proteomic study was performed with two sugarcane cultivars inoculated with X. albilineans: one resistant (LCP 85-384) and one susceptible (ROC20) to leaf scald. The iTRAQ (isobaric tags for relative and absolute quantification) approach at 0 and 48 h post-inoculation (hpi) was used to identify and annotate differentially expressed proteins (DEPs). A total of 4295 proteins were associated with 1099 gene ontology (GO) terms by GO analysis. Among those, 285 were DEPs during X. albilineans infection in cultivars LCP 85-384 and ROC20. One hundred seventy-two DEPs were identified in resistant cultivar LCP 85-384, and 113 of these proteins were upregulated and 59 were downregulated. One hundred ninety-two DEPs were found in susceptible cultivar ROC20 and half of these (92) were upregulated, whereas the other half corresponded to downregulated proteins. The significantly upregulated DEPs in LCP 85-384 were involved in metabolic pathways, the biosynthesis of secondary metabolites, and the phenylpropanoid biosynthesis pathway. Additionally, the expression of seven candidate genes related to photosynthesis and glycolytic pathways, plant innate immune system, glycosylation process, plant cytochrome P450, and non-specific lipid transfer protein was verified based on transcription levels in sugarcane during infection by X. albilineans. Our findings shed new light on the differential expression of proteins in sugarcane cultivars in response to infection by X. albilineans. The identification of these genes provides important information for sugarcane variety improvement programs using molecular breeding strategies.

Global Gene Responses of Resistant and Susceptible Sugarcane Cultivars to Acidovorax avenae subsp. avenae Identified Using Comparative Transcriptome Analysis

Red stripe disease in sugarcane caused by Acidovorax avenae subsp. avenae (Aaa) is related to serious global losses in yield. However, the underlying molecular mechanisms associated with responses of sugarcane plants to infection by this pathogen remain largely unknown. Here, we used Illumina RNA-sequencing (RNA-seq) to perform large-scale transcriptome sequencing of two sugarcane cultivars to contrast gene expression patterns of plants between Aaa and mock inoculations, and identify key genes and pathways involved in sugarcane defense responses to Aaa infection. At 0–72 hours post-inoculation (hpi) of the red stripe disease-resistant cultivar ROC22, a total of 18,689 genes were differentially expressed between Aaa-inoculated and mock-inoculated samples. Of these, 8498 and 10,196 genes were up- and downregulated, respectively. In MT11-610, which is susceptible to red stripe disease, 15,782 genes were differentially expressed between Aaa-inoculated and mock-inoculated samples and 8807 and 6984 genes were up- and downregulated, respectively. The genes that were differentially expressed following Aaa inoculation were mainly involved in photosynthesis and carbon metabolism, phenylpropanoid biosynthesis, plant hormone signal transduction, and plant–pathogen interaction pathways. Further, qRT-PCR and RNA-seq used for additional validation of 12 differentially expressed genes (DEGs) showed that eight genes in particular were highly expressed in ROC22. These eight genes participated in the biosynthesis of lignin and coumarin, as well as signal transduction by salicylic acid, jasmonic acid, ethylene, and mitogen-activated protein kinase (MAPK), suggesting that they play essential roles in sugarcane resistance to Aaa. Collectively, our results characterized the sugarcane transcriptome during early infection with Aaa, thereby providing insights into the molecular mechanisms responsible for bacterial tolerance.

Cucurbit Chlorotic Yellows Virus p22 Protein Interacts with Cucumber SKP1LB1 and Its F-Box-Like Motif Is Crucial for Silencing Suppressor Activity

Plants use RNA silencing as a defense against viruses. In response, viruses encode various RNA silencing suppressors to counteract the antiviral silencing. Here, we identified p22 as a silencing suppressor of cucurbit chlorotic yellows crinivirus and showed that p22 interacts with CsSKP1LB1, a Cucumis sativus ortholog of S-phase kinase-associated protein 1 (SKP1). The F-box-like motif of p22 was identified through sequence analysis and found to be necessary for the interaction using a yeast two-hybrid assay. The involvement of the F-box-like motif in p22 silencing suppressor activity was determined. Proteomics analysis of Nicotiana benthamiana leaves expressing p22, and its F-box-like motif deletion mutant showed 228 differentially expressed proteins and five enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways: ABC transporters, sesquiterpenoid and triterpenoid biosynthesis, ubiquitin-mediated proteolysis, riboflavin metabolism, and cysteine and methionine metabolism. Collectively, our results demonstrate the interaction between p22 and CsSKP1LB1 and show that the deletion of F-box-like motif inhibits p22 silencing suppressor activity. The possible pathways regulated by the p22 through the F-box-like motif were identified using proteomics analysis.