Dissection of the Genetic Basis of Yield-Related Traits in the Chinese Peanut Mini-Core Collection Through Genome-Wide Association Studies

Fine mapping and candidate gene analysis of major QTLs for number of seeds per pod in Arachis hypogaea L

BACKGROUND

Peanut (Arachis hypogaea L., 2n = 2x = 20) is an important industrial and oil crop that is widely grown in more than 100 countries. In recent years, breeders have focused on increasing the seed number per pod to improve their yield in addition to other breeding for other key components of yield, including the pod number, seeds per pod, and 100-seed weight.

RESULTS

In this study, a secondary population of 1,114 BC1F2 lines was derived from a cross between the parents R45 and JNH3. Two stable major-effect quantitative trait loci of qRMPA09.1 and qRMPA09.2 were detected simultaneously and mapped within chromosomal intervals of approximately 400 Kb and 600 Kb on chromosome A09. Additionally, combined whole-genome and RNA sequencing analyses showed the differential expression of the Arahy.04JNDX gene that belongs to a MYB transcription factor (TF) between the two parents. The AhMYB51 gene was also inferred to influence the number of seeds per pod in peanuts. An examination of the backcross lines L2/L4 showed that AhMYB51 increases the rate of multiple pods per plant (RMSP) primarily by affecting brassinosteroids in the flowers, while its overexpression promotes the length of siliques in Arabidopsis thaliana.

CONCLUSIONS

Our findings provide valuable insights for the cloning of favorable alleles for RMSP in peanuts. The qRMSPA09.1 and qRMSPA09.2 are two novel QTL associated with the RMSP trait, with AhMYB51 predicted as its candidate gene. Moreover, the closely linked polymorphic SNP markers for loci of two significant QTLs may be useful in accelerating marker-assisted breeding in peanuts.

Integrated RNA-Seq and Metabolomics Analyses of Biological Processes and Metabolic Pathways Involved in Seed Development in Arachis hypogaea L

In peanut cultivation, fertility and seed development are essential for fruit quality and yield, while pod number per plant, seed number per pod, kernel weight, and seed size are indicators of peanut yield. In this study, metabolomic and RNA-seq analyses were conducted on the flowers and aerial pegs (aerpegs) of two peanut cultivars JNH3 (Jinonghei) and SLH (Silihong), respectively. Compared with SLH, JNH3 had 3840 up-regulated flower-specific differentially expressed genes (DEGs) and 5890 up-regulated aerpeg-specific DEGs. Compared with the JNH3 aerpegs, there were 4079 up-regulated variety-specific DEGs and 18 up-regulated differentially accumulated metabolites (DAMs) of JNH3 flowers, while there were 3732 up-regulated variety-specific DEGs and 48 up-regulated DAMs in SLH flowers. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the DEGs of JNH3 were associated with pollen germination and phenylalanine metabolism in flower and aerpeg tissues, respectively. In contrast, the DEGs of SLH were associated with protein degradation, amino acid metabolism, and DNA repair. However, there were significant differences in the lipids and lipid-like molecules between JNH3 flowers and SLH flowers. This investigation provides candidate genes and an experimental basis for the further improvement of high-quality and high-yield peanut varieties.

Using cross-country datasets for association mapping in Arachis hypogaea L

Groundnut (Arachis hypogaea L.) is one of the most important climate-resilient oil crops in sub-Saharan Africa. There is a significant yield gap for groundnut in Africa because of poor soil fertility, low agricultural inputs, biotic and abiotic stresses. Cross-country evaluations of promising breeding lines can facilitate the varietal development process. The objective of our study was to characterize popular test environments in Uganda (Serere and Nakabango) and Malawi (Chitala and Chitedze) and identify genotypes with stable superior yields for potential future release. Phenotypic data were generated for 192 breeding lines for yield-related traits, while genotypic data were generated using skim-sequencing. We observed significant variation (p < 0.001; p < 0.01; p < 0.05) across genotypes for all yield-related traits: days to flowering (DTF), pod yield (PY), shelling percentage, 100-seed weight, and grain yield within and across locations. Nakabango, Chitedze, and Serere were clustered as one mega-environment with the top five most stable genotypes being ICGV-SM 01709, ICGV-SM 15575, ICGV-SM 90704, ICGV-SM 15576, and ICGV-SM 03710, all Virginia types. Population structure analysis clustered the genotypes in three distinct groups based on market classes. Eight and four marker-trait associations (MTAs) were recorded for DTF and PY, respectively. One of the MTAs for DTF was co-localized within an uncharacterized protein on chromosome 13, while another one (TRv2Chr.11_3476885) was consistent across the two countries. Future studies will need to further characterize the candidate genes as well as confirm the stability of superior genotypes across seasons before recommending them for release.

Omics-driven advances in the understanding of regulatory landscape of peanut seed development

Peanuts (Arachis hypogaea) are an essential oilseed crop known for their unique developmental process, characterized by aerial flowering followed by subterranean fruit development. This crop is polyploid, consisting of A and B subgenomes, which complicates its genetic analysis. The advent and progression of omics technologies-encompassing genomics, transcriptomics, proteomics, epigenomics, and metabolomics-have significantly advanced our understanding of peanut biology, particularly in the context of seed development and the regulation of seed-associated traits. Following the completion of the peanut reference genome, research has utilized omics data to elucidate the quantitative trait loci (QTL) associated with seed weight, oil content, protein content, fatty acid composition, sucrose content, and seed coat color as well as the regulatory mechanisms governing seed development. This review aims to summarize the advancements in peanut seed development regulation and trait analysis based on reference genome-guided omics studies. It provides an overview of the significant progress made in understanding the molecular basis of peanut seed development, offering insights into the complex genetic and epigenetic mechanisms that influence key agronomic traits. These studies highlight the significance of omics data in profoundly elucidating the regulatory mechanisms of peanut seed development. Furthermore, they lay a foundational basis for future research on trait-related functional genes, highlighting the pivotal role of comprehensive genomic analysis in advancing our understanding of plant biology.

Genome-wide association study and development of molecular markers for yield and quality traits in peanut (Arachis hypogaea L.)

BACKGROUND

This study aims to decipher the genetic basis governing yield components and quality attributes of peanuts, a critical aspect for advancing molecular breeding techniques. Integrating genotype re-sequencing and phenotypic evaluations of seven yield components and two grain quality traits across four distinct environments allowed for the execution of a genome-wide association study (GWAS).

RESULTS

The nine phenotypic traits were all continuous and followed a normal distribution. The broad heritability ranged from 88.09 to 98.08%, and the genotype-environment interaction effects were all significant. There was a highly significant negative correlation between protein content (PC) and oil content (OC). The 10× genome re-sequencing of 199 peanut accessions yielded a total of 631,988 high-quality single nucleotide polymorphisms (SNPs), with 374 significant SNP loci identified in association with the nine traits of interest. Notably, 66 of these pertinent SNPs were detected in multiple environments, and 48 of them were linked to multiple traits of interest. Five loci situated on chromosome 16 (Chr16) exhibited pleiotropic effects on yield traits, accounting for 17.64-32.61% of the observed phenotypic variation. Two loci on Chr08 were found to be strongly associated with protein and oil contents, accounting for 12.86% and 14.06% of their respective phenotypic variations, respectively. Linkage disequilibrium (LD) block analysis of these seven loci unraveled five nonsynonymous variants, leading to the identification of one yield-related candidate gene and two quality-related candidate genes. The correlation between phenotypic variation and SNP loci in these candidate genes was validated by Kompetitive allele-specific PCR (KASP) marker analysis.

CONCLUSIONS

Overall, molecular markers were developed for genetic loci associated with yield and quality traits through a GWAS investigation of 199 peanut accessions across four distinct environments. These molecular tools can aid in the development of desirable peanut germplasm with an equilibrium of yield and quality through marker-assisted breeding.

Genome-wide association study for yield-related traits in faba bean (Vicia faba L.)

Yield is the most complex trait to improve crop production, and identifying the genetic determinants for high yield is a major issue in breeding new varieties. In faba bean (Vicia faba L.), quantitative trait loci (QTLs) have previously been detected in studies of biparental mapping populations, but the genes controlling the main trait components remain largely unknown. In this study, we investigated for the first time the genetic control of six faba bean yield-related traits: shattering (SH), pods per plant (PP), seeds per pod (SP), seeds per plant (SPL), 100-seed weight (HSW), and plot yield (PY), using a genome-wide association study (GWAS) on a worldwide collection of 352 homozygous faba bean accessions with the aim of identifying markers associated with them. Phenotyping was carried out in field trials at three locations (Spain, United Kingdom, and Serbia) over 2 years. The faba bean panel was genotyped with the Affymetrix faba bean SNP-chip yielding 22,867 SNP markers. The GWAS analysis identified 112 marker-trait associations (MTAs) in 97 candidate genes, distributed over the six faba bean chromosomes. Eight MTAs were detected in at least two environments, and five were associated with multiple traits. The next step will be to validate these candidates in different genetic backgrounds to provide resources for marker-assisted breeding of faba bean yield.

Designing future peanut: the power of genomics-assisted breeding

KEY MESSAGE

Integrating GAB methods with high-throughput phenotyping, genome editing, and speed breeding hold great potential in designing future smart peanut cultivars to meet market and food supply demands. Cultivated peanut (Arachis hypogaea L.), a legume crop greatly valued for its nourishing food, cooking oil, and fodder, is extensively grown worldwide. Despite decades of classical breeding efforts, the actual on-farm yield of peanut remains below its potential productivity due to the complicated interplay of genotype, environment, and management factors, as well as their intricate interactions. Integrating modern genomics tools into crop breeding is necessary to fast-track breeding efficiency and rapid progress. When combined with speed breeding methods, this integration can substantially accelerate the breeding process, leading to faster access of improved varieties to farmers. Availability of high-quality reference genomes for wild diploid progenitors and cultivated peanuts has accelerated the process of gene/quantitative locus discovery, developing markers and genotyping assays as well as a few molecular breeding products with improved resistance and oil quality. The use of new breeding tools, e.g., genomic selection, haplotype-based breeding, speed breeding, high-throughput phenotyping, and genome editing, is probable to boost genetic gains in peanut. Moreover, renewed attention to efficient selection and exploitation of targeted genetic resources is also needed to design high-quality and high-yielding peanut cultivars with main adaptation attributes. In this context, the combination of genomics-assisted breeding (GAB), genome editing, and speed breeding hold great potential in designing future improved peanut cultivars to meet market and food supply demands.

The groundnut improvement network for Africa (GINA) germplasm collection: a unique genetic resource for breeding and gene discovery

Cultivated peanut or groundnut (Arachis hypogaea L.) is a grain legume grown in many developing countries by smallholder farmers for food, feed, and/or income. The speciation of the cultivated species, that involved polyploidization followed by domestication, greatly reduced its variability at the DNA level. Mobilizing peanut diversity is a prerequisite for any breeding program for overcoming the main constraints that plague production and for increasing yield in farmer fields. In this study, the Groundnut Improvement Network for Africa assembled a collection of 1,049 peanut breeding lines, varieties, and landraces from 9 countries in Africa. The collection was genotyped with the Axiom_Arachis2 48K SNP array and 8,229 polymorphic single nucleotide polymorphism (SNP) markers were used to analyze the genetic structure of this collection and quantify the level of genetic diversity in each breeding program. A supervised model was developed using dapc to unambiguously assign 542, 35, and 172 genotypes to the Spanish, Valencia, and Virginia market types, respectively. Distance-based clustering of the collection showed a clear grouping structure according to subspecies and market types, with 73% of the genotypes classified as fastigiata and 27% as hypogaea subspecies. Using STRUCTURE, the global structuration was confirmed and showed that, at a minimum membership of 0.8, 76% of the varieties that were not assigned by dapc were actually admixed. This was particularly the case of most of the genotype of the Valencia subgroup that exhibited admixed genetic heritage. The results also showed that the geographic origin (i.e. East, Southern, and West Africa) did not strongly explain the genetic structure. The gene diversity managed by each breeding program, measured by the expected heterozygosity, ranged from 0.25 to 0.39, with the Niger breeding program having the lowest diversity mainly because only lines that belong to the fastigiata subspecies are used in this program. Finally, we developed a core collection composed of 300 accessions based on breeding traits and genetic diversity. This collection, which is composed of 205 genotypes of fastigiata subspecies (158 Spanish and 47 Valencia) and 95 genotypes of hypogaea subspecies (all Virginia), improves the genetic diversity of each individual breeding program and is, therefore, a unique resource for allele mining and breeding.

Identification of QTL for kernel weight and size and analysis of the pentatricopeptide repeat (PPR) gene family in cultivated peanut (Arachis hypogaea L.)

Peanut (Arachis hypogaea L.) is an important oilseed crop worldwide. Improving its yield is crucial for sustainable peanut production to meet increasing food and industrial requirements. Deciphering the genetic control underlying peanut kernel weight and size, which are essential components of peanut yield, would facilitate high-yield breeding. A high-density single nucleotide polymorphism (SNP)-based linkage map was constructed using a recombinant inbred lines (RIL) population derived from a cross between the variety Yuanza9102 and a germplasm accession wt09-0023. Kernel weight and size quantitative trait loci (QTLs) were co-localized to a 0.16 Mb interval on Arahy07 using inclusive composite interval mapping (ICIM). Analysis of SNP, and Insertion or Deletion (INDEL) markers in the QTL interval revealed a gene encoding a pentatricopeptide repeat (PPR) superfamily protein as a candidate closely linked with kernel weight and size in cultivated peanut. Examination of the PPR gene family indicated a high degree of collinearity of PPR genes between A. hypogaea and its diploid progenitors, Arachis duranensis and Arachis ipaensis. The candidate PPR gene, Arahy.JX1V6X, displayed a constitutive expression pattern in developing seeds. These findings lay a foundation for further fine mapping of QTLs related to kernel weight and size, as well as validation of candidate genes in cultivated peanut.

Dissection of the Genetic Basis of Resistance to Stem Rot in Cultivated Peanuts (Arachis hypogaea L.) through Genome-Wide Association Study

Peanut (Arachis hypogaea) is an important oilseed and cash crop worldwide, contributing an important source of edible oil and protein for human nutrition. However, the incidence of stem rot disease caused by Athelia rolfsii poses a major challenge to peanut cultivation, resulting in significant yield losses. In this study, a panel of 202 peanut accessions was evaluated for their resistance to stem rot by inoculating plants in the field with A. rolfsii-infested oat grains in three environments. The mean disease index value of each environment for accessions in subsp. fasitigiate and subsp. hypogaea showed no significant difference. Accessions from southern China displayed the lowest disease index value compared to those from other ecological regions. We used whole-genome resequencing to analyze the genotypes of the accessions and to identify significant SNPs associated with stem rot resistance through genome-wide association study (GWAS). A total of 121 significant SNPs associated with stem rot resistance in peanut were identified, with phenotypic variation explained (PVE) ranging from 12.23% to 15.51%. A total of 27 candidate genes within 100 kb upstream and downstream of 23 significant SNPs were annotated, which have functions related to recognition, signal transduction, and defense response. These significant SNPs and candidate genes provide valuable information for further validation and molecular breeding to improve stem rot resistance in peanut.

An Overview of Mapping Quantitative Trait Loci in Peanut (Arachis hypogaea L.)

Quantitative Trait Loci (QTL) mapping has been thoroughly used in peanut genetics and breeding in spite of the narrow genetic diversity and the segmental tetraploid nature of the cultivated species. QTL mapping is helpful for identifying the genomic regions that contribute to traits, for estimating the extent of variation and the genetic action (i.e., additive, dominant, or epistatic) underlying this variation, and for pinpointing genetic correlations between traits. The aim of this paper is to review the recently published studies on QTL mapping with a particular emphasis on mapping populations used as well as traits related to kernel quality. We found that several populations have been used for QTL mapping including interspecific populations developed from crosses between synthetic tetraploids and elite varieties. Those populations allowed the broadening of the genetic base of cultivated peanut and helped with the mapping of QTL and identifying beneficial wild alleles for economically important traits. Furthermore, only a few studies reported QTL related to kernel quality. The main quality traits for which QTL have been mapped include oil and protein content as well as fatty acid compositions. QTL for other agronomic traits have also been reported. Among the 1261 QTL reported in this review, and extracted from the most relevant studies on QTL mapping in peanut, 413 (~33%) were related to kernel quality showing the importance of quality in peanut genetics and breeding. Exploiting the QTL information could accelerate breeding to develop highly nutritious superior cultivars in the face of climate change.

Genetic mapping identified three hotspot genomic regions and candidate genes controlling heat tolerance-related traits in groundnut

Groundnut productivity and quality have been impeded by rising temperatures in semi-arid environments. Hence, understanding the effects and molecular mechanisms of heat stress tolerance will aid in tackling yield losses. In this context, a recombinant inbred line (RIL) population was developed and phenotyped for eight seasons at three locations for agronomic, phenological, and physiological traits under heat stress. A genetic map was constructed using genotyping-by-sequencing with 478 single-nucleotide polymorphism (SNP) loci spanning a map distance of 1,961.39 cM. Quantitative trait locus (QTL) analysis using phenotypic and genotypic data identified 45 major main-effect QTLs for 21 traits. Intriguingly, three QTL clusters (Cluster-1-Ah03, Cluster-2-Ah12, and Cluster-3-Ah20) harbor more than half of the major QTLs (30/45, 66.6%) for various heat tolerant traits, explaining 10.4%-38.6%, 10.6%-44.6%, and 10.1%-49.5% of phenotypic variance, respectively. Furthermore, important candidate genes encoding DHHC-type zinc finger family protein (arahy.J0Y6Y5), peptide transporter 1 (arahy.8ZMT0C), pentatricopeptide repeat-containing protein (arahy.4A4JE9), Ulp1 protease family (arahy.X568GS), Kelch repeat F-box protein (arahy.I7X4PC), FRIGIDA-like protein (arahy.0C3V8Z), and post-illumination chlorophyll fluorescence increase (arahy.92ZGJC) were the underlying three QTL clusters. The putative functions of these genes suggested their involvement in seed development, regulating plant architecture, yield, genesis and growth of plants, flowering time regulation, and photosynthesis. Our results could provide a platform for further fine mapping, gene discovery, and developing markers for genomics-assisted breeding to develop heat-tolerant groundnut varieties.

Genome-wide association studies reveal novel loci for resistance to groundnut rosette disease in the African core groundnut collection

KEY MESSAGE

We identified markers associated with GRD resistance after screening an Africa-wide core collection across three seasons in Uganda Groundnut is cultivated in several African countries where it is a major source of food, feed and income. One of the major constraints to groundnut production in Africa is groundnut rosette disease (GRD), which is caused by a complex of three agents: groundnut rosette assistor luteovirus, groundnut rosette umbravirus and its satellite RNA. Despite several years of breeding for GRD resistance, the genetics of the disease is not fully understood. The objective of the current study was to use the African core collection to establish the level of genetic variation in their response to GRD, and to map genomic regions responsible for the observed resistance. The African groundnut core genotypes were screened across two GRD hotspot locations in Uganda (Nakabango and Serere) for 3 seasons. The Area Under Disease Progress Curve combined with 7523 high quality SNPs were analyzed to establish marker-trait associations (MTAs). Genome-Wide Association Studies based on Enriched Compressed Mixed Linear Model detected 32 MTAs at Nakabango: 21 on chromosome A04, 10 on B04 and 1 on B08. Two of the significant markers were localised on the exons of a putative TIR-NBS-LRR disease resistance gene on chromosome A04. Our results suggest the likely involvement of major genes in the resistance to GRD but will need to be further validated with more comprehensive phenotypic and genotypic datasets. The markers identified in the current study will be developed into routine assays and validated for future genomics-assisted selection for GRD resistance in groundnut.

Genetic dissection of fatty acid components in the Chinese peanut (Arachis hypogaea L.) mini-core collection under multi-environments

Peanut (Arachis hypogaea L.) is an important source of edible oil and protein for human nutrition. The quality of peanut seed oil is mainly determined by the composition of fatty acids, especially the contents of oleic acid and linoleic acid. Improving the composition of fatty acids in the seed oil is one of the main objectives for peanut breeding globally. To uncover the genetic basis of fatty acids and broaden the genetic variation in future peanut breeding programs, this study used genome-wide association studies (GWAS) to identify loci associated with target traits and developed diagnostic marker. The contents of eight fatty acid components of the Chinese peanut mini-core collection were measured under four environments. Using the phenotypic information and over one hundred thousand single nucleotide polymorphisms (SNPs), GWAS were conducted to investigate the genetics basis of fatty acids under multi-environments. Overall, 75 SNPs were identified significant trait associations with fatty acid components. Nineteen associations were repeatedly identified in multiple environments, and 13 loci were co-associated with two or three traits. Three stable major associated loci were identified, including two loci for oleic acid and linoleic acid on chromosome A09 [mean phenotypic variation explained (PVE): 38.5%, 10.35%] and one for stearic acid on B06 (mean PVE: 23%). According to functional annotations, 21 putative candidate genes related to fatty acid biosynthesis were found underlying the three associations. The allelic effect of SNP A09-114690064 showed that the base variation was highly correlated with the phenotypic variation of oleic acid and linoleic acid contents, and a cost-effective Kompetitive allele-Specific PCR (KASP) diagnostic marker was developed. Furthermore, the SNP A09-114690064 was found to change the cis-element CAAT (-) in the promoter of ahFAD2A to YACT (+), leading dozens of times higher expression level. The enhancer-like activity of ahFAD2A promoter was identified that was valuable for enriching the regulation mechanism of ahFAD2A. This study improved our understanding on the genetic architecture of fatty acid components in peanut, and the new effective diagnostic marker would be useful for marker-assisted selection of high-oleic peanut breeding.

A first insight into the genetics of maturity trait in Runner × Virginia types peanut background

'Runner' and 'Virginia', the two main market types of Arachis hypogaea subspecies hypogaea, differ in several agricultural and industrial characteristics. One such trait is time to maturation (TTM), contributing to the specific environmental adaptability of each subspecies. However, little is known regarding TTM's genetic and molecular control in peanut in general, and particularly in the Runner/Virginia background. Here, a recombinant inbred line population, originating from a cross between an early-maturing Virginia and a late-maturing Runner type, was used to detect quantitative trait loci (QTL) for maturity. An Arachis SNP-array was used for genotyping, and a genetic map with 1425 SNP loci spanning 24 linkage groups was constructed. Six significant QTLs were identified for the maturity index (MI) trait on chromosomes A04, A08, B02 and B04. Two sets of stable QTLs in the same loci were identified, namely qMIA04a,b and qMIA08_2a,b with 11.5%, 8.1% and 7.3%, 8.2% of phenotypic variation explained respectively in two environments. Interestingly, one consistent QTL, qMIA04a,b, overlapped with the previously reported QTL in a Virginia × Virginia population having the same early-maturing parent ('Harari') in common. The information and materials generated here can promote informed targeting of peanut idiotypes by indirect marker-assisted selection.

High-Density Genetic Variation Map Reveals Key Candidate Loci and Genes Associated With Important Agronomic Traits in Peanut

Peanut is one of the most important cash crops with high quality oil, high protein content, and many other nutritional elements, and grown globally. Cultivated peanut (Arachis hypogaea L.) is allotetraploid with a narrow genetic base, and its genetics and molecular mechanisms controlling the agronomic traits are poorly understood. Here, we report a comprehensive genome variation map based on the genotyping of a panel of 178 peanut cultivars using Axiom_Arachis2 SNP array, including 163 representative varieties of different provinces in China, and 15 cultivars from 9 other countries. According to principal component analysis (PCA) and phylogenetic analysis, the peanut varieties were divided into 7 groups, notable genetic divergences between the different areas were shaped by environment and domestication. Using genome-wide association study (GWAS) analysis, we identified several marker-trait associations (MTAs) and candidate genes potentially involved in regulating several agronomic traits of peanut, including one MTA related with hundred seed weight, one MTA related with total number of branches, and 14 MTAs related with pod shape. This study outlines the genetic basis of these peanut cultivars and provides 13,125 polymorphic SNP markers for further distinguishing and utility of these elite cultivars. In addition, the candidate loci and genes provide valuable information for further fine mapping of QTLs and improving the quality and yield of peanut using a genomic-assisted breeding method.

Genotyping-by-Sequencing Based Genetic Mapping Identified Major and Consistent Genomic Regions for Productivity and Quality Traits in Peanut

With an objective of identifying the genomic regions for productivity and quality traits in peanut, a recombinant inbred line (RIL) population developed from an elite variety, TMV 2 and its ethyl methane sulfonate (EMS)-derived mutant was phenotyped over six seasons and genotyped with genotyping-by-sequencing (GBS), Arachis hypogaea transposable element (AhTE) and simple sequence repeats (SSR) markers. The genetic map with 700 markers spanning 2,438.1 cM was employed for quantitative trait loci (QTL) analysis which identified a total of 47 main-effect QTLs for the productivity and oil quality traits with the phenotypic variance explained (PVE) of 10-52% over the seasons. A common QTL region (46.7-50.1 cM) on Ah02 was identified for the multiple traits, such as a number of pods per plant (NPPP), pod weight per plant (PWPP), shelling percentage (SP), and test weight (TW). Similarly, a QTL (7.1-18.0 cM) on Ah16 was identified for both SP and protein content (PC). Epistatic QTL (epiQTL) analysis revealed intra- and inter-chromosomal interactions for the main-effect QTLs and other genomic regions governing these productivity traits. The markers identified by a single marker analysis (SMA) mapped to the QTL regions for most of the traits. Among the five potential candidate genes identified for PC, SP and oil quality, two genes (Arahy.7A57YA and Arahy.CH9B83) were affected by AhMITE1 transposition, and three genes (Arahy.J5SZ1I, Arahy.MZJT69, and Arahy.X7PJ8H) involved functional single nucleotide polymorphisms (SNPs). With major and consistent effects, the genomic regions, candidate genes, and the associated markers identified in this study would provide an opportunity for gene cloning and genomics-assisted breeding for increasing the productivity and enhancing the quality of peanut.