Complex relationship between DNA methylation and gene expression due to Lr28 in wheat-leaf rust pathosystem

5-mC methylation study of sORFs in 3'UTR of transcription factor JUNGBRUNNEN 1-like during leaf rust pathogenesis in wheat

BACKGROUND

JUB1, a NAC domain containing hydrogen peroxide-induced transcription factor, plays a critical role in plant immunity. Little is known about how JUB1 responds to leaf rust disease in wheat. Recent discoveries in genomics have also unveiled a multitude of sORFs often assumed to be non-functional, to argue for the necessity of including them as potential regulatory players of translation. However, whether methylation on sORFs spanning the 3'UTR of regulatory genes like JUB1 modulate gene expression, remains unclear.

METHODS AND RESULTS

In this study, we identified the methylation states of two sORFs in 3'UTR of a homologous gene of JUB1 in wheat, TaJUB1-L, at cytosine residues in CpG, CHH and CHG sites at different time points of disease progression in two near-isogenic lines of wheat (HD2329), with and without Lr24 gene during leaf rust pathogenesis. Here, we report a significant demethylation of the CpG dinucleotides occurring in the sORFs of the 3'UTR in the resistant isolines after 24 h post-infection. Also, the up-regulated gene expression observed through RT-qPCR was directly proportional to the demethylation of the CpG sites in the sORFs.

CONCLUSIONS

Our findings indicate that TaJUB1-L might be a positive regulator in providing tolerance during leaf rust pathogenesis and cytosine methylation at 3'UTR might act as a switch for its expression control. These results enrich the potential benefit of conventional methylation assay techniques for unraveling the unexplored enigma in epigenetics during plant-pathogen interaction in a cost-effective and confidentially conclusive manner.

The Sink-Source Relationship in Cucumber (Cucumis sativus L.) Is Modulated by DNA Methylation

The optimization of the sink-source relationship is of great importance for crop yield regulation. Cucumber is a typical raffinose family oligosaccharide (RFO)-transporting crop. DNA methylation is a common epigenetic modification in plants, but its role in sink-source regulation has not been demonstrated in RFO-translocating species. Here, whole-genome bisulfite sequencing (WGBS-seq) was conducted to compare the nonfruiting-node leaves (NFNLs) and leaves of fruit setting (FNLs) at the 12th node by removing all female flowers in other nodes of the two treatments. We found considerable differentially methylated genes enriched in photosynthesis and carbohydrate metabolic processes. Comparative transcriptome analysis between FNLs and NFNLs indicated that many differentially expressed genes (DEGs) with differentially methylated regions were involved in auxin, ethylene and brassinolide metabolism; sucrose metabolism; and RFO synthesis pathways related to sink-source regulation. Moreover, DNA methylation levels of six sink-source-related genes in the pathways mentioned above decreased in leaves after 5-aza-dC-2'-deoxycytidine (5-Aza-dC, a DNA methyltransferase inhibitor) treatment on FNLs, and stachyose synthase (CsSTS) gene expression, enzyme activity and stachyose content in RFO synthesis pathway were upregulated, thereby increasing fruit length and dry weight. Taken together, our findings proposed an up-to-date inference for the potential role of DNA methylation in the sink-source relationship, which will provide important references for further exploring the molecular mechanism of DNA methylation in improving the yield of RFO transport plants.

Leaf rust responsive miRNA and their target genes in wheat

Small RNA sequencing (sRNA-seq) and degradome analysis were used for the identification of miRNAs and their target host genes in a pair of near-isogenic lines (NILs), which differed for the presence of leaf rust resistance gene Lr28. The study led to identification of (i) 506 known and 346 novel miRNAs; and (ii) 5054 target genes including 4557 in silico predicted and 497 degradome-based genes using 105 differentially expressed (DE) miRNAs. A subset of 128 targets (67 in silico + 61 degradome-based) was differentially expressed in RNA-seq data that was generated by us earlier using the same pair of NILs; among these 128 targets, 58 target genes exhibited an inverse relationship with the DE miRNAs (expression of miRNAs and activation/suppression of target genes). Eight miRNAs which belonged to the conserved miRNA families and were known to be induced in response to fungal diseases in plants included the following: miR156, miR158, miR159, miR168, miR169, miR172, miR319, miR396. The target genes belonged to the following classes of genes known to be involved in downstream disease resistance pathways; peroxidases, sugar transporters, auxin response signaling, oxidation-reduction, etc. It was also noticed that although a majority of miRNAs and target genes followed the above classical inverse relationship, there were also examples, where no such relationship was observed. Among the target genes, there were also 51 genes that were not only regulated by miRNAs, but were also differentially methylated at sequences including the following segments: promotors, introns, TSS, exons. The results of the present study suggest a complex interplay among miRNA genes, target genes, and various epigenetic controls, which regulate the expression of genes involved in downstream pathways for disease resistance.

An Update on Resistance Genes and Their Use in the Development of Leaf Rust Resistant Cultivars in Wheat

Wheat is one of the most important cereal crops in the world. The production and productivity of wheat is adversely affected by several diseases including leaf rust, which can cause yield losses, sometimes approaching >50%. In the present mini-review, we provide updated information on (i) all Lr genes including those derived from alien sources and 14 other novel resistance genes; (ii) a list of QTLs identified using interval mapping and MTAs identified using GWAS (particular those reported recently i.e., after 2018) and their association with known Lr genes; (iii) introgression/pyramiding of individual Lr genes in commercial/prominent cultivars from 18 different countries including India. Challenges and future perspectives of breeding for leaf rust resistance are also provided at the end of this mini-review. We believe that the information in this review will prove useful for wheat geneticists/breeders, not only in the development of leaf rust-resistant wheat cultivars, but also in the study of molecular mechanism of leaf rust resistance in wheat.

BS-Seq reveals major role of differential CHH methylation during leaf rust resistance in wheat (Triticum aestivum L.)

Epigenetic regulation of the activity of defense genes during onset of diseases or resistance against diseases in plants is an active area of research. In the present study, a pair of wheat NILs for leaf rust resistance gene Lr28 (R) in the background of an Indian cultivar HD2329 (S) was used for a study of DNA methylation mediated regulation of gene expression. Leaf samples were collected at 0 h before (S0 and R0) and 96 h after inoculation (S96 and R96). The DNA samples were subjected to BS-Seq and sequencing data were used for identification of differentially methylated/demethylated regions/genes (DMRs and DMGs). Following four pairs of comparisons were used for this purpose: S0 vs S96; S0 vs R0; R0 vs R96; S96 vs R96. Major role of CHH methylation relative to that of CG and CHG methylation was observed. Some important observations include the following: (i) abundance of CHH methylation among DMRs; (ii) predominance of DMRs in intergenic region, relative to other genomic regions (promoters, exons, introns, TSS and TTS); (iii) abundance of transposable elements (TEs) in DMRs with CHH context; (iv) demethylation mediated high expression of genes during susceptible reaction (S0 vs S96) and methylation mediated low expression of genes during resistant reaction (R0 vs R96 and S96 vs R96); (v) major genes under regulation encode proteins, which differ from those encoded by genes regulated during susceptible reaction and (vi) ~ 500 DMGs carried differential binding sites for H3K4/K27me3 marks suggesting joint involvement of DNA and H3 methylation. Thus, CHH methylation either alone or in combination with histone methylation plays a major role in regulating the expression of genes involved in wheat-leaf rust interaction.

Genome-Wide DNA Methylation and Its Effect on Gene Expression During Subclinical Mastitis in Water Buffalo

Subclinical mastitis (SCM) in buffalo is one of the most challenging paradoxes for the dairy sector with very significant milk production losses and poses an imminent danger to milch animal’s milk-producing ability. We present here the genome-wide methylation specific to SCM in water buffalo and its consequential effect on the gene expression landscape for the first time. Whole-genome DNA methylation profiles from peripheral blood lymphocytes and gene expression profiles from milk somatic cells of healthy and SCM cases were catalogued from the MeDIP-Seq and RNA-Seq data. The average methylation in healthy buffaloes was found to be higher than that in the SCM-infected buffaloes. DNA methylation was abundant in the intergenic region followed by the intronic region in both healthy control and SCM groups. A total of 3,950 differentially methylated regions (DMRs) were identified and annotated to 370 differentially methylated genes (DMGs), most of which were enriched in the promoter region. Several important pathways were activated due to hypomethylation and belonged to the Staphylococcus aureus infection, Th17 cell differentiation, and antigen processing and presentation pathways along with others of defense responses. DNA methylome was compared with transcriptome to understand the regulatory role of DNA methylation on gene expression specific to SCM in buffaloes. A total of 4,778 significant differentially expressed genes (DEGs) were extracted in response to SCM, out of which 67 DMGs were also found to be differentially expressed, suggesting that during SCM, DNA methylation could be one of the epigenetic regulatory mechanisms of gene expression. Genes like CSF2RB, LOC102408349, C3 and PZP like, and CPAMD8 were found to be downregulated in our study, which are known to be involved in the immune response to SCM. Association of DNA methylation with transposable elements, miRNAs, and lncRNAs was also studied. The present study reports a buffalo SCM web resource (BSCM2TDb) available at http://webtom.cabgrid.res.in/BSCM2TDb that catalogues all the mastitis-related information of the analyses results of this study in a single place. This will be of immense use to buffalo researchers to understand the host–pathogen interaction involving SCM, which is required in endeavors of mastitis control and management.

Epigenetics of wheat-rust interaction: an update

The outcome of different host-pathogen interactions is influenced by both genetic and epigenetic systems, which determine the response of plants to pathogens and vice versa. This review highlights key molecular mechanisms and conceptual advances involved in epigenetic research and the progress made in epigenetics of wheat-rust interactions. Epigenetics implies the heritable changes in the way of gene expression as a consequence of the modification of DNA bases, histone proteins, and/or non-coding-RNA biogenesis without disturbing the underlying nucleotide sequence. The changes occurring between DNA and its surrounding chromatin without altering its DNA sequence and leading to significant changes in the genome of any organism are called epigenetic changes. Epigenetics has already been used successfully to explain the mechanism of human pathogens and in the identification of pathogen-induced modifications within various host plants. Wheat rusts are one of the most vital fungal diseases throughout the major wheat-growing areas of the world. The epigenome in plant pathogens causing diseases such as wheat rusts is mysterious. The investigations of host and pathogen epigenetics in the wheat rusts system can offer a piece of suitable evidence for elucidation of the molecular basis of host-pathogen interaction. Besides, the information on the epigenetic regulation of the genes involved in resistance or pathogenicity will provide better insights into the complex resistance signaling pathways and could provide answers to certain key questions, such as whether epigenetic regulation of certain genes is imparting resistance to host in response of certain pathogen elicitors or not. In the last few years, there has been an upsurge in research on the host as well as pathogen epigenetics and its outcome in plant-pathogen interactions. This review summarizes the progress made in the areas related to the epigenetic control of host-pathogen interaction with particular emphasis on wheat rusts.

Crosstalk of DNA Methylation Triggered by Pathogen in Poplars With Different Resistances

DNA methylation plays crucial roles in responses to environmental stimuli. Modification of DNA methylation during development and abiotic stress responses has been confirmed in increasing numbers of plants, mainly annual plants. However, the epigenetic regulation mechanism underlying the immune response to pathogens remains largely unknown in plants, especially trees. To investigate whether DNA methylation is involved in the response to infection process or is related to the resistance differences among poplars, we performed comprehensive whole-genome bisulfite sequencing of the infected stem of the susceptible type Populus × euramerican ‘74/76’ and resistant type Populus tomentosa ‘henan’ upon Lonsdalea populi infection. The results revealed that DNA methylation changed dynamically in poplars during the infection process with a remarkable decrease seen in the DNA methylation ratio. Intriguingly, the resistant P. tomentosa ‘henan’ had a much lower basal DNA methylation ratio than the susceptible P. × euramerican ‘74/76’. Compared to mock-inoculation, both poplar types underwent post-inoculation CHH hypomethylation; however, significant decreases in mC and mCHH proportions were found in resistant poplar. In addition, most differentially CHH-hypomethylated regions were distributed in repeat and promoter regions. Based on comparison of DNA methylation modification with the expression profiles of genes, DNA methylation occurred in resistance genes, pathogenesis-related genes, and phytohormone genes in poplars during pathogen infection. Additionally, transcript levels of genes encoding methylation-related enzymes changed during pathogen infection. Interestingly, small-regulator miRNAs were subject to DNA methylation in poplars experiencing pathogen infection. This investigation highlights the critical role of DNA methylation in the poplar immune response to pathogen infection and provides new insights into epigenetic regulation in perennial plants in response to biotic stress.

Fungi Inhabiting the Wheat Endosphere

Wheat production is influenced by changing environmental conditions, including climatic conditions, which results in the changing composition of microorganisms interacting with this cereal. The group of these microorganisms includes not only endophytic fungi associated with the wheat endosphere, both pathogenic and symbiotic, but also those with yet unrecognized functions and consequences for wheat. This paper reviews the literature in the context of the general characteristics of endophytic fungi inhabiting the internal tissues of wheat. In addition, the importance of epigenetic regulation in wheat-fungus interactions is recognized and the current state of knowledge is demonstrated. The possibilities of using symbiotic endophytic fungi in modern agronomy and wheat cultivation are also proposed. The fact that the current understanding of fungal endophytes in wheat is based on a rather small set of experimental conditions, including wheat genotypes, plant organs, plant tissues, plant development stage, or environmental conditions, is recognized. In addition, most of the research to date has been based on culture-dependent methods that exclude biotrophic and slow-growing species and favor the detection of fast-growing fungi. Additionally, only a few reports of studies on the entire wheat microbiome using high-throughput sequencing techniques exist. Conducting comprehensive research on the mycobiome of the endosphere of wheat, mainly in the context of the possibility of using this knowledge to improve the methods of wheat management, mainly the productivity and health of this cereal, is needed.

Transcriptomic and Epigenomic Dynamics of Honey Bees in Response to Lethal Viral Infection

Honey bees (Apis mellifera L.) suffer from many brood pathogens, including viruses. Despite considerable research, the molecular responses and dynamics of honey bee pupae to viral pathogens remain poorly understood. Israeli Acute Paralysis Virus (IAPV) is emerging as a model virus since its association with severe colony losses. Using worker pupae, we studied the transcriptomic and methylomic consequences of IAPV infection over three distinct time points after inoculation. Contrasts of gene expression and 5 mC DNA methylation profiles between IAPV-infected and control individuals at these time points - corresponding to the pre-replicative (5 h), replicative (20 h), and terminal (48 h) phase of infection - indicate that profound immune responses and distinct manipulation of host molecular processes accompany the lethal progression of this virus. We identify the temporal dynamics of the transcriptomic response to with more genes differentially expressed in the replicative and terminal phases than in the pre-replicative phase. However, the number of differentially methylated regions decreased dramatically from the pre-replicative to the replicative and terminal phase. Several cellular pathways experienced hyper- and hypo-methylation in the pre-replicative phase and later dramatically increased in gene expression at the terminal phase, including the MAPK, Jak-STAT, Hippo, mTOR, TGF-beta signaling pathways, ubiquitin mediated proteolysis, and spliceosome. These affected biological functions suggest that adaptive host responses to combat the virus are mixed with viral manipulations of the host to increase its own reproduction, all of which are involved in anti-viral immune response, cell growth, and proliferation. Comparative genomic analyses with other studies of viral infections of honey bees and fruit flies indicated that similar immune pathways are shared. Our results further suggest that dynamic DNA methylation responds to viral infections quickly, regulating subsequent gene activities. Our study provides new insights of molecular mechanisms involved in epigenetic that can serve as foundation for the long-term goal to develop anti-viral strategies for honey bees, the most important commercial pollinator.

Identification and characterization of SET domain family genes in bread wheat (Triticum aestivum L.)

SET domain genes (SDGs) that are involved in histone methylation have been examined in many plant species, but have never been examined in bread wheat; the histone methylation caused due to SDGs is associated with regulation of gene expression at the transcription level. We identified a total of 166 bread wheat TaSDGs, which carry some interesting features including the occurrence of tandem/interspersed duplications, SSRs (simple sequence repeats), transposable elements, lncRNAs and targets for miRNAs along their lengths and transcription factor binding sites (TFBS) in the promoter regions. Only 130 TaSDGs encoded proteins with complete SET domain, the remaining 36 proteins had truncated SET domain. The TaSDG encoded proteins were classified into six classes (I–V and VII). In silico expression analysis indicated relatively higher expression (FPKM > 20) of eight of the 130 TaSDGs in different tissues, and downregulation of 30 TaSDGs under heat and drought at the seedling stage. qRT-PCR was also conducted to validate the expression of seven genes at the seedling stage in pairs of contrasting genotypes in response to abiotic stresses (water and heat) and biotic stress (leaf rust). These genes were generally downregulated in response to the three stresses examined.

Genome-wide analysis of H3K4me3 and H3K27me3 modifications due to Lr28 for leaf rust resistance in bread wheat (Triticum aestivum)

Present study revealed a complex relationship among histone H3 methylation (examined using H3K4/K27me3 marks), cytosine DNA methylation and differential gene expression during Lr28 mediated leaf rust resistance in wheat. During the present study, genome-wide histone modifications were examined in a pair of near isogenic lines (NILs) (with and without Lr28 in the background of cv. HD2329). The two histone marks used included H3K4me3 (an activation mark) and H3K27me3 (a repression mark). The results were compared with levels of expression (using RNA-seq) and DNA methylation (MeDIP) data obtained using the same pair of NILs. Some of the salient features of the present study include the following: (i) large scale differential binding sites (DBS) were available for only H3K4me3 in the susceptible cultivar, but for both H3K4me3 and H3K27me3 in its resistant NIL; (ii) DBSs for H3K27me3 mark were more abundant (> 80%) in intergenic regions, whereas DBSs for H3K4me3 were distributed in all genomic regions including exons, introns, intergenic, TTS (transcription termination sites) and promoters; (iii) fourteen (14) genes associated with DBSs showed co-localization for both the marks; (iv) only a small fraction (7% for H3K4me3 and 12% for H3K27me3) of genes associated with DBSs matched with the levels of gene expression inferred from RNA-seq data; (v) validation studies using qRT-PCR were conducted on 26 selected representative genes; results for only 11 genes could be validated. The proteins encoded by important genes involved in promoting infection included domains generally carried by R gene proteins such as Mlo like protein, protein kinases and purple acid phosphatase. Similarly, proteins encoded by genes involved in resistance included those carrying domains for lectin kinase, R gene, aspartyl protease, etc. Overall, the results suggest a very complex network of downstream genes that are expressed during compatible and incompatible interactions; some of the genes identified during the present study may be used in future validation studies involving RNAi/overexpression approaches.

Genome-wide identification and characterization of novel non-coding RNA-derived SSRs in wheat

Expression of eukaryotic genes is largely regulated by non-coding RNAs (ncRNA). Sequence variations in the regulatory RNAs may have critical biological consequences including transcriptional and post-transcriptional gene regulation. ncRNA-derived markers thus can be proved useful in molecular breeding, QTL mapping and association studies for trait dissection. In present study, we identified a total of 661 SSRs dwelling in pre-miRNA (15), small nuclear RNA (25) and lncRNA (621). Of these, 46 were validated and 100% amplification success was observed in selected wheat genotypes. A set of 36 ncRNA-SSRs markers was utilized for genetic variability assessment in forty-eight Indian wheat genotypes (which includes bread wheat, durum wheat and relatives). Number of alleles ranged from 1 to 4 with an average of two alleles per SSR locus. Mean PIC, observed heterozygosity and Shannon information index were found to be 0.258, 0.37 and 0.476 which suggests ncRNA-SSRs show higher polymorphism compared to genic SSRs but lower polymorphism compared to genomic SSRs. Thirty-six ncRNA-SSRs showed transferability ranging from 42.1% to 100%. Average genetic dissimilarity among wheat genotypes was found to be 0.29 based on Jaccard's dissimilarity. This is the first report of ncRNA-SSRs in wheat which will be useful for molecular breeding and genetic improvement of wheat.

The progress of leaf rust research in wheat

Leaf rust (also called brown rust) in wheat, caused by fungal pathogen Puccinia triticina Erikss. (Pt) is one of the major constraints in wheat production worldwide. Pt is widespread with diverse population structure and undergoes rapid evolution to produce new virulent races against resistant cultivars that are regularly developed to provide resistance against the prevailing races of the pathogen. Occasionally, the disease may also take the shape of an epidemic in some wheat-growing areas causing major economic losses. In the recent past, substantial progress has been made in characterizing the sources of leaf rust resistance including non-host resistance (NHR). Progress has also been made in elucidating the population biology of Pt and the mechanisms of wheat-Pt interaction. So far, ∼80 leaf rust resistance genes (Lr genes) have been identified and characterized; some of them have also been used for the development of resistant wheat cultivars. It has also been shown that a gene-for-gene relationship exists between individual wheat Lr genes and the corresponding Pt Avr genes so that no Lr gene can provide resistance unless the prevailing race of the pathogen carries the corresponding Avr gene. Several Lr genes have also been cloned and their products characterized, although no Avr gene corresponding a specific Lr gene has so far been identified. However, several candidate effectors for Pt have been identified and functionally characterized using genome-wide analyses, transcriptomics, RNA sequencing, bimolecular fluorescence complementation (BiFC), virus-induced gene silencing (VIGS), transient expression and other approaches. This review summarizes available information on different aspects of the pathogen Pt, genetics/genomics of leaf rust resistance in wheat including cloning and characterization of Lr genes and epigenetic regulation of disease resistance.