Identification of Differentially Methylated miRNA Genes During Compatible and Incompatible Interactions Between Soybean and Soybean Cyst Nematode.

DNA methylation is a widespread epigenetic mark that affects gene expression and transposon mobility during plant development and stress responses. However, the role of DNA methylation in regulating the expression of microRNA (miRNA) genes remains largely unexplored. Here, we analyzed DNA methylation changes of miRNA genes using a pair of soybean (Glycine max) near-isogenic lines (NILs) differing in their response to soybean cyst nematode (SCN; Heterodera glycines). Differences in global DNA methylation levels over miRNA genes in response to SCN infection were observed between the isogenic lines. miRNA genes with significant changes in DNA methylation levels in the promoter and primary transcript-coding regions were detected in both lines. In the susceptible isogenic line (NIL-S), 82 differentially methylated miRNAs were identified in response to SCN infection whereas, in the resistant isogenic line (NIL-R), only 16 differentially methylated miRNAs were identified. Interestingly, gma-miR5032, gma-miR5043, gma-miR1520b, and gma-2107-ch16 showed opposite methylation patterns in the isogenic lines. In addition, the miRNA paralogs gma-miR5770a and gma-miR5770b showed hypermethylation and hypomethylation in NIL-S and NIL-R, respectively. Gene expression quantification of gma-miR5032, gma-miR5043, gma-miR1520b, and gma-miR5770a/b and their confirmed targets indicated a role of DNA methylation in regulating miRNA expression and, thus, their targets upon SCN infection. Furthermore, overexpression of these four miRNAs in NIL-S using transgenic hairy root system enhanced plant resistance to SCN to various degrees with a key role observed for miR5032. Together, our results provide new insights into the role of epigenetic mechanisms in controlling miRNA regulatory function during SCN-soybean interactions.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Identification of introduced and stably inherited DNA methylation variants in soybean associated with soybean cyst nematode parasitism.

DNA methylation is a widespread epigenetic mark that contributes to transcriptome reprogramming during plant-pathogen interactions. However, the distinct role of DNA methylation in establishing resistant and susceptible responses remains largely unexplored. Here, we developed and used a pair of near-isogenic lines (NILs) to characterize DNA methylome landscapes of soybean roots during the susceptible and resistant interactions with soybean cyst nematode (SCN; Heterodera glycines). We also compared the methylomes of the NILs and their parents to identify introduced and stably inherited methylation variants. The genomes of the NILs were substantially differentially methylated under uninfected conditions. This difference was associated with differential gene expression that may prime the NIL responses to SCN infection. In response to SCN infection, the susceptible line exhibited reduced global methylation levels in both protein-coding genes and transposable elements, whereas the resistant line showed the opposite response, increased global methylation levels. Heritable and novel nonparental differentially methylated regions overlapping with genes associated with soybean response to SCN infection were identified and validated using transgenic hairy root system. Our analyses indicate that DNA methylation patterns associated with the susceptible and resistant interactions are highly specific and that novel and stably inherited methylation variants are of biological significance.

Whole-genome re-sequencing reveals the impact of the interaction of copy number variants of the rhg1 and Rhg4 genes on broad-based resistance to soybean cyst nematode.

Soybean cyst nematode (SCN) is the most devastating plant-parasitic nematode. Most commercial soybean varieties with SCN resistance are derived from PI88788. Resistance derived from PI88788 is breaking down due to narrow genetic background and SCN population shift. PI88788 requires mainly the rhg1-b locus, while 'Peking' requires rhg1-a and Rhg4 for SCN resistance. In the present study, whole genome re-sequencing of 106 soybean lines was used to define the Rhg haplotypes and investigate their responses to the SCN HG-Types. The analysis showed a comprehensive profile of SNPs and copy number variations (CNV) at these loci. CNV of rhg1 (GmSNAP18) only contributed towards resistance in lines derived from PI88788 and 'Cloud'. At least 5.6 copies of the PI88788-type rhg1 were required to confer SCN resistance, regardless of the Rhg4 (GmSHMT08) haplotype. However, when the GmSNAP18 copies dropped below 5.6, a 'Peking'-type GmSHMT08 haplotype was required to ensure SCN resistance. This points to a novel mechanism of epistasis between GmSNAP18 and GmSHMT08 involving minimum requirements for copy number. The presence of more Rhg4 copies confers resistance to multiple SCN races. Moreover, transcript abundance of the GmSHMT08 in root tissue correlates with more copies of the Rhg4 locus, reinforcing SCN resistance. Finally, haplotype analysis of the GmSHMT08 and GmSNAP18 promoters inferred additional levels of the resistance mechanism. This is the first report revealing the genetic basis of broad-based resistance to SCN and providing new insight into epistasis, haplotype-compatibility, CNV, promoter variation and its impact on broad-based disease resistance in plants.

Cyst Nematode Parasitism Induces Dynamic Changes in the Root Epigenome.

A growing body of evidence indicates that epigenetic modifications can provide efficient, dynamic, and reversible cellular responses to a wide range of environmental stimuli. However, the significance of epigenetic modifications in plant-pathogen interactions remains largely unexplored. In this study, we provide a comprehensive analysis of epigenome changes during the compatible interaction between the beet cyst nematode Heterodera schachtii and Arabidopsis (Arabidopsis thaliana). Whole-genome bisulfite sequencing was conducted to assess the dynamic changes in the methylome of Arabidopsis roots in response to H. schachtii infection. H. schachtii induced widespread hypomethylation of protein-coding genes and transposable elements (TEs), preferentially those adjacent to protein-coding genes. The abundance of 24-nt siRNAs was associated with hypermethylation of TEs and gene promoters, with influence observed for methylation context and infection time points. mRNA sequencing revealed a significant enrichment for the differentially methylated genes among the differentially expressed genes, specifically those with functions corresponding to primary metabolic processes and responses to stimuli. The differentially methylated genes overlapped with more than one-fourth of the syncytium differentially expressed genes and are of functional significance. Together, our results provide intriguing insights into the potential regulatory role of differential DNA methylation in shaping the biological interplay between cyst nematodes and host plants.

Cooperative Regulatory Functions of miR858 and MYB83 during Cyst Nematode Parasitism.

MicroRNAs (miRNAs) recently have been established as key regulators of transcriptome reprogramming that define cell function and identity. Nevertheless, the molecular functions of the greatest number of miRNA genes remain to be determined. Here, we report cooperative regulatory functions of miR858 and its MYB83 transcription factor target gene in transcriptome reprogramming during Heterodera cyst nematode parasitism of Arabidopsis (Arabidopsis thaliana). Gene expression analyses and promoter-GUS fusion assays documented a role of miR858 in posttranscriptional regulation of MYB83 in the Heterodera schachtii-induced feeding sites, the syncytia. Constitutive overexpression of miR858 interfered with H. schachtii parasitism of Arabidopsis, leading to reduced susceptibility, while reduced miR858 abundance enhanced plant susceptibility. Similarly, MYB83 expression increases were conducive to nematode infection because overexpression of a noncleavable coding sequence of MYB83 significantly increased plant susceptibility, whereas a myb83 mutation rendered the plants less susceptible. In addition, RNA-seq analysis revealed that genes involved in hormone signaling pathways, defense response, glucosinolate biosynthesis, cell wall modification, sugar transport, and transcriptional control are the key etiological factors by which MYB83 facilitates nematode parasitism of Arabidopsis. Furthermore, we discovered that miR858-mediated silencing of MYB83 is tightly regulated through a feedback loop that might contribute to fine-tuning the expression of more than a thousand of MYB83-regulated genes in the H. schachtii-induced syncytium. Together, our results suggest a role of the miR858-MYB83 regulatory system in finely balancing gene expression patterns during H. schachtii parasitism of Arabidopsis to ensure optimal cellular function.

The cyst nematode effector protein 10A07 targets and recruits host posttranslational machinery to mediate its nuclear trafficking and to promote parasitism in Arabidopsis.

Plant-parasitic cyst nematodes synthesize and secrete effector proteins that are essential for parasitism. One such protein is the 10A07 effector from the sugar beet cyst nematode, Heterodera schachtii, which is exclusively expressed in the nematode dorsal gland cell during all nematode parasitic stages. Overexpression of H. schachtii 10A07 in Arabidopsis thaliana produced a hypersusceptible phenotype in response to H. schachtii infection along with developmental changes reminiscent of auxin effects. The 10A07 protein physically associates with a plant kinase and the IAA16 transcription factor in the cytoplasm and nucleus, respectively. The interacting plant kinase (IPK) phosphorylates 10A07 at Ser-144 and Ser-231 and mediates its trafficking from the cytoplasm to the nucleus. Translocation to the nucleus is phosphorylation dependent since substitution of Ser-144 and Ser-231 by alanine resulted in exclusive cytoplasmic accumulation of 10A07. IPK and IAA16 are highly upregulated in the nematode-induced syncytium (feeding cells), and deliberate manipulations of their expression significantly alter plant susceptibility to H. schachtii in an additive fashion. An inactive variant of IPK functioned antagonistically to the wild-type IPK and caused a dominant-negative phenotype of reduced plant susceptibility. Thus, exploitation of host processes to the advantage of the parasites is one mechanism by which cyst nematodes promote parasitism of host plants.

Arabidopsis miR827 mediates post-transcriptional gene silencing of its ubiquitin E3 ligase target gene in the syncytium of the cyst nematode Heterodera schachtii to enhance susceptibility.

MicroRNAs (miRNAs) are a major class of small non-coding RNAs with emerging functions in biotic and abiotic interactions. Here, we report on a new functional role of Arabidopsis miR827 and its NITROGEN LIMITATION ADAPTATION (NLA) target gene in mediating plant susceptibility to the beet cyst nematode Heterodera schachtii. Cyst nematodes are sedentary endoparasites that induce the formation of multinucleated feeding structures termed syncytia in the roots of host plants. Using promoter:GUS fusion assays we established that miR827 was activated in the initial feeding cells and this activation was maintained in the syncytium during all sedentary stages of nematode development. Meanwhile, the NLA target gene, which encodes an ubiquitin E3 ligase enzyme, was post-transcriptionally silenced in the syncytium to permanently suppress its activity during all nematode parasitic stages. Overexpression of miR827 in Arabidopsis resulted in hyper-susceptibility to H. schachtii. In contrast, inactivation of miR827 activity through target mimicry or by overexpression a miR827-resistant cDNA of NLA produced the opposite phenotype of reduced plant susceptibility to H. schachtii. Gene expression analysis of several pathogenesis-related genes together with Agrobacterium-mediated transient expression in Nicotiana benthamiana provided strong evidence that miR827-mediated downregulation of NLA to suppress basal defense pathways. In addition, using yeast two-hybrid screens we identified several candidates of NLA-interacting proteins that are involved in a wide range of biological processes and molecular functions, including three pathogenesis-related proteins. Taken together, we conclude that nematode-activated miR827 in the syncytium is necessary to suppress immune responses in order to establish infection and cause disease.

Altered Expression of a Chloroplast Protein Affects the Outcome of Virus and Nematode Infection.

The chloroplast-resident RNA helicase ISE2 (INCREASED SIZE EXCLUSION LIMIT2) can modulate the formation and distribution of plasmodesmata and intercellular trafficking. We have determined that ISE2 expression is induced by viral infection. Therefore, the responses of Nicotiana benthamiana plants with varying levels of ISE2 expression to infection by Tobacco mosaic virus and Turnip mosaic virus were examined. Surprisingly, increased or decreased ISE2 expression led to faster viral systemic spread and, in some cases, enhanced systemic necrosis. The contributions of RNA silencing and hormone-mediated immune responses to the increased viral susceptibility of these plants were assessed. In addition, Arabidopsis thaliana plants with increased ISE2 expression were found to be more susceptible to infection by the beet cyst nematode Heterodera schachtii. Our analyses provide intriguing insights into unexpected functional roles of a chloroplast protein in mediating plant-pathogen interactions. The possible roles of plasmodesmata in determining the outcomes of these interactions are also discussed.

The Role of Cytokinin During Infection of Arabidopsis thaliana by the Cyst Nematode Heterodera schachtii.

Plant-parasitic cyst nematodes induce the formation of hypermetabolic feeding sites, termed syncytia, as their sole source of nutrients. The formation of the syncytium is orchestrated by the nematode, in part, by modulation of phytohormone responses, including cytokinin. In response to infection by the nematode Heterodera schachtii, cytokinin signaling is transiently induced at the site of infection and in the developing syncytium. Arabidopsis lines with reduced cytokinin sensitivity show reduced susceptibility to nematode infection, indicating that cytokinin signaling is required for optimal nematode development. Furthermore, lines with increased cytokinin sensitivity also exhibit reduced nematode susceptibility. To ascertain why cytokinin hypersensitivity reduces nematode parasitism, we examined the transcriptomes in wild type and a cytokinin-hypersensitive type-A arr Arabidopsis mutant in response to H. schachtii infection. Genes involved in the response to biotic stress and defense response were elevated in the type-A arr mutant in the absence of nematodes and were hyperinduced following H. schachtii infection, which suggests that the Arabidopsis type-A arr mutants impede nematode development because they are primed to respond to pathogen infection. These results suggest that cytokinin signaling is required for optimal H. schachtii parasitism of Arabidopsis but that elevated cytokinin signaling triggers a heightened immune response to nematode infection.

The Methylome of Soybean Roots during the Compatible Interaction with the Soybean Cyst Nematode.

The soybean cyst nematode (SCN; Heterodera glycines) induces the formation of a multinucleated feeding site, or syncytium, whose etiology includes massive gene expression changes. Nevertheless, the genetic networks underlying gene expression control in the syncytium are poorly understood. DNA methylation is a critical epigenetic mark that plays a key role in regulating gene expression. To determine the extent to which DNA methylation is altered in soybean (Glycine max) roots during the susceptible interaction with SCN, we generated whole-genome cytosine methylation maps at single-nucleotide resolution. The methylome analysis revealed that SCN induces hypomethylation to a much higher extent than hypermethylation. We identified 2,465 differentially hypermethylated regions and 4,692 hypomethylated regions in the infected roots compared with the noninfected control. In addition, 703 and 1,346 unique genes were identified as overlapping with hyper- or hypomethylated regions, respectively. The differential methylation in genes apparently occurs independently of gene size and GC content but exhibits strong preference for recently duplicated paralogs. Furthermore, a set of 278 genes was identified as specifically syncytium differentially methylated genes. Of these, we found genes associated with epigenetic regulation, phytohormone signaling, cell wall architecture, signal transduction, and ubiquitination. This study provides, to our knowledge, new evidence that differential methylation is part of the regulatory mechanisms controlling gene expression changes in the nematode-induced syncytium.

Assessing the bioconfinement potential of a Nicotiana hybrid platform for use in plant molecular farming applications.

BACKGROUND: The introduction of pharmaceutical traits in tobacco for commercial production could benefit from the utilization of a transgene bioconfinement system. It has been observed that interspecific F1Nicotiana hybrids (Nicotiana tabacum × Nicotiana glauca) are sterile and thus proposed that hybrids could be suitable bioconfined hosts for biomanufacturing. We genetically tagged hybrids with green fluorescent protein (GFP), which was used as a visual marker to enable gene flow tracking and quantification for field and greenhouse studies. GFP was used as a useful proxy for pharmaceutical transgenes. RESULTS: Analysis of DNA content revealed significant genomic downsizing of the hybrid relative to that of N. tabacum. Hybrid pollen was capable of germination in vitro, albeit with a very low frequency and with significant differences between plants. In two field experiments, one each in Tennessee and Kentucky, we detected outcrossing at only one location (Tennessee) at 1.4%. Additionally, from 50 hybrid plants at each field site, formation of 84 and 16 seed was observed, respectively. Similar conclusions about hybrid fertility were drawn from greenhouse crosses. In terms of above-ground biomass, the hybrid yield was not significantly different than that of N. tabacum in the field. CONCLUSION: N. tabacum × N. glauca hybrids show potential to contribute to a bioconfinement- and biomanufacturing host system. Hybrids exhibit extremely low fertility with no difference of green biomass yields relative to N. tabacum. In addition, hybrids are morphologically distinguishable from tobacco allowing for identity preservation. This hybrid system for biomanufacturing would optimally be used where N. glauca is not present and in physical isolation of N. tabacum production to provide total bioconfinement.

Establishment and maintenance of DNA methylation in nematode feeding sites.

A growing body of evidence indicates that epigenetic mechanisms, particularly DNA methylation, play key regulatory roles in plant-nematode interactions. Nevertheless, the transcriptional activity of key genes mediating DNA methylation and active demethylation in the nematode feeding sites remains largely unknown. Here, we profiled the promoter activity of 12 genes involved in maintenance and de novo establishment of DNA methylation and active demethylation in the syncytia and galls induced respectively by the cyst nematode Heterodera schachtii and the root-knot nematode Meloidogyne incognita in Arabidopsis roots. The promoter activity assays revealed that expression of the CG-context methyltransferases is restricted to feeding site formation and development stages. Chromomethylase1 (CMT1), CMT2, and CMT3 and Domains Rearranged Methyltransferase2 (DRM2) and DRM3, which mediate non-CG methylation, showed similar and distinct expression patterns in the syncytia and galls at various time points. Notably, the promoters of various DNA demethylases were more active in galls as compared with the syncytia, particularly during the early stage of infection. Mutants impaired in CG or CHH methylation similarly enhanced plant susceptibility to H. schachtii and M. incognita, whereas mutants impaired in CHG methylation reduced plant susceptibility only to M. incognita. Interestingly, hypermethylated mutants defective in active DNA demethylation exhibited contrasting responses to infection by H. schachtii and M. incognita, a finding most likely associated with differential regulation of defense-related genes in these mutants upon nematode infection. Our results point to methylation-dependent mechanisms regulating plant responses to infection by cyst and root-knot nematodes.

Spatial and temporal expression patterns of auxin response transcription factors in the syncytium induced by the beet cyst nematode Heterodera schachtii in Arabidopsis.

Plant-parasitic cyst nematodes induce the formation of a multinucleated feeding site in the infected root, termed the syncytium. Recent studies point to key roles of the phytohormone auxin in the regulation of gene expression and establishment of the syncytium. Nevertheless, information about the spatiotemporal expression patterns of the transcription factors that mediate auxin transcriptional responses during syncytium formation is limited. Here, we provide a gene expression map of 22 auxin response factors (ARFs) during the initiation, formation and maintenance stages of the syncytium induced by the cyst nematode Heterodera schachtii in Arabidopsis. We observed distinct and overlapping expression patterns of ARFs throughout syncytium development phases. We identified a set of ARFs whose expression is predominantly located inside the developing syncytium, whereas others are expressed in the neighbouring cells, presumably to initiate specific transcriptional programmes required for their incorporation within the developing syncytium. Our analyses also point to a role of certain ARFs in determining the maximum size of the syncytium. In addition, several ARFs were found to be highly expressed in fully developed syncytia, suggesting a role in maintaining the functional phenotype of mature syncytia. The dynamic distribution and overlapping expression patterns of various ARFs seem to be essential characteristics of ARF activity during syncytium development.

Synchronization of developmental processes and defense signaling by growth regulating transcription factors.

Growth regulating factors (GRFs) are a conserved class of transcription factor in seed plants. GRFs are involved in various aspects of tissue differentiation and organ development. The implication of GRFs in biotic stress response has also been recently reported, suggesting a role of these transcription factors in coordinating the interaction between developmental processes and defense dynamics. However, the molecular mechanisms by which GRFs mediate the overlaps between defense signaling and developmental pathways are elusive. Here, we report large scale identification of putative target candidates of Arabidopsis GRF1 and GRF3 by comparing mRNA profiles of the grf1/grf2/grf3 triple mutant and those of the transgenic plants overexpressing miR396-resistant version of GRF1 or GRF3. We identified 1,098 and 600 genes as putative targets of GRF1 and GRF3, respectively. Functional classification of the potential target candidates revealed that GRF1 and GRF3 contribute to the regulation of various biological processes associated with defense response and disease resistance. GRF1 and GRF3 participate specifically in the regulation of defense-related transcription factors, cell-wall modifications, cytokinin biosynthesis and signaling, and secondary metabolites accumulation. GRF1 and GRF3 seem to fine-tune the crosstalk between miRNA signaling networks by regulating the expression of several miRNA target genes. In addition, our data suggest that GRF1 and GRF3 may function as negative regulators of gene expression through their association with other transcription factors. Collectively, our data provide new insights into how GRF1 and GRF3 might coordinate the interactions between defense signaling and plant growth and developmental pathways.

An orange fluorescent protein tagging system for real-time pollen tracking.

BACKGROUND: Monitoring gene flow could be important for future transgenic crops, such as those producing plant-made-pharmaceuticals (PMPs) in open field production. A Nicotiana hybrid (Nicotiana. tabacum × Nicotiana glauca) shows limited male fertility and could be used as a bioconfined PMP platform. Effective assessment of gene flow from these plants is augmented with methods that utilize fluorescent proteins for transgenic pollen identification. RESULTS: We report the generation of a pollen tagging system utilizing an orange fluorescent protein to monitor pollen flow and as a visual assessment of transgene zygosity of the parent plant. This system was created to generate a tagged Nicotiana hybrid that could be used for the incidence of gene flow. Nicotiana tabacum 'TN 90' and Nicotiana glauca were successfully transformed via Agrobacterium tumefaciens to express the orange fluorescent protein gene, tdTomato-ER, in pollen and a green fluorescent protein gene, mgfp5-er, was expressed in vegetative structures of the plant. Hybrids were created that utilized the fluorescent proteins as a research tool for monitoring pollen movement and gene flow. Manual greenhouse crosses were used to assess hybrid sexual compatibility with N. tabacum, resulting in seed formation from hybrid pollination in 2% of crosses, which yielded non-viable seed. Pollen transfer to the hybrid formed seed in 19% of crosses and 10 out of 12 viable progeny showed GFP expression. CONCLUSION: The orange fluorescent protein is visible when expressed in the pollen of N. glauca, N. tabacum, and the Nicotiana hybrid, although hybrid pollen did not appear as bright as the parent lines. The hybrid plants, which show limited ability to outcross, could provide bioconfinement with the benefit of detectable pollen using this system. Fluorescent protein-tagging could be a valuable tool for breeding and in vivo ecological monitoring.