Integrated single-base resolution maps of transcriptome, sRNAome and methylome of Tomato yellow leaf curl virus (TYLCV) in tomato

In-depth analysis of 17,115 rice transcriptomes reveals extensive viral diversity in rice plants

Rice viruses seriously threaten rice cultivation and cause significant economic losses, but they have not yet been systematically identified, with only 20 rice-infecting viruses reported. Here, we perform a large-scale analysis of 17,115 RNA-seq libraries spanning 24 Oryza species across 51 countries. Using de novo assembly and homology-based methods, we identify 810 complete or near-complete viruses, including 276 known viruses and 534 novel viruses. Given the high divergence and atypical genome organizations of novel viruses, more than a half of them are tentatively assigned to 1 new order, 61 new families, and at least 104 new genera. Utilizing homology-independent approaches, we additionally identify 49 divergent RNA-dependent RNA polymerases (RdRPs), which are confirmed by protein structural alignment. Furthermore, we analyze the metadata of related Sequence Read Archive (SRA) libraries and estimated viral abundance in each library, leading to the screening of 427 viruses closely associated with rice plants. Overall, our study vastly expands the viral diversity in rice plants, providing insights for the prevention and control of viral disease.

Viral Recognition and Evasion in Plants

Viruses, causal agents of devastating diseases in plants, are obligate intracellular pathogens composed of a nucleic acid genome and a limited number of viral proteins. The diversity of plant viruses, their diminutive molecular nature, and their symplastic localization pose challenges to understanding the interplay between these pathogens and their hosts in the currently accepted framework of plant innate immunity. It is clear, nevertheless, that plants can recognize the presence of a virus and activate antiviral immune responses, although our knowledge of the breadth of invasion signals and the underpinning sensing events is far from complete. Below, I discuss some of the demonstrated or hypothesized mechanisms enabling viral recognition in plants, the step preceding the onset of antiviral immunity, as well as the strategies viruses have evolved to evade or suppress their detection.

Transcriptome and small RNAome profiling uncovers how a recombinant begomovirus evades RDRγ-mediated silencing of viral genes and outcompetes its parental virus in mixed infection

Tomato yellow leaf curl virus (TYLCV, genus Begomovirus, family Geminiviridae) causes severe disease of cultivated tomatoes. Geminiviruses replicate circular single-stranded genomic DNA via rolling-circle and recombination-dependent mechanisms, frequently generating recombinants in mixed infections. Circular double-stranded intermediates of replication also serve as templates for Pol II bidirectional transcription. IS76, a recombinant derivative of TYLCV with a short sequence in the bidirectional promoter/origin-of-replication region acquired from a related begomovirus, outcompetes TYLCV in mixed infection and breaks disease resistance in tomato Ty-1 cultivars. Ty-1 encodes a γ-clade RNA-dependent RNA polymerase (RDRγ) implicated in Dicer-like (DCL)-mediated biogenesis of small interfering (si)RNAs directing gene silencing. Here, we profiled transcriptome and small RNAome of Ty-1 resistant and control susceptible plants infected with TYLCV, IS76 or their combination at early and late infection stages. We found that RDRγ boosts production rates of 21, 22 and 24 nt siRNAs from entire genomes of both viruses and modulates DCL activities in favour of 22 and 24 nt siRNAs. Compared to parental TYLCV, IS76 undergoes faster transition to the infection stage favouring rightward transcription of silencing suppressor and coat protein genes, thereby evading RDRγ activity and facilitating its DNA accumulation in both single and mixed infections. In coinfected Ty-1 plants, IS76 efficiently competes for host replication and transcription machineries, thereby impairing TYLCV replication and transcription and forcing its elimination associated with further increased siRNA production. RDRγ is constitutively overexpressed in Ty-1 plants, which correlates with begomovirus resistance, while siRNA-generating DCLs (DCL2b/d, DCL3, DCL4) and genes implicated in siRNA amplification (α-clade RDR1) and function (Argonaute2) are upregulated to similar levels in TYLCV- and IS76-infected susceptible plants. Collectively, IS76 recombination facilitates replication and promotes expression of silencing suppressor and coat proteins, which allows the recombinant virus to evade the negative impact of RDRγ-boosted production of viral siRNAs directing transcriptional and posttranscriptional silencing.

Transcriptional and epigenetic changes during tomato yellow leaf curl virus infection in tomato

BACKGROUND

Geminiviruses are DNA plant viruses that cause highly damaging diseases affecting crops worldwide. During the infection, geminiviruses hijack cellular processes, suppress plant defenses, and cause a massive reprogramming of the infected cells leading to major changes in the whole plant homeostasis. The advances in sequencing technologies allow the simultaneous analysis of multiple aspects of viral infection at a large scale, generating new insights into the molecular mechanisms underlying plant-virus interactions. However, an integrative study of the changes in the host transcriptome, small RNA profile and methylome during a geminivirus infection has not been performed yet. Using a time-scale approach, we aim to decipher the gene regulation in tomato in response to the infection with the geminivirus, tomato yellow leaf curl virus (TYLCV).

RESULTS

We showed that tomato undergoes substantial transcriptional and post-transcriptional changes upon TYLCV infection and identified the main altered regulatory pathways. Interestingly, although the principal plant defense-related processes, gene silencing and the immune response were induced, this cannot prevent the establishment of the infection. Moreover, we identified extra- and intracellular immune receptors as targets for the deregulated microRNAs (miRNAs) and established a network for those that also produced phased secondary small interfering RNAs (phasiRNAs). On the other hand, there were no significant genome-wide changes in tomato methylome at 14 days post infection, the time point at which the symptoms were general, and the amount of viral DNA had reached its maximum level, but we were able to identify differentially methylated regions that could be involved in the transcriptional regulation of some of the differentially expressed genes.

CONCLUSION

We have conducted a comprehensive and reliable study on the changes at transcriptional, post-transcriptional and epigenetic levels in tomato throughout TYLCV infection. The generated genomic information is substantial for understanding the genetic, molecular and physiological changes caused by TYLCV infection in tomato.

Enhanced Susceptibility to Tomato Chlorosis Virus (ToCV) in Hsp90- and Sgt1-Silenced Plants: Insights from Gene Expression Dynamics

The emerging whitefly-transmitted crinivirus tomato chlorosis virus (ToCV) causes substantial economic losses by inducing yellow leaf disorder in tomato crops. This study explores potential resistance mechanisms by examining early-stage molecular responses to ToCV. A time-course transcriptome analysis compared naïve, mock, and ToCV-infected plants at 2, 7, and 14 days post-infection (dpi). Gene expression changes were most notable at 2 and 14 dpi, likely corresponding to whitefly feeding and viral infection. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed key genes and pathways associated with ToCV infection, including those related to plant immunity, flavonoid and steroid biosynthesis, photosynthesis, and hormone signaling. Additionally, virus-derived small interfering RNAs (vsRNAs) originating from ToCV predominantly came from RNA2 and were 22 nucleotides in length. Furthermore, two genes involved in plant immunity, Hsp90 (heat shock protein 90) and its co-chaperone Sgt1 (suppressor of the G2 allele of Skp1) were targeted through viral-induced gene silencing (VIGS), showing a potential contribution to basal resistance against viral infections since their reduction correlated with increased ToCV accumulation. This study provides insights into tomato plant responses to ToCV, with potential implications for developing effective disease control strategies.

Replication mechanisms of circular ssDNA plant viruses and their potential implication in viral gene expression regulation

The replication of members of the two circular single-stranded DNA (ssDNA) virus families Geminiviridae and Nanoviridae, the only ssDNA viruses infecting plants, is believed to be processed by rolling-circle replication (RCR) and recombination-dependent replication (RDR) mechanisms. RCR is a ubiquitous replication mode for circular ssDNA viruses and involves a virus-encoded Replication-associated protein (Rep) which fulfills multiple functions in the replication mechanism. Two key genomic elements have been identified for RCR in Geminiviridae and Nanoviridae: (i) short iterative sequences called iterons which determine the specific recognition of the viral DNA by the Rep and (ii) a sequence enabling the formation of a stem-loop structure which contains a conserved motif and constitutes the origin of replication. In addition, studies in Geminiviridae provided evidence for a second replication mode, RDR, which has also been documented in some double-stranded DNA viruses. Here, we provide a synthesis of the current understanding of the two presumed replication modes of Geminiviridae and Nanoviridae, and we identify knowledge gaps and discuss the possibility that these replication mechanisms could regulate viral gene expression through modulation of gene copy number.

Infection Defects of RNA and DNA Viruses Induced by Antiviral RNA Interference

Immune recognition of viral genome-derived double-stranded RNA (dsRNA) molecules and their subsequent processing into small interfering RNAs (siRNAs) in plants, invertebrates, and mammals trigger specific antiviral immunity known as antiviral RNA interference (RNAi). Immune sensing of viral dsRNA is sequence-independent, and most regions of viral RNAs are targeted by virus-derived siRNAs which extensively overlap in sequence. Thus, the high mutation rates of viruses do not drive immune escape from antiviral RNAi, in contrast to other mechanisms involving specific virus recognition by host immune proteins such as antibodies and resistance (R) proteins in mammals and plants, respectively. Instead, viruses actively suppress antiviral RNAi at various key steps with a group of proteins known as viral suppressors of RNAi (VSRs). Some VSRs are so effective in virus counter-defense that potent inhibition of virus infection by antiviral RNAi is undetectable unless the cognate VSR is rendered nonexpressing or nonfunctional. Since viral proteins are often multifunctional, resistance phenotypes of antiviral RNAi are accurately defined by those infection defects of VSR-deletion mutant viruses that are efficiently rescued by host deficiency in antiviral RNAi. Here, we review and discuss in vivo infection defects of VSR-deficient RNA and DNA viruses resulting from the actions of host antiviral RNAi in model systems.

The identification of the methylation patterns of tomato curly stunt virus in resistant and susceptible tomato lines

Tomato curly stunt virus (ToCSV) is a monopartite begomovirus infecting tomatoes in South Africa, with sequence similarity to tomato yellow leaf curl virus (TYLCV). While there are numerous reports on the mechanism of TYLCV resistance in tomato, the underlying mechanisms in the tomato-ToCSV pathosystem is still relatively unknown. The main aim of this study was to investigate and compare the global methylation profile of ToCSV in two near-isogenic tomato lines, one with a tolerant phenotype (T, NIL396) and one with a susceptible phenotype (S, NIL395). Bisulfite conversion and PCR amplification, coupled with a next-generation sequencing approach, were used to elucidate the global pattern of methylation of ToCSV cytosine residues in T and S leave tissue at 35 days post-infection (dpi). The extent of methylation was more pronounced in tolerant plants compared to susceptible plants in all sequence (CG, CHG and CHH) contexts, however, the overall methylation levels were relatively low (<3%). Notably, a significant interaction (p < 0.05) was observed between the viral genomic region and susceptible vs. tolerant status for CG methylated regions where it was observed that the 3'IR CG methylation was significantly (p < 0.05) higher than CG methylation of other genomic regions in tolerant and susceptible plants. Additionally, statistically significant (EdgeR p < 0.05) differentially methylated cytosines were located primarily in the genomic regions V2/V1 and C4/C1 of ToCSV. The relative expression, using RT-qPCR, was also employed in order to quantify the expression of various key methylation-related genes, MET1, CMT2, KYP4/SUVH4, DML2, RDM1, AGO4 and AGO6 in T vs. S plants at 35dpi. The differential expression between T and S was significant for MET1, KYP4/SUVH4 and RDM1 at p<0.05 which further supports more pronounced methylation observed in ToCSV from T plants vs. S plants. While this study provides new insights into the differences in methylation profiles of ToCSV in S vs. T tomato plants, further research is required to link tolerance and susceptibility to ToCSV.

Differences in the 3' intergenic region and the V2 protein of two sequence variants of tomato curly stunt virus play an important role in disease pathology in Nicotiana benthamiana

Tomato production in South Africa is threatened by the emergence of tomato curly stunt virus (ToCSV), a monopartite Begomovirus transmitted by the whitefly vector Bemisia tabaci (Genn.). We investigated the role of sequence differences present in the 3' intergenic region (IR) and the V2 coding region on the differing infectivity of ToCSV sequence variant isolates V30 and V22 in the model host Nicotiana benthamiana. Using virus mutant chimeras, we determined that the development of the upward leaf roll symptom phenotype is mediated by sequence differences present in the 3' IR containing the TATA-associated composite element. Sequence differences present in the V2 coding region are responsible for modulating disease severity and symptom recovery in V22-infected plants. Serine substitution of V22 V2 Val27 resulted in a significant increase in disease severity with reduced recovery, the first study to demonstrate the importance of this V2 residue in disease development. Two putative ORFs, C5 and C6, were identified using in silico analysis and detection of an RNA transcript spanning their coding region suggests that these ORFs may be transcribed during infection. Additional virus-derived RNA transcripts spanning multiple ORFs and crossing the boundaries of recognised polycistronic transcripts, as well as the origin of replication within the IR, were detected in ToCSV-infected plants providing evidence of bidirectional readthrough transcription. From our results, we conclude that the diverse responses of the model host to ToCSV infection is influenced by select sequence differences and our findings provide several avenues for further investigation into the mechanisms behind these responses to infection.

Differential Effects of RNA-Dependent RNA Polymerase 6 (RDR6) Silencing on New and Old World Begomoviruses in Nicotiana benthamiana

RNA-dependent RNA polymerases (RDRs) are key players in the antiviral defence mediated by RNA silencing in plants. RDR6 is one of the major components of the process, regulating the infection of certain RNA viruses. To better clarify its function against DNA viruses, we analyzed the effect of RDR6 inactivation (RDR6i) in N. benthamiana plants on two phloem-limited begomoviruses, the bipartite Abutilon mosaic virus (AbMV) and the monopartite tomato yellow leaf curl Sardinia virus (TYLCSV). We observed exacerbated symptoms and DNA accumulation for the New World virus AbMV in RDR6i plants, varying with the plant growth temperature (ranging from 16 °C to 33 °C). However, for the TYLCSV of Old World origin, RDR6 depletion only affected symptom expression at elevated temperatures and to a minor extent; it did not affect the viral titre. The accumulation of viral siRNA differed between the two begomoviruses, being increased in RDR6i plants infected by AbMV but decreased in those infected by TYLCSV compared to wild-type plants. In situ hybridization revealed a 6.5-fold increase in the number of AbMV-infected nuclei in RDR6i plants but without egress from the phloem tissues. These results support the concept that begomoviruses adopt different strategies to counteract plant defences and that TYLCSV evades the functions exerted by RDR6 in this host.

Soaud S, Huang Q, M. Saleh A, A. S. Abourehab M, Wan L, Cheng GT, Liu J, Ihtisham M, Noor Z, Rouf Mir R, Zhao X, Yan K, Abbas M, Li J. Natural resistance of tomato plants to Tomato yellow leaf curl virus

Tomato yellow leaf curl virus (TYLCV) is one of the most harmful afflictions in the world that affects tomato growth and production. Six regular antagonistic genes (Ty-1, Ty-2, Ty-3, Ty-4, ty-5, and Ty-6) have been transferred from wild germplasms to commercial cultivars as TYLCV protections. With Ty-1 serving as an appropriate source of TYLCV resistance, only Ty-1, Ty-2, and Ty-3 displayed substantial levels of opposition in a few strains. It has been possible to clone three TYLCV opposition genes (Ty-1/Ty-3, Ty-2, and ty-5) that target three antiviral safety mechanisms. However, it significantly impacts obtaining permanent resistance to TYLCV, trying to maintain opposition whenever possible, and spreading opposition globally. Utilizing novel methods, such as using resistance genes and identifying new resistance resources, protects against TYLCV in tomato production. To facilitate the breeders make an informed decision and testing methods for TYLCV blockage, this study highlights the portrayal of typical obstruction genes, common opposition sources, and subatomic indicators. The main goal is to provide a fictitious starting point for the identification and application of resistance genes as well as the maturation of tomato varieties that are TYLCV-resistant.

DNA methylation dynamics in response to abiotic and pathogen stress in plants

DNA methylation is a dynamic epigenetic mechanism that plays a significant role in gene expression and also maintains chromatin stability. The process is conserved in both plants and animals, and crucial for development and stress responses. Differential DNA methylation during adverse environmental conditions or pathogen attack facilitates the selective expression of defense-related genes. Both stress-induced DNA hypomethylation and hypermethylation play beneficial roles in activating the defense response. These DNA marks may be carried to the next generation making the progenies 'primed' for abiotic and biotic stress responses. Over the recent years, rapid advancements in the area of high throughput sequencing have enabled the detection of methylation status at genome levels in several plant species. Epigenotyping offers an alternative tool to plant breeders in addition to conventional markers for the selection of the desired offspring. In this review, we briefly discuss the mechanism of DNA methylation, recent understanding of DNA methylation-mediated gene regulation during abiotic and biotic stress responses, and stress memory in plants.

Role of main RNA modifications in cancer: N6-methyladenosine, 5-methylcytosine, and pseudouridine

Cancer is one of the major diseases threatening human life and health worldwide. Epigenetic modification refers to heritable changes in the genetic material without any changes in the nucleic acid sequence and results in heritable phenotypic changes. Epigenetic modifications regulate many biological processes, such as growth, aging, and various diseases, including cancer. With the advancement of next-generation sequencing technology, the role of RNA modifications in cancer progression has become increasingly prominent and is a hot spot in scientific research. This review studied several common RNA modifications, such as N⁶-methyladenosine, 5-methylcytosine, and pseudouridine. The deposition and roles of these modifications in coding and noncoding RNAs are summarized in detail. Based on the RNA modification background, this review summarized the expression, function, and underlying molecular mechanism of these modifications and their regulators in cancer and further discussed the role of some existing small-molecule inhibitors. More in-depth studies on RNA modification and cancer are needed to broaden the understanding of epigenetics and cancer diagnosis, treatment, and prognosis.

Geminiviruses employ host DNA glycosylases to subvert DNA methylation-mediated defense

DNA methylation is an epigenetic mechanism that plays important roles in gene regulation and transposon silencing. Active DNA demethylation has evolved to counterbalance DNA methylation at many endogenous loci. Here, we report that active DNA demethylation also targets viral DNAs, tomato yellow leaf curl China virus (TYLCCNV) and its satellite tomato yellow leaf curl China betasatellite (TYLCCNB), to promote their virulence. We demonstrate that the βC1 protein, encoded by TYLCCNB, interacts with a ROS1-like DNA glycosylase in Nicotiana benthamiana and with the DEMETER (DME) DNA glycosylase in Arabidopsis thaliana. The interaction between βC1 and DME facilitates the DNA glycosylase activity to decrease viral DNA methylation and promote viral virulence. These findings reveal that active DNA demethylation can be regulated by a viral protein to subvert DNA methylation-mediated defense.

Small RNA Profiling of Susceptible and Resistant Ty-1 Encoding Tomato Plants Upon Tomato Yellow Leaf Curl Virus Infection

Ty-1 presents an atypical dominant resistance gene that codes for an RNA-dependent RNA polymerase (RDR) of the gamma class and confers resistance to tomato yellow leaf curl virus (TYLCV) and other geminiviruses. Tomato lines bearing Ty-1 not only produce relatively higher amounts of viral small interfering (vsi)RNAs, but viral DNA also exhibits a higher amount of cytosine methylation. Whether Ty-1 specifically enhances posttranscriptional gene silencing (PTGS), leading to a degradation of RNA target molecules and primarily relying on 21-22 nucleotides (nts) siRNAs, and/or transcriptional gene silencing (TGS), leading to the methylation of cytosines within DNA target sequences and relying on 24-nts siRNAs, was unknown. In this study, small RNAs were isolated from systemically TYLCV-infected leaves of Ty-1 encoding tomato plants and susceptible tomato Moneymaker (MM) and sequence analyzed. While in susceptible tomato plants vsiRNAs of the 21-nt size class were predominant, their amount was drastically reduced in tomato containing Ty-1. The latter, instead, revealed elevated levels of vsiRNAs of the 22- and 24-nt size classes. In addition, the genomic distribution profiles of the vsiRNAs were changed in Ty-1 plants compared with those from susceptible MM. In MM three clear hotspots were seen, but these were less pronounced in Ty-1 plants, likely due to enhanced transitive silencing to neighboring viral genomic sequences. The largest increase in the amount of vsiRNAs was observed in the intergenic region and the V1 viral gene. The results suggest that Ty-1 enhances an antiviral TGS response. Whether the elevated levels of 22 nts vsiRNAs contribute to an enhanced PTGS response or an additional TGS response involving a noncanonical pathway of RNA dependent DNA methylation remains to be investigated.

Epigenetic regulation of geminivirus pathogenesis: a case of relentless recalibration of defence responses in plants

Geminiviruses constitute one of the largest families of plant viruses and they infect many economically important crops. The proteins encoded by the single-stranded DNA genome of these viruses interact with a wide range of host proteins to cause global dysregulation of cellular processes and help establish infection in the host. Geminiviruses have evolved numerous mechanisms to exploit host epigenetic processes to ensure the replication and survival of the viral genome. Here, we review our current knowledge of diverse epigenetic processes that have been implicated in the regulation of geminivirus pathogenesis, including DNA methylation, histone post-transcriptional modification, chromatin remodelling, and nucleosome repositioning. In addition, we discuss the currently limited evidence of host epigenetic defence responses that are aimed at counteracting geminivirus infection, and the potential for exploiting these responses for the generation of resistance against geminiviruses in crop species.

A virus-encoded protein suppresses methylation of the viral genome through its interaction with AGO4 in the Cajal body

In plants, establishment of de novo DNA methylation is regulated by the RNA-directed DNA methylation (RdDM) pathway. RdDM machinery is known to concentrate in the Cajal body, but the biological significance of this localization has remained elusive. Here, we show that the antiviral methylation of the Tomato yellow leaf curl virus (TYLCV) genome requires the Cajal body in Nicotiana benthamiana cells. Methylation of the viral genome is countered by a virus-encoded protein, V2, which interacts with the central RdDM component AGO4, interfering with its binding to the viral DNA; Cajal body localization of the V2-AGO4 interaction is necessary for the viral protein to exert this function. Taken together, our results draw a long sought-after functional connection between RdDM, the Cajal body, and antiviral DNA methylation, paving the way for a deeper understanding of DNA methylation and antiviral defences in plants.

Recent advances on the plant manipulation by geminiviruses

As intracellular parasites, viruses co-opt the molecular machinery of the cells they infect in order to multiply and spread, and the extensiveness and effectiveness of this manipulation ultimately determine the outcome of the interaction between virus and host. Members of the Geminiviridae family, causal agents of devastating diseases in crops, encode only a handful of multifunctional, fast-evolving proteins, which efficiently target host proteins to re-wire plant development and physiology and enable replication and spread of the viral genome. In this review, we offer an overview of the different steps in the geminiviral invasion of the host plant, and explore the knowns and unknowns in geminivirus biology.

The genome of the Cauliflower mosaic virus, a plant pararetrovirus, is highly methylated in the nucleus

Cytosine methylation is an important defense against invasive DNAs. Here, cytosine methylation profiles of a plant pararetrovirus, Cauliflower mosaic virus (CaMV), were investigated. Nuclear CaMV DNA is highly methylated throughout the genome including at transcription regulatory regions, but the virion DNA is unmethylated. In vitro CG methylation of the viral 35S promoter reduces transcription from the downstream gene. Although nuclear CaMV DNA is highly methylated, its transcripts are accumulated in the nucleus. The data suggest that a small population of unmethylated viral genomes produced through reverse transcription are constantly delivered back to the nucleus. Small RNA profiles suggest that methylation of the CaMV DNA may be due to de novo methylation through 21-, 22-, and 24-nt small RNAs with adenines at their 5' terminus.

Towards Exploitation of Adaptive Traits for Climate-Resilient Smart Pulses

Pulses are the main source of protein and minerals in the vegetarian diet. These are primarily cultivated on marginal lands with few inputs in several resource-poor countries of the world, including several in South Asia. Their cultivation in resource-scarce conditions exposes them to various abiotic and biotic stresses, leading to significant yield losses. Furthermore, climate change due to global warming has increased their vulnerability to emerging new insect pests and abiotic stresses that can become even more serious in the coming years. The changing climate scenario has made it more challenging to breed and develop climate-resilient smart pulses. Although pulses are climate smart, as they simultaneously adapt to and mitigate the effects of climate change, their narrow genetic diversity has always been a major constraint to their improvement for adaptability. However, existing genetic diversity still provides opportunities to exploit novel attributes for developing climate-resilient cultivars. The mining and exploitation of adaptive traits imparting tolerance/resistance to climate-smart pulses can be accelerated further by using cutting-edge approaches of biotechnology such as transgenics, genome editing, and epigenetics. This review discusses various classical and molecular approaches and strategies to exploit adaptive traits for breeding climate-smart pulses.

DNA Methylation Analysis of the Citrullus lanatus Response to Cucumber Green Mottle Mosaic Virus Infection by Whole-Genome Bisulfite Sequencing

DNA methylation is an important epigenetic mark associated with plant immunity, butlittle is known about its roles in viral infection of watermelon. We carried out whole-genomebisulfite sequencing of watermelon leaves at 0 h (ck), 48 h, and 25 days post-inoculation withCucumber green mottle mosaic virus (CGMMV). The number of differentially methylated regions(DMRs) increased during CGMMV infection and 2788 DMR-associated genes (DMGs) werescreened out among three libraries. Most DMRs and DMGs were obtained under the CHH context.These DMGs were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG)pathways of secondary biosynthesis and metabolism, plant-pathogen interactions, Toll-likereceptor signaling, and ABC transporters. Additionally, DMGs encoding PR1a, CaMs, calciumbindingprotein, RIN4, BAK1, WRKYs, RBOHs, STKs, and RLPs/RLKs were involved in thewatermelon-CGMMV interaction and signaling. The association between DNA methylation andgene expression was analyzed by RNA-seq and no clear relationship was detected. Moreover,downregulation of genes in the RdDM pathway suggested the reduced RdDM-directed CHHmethylation plays an important role in antiviral defense in watermelon. Our findings providegenome-wide DNA methylation profiles of watermelon and will aid in revealing the molecularmechanism in response to CGMMV infection at the methylation level.