- Research
- Open access
- Published:
Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass
BMC Genomics volume 25, Article number: 223 (2024)
Abstract
Background
Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress.
Results
In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls.
Conclusion
Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.
Introduction
Switchgrass is cultivated across a vast geographical area, and it is economically important as it produces relatively higher cellulosic biomass with little water and nutrient input for bioenergy production. Growing biomass crops such as switchgrass on marginal lands is sustainable, and it helps in mitigating climate change, as they reduce CO2 emissions over four times more effectively than forest and grassland ecosystems. Since the US Department of Energy (DOE) identified switchgrass as a model biofuel crop, molecular and genomic breeding studies have gained momentum [1]. The ecotype- or genotype-specific genomic/epigenomic variation linked with stress tolerance and cellulosic biomass production can be exploited using biotechnological, and Next Generation Sequencing (NGS) approaches [2]. The polyploid (1C = ~ 1,130 Mb) genome of the lowland ecotype of switchgrass, AP13, that is used in this study is large and repetitive (~ 60%) [2]. Climate models predicted that the United States would be affected by moderate-to-severe drought and excessive heat in the next two decades due to climate change and global warming [3], which may decrease the productivity of agronomically important crops, including C3 and C4 grasses.
Plants are sessile and often are exposed to various abiotic stresses, such as heat and drought, either individually or in combination, adversely affecting plant growth, development, and productivity [4]. Drought is a meteorological condition where the plants receive inadequate rainfall and often show a high evapotranspiration rate and low water-holding capacity around the rhizosphere [4, 5]. Drought can influence plants' molecular, biochemical, and physiological processes, resulting in yield loss [4, 5]. Heat stress is often referred to as a condition when plants are exposed to excessive air temperatures for prolonged periods [6, 7]. Heat stress can influence physiological processes, such as photosynthesis, respiration, transpiration, membrane transport, homeostasis, and osmotic regulation, resulting in plant yield loss [6, 7]. The essential genes, pathways, regulatory networks, underlying molecular mechanisms, and interplay associated with plant drought and heat stress responses are not fully understood. Though drought and heat can independently induce stress-responsive genes and their associated pathways, it is quite possible that plants are exposed to multiple stresses [7], often a combination of drought and heat. When multiple stresses co-occur, plants respond simultaneously by triggering a few transcription factors (TFs), signaling molecules, and pathways. Complex crosstalk or interaction between key, intermediate, and regulatory biomolecules happens when plant stresses co-occur [8].
The concurrent effect of water deficit and high temperature on the growth and yield of crops under drought and heat stresses have been reported [9]. However, we are still trying to understand to what extent the interactive effects of heat and drought would affect the physiological responses of crops. Also, there is no clear evidence to support how crops would recover from heat, drought, and combined stresses. The consequence of drought and heat stress on carbon fluxes and storage was studied in switchgrass [10]. A previous report suggested that drought significantly affects various stages of switchgrass, and yield is dramatically reduced after three consecutive years of drought [11]. The authors further suggested that the early stages of switchgrass growth are more important to ensure biomass yield during subsequent years [11]. The physiological and gene expression variation on drought responses has been studied in switchgrass [12]. The role of small RNAs (microRNAs) has been reported in drought and heat stresses in switchgrass [13]. In another study, deep sequencing identified the regulatory role of miRNAs in drought and salinity stresses in switchgrass [14]. Transcriptome analyses identified over 16 heat-responsive genes in switchgrass [15]. Despite significant progress in classical genetics, molecular breeding, and modern biotechnological approaches, developing drought- and heat-tolerant crops are still challenging because of this constantly changing environment due to anthropogenic factors. Moreover, developing tolerant crops for combined stresses by gene stacking or pyramiding is complex and requires understanding genotype- or species-specific gene expression and regulation.
Chromatin immunoprecipitation-based sequencing (ChIP-Seq) is remarkable due to its high sensitivity and specificity in identifying protein-DNA interactions across the genome and provides high-resolution epigenomes. This technology has proven to be an efficient tool for generating genome-wide histone marks in Arabidopsis [16], rice [17], maize [18], brassica [19], and poplar [20]. The methylation and deacetylation of histone H3 lysine 9 (H3K9) and H3K27 are often associated with gene repression. While the acetylation and methylation of H3K4 and H3K36 are often associated with transcriptional activation in plant stress responses [21]. The positive correlation of H3K4me3 levels with the change in gene expression was first reported in Arabidopsis when subjected to dehydration stress [21]. Also, the genome-wide H3K4me3 modifications associated with genes responsive to drought stress have been reported in rice [22]. Using multiple NGS approaches, it is demonstrated that various epigenetic factors (H3K4me3, DNA methylation, and small RNAs) interact in a coordinated fashion to regulate the expression of heat stress-responsive genes in Arabidopsis [23]. The expression of heat shock proteins (HSPs) was altered in Arabidopsis when subjected to multiple abiotic stresses [24]. However, epigenomic modifications, including H3K4me3 in combined drought and heat stresses, have not been reported in switchgrass. As the transcriptionally active chromatin is marked by H3K4me3 around transcription start sites [25], genome-wide changes in gene expression in response to stress can be studied using ChIP-Seq. Therefore, this study evaluated genome-wide marking of H3K4me3 in switchgrass in response to drought, heat, and combined stresses using ChIP-Seq analysis.
Materials and methods
Experimental design and assessment of data quality
We adopted the cost-effective experimental design reported in our previous study on switchgrass [26]. This study used Alamo, AP13 genotype, a lowland ecotype. The selection of the AP13 from the switchgrass cultivar 'Alamo' was initially made at the University of Georgia. Then the genotype was relocated, and clonal copies were maintained at the greenhouse of Noble Research Institute, LLC, Ardmore, OK. The ramets of AP13 were then shifted into the 3-gallon nursery pots in the greenhouse, maintained under optimum conditions for 40 days, and then transported to growth chambers at the Noble Research Institute, LLC. Then the experiment was conducted with three biological replicates in a randomized complete block design starting five days after transfer to the growth chamber.
We randomly chose six pots for control (C), nine for DT, and nine for combined DT and HT (DT followed by HT) treatments during the transfer. The pots were separated into three treatments, randomly split into three groups, and maintained a minimum of three replicates per treatment. A diagram that shows how the growth chamber was partitioned for DT and HT treatments is depicted (Fig. 1). The control and DT treatments were maintained in a growth chamber. The DT inflicted with HT treatment was kept in another growth chamber, similar brand and model. The replicates for control and treatments were randomly distributed in the chamber (Fig. 1). In addition, for each condition, we maintained a positive/input control PC to compare it against the conditions and time points. The leaf samples were collected around the same time (~ 2:00 PM) and in the same position for all the treatments, with five leaves of the same age from major tillers (that had > 10 fully developed leaves) pooled per replicate. The plant tissues thus collected were flash-frozen using liquid nitrogen (LN2) and later stored at -80 °C for further processing.
Isolation and immunoprecipitation of chromatin
ChIP assay was carried out as described previously [27] and modified for Switchgrass. First, the finely powdered samples were settled in cold nuclear isolation buffer consisting of 1% formaldehyde with 20 µl of protease inhibitor (Thermo Scientific, Waltham, MA; Product #87786) at room temperature. Next, the sonicated-nuclear lysate chromatin samples were incubated with five µg of rabbit polyclonal antibody raised against a synthetic peptide corresponding to trimethyl-lysine 4 of histone H3 (Active Motif, Carlsbad, CA; 39915) along with 20 µl of Pierce protein A/G magnetic beads (Thermo Scientific, Waltham, MA) and incubated at 4 °C overnight. A negative control sample was processed similarly without adding the antibody. Additionally, input (positive control) samples were maintained with 50 µl of sonicated chromatin in the elution buffer. Then, the antibody-chromatin complex was washed, eluted, and de-cross-linked with 20 µl of 5 M NaCl and placed overnight at 65 °C. Then eluted ChIP DNA was recovered after Proteinase-K digestion, purified by phenol–chloroform extraction followed by ethanol precipitation using glycogen, and re-suspended in 20 µl of TE buffer. Finally, 100 ng of purified DNA generated a ChIP-Seq library for each treatment.
ChIP-Seq and data analysis
Data collection
In total, 63 libraries were generated from C, DT, DTHT, and PC samples and sequenced on Illumina HiSeqTM 2000 (Illumina, San Diego, CA) at the Delaware Biotechnology Institute (DBI, Newark, DE). Among these, 54 (3 conditions × 6 time-points × 3 replicates) libraries were from C, DT, and DTHT, and 9 (3 conditions × 1 positive/input control × 3 replicates) libraries were from PC. The ChIP-Seq data generated in this study is submitted to the GSE (Genomic Spatial Event database) section of the NCBI under bioproject number GSE196295 for the ChIP-Seq experiment.
Quality assessment
Read quality was evaluated using FastQC (v 0.11.2) for all samples, and quality trimming was performed with cutadapt (v 1.13) to remove bases with Phred33 score less than 30, and resulting reads of at least 50 bases were retained. The trimmed and filtered reads were mapped against the bowtie2-indexed Panicum virgatum reference genome (v4.1) using bowtie2 (v 2.3.1) [28] with default parameters other than the number of mismatches set to 1. The summary statistics, including the number of sequenced, trimmed, and % mapped reads for C, DT, and DTHT immunoprecipitated and PC samples, were given in Table 1.
Peak calling
Histone modification peaks were detected using SICER (Spatial Clustering for Identification of ChIP-Enriched Regions) (Version 1.1) [29] with the parameters as follows – redundancy threshold of 1, window size of 200, gap size of 200, the effective genome size of 0.7 (taken as a fraction of reference genome of Switchgrass), and FDR threshold controlling significance as 0.05 for each sample relative to their control sample. BAM files generated from Bowtie2 were converted to bed files using BEDTools (Version 2.21.0) [30] and used as input for SICER.
Identification of stress-responsive peaks
We prepared our experiment in such a way that additional heat treatment is not required when the heat-responsive peaks can be identified bioinformatically by comparing un-stressed (C), DT, DTHT (drought followed by heat), and input controls, especially when we included different collection time-points (0, 72, 96, 120, 144 or 168 h) in our study. Therefore, differential peaks were identified using different replications with a P-value cutoff of 0.01, window size 200, step size 100, gap size 0, and fragment size 300. To identify DT-responsive peaks, we required a region: 1) to be not identified as a differential peak between the drought treatment group and control group at 0 h; and 2) to be identified as a differential peak between the drought treatment group and control group in at least one of the time-points as follows (72, 96, 120, 144 or 168 h). To identify DTHT-responsive peaks, we required a region: 1) to be not identified as a differential peak between the group with a combination of drought and heat treatment and the control group at 0 h; and 2) to be identified as a differential peak between the group with a combination of drought and heat treatment and control group (C) in at least one of the following time-points (72, 96, 120, 144 or 168 h). To identify HT-responsive peaks, we required a region: 1) to be not identified as a differential peak between the group with a combination of drought and heat treatment and drought group at 0 h and 72 h; and 2) to be identified as a differential peak between the group with a combination of drought and heat treatment and control group in at least one of the time-points as follows (96, 120, 144 or 168 h).
Functional analysis of stress-responsive peaks
To investigate how DT and DTHT responsive peaks identified in this study correlated with the corresponding responsive genes in RNA-Seq data reported [25] (NCBI GEO accession # GSE174278), mapped files were screened for overlapping regions (with flanking 2 kb regions included) using intersect in BEDTools. Then GO enrichment analysis was conducted using agriGO for the overlapping responsive genes [31]. Overrepresented GO categories of P. virgatum genes were found using the PANTHER Classification System Version 13.1 [32] with default settings (including FDR < 0.05). PANTHER Pathways and the PANTHER classes of proteins were also identified.
Further, to determine functional gene enrichment and differentially acetylated gene interaction networks in Switchgrass, we used the genes identified in SICER for functional annotation. To analyze functionally enriched genes in the SICER analyses, we used the Database for Annotation, Visualization, and Integrated Discovery, DAVID [33]. As in the functional annotation charts, the fold enrichment score obtained from DAVID is defined as the percentage of genes in a given genome class. This study considered clusters with Benjamini factor < 0.05 and fold enrichment > 1.3 as significant and non-informative classifications were eliminated [34]. Pathway analysis was performed using MapMan with customized input files generated explicitly for Switchgrass using the Mercator tool to study the gene functions associated with responsive peaks.
The Mercator tool can batch classify sequences (gene/protein) into functional plant categories and make a draft metabolic network that can be used directly in MapMan software [35]. The experiment file has three possible values: 0, 1, and -1. i) "0" means a given gene was not identified as responsive in a particular condition. ii) "1" means the gene was identified as responsive and showed upregulated in at least one of the comparisons in a specific condition. iii) "-1" means the gene was identified as responsive and showed down-regulated in the comparisons. To make the color scale suitable for visualization, "scale" was adjusted to 1 in MapMan.
Visualization of data
Plots were generated using Circos [36], color keys were chosen from ColorBrewer (http://colorbrewer2.org) [37], and figure legends were included in the Circos plot by Inkscape (http://inkscape.org). Then ChIP-Seq data were visualized by Integrative Genome Viewer (IGV; http://software.broadinstitute.org/software/igv/).
RT-PCR and real-time quantitative RT-PCR (qPCR) validation
The RT-PCR and real-time quantitative RT-PCR (qPCR) validations were performed to assess the amplification of DNA/cDNA by using MyCycler thermocycler (Bio-Rad Laboratories, Hercules, CA) and ABI 7500 real-time PCR (Applied Biosystems, Foster City, CA), respectively. Our ChIP-Seq analysis selected six significant peaks associated with DT- and DTHT-responsive genes also reported in Switchgrass [38]. The details of genes and respective primers are provided in Supplemental Table 1. The selected genes' primers were designed using a TaqMan primer design tool for real-time PCR (GenScript USA Inc., Piscataway, NJ). The cDNA (10 ng) extracted from switchgrass leaves was used as a template for 25 μl qPCR reactions in triplicates, in which 10 μM of primer pairs (Forward and Reverse) and 12.5 μl of SYBR Green (Germantown, MD) PCR Master Mix. PCR conditions for qPCR were as follows: 95 °C for 10 min, 40 iterations of 95 °C for 15 s, and 65 °C for 60 s. To normalize the results, cons7 was used as a constitutive control of expression for all tissue samples. The efficiency of primers was tested and analyzed by using the previously reported 2-ΔΔCT method [39], where ΔΔCT = (CT of gene—CT of cons7) tissue to be observed—(CT of Genex—CT of cons7) leaf tissue. The normalized CT values (ΔΔCT) from qPCR analysis were collected and analyzed using Minitab 17, and the expression results were presented as mean ± SE. One-way ANOVA was done on qPCR experiments for multiple comparisons between the mean of samples.
Results
Analysis of ChIP-Seq peaks
ChIP-Seq analysis of 63 frozen leaf samples generated 5,361 million reads (300 bp fragments) at 85 million reads per sequenced sample (Table 1). For each condition, we maintained a positive control, namely 0PC, DTPC, DTHTPC. The total number of reads in C, DT, DTHT, and PC were approximately 1.6, 1.4, 1.6, and 0.7 billion reads, respectively. Deconseq [40] analysis revealed that the data collected was of high quality and devoid of contamination. Mapping of quality ChIP-Seq reads to reference genome (v) revealed that about 90% of the reads from each library were aligned and used for downstream analysis (Table 1).
Histone modification peaks were detected using the SICER program [29] with a false discovery cutoff of 0.05. In total, 1,374,515 H3K4me3 peaks were identified from all samples (Supplemental Table 2). Among these, DT-H3K4me3 and DTHT-H3K4me3 peaks were 754,409 and 620,106, respectively. On an average, each DT and DTHT sample had 41,900 and 39,481 peaks, respectively (Supplemental Table 3). The median width of a peak in the samples ranged from 1,200 bp to 1,600 bp, and the mean ranged from 1,356 bp to 2,119 bp. On average, the peaks identified on the mapped genome were about 68 Mb per sample. The sites of the enriched peaks on the switchgrass genome are shown (Supplemental Table 2). H3K4me3, an active histone modification that is associated with a transcriptionally active state of chromatin, exhibits conventional and non-conventional patterns, which are generally conserved in most eukaryotes. The canonical H3K4me3 pattern is narrow and associated with actively transcribed promoters and CpG islands. In contrast, the non-canonical H3K4me3 pattern is broad and associated with transcribed gene bodies. In our study, the distribution of enriched H3K4me3 peaks was higher downstream of the transcription start sites (TSS) of transcriptionally active genes (Figs. 2 and 3). While a few peaks (< 50%) were found in genic regions (gene body and flanking 2 kb areas), including CDS, downstream (2 kb from the downstream of the stop codon), 5’-UTR, and upstream (from 2 kb upstream of TSS) elements (Supplemental Table 2; Supplemental Fig. 1).
Identification of stress-responsive peaks
H3K4me3 is a near-universal chromatin modification at the TSS of active genes in eukaryotes and its levels reflect the rate of transcription. Our results indicate that many transcribed genes are not annotated. Reads per million (RPM) of the differential peaks were presented in Supplemental Table 4, and the distribution of stress-responsive peaks based on genomic features was extracted and presented (Fig. 3). The peaks in genic regions were enriched in 5’-UTR and CDS (Fig. 3). The number of DT, DTHT, and HT stress-responsive peaks identified were 7,342; 6,510; and 8,536, respectively (Table 2). We further studied how the DT- and DTHT-responsive peaks correlate with the corresponding candidate genes identified in previously reported RNA-Seq analyses [25]. It is interesting to report that 155 DT-responsive peaks overlapped with 118 DT-responsive genes (Supplemental Table 5). Similarly, 121 DTHT-responsive peaks overlapped with 110 DTHT-responsive genes (Supplemental Table 5). The overlap of the epigenomic peaks and genes could be seen as an interplay between the master regulators of the drought and heat stress-responsive genes. The nearest gene distribution was determined by BEDTools (v2.21.0) [29]; over 6000 H3K4me3 DT and DTHT peaks were spanned between 0 and + 2000 bases from TSS (Supplemental Table 5). Among the peaks identified, ~ 25% of peaks were enriched above 2000 bp from the TSS. The distribution of enriched H3K4me3 levels 2 Kb upstream and downstream of TSS decreased in both treatments gradually. We provided the GO term enrichment results of the genes with the stress-responsive peaks in Supplemental Table 6. We performed a Circos plot analysis to study the binding patterns of H3K4me3 during heat and drought stress in switchgrass. In the circos plot, we presented peak densities of DT and DTHT treatments in the respective chromosomes by comparing them with the genome-wide peak density (Fig. 4). Chr01K has the lowest peak density among them, and Chr07N has the highest peak density.
Functional (GO) analysis of stress-responsive peaks
Functional (GO) analysis was performed using agriGO for the overlapping responsive genes (http://bioinfo.cau.edu.cn/agriGO/analysis.php). GO terms for the genes were extracted from “Pvirgatum_450_v4.1.annotation_info.txt”. We identified all three GO categories with significantly over-represented members: Biological Process (45), Molecular Function (28), and Cellular Component (8) when DTHT was compared against DT ChIP-Seq. Similarly, the GO categories were identified between DT and C, and DTHT and C ChIP-Seq samples of switchgrass (Supplemental Table 6). The functional (GO) categories significantly (FDR < 0.05 and p-values < 0.01) enriched in biological processes were: metabolic processes (cellular, primary, nitrogen compound, and macromolecular), cellular processes, and response to stimuli (Fig. 5). Four members of molecular function category that significantly enriched were: transmembrane transporter, transferase, catalytic, and binding activities (Supplemental Table 6). Mapman analysis was performed to study the relationship between heat tolerance and fatty acid metabolism. Our results showed genes involved in fatty acid synthesis were highly marked during the HT and DT stresses (Fig. 6). We performed pathway analysis of the enriched differentially marked genes using the PANTHER database. Most of the differentially marked genes in cellular component category were classified as cell parts (41%), organelle (30%), membrane activity (14%), and macromolecular complex (13%) (Fig. 7).
Transcription factors (TFs) associated with stress-responsive peaks
This study identified various TFs that were involved in DT, DTHT, and HT stresses. The list of TFs involved in DT treatment was given in Tables 3, 4 and 5. There were five TFs viz., ABA inducible bHLH type (AtAIB), MADS-box protein Soc 1, growth-regulating factor, chaperone-protein dnaj-related, and cleavage and polyadenylation specificity factor 100. The TFs that responded during HT stress were KUA1, MADS-box protein Soc 1, CW-type Zinc finger, and NAC-domain containing protein 86. The TFs that were enriched during DTHT stress were auxin-response factor 16 (ARF 16), ABA-inducible bHLH types A1B, BR1-EMS suppressor (BES), chromatin remodeling factor CHD3, LEAFY (LFY), and myeloblastosis (MYB) (Tables 3, 4 and 5). This study also revealed various genes that responded to the stimulus. The commonly identified gene in DT and DTHT was L-type lectin domain-containing receptor-like protein kinase IV.1 (LFCRK41) (Tables 6, 7 and 8). The other stimulus-responsive genes identified belonged to chloroplastic exonuclease V, cysteine-rich receptor-like kinase 8 (CRK8), and 17.4 kDa class III heat-shock proteins, HSPs. We also identified various transporter genes that were marked in DT and DTHT treatments. The important transporter genes differentially marked in DT and DTHT stresses were ABC transporter C family member 3 and calcium-transporting ATPase-8 (Tables 6, 7 and 8).
Validation of ChIP-Seq results
To validate our results from ChIP-Seq analysis, ChIP-qPCR was conducted using triplicates for control and stressed leaf samples. Our comprehensive ChIP-Seq analysis revealed 162, 185, and 125 genes correlated with peaks between DT vs. C, DTHT vs. C, and DTHT vs. DT, respectively, for H3K4me3 modification. Among these, six essential genes (Supplemental Table 5) involved in abiotic stress tolerance have been selected to validate the expression level by qPCR. The relative expressions of all six genes indicated significant (P < 0.05) marks under drought and heat stresses compared to combined stress treatment. Further, for all the treatments, the expression levels were higher for transmembrane protein 14C, followed by fatty acid desaturase A, and chromatin-remodeling factor CHD3 (Fig. 8). The expression of the above-validated genes by ChIP-qPCR correlated with our ChIP-Seq analysis. Furthermore, the distribution of histone H3K4me3 peaks in the Auxin responsive factor 16 gene on chromosome 09 K was visualized using IGV genome browser (Fig. 9).
Discussion
Analysis of ChIP-Seq peaks
Histone H3 lysine 4 trimethylation (H3K4me3) has been used as an active mark to identify histone modifications associated with drought-responsive genes in Arabidopsis [21, 41]. A substantial increase in H3K4me3 activity in response to heat stress has been reported in Arabidopsis [21]. Similarly, we utilized H3K4me3 mark to detect transcriptionally active genes that respond to drought (DT), heat (HT), and combined drought and heat (DTHT) stresses. This is the first ChIP-Seq report on DT, HT, and DTHT stresses in switchgrass.
This study presents a modified ChIP-Seq pipeline for preprocessing and analyzing H3K4me3 data in AP13, a lowland switchgrass ecotype. This helps us to investigate the role of unfathomed epigenetic mechanisms in switchgrass that are likely linked with habitat preference. Several lines of evidence evaluated the efficacy of our pipeline: 1) ~ 89.21 – 94.66% of reads were aligned with the reference genome, which is comparable to the results obtained from other contemporary studies in plants, for example, Arabidopsis (95%) [42] brassica (93.0–98.0%) [43], sorghum (97.0%) [44], rice (90.0%–99.0%) [45], and maize (85.0%–92.0%) [46]. 2) The number of uniquely mapped H3K4me3 peaks (7,342, 6,510, and 8,536 peaks were responded to DT, DTHT, and HT, respectively) identified in our study were smaller when compared with those in Arabidopsis [42, 47], brassica [43, 48], rice [45], sorghum [44, 49], and maize [46, 49], yet the genome size of switchgrass (1129.9 Mb) is larger than Arabidopsis (~ 140 Mb), brassica (~ 920 Mb), rice (~ 430 Mb), and sorghum (~ 730 Mb), which may be due to lack of high-quality annotations in switchgrass. 3) We found diverse H3K4me3 distribution patterns in switchgrass, predominantly enriched in the genic regions and TSS, consistent with previous studies in other organisms [42,43,44,45,46,47,48,49,50,51,52,53]. 4) Our ChIP-Seq analysis was compared with RNA-seq data to identify the relationship between actively transcribed genes and genes altered by H3K4me3. However, the mean expression level of unmarked genes is lower than H3K4me3 marked genes. 5) Gene annotation and ontology analysis of H3K4me3 showed the enrichment of genes in cellular and metabolic processes and signaling. In addition, the H3K4me3 distribution patterns and functional annotations found in this study were comparable with previous reports in other plant species, indicating the high dependability of our ChIP-Seq data.
Identification and Functional analyses of stress-responsive peaks
After alignment, we annotated the peaks for genic and non-genic regions. We identified that the peaks in the genic regions were mainly enriched in 5'-UTR and CDS regions in switchgrass. Our study is in accordance with a recent report that showed the enrichment of peaks in the SUMO-associated genes, primarily around TSS, promoter-proximal region, and putative 5' UTR regions in rice [54]. In another report, enriched H3K4me3 histone modification correlated with the transcription of 3' upstream and TSS regions [41]. Most of the H3K4me3 modifications have been identified at the 5'-ends under drought treatment in Arabidopsis [41]. In addition, an increase in H3K4me3 trimethylation in the promoter region has been reported in Arabidopsis [21].
Our functional analyses showed the enrichment of H3K4me3 around the TSS and 5’end of the transcriptionally active genes, which is consistent with the distribution of histone modifications in other plants (Arabidopsis [55], rice [44, 51], sorghum [49], apple, [56], soybean [57], eucalyptus [58]), implying a conserved regulatory role of H3K4me3. In addition, we assigned the functions associated with the genes marked by H3K4me3. The H3K4me3 marked peaks that overlapped with annotated genes mostly belonged to molecular regulation, physiological processes, and energy metabolism, indicating the involvement of H3K4me3 in controlling housekeeping genes [59] and stress-induced signaling pathways [60, 61]. The H3K4me3 peaks corresponded with the annotated genes and were more likely to be activated, including the highly expressed genes in response to abiotic stresses, for example, 17.4 kDa class III HSP, chloroplast exonuclease V, ABA Inducible bHLH-Type (AtAIB), MADS-box protein Soc 1, NAC domain-containing protein 86, Auxin response factor 16 (ARF16), and Chromatin remodeling factor (CHD3). Most of these genes linked with stress-related processes altered by H3K4me3 have also been reported in other plant species [61,62,63,64]. These findings suggest that H3K4me3 probably modulates active transcription in switchgrass but with a distinct pattern of genomic distributions and regulatory roles. However, in contrast to the previous findings [65,66,67], H3K4me3 marks were widely distributed in switchgrass. We identified 22% (1950 out of 8536) of all the peaks were enriched around TSS and 34% (3015 out of 8536) in genic regions in switchgrass. This is compatible with recent findings in Arabidopsis that the H3K4me3 mark was significant in TSS and genic regions [60, 65, 68, 69], indicating that domains might be conserved in organisms [45, 70, 71]. Species-specific H3K4me3 distribution patterns were reported in plants, while some genomes were highly methylated, and a few were sparsely methylated [72,73,74]. The epigenetic mechanisms often act in a coordinated fashion to influence plant functions. Our results indicate that most H3K4me3 distribution and regulation were conserved, while a few are phylogenetically unique.
Our results indicated that most of the H3K4me3 peaks responded to DT, DTHT, and HT were enriched within 2000 bp from the TSS of protein-coding genes. A previous study reported over 40% of the HSFA1b-bound genes spanned around 250 bp from the TSS under heat stress in Arabidopsis [75]. In another study, the enrichment of OSABF1V has been found close to 200 bp of the TSS site in response to drought stress in rice [76]. The H3K4me3 mark is prevalent in actively transcribing genes in Arabidopsis [77]. Our study identified 155, 121, and 175 genes marked by H3K4me3 modification in response to DT, DTHT, and HT, respectively. A previous rice study showed a direct correlation between the enrichment of the H3K4me3 and 423 drought stress-responsive genes [22].
RNA-dependent RNA polymerase (RDR) family protein is involved in various abiotic stresses. Our study identified the multiple copies of the RDR family protein on chromosome 07. Among them, RDR6 is the most responsive to various stresses. The ChIP-Seq analysis of 15-day-old seedlings identified H3K4me3 enrichment in four drought-responsive genes including RDR6 under the dehydration stress in plants [78]. The positive correlation between H3K4me3 abundance and RDR6 in dehydration stress has been reported in 4-week-old rosette Arabidopsis leaves [79]. Histone-mediated regulation of the RDR6 gene under high temperatures has been observed in Arabidopsis [80]. A recent study reported an essential role for sRDR6 in PEG-induced drought stress in Saccharum spontaneum [81]. In another study, HvRDR2 has been significantly induced in response to heat stress in barley [82].
Gene Ontology analysis revealed various GO categories involved in biological, molecular, and cellular processes of switchgrass. Most biological processes identified here belonged to transmembrane transporter, transferase, catalytic, and binding activities. Similar to a previous drought stress study in rice [83], we found a significant number of drought-responsive genes that belong to transmembrane/transporter activity, plant hormones, and carbohydrate metabolism categories in switchgrass. We found ten genes commonly expressed or marked between two datasets by comparing ChIP-Seq and RNA-Seq analyses of DT, DTHT, and HT stresses. Among these, the essential genes that triggered in response to abiotic and biotic stresses with known functions were: 1) Aldolase-type TIM barrel family protein, also known as pyridoxal 5'-phosphate synthase pdxS subunit (pdxS and pdx1), which is involved in protein binding and cofactor biosynthesis that plays a role in protecting cellular membranes from lipid peroxidation [84]; 2) Condensin, subunit H, is an essential gene required for chromosome stability, condensation, and segregation [85]; 3) NB-ARC domain-containing or LRR containing proteins that are known to be involved in resistance to plant stresses [86]; 4) Pectin acetylesterase family protein or Notum-related proteins that play a role in structural stability, binding, and catalytic activity [87]; 5) Putative Mg-protoporphyrin IX chelatase subunit CHLH that is involved in chlorophyll biosynthesis [88]; 6) Receptor lectin kinases that are involved in plant development and stresses [89]; and 7) rRNA N-glycosylase/N-glycosidase enzymes that are induced during abiotic stresses [90] (Table 9).
Transcription factors (TFs) associated with stress-responsive peaks
TFs regulate the expression of gene/s or gene sets or networks that often respond to various abiotic and biotic stimuli. For example, MADS-Box proteins are crucial in triggering gene regulatory webs that impart drought stress tolerance in plants [91]. This study identified the MADS-Box protein Soc 1, a highly conserved DNA binding domain-containing transcriptional factor in both drought and combined (drought and heat) stresses. Often MADS-domain proteins interact with DNA binding sites at the CArG box in the promoter region to initiate transcription by recruiting co-factors and chromatin remodeling proteins [92].
The activation mark H3K4me3 recruits TFs, co-factors, and CRFs to modulate gene expression. Another TF, bHLH (basic helix-loop-helix), was found to be a positive regulator for root development and drought tolerance in maize [93], peanut [94], and Arabidopsis [95]. We found two bHLH TFs, ABA inducible bHLH-type (AtAIB) and bHLH DNA binding family proteins that were differentially marked by H3K4me3 under drought and combined stresses. In a previous report, AtAIB positively regulated the phytohormone ABA, which activates a significant set of genes involved in signaling when subjected to drought stress in Arabidopsis [96, 97]. We also identified the MYB, NAC, LFY, and zinc finger families of TFs that play a critical role in imparting abiotic stress tolerance. These families of TFs are not direct components of the ABA-induced signaling pathway. However, the direct correlation of ABA sensitivity with the MYB activity resulted in ABA-induced stomatal closure and increased drought tolerance in plants [98].
In rice, OsbZIP23 is identified as a key regulator in ABA signaling and drought resistance [99]. In a recent study, transcription factor OsbZIP23 has been reported to regulate the expression of dehydrin genes directly and positively and H3K4me3 modification levels under drought stress. The use of dehydrin gene as a drought-sensitive marker in investigating gene regulation at both the transcriptional and epigenetic levels has been reported [100]. Another study in rice reported that wide grain 7, a cysteine-tryptophan (CW) domain-containing transcriptional activator may increase the expression of grain size gene (OsMADS1), by binding directly to its promoter, thus increasing the histone H3K4me3 enrichment, and ultimately resulting in increased grain width [101].
The critical role of DNAJ proteins in enhancing drought tolerance and resistance to Pseudomonas solanacearum in transgenic tobacco has been reported [102]. In Arabidopsis, the interaction of chaperone DNAJ-like zinc finger with HSP70 in protein homeostasis under heat stress has been reported [103]. In a recent study in rice, the elevated expression levels of DNAJ-proteins at 3 h and 6 h time-points under drought stress has been reported [104]. The NtDnaJ1 gene overexpression in A. thaliana exhibited higher DNAJ protein ortholog levels and tolerance to drought stress [105]. The abnormal recruitment of DNAJ1-mediated transcription activation and its role in crosstalk between other epigenetic factors in plants, yeast, and mammals has been reported [106]. The SET- and RING-associated (SRA) domain of SUVH1/SUVH3 that binds all three contexts (CG, CHG, and CHH) of methylated DNA with a preference towards the CHH context was reported in Arabidopsis [106]. Further, these methylated regions recruit DNAJ domain-containing proteins to promote the expression of adjacent genes [106]. Similarly, in rice, when subjected to salt stress, the SRA domain of SUVH7 binds to CHH methylated regions (MITE transposable elements) within the promoter of a Na + /K + transporter (HKT15) to repress the gene expression [107]. While SUVH7 recruits a BAG4 chaperone or an MYB106 TF to promote the activation of genes in the vicinity of HKT15 transporters. The SRA domain of SUVH homologs has a role in evading transposon- or methylation-mediated gene silencing, thus activating the neighboring genes by recruiting chaperone regulators or TFs. The DNAJ-related chaperones enriched in drought, heat, and combined stresses in this study may have similar roles in promoting gene expression.
The auxin response factors (ARFs) are primarily essential for plant growth and development. Also, several reports corroborated the role of ARFs in tolerance to diverse abiotic stresses [108,109,110]. The extent of changes in auxin distribution and signaling has been investigated in drought stress in tomatoes [108]. They identified many ARFs that potentially mediate the auxin-responsive genes in drought and heat stresses. Similarly, we identified multiple copies of auxin response factor 16 under drought, heat, and combined stresses in switchgrass.
Stress-induced genes associated with stress-responsive peaks
We identified a plasma membrane protein, receptor-like protein kinase (RLK), that triggers Reactive Oxygen Species (ROS) signaling in abiotic stresses. The cysteine-rich protein kinase (CRK) receptor identified here may act as a sensor to recognize environmental stimuli, and its role in ROS signaling has been suggested in overcoming stress responses. In an earlier study, the overexpression of subfamily CRK45 showed increased tolerance against drought in Arabidopsis 110]. Our study revealed H3K4me3-mediated activation of CRK8 under drought stress, as identified in Arabidopsis [109]. Also, histone-mediated activation was seen in proteins that responded to stimuli, such as LECRK41 protein kinase, chloroplastic exonuclease V, and 17.4 kDa class III HSP.
Transporters associated with stress-responsive peaks
ATP-binding cassette (ABC) transporter plays a vital role in transporting abscisic acid (ABA) to guard cells to prevent water loss during drought stress [111]. In Arabidopsis, the ABCG40 transporter is localized in guard cells, and it is critical for the proper plant response to ABA. Furthermore, the functional role of the ABCG40 transporter is implicated in importing ABA to stomatal cells and drought stress [112]. In Medicago, the overexpression of the Chr6g0452171 gene that encodes ABC transporter B family protein enriched the histone mark under phosphorus deficiency (PD) and their possible role in PD-induced root system architecture (RSA) remodeling has been reported [113].
Also, the ABC transporter B family proteins are involved in auxin transport and iron homeostasis under PD in Arabidopsis [114]. Similarly, H3K4me3 marked over ten ABC transporter C family member genes in our study in response to drought, heat, and combined stresses, suggesting a functional role in RSA remodeling by modulating ion homeostasis. In poplar, the possibility of a change in phosphorylation is one of the most critical regulatory ways to modulate stress responses has been reported [115]. Also, they identified amino acid transporter (AVT6A-like) as the most interconnected gene with protein phosphatase 4 mediated phosphorylation when subjected to drought stress [115]. Similarly, we found enriched AVT6A-like transporter levels during drought, heat, and combined stresses in switchgrass.
Conclusion and future perspectives
This is the first genome-wide study investigating heat, drought, and combined stress in switchgrass. We observed an enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes. Circos plot indicates the peak density of DT and DTHT treatments across the respective chromosomes relative to a genome-wide peak density. GO categories showed significant enrichment primarily in four functional categories: transmembrane transporter, transferase, catalytic, and binding activities. Mapman analysis revealed that fatty acid synthesis genes were highly marked in response to DT and HT stresses. The PANTHER analysis revealed that many genes can be classified into cellular components, organelles, membrane activity, and macromolecular complexes. During combined drought and heat stresses, TFs auxin-response factor 16 (ARF 16), ABA-inducible bHLH type A1B, BR1-EMS suppressor (BES), chromatin remodeling factor CHD3, LEAFY (LFY), and myeloblastosis (MYB) were significantly enriched. ChIP-qPCR analysis showed that transmembrane protein 14C had the highest levels of expression, followed by fatty acid desaturase A, and chromatin-remodeling factor CHD3. Only a few studies are available investigating how combined heat and drought stress affect plant histone modifications. Therefore, results from this study may be helpful in understanding the epigenome regulation in plants exposed to the combination of abiotic stresses.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in the GSE (Genomic Spatial Event database) repository, GSE196295 [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE196295].
Abbreviations
- DT:
-
Drought
- HT:
-
Heat
- DTHT:
-
Drought and heat stress
- GO:
-
Gene Ontology
- IGV:
-
Integrative Genomics Viewer
- qRT-PCR:
-
Quantitative Real Time PCR
- RPM:
-
Reads Per Million
References
McLaughlin SB, Kszos LA. Development of switchgrass (panicum virgatum) as a bioenergy feedstock in the United States. Biomass Bioenerg. 2005;28:515–35. https://doi.org/10.1016/j.biombioe.2004.05.006.
Lovell JT, et al. Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass. Nature. 2021;590(7846):438–44. https://doi.org/10.1038/s41586-020-03127-1.
Norton MR, et al. Plant drought survival under climate change and strategies to improve perennial grasses. A review. Agron Sustain Dev. 2016;36:29. https://doi.org/10.1007/s13593-016-0362-1.
Maqbool N, Sadiq R. Allelochemicals as growth stimulators for drought stressed maize. Am J Plant Sci. 2017;8:985–97. https://doi.org/10.4236/ajps.2017.85065.
Pinheiro C, Chaves MM. Photosynthesis and drought: can we make metabolic connections from available data? J Exp Bot. 2011;62(3):869–82. https://doi.org/10.1093/jxb/erq340.
Lipiec J, et al. Effect of drought and heat stresses on plant growth and yield: a review. Int Agrophys. 2013;27:463–77. https://doi.org/10.2478/intag-2013-0017.
Atkinson NJ, Urwin PE. The interaction of plant biotic and abiotic stresses: from genes to the field. J Exp Bot. 2012;63(10):3523–43. https://doi.org/10.1093/jxb/ers100.
Khan A, et al. Coping with drought: stress and adaptive mechanisms, and management through cultural and molecular alternatives in cotton as vital constituents for plant stress resilience and fitness. Biol Res. 2018;51:47. https://doi.org/10.1186/s40659-018-0198-z.
Lamaoui M, et al. Heat and drought stresses in crops and approaches for their mitigation. Front Chem. 2018;6:26. https://doi.org/10.3389/fchem.2018.00026.
Joo E, et al. The influence of drought and heat stress on long-term carbon fluxes of bioenergy crops grown in the Midwestern USA. Plant Cell Environ. 2016;39(9):1928–40. https://doi.org/10.1111/pce.12751.
Liu Y, et al. Assessment of drought tolerance of 49 switchgrass (Panicum virgatum) genotypes using physiological and morphological parameters. Biotechnol Biofuels. 2015;8:152. https://doi.org/10.1186/s13068-015-0342-8.
Meyer E, et al. Integrating transcriptional, metabolomic, and physiological responses to drought stress and recovery in switchgrass (Panicum virgatum L.). BMC Genomics. 2014;15(1):527. https://doi.org/10.1186/1471-2164-15-527.
Hivrale V, et al. Characterization of drought- and heat-responsive microRNAs in switchgrass. Plant Sci. 2016;242:214–23. https://doi.org/10.1016/j.plantsci.2015.07.018.
Xie F, et al. High-throughput deep sequencing shows that microRNAs play important roles in switchgrass responses to drought and salinity stress. Plant Biotechnol J. 2014;12(3):354–66. https://doi.org/10.1111/pbi.12142.
Li YF, et al. Transcriptome analysis of heat stress response in switchgrass (Panicum virgatum L.). BMC Plant Biol. 2013;13:153. https://doi.org/10.1186/1471-2229-13-153.
Cortijo S, et al. Chromatin Immunoprecipitation Sequencing (ChIP-Seq) for transcription factors and chromatin factors in Arabidopsis thaliana roots: from material collection to data analysis. Methods Mol Biol. 2018;1761:231–48. https://doi.org/10.1007/978-1-4939-7747-5_18.
Malone BM, et al. Comparison of four ChIP-Seq analytical algorithms using rice endosperm H3K27 trimethylation profiling data. PLoS One. 2011;6(9):e25260. https://doi.org/10.1371/journal.pone.0025260.
Tu X, et al. Reconstructing the maize leaf regulatory network using ChIP-seq data of 104 transcription factors. Nat Commun. 2020;11:5089. https://doi.org/10.1038/s41467-020-18832-8.
Shen X, et al. Comparative transcriptome analysis and ChIP-sequencing reveals stage-specific gene expression and regulation profiles associated with pollen wall formation in Brassica rapa. BMC Genomics. 2019;20(1):264. https://doi.org/10.1186/s12864-019-5637-x.
Liu L, et al. Evaluation of experimental design and computational parameter choices affecting analyses of ChIP-seq and RNA-seq data in undomesticated poplar trees. BMC Genomics. 2014;15(Suppl. 5):S3. https://doi.org/10.1186/1471-2164-15-S5-S3.
van Dijk K, et al. Dynamic changes in genome-wide histone H3 lysine 4 methylation patterns in response to dehydration stress in Arabidopsis thaliana. BMC Plant Biol. 2010;10:238. https://doi.org/10.1186/1471-2229-10-238.
Zong W, et al. Genome-wide profiling of histone H3K4-tri-methylation and gene expression in rice under drought stress. Plant Mol Biol. 2013;81(1–2):175–88. https://doi.org/10.1007/s11103-012-9990-2.
Song ZT, et al. Histone H3K4 methyltransferases SDG25 and ATX1 maintain heat-stress gene expression during recovery in Arabidopsis. Plant J. 2021;105(5):1326–38. https://doi.org/10.1111/tpj.15114.
Georgii E, et al. Relationships between drought, heat and air humidity responses revealed by transcriptome-metabolome co-analysis. BMC Plant Biol. 2017;17:1–23. https://doi.org/10.1186/s12870-017-1062-y.
Beacon TH, Delcuve GP, López C, et al. The dynamic broad epigenetic (H3K4me3, H3K27ac) domain as a mark of essential genes. Clin Epigenet. 2021;13:138. https://doi.org/10.1186/s13148-021-01126-1.
Hayford R, et al. Global analysis of switchgrass transcriptomes in response to combined drought and heat stress. BMC Plant Biol. 2022;22(107):1–23. https://doi.org/10.6084/m9.figshare.c.5883880.v1.
Ayyappan V, et al. Genome-wide profiling of histone modifications (H3K9me2 and H4K12ac) and gene expression in rust (Uromyces appendiculatus) inoculated common bean (Phaseolus vulgaris L.). PLoS ONE. 2015;10(7):e0132176. https://doi.org/10.1371/journal.pone.0132176.
Langmead B, Salzberg S. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9. https://doi.org/10.1038/nmeth.1923.
Xu S, et al. Spatial clustering for identification of ChIP-enriched regions (SICER) to map regions of histone methylation patterns in embryonic stem cells. Methods Mol Biol. 2014;1150:97–111. https://doi.org/10.1007/978-1-4939-0512-6_5.
Aaron R, et al. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2. https://doi.org/10.1093/bioinformatics/btq033.
Zhou D, et al. agriGO: a GO analysis toolkit for the agricultural community nucleic acids research advance access. Nucl Acids Res. 2010;2010(38):W64–70. https://doi.org/10.1093/nar/gkq310.
Mi K, et al. PANTHER version 16: a revised family classification, tree-based classification tool, enhancer regions and extensive API. Nucleic Acids Res. 2021;49(D1):D394–403. https://doi.org/10.1093/nar/gkaa1106.
Huang DW, et al. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Xiao Y, et al. A novel significance score for gene selection and ranking. Bioinformatics (Oxford, England). 2014;30(6):801–7. https://doi.org/10.1093/bioinformatics/btr671.
Schwacke R, et al. MapMan4: a refined protein classification and annotation framework applicable to multi-omics data analysis. Mol Plant. 2019;12(6):879–92. https://doi.org/10.1016/j.molp.2019.01.003.
Krzywinski M, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45. https://doi.org/10.1101/gr.092759.109.
Mark H, Brewer CA. ColorBrewer.org: an online tool for selecting colour schemes for maps. Cartographic J. 2003;40(1):27–37. https://doi.org/10.1179/000870403235002042.
Sharma MK, et al. A genome-wide survey of Switchgrass genome structure and organization. PLoS ONE. 2012;7(4):e33892. https://doi.org/10.1371/journal.pone.0033892.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.
Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4. https://doi.org/10.1093/bioinformatics/btr026.
Kim JM, et al. Alterations of lysine modifications on the histone H3 N-tail under drought stress conditions in Arabidopsis thaliana. Plant Cell Physiol. 2008;49(10):1580–8. https://doi.org/10.1093/pcp/pcn133.
Zhang Y, et al. Functional analysis of tomato CHIP ubiquitin E3 ligase in heat tolerance. Sci Rep. 2021;11(1):1713. https://doi.org/10.1038/s41598-021-81372-8.
Ma M, et al. Genome-wide analysis of transcriptome and histone modifications in Brassica napus hybrid. Front Plant Sci. 2023;27(14):1123729. https://doi.org/10.3389/fpls.2023.1123729.
Tiley GP, Burleigh JG. The relationship of recombination rate, genome structure, and patterns of molecular evolution across angiosperms. BMC Evol Biol. 2015;15:114. https://doi.org/10.1186/s12862-015-0473-3.
Zhao L, et al. Integrative analysis of reference epigenomes in 20 rice varieties. Nat Commun. 2020;11:2658. https://doi.org/10.1038/s41467-020-16457-5.
Peng Y, et al. Chromatin interaction maps reveal genetic regulation for quantitative traits in maize. Nat Commun. 2019;10:2632. https://doi.org/10.1038/s41467-019-10602-5.
Yan WH, et al. Dynamic control of enhancer activity drives stage-specific gene expression during flower morphogenesis. Nat Commun. 2019;10:16. https://doi.org/10.1038/s41467-019-09513-2.
Mehraj H, et al. Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa. PLoS ONE. 2021;16(3):e0242530. https://doi.org/10.1371/journal.pone.0242530.
Zhang Y, et al. Differentially Regulated Orthologs in Sorghum and the Subgenomes of Maize. Plant Cell. 2017;29(8):1938–51. https://doi.org/10.1105/tpc.17.00354.
Kim T, Buratowski S. Dimethylation of H3K4 by Set1 recruits the Set3 histone deacetylase complex to 5’ transcribed regions. Cell. 2009;137(2):259–72. https://doi.org/10.1016/j.cell.2009.02.045.
Liu X, et al. Regulation of histone methylation and reprogramming of gene expression in the rice inflorescence meristem. Plant Cell. 2015;27(5):1428–44. https://doi.org/10.1105/tpc.15.00201.
Wang W, et al. SET-9 and SET-26 are H3K4me3 readers and play critical roles in germline development and longevity. Elife. 2018;1(7):e34970. https://doi.org/10.7554/eLife.34970.
Benayoun BA, et al. H3K4me3 breadth is linked to cell identity and transcriptional consistency. Cell. 2014;158(3):673–88. https://doi.org/10.1016/j.cell.2014.06.027.
Elrouby N. Analysis of Small Ubiquitin-Like Modifier (SUMO) targets reflects the essential nature of protein SUMOylation and provides Insight to elucidate the role of SUMO in plant development. Plant Physiol. 2015;169(2):1006–17. https://doi.org/10.1104/pp.15.01014.
Wang M, et al. Arabidopsis TRB proteins function in H3K4me3 demethylation by recruiting JMJ14. Nat Commun. 2023;14:1736. https://doi.org/10.1038/s41467-023-37263-9.
Chen W, et al. H3K4me3 plays a key role in establishing permissive chromatin states during bud dormancy and bud break in apple. Plant J. 2022;111(4):1015–31. https://doi.org/10.1111/tpj.15868.
Huang MK, et al. Genomic features of Open Chromatin Regions (OCRs) in wild soybean and their effects on gene expressions. Genes (Basel). 2021;12(5):640. https://doi.org/10.3390/genes12050640.
Hussey SG, et al. Integrated analysis and transcript abundance modelling of H3K4me3 and H3K27me3 in developing secondary xylem. Sci Rep. 2017;7(1):3370. https://doi.org/10.1038/s41598-017-03665-1.
Fiorucci AS, et al. Arabidopsis S2Lb links AtCOMPASS-like and SDG2 activity in H3K4me3 independently from histone H2B monoubiquitination. Genome Biol. 2019;20(1):100. https://doi.org/10.1186/s13059-019-1705-4.
Yan H, et al. Chromatin state-based analysis of epigenetic H3K4me3 marks of Arabidopsis in response to dark stress. Front Genet. 2019;3(10):306. https://doi.org/10.3389/fgene.2019.00306.
Dong Y, et al. TOR represses stress responses through global regulation of H3K27 trimethylation in plants. J Exp Bot. 2023;74(5):1420–31. https://doi.org/10.1093/jxb/erac486.
Godbole AA, et al. S-adenosylmethionine synthases specify distinct H3K4me3 populations and gene expression patterns during heat stress. Elife. 2023;9(12):e79511. https://doi.org/10.7554/eLife.79511.
Zeng Z, et al. Cold stress induces enhanced chromatin accessibility and bivalent histone modifications H3K4me3 and H3K27me3 of active genes in potato. Genome Biol. 2019;20(1):123. https://doi.org/10.1186/s13059-019-1731-2.
Shen Y, et al. Involvement of JMJ15 in the dynamic change of genome-wide H3K4me3 in response to salt stress. Front Plant Sci. 2022;26(13):1009723. https://doi.org/10.3389/fpls.2022.1009723.
Moreno-Pérez AJ, et al. Genome-wide mapping of H3 Lysine 4 Trimethylation (H3K4me3) and its involvement in fatty acid biosynthesis in sunflower developing seeds. Plants (Basel). 2021;10(4):706. https://doi.org/10.3390/plants10040706.
Chen HS, et al. Genome-wide analyses of histone modifications and chromatin accessibility reveal the distinct genomic compartments in Histone the Irish potato famine pathogen Phytophthorainfestans.Han bioRxiv. 2022:480484. https://doi.org/10.1101/2022.02.18.480484.
Engelhorn J, et al. Dynamics of H3K4me3 chromatin marks prevails over H3K27me3 for gene regulation during flower morphogenesis in Arabidopsis thaliana. Epigenomes. 2017;1:8. https://doi.org/10.3390/epigenomes1020008.
Tan LM, et al. Dual recognition of H3K4me3 and DNA by the ISWI component ARID5 regulates the floral transition in Arabidopsis. Plant Cell. 2020;32(7):2178–95. https://doi.org/10.1105/tpc.19.00944.
Wang Y, et al. The Arabidopsis DREAM complex antagonizes WDR5A to modulate histone H3K4me2/3 deposition for a subset of genome repression. Proc Natl Acad Sci U S A. 2022;119(27):e2206075119. https://doi.org/10.1073/pnas.2206075119.
Moreno-Pérez AJ, et al. Genome-wide mapping of H3 Lysine 4 Trimethylation (H3K4me3) and its involvement in fatty acid biosynthesis in sunflower developing seeds. Plants (Basel). 2021;10(4):706. https://doi.org/10.3390/plants10040706.
Qian S, et al. Dual recognition of H3K4me3 and H3K27me3 by a plant histone reader SHL. Nat Commun. 2018;9:2425. https://doi.org/10.1038/s41467-018-04836-y.
Miller JL, Grant PA. The role of DNA methylation and histone modifications in transcriptional regulation in humans. Subcell Biochem. 2013;61:289–317. https://doi.org/10.1007/978-94-007-4525-4_13.
Richetto J, et al. Genome-wide DNA methylation changes in a mouse model of infection-mediated neurodevelopmental disorders. Biol Psychiatry. 2017;81(3):265–76. https://doi.org/10.1016/j.biopsych.2016.08.010.
Nawa K, et al. Histone modifications and DNA methylation act cooperatively in regulating symbiosis genes in the sea anemone Aiptasia. BMC Biol. 2022;20:265. https://doi.org/10.1186/s12915-022-01469-y.
Albihlal WS, et al. Arabidopsis HEAT SHOCK TRANSCRIPTION FACTORA1b regulates multiple developmental genes under benign and stress conditions. J Exp Bot. 2018;69(11):2847–62. https://doi.org/10.1093/jxb/ery142.
Zhang C, et al. The OsABF1 transcription factor improves drought tolerance by activating the transcription of COR413-TM1 in rice. J Exp Bot. 2017;68(16):4695–707. https://doi.org/10.1093/jxb/erx260.
Zeng Z, et al. Cold stress induces enhanced chromatin accessibility and bivalent histone modifications H3K4me3 and H3K27me3 of active genes in potato. Genome Biol. 2019;20(1):123. https://doi.org/10.1186/s13059-019-1731-2.
Dhar MK, et al. Epigenetic dynamics: role of epimarks and underlying machinery in plants exposed to abiotic stress. Int J Genomics. 2014;2014:187146. https://doi.org/10.1155/2014/187146.
Liu P, et al. The histone H3K4 demethylase JMJ16 represses leaf senescence in Arabidopsis. Plant Cell. 2018;31(2):430–43. https://doi.org/10.1105/tpc.18.00693.
Popova OV, et al. The RdDM pathway is required for basal heat tolerance in Arabidopsis. Mol plant. 2013;6(2):396–410. https://doi.org/10.1093/mp/sst023.
Agrawal L, et al. Elucidation of complex nature of PEG induced drought-stress response in rice root using comparative proteomics approach. Front Plant Sci. 2016;7:1466. https://doi.org/10.3389/fpls.2016.01466.
Hamar E. Genome-wide identification of RNA silencing-related genes and their expressional analysis in response to heat stress in barley (Hordeum vulgare L.). Biomolecules. 2020;10(6):929. https://doi.org/10.3390/biom10060929.
Chung PJ, et al. Genome-wide analyses of direct target genes of four rice NAC-domain transcription factors involved in drought tolerance. BMC Genomics. 2018;19(1):40. https://doi.org/10.1186/s12864-017-4367-1.
Richts B, Commichau FM. Underground metabolism facilitates the evolution of novel pathways for vitamin B6 biosynthesis. Appl Microbiol Biotechnol. 2021;105(6):2297–305. https://doi.org/10.1007/s00253-021-11199-w.
Kim JH, et al. Non-SMC condensin I complex subunit H mediates mature chromosome condensation and DNA damage in pancreatic cancer cells. Sci Rep. 2019;9:17889. https://doi.org/10.1038/s41598-019-54478-3.
Wen Z, et al. Over expression of a Chitinase 2 gene from Chinese Wild Strawberry improves resistance to anthracnose disease in transgenic Arabidopsis thaliana. Plant Biotechnol Rep. 2020;14:725–36. https://doi.org/10.1007/s11816-020-00648-z.
Philippe F, et al. Plant pectin acetylesterase structure and function: new insights from bioinformatic analysis. BMC Genomics. 2017;18(1):456. https://doi.org/10.1186/s12864-017-3833-0.
Zhang J, et al. Transcriptomic and proteomic analyses reveal new insight into chlorophyll synthesis and chloroplast structure of maize leaves under zinc deficiency stress. J Proteomics. 2019;199:123–34. https://doi.org/10.1016/j.jprot.2019.03.001.
Sun Y, Lectin Receptor-Like Kinases: the sensor and mediator at the plant cell surface. Switzerland: N. p., 2020. Web. https://doi.org/10.3389/fpls.2020.596301.
Glesias R, et al. Ribosomal RNA N-glycosylase activity assay of ribosome-inactivating proteins. Bio Protoc. 2017;7(6):e2180. https://doi.org/10.21769/BioProtoc.2180.
Castelán-Muñoz N, et al. MADS-box genes are key components of genetic regulatory networks involved in abiotic stress and plastic developmental responses in plants. Front Plant Sci. 2019;10:853. https://doi.org/10.3389/fpls.2019.00853.
Smaczniak C, et al. Characterization of MADS-domain transcription factor complexes in Arabidopsis flower development. Proc Natl Acad Sci USA. 2012;109(5):1560–5. https://doi.org/10.1073/pnas.1112871109.
Li Z, et al. The bHLH family member ZmPTF1 regulates drought tolerance in maize by promoting root development and abscisic acid synthesis. J Exp Bot. 2019;70(19):5471–86. https://doi.org/10.1093/jxb/erz307.
Li C, et al. The bHLH transcription factor AhbHLH112 improves the drought tolerance of peanut. BMC Plant Biol. 2021;21:540. https://doi.org/10.1186/s12870-021-03318-6.
Hao Y. Basic Helix-Loop-Helix (bHLH) transcription factors regulate a wide range of functions in Arabidopsis. International J of Mol Sci. 2021;22(13):7152. https://doi.org/10.3390/ijms22137152.
Li H, et al. The bHLH-type transcription factor AtAIB positively regulates ABA response in Arabidopsis. Plant Mol Biol. 2007;65(5):655–65. https://doi.org/10.1007/s11103-007-9230-3.
Zou J, et al. Arabidopsis LSH8 positively regulates ABA signaling by changing the expression pattern of aba-responsive proteins. International J Mol Sci. 2021;22(19):10314. https://doi.org/10.3390/ijms221910314.
Bharath P. Abscisic acid-induced stomatal closure: an important component of plant defense against abiotic and biotic stress. Front Plant Sci. 2021;12:615114. https://doi.org/10.3389/fpls.2021.615114.
Zong W, et al. Feedback regulation of ABA signaling and biosynthesis by a bZIP transcription factor targets drought-resistance-related genes. Plant Physiol. 2016;171:2810–25. https://doi.org/10.1104/pp.16.00469.
Zong W, et al. Synergistic regulation of drought-responsive genes by transcription factor OsbZIP23 and histone modification in rice. J Integr Plant Biol. 2020;62(6):723–9. https://doi.org/10.1111/jipb.12850.
Huang Y, et al. Wide grain 7 increases grain width by enhancing H3K4me3 enrichment in the OsMADS1 promoter in rice (Oryza sativa L.). Plant J. 2022;102:517–28. https://doi.org/10.1111/tpj.14646.
Wang G, et al. Overexpression of tomato chloroplast-targeted DnaJ protein enhances tolerance to drought stress and resistance to Pseudomonas solanacearum in transgenic tobacco. Plant Physiol Biochem. 2014;82:95–104. https://doi.org/10.1016/j.plaphy.2014.05.011.
Pulido P, et al. Novel DNAJ-related proteins in Arabidopsis thaliana. New Phytol. 2018;217(2):480–90. https://doi.org/10.1111/nph.14827.
Luo Y, et al. Genome-wide analysis of the rice J-protein family: identification, genomic organization, and expression profiles under multiple stresses. 3 Biotech. 2019;9(10):358. https://doi.org/10.1007/s13205-019-1880-8.
Verma, et al. The expanding world of plant J-domain proteins. Crit Rev Plant Sci. 2019;38(5–6):382–400. https://doi.org/10.1080/07352689.2019.1693716.
Harris CJ, et al. A DNA methylation reader complex that enhances gene transcription. Science. 2018;362(6419):1182–6.
Scheid R, Chen J, Zhong X. Biological role and mechanism of chromatin readers in plants. Curr Opin Plant Biol. 2021;61:102008.
Bouzroud S, et al. Auxin Response Factors (ARFs) are potential mediators of auxin action in tomato response to biotic and abiotic stress (Solanum Lycopersicum). PLoS One. 2018;13(2):e0193517. https://doi.org/10.1371/journal.pone.0193517.
Lu K, et al. Overexpression of an Arabidopsis cysteine-rich receptor-like protein kinase, CRK5, enhances abscisic acid sensitivity and confers drought tolerance. J Exp Bot. 2016;67(17):5009–27. https://doi.org/10.1093/jxb/erw266.
Tripathy BC, Oelmüller R. Reactive oxygen species generation and signaling in plants. Plant Signal Behav. 2012;7(12):1621–33. https://doi.org/10.4161/psb.22455.
Kang J, et al. PDR-type ABC transporter mediates cellular uptake of the phytohormone abscisic acid. Proc Natl Acad Sci USA. 2010;107(5):2355–60. https://doi.org/10.1073/pnas.0909222107.
Jarzyniak KM, Jasiński M. Membrane transporters and drought resistance - a complex issue. Front Plant Sci. 2014;5:687. https://doi.org/10.3389/fpls.2014.00687.
Wang T, et al. Reorganization of three-dimensional chromatin architecture in Medicago truncatula under phosphorus deficiency. J Exp Bot. 2022;74(6):2005–15. https://doi.org/10.1093/jxb/erac517.
Dong J, et al. An Arabidopsis ABC transporter mediates phosphate deficiency-induced remodeling of root architecture by modulating iron homeostasis in roots. Mol Plant. 2017;10(2):244–59. https://doi.org/10.1016/j.molp.2016.11.001.
Tahmasebi A, et al. Integration of meta-analysis, machine learning and systems biology approach for investigating the transcriptomic response to drought stress in Populus species. Sci Rep. 2023;13(1):847. https://doi.org/10.1038/s41598-023-27746-6.
Acknowledgements
We thank the Department of Energy Joint Genome Institute for pre-publication access and permission to publish this analysis using the reference genome of Panicum virgatum. The authors acknowledge the assistance of Sowndarya Karapareddy, Varsha C. Anche, and Shalini P. Etukuri at Alabama A&M University in reviewing the manuscript.
Funding
This work was supported by Delaware EPSCoR with funds from the National Science Foundation Grant EPS-0814251 and the State of Delaware to VK. These funds have been utilized in the design of the study, data collection, analysis, interpretation, and in writing the manuscript.
Author information
Authors and Affiliations
Contributions
VK designed, planned the research, and helped to edit the manuscript; MCS designed, planned the experiment and helped to edit the manuscript; VA conducted the experiments, wrote the manuscript and helped edit the manuscript; SX analyzed the data; RH designed and helped to conduct the experiments; VRS helped to design the experiment, edited, and revised the manuscript; JT helped to analyze the data and edit the manuscript; DDS designed the experiments and helped to edit the manuscript; MS helped to edit the manuscript; AT helped to edit the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Study complied with local and national regulations for using plants.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Supplemental Table 1.
List of primers selected for real-time quantitative PCR (q-RT-PCR) from ChIP-Seq.
Additional file 2: Supplemental Table 2.
ChIP-Seq analysis of differentially marked histone H3K4me3 peaks in response to DT and DTHT.
Additional file 3: Supplemental Table 3.
Average number of H3K4me3 peaks in DT and DTHT-treated samples.
Additional file 4: Supplemental Table 4.
Reads per million (RPM) of the differential peaks identified between DT vs Ctrl, DTHT vs Ctrl, and DTHT vs DT.
Additional file 5: Supplemental Table 5.
DT and DTHT-responsive peaks correlated with corresponding protein coding genes identified by RNA-Seq analysis.
Additional file 6: Supplemental Table 6.
AgriGO analysis of genes associated with differential peaks involved in biological, molecular, and cellular functions in response to DT and DTHT.
Additional file 7: Supplemental Table 7.
List of DT responsive genes overlapping DT responsive peaks (for MapMan visualization).
Additional file 8: Supplemental Table 8.
List of DTHT responsive genes overlapping DTHT responsive peaks (for MapMan visualization).
Additional file 9: Supplemental Table 9.
List of HT responsive genes overlapping HT responsive peaks (for MapMan visualization).
Additional file 10: Supplemental Figure 1.
Distribution of stress-responsive peaks based on genomic features.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ayyappan, V., Sripathi, V.R., Xie, S. et al. Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass. BMC Genomics 25, 223 (2024). https://doi.org/10.1186/s12864-024-10068-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-024-10068-w