Histone modifications such as methylation and acetylation play a significant role in controlling gene expression in unstressed and stressed plants. Genome-wide analysis of such stress-responsive modifications and genes in non-model crops is limited. We report the genome-wide profiling of histone methylation (H3K9me2) and acetylation (H4K12ac) in common bean (Phaseolus vulgaris L.) under rust (Uromyces appendiculatus) stress using two high-throughput approaches, chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq). ChIP-Seq analysis revealed 1,235 and 556 histone methylation and acetylation responsive genes from common bean leaves treated with the rust pathogen at 0, 12 and 84 hour-after-inoculation (hai), while RNA-Seq analysis identified 145 and 1,763 genes differentially expressed between mock-inoculated and inoculated plants. The combined ChIP-Seq and RNA-Seq analyses identified some key defense responsive genes (calmodulin, cytochrome p450, chitinase, DNA Pol II, and LRR) and transcription factors (WRKY, bZIP, MYB, HSFB3, GRAS, NAC, and NMRA) in bean-rust interaction. Differential methylation and acetylation affected a large proportion of stress-responsive genes including resistant (R) proteins, detoxifying enzymes, and genes involved in ion flux and cell death. The genes identified were functionally classified using Gene Ontology (GO) and EuKaryotic Orthologous Groups (KOGs). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified a putative pathway with ten key genes involved in plant-pathogen interactions. This first report of an integrated analysis of histone modifications and gene expression involved in the bean-rust interaction as reported here provides a comprehensive resource for other epigenomic regulation studies in non-model species under stress.

BackgroundCommon bean (Phaseolus vulgaris L.) is one of the most important legumes in the world. Several diseases severely reduce bean production and quality; therefore, it is very important to better understand disease resistance in common bean in order to prevent these losses. More than 70 resistance (R) genes which confer resistance against various pathogens have been cloned from diverse plant species. Most R genes share highly conserved domains which facilitates the identification of new candidate R genes from the same species or other species. The goals of this study were to isolate expressed R gene-like sequences (RGLs) from 454-derived transcriptomic sequences and expressed sequence tags (ESTs) of common bean, and to develop RGL-tagged molecular markers.ResultsA data-mining approach was used to identify tentative P. vulgaris R gene-like sequences from approximately 1.69 million 454-derived sequences and 116,716 ESTs deposited in GenBank. A total of 365 non-redundant sequences were identified and named as common bean (P. vulgaris = Pv) resistance gene-like sequences (PvRGLs). Among the identified PvRGLs, about 60% (218 PvRGLs) were from 454-derived sequences. Reverse transcriptase-polymerase chain reaction (RT-PCR) analysis confirmed that PvRGLs were actually expressed in the leaves of common bean. Upon comparison to P. vulgaris genomic sequences, 105 (28.77%) of the 365 tentative PvRGLs could be integrated into the existing common bean physical map. Based on the syntenic blocks between common bean and soybean, 237 (64.93%) PvRGLs were anchored on the P. vulgaris genetic map and will need to be mapped to determine order. In addition, 11 sequence-tagged-site (STS) and 19 cleaved amplified polymorphic sequence (CAPS) molecular markers were developed for 25 unique PvRGLs.ConclusionsIn total, 365 PvRGLs were successfully identified from 454-derived transcriptomic sequences and ESTs available in GenBank and about 65% of PvRGLs were integrated into the common bean genetic map. A total of 30 RGL-tagged markers were developed for 25 unique PvRGLs, including 11 STS and 19 CAPS markers. The expressed PvRGLs identified in this study provide a large sequence resource for development of RGL-tagged markers that could be used further for genetic mapping of disease resistant candidate genes and quantitative trait locus/loci (QTLs). This work also represents an additional method for identifying expressed RGLs from next generation sequencing data.

Common bean (Phaseolus vulgaris L.) is an important legume, useful for its high protein and dietary fiber. The fungal pathogen Uromyces appendiculatus (Pers.) Unger can cause major loss in susceptible varieties of the common bean. The Ur-3 locus provides race specific resistance to virulent strains or races of the bean rust pathogen along with Crg, (Complements resistance gene), which is required for Ur-3-mediated rust resistance. In this study, we inoculated two common bean genotypes (resistant “Sierra” and susceptible crg) with rust race 53 of U. appendiculatus, isolated leaf RNA at specific time points, and sequenced their transcriptomes. First, molecular markers were used to locate and identify a 250 kb deletion on chromosome 10 in mutant crg (which carries a deletion at the Crg locus). Next, we identified differential expression of several disease resistance genes between Mock Inoculated (MI) and Inoculated (I) samples of “Sierra” leaf RNA within the 250 kb delineated region. Both marker assisted molecular profiling and RNA-seq were used to identify possible transcriptomic locations of interest regarding the resistance in the common bean to race 53. Identification of differential expression among samples in disease resistance clusters in the bean genome may elucidate significant genes underlying rust resistance. Along with preserving favorable traits in the crop, the current research may also aid in global sustainability of food stocks necessary for many populations.

Histone deacetylases (HDACs) play an important role in plant growth, development, and defense processes and are one of the primary causes of epigenetic modifications in a genome. There was only one study reported on epigenetic modifications of the important legume crop, common bean, and its interaction with the fungal rust pathogen Uromyces appendiculatus prior to this project. We measured the total active HDACs levels in leaf tissues and observed expression patterns for the selected HDAC genes at 0, 12, and 84 hours after inoculation in mock inoculated and inoculated plants. Colorimetric analysis showed that the total amount of HDACs present in the leaf tissue decreased at 12 hours in inoculated plants compared to mock inoculated control plants. Gene expression analyses indicated that the expression pattern of gene PvSRT1 is similar to the trend of total active HDACs in this time course experiment. Gene PvHDA6 showed increased expression in the inoculated plants during the time points measured. This is one of the first attempts to study expression levels of HDACs in economically important legumes in the context of plant pathogen interactions. Findings from our study will be helpful to understand trends of total active HDACs and expression patterns of these genes under study during biotic stress.