Skip to main content

Transcriptomic analyses reveal the expression and regulation of genes associated with resistance to early leaf spot in peanut

Abstract

Objective

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.

Results

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.

Introduction

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 [1]. 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 [2] 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 [5]. Dominant gene action controls resistance components and epistasis played a major role in the inheritance of resistance components [6]. 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 [7].

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 [10], pod development [11, 12], aflatoxin resistance [13, 14], defense related genes to Ralstonia solanacearum [15], and resistance to LLS [16] 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.

Main text

Materials and methods

Plant materials and inoculation

Two peanut lines 904 and 1006 selected from a F9 RIL population [17] 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 [18] 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.

Results

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.

Fig. 1
figure 1

Comparison of DEGs after inoculation. Scatter plot a T1 (4 h) vs T3 (44 h) in the line 904 (R); b T1 vs T3 in the 1006 line (S); c lines 904 (R) vs 1006 (S)

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).

Fig. 2
figure 2

Comparisons of enriched KEGG pathway of DEGs. a T1 (4 h) vs T3 (44 h) in 904 R line, b T1 vs T3 in 1006 S line, and (c) 906 R vs 1006 S line. Dot size and color refer to gene number and the corrected P value, respectively

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.

Table 1 Key DEGs related with the resistance to early leaf spot

Discussion

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 [23] and plays the potential role of defense signaling genes in quantitative disease resistance [22]. Several reports have demonstrated the increasing resistance against pathogens via expressing PAD4, such as in Arabidopsis [24]. 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 [25]. 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. [26] 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 [27]. 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 [28]. 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 [34] and to insects in tomato and cotton [35, 36]. [37] 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. [16] 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.

Limitations

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).

Abbreviations

DEG:

Differentially expressed genes

FDR:

False discovery rate

FPKM:

Fragments per kilobase of exon per million mapped reads

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

MYB:

Myeloblastosis

PAD4:

Phytoalexin deficient 4

PDA:

Potato dextrose agar

PPEE:

Posterior probability of being equivalent expression

PPO:

Polyphenol oxidase

ROS:

Reactive oxygen species

References

  1. FAO. Foreign agricultural service: official USDA estimates. FAO 2016; http://apps.usda.gov/psdonline/psdQuery.aspx.

  2. 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.

    Article  Google Scholar 

  3. 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.

    Article  Google Scholar 

  4. Wynne JC, Beute MK, Nigam SN. Breeding for disease resistance in peanut (Arachis hypogaea L.). Annu Rev Phytopathol. 1991;29:279–303.

    Article  Google Scholar 

  5. 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.

    Article  Google Scholar 

  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.

    Article  CAS  Google Scholar 

  7. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 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.

    Article  Google Scholar 

  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.

    Article  CAS  Google Scholar 

  10. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 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.

    Article  CAS  Google Scholar 

  15. 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.

    Article  Google Scholar 

  16. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. 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.

    Article  Google Scholar 

  18. Thomma BP, Nurnberger T, Joosten MH. Of PAMPs and effectors: the blurred PTI-ETI dichotomy. Plant Cell. 2011;23:4–15.

    Article  CAS  Google Scholar 

  19. 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.

    Article  CAS  Google Scholar 

  20. 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.

    Article  CAS  Google Scholar 

  21. 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.

    Article  CAS  Google Scholar 

  22. 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.

    Article  CAS  Google Scholar 

  23. 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.

    Google Scholar 

  24. Juge N. Plant protein inhibitors of cell wall degrading enzymes. Trend Plant Sci. 2006;11(7):359–67.

    Article  CAS  Google Scholar 

  25. 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.

    Article  CAS  Google Scholar 

  26. 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.

    Article  PubMed  Google Scholar 

  27. Constabel CP, Ryan CA. A survey of wound- and methyl jasmonate-induced leaf polyphenol oxidase in crop plants. Phytochemistry. 1998;47:507–11.

    Article  CAS  Google Scholar 

  28. 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.

    Article  CAS  Google Scholar 

  29. 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.

    Article  CAS  Google Scholar 

  30. 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.

    CAS  Google Scholar 

  31. 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.

    CAS  Google Scholar 

  32. Li L, Steffens JC. Overexpression of polyphenol oxidase in transgenic tomato plants results in enhanced bacterial disease resistance. Planta. 2002;2:239–47.

    Article  Google Scholar 

  33. 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.

    Article  CAS  Google Scholar 

  34. Thipyapong P, Hunt MD, Steffens JC. Antisense downregulation of polyphenol oxidase results in enhanced disease susceptibility. Planta. 2004;220:105–17.

    Article  CAS  Google Scholar 

  35. Mahanil S, Attajarusit J, Stout MJ, Thipyapong P. Overexpression of tomato polyphenol oxidase increases resistance to common cutworm. Plant Sci. 2008;174:456–66.

    Article  CAS  Google Scholar 

  36. 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.

    Article  CAS  Google Scholar 

  37. 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.

    Article  Google Scholar 

Download references

Acknowledgements

Authors appreciate Drs. T. Isleib, C. Holbrook, P. Ozias-Akins, and Y. Chu developed the RIL population and provided materials in this study.

Funding

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).

Author information

Authors and Affiliations

Authors

Contributions

GH Conceptualized the idea for this study. LG performed the fungus culture and inoculation. AH conducted phenotyping in field. LG, SH, MY, and XM analyzed the date. GH, LG, and SH prepared original draft. GH, LG, SH, MY, XM, and AH reviewed and edited the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Guohao He.

Ethics declarations

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

Not applicable.

Competing interests

The authors declare that they have 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: Figure S1.

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.

Additional file 2: Table S1.

List of 133 differentially expressed genes identified between the R line and the S line.

Additional file 3: Figure S2.

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.

Additional file 4: Table S2.

DNA sequences of resistance genes to early leaf spot disease in peanut.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13104-020-05225-9

Keywords