Skip to main content

High-throughput diagnostic markers for foliar fungal disease resistance and high oleic acid content in groundnut

A Correction to this article was published on 25 April 2024

This article has been updated

Abstract

Background

Foliar diseases namely late leaf spot (LLS) and leaf rust (LR) reduce yield and deteriorate fodder quality in groundnut. Also the high oleic acid content has emerged as one of the most important traits for industries and consumers due to its increased shelf life and health benefits.

Results

Genetic mapping combined with pooled sequencing approaches identified candidate resistance genes (LLSR1 and LLSR2 for LLS and LR1 for LR) for both foliar fungal diseases. The LLS-A02 locus housed LLSR1 gene for LLS resistance, while, LLS-A03 housed LLSR2 and LR1 genes for LLS and LR resistance, respectively. A total of 49 KASPs markers were developed from the genomic regions of important disease resistance genes, such as NBS-LRR, purple acid phosphatase, pentatricopeptide repeat-containing protein, and serine/threonine-protein phosphatase. Among the 49 KASP markers, 41 KASPs were validated successfully on a validation panel of contrasting germplasm and breeding lines. Of the 41 validated KASPs, 39 KASPs were designed for rust and LLS resistance, while two KASPs were developed using fatty acid desaturase (FAD) genes to control high oleic acid levels. These validated KASP markers have been extensively used by various groundnut breeding programs across the world which led to development of thousands of advanced breeding lines and few of them also released for commercial cultivation.

Conclusion

In this study, high-throughput and cost-effective KASP assays were developed, validated and successfully deployed to improve the resistance against foliar fungal diseases and oleic acid in groundnut. So far deployment of allele-specific and KASP diagnostic markers facilitated development and release of two rust- and LLS-resistant varieties and five high-oleic acid groundnut varieties in India. These validated markers provide opportunities for routine deployment in groundnut breeding programs.

Peer Review reports

Introduction

Groundnut (Arachis hypogaea L.), or peanut, is cultivated in more than 100 countries worldwide covering33.2 million hectares with an annual production of 72.3 million tons [1]. A total of 100 g of groundnut oil contained 17.7 g of saturated fat, 48.3 g of monounsaturated fat (primarily oleic acid), and 33.4 g of polyunsaturated fat (comprising linoleic and linolenic acid). High oleic acid increases the shelf life of groundnut oil and has important health benefits [2]. The popularity of groundnut oil in Chinese, South Asian and Southeast Asian countries is because of its high smoke point, making it ideal for preparation of fried food delicacies.

Foliar fungal diseases, namely, rust caused by Puccinia arachidis and late leaf spot (LLS) caused by Nothopassalora personata, are highly destructive when they occur together, leading to substantial yield losses of 50–70%, as well as deterioration of fodder quality [3]. The rate of photosynthesis decreases as the healthy leaf area decreases due to the large number of leaf spots and defoliation. In addition to yield loss and reduced fodder quality, farmers bear the extra cost of fungicide application which impacts both their income and environmental health. As a result, developing groundnut varieties with rust and LLS resistance is a key objective for groundnut breeding programs worldwide to sustain productivity .

The marker-assisted backcrossing (MABC) approach has been successfully used to improve resistance to rust and late leaf spot as well as to increase the oleic acid content in groundnut varieties [4]. There has been successful marker deployment in groundnut breeding using low-throughput genotyping assays such as simple sequence repeats (SSRs) and allele-specific markers for improving foliar disease resistance [2, 5,6,7,8] and high oleic acid content [2, 8,9,10,11,12]. The deployment of available allele-specific markers and KASP (Kompetitive Allele Specific PCR) assays developed under this study have already resulted in the commercial cultivation of five high-oleic acid varieties in five Indian states (Gujarat, Rajasthan, Karnataka, Andhra Pradesh, Telangana and Tamil Nadu). The two Virginia Bunch varieties, Girnar 4 (ICGV 15083) and Girnar 5 (ICGV 15090) were developed from the cross ICGV 06420 × (ICGV 06420 × Sun Oleic 95-R) while the three Spanish Bunch varieties, GG 39 (ICGV 16697), ICRC-1 (ICGV 16690) and GG 40 (ICGV 16668), were developed from the cross ICGV 06110 × (ICGV 061l0 × Sun Oleic 95-R); and the high oleic acid donor SunOleic 95R [9]. In addition, there are several promising marker-assisted bred lines (ICGV 181023, ICGV 181025 and ICGV 171025) in the genetic background of three elite varieties, ICGV 06142, ICGV 06420 and ICGV 061l0 [9] which are under third-year testing in the All India Coordinated Research Project on Groundnut (AICRP-G). Furthermore, two foliar disease-resistant varieties, namely, Improved JL 24 (DBG 3) and Super TMV 2 (DBG 4) have also been developed using the resistant donor GPBD 4 [6, 7] and released in the state of Karnataka, India for commercial cultivation. The low-throughput markers such as allele-specific and SSR markers are good for genotyping small number of breeding lines, however, genotyping and selection in large scale breeding samples require high-throughput genotyping technologies to perform selection in time. Further, the availability of high-throughput markers such as KASP assays provide an opportunity to perform selection for multiple traits which benefit breeding teams in saving resources and to stack multiple traits.

With recent advancements in sequencing technologies, high-throughput genotyping has become more feasible, facilitating the identification of associated genomic regions and candidate genes and the development of diagnostic markers for use in early-generation selection in several crops, including groundnut [4, 13, 14]. In addition to the already available reference genomes of diploids [15,16,17], the further availability of the high-density genotyping assay “Axiom_Arachis” with 58 K SNPs [18, 19] and high-quality reference genomes [20,21,22] for both subspecies of cultivated tetraploid groundnut further enhances trait dissection and marker development for target traits in groundnut.

Previous studies using a recombinant inbred line (RIL) population developed from the cross TAG 24 and GPBD 4 identified major quantitative trait loci (QTLs) for rust resistance on chromosome A03 and two major QTLs for LLS resistance on chromosomes A03 and A02 [23,24,25]. The QTL for rust explained up to 82.6% of the phenotypic variance (PVE), while both QTLs for LLS resistance explained 40–60% PVE; the A03 genomic region is common for both diseases, rust and LLS [23]. SSR markers linked to these traits have been developed and used for the introgression of rust and LLS resistance [2, 5]. However, genotyping with the currently available linked SSR markers for these traits is time consuming and laborious, hindering their large-scale adoption in breeding programs. Marker-assisted early generation selection (MEGS) has emerged as an effective approach for selecting desired alleles in early generations, such as F2 in groundnut [4, 26]. However, MEGS requires diagnostic markers with high selection efficiency so that decisions can be made at very early stages of the breeding program. Therefore, this study reports the development of low-density and high-throughput single nucleotide polymorphism (SNP)-based KASP genotyping assays for two foliar fungal diseases and high oleic acid content in groundnut. These markers are now routinely used for tracking desirable alleles for these traits in breeding populations and germplasms.

Materials and methods

Genetic mapping population and validation panel

The RIL mapping population (TAG 24 × GPBD 4) comprising 266 individuals was developed at the University of Agricultural Sciences-Dharwad (UAS-Dharwad), India. The resistant parent, GPBD 4, is an elite groundnut variety used as a national check, derived from the cross KRG 1 × ICGV 86855 (CS 16). Notably, ICGV 86855 (CS 16), an interspecific derivative of A. cardenasii, is a source of rust and LLS resistance [27]. The RIL population (TAG 24 × GPBD 4) was phenotyped extensively for both foliar fungal diseases, i.e., leaf rust and LLS, at the University of Agricultural Sciences, Dharwad, India, and the details are available in our previous studies [23,24,25]. Multiple seasons of phenotyping data have been generated on the RIL population TAG 24 × GPBD 4 for leaf rust and LLS resistance for seven years/seasons (2004 to 2010) and reported in previous studies with SSR [23] and RAD-Seq to identify the QTLs associated with rust and LLS resistance. The disease rating of the RIL population was performed at 80 days after sowing (DAS) and 90 DAS for rust resistance and at 70 DAS and 90 DAS for LLS resistance [23]. Further details on the phenotyping procedure for rust and LLS have been provided by Sujay et al. [23]. The extreme resistant and susceptible individuals from the RIL population (TAG 24 × GPBD 4) based on phenotyping data were used for developing two DNA bulks for sequencing and QTL-seq analysis to identify the genomic regions and candidate genes for leaf rust and late leaf spot resistance in groundnut [22, 24].

We used separate validation panels individually for each LLS, rust and high oleic acid content. The validation panel for LLS included 20 lines (10 resistant and 10 susceptible), 43 for rust (21 resistant and 22 susceptible) and 51 for oleic acid (25 high oleic acid and 26 with low oleic acid content). The sources of high oleic acid, SunOleic 95R, as well as the derived high oleic acid varieties Girnar 4 (ICGV 15083), Girnar 5 (ICGV 15090), GG 40 (ICGV 16688) and ICRC-1 (ICGV 16690) from SunOleic 95R, were also included in the validation panel. The abovementioned genotype sets were used for the validation of KASP markers and correlated with preliminary phenotypic data. Furthermore, based on accuracy, highly efficient KASP markers were identified for breeding purposes.

Identification of QTLs for rust and LLS resistance

Genotyping by sequencing of the RIL population

DNA from 217 RILs and the two parents was extracted using a Nucleospin Plant II kit (Macharey-Nigel, Duren, Germany). The DNA quality and quantity were checked on 0.8% agarose and then with a Qubit 2.0 fluorometer (Thermo Fisher Scientific Inc., USA). Low-coverage sequencing (GBS: genotyping-by-sequencing) was performed for the entire RIL population for simultaneous SNP discovery. A total of 10 ng of DNA from each RIL was digested using the restriction enzyme ApeKI endonuclease, which recognizes the G/CWCG site. The ligation enzyme T4 ligase was used to ligate the digested products with uniquely barcoded adapters. Such digestion and ligation were performed for each RIL, and then equal proportions from each sample were mixed to construct the libraries. These libraries were then amplified and purified to remove excess adapters. These DNA libraries were then sequenced on the HiSeq 2500 platform (Illumina Inc., San Diego, CA, USA) to generate genome-wide sequence reads at Center of Excellence in Genomics and Systems Biology (CEGSB), ICRISAT, Hyderabad, India. The sequence raw reads in FASTQ files generated for the RIL population and parental genotypes were used for SNP discovery using TASSEL v4.0 [28]. The groundnut progenitor (A. ipaensis and A. duranensis) genomes were used as reference assemblies for SNP calling [15].

Genetic mapping using a GBS-based genetic map

The GBS data generated from 217 RILs and two parental genotypes (TAG 24 and GPBD 4) were used for SNP discovery. Diploid progenitor genome assemblies (A- and B- subgenome) were used as reference assemblies for SNP calling. The chi-square (χ2) values calculated for each SNP marker were used to determine the goodness of fit to the expected 1:1 segregation ratio; highly distorted markers were filtered out and not included in the genetic map. The dense genetic map was constructed using JoinMap version 4 [29] and was redrawn using MapChart [30]. The Kosambi map function [31] with a recombination frequency of 0.45 was used to determine the map order for these new SNP markers by keeping the order fixed for earlier mapped SSR marker loci [23]. QTL analysis was performed using the composite interval mapping model implemented in Win-QTL cartographer 2.5 software [32].

Identification of genomic regions for rust and LLS resistance using QTL-Seq

For QTL-Seq analysis, DNA from 25 RILs with high susceptibility scores for rust (Rust_Sbulk) and LLS (LLS_Sbulk) but high resistance scores for rust (Rust_Rbulk) and LLS (LLS_Rbulk) were pooled separately to construct sequencing bulks. QTL-Seq analysis was performed using diploid [24] and tetraploid [22] genomes (Supplementary Table 1). These libraries generated 250-base paired-end reads on the Illumina HiSeq 2500 platform and the analysis was performed using the QTL-Seq pipeline (http://genome-e.ibrc.or.jp/home/bioinformatics-team/QTL-seq).

Candidate gene discovery and marker development

The discovery of candidate genes and allele-specific markers was performed from the QTL/genomic regions identified in three different analyses, namely, (a) GBS-based genetic mapping and QTL analysis; (b) QTL-Seq analysis using a synthetic tetraploid reference genome built from individual diploid genomes for the A- and B-sub-genomes (genomes reported by Bertioli et al. [15]; and (c) QTL-Seq analysis using a tetraploid reference genome for the fastigiata subspecies of the cultivated groundnut genome [22]. Initially, low-throughput and easy-to-use allele-specific markers were developed and validated. The validated markers were converted into low-cost and high-throughput genotyping assays, also called KASP assays, for large-scale and easy deployment in breeding programs. To make this assay more valuable, KASP markers for high oleic acid content from fatty acid desaturase genes (FAD2B and FAD2A) have also been included in the genotyping assay for pyramiding high oleic acid content and foliar disease resistance. KASP assays were validated on diverse sets of resistant and susceptible genotypes and deployed successfully in developing marker-assisted backcrossed introgression lines [2, 33].

For high oleic acid, the high-oleic acid varieties Girnar 4 and Girnar 5 were sequenced at 10X coverage using an Illumina HiSeq2500 platform. Subsequently, the obtained sequences were aligned to the reference genome to identify non-synonymous SNPs in the FAD2A and FAD2B genes. Next, we developed allele-specific genotyping assays to accurately select both mutant alleles.

Development and validation of a high-throughput genotyping assay for foliar fungal disease (rust and late leaf spot) resistance and high oleic acid content

SNPs from QTL regions, specifically from coding regions and affecting the function of important disease resistance-related candidate genes, were used for the development of KASP assays. These SNPs were selected from the regions 1 kb upstream or downstream or ideally within candidate genes present in the QTL regions on chromosomes A02 (LLS resistance) and A03 (rust and LLS resistance) using mapping information from GBS analysis and B03 (rust and LLS resistance) using information from QTL-seq analysis. Additionally, SNPs from genic regions of the fatty acid desaturase (FAD) gene from both sub-genomes, A. ipaensis (fad2a) on chromosome Ah09 and A. duranensis (fad2b) on Ah19, were targeted to design KASP markers for high oleic acid content. The SNPs were converted into KASP markers by using 50 bp upstream and 50 bp downstream sequences to develop user-friendly and cost-effective markers. KASP assays for each SNP marker were designed with two allele-specific forward primers and a common reverse primer. All the KASP primers used in this study are listed in Supplementary Table 2. Following development, the KASP markers were further validated on a diverse validation panel comprising both resistant and susceptible genotypes for rust and LLS as well as low- and high-oleic acid breeding lines.

Results

Genomic regions for leaf rust and LLS resistance

A total of 50 Gb of clean read data were generated for 217 RILs and two parental genotypes (TAG 24 and GPBD 4) of the mapping population. Additionally, resistant and susceptible bulks were constructed by mixing pooled DNA from 25 RILs exhibiting extreme phenotypes, i.e. resistant and susceptible for leaf rust (mean disease score of 3.7 for the resistant bulk and 7.7 for the susceptible bulk) and LLS (mean disease score of 4.4 for the resistant bulk and 8.1 for the susceptible bulk, as explained in [24]). The susceptible parent, TAG 24, had mean disease scores of 7.5 and 8.4 for rust and LLS, respectively, while the resistant parent, GPBD 4, had scores of 3.0 and 3.7, respectively.

A total of 53,103 single nucleotide polymorphisms (SNPs) were identified across the RIL population using genotyping-by-sequencing. After filtration, 1,529 SNPs were found to be polymorphic between the parents. The chi-square (χ2) values calculated for each SNP marker were used to determine the goodness of fit to the expected 1:1 segregation ratio; highly distorted and unlinked markers were filtered out and not considered for linkage map construction. A dense genetic map with 1,119 SNP loci, including 831 SNPs and 288 simple sequence repeats (SSRs), was constructed with a map density of 1.48 cM/loci and a map length of 1,660.5 cM (Fig. 1A; Table 1).

Table 1 Major QTLs identified for leaf rust and late leaf spot resistance in groundnuts using the RIL population TAG 24 × GPBD 4
Fig. 1
figure 1

Genetic map and major QTLs identified for leaf rust and late leaf spot resistance. (A) Genetic map constructed using genotyping-by-sequencing data generated on the RIL population TAG 24 × GPBD 4. (B) Consistent major effect QTL identified on chromosome A02. (C) Consistent major effect QTL identified on chromosome A03

Genome-wide QTL discovery revealed 10 main effect QTLs, with LODs ranging from 3.5 to 71.6 and PVEs ranging from 5.1 to 87.1%. Of these 10 QTLs, three were associated with leaf rust resistance, and seven were associated with LLS resistance. For rust resistance, 3 major main effect QTLs were detected, with LODs ranging from 6.7 to 56.8 and PVEs ranging from 15.9 to 80.9%. Similarly, 7 main effect QTLs were identified for LLS resistance, with LODs ranging from 3.5 to 71.6 and PVEs ranging from 5.6 to 87.1%. Furthermore, all three QTLs were identified as major main effect QTLs for leaf rust. However, of the 7 QTLs for LLS resistance, 5 had major effects (Table 1).

Of the 3 major main effect QTLs detected for leaf rust resistance, qRust-A02.1 was mapped at 1.2 cM (S2_464902 - S2_768713) on A02, with a LOD value of 11.3 and a PVE of 30.6%, showing a negative additive effect, indicating that a high level of resistance was contributed by GPBD 4. Similarly, a major main effect QTL, qRust-A03.1, was mapped on chromosome A03 at 87.8 cM (S3_133697495-S3_134897269), with LODs of 6.7 and 15.9 PVE%. The third major QTL for rust resistance, qRust-B02.1 (S12_1692764-S12_1687343), was mapped at 10.5 cM with an LOD of 56.8 and 80.9% PVE, and a negative additive effect indicated a high level of the rust-resistant segment from the donor GPBD 4. Although the major QTL detected in B02 for rust resistance was not stable, its high PVE and LOD indicate the necessity for further investigation using a larger genetic population.

In the case of late leaf spot resistance, of the 5 major main effect QTLs detected for LLS resistance, 4 QTLs were mapped to the 5 cM genomic region at the start of chromosome A02. Of these four major main effect QTLs mapped on chromosome A02, qLLS-A02.1 was mapped at 0.7 cM (S2_880921 - S2_99621) with (14.6–18.3 LOD and 23.3–30.1% PVE), qLLS-A02.2 was mapped at 2.0 cM (S2_976663 - S2_454537) with (11.7–30.0 LOD to 16.3% PVE), qLLS-A02.3 was mapped at 5.7 cM (S2_1228446 - S2_1281741) with (9.3–14.0 LOD and 17.3–23.9% PVE), and qLLS-A02.4 was mapped at 0.1 cM (S2_1240859 - S2_1209628) with (8.7–12.4 LOD and 14.7–18.2% PVE). One main effect QTL on chromosome A02 at 40.5 cM (S2_2024918-S2_4433397) was consistent across the four environments for LLS (Fig. 1B; Table 1). A major main effect QTL, qLLS-A03.1, on chromosome A03 was mapped at 87.8 cM (S3_133697495-S3_134897269), with 13.1–33.1 LOD and 25.4–69.0% PVE. The major main effect QTLs located on A02 spanning a physical map distance of 2.08 Mb (0.7–2.8 Mb) were the major QTLs for both rust and LLS, with up to 30.6 and 44.2% PVE, respectively. However, a common genomic region was identified on chromosome A03 for both rust and LLS resistance, with up to 87.1% and 69.0% PVE, respectively (Fig. 1C).

Initially, QTL-seq analyses based on WGRS data for six samples were performed using available genome assemblies for diploid progenitors [24]. Later, the analyses were repeated using a cultivated tetraploid genome assembly for the subspecies fastigiata [22]. The reference-guided assembly for the resistant parent GPBD 4 was developed with diploid progenitors and a cultivated tetraploid genome. We observed that the mapping percentage of the tetraploid reference genome was 10% greater than that of the diploid progenitor genome in the resistant parent, GPBD 4 (85.07% and 96.58%), Rust_Rbulk (85.24% and 96.68%), Rust_Sbulk (85.34% and 96.77%), LLS_Rbulk (85.01% and 96.34%) and LLS_Sbulk (84.97% and 96.78%). However, there was no significant difference in average mapping depth between the diploid and tetraploid reference genomes, as shown in Supplementary Table 1.

The initial QTL-Seq analysis of the diploid genome-based GPBD 4 assembly revealed a significant colocalized genomic region on chromosome A03 for leaf rust (131.60-134.66 Mb) and LLS resistance (131.67-134.65 Mb) [24]. After analyzing the resistant and susceptible bulks from the QTL region identified in A03, a total of 66 SNPs were effective SNPs identified at a read depth of ≥ 7 for leaf rust and LLS resistance, respectively (Supplementary Table 1). Similarly, QTL-seq analysis of the cultivated tetraploid genome-based GPBD 4 assembly also revealed a colocalised genomic region on the pseudomolecule Chr13 for rust (140.405–144.882 Mb) and LLS (140.808–144.705 Mb) resistance [22]. A total of 3,270 SNPs and 1,620 effective SNPs were identified at a read depth of ≥ 7 for leaf rust and LLS resistance, respectively, after analysing the resistant and susceptible bulks from the QTL region identified on Chr13 (Supplementary Table 1). It is important to note that the source for leaf rust and LLS resistance in GPBD 4 was traced back to A. cardenasii while performing analysis with the diploid genome; in contrast, mapping to Chr13 with the tetraploid genome indicated possible translocation of the genomic region controlling resistance from Chr03 to Chr13 after tetraploidization or due to high sequence similarity between the two sub-genomes. This colocalized genomic region yielded more effective SNPs for both diseases in the tetraploid reference genome than in the diploid reference genome.

Candidate genes for resistance to leaf rust and LLS

In GBS-based QTL analysis, major main effect QTLs for leaf rust and late leaf spot were identified on chromosomes A02 and A03. The genomic region of length 2.08 Mb (0.7–2.8 Mb) on chromosome A02 comprised 6 adjacent QTL regions (1 for rust resistance and 5 for LLS resistance) detected in various environments (Fig. 2a). Within this region, a total of 245 genes were located (www.peanutbase.org). Of these, 38 candidate genes were found to be affected by 72 GBS-SNPs. The 38 candidate genes included important disease resistance genes, such as LRR and NB-ARC domain disease resistance protein, calmodulin-binding heat-shock protein (Aradu. MS73B), DEAD-box ATP-dependent RNA helicase (Aradu. JB556), disease resistance LRR family protein (Aradu. KFQ7S), ubiquitin-protein ligase (Aradu. MMY8Q) and a few unknown proteins (Fig. 2b).

Fig. 2
figure 2

Fine mapping, gene discovery (LLSR1 gene for LLS resistance), marker development and validation of major main-effect QTLs controlling LLS resistance. (a) The 2.08 Mb genomic region on chromosome A02 wherein QTLs were detected across different seasons, with eight candidate resistance genes targeted for marker development. The gene IDs for these candidate resistance genes and KASP assay information are also provided. (b) Scatter plot showing the validation results for selected KASP assays (snpAH0003, snpAH0004, snpAH0005, snpAH0003, snpAH0010, snpAH0011, snpAH0014, snpAH0120, and snpAH0121). There are three clusters on each scatter plot. The red cluster indicates the resistant lines, the blue cluster indicates the susceptible lines, and the green clusters are the heterozygotes derived from the crosses between the resistant and susceptible lines

Interestingly, multiple SNPs were identified in the genomic regions of a few candidate genes. For instance, 3 SNPs were identified in the region of the disease resistance gene LRR and NB-ARC domain disease resistance protein (Aradu.0G2IC), one SNP (snpAH0014) in the disease resistance LRR family protein (Aradu. KFQ7S) (Fig. 2b).

In GBS-based genetic mapping, a common genomic region of 1.19 Mb (133.7-134.9 Mb) on A03 was identified, which overlapped with the regions identified via QTL-seq on chromosome A03 and Chr13 in diploid and tetraploid genomes, respectively. In the QTL-seq analysis, 3,139 SNPs, including 30 nonsynonymous SNPs affecting 25 candidate genes related to plant growth and defense, were detected in the genomic region of 3.06 Mb (131.60-134.66 Mb) on A03 identified for rust resistance using the diploid genome (Fig. 3a).

Fig. 3
figure 3

Fine mapping, gene discovery (LLSR1 gene for LLS resistance and LR1 gene for LR resistance), marker development and validation for major main-effect colocalised QTLs for rust (qRust-A03.0) and LLS (qLLS-A03.0) resistance. Figure (a) shows the 1.19 Mb genomic region on chromosome A03, wherein QTLs were detected across different seasons and 13 candidate resistance genes were targeted for marker development. The gene IDs for these candidate resistance genes and KASP assay information are also provided. (b) Scatter plot showing the validation results for selected KASP assays (snpAH0017, snpAH0020, snpAH0022, snpAH0024, snpAH0026, snpAH0135, snpAH0136, snpAH0155, snpAH0130, snpAH0139, snpAH0156, and snpAH0158). There are three clusters on each scatter plot. The red cluster indicates the resistant lines, the blue cluster indicates the susceptible lines, and the green clusters are the heterozygotes derived from the crosses between the resistant and susceptible lines

Similarly, although the genomic region of 2.98 Mb (131.67-134.65 Mb) on A03 was identified for rust resistance using a diploid genome that had 66 SNPs, none were nonsynonymous. Upon repeating the QTL-seq analysis using the tetraploid reference genome, this colocalized genomic region possessed 216 and 171 candidate genes for leaf rust and LLS resistance, respectively. The colocalized genomic region on pseudomolecule Chr13 for rust (140.405–144.882 Mb) and LLS (140.808–144.705 Mb) resistance had 121 and 51 nonsynonymous SNPs, respectively. The number of candidate genes detected through analysis of the tetraploid reference genome was much greater than that detected through analysis of the diploid reference genome. These genes included TIR-NBS-LRR, PPR proteins (AH13G54010), C2H2-like zinc finger protein (Aradu. AB2YQ), the disease resistance protein TIR-NBS-LRR class (Aradu. Z87JB, AH13G54010), pentatricopeptide repeats (AH13G53720), leucine-rich repeat kinase family protein (AH13G53800), serine/threonine-protein kinase (AH13G52650), RNA 3-terminal phosphate cyclase-like protein (Aradu. HH10J), C2H2-like zinc finger protein (Aradu. AB2YQ), and ubiquitin-associated (UBA)/TS-N domain-containing (AH13G54110) (Fig. 3b).

Development, validation and deployment of genotyping assays for foliar fungal diseases and high oleic acid

From the genomic regions identified via GBS-based genetic mapping and QTL-seq analysis [13, 22], 120 primer pairs were designed and synthesized. Among these primer pairs, 90 showed precise amplification, and 80 exhibited polymorphisms between contrasting genotypes for foliar disease resistance and oleic acid content. Further validation was performed on a panel of diverse genotypes comprising both susceptible and resistant genotypes. The susceptible genotypes included GJ 9, GJ 20, GJGHPS 1, SunOleic 95R, ICGV 07368, ICGV 06420, TMV 2, DH 86, TAG 24, TG 26, ICGV 91114 and JL 24, and the resistant genotypes included GPBD 4 and 11 introgression lines in the genetic background of ICGV 91114, JL 24 and TAG 24 developed through a marker-assisted backcrossing (MABC) approach.

In parallel, the KASP assay was designed for 49 SNPs, for which primer validation was successful among the parental genotypes of the RIL population (Table 2; Supplementary Table 2). These KASP markers were validated on a panel of 96 genotypes, which included introgression lines carrying resistance from GPBD 4. Of the KASP markers tested, the eight best performing KASP markers were finalized for further use in breeding. In addition to rust and LLS, we also developed a KASP assay in which the FAD2A and FAD2B mutant alleles regulate the oleic acid trait (Fig. 4).

Table 2 Major QTLs targeted for candidate gene discovery and KASP assay development for late leaf spot (LLS), leaf rust and high oleic acid content
Fig. 4
figure 4

Development and validation of KASP markers for fatty acid desaturase genes controlling high oleic acid. The figure shows the genomic position of FAD2A on chromosome 9 and FAD2B on chromosome 19 as well as the two validated KASP assays, namely, snpAH0116 (FAD2A) and snpAH0002 (FAD2B), for high oleic acid. There are three clusters on each scatter plot. The red cluster indicates the resistant lines, the blue cluster indicates the susceptible lines, and the green clusters are the heterozygotes derived from the crosses between the resistant and susceptible lines

The KASP marker successfully distinguished the FAD2B mutant allele on the HTPG platform, while the KASP marker for the FAD2A mutant allele failed to perform as expected. Therefore, we sequenced two newly developed high-oleic acid varieties (Girnar 4 and Girnar 5) and completed primer development and validation after sequence analysis of these FAD genes from Girnar 4 and Girnar 5. Finally, we developed a high-throughput genotyping assay platform that is easily accessible through genotyping platforms in India, Sweden and Australia. The transition from genotyping with low-throughput allele-specific markers to high-throughput KASP genotyping has significantly reduced the cost of selection for all three traits to a mere 2.5 USD from 13 USD, bringing more precision and time savings in performing selection and decision making. Encouragingly, KASP assays are now being used regularly in several breeding programs across different countries, including India, Malawi, Mali, Burkina Faso, Nigeria and Uganda.

Discussion

Resistance to foliar fungal diseases such as rust and late leaf spot are the preferred traits of farmers because of their devastating nature, while a high oleic acid content is an industry- and consumer-preferred trait due to increased shelf life and health benefits. Current groundnut breeding programs have identified resistance to foliar fungal diseases (LLS & rust) and high oleic acid as must have traits. Therefore, in this study, we combined KASP markers for foliar diseases and a high oleic acid content. Leaf rust and late leaf spot are two devastating foliar fungal diseases of groundnuts that cause significant yield loss in all groundnut-growing areas worldwide [23, 34]. Potential resistant sources derived from Arachis cardenasii, such as GPBD4 in India and SPT06-06 (GP-NC WS16) in the USA, have been extensively used to improve late leaf spot resistance through conventional breeding approaches. To further accelerate conventional breeding, the development of genetic markers for the introgression of late leaf spot resistance using marker-assisted backcrossing (MABC) has become a focus of groundnut breeding programs. Linked simple sequence repeat (SSR) markers were successfully used for foreground selection in groundnuts [2, 5, 33, 35]. However, genotyping with SSR markers is laborious and time consuming, limiting its utility in molecular breeding. Advances in sequencing and genotyping technologies resulted in the discovery of new genotyping platforms. For instance, diagnostic KASP markers have proven to be a rapid genotyping platform for foreground selection, thereby accelerating marker-assisted introgressions. In this study, we reported the development and scaling up of newly developed and validated high-quality KASP markers for high oleic acid content, leaf rust and late leaf spot.

In groundnut, QTL-seq has been successfully used for the discovery of genomic regions and candidate genes associated with LLS and rust resistance. To identify the genomic regions for LLS and rust resistance, diploid reference genomes (A. ipaensis and A. duranensis) were initially used [13], and the genomic region on chromosome A03 was identified; subsequently, cultivated tetraploid reference genomes (Shitoqui) were used [22], and the genomic region on chromosome B03 (chr13) in the tetraploid genome was identified. Furthermore, to confirm the genomic regions identified in QTL-Seq studies, we used a genotyping-by-sequencing (GBS)-based genetic mapping approach to identify the QTLs linked with late leaf spot and rust resistance. The overlapping genomic regions from QTL-Seq studies and the QTLs identified via GBS-based genetic mapping were used for the development of KASP markers. A total of 47 SNPs with a Δ-SNP index of -1 from both studies were used for the development of KASP markers for LLS and rust. Among these 47 SNPs, 39 SNPs were validated on a panel of resistant and susceptible lines for LLS and rust, and 2 KASPs for high oleic acid from the A- and B-genomes were validated successfully on contrasting oleic acid-containing lines.

All the targeted SNPs were selected from the genomic regions of potential disease resistance genes reported earlier in various crop plants. According to the pedigrees of the resistance donor and reference genomes used for analysis, the resistance QTLs for rust and LLS resistance have been reported either on chromosomes A02 or B02 or on chromosomes A03 or B03 of both sub-genomes in various studies. For instance, a common region was identified for rust and LLS resistance on chromosome A03 by [24] using a QTL-seq approach with diploid reference genomes. Using the same dataset, QTL-Seq analysis of the tetraploid genome was used to identify the region on chromosome chr13 [22]. There are also reports on the identification of resistant QTLs on chromosome A05 in the Tifrunner genetic background using the QTL-Seq approach [18]. However, in this study, potential markers for rust and LLS were developed and successfully validated for chromosome A03 [24] and chr13 (B03; [22]). Interestingly, in this study, which used GBS-based genetic mapping, an additional major QTL (qRust-B02.1) identified on chromosome B02 for rust resistance with a PVE of 80.9% is also a promising genomic region because of its high phenotypic variance. In the future, this region can be further explored using fine mapping to identify potential resistance alleles for rust and LLS resistance.

Fatty acid desaturase enzyme catalyzes desaturation of oleic to linoleic acid, and it has two copies in both sub-genomes of cultivated groundnut. Here, we developed two markers for high oleic acid content (snpAH0116 from the A sub-genome and snpAH002 from the B sub-genome) located in the genic region of fatty acid desaturase (FAD) from the A and B sub-genomes, respectively (Fig. 5a). Similarly, KASP markers for LLS (snpAH0010, snpAH0011, and snpAH00121) were located in the genic regions of the LRR and NB-ARC domain disease resistance proteins (Aradu.0G2IC). KASP snpAH0014 was also developed from a protein/LRR family protein in the disease resistance family (Aradu.KFQ7S). snpAH0022 is a disease resistance protein (TIR-NBS-LRR class) (Aradu.Z87JB). NBS-LRR proteins recognize effectors secreted by pathogens directly or indirectly that in turn activate downstream signaling pathways, leading to activation of the plant defense response [36]. The KASP marker snpAH0120 represents an important gene, SKP1-like protein (Aradu.D80L9). The SKP1 protein is a component of the SKP1 (S-phase kinase-associated protein 1) complex and regulates the ubiquitination of proteins targeted for proteasomal degradation. Ubiquitination is implicated in many cellular processes, including plant defense responses [37]. Moreover, it plays a role in plant hormone signaling [38] and the accumulation of nucleotide-binding leucine-rich repeat-type immune receptors [39]. Induced expression of SKP1 was observed in N. benthamiana after infection with potato virus X, which is an indication of the involvement of ubiquitination in plant defense [40], and the upregulation of Arabidopsis SKP1 (ASK1 and ASK2) was shown to be required for successful Agrobacterium-mediated plant transformation [41]. snpAH0004 is located in the genic region of oxoglutarate (2OG) and Fe(II)-dependent oxygenase superfamily protein Aradu.2U1V2. 2OGs are reportedly involved in melatonin metabolism and subsequently affect plant responses to cold, heat, salt, drought, heavy metal stress, and pathogen invasion [42]. SnpAH0017 was developed from the genic region of purple acid phosphatase (Aradu.H1HIG). It has been reported that an optimal level of PAP maintains resistance against Pseudomonas syringae. It has also been reported that a mutation in PAP results in susceptibility to Pseudomonas syringae [43]. SnpAH0018 was located in the genic region of the P-loop containing nucleoside triphosphate hydrolase (Aradu.14 × 1M). The P-loop-containing nucleotide triphosphate hydrolase superfamily contains several kinases that are known to be involved in disease resistance [44]. C2H2-like zinc finger proteins positively modulate the expression levels of stress-related genes by directly binding to TACAAT motifs in the promoter regions of pathogen-related genes such as ENHANCED DISEASE SUSCEPTIBILITY1, PHYTOALEXIN DEFICIENT4, and PATHOGENESIS-RELATED GENE1 (PR1, PR2, and PR5) [45]. KASP was developed from the transcription factor TCP2 (AH13G52920). TCPs regulate plant development and defense responses by stimulating the biosynthetic pathways of bioactive metabolites, such as brassinosteroid (BR), jasmonic acid (JA) and flavonoids [46]. SnpAH0026 was developed from the major intrinsic protein (MIP) family transporter Aradu.HH10J. MIP isoforms are transporters of hydrogen peroxide. H2O2 experiences cross talk with immune pathways, such as those involved in systemic acquired resistance (SAR) and pathogen-associated molecular pattern-triggered immunity (PTI), to regulate plant disease resistance [47]. Topless-related protein 1 (AH13G53060) was also targeted for the development of the KASP marker snpAH0151. TOPLESS (TPL)/TOPLESS-RELATED (TPR), TPL has been linked to plant stress responses through several phytohormones [48] and is also involved in plant responses to pathogen attack. Arabidopsis plants with depleted TPL/TPR activity were reported to be more susceptible to infection due to disruption of SNC1-mediated immune responses [49]. Additionally, TPL is recruited by TGA family bZIP transcription factors that control development and stress responses via ROXY adaptor proteins to repress TGA targets [50]. SnpAH0159 was targeted from the genic region of CBL-interacting serine/threonine-protein kinase (AH13G54550). Calcineurin B-like proteins (CBLs) are a family of calcium sensor proteins that interact with a group of serine/threonine kinases designated CBL-interacting protein kinases (CIPKs). CBL-CIPK complexes are crucially involved in relaying plant responses to many environmental signals and in regulating ion fluxes [51] (Table 3).

Table 3 Details of validated KASP assays for resistance to leaf rust and LLS and high oleic acid content tables
Fig. 5
figure 5

Lab-to-field applications and impact of KASP markers in groundnuts. (a) Illustration of the roadmap for developing diagnostic KASP markers; (b) Molecular breeding products using diagnostic markers for foliar fungal diseases and a high oleic acid content

These newly developed KASP markers were made publicly available (https://excellenceinbreeding.org/module3/kasp) before this publication, with genotyping options available in India, Sweden and Australia. To date, these markers have been successfully used in breeding programs in multiple countries, such as India, Australia, Kenya, Malawi, Uganda, Ghana, Senegal, Mali and the United States. The deployment of low-throughput allele-specific markers and these newly developed high-throughput KASP assays has resulted in the development and release of commercial cultivation of five high-oleic acid groundnut varieties, namely, Girnar 4 (ICGV 15083), Girnar 5 (ICGV 15090), GG 39 (ICGV 16697), GG 40 (ICGV 16668) and ICRC-1 (ICGV 16690), in India for cultivation in Gujarat, Rajasthan, Karnataka, Andhra Pradesh and Tamil Nadu. Low-throughput genotyping assays, such as allele-specific and simple sequence repeat (SSR) genotyping, are being used in developing two new foliar disease-resistant varieties, namely, Improved JL 24 (DBG 3) and Super TMV 2 (DBG 4), at the University of Agricultural Sciences (UAS), Dharwad, India. National Agricultural Research Stations (NARSs), such as the Indian Council of Agricultural Research-Directorate of Groundnut Research (ICAR-DGR), Junagadh, India [10]; Acharya N.G. Ranga Agricultural University, Andhra Pradesh, India; Professor Jayashankar Telangana State Agricultural University, Telangana, India [33]; and University of Agricultural Sciences Dharwad, India [6, 12], used allele-specific and KASP markers, leading to the development of improved groundnut varieties with enhanced LLS and rust resistance and oleic acid content. The newly developed diagnostic kit for LLS, rust and oleic acid has 10 KASP markers, including two KASPs for the selection of alleles with high oleic acid content (Fig. 3). These high-throughput genotyping assays are also being developed for other traits, such as fresh seed dormancy [52, 53] and seed weight [54, 55], and for identifying true allotetraploids derived from crosses of the A- and B/K-genomes of Arachis diploid species [56].

Conclusion

Although the contributions of homologous chromosomes (A02/B02 and A03/B03) from the two sub-genomes are debated, these validated markers provide opportunities for routine deployment in groundnut breeding programs. The reason for this confusion is mostly because of the high sequence similarity between the two sub-genomes, which poses serious challenges to existing software for mapping such associated genomic regions. Nevertheless, the authors still vouch that the contribution of resistance comes from the A sub-genome (A02 and A03) due to the resistance source being A. cardenasii, which belongs to the A genome. Most importantly, a diagnostic kit for LLS, rust and high oleic acid content will be an important genomic tool for improving foliar disease resistance and high oleic acid content in groundnuts. Moreover, the successful deployment of genomic resources underpinned a significant milestone in our research, effectively advancing the development and release of improved groundnut varieties (Fig. 5b). Selection in breeding progenies can be performed for four traits (resistance to rust and LLS, high oleic acid content and fresh seed dormancy) using only 10 SNPs, making large-scale deployment most affordable and cost-effective. Recently, we also optimized a protocol to perform single-seed-based genotyping in groundnuts and started processing a large number of samples in groundnut breeding programs for selection even before seeds went to the field [26], thereby saving resources and time. The integration of such technologies is highly important for reducing the duration of breeding cycles.

Data availability

The data presented in the study are deposited in NCBI, accession number PRJNA999952 and PRJNA992198.

Change history

Abbreviations

LLS:

Late Leaf Spot

MEGS:

Marker-Based Early Generation Advancement

QTL:

Quantitative trait locus

KASP:

Kompetitive allele-specific PCR

MABC:

Marker-assisted back-cross

SSR:

Simple sequence repeat

RIL:

Recombinant inbred line population

PVE:

Phenotypic variance explained

AICRP-G:

All India Coordinated Research Project on Groundnut

DAS:

Days after sowing

FAD:

Fatty acid desaturase

GBS:

Genotyping-by-sequencing

cM:

Centi-Morgan

SNP:

Single nucleotide polymorphism

References

  1. FAOSTAT. (2020), http://faostat3.fao.org/home/index.html.

  2. Shasidhar Y, Variath MT, Vishwakarma MK, et al. Improvement of three popular Indian groundnut varieties for foliar disease resistance and high oleic acid using SSR markers and SNP array in marker-assisted backcrossing. Crop J. 2020;8:1–15.

    Article  Google Scholar 

  3. Subrahmanyam P, Williams JH, McDonald D, Gibbons RW. The influence of foliar diseases and their control by selective fungicides on a range of groundnut (Arachis hypogaea L.) genotypes. Ann Appl Biol. 1984;104:467–76.

    Article  CAS  Google Scholar 

  4. Pandey MK, Pandey AK, Kumar R, et al. Translational genomics for achieving higher genetic gains in groundnut. Theor Appl Genet. 2020;133:1679–702.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Varshney RK, Pandey MK, Janila P, et al. Marker-assisted introgression of a QTL region to improve rust resistance in three elite and popular varieties of peanut (Arachis hypogaea L). Theor Appl Genet. 2014;127:1771–81.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Sharanabasappa BY, Bhat RS. (2016) Development of late leaf spot and rust resistant backcross lines in Jl 24 variety of groundnut (Arachis hypogaea L.).

  7. Kolekar RM, Sukruth M, Shirasawa K, et al. Marker-assisted backcrossing to develop foliar disease‐resistant genotypes in TMV 2 variety of peanut (Arachis hypogaea L). Plant Breed. 2017;136:948–53.

    Article  CAS  Google Scholar 

  8. Nawade B, Mishra GP, Radhakrishnan T, et al. Development of high oleic peanut lines through marker-assisted introgression of mutant ahFAD2 alleles and its fatty acid profiles under open-field and controlled conditions. 3 Biotech. 2019;9:1–16.

    Article  Google Scholar 

  9. Janila P, Pandey MK, Shasidhar Y, et al. Molecular breeding for introgression of fatty acid desaturase mutant alleles (ahFAD2A and ahFAD2B) enhances oil quality in high and low oil containing peanut genotypes. Plant Sci. 2016;242:203–13.

    Article  CAS  PubMed  Google Scholar 

  10. Bera SK, Kamdar JH, Kasundra SV, et al. Improving oil quality by altering levels of fatty acids through marker-assisted selection of ahfad2 alleles in peanut (Arachis hypogaea L). Euphytica. 2018;214:1–15.

    Article  CAS  Google Scholar 

  11. Huang B, Qi F, Sun Z et al. (2019) Marker-assisted backcrossing to improve seed oleic acid content in four elite and popular peanut (Arachis hypogaea L.) cultivars with high oil content. Breed Sci 18107.

  12. Jadhav MP, Gangurde SS, Hake AA et al. (2021) Genotyping-by-sequencing based genetic mapping identified major and consistent genomic regions for productivity and quality traits in peanut. Front Plant Sci 2034.

  13. Pandey MK, Roorkiwal M, Singh VK, et al. Emerging genomic tools for legume breeding: current status and future prospects. Front Plant Sci. 2016;7:455.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Varshney RK, Singh VK, Kumar A, et al. Can genomics deliver climate-change ready crops? Curr Opin Plant Biol. 2018;45:205–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Bertioli DJ, Cannon SB, Froenicke L, et al. The genome sequences of Arachis duranensis and Arachis Ipaensis, the diploid ancestors of cultivated peanut. Nat Genet. 2016;48:438–46.

    Article  CAS  PubMed  Google Scholar 

  16. Chen X, Li H, Pandey MK, et al. Draft genome of the peanut A-genome progenitor (Arachis duranensis) provides insights into geocarpy, oil biosynthesis, and allergens. Proc Natl Acad Sci. 2016;113:6785–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Lu Q, Li H, Hong Y, Zhang G, Wen S, Li X, Zhou G, Li S, Liu H, Liu H, Liu ZJ. Genome sequencing and analysis of the peanut B-genome progenitor (Arachis Ipaensis). Front Plant Sci. 2018;9:604.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Clevenger J, Chu Y, Chavarro C, et al. Genome-wide SNP genotyping resolves signatures of selection and tetrasomic recombination in peanut. Mol Plant. 2017;10:309–22.

    Article  CAS  PubMed  Google Scholar 

  19. Pandey MK, Agarwal G, Kale SM, et al. Development and evaluation of a high density genotyping ‘Axiom_Arachis’ array with 58 K SNPs for accelerating genetics and breeding in groundnut. Sci Rep. 2017a;7:1–10.

    Article  Google Scholar 

  20. Bertioli DJ, Jenkins J, Clevenger J, et al. The genome sequence of segmental allotetraploid peanut Arachis hypogaea. Nat Genet. 2019;51:877–84.

    Article  CAS  PubMed  Google Scholar 

  21. Chen X, Lu Q, Liu H, et al. Sequencing of cultivated peanut, Arachis hypogaea, yields insights into genome evolution and oil improvement. Mol Plant. 2019;12:920–34.

    Article  CAS  PubMed  Google Scholar 

  22. Zhuang W, Chen H, Yang M, et al. The genome of cultivated peanut provides insight into legume karyotypes, polyploid evolution and crop domestication. Nat Genet. 2019;51:865–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Sujay V, Gowda MVC, Pandey MK, et al. QTL analysis and construction of consensus genetic map for foliar diseases resistance based on two RIL populations in cultivated groundnut (Arachis hypogaea L). Mol Breed. 2012;32:773–88.

    Article  Google Scholar 

  24. Pandey MK, Khan AW, Singh VK, et al. QTL-seq approach identified genomic regions and diagnostic markers for rust and late leaf spot resistance in groundnut (Arachis hypogaea L). Plant Biotechnol J. 2017b;15:927–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Shirasawa K, Bhat RS, Khedikar YP, et al. Sequencing analysis of genetic loci for resistance for late leaf spot and rust in peanut (Arachis hypogaea L). Front Plant Sci. 2018;9:1727.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Parmar S, Deshmukh DB, Kumar R, et al. Single seed-based high-throughput genotyping and Rapid Generation Advancement for Accelerated Groundnut Genetics and Breeding Research. Agronomy. 2021;11:1226.

    Article  CAS  Google Scholar 

  27. Gowda MVC, Motagi BN, Naidu GK et al. (2002) GPBD 4: a Spanish bunch groundnut genotype resistant to rust and late leaf spot.

  28. Bradbury PJ, Zhang Z, Kroon DE, et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23:2633–5.

    Article  CAS  PubMed  Google Scholar 

  29. Van Ooijen JW. (2006) JoinMap® 4, Software for the calculation of genetic linkage maps in experimental populations. Kyazma BV, Wageningen 33.

  30. Voorrips R. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93:77–8.

    Article  CAS  PubMed  Google Scholar 

  31. Kosambi DD. (1944) The estimation of map distances from recombination values. DD Kosambi Sel Work Math Stat 125–30.

  32. Wang S, Basten CJ, Zeng ZB. (2007) Windows QTL cartographer 2.5. Department of Statistics, North Carolina State University, Raleigh, NC. http://statgen.ncsu.edu/qtlcart/WQTLCart.htm (accessed 12 Jan. 2023).

  33. Deshmukh DB, Marathi B, Sudini HK, et al. Combining high oleic acid trait and resistance to late leaf spot and rust diseases in groundnut (Arachis hypogaea L). Front Genet. 2020;11:514.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Gangurde SS, Kumar R, Pandey AK et al. (2019) Climate-smart groundnuts for achieving high productivity and improved quality: current status, challenges, and opportunities. Genomic Des Clim Oilseed Crop 133–72.

  35. Chu Y, Clevenger JP, Holbrook CC, et al. Registration of two peanut recombinant inbred lines (TifGP-5 and TifGP‐6) resistant to late leaf spot disease. J Plant Regist. 2022;16:635–40.

    Article  Google Scholar 

  36. DeYoung BJ, Innes RW. Plant NBS-LRR proteins in pathogen sensing and host defense. Nat Immunol. 2006;7:1243–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Furlan G, Klinkenberg J, Trujillo M. Regulation of plant immune receptors by ubiquitination. Front Plant Sci. 2012;3:238.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Santner A, Estelle M. Recent advances and emerging trends in plant hormone signaling. Nature. 2009;459:1071–8.

    Article  CAS  PubMed  Google Scholar 

  39. Cheng YT, Li Y, Huang S, et al. Stability of plant immune-receptor resistance proteins is controlled by SKP1-CULLIN1-F-box (SCF)-mediated protein degradation. Proc Natl Acad Sci. 2011;108:14694–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Ye C, Chen S, Payton M, et al. TGB p3 triggers the unfolded protein response and SKP 1-dependent programmed cell death. Mol Plant Pathol. 2013;14:241–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Anand A, Rojas CM, Tang Y, Mysore KS. Several components of SKP1/Cullin/F-box E3 ubiquitin ligase complex and associated factors play a role in Agrobacterium‐mediated plant transformation. New Phytol. 2012;195:203–16.

    Article  CAS  PubMed  Google Scholar 

  42. Wei S, Zhang W, Fu R, Zhang Y. Genome-wide characterization of 2-oxoglutarate and fe (II)-dependent dioxygenase family genes in tomato during growth cycle and their roles in metabolism. BMC Genomics. 2021;22:1–14.

    Article  Google Scholar 

  43. Ravichandran S, Stone SL, Benkel B, et al. Optimal level of purple acid phosphatase5 is required for maintaining complete resistance to Pseudomonas syringae. Front Plant Sci. 2015;6:568.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Arya P, Acharya V. Plant STAND P-loop NTPases: a current perspective of genome distribution, evolution, and function: plant STAND P-loop NTPases: genomic organization, evolution, and molecular mechanism models contribute broadly to plant pathogen defense. Mol Genet Genomics. 2018;293:17–31.

    Article  CAS  PubMed  Google Scholar 

  45. Shi H, Chan Z. The cysteine2/histidine2-type transcription factor ZINC FINGER OF ARABIDOPSIS THALIANA 6‐activated C‐REPEAT‐BINDING FACTOR pathway is essential for melatonin‐mediated freezing stress resistance in Arabidopsis. J Pineal Res. 2014;57:185–91.

    Article  CAS  PubMed  Google Scholar 

  46. Li S. The Arabidopsis thaliana TCP transcription factors: a broadening horizon beyond development. Plant Signal Behav. 2015;10:e1044192.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Bienert GP, Chaumont F. Aquaporin-facilitated transmembrane diffusion of hydrogen peroxide. Biochim Biophys Acta (BBA)-General Subj. 2014;1840:1596–604.

    Article  CAS  Google Scholar 

  48. Causier B, Lloyd J, Stevens L, Davies B. TOPLESS corepressor interactions and their evolutionary conservation in plants. Plant Signal Behav. 2012;7:325–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Zhu Z, Xu F, Zhang Y, et al. Arabidopsis resistance protein SNC1 activates immune responses through association with a transcriptional corepressor. Proc Natl Acad Sci. 2010;107:13960–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Uhrig JF, Huang L-J, Barghahn S, et al. CC-type glutaredoxins recruit the transcriptional corepressor TOPLESS to TGA-dependent target promoters in Arabidopsis thaliana. Biochim Biophys Acta (BBA)-Gene Regul Mech. 2017;1860:218–26.

    Article  CAS  Google Scholar 

  51. Hashimoto K, Eckert C, Anschütz U, et al. Phosphorylation of calcineurin B-like (CBL) calcium sensor proteins by their CBL-interacting protein kinases (CIPKs) is required for full activity of CBL-CIPK complexes toward their target proteins. J Biol Chem. 2012;287:7956–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Kumar R, Janila P, Vishwakarma MK, et al. Whole-genome resequencing‐based QTL‐seq identified candidate genes and molecular markers for fresh seed dormancy in groundnut. Plant Biotechnol J. 2020;18:992–1003.

    Article  CAS  PubMed  Google Scholar 

  53. Bomireddy D, Gangurde SS, Variath MT, et al. Discovery of major quantitative trait loci and candidate genes for fresh seed dormancy in groundnut. Agronomy. 2022;12:404.

    Article  CAS  Google Scholar 

  54. Gangurde SS, Khan AW, Janila P et al. (2022) Whole-genome sequencing based discovery of candidate genes and diagnostic markers for seed weight in groundnut. Plant Genome e20265.

  55. Gangurde SS, Pasupuleti J, Parmar S, et al. Genetic mapping identifies genomic regions and candidate genes for seed weight and shelling percentage in groundnut. Front Genet. 2023;14:1128182.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Levinson CM, Bertioli D, Chu Y, et al. Development and applications of KASP markers distinguishing A-and B/K-genomes of Arachis. Euphytica. 2021;217:1–16.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to acknowledge MARS-Wrigley, USA; the Government of Karnataka (GoK), India; the Bill & Melinda Gates Foundation (BMGF), USA; and the Indian Council of Agricultural Research (ICAR), India, for financial support for this research. We are grateful to Richard Whitehead for editing the language of the manuscript.

Funding

Financial assistance was received for this work from the ICAR-ICRISAT collaborative project, India; MARS-Wrigley, USA; the World Bank-assisted Karnataka Watershed Development Project-II (KWDP-II) funded by the Government of Karnataka (GoK), India; and the Bill & Melinda Gates Foundation (BMGF), USA.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization; data curation; funding acquisition; investigation; methodology; project administration; resources; software; supervision; validation; writing-original draft; writing-review & editing: Manish K. Pandey; Investigation; methodology; validation; writing original draft; writing-review & editing: Sunil S. Gangurde; Data curation; investigation; methodology; formal analysis; writing review & editing: Yaduru Shasidhar; Data curation; investigation; validation; methodology; writing-review & editing: Vinay Sharma; Data curation; validation; formal analysis: Sandeep Kale; Data curation; formal analysis: Aamir W. Khan; Investigation; validation; writing-review & editing: Priya Shah; Investigation; validation; writing-review & editing: Pushpesh Joshi; Investigation; validation; writing-review & editing: Ramesh Bhat; Resources; validation; writing-review & editing: Pasupuleti Janila; Resources; validation; writing-review & editing: Sandip K. Bera; Conceptualization; funding acquisition; methodology; resources; writing review editing: Rajeev K. Varshney.

Corresponding authors

Correspondence to Manish K. Pandey or Rajeev K. Varshney.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

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.

The original online version of this article was revised: errors were spotted in Introduction and Discussion sections of the article.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

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

Pandey, M.K., Gangurde, S.S., Shasidhar, Y. et al. High-throughput diagnostic markers for foliar fungal disease resistance and high oleic acid content in groundnut. BMC Plant Biol 24, 262 (2024). https://doi.org/10.1186/s12870-024-04987-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-04987-9

Keywords