- Research note
- Open Access
Transcriptomic analyses reveal the expression and regulation of genes associated with resistance to early leaf spot in peanut
BMC Research Notes volume 13, Article number: 381 (2020)
Early leaf spot (ELS) caused by Cercospora arachidicola (Hori) is a serious foliar disease in peanut worldwide, which causes considerable reduction of yield. Identification of resistance genes is important for both conventional and molecular breeding. Few resistance genes have been identified and the mechanism of defense responses to this pathogen remains unknown.
We detected several genes involved in disease resistance to ELS through transcriptome analysis. Using RNA-seq technology, one hundred thirty-three differentially expressed genes (DEGs) were identified between resistant and susceptible lines. Among these DEGs, coiled coil-nucleotide binding-leucine rich repeat (NLR) type resistance genes were identified as duplicated R genes on the chromosome B2. Peanut phytoalexin deficient 4 (PAD4) regulator of effector-triggered immunity mediated by NLR resistance proteins and polyphenol oxidase (PPO) genes play important roles in early leaf spot resistance. Our study provides the useful information on plant response to C. arachidicola infection in peanut. The results suggest that a few major genes and several factors mediate the resistance to ELS disease, showing the characteristics of quantitative trait in defense responses.
Peanut (Arachis hypogaea L.), originated in South America, is a major oilseed crop and is cultivated in tropical and subtropical regions of the world. Peanut production around the world totaled ~ 38 million tons in 2015 . However, peanut productivity is severely constrained by several biotic and abiotic stresses. Amongst those stresses, early leaf spot (ELS) caused by Cercospora arachidicola Hori and late leaf spot (LLS) caused by Cercosporidium personatum (Berk & Curtis) are serious foliar diseases in peanut worldwide. Peanut yield losses may reach up to 50% worldwide due to these two diseases  and their management is highly dependent on multiple fungicide applications. The best strategy to avoid the significant economic loss is to develop resistant cultivars.
The genetic nature of resistance to early and late leaf spot diseases is quantitative and inherited independently [3, 4]. Resistance to ELS is associated with resistance components, such as longer latent period, reduced sporulation, smaller lesion diameter, and lower infection frequencies . Dominant gene action controls resistance components and epistasis played a major role in the inheritance of resistance components . Sequencing of the genome of C. arachidicola provides relevant information for the advancement of resistant cultivars and is a useful resource to aid in the selection of target resistance genes for enhanced disease control .
RNA sequencing (RNA-seq) is an effective technology for understanding of metabolic processes, genome-wide quantification of gene expressions, and identification of key genes associated with the traits of interest. Recently, RNA-seq analysis provides insights into transcriptome profiles associated with several traits in peanut. Transcriptomes have revealed the mechanism of salinity resistance [8, 9], nodulation , pod development [11, 12], aflatoxin resistance [13, 14], defense related genes to Ralstonia solanacearum , and resistance to LLS  in peanut. However, transcriptome profiling analysis of resistance to ELS has been limited. To explore the mechanism of resistance to ELS in peanut, the transcriptional profiles of two sister inbred lines, resistant (R) and susceptible (S) lines to ELS, from a F9 RIL population were investigated under C. arachidicola infection. The objectives of this study were to (1) compare the differentially expressed genes (DEGs) between R and S lines, (2) understand the possible roles of DEGs in plant defense response, (3) identify resistance genes to ELS disease.
Materials and methods
Plant materials and inoculation
Two peanut lines 904 and 1006 selected from a F9 RIL population  derived from a cross of Florida-07 × GP-NC WS 16 were used in this study. The parent GP-NC WS 16 is resistant and Florida-07 is susceptible to ELS, respectively. Early leaf spot disease was evaluated in the E.V. Smith Research Center, Tallassee, Auburn University (2016 and 2017). The inbred line 904 was selected due to its higher resistance level similar to the resistant parent and the sister line1006 demonstrated higher susceptibility to ELS compared to the other lines in the F9 RIL population. Nine seeds of each line were individually planted in 40 Magenta GA-7 plant culture boxes containing a potting mix.
Fungus culture and inoculation
A total volume of 0.5 ml C. arachidicola conidial suspension was applied to the foliage of each plant by spraying. The inoculated plants were kept at room temperature with 100% humidity in the closed culture box. Four leaves from each plant were collected at three time points, T1 (4 h), T2 (24 h), and T3 (44 h) after inoculation, with three biological replicates. RNA isolation from 18 leaf samples (2 lines × 3 time points × 3 replicates), library construction and RNA sequencing were performed in the Beijing Genomics Institute (BGI).
Identification of differentially expressed genes (DEGs), hierarchical clustering analysis, gene ontology (GO) enrichment analysis, and KEGG enrichment analysis
The differentially expressed genes (DEGs) were identified using EBSeq algorithms  based on the gene expression level with cutoff: fold change ≥ 2.00 and posterior probability of being equivalent expression (PPEE) ≤ 0.05. DEGs were compared between lines and within a line at different time points by Heatmap plot and Hierarchical clustering analysis. Gene ontology (GO) and KEGG annotation were performed using phyper, a function of R, with pvalue calculating formula in hypergenometric test (resource/hypergenometric.jpg), then calculate false discovery rate (FDR) for each p-value. The FDR < 0.01 was defined as significant enriched.
RNA-seq and transcriptome profiles between R and S lines
To compare gene expression of plants infected by the C. arachidicola, the resistant and susceptible lines were used and RNA sequencing of 18 samples, including samples from the R line and the S line at T1, T2, T3 time points with three biological replicates, were used for transcriptome analysis. All clean reads were subjected to genome mapping. On average, 83.42% reads were mapped to the reference genome and 66.48% reads of each sample were mapped to only one location of the reference genome. The uniformity of the mapping results for each sample suggested that the samples were comparable.
Identification of differentially expressed genes (DEGs)
A total of 1711 DEGs were detected from all three time points between two lines. Among these DEGs, 595 were up-regulated and 1116 were down-regulated. The resultant data showed that the number of down-regulated DEGs was significantly higher than those of up-regulated DEGs in both lines (Additional file 1: Figure S1). Due to gene expression dynamically changing between T1 (4 h) and T3 (44 h), the comparison of DEGs focused on T1 vs T3. When comparing R and S lines at T3, 133 DEGs including 52 up and 81 down-regulated were identified according to gene expression level with fold change ≥ 2 (Additional file 2: Table S1). The summary of DEGs was illustrated in the Fig. 1.
Gene ontology and KEGG enrichment analysis of DEGs
Gene ontology (GO) enrichment analysis was performed to analyze the functions of DEGs. Metabolic process with 137 and 25 genes were dominant in the biological process category in R and S lines, respectively. Catalytic activity (181 genes) and binding (102 genes) were significant active responses in the molecular function category in the R line, whereas only 33 and 25 genes with the same function in the S line, respectively, indicating there was a higher number of expressed genes in R line than in S line between T1 and T3 (Additional file 3: Figure S2).
To further explore biological functions and gene interaction, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was used. For the R line, comparison of DEGs with T1 vs T3 showed 11 KEGG pathways were significantly enriched, including biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, MAPK signaling pathway-plant, flavonoid biosynthesis, plant-pathogen interaction, circadian rhythm-plant, isoflavonoid biosynthesis, etc. While in the S line, 4 KEGG pathways including photosynthesis, plant-pathogen interaction, oxidative phosphorylation, MAPK signaling pathway-plant were significantly enriched for DEGs between T1 and T3 (Fig. 2).
DEGs related the resistance to ELS
There were 133 DEGs with 52 up and 81 down-regulated genes identified between R and S lines (Additional file 2: Table S1). Some DEGs were crucially involved in plant responses to the pathogen infection. For instance, nine of such genes were significantly up-regulated in the S line while down-regulated in the R line based on the log2 fold change, including two nucleotide binding-leucine rich repeats (NLRs), Interleukin-1 receptor-associated kinase 4, ATP-dependent RNA helicase, Nuclear pore complex protein Nup85, and E3 ubiquitin-protein ligase. On the other hand, 10 different DEGs including phytoalexin deficient 4 (PAD4), Interleukin-1 receptor-associated kinase 4, ATP-dependent RNA helicase, ATP-binding cassette, Serine/threonine-protein kinase, Pectinesterase, Auxin responsive GH3 gene family, and Polyphenol oxidase were significantly up-regulated in the R line but down-regulated in the S line (Table 1). The sequence of putative NLRs were showed in Additional file 4: Table S2.
The development of disease-resistant plants is important in mitigating the losses of productions. Understanding the mechanism of plant-pathogen interactions and searching for resistant resources is a prerequisite in plant breeding. Plant responses to pathogen attack are involved in several defense strategies against pathogens. The first barrier for pathogen invasion is the uses of plant reprogramming the cell wall damaged by virulent pathogens, resulting in pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) [19, 20]. This basal resistance layer is controlled by lipase-like protein phytoalexin deficient 4 (PAD4) either in combination with enhanced disease susceptibility 1 (EDS1) or independent of EDS1 [21, 22]. PAD4 is required for activation of responses during some other gene-for-gene resistance responses  and plays the potential role of defense signaling genes in quantitative disease resistance . Several reports have demonstrated the increasing resistance against pathogens via expressing PAD4, such as in Arabidopsis . In this study, PAD4 was significantly down-regulated in the S line compared to the R line, suggesting that the lower expression level of PAD4 might affect the resistance to ELS in the susceptible genotype.
Because PAD4 is an important component of the downstream signaling in a different immune pathway, the reduced expression of PAD4 may also have a negative effect on resistance pathways. Several other genes participated in various pathways responding to ELS pathogen challenges were identified in this study. These genes encoded various enzymes to conspire in biological processes. For instance, expression of pectinesterase gene was significantly reduced in the S line that could result in weak cell walls due to degradation of pectin, a critical cell wall component. In plants, pectinesterase plays a role in plant response to pathogen attack by influencing numerous physiological processes . However, the repressed expression of pectinesterase gene led to degrade cell wall integrity, thus increase C. arachidicola colonization of leaf tissue and enhance disease development.  described that the impact of pectin methyl esterification on plant-pathogen interactions and on the dynamic role of its alteration during pathogenesis. If these passive defenses were breached, the plant innate immune system would launch active defenses . Nevertheless, due to the down-regulation of PAD4 in the S line, it lost the functions of the signaling accumulation for downstream activities.
Besides PAD4, polyphenol oxidase (PPO) genes also significantly down-regulated in the S line. Normally, PPO is up-regulated by abiotic and biotic stresses though the responses to stresses varies within PPO gene families by plant species . Overexpression of PPO displayed enhanced resistance to pathogen in wheat, chickpea, tomato, and populus [29,30,31,32,33]. Conversely, down-regulation of PPO resulted in increased disease susceptibility. For instance, silencing of PPOs leads to increased susceptibility to disease in tomato  and to insects in tomato and cotton [35, 36].  studied the interactions of host with fungi in pearl millet and reported that the resistant genotypes showed high, rapid accumulation of localized levels of PPO while susceptible genotypes failed to accumulate PPO even after consideration time. In this study, significant down-regulation of PPO in the S line could enhance susceptibility to ELS disease.
In addition to PTI, another type of defense response is activated by resistance proteins as sensors that indirectly or directly recognize specific effectors from many different pathogens, leading to effector-triggered immunity (ETI). The resistance genes, NLR with N-terminal Toll-interleukin-1-receptor domain or coiled coil domain encoding resistance proteins, are a major class of resistance genes that evolved to intercept these effectors. Two putative NLRs were identified as DEGs between the R line and S line in this study. Two R genes possessed the same DNA sequences with several SNPs and were closely located in the B2 chromosome. Higher R gene expression at the earlier reaction period in the susceptible genotype implied that resistance genes trigged by pathogen infection responded quicker in the susceptible than the resistant genotype.  observed the same scenario in a different peanut population responding to LLS disease. The susceptible genotype was eventually failed in defense response might be due to the lower expression of PAD4. Overexpression of PAD4 would allow us to have further insight into the plant response in the susceptible genotype.
The resistance genes identified from differentially expressed genes still need to validate by functional study, such as loss of function mutation targeting on these genes.
Availability of data and materials
All raw fastq sequences for these eighteen libraries were deposited in the National Center of Biotechnology Information (NCBI). The BioProject number is PRJNA494675 (https://www.ncbi.nlm.nih.gov/sra/PRJNA494675).
Differentially expressed genes
False discovery rate
Fragments per kilobase of exon per million mapped reads
Kyoto encyclopedia of genes and genomes
Phytoalexin deficient 4
Potato dextrose agar
Posterior probability of being equivalent expression
Reactive oxygen species
FAO. Foreign agricultural service: official USDA estimates. FAO 2016; http://apps.usda.gov/psdonline/psdQuery.aspx.
Tshilenge LL, Nkongolo KKC, Kalonji MA, Kizungu RV. Epidemiology of the groundnut (Arachis hypogaea L.) leaf spot disease: genetic analysis and developmental cycles. Am J Plant Sci. 2012;3:582–8.
Anderson WF, Wynne JC, Green CC. Potential for incorporation of early and late leaf spot resistance in peanut. Plant Breeding. 1986;97(2):163–70.
Wynne JC, Beute MK, Nigam SN. Breeding for disease resistance in peanut (Arachis hypogaea L.). Annu Rev Phytopathol. 1991;29:279–303.
Waliyar F, McDonald D, Subba Rao PV, Reddy PM. Components of resistance to an India source of Cercospora arachidicola in selected peanut lines. Peanut Sci. 1993;20(2):93–6.
Green CC, Wynne JC. Diallel and generation means analyses for the components of resistance to Cercospora arachidicola in peanut. Theor Appl Genet. 1986;73(2):228–35.
Orner VA, Cantonwine EG, Wang XM, Abouelleil A, Bochicchio J, Nusbaum C, et al. Draft genome sequence of Cercospora arachidicola, causal agent of early leaf spot in peanuts. Genome Announc. 2015;3(6):e01281–15. https://doi.org/10.1128/genomeA.01281-15.
Sui JM, Li G, Chen GX, Zhao CM, Kong XY, Hou XM, et al. RNA-seq analysis reveals the role of a small GTP-binding protein, Rab7, in regulating clathrin-mediated endocytosis and salinity-stress resistance in peanut. Plant Biotechnol Rep. 2017;11:43–52. https://doi.org/10.1007/s11816-017-0428-9.
Sui JM, Jiang PP, Qin GL, Gai SP, Zhu D, Qiao LX, Wang JS. Transcriptome profiling and digital gene expression analysis of genes associated with salinity resistance in peanut. Electronic J Biotechnol. 2018;32:19–25. https://doi.org/10.1016/j.ejbt.2017.12.002.
Peng Z, Liu FX, Wang LP, Zhou H, Paudel D, Tan LB, et al. Transcriptome profiles reveal gene regulation of peanut (Arachis hypogaea L.) nodulation. Sci Rep. 2017;7:40066. https://doi.org/10.1038/srep40066.
Yang S, Li L, Zhang JL, Geng Y, Guo F, Wang JG, et al. Transcriptome and differential expression profiling analysis of the mechanism of Ca2+ regulation in peanut (Arachis hypogaea) pod development. Front Plant Sci. 2017;8:1609. https://doi.org/10.3389/fpls.2017.01609.
Wan LY, Li B, Lei Y, Yan LY, Ren XP, Chen YN, et al. Mutant transcriptome sequencing provides insights into pod development in peanut (Arachis hypogaea L.). Front Plant Sci. 2017;8:1900. https://doi.org/10.3389/fpls.2017.01900.
Wang HM, Lei Y, Yan LY, Wan LY, Ren XP, Chen SL, et al. Functional genomic analysis of Aspergillus flavus interacting with resistant and susceptible peanut. Toxins. 2016;8:46. https://doi.org/10.3390/toxins8020046.
Korani W, Chu Y, Holbrook CC, Ozias-Akins P. Insight into genes regulating post-harvest aflatoxin contamination of tetraploid peanut from transcriptional profiling. Genetics. 2018;209(1):143–56.
Chen YN, Ren XP, Zhou XJ, Huang L, Yan LY, Lei Y, et al. Dynamics in the resistant and susceptible peanut (Arachis hypogaea L.) root transcriptome on infection with the Ralstonia solanacearum. BMC Genomics. 2014;15:1078.
Han SY, Liu H, Yan M, Qi FY, Wang YQ, Sun ZQ, et al. Differential gene expression in leaf tissues between mutant and wild-type genotypes response to late leaf spot in peanut (Arachis hypogaea L.). PLoS ONE. 2017;12(8):e0183428. https://doi.org/10.1371/journal.pone.0183428.
Holbrook CC, Isleib TG, Ozias-Akins P, Chu Y, Knapp SJ, Tillman B, et al. Development and pehnotyping of recombinant inbred line (RIL) populations for peanut (Arachis hypogaea). Peanut Sci. 2013;40:89–94.
Thomma BP, Nurnberger T, Joosten MH. Of PAMPs and effectors: the blurred PTI-ETI dichotomy. Plant Cell. 2011;23:4–15.
Silva MS, Arraes FBM, de Araujo Campos M, Grossi-de-Sa M, Fernandez D, de Souza Candido E, Cardoso MH, Franco OL, Grossi-de-Sa MF. Review: potential biotechnological assets related to plant immunity modulation applicable in engineering disease-resistant crops. Plant Sci. 2018;270:72–84.
Wagner S, Stuttmann J, Rietz S, Guerois R, Brunstein E, Bautor J, Niefind K, Parker JE. Structural basis for signaling by exclusive EDS1 heteromeric complexes with SAG101 or PAD4 in plant innate immunity. Cell Host Microbe. 2013;14:619–30.
Youssef RM, MacDonald MH, Brewer EP, Bauchan GR, Kim KH, Matthews BF. Ectopic expression of AtPAD4 broadens resistance of soybean to soybean cyst and root-knot nematodes. BMC Plant Biol. 2013;13:67.
Zhou N, Tootle TL, Tsui F, Klessig DF, Glazebrook J. PAD4 functions upstream from salicylic acid to control defense responses in Arabidopsis. Plant Cell. 1998;10:1021–30.
Harun-Or-Rashid M, Khan A, Hossain MT, Chung YR. Induction of systemic resistance against aphids by endophytic Bacillus velezensis YC7010 via expressing phytoalexin deficient 4 in Arabidopsis. Front Plant Sci. 2017;8:211.
Juge N. Plant protein inhibitors of cell wall degrading enzymes. Trend Plant Sci. 2006;11(7):359–67.
Lionetti V, Cervone F, Bellincampi D. Methyl esterification of pectin plays a role during plant-pathogen interactions and affects plant resistance to diseases. J Plant Physiol. 2012;169(16):1623–30.
Klessig DF, Choi HW, Dempsey DA. Systemic acquired resistance and salicylic acid: past, present, and future. MPMI. 2018. https://doi.org/10.1094/MPMI-03-19-0067-CR.
Constabel CP, Ryan CA. A survey of wound- and methyl jasmonate-induced leaf polyphenol oxidase in crop plants. Phytochemistry. 1998;47:507–11.
Araji S, Grammer TA, Gertzen R, Anderson SD, Mikulic-Petkovsek M, Veberic R, Phu ML, Solar A, Leslie CA, Dandekar AM, Escobar MA. Novel roles for the polyphenol oxidase enzyme in secondary metabolism and the regulation of cell death in walnut. Plant Physiol. 2014;164:1191–203.
Mohammadi M, Kazemi H. Changes in peroxidase and polyphenol oxidase activities in susceptible and resistant wheat heads inoculated with Fusarium graminearum and induced resistance. Plant Sci. 2002;162:491–8.
Raju S, Jayalakshmi SK, Sreeramulu K. Comparative study on the induction of defense related enzymes in two different cultivars of chickpea (Cicer arietinum L.) genotypes by salicylic acid, spermin and Fusarium oxysporum f. sp. ciceri. Aust J Crop Sci. 2008;2:121–40.
Jayalakshmi SK, Raju S, Rani SU, Benagi VI, Sreeramulu K. Trichoderma harzianum L1 as a potential source for lytic enzymes and elicitor of defense responses in chickpea (Cicer arietinum L.) against wilt disease caused by Fusarium oxysporum f. sp. ciceri. Aust J Crop Sci. 2009;3:44–52.
Li L, Steffens JC. Overexpression of polyphenol oxidase in transgenic tomato plants results in enhanced bacterial disease resistance. Planta. 2002;2:239–47.
Wang J, Constabel CP. Polyphenol oxidase overexpression in transgenic Populus enhances resistance to herbivory by forest tent caterpillar (Malacosoma disstria). Planta. 2004;220:87–96.
Thipyapong P, Hunt MD, Steffens JC. Antisense downregulation of polyphenol oxidase results in enhanced disease susceptibility. Planta. 2004;220:105–17.
Mahanil S, Attajarusit J, Stout MJ, Thipyapong P. Overexpression of tomato polyphenol oxidase increases resistance to common cutworm. Plant Sci. 2008;174:456–66.
Bhonwong A, Stout MJ, Attajarusit J, Tantasawat P. Defensive role of tomato polyphenol oxidases against cotton bollworm (Helicoverpa armigera) and beet armyworm (Spodoptera exigua). J Chem Ecol. 2009;35:28–38.
Niranjan RS, Sarosh BR, Shetty HS. Induction and accumulation of polyphenol oxidase activities as implicated in development of resistance against pearl millet downy mildew disease. Funct Plant Biol. 2006;33:563–71.
Authors appreciate Drs. T. Isleib, C. Holbrook, P. Ozias-Akins, and Y. Chu developed the RIL population and provided materials in this study.
This work was support by USDA/NIFA (award number 2014-38821-22453) and the George Washington Carver Agricultural Experiment Station at Tuskegee University; National Natural Science Foundation of China (No. 31871663); and Major R&D and Promotion Projects in Henan, China (No. 182102110137).
Ethics approval and consent to participate
Peanut seeds of two lines were obtained from Auburn University. After germinated, only seedlings were used in the lab for this research.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of DEGs among different groups. X-axis represents comparison of DEGs between each group and Y-axis refers to the number of DEG. Red color indicates up-regulated DEGs and blue color for down-regulated DEGs.
List of 133 differentially expressed genes identified between the R line and the S line.
Gene ontology (GO) enrichment analysis showed DEGs involved in three function categories. A. T1 (4 h) vs T3 (44 h) in R line; B. T1 vs T3 in S line; C. R vs S lines.
DNA sequences of resistance genes to early leaf spot disease in peanut.
About this article
Cite this article
Gong, L., Han, S., Yuan, M. et al. Transcriptomic analyses reveal the expression and regulation of genes associated with resistance to early leaf spot in peanut. BMC Res Notes 13, 381 (2020). https://doi.org/10.1186/s13104-020-05225-9
- Defense response
- Differentially expressed gene
- Gene regulation