Next Article in Journal
Peanut Pickup Combine Harvester Seedling Vines Crushing Mechanism Performance Enhancement and Pilot Studies
Previous Article in Journal
Design and Optimization of Key Parameters for a Machine Vision-Based Walnut Shell–Kernel Separation Device
Previous Article in Special Issue
Analysis of BnGPAT9 Gene Expression Patterns in Brassica napus and Its Impact on Seed Oil Content
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Hotspot Regions for Candidate Genes Associated with Peanut (Arachis hypogaea L.) Pod and Seed Size on Chromosome A05

1
College of Agronomy, Shanxi Agricultural University, Taigu, Jinzhong 030801, China
2
Institute of Industrial Crops, Shanxi Agricultural University, Taiyuan 030031, China
3
College of Plant Protection, Shanxi Agricultural University, Taigu, Jinzhong 030801, China
*
Authors to whom correspondence should be addressed.
Agriculture 2024, 14(9), 1634; https://doi.org/10.3390/agriculture14091634
Submission received: 27 August 2024 / Revised: 10 September 2024 / Accepted: 13 September 2024 / Published: 18 September 2024
(This article belongs to the Special Issue Molecular Breeding and Genetic Improvement of Oilseed Crops)

Abstract

:
The size of peanut pods and seeds, which directly affects yield and quality, also has significant implications for mechanized production and market efficiency. Identifying relevant loci and mining candidate genes is crucial for cultivating high-yield peanut varieties. In this study, we employed advanced generation recombinant inbred lines developed by crossbreeding Huayu 44 and DF12 as the experimental material. Quantitative trait locus (QTL) mapping for traits related to pod and seed size was conducted across six environments. A total of 44 QTLs were detected, distributed on chromosomes A02, A05, B04, B08, and B10. An enrichment region for multiple QTLs was also identified on chromosome A05 (19.28~52.32 cm). In this region, 10 KASP markers were developed, narrowing the enrichment area to two candidate gene hotspot regions of 600.9 kb and 721.2 kb. By combining gene prediction and functional annotation within the intervals, 10 candidate genes, including those encoding cytochrome P450 protein, polyamine synthase, mannose-1-phosphate guanylyltransferase, pentatricopeptide repeat protein, and E2F transcription factor, were identified as regulators of pod and seed size. This study provides technical support for the genetic improvement and key gene identification of peanut pod and seed size.

1. Introduction

Peanut (Arachis hypogaea L.) is a globally significant oilseed and economic crop. Owing to its rich content of protein, fat, and other nutrients, it is widely used in food and oil production, making it an important source of plant-based protein for human consumption. Peanuts are primarily cultivated in the vast region between 40° S latitude and 40° N latitude, with major producing countries including China, India, Nigeria, the United States, and Indonesia. As of 2021, the global peanut cultivation area was 32.72 million hectares, with a total production of 53.93 million tons [1].
Cultivated peanuts (2n = 4x = 40) originated in South America. They evolved gradually through distant hybridization between two wild diploid species, A. duranensis (AA genome) and A. ipaensis (BB genome), and spontaneous chromosome doubling [2,3,4,5,6]. The seeds of the wild diploid species are very small, weighing only 0.1 to 0.3 g [7]. During the early domestication of peanuts, the size of the pods and seeds directly influenced yield and harvesting efficiency. People tend to select individuals with larger seeds and pods, increasing the size of the seeds and pods in cultivated peanuts. Larger pods not only provide more space and nutrients for the seeds, but also enhance their adaptability and resistance to adverse conditions. Therefore, they are considered advantageous for increasing yield [8]. In modern agriculture, larger pods facilitate mechanized harvesting, thereby improving efficiency. Larger seeds typically accumulate more nutrients, and during seed maturation, the accumulation of proteins and fats within the seeds increases. The pod-filling process is considered complete when the seeds accumulate sufficient nutrients and reach their maximum size [9]. Additionally, when the size of the pods and seeds reaches an appropriate ratio, shelling becomes easier and the shelling percentage is higher. Therefore, effectively regulating this trait can not only enhance the yield and quality of peanuts but also optimize processing and utilization, expanding the application range of peanut products; this maximizes the value of peanuts as an economically significant crop.
Peanut pod and seed sizes are complex quantitative traits regulated by multiple genes and influenced by genetic factors and environmental conditions [10]. Marker-assisted selection (MAS) technology, an important tool for quantitative trait locus (QTL) mapping, has played a crucial role in the screening and detection of QTLs related to peanut pod and seed size. Hopkins et al. (1999) discovered many CT and GT repeat sequences in the cultivated peanut genome. Owing to their high polymorphism, simple sequence repeats (SSR) are particularly well suited for developing molecular markers [11]. Subsequently, Varshney et al. integrated 135 SSR loci to establish the first relatively comprehensive SSR-based molecular linkage map of peanuts [12], significantly advancing the research progress in QTL mapping. With extensive development and application of many SSR markers, researchers have successfully identified numerous QTLs associated with pod and seed size in peanuts [13,14,15,16,17,18]. However, SSR markers are limited by factors such as low polymorphism and marker density, resulting in a lack of precision and efficiency in the localization of candidate genes.
With the rapid advancement of high-throughput sequencing technology, the genome sequences of the tetraploid peanut cultivars A. hypogaea ‘Shitouqi’ [19], A. hypogaea ‘Fuhuasheng’ [20], and A. hypogaea ‘Tifrinner’ [21] have been published to ensure the accuracy of genetic map construction and candidate gene identification. Applying whole-genome re-sequencing and reduced representation sequencing technologies (e.g., RAD-seq and SLAF-seq) has significantly promoted the development of single nucleotide polymorphism (SNP) markers in the peanut genome. High-density SNP markers provide a more refined genomic map and allow for the accurate identification and localization of numerous genomic regions related to pod and seed size [9,22,23,24,25].
Recent studies have identified several candidate genes associated with peanut pod and seed size. Yang et al. finely mapped the major quantitative trait locus qAHPS07 related to pod size to a 36.46 kb interval on chromosome A07. They identified AhRUVBL2 as the most likely candidate gene within this interval [26]. Zhao et al. located the major QTL controlling pod size within an approximately 1.17 Mb interval on chromosome 07 and fine-mapped it to a 20.4 kb interval using F6:7 and F6:8 populations, which contained only one predicted nonsynonymous gene mutation [27].
However, previous QTL mapping studies on traits related to peanut pod and seed size primarily relied on SNP chips and reduced representation sequencing technologies. In contrast, the present study employed whole-genome re-sequencing technology, a more precise and comprehensive approach, which resulted in a genetic map that is more detailed and accurate. This high-resolution genetic map not only aids in the precise localization of key QTL, but also provides a deeper understanding of the underlying genetic mechanisms. The candidate genes linked to peanut pod and seed size remain incompletely discovered, and existing research offers only a restricted comprehension of the genetic foundation of these characteristics. Although several QTL mapping investigations have been carried out, only a limited number of candidate genes have been identified due to the use of somewhat restricted techniques. The complexity of these qualities has posed a challenge in fully elucidating the genetic mechanisms that underlie them. The objective of this work is to conduct high-resolution mapping of quantitative trait loci (QTLs) associated with pod and seed size. The goal is to find more potential genes and enhance the knowledge of the genetic control mechanisms that govern these characteristics.
This study used a high-generation recombinant inbred line (RIL) population derived from a cross between Huayu 44 and DF12. A high-density genetic map was constructed based on whole-genome re-sequencing, and QTL analysis was conducted for hundred pod weight (HPW), hundred seed weight (HSW), pod length (PL), pod width (PW), seed length (SL), and seed width (SW) across six environments. A total of 44 QTLs related to pod and seed size were detected. Notably, a QTL enrichment region was identified at the distal end of chromosome A05. This enrichment region was narrowed down by developing molecular markers for two candidate gene hotspot regions of 600.9 kb and 721.2 kb. Ten candidate genes that potentially regulate pod and seed size were identified from these regions. The experimental results provide a theoretical basis for the genetic improvement of peanut pod and seed size traits through MAS and identifying genes associated with these traits.

2. Materials and Methods

2.1. Plant Materials and Development of RIL Population

By selecting Huayu 44 (a large-seed variety) as the female parent and DF12 (a small-seed variety) as the male parent, a high-generation RIL population comprising 807 lines was constructed using the Single-Seed descent method. Huayu 44 was created by the Shandong Academy of Agricultural Sciences using interspecific hybridization of the farmed species Arachis hypogaea with the biologically incompatible wild species Arachis glabrata Benth. The pods of this cultivar are somewhat huge and it is classified as the ordinary type. DF12 was created by the Industrial Crops Research Institute of the Henan Academy of Agricultural Sciences through a cross between Kaifeng Xuan 01-6 and Baisha 1016. It is an ordinary small-podded kind (Figure 1). The 807 lines and their parents were planted in six different environments. Those grown at the Fen-Yang (FY) experimental station (37°42′ N, 111°79′ E) in 2021, 2022, and 2023 are referred to as 21FY, 22FY, and 23FY, respectively. Those used for the experiment conducted at the Le-Dong (LD) experimental station (18°43′ N, 109°00′ E) in 2021 and 2023 are referred to as 21LD and 23LD, respectively. Experiments were also conducted at the Nan-Bin (NB) experimental station (18°38′ N, 109°21′ E) in 2021, and those grown there are referred to as 21NB (Table S1). Experimental planting in each environment followed a replicate randomized block design, with each row containing 15–20 plants spaced 15 cm apart within rows and 25 cm between rows. Standard agricultural practices were followed for field management during planting.

2.2. Trait Measurements and Data Analysis

After physiological maturity, RIL populations and their parents were harvested individually in rows. After drying in the sun, ten full pods and 20 intact seeds were selected from each family. The traits of the parents and families across the six environments were measured using a Wanshen SC-G Automatic Seed Analysis Instrument (Model: SC-G). The equipment was manufactured by Hangzhou Wanshen Testing Technology Co., Ltd., located in Hangzhou, China.
Based on image recognition principles, the six traits measured were hundred pod weight, pod length, pod width, hundred seed weight, seed length, and seed width. Basic statistical analyses of both the pod and seed traits of the parents and the RIL population, including mean values, standard deviations, and coefficients of variation, were performed using SPSS 20.0. Additionally, correlation coefficients between the six traits were calculated. Normal distribution graphs were plotted using Origin 7.0.

2.3. QTL Analysis

We used the high-density genetic linkage map previously constructed in our laboratory [28]. Combining the phenotypic data of six traits related to pod and seed size from 200 re-sequenced peanuts of the RIL population and the two parents across six planting environments, QTL identification was performed using the composite interval mapping (CIM) model in the R/qtl package [29], with 1000 permutation tests. A log of the odds (LODs) threshold of 3.0 at the 95% significance level was used to declare the presence of a QTL. Positive and negative additive effects indicated favorable alleles in Huayu 44 and DF12, respectively. QTLs with a variance contribution rate of ≥10% were considered major QTLs, while others were considered minor QTLs [18]. A QTL detected in two or more environments was considered stable. QTL naming followed previous conventions, starting with the letter “q” followed by the trait name and chromosome number [30]. If multiple QTLs were detected for the same trait and chromosome, they were distinguished by numbers, such as qHPWA05.1 and qHPWA05.2 for two QTLs related to HPW on chromosome A05.

2.4. Further QTL Mapping Based on KASP Genotyping

Twelve KASP markers were developed based on the identified QTL enrichment regions (Table S2), covering the main QTL enrichment area and its flanking regions. For KASP genotyping, 607 families from the RIL population were selected after excluding 200 re-sequenced families. Ten of the twelve KASP molecular markers were successfully genotyped (Figure S1). The genotyping results and corresponding phenotypic data were used to construct a local genetic linkage map using the MAP module in IciMapping v4.2. QTL mapping was performed using the BIP module with the parameters set to step (cm) = 0.1000.

2.5. Prediction of Candidate Genes in Hotspot Regions

Narrowed intervals on chromosome A05 were defined as hotspot regions for candidate genes regulating peanut pod and seed size. Genes within these hotspot regions were predicted based on the “Tifrunner.gnm2.J5K5” genome sequence. The Candidate genes were annotated using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Non-Redundant Protein (NR) databases. The potential functions of these genes were categorized and analyzed to identify potential candidate genes.

3. Results

3.1. Agronomic Performance of the Parents and RIL Individuals

Using Huayu 44 and DF12 as parental lines, an RIL population consisting of 807 families was constructed. The female parent, Huayu 44, has larger pods and seeds, whereas the male parent, DF12, is a small-seeded variety (Figure 1a). Across the six environments, Huayu 44 showed higher values for HPW, PL, PW, HSW, SL, and SW than DF12 (Figure 1b). The six traits related to pod and seed size in the RIL population exceeded the mean values of the two parents, indicating a wide range of variation and transgressive segregation (Table 1). The population exhibited a continuous distribution and rich genetic diversity. Most phenotypic data had absolute kurtosis and skewness of <1. Combined with the frequency distribution histograms (Figure 2), all traits displayed a normal distribution, characteristic of typical quantitative traits. Significant positive correlations were observed among most traits across the six environments (Table 2), particularly between HPW and HSW, which showed strong correlations in all six environments (0.36~0.95, p < 0.001). The highest correlation was 0.95 (p < 0.001) in the 23LD environment, with HSW and SW, PL and SL, and HPW and SW also highly correlated in this environment with correlation coefficients of 0.87, 0.82, and 0.79, respectively (p < 0.001).

3.2. QTL Mapping for Traits Related to Pod and Seed Size

Our laboratory previously performed whole-genome re-sequencing on 200 of the RIL population and the two parents using the genome of the cultivated peanut variety Tifrunner.gnm2.J5K5 as a reference to construct a genetic linkage map. The total length of this linkage map was 1216.32 cm, comprising 20 linkage groups and 2494 bin markers, which included 447,528 SNPs and 220,443 indels. The average distance between adjacent markers was 0.49 cm. We conducted QTL mapping analysis for traits related to pod and seed sizes. A total of 44 QTLs associated with six traits were detected (Table 3, Figure 3) and distributed across chromosomes A02, A05, B04, B08, and B10. Among these, 40 were major QTLs and four were minor QTLs. The additive effects of these QTLs ranged from −5.94% to 12.8%, LOD values ranged from 4.17 to 22.72, and the phenotypic variance explained (PVE) ranged from 9.16% to 40.73%. These relatively high values reflect the strong stability and reliability of these QTLs. On chromosomes A05, B08, and B10, seven QTLs related to HPW were detected. Among them, qHPWA05.1, which was repeatedly detected in two environments (21FY and 21LD), showed high positive additive effects (12.8% and 11.77%) and explained 26.63% and 19.39% of the phenotypic variance, respectively. Therefore, qHPWA05.1, considered a major and stable QTL, is a key target region for breeding efforts to increase HPW. Eight QTLs related to PL were detected on chromosomes A02, A05, and B04, with all but qPLB04 showing positive additive effects. Among these, qPLA05.2 had the highest PVE for PL at 35.85%, making it the most significant contributor to PL variation. For PW, six QTLs were detected, all located on chromosome A05. qPWA05.2 exhibited the highest positive additive effect (1.37%), explaining 30.65% of the phenotypic variance, thus making the most significant contribution to PW. Eight QTLs related to HSW were detected on chromosomes A05 and B10, all showing positive additive effects (2.48–4.97%). qHSWA05.2 had the highest PVE at 32.20%, whereas qHSWB10.2 had an explained variance of 9.61%, classifying it as a minor QTL. qHSWA05.1, qHSWA05.3, qHSWA05.4, qHSWA05.5, qHSWB10.1, and qHSWB10.3 explained 14.08%, 17.75%, 21.66%, 16.68%, 12.82%, and 10.52% of the phenotypic variance, respectively, all being major QTLs. Six QTLs were detected for SL, located on chromosome A05 and showing positive additive effects. qSLA05.2 had the highest phenotypic variance (40.73%), the highest among all the QTLs for SL. qSLA05.1, qSLA05.3, qSLA05.4, qSLA05.5, and qSLA05.6 explained 22.08%, 37.86%, 25.39%, 24.31%, and 10.54% of the phenotypic variance, respectively. Nine QTLs related to seed width SW were detected on chromosomes A02, A05, B04, B08, and B10, with PVE ranging from 10.54% to 16.00%, all of which were major QTLs. Among these, qSWA05.1 had the highest contribution, explaining 16.00% of the phenotypic variance.

3.3. Further Identification of QTL Enrichment Regions on Chromosome A05

The preliminary QTL mapping results indicate an enriched region of multiple QTLs on chromosome A05 (Figure 4). To further confirm and narrow down the QTL on chromosome A05 related to HPW, PL, PW, HSW, SL, and SW, we constructed a local genetic linkage map based on 10 successfully genotyped KASP molecular markers and further refined the QTL enrichment region on chromosome A05 (Figure 5). Six QTLs were located (Table 4), with additive effect values ranging from −9.89 to −0.20. The negative additive effect values indicated that the favorable alleles for these traits originate from the male parent DF12. High LOD values (7.08 to 20.26) indicated a strong association between the QTL and the target traits. At the same time, the PVE ranged from 10.75% to 27.92%, indicating significant contributions of these major QTLs to the target traits. qHPWA05.6 explained 16.69% of the phenotypic variance for HPW, qPLA05.7 explained 23.90% of the phenotypic variance for PL, qPWA05.7 explained 18.06% of the phenotypic variance for PW, qSLA05.7 explained 18.36% of the phenotypic variance for SL, and qSWA05.5 explained 10.75% of the phenotypic variance for SW. These QTLs were co-located between the markers Ah900023 and Ah900027. The co-located interval was named after the QTL with the highest contribution rate (qPLA05.7). Genes within this interval regulated multiple traits, including HPW, PL, PW, SL, and SW, indicating pleiotropic effects. qHSWA05.6, located between the markers Ah900027 and Ah900026, explained 27.92% of the phenotypic variance in HSW.

3.4. Confirmation and Identification of Candidate Gene Hotspot Regions on Chromosome A05

Under the QTL enrichment region on chromosome A05, we developed ten KASP markers, refining the enrichment area to two QTL regions: qPLA05.7 (721.2 kb) and qHSWA05.6 (600.9 kb). These regions have been identified as hotspots for candidate genes associated with pod and seed size. We annotated 22 candidate genes within the qPLA05.7 region and 27 candidate genes within the qHSWA05.6 region using the GO, KEGG, and NR databases (Table S3). Through GO annotation, 36 of the 49 candidate genes were assigned at least one GO term. The 36 genes were classified into three GO categories: “cellular component”, “molecular function”, and “biological process.” The biological processes these genes are involved in include the methionine biosynthetic process (GO:0009086), threonine biosynthetic process (GO:0009088), isoleucine biosynthetic process (GO:0009097), embryo development ending in seed dormancy (GO:0009793), regulation of transcription, DNA-templated (GO:0006355), sugar-mediated signaling pathway (GO:0010182), protein ubiquitination (GO:0016567), ubiquitin-dependent protein catabolic process (GO:0042787), response to auxin (GO:0009733), response to jasmonic acid (GO:0009753), and polyamine biosynthetic process (GO:0006596). These biological processes involve amino acid synthesis, embryo development, transcription regulation, sugar signal transduction, protein degradation, and plant hormone response. KEGG analysis revealed that these genes encode proteins including chloroplast-like ferredoxin-thioredoxin reductase catalytic chain, 40S ribosomal protein S10, mannose-1-phosphate guanylyltransferase, E3 ubiquitin-protein ligase, cyclic nucleotide-gated ion channel, spermidine synthase, Cullin-1, E2F transcription factor, and lysine-specific histone demethylase. These genes are involved in protein synthesis, modification of cell signaling, and metabolic regulation, highlighting their potential roles in pod and seed development. Using the NR database, we identified the following homologies for the candidate genes: 18 candidate genes were homologous to Glycine max, nine to Cicer arietinum, eight to Medicago truncatula, and seven to Phaseolus vulgaris. Additionally, one gene was homologous to Fragaria vesca subsp. vesca, Pisum sativum, Glycine soja, and Lotus japonicus.

4. Discussion

The present study constructed a new RIL population through the cross between Huayu 44 and DF12. The extensive variation among offspring provided a rich diverse foundation for QTL mapping, facilitating the identification of different genes affecting these traits. Some offspring exhibit trait values beyond the extremes of their parents, which helps in identifying potentially favorable gene combinations and interactions. The six traits also exhibited a normal distribution across the six different environments, providing a solid population foundation for accurately mapping the QTL-controlling pod and seed size. Whole-genome re-sequencing technology covers the entire genome, avoiding gene loss or bias and ensuring the completeness and accuracy of the constructed map. We constructed a high-density genetic map based on whole-genome re-sequencing in the present study. We conducted QTL mapping to ensure the precision and accuracy of the identified QTL.
A total of 44 QTLs were detected across six environments, with 40 major QTLs and four minor QTLs with PVE values close to 10%, indicating a high contribution to the target traits. The high additive effects and PVE values of these QTLs are important for breeding practices because selecting these QTLs can substantially improve target traits. High LOD scores provided reliable statistical evidence for the significance of the detected QTL, laying the foundation for further functional validation and application. Previous studies have identified numerous QTLs related to pod and seed size on chromosome A05. We identified 32 additive QTLs with PVE values up to 40.73% on chromosome A05, all of which were major QTLs, emphasizing the central role of chromosome A05 in the genetic control of pod and seed size. These QTLs highly overlapped, with peaks at the same or adjacent positions, suggesting the presence of one or more gene clusters on chromosome A05 that collectively influence pod and seed sizes through complex genetic mechanisms. Several stable and consistent QTLs detected in the initial mapping had widespread effects on multiple traits of peanut pod and seed size. This phenomenon was also observed in the subsequent validation and fine mapping, indicating that the strong correlations between HPW, PL, PW, HSW, SL, and SW were the main reasons for this phenomenon.
By developing 10 KASP markers on chromosome A05 to cover the overlapping major QTL enrichment region and its flanking areas, we narrowed down the overlapping section to intervals of approximately 600.9 kb between markers Ah900026 and Ah900027 (106,625,271–107,226,203 bp) and approximately 721.2 kb between markers Ah900027 and Ah900023 (107,226,203–107,947,402 bp). Comparatively, these two narrowed intervals overlap with the QTL region of 101.70–111.64 Mb on chromosome A05 identified by Wang [31] and the QTL reported by Luo et al. [18,32]. Identifying similar or overlapping QTL regions in populations with different genetic backgrounds demonstrates that these QTLs generally influence pod and seed size traits rather than being specific to a particular parental genetic background. This further suggests that these regions may contain genes crucial for pod and seed size that are conserved across different parents; that is, the position and function of these genes did not change significantly among the parents. This consistency in gene position and function across different genetic backgrounds indicates that these genes and their effects have been preserved throughout evolution.
Numerous QTLs associated with peanut pod and seed size have been reported on chromosomes A02, A03, A05, A07, A09, A10, B02, B03, B06, B07, and B08 [14,16,22,33,34]. However, in our recent study, five new QTLs were detected on chromosome B10. These findings are significant because QTLs associated with pod and seed traits on chromosome B10 have rarely been reported. Therefore, chromosome B10 is crucial for studying peanuts’ pod- and seed-related traits. This warrants further research with an increased density of molecular markers for precise localization.
Using GO, KEGG, and NR databases annotated 49 candidate genes within the identified hotspot regions. By comparing these genes with homologous genes in other species with related functions, we identified 10 candidate genes that may regulate pod and seed size. These genes encode cytochrome P450, spermidine synthase, mannose-1-phosphate guanylyltransferase, pentatricopeptide repeat proteins, and E2F transcription factors. The genes Arahy.QW591P, Arahy.Q9F5SX, Arahy.Z6EAYS, Arahy.4A049R, and Arahy.TNFT39 encode proteins belonging to the cytochrome P450 (CYP450) family. As one of the largest plant protein families, CYP450 proteins have been found to play a significant role in seed development in various crops, such as wheat [35], tomato [36], and rapeseed [37]. Adamski et al. found that the CYP450 protein encoded by KLUH determines the growth potential of the seed coat by promoting cell proliferation in both the inner and outer integuments, thereby affecting final seed size [38]. Fang et al. discovered that EOD3 encodes a CYP450 protein in Arabidopsis, and its overexpression significantly increases seed size [39]. Ma et al. found that the activity of wheat TaCYP78A3 is positively correlated with the final seed size, as it influences the number of cells in the seed coat. Silencing TaCYP78A3 reduces seed coat cell numbers, resulting in an 11% decrease in seed size, whereas overexpression of TaCYP78A3 in Arabidopsis induces more seed coat cells, increasing seed size by 11–48% [35]. Liu et al. utilized multi-omics analysis to identify a candidate gene encoding the CYP450 protein associated with increased peanut seed size. They hypothesized that the significant upregulation of CYP450 may contribute to larger seed size [40]. Based on existing research evidence, we hypothesize that the genes mentioned above encoding CYP450 proteins (Arahy.QW591P, Arahy.Q9F5SX, Arahy.Z6EAYS, Arahy.4A049R, and Arahy.TNFT39) may significantly influence seed size by regulating the number and size of cells within the seed coat and other related organs during seed development.
Arahy.DY91G5 encodes spermidine synthase (SPDS). Studies have shown that the OsSPMS1 gene in rice plays a crucial role in converting spermidine (SPD) to spermine. Overexpression of OsSPMS1 leads to the excessive consumption of SPD, thereby negatively regulating seed germination, grain size, and yield of rice [41]. Recent studies have also found that the overexpression of the TaSPDS-7D1 gene, which encodes SPDS in wheat, significantly increases the thousand-grain weight of wheat seeds. Therefore, the SPD expression level has a crucial effect on seed size [42].
Arahy.Q67Z2R encodes an E2F transcription factor that plays a vital role in cell-cycle regulation. Studies have shown that the activators E2FA and E2FB primarily regulate cell cycle progression. These factors also inhibit the premature expression of seed maturation-related genes and coordinate the orderly progression of cell proliferation and maturation during seed development [43]. Therefore, E2F transcription factors are critical regulators of seed development.
Arahy.SEJH4Y encodes mannose-1-phosphate guanylyltransferase, which mediates N-glycosylation during cellulose biosynthesis. Cellulose is a crucial component of peanut seed coats and pods [44]. Thus, mannose-1-phosphate guanylyltransferase may influence the development of seed coats and pods by regulating cellulose synthesis and affecting peanut pod size.
Arahy.7TY6PE and Arahy.X5R2L4 encode pentatricopeptide repeat (PPR) proteins, whose functional impairments can lead to abnormal phenotypes related to embryonic development, photosynthetic capacity, and seed pigment deposition [45]. In maize, qKW9, which encodes a PPR protein, regulates kernel size and weight by influencing photosynthesis and grain filling [46]. Recent studies have identified several PPR-encoding genes as candidates for regulating seed size in peanuts. Gangurde et al. identified several candidate genes related to peanut SW using linkage and association analyses, one of which was predicted to encode a PPR protein [44]. Using transcriptome analysis combined with previously reported QTL data, Li et al. predicted the PPR protein-encoding gene EVM0025654 to be a candidate gene associated with peanut seed size variation [45]. Therefore, the PPR protein-encoding genes rahy.7TY6PE and Arahy.X5R2L4 identified in this study may be potential candidate genes associated with peanut pod and seed size variations.
Notably, several key candidate genes identified in this study, including those encoding SPDS, mannose-1-phosphate guanylyltransferase, PPR proteins, and E2F transcription factors, were also mentioned by Gangurde et al. in their identification of candidate genes related to peanut pod and seed weight [44]. Although we have identified several candidate genes within the hotspot regions that may regulate peanut pod and seed size, there remains some uncertainty due to the lack of functional validation. Functional validation of these candidate genes is crucial in confirming their role in trait regulation. It is possible to further elucidate whether and how these genes influence pod and seed size through functional experiments such as gene knockout, overexpression, or gene editing.

5. Conclusions

This study constructed a high-density genetic map through whole-genome re-sequencing using a high-generation RIL population derived from a Huayu 44 and DF12 cross. QTL analysis was then employed to identify potential genomic regions associated with pod and seed size traits found on five chromosomes (A02, A05, B04, B08, and B10). In addition, a QTL-rich region was discovered on chromosome A05. To further narrow the interval, ten molecular markers were developed to cover the QTL-rich region and its flanking regions, reducing the interval lengths to 600.9 kb and 721.2 kb. These two QTL regions were identified as candidate gene hotspots regulating pod and seed size. Through prediction and annotation of candidate genes within these hotspot regions, 10 genes were proposed as candidate genes for regulating pod and seed size. These include those encoding CYP450 protein (Arahy.QW591P, Arahy.Q9F5SX, Arahy.Z6EAYS, Arahy.4A049R, and Arahy.TNFT39), polyamine synthase (Arahy.DY91G5), mannose-1-phosphate guanylyltransferase (Arahy.SEJH4Y), pentapeptide repeat proteins (Arahy.7TY6PE and Arahy.X5R2L4), and E2F transcription factor (Arahy.Q67Z2R). Identifying and characterizing these candidate genes will enhance our understanding of peanuts’ molecular mechanisms underlying pod and seed size traits. Additionally, the molecular markers developed in this study can serve as potential selection markers, laying the foundation for marker-assisted breeding.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture14091634/s1, Table S1: annual average temperature and annual average precipitation in six environments; Table S2: The primers of 12 SNP-based KASP markers; Table S3: Annotation of candidate genes on spot region on chromosome A05; Figure S1: Genotyping profiles of the 12 developed KASP markers.

Author Contributions

Conceptualization, supervision, writing—review and editing, project administration, funding acquisition, D.B and X.Z. (Xin Zhang).; data curation, formal analysis, writing—original draft preparation, visualization, X.Z. (Xiaoji Zhang); methodology, H.Z. and N.L.; validation, Y.T. and Y.X.; investigation, X.Z. (Xiaoji Zhang), L.W., Q.L., and X.Z. (Xiaoyu Zhang); resources, D.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Major Project of Shanxi Province (No.: 202201140601025), the National Peanut Industry Technology System Construction (No.: CARS-13), The Biological Breeding Engineering Project of Shanxi Agricultural University (No.: YZGC049), the earmarked fund for Modern Agro-industry Technology Research System (No.: 2024CYJSTX05), and the Fundamental Research Program of Shanxi Province (No.: 202203021221178).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviation

AbbreviationFull Term
CIMComposite Interval Mapping
KASPKompetitive Allele-Specific PCR
PVEPhenotypic Variance Explained
QTLQuantitative Trait Loci
RAD-seqRestriction site-Associated DNA sequencing
RILRecombinant Inbred Line
SLAF-seqSpecific-Locus Amplified Fragment sequencing
SNPSingle Nucleotide Polymorphism

References

  1. Guo, M.; Deng, L.; Gu, J.; Miao, J.; Yin, J.; Li, Y.; Fang, Y.; Huang, B.; Sun, Z.; Qi, F. Genome-wide association study and development of molecular markers for yield and quality traits in peanut (Arachis hypogaea L.). BMC Plant Biol. 2024, 24, 244. [Google Scholar] [CrossRef] [PubMed]
  2. Bertioli, D.J.; Cannon, S.B.; Froenicke, L.; Huang, G.; Farmer, A.D.; Cannon, E.K.; Liu, X.; Gao, D.; Clevenger, J.; Dash, S. The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat. Genet. 2016, 48, 438–446. [Google Scholar] [CrossRef] [PubMed]
  3. Cuc, L.M.; Mace, E.S.; Crouch, J.H.; Quang, V.D.; Long, T.D.; Varshney, R.K. Isolation and characterization of novel microsatellite markers and their application for diversity assessment in cultivated groundnut (Arachis hypogaea). BMC Plant Biol. 2008, 8, 55. [Google Scholar] [CrossRef] [PubMed]
  4. Kochert, G.; Stalker, H.T.; Gimenes, M.; Galgaro, L.; Lopes, C.R.; Moore, K. RFLP and cytogenetic evidence on the origin and evolution of allotetraploid domesticated peanut, Arachis hypogaea (Leguminosae). Am. J. Bot. 1996, 83, 1282–1291. [Google Scholar] [CrossRef]
  5. Seijo, J.G.; Lavia, G.I.; Fernández, A.; Krapovickas, A.; Ducasse, D.; Moscone, E.A. Physical mapping of the 5S and 18S–25S rRNA genes by FISH as evidence that Arachis duranensis and A. ipaensis are the wild diploid progenitors of A. hypogaea (Leguminosae). Am. J. Bot. 2004, 91, 1294–1303. [Google Scholar] [CrossRef]
  6. Robledo, G.; Lavia, G.; Seijo, G. Species relations among wild Arachis species with the A genome as revealed by FISH mapping of rDNA loci and heterochromatin detection. Theor. Appl. Genet. 2009, 118, 1295–1307. [Google Scholar] [CrossRef]
  7. Chu, Y.; Chee, P.; Isleib, T.G.; Holbrook, C.C.; Ozias-Akins, P. Major seed size QTL on chromosome A05 of peanut (Arachis hypogaea) is conserved in the US mini core germplasm collection. Mol. Breed. 2020, 40, 6. [Google Scholar] [CrossRef]
  8. Wu, Y.; Sun, Z.; Qi, F.; Tian, M.; Wang, J.; Zhao, R.; Wang, X.; Wu, X.; Shi, X.; Liu, H.; et al. Comparative transcriptomics analysis of developing peanut (Arachis hypogaea L.) pods reveals candidate genes affecting peanut seed size. Front. Plant Sci. 2022, 13, 958808. [Google Scholar] [CrossRef] [PubMed]
  9. Chavarro, C.; Chu, Y.; Holbrook, C.; Isleib, T.; Bertioli, D.; Hovav, R.; Butts, C.; Lamb, M.; Sorensen, R.; Jackson, S.A. Pod and seed trait QTL identification to assist breeding for peanut market preferences. G3 Genes Genomes Genet. 2020, 10, 2297–2315. [Google Scholar] [CrossRef]
  10. Xiang, D.-Q.; Cao, H.-H.; Cao, Y.-G.; Yang, J.-P.; Huang, L.-J.; Wang, S.-C.; Dai, J. Construction of a genetic map and location of quantitative trait loci for yield component traits in maize by SSR markers. Yi Chuan Xue Bao = Acta Genet. Sin. 2001, 28, 778–784. [Google Scholar]
  11. Hopkins, M.; Casa, A.; Wang, T.; Mitchell, S.; Dean, R.; Kochert, G.; Kresovich, S. Discovery and characterization of polymorphic simple sequence repeats (SSRs) in peanut. Crop Sci. 1999, 39, 1243–1247. [Google Scholar] [CrossRef]
  12. Varshney, R.; Bertioli, D.; Moretzsohn, M.d.C.; Vadez, V.; Krishnamurthy, L.; Aruna, R.; Nigam, S.; Moss, B.; Seetha, K.; Ravi, K. The first SSR-based genetic linkage map for cultivated groundnut (Arachis hypogaea L.). Theor. Appl. Genet. 2009, 118, 729–739. [Google Scholar] [CrossRef] [PubMed]
  13. Ravi, K.; Vadez, V.; Isobe, S.; Mir, R.; Guo, Y.; Nigam, S.; Gowda, M.; Radhakrishnan, T.; Bertioli, D.; Knapp, S. Identification of several small main-effect QTL and a large number of epistatic QTL for drought tolerance related traits in groundnut (Arachis hypogaea L.). Theor. Appl. Genet. 2011, 122, 1119–1132. [Google Scholar] [CrossRef]
  14. Shirasawa, K.; Koilkonda, P.; Aoki, K.; Hirakawa, H.; Tabata, S.; Watanabe, M.; Hasegawa, M.; Kiyoshima, H.; Suzuki, S.; Kuwata, C. In silico polymorphism analysis for the development of simple sequence repeat and transposon markers and construction of linkage map in cultivated peanut. BMC Plant Biol. 2012, 12, 80. [Google Scholar] [CrossRef] [PubMed]
  15. Fonceka, D.; Tossim, H.-A.; Rivallan, R.; Vignes, H.; Faye, I.; Ndoye, O.; Moretzsohn, M.C.; Bertioli, D.J.; Glaszmann, J.-C.; Courtois, B. Fostered and left behind alleles in peanut: Interspecific QTL mapping reveals footprints of domestication and useful natural variation for breeding. BMC Plant Biol. 2012, 12, 26. [Google Scholar] [CrossRef]
  16. Huang, L.; He, H.; Chen, W.; Ren, X.; Chen, Y.; Zhou, X.; Xia, Y.; Wang, X.; Jiang, X.; Liao, B. Quantitative trait locus analysis of agronomic and quality-related traits in cultivated peanut (Arachis hypogaea L.). Theor. Appl. Genet. 2015, 128, 1103–1115. [Google Scholar] [CrossRef]
  17. Chen, W.; Jiao, Y.; Cheng, L.; Huang, L.; Liao, B.; Tang, M.; Ren, X.; Zhou, X.; Chen, Y.; Jiang, H. Quantitative trait locus analysis for pod-and kernel-related traits in the cultivated peanut (Arachis hypogaea L.). BMC Genet. 2016, 17, 25. [Google Scholar] [CrossRef]
  18. Luo, H.; Ren, X.; Li, Z.; Xu, Z.; Li, X.; Huang, L.; Zhou, X.; Chen, Y.; Chen, W.; Lei, Y. Co-localization of major quantitative trait loci for pod size and weight to a 3.7 cM interval on chromosome A05 in cultivated peanut (Arachis hypogaea L.). BMC Genom. 2017, 18, 58. [Google Scholar] [CrossRef]
  19. Zhuang, W.; Chen, H.; Yang, M.; Wang, J.; Pandey, M.K.; Zhang, C.; Chang, W.-C.; Zhang, L.; Zhang, X.; Tang, R. The genome of cultivated peanut provides insight into legume karyotypes, polyploid evolution and crop domestication. Nat. Genet. 2019, 51, 865–876. [Google Scholar] [CrossRef]
  20. Chen, X.; Lu, Q.; Liu, H.; Zhang, J.; Hong, Y.; Lan, H.; Li, H.; Wang, J.; Liu, H.; Li, S. Sequencing of cultivated peanut, Arachis hypogaea, yields insights into genome evolution and oil improvement. Mol. Plant 2019, 12, 920–934. [Google Scholar] [CrossRef]
  21. Bertioli, D.J.; Jenkins, J.; Clevenger, J.; Dudchenko, O.; Gao, D.; Seijo, G.; Leal-Bertioli, S.C.; Ren, L.; Farmer, A.D.; Pandey, M.K. The genome sequence of segmental allotetraploid peanut Arachis hypogaea. Nat. Genet. 2019, 51, 877–884. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, S.; Hu, X.; Miao, H.; Chu, Y.; Cui, F.; Yang, W.; Wang, C.; Shen, Y.; Xu, T.; Zhao, L. QTL identification for seed weight and size based on a high-density SLAF-seq genetic map in peanut (Arachis hypogaea L.). BMC Plant Biol. 2019, 19, 537. [Google Scholar] [CrossRef]
  23. Alyr, M.H.; Pallu, J.; Sambou, A.; Nguepjop, J.R.; Seye, M.; Tossim, H.-A.; Djiboune, Y.R.; Sane, D.; Rami, J.-F.; Fonceka, D. Fine-mapping of a wild genomic region involved in pod and seed size reduction on chromosome A07 in peanut (Arachis hypogaea L.). Genes 2020, 11, 1402. [Google Scholar] [CrossRef]
  24. Yang, Y.; Su, Q.; Li, Y.; Cheng, Z.; Song, Y.; Jin, X.; Wang, J. Fine mapping of a major QTL qHYF_B06 for peanut yield. Crop J. 2023, 11, 1533–1540. [Google Scholar] [CrossRef]
  25. Lv, Z.; Lan, G.; Bai, B.; Yu, P.; Wang, C.; Zhang, H.; Zhong, C.; Zhao, X.; Yu, H. Identification of candidate genes associated with peanut pod length by combined analysis of QTL-seq and RNA-seq. Genomics 2024, 116, 110835. [Google Scholar] [CrossRef]
  26. Yang, H.; Luo, L.; Li, Y.; Li, H.; Zhang, X.; Zhang, K.; Zhu, S.; Li, X.; Li, Y.; Wan, Y. Fine mapping of qAHPS07 and functional studies of AhRUVBL2 controlling pod size in peanut (Arachis hypogaea L.). Plant Biotechnol. J. 2023, 21, 1785–1798. [Google Scholar] [CrossRef] [PubMed]
  27. Zhao, K.; Wang, L.; Qiu, D.; Cao, Z.; Wang, K.; Li, Z.; Wang, X.; Wang, J.; Ma, Q.; Cao, D. PSW1, an LRR receptor kinase, regulates pod size in peanut. Plant Biotechnol. J. 2023, 21, 2113–2124. [Google Scholar] [CrossRef]
  28. Zhang, X.; Zhang, X.; Wang, L.; Liu, Q.; Liang, Y.; Zhang, J.; Xue, Y.; Tian, Y.; Zhang, H.; Li, N. Fine mapping of a QTL and identification of candidate genes associated with cold tolerance during germination in peanut (Arachis hypogaea L.) on chromosome B09 using whole genome re-sequencing. Front. Plant Sci. 2023, 14, 1153293. [Google Scholar] [CrossRef]
  29. Arends, D.; Prins, P.; Jansen, R.C.; Broman, K.W. R/qtl: High-throughput multiple QTL mapping. Bioinformatics 2010, 26, 2990–2992. [Google Scholar] [CrossRef]
  30. Li, L.; Yang, X.; Cui, S.; Meng, X.; Mu, G.; Hou, M.; He, M.; Zhang, H.; Liu, L.; Chen, C.Y. Construction of high-density genetic map and mapping quantitative trait loci for growth habit-related traits of peanut (Arachis hypogaea L.). Front. Plant Sci. 2019, 10, 745. [Google Scholar] [CrossRef]
  31. Wang, Z.; Yan, L.; Chen, Y.; Wang, X.; Huai, D.; Kang, Y.; Jiang, H.; Liu, K.; Lei, Y.; Liao, B. Detection of a major QTL and development of KASP markers for seed weight by combining QTL-seq, QTL-mapping and RNA-seq in peanut. Theor. Appl. Genet. 2022, 135, 1779–1795. [Google Scholar] [CrossRef] [PubMed]
  32. Luo, H.; Guo, J.; Ren, X.; Chen, W.; Huang, L.; Zhou, X.; Chen, Y.; Liu, N.; Xiong, F.; Lei, Y. Chromosomes A07 and A05 associated with stable and major QTL for pod weight and size in cultivated peanut (Arachis hypogaea L.). Theor. Appl. Genet. 2018, 131, 267–282. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, Z.; Huai, D.; Zhang, Z.; Cheng, K.; Kang, Y.; Wan, L.; Yan, L.; Jiang, H.; Lei, Y.; Liao, B. Development of a high-density genetic map based on specific length amplified fragment sequencing and its application in quantitative trait loci analysis for yield-related traits in cultivated peanut. Front. Plant Sci. 2018, 9, 827. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, Y.; Ren, X.; Zheng, Y.; Zhou, X.; Huang, L.; Yan, L.; Jiao, Y.; Chen, W.; Huang, S.; Wan, L. Genetic mapping of yield traits using RIL population derived from Fuchuan Dahuasheng and ICG6375 of peanut (Arachis hypogaea L.). Mol. Breed. 2017, 37, 17. [Google Scholar] [CrossRef]
  35. Ma, M.; Wang, Q.; Li, Z.; Cheng, H.; Li, Z.; Liu, X.; Song, W.; Appels, R.; Zhao, H. Expression of Ta CYP 78A3, a gene encoding cytochrome P450 CYP 78A3 protein in wheat (Triticum aestivum L.), affects seed size. Plant J. 2015, 83, 312–325. [Google Scholar] [CrossRef]
  36. Chakrabarti, M.; Zhang, N.; Sauvage, C.; Munos, S.; Blanca, J.; Cañizares, J.; Diez, M.J.; Schneider, R.; Mazourek, M.; McClead, J. A cytochrome P450 regulates a domestication trait in cultivated tomato. Proc. Natl. Acad. Sci. USA 2013, 110, 17125–17130. [Google Scholar] [CrossRef]
  37. Zhou, X.; Zhang, H.; Wang, P.; Liu, Y.; Zhang, X.; Song, Y.; Wang, Z.; Ali, A.; Wan, L.; Yang, G. BnaC7. ROT3, the causal gene of cqSL-C7, mediates silique length by affecting cell elongation in Brassica napus. J. Exp. Bot. 2022, 73, 154–167. [Google Scholar] [CrossRef]
  38. Adamski, N.M.; Anastasiou, E.; Eriksson, S.; O’Neill, C.M.; Lenhard, M. Local maternal control of seed size by KLUH/CYP78A5-dependent growth signaling. Proc. Natl. Acad. Sci. USA 2009, 106, 20115–20120. [Google Scholar] [CrossRef]
  39. Fang, W.; Wang, Z.; Cui, R.; Li, J.; Li, Y. Maternal control of seed size by EOD3/CYP78A6 in Arabidopsis thaliana. Plant J. 2012, 70, 929–939. [Google Scholar] [CrossRef]
  40. Liu, Y.; Yi, C.; Liu, Q.; Wang, C.; Wang, W.; Han, F.; Hu, X. Multi-omics profiling identifies candidate genes controlling seed size in peanut. Plants 2022, 11, 3276. [Google Scholar] [CrossRef]
  41. Tao, Y.; Wang, J.; Miao, J.; Chen, J.; Wu, S.; Zhu, J.; Zhang, D.; Gu, H.; Cui, H.; Shi, S. The spermine synthase OsSPMS1 regulates seed germination, grain size, and yield. Plant Physiol. 2018, 178, 1522–1536. [Google Scholar] [CrossRef] [PubMed]
  42. Ren, J.; Li, C.; Xiu, Q.; Xu, M.; Liu, H. Overexpression of wheat spermidine synthase gene enhances wheat resistance to Fusarium head blight. Phytopathol. Res. 2024, 6, 24. [Google Scholar] [CrossRef]
  43. Leviczky, T.; Molnár, E.; Papdi, C.; Őszi, E.; Horváth, G.V.; Vizler, C.; Nagy, V.; Pauk, J.; Bögre, L.; Magyar, Z. E2FA and E2FB transcription factors coordinate cell proliferation with seed maturation. Development 2019, 146, dev179333. [Google Scholar] [CrossRef] [PubMed]
  44. Gangurde, S.S.; Wang, H.; Yaduru, S.; Pandey, M.K.; Fountain, J.C.; Chu, Y.; Isleib, T.; Holbrook, C.C.; Xavier, A.; Culbreath, A.K. Nested-association mapping (NAM)-based genetic dissection uncovers candidate genes for seed and pod weights in peanut (Arachis hypogaea). Plant Biotechnol. J. 2020, 18, 1457–1471. [Google Scholar] [CrossRef] [PubMed]
  45. Li, Z.; Zhang, X.; Zhao, K.; Zhao, K.; Qu, C.; Gao, G.; Gong, F.; Ma, X.; Yin, D. Comprehensive transcriptome analyses reveal candidate genes for variation in seed size/weight during peanut (Arachis hypogaea L.) domestication. Front. Plant Sci. 2021, 12, 666483. [Google Scholar] [CrossRef]
  46. Huang, J.; Lu, G.; Liu, L.; Raihan, M.S.; Xu, J.; Jian, L.; Zhao, L.; Tran, T.M.; Zhang, Q.; Liu, J. The kernel size-related quantitative trait locus qKW9 encodes a pentatricopeptide repeat protein that affects photosynthesis and grain filling. Plant Physiol. 2020, 183, 1696–1709. [Google Scholar] [CrossRef]
Figure 1. Phenotypic characterization of seeds from the female parent Huayu 44 and the male parent DF12: (a) Pod and seed morphologies of the two parents. (b) Comparisons of hundred pod weight (HPW), pod length (PL), pod width (PW), hundred seed weight (HSW), seed length (SL), and seed width (SW) between Huayu 44 and DF12. Data are shown as mean ± standard error. p values were generated using Student’s t-test.
Figure 1. Phenotypic characterization of seeds from the female parent Huayu 44 and the male parent DF12: (a) Pod and seed morphologies of the two parents. (b) Comparisons of hundred pod weight (HPW), pod length (PL), pod width (PW), hundred seed weight (HSW), seed length (SL), and seed width (SW) between Huayu 44 and DF12. Data are shown as mean ± standard error. p values were generated using Student’s t-test.
Agriculture 14 01634 g001
Figure 2. Phenotypic distributions of six traits related to pod and seed size in the RIL population under six environments. P1 represents the phenotypic distribution of DF12, and P2 represents the phenotypic distribution of Huayu 44.
Figure 2. Phenotypic distributions of six traits related to pod and seed size in the RIL population under six environments. P1 represents the phenotypic distribution of DF12, and P2 represents the phenotypic distribution of Huayu 44.
Agriculture 14 01634 g002
Figure 3. QTL mapping for traits related to pod and seed size.
Figure 3. QTL mapping for traits related to pod and seed size.
Agriculture 14 01634 g003
Figure 4. Distribution of QTLs related to pod and seed size on the genetic linkage map. The distributions of recombination hotspots and coldspots on the genetic map are shown by the color blocks along the horizontal axis. Recombination density increases with color; conversely, it decreases with color. Redder colors indicate higher recombination densities.
Figure 4. Distribution of QTLs related to pod and seed size on the genetic linkage map. The distributions of recombination hotspots and coldspots on the genetic map are shown by the color blocks along the horizontal axis. Recombination density increases with color; conversely, it decreases with color. Redder colors indicate higher recombination densities.
Agriculture 14 01634 g004
Figure 5. Further narrowing of the QTL related to pod and seed size in the enriched region. The red region represents qPLA05.7, and the green region represents qHSWA05.6.
Figure 5. Further narrowing of the QTL related to pod and seed size in the enriched region. The red region represents qPLA05.7, and the green region represents qHSWA05.6.
Agriculture 14 01634 g005
Table 1. Phenotypic variations for six traits of the RIL population and their parents in six environments.
Table 1. Phenotypic variations for six traits of the RIL population and their parents in six environments.
TraitEnv.ParentsRIL Population
DF12HY44MaxMinMean ± SDCV/%KurtSkew
HPW20FY132.64158.65356.1577.80140.62 ± 35.5125.2514.933.24
(g)21FY162.97192.86203.2365.50149.30 ± 24.6916.540.02−0.36
21LD127.86154.54220.9236.10126.57 ± 26.2420.731.48−0.02
21NB112.37140.68180.2462.71123.27 ± 21.1017.110.03−0.12
22FY94.01126.99179.8177.10126.52 ± 19.5115.430.330.06
23LD135.66213.91236.88111.26166.74 ± 23.7714.250.020.32
PL20FY29.8133.4443.8526.6234.40 ± 3.088.96−0.170.10
(mm)21FY31.2234.6539.4027.0633.44 ± 2.808.36−0.73−0.14
21LD29.8534.7841.9915.7032.15 ± 3.5811.121.98−0.63
21NB29.7733.6941.7722.4033.10 ± 3.4610.44−0.31−0.18
22FY29.4934.7540.0425.2731.94 ± 2.648.26−0.100.09
23LD29.7833.8737.7225.4932.25 ± 2.527.83−0.29−0.26
PW20FY12.4814.8519.0411.7014.91 ± 1.157.740.940.48
(mm)21FY15.9816.2918.0712.2514.98 ± 1.348.95−0.810.21
21LD14.4515.4617.9010.7014.54 ± 1.268.690.060.23
21NB12.2514.7717.4611.8914.59 ± 1.077.310.230.42
22FY12.5715.6717.8411.5614.48 ± 1.238.49−0.230.26
23LD13.3715.4317.2411.8814.59 ± 0.976.62−0.100.19
HSW20FY54.1574.3677.4633.1556.45 ± 8.5815.20−0.40−0.05
(g)21FY60.3475.8785.8233.2665.07 ± 8.6913.360.22−0.20
21LD48.1285.8592.3032.3061.50 ± 12.2619.93−0.390.04
21NB62.5776.4581.9329.8055.49 ± 9.4317.000.02−0.04
22FY33.1944.6966.0910.2946.48 ± 8.1217.481.46−0.38
23LD55.7089.9498.2444.3868.54 ± 10.0014.60−0.100.29
SL20FY14.9216.2418.5510.8915.27 ± 1.278.300.05−0.15
(mm)21FY15.4716.6018.7512.3415.81 ± 1.217.65−0.43−0.17
21LD13.7918.1118.4110.5815.47 ± 1.469.44−0.15−0.30
21NB13.9117.7918.1311.6815.22 ± 1.308.53−0.53−0.23
22FY12.9015.3517.3711.8314.48 ± 1.147.85−0.370.04
23LD14.5216.7918.0112.8515.50 ± 1.107.10−0.47−0.23
SW20FY8.4310.019.756.328.62 ± 0.546.281.08−0.56
(mm)21FY8.8110.3210.817.549.31 ± 0.515.440.490.04
21LD8.1810.2210.496.558.78 ± 0.707.990.14−0.22
21NB8.129.3510.926.768.73 ± 0.647.350.53−0.17
22FY6.967.809.356.478.06 ± 0.556.830.05−0.24
23LD8.7711.0411.107.879.54 ± 0.636.63−0.300.18
The female parent is Huayu 44 and the male parent is DF12. Env. is the environment. The unit of measurement for HPW and HSW is grams. The unit of measurement for PL, PW, SL, and SW is centimeters. Max, maximum; Min, minimum; SD, standard deviation; cv is the coefficient of variation; Skew, skewness; Kurt, kurtosis.
Table 2. Correlation analysis of six traits in the RIL population under six environments.
Table 2. Correlation analysis of six traits in the RIL population under six environments.
EnvironmentTraitHPWPLPWHSWSLSW
20FYHPW1
PL0.20 **1
PW0.37 ***0.64 ***1
HSW0.36 ***0.39 ***0.44 ***1
SL0.28 ***0.70 ***0.52 ***0.76 ***1
SW0.28 ***0.0450.20 **0.83 ***0.38 ***1
21FYHPW1
PL0.66 ***1
PW0.65 ***0.65 ***1
HSW0.85 ***0.62 ***0.61 ***1
SL0.70 ***0.86 ***0.65 ***0.77 ***1
SW0.68 ***0.26 ***0.44 ***0.83 ***0.41 ***1
21LDHPW1
PL0.71 ***1
PW0.68 ***0.68 ***1
HSW0.83 ***0.55 ***0.57 ***1
SL0.69 ***0.77 ***0.59 ***0.78 ***1
SW0.65 ***0.19 **0.42 ***0.84 ***0.42 ***1
21NBHPW1
PL0.58 ***1
PW0.64 ***0.65 ***1
HSW0.82 ***0.51 ***0.57 ***1
SL0.44 ***0.63 ***0.44 ***0.52 ***1
SW0.46 ***0.070.27 ***0.61 ***0.25 ***1
22FYHPW1
PL0.65 ***1
PW0.62 ***0.72 ***1
HSW0.90 ***0.53 ***0.48 ***1
SL0.74 ***0.85 ***0.64 ***0.71 ***1
SW0.77 ***0.23 ***0.35 ***0.81 ***0.43 ***1
23LDHPW1
PL0.54 ***1
PW0.63 ***0.69 ***1
HSW0.95 ***0.44 ***0.55 ***1
SL0.63 ***0.82 ***0.63 ***0.61 ***1
SW0.79 ***0.0930.37 ***0.87 ***0.29 ***1
The correlations were calculated with Pearson correlation coefficients, and the p values are indicated as follows: **, p < 0.01; ***, p < 0.001.
Table 3. QTL information for traits related to pod and seed size.
Table 3. QTL information for traits related to pod and seed size.
TraitQTLChr.Range (cm)LODEnv.AdditivePVE%
HPWqHPWA05.1A0529.91–49.5213.4521FY12.8026.63
qHPWA05.2A0534.81–49.5211.7122FY9.5223.64
qHPWA05.3A0529.91–49.429.3621LD11.7719.39
qHPWA05.4A0530.17–48.9911.1623LD11.4122.66
qHPWA05.5A0535.34–48.725.9521NB7.6812.81
qHPWB08B080–0.804.1822FY−5.949.18
qHPWB10B1023.66–24.764.1721FY7.509.16
PLqPLA02A0267.60–68.374.2121NB1.059.24
qPLA05.1A0524.34–48.197.8820FY1.2616.59
qPLA05.2A0525.14–52.3219.2821FY1.6835.85
qPLA05.3A0521.43–49.5213.4822FY1.3726.82
qPLA05.4A0527.78–48.199.4121LD1.6119.49
qPLA05.5A0527.78–49.5212.6223LD1.2825.22
qPLA05.6A0526.73–52.329.3621NB1.5519.38
qPLB04B0435.65–37.534.8321FY−0.9110.52
PWqPWA05.1A0534.02–48.197.8820FY0.4716.60
qPWA05.2A0529.91–49.5215.8921FY1.3730.65
qPWA05.3A0530.96–49.5213.2622FY0.6326.31
qPWA05.4A0525.14–52.3215.5121LD0.7030.04
qPWA05.5A0535.87–49.5210.9123LD0.4622.22
qPWA05.6A0534.02–52.3210.2621NB0.4921.04
HSWqHSWA05.1A0534.81–46.336.5920FY3.2914.08
qHSWA05.2A0527.8–49.5216.8821FY4.9732.20
qHSWA05.3A0535.34–47.398.4922FY3.4417.75
qHSWA05.4A0529.91–49.5210.6023LD4.7021.66
qHSWA05.5A0534.02–49.527.9221NB3.9016.68
qHSWB10.1B1020.64–35.725.9521FY3.1312.82
qHSWB10.2B1023.12–26.914.3922FY2.489.61
qHSWB10.3B1020.64–30.134.8321NB3.1010.52
SLqSLA05.1A0530.17–52.3210.8320FY0.6022.08
qSLA05.2A0521.17–52.3222.7221FY0.7740.73
qSLA05.3A0519.84–52.3220.6622FY0.7037.86
qSLA05.4A0523.82–49.5212.7221LD0.7425.39
qSLA05.5A0521.17–48.1912.1023LD0.5424.31
qSLA05.6A0538.22–44.254.8421NB0.4210.54
SWqSWA02A0266.49–68.375.3523LD−0.2111.60
qSWA05.1A0530.70–46.597.5721FY0.2016.00
qSWA05.2A0536.92–43.735.1822FY0.1811.25
qSWA05.3A0537.18–48.995.8821LD0.2512.67
qSWA05.4A0538.75–45.555.9023LD0.2312.69
qSWB04.1B0447.79–51.784.7322FY0.1810.31
qSWB04.2B0432.43–51.255.4021LD0.2411.70
qSWB08B080–1.074.9222FY−0.1810.71
qSWB10B1020.64–35.454.8421FY0.1710.54
Table 4. Information on QTL related to pod and seed size after further narrowing of the enriched region.
Table 4. Information on QTL related to pod and seed size after further narrowing of the enriched region.
TraitChromosomeMarker IntervalLODAddPVE(%)
qHPWA05.6A05Ah900023–Ah90002711.31−9.8916.69
qPLA05.7A05Ah900023–Ah90002716.82−1.3223.90
qPWA05.7A05Ah900023–Ah90002712.43−0.5618.06
qHSWA05.6A05Ah900027–Ah90002620.26−0.6527.92
qSLA05.7A05Ah900023–Ah90002712.15−4.1018.36
qSWA05.5A05Ah900023–Ah9000277.08−0.2010.75
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, X.; Wang, L.; Liu, Q.; Zhang, X.; Tian, Y.; Xue, Y.; Zhang, H.; Li, N.; Zhang, X.; Bai, D. Identification of Hotspot Regions for Candidate Genes Associated with Peanut (Arachis hypogaea L.) Pod and Seed Size on Chromosome A05. Agriculture 2024, 14, 1634. https://doi.org/10.3390/agriculture14091634

AMA Style

Zhang X, Wang L, Liu Q, Zhang X, Tian Y, Xue Y, Zhang H, Li N, Zhang X, Bai D. Identification of Hotspot Regions for Candidate Genes Associated with Peanut (Arachis hypogaea L.) Pod and Seed Size on Chromosome A05. Agriculture. 2024; 14(9):1634. https://doi.org/10.3390/agriculture14091634

Chicago/Turabian Style

Zhang, Xiaoji, Luhuan Wang, Qimei Liu, Xiaoyu Zhang, Yuexia Tian, Yunyun Xue, Huiqi Zhang, Na Li, Xin Zhang, and Dongmei Bai. 2024. "Identification of Hotspot Regions for Candidate Genes Associated with Peanut (Arachis hypogaea L.) Pod and Seed Size on Chromosome A05" Agriculture 14, no. 9: 1634. https://doi.org/10.3390/agriculture14091634

APA Style

Zhang, X., Wang, L., Liu, Q., Zhang, X., Tian, Y., Xue, Y., Zhang, H., Li, N., Zhang, X., & Bai, D. (2024). Identification of Hotspot Regions for Candidate Genes Associated with Peanut (Arachis hypogaea L.) Pod and Seed Size on Chromosome A05. Agriculture, 14(9), 1634. https://doi.org/10.3390/agriculture14091634

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop