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

Consensus genetic linkage map and QTL mapping allow to capture the genomic regions associated with agronomic traits in pearl millet

MAIN CONCLUSION

A genetic linkage map representing the pearl millet genome was constructed with SNP markers. Major and stable QTL associated with flowering, number of productive tillers, ear head length, and test weight were mapped on chromosomes 1 and 3. Pearl millet (Pennisetum glaucum) is a major cereal and fodder crop in arid and semi-arid regions of Asia and Africa. Agronomic traits are important traits in pearl millet breeding and genetic and environmental factors highly influence them. In the present study, an F9 recombinant inbred line (RIL) population derived from a cross between PT6029 and PT6129 was evaluated for agronomic traits in three environments. Utilizing a genotyping by sequencing approach, a dense genetic map with 993 single nucleotide polymorphism markers covering a total genetic distance of 1035.4 cM was constructed. The average interval between the markers was 1.04 cM, and the seven chromosomes varied from 115.39 to 206.72 cM. Quantitative trait loci (QTL) mapping revealed 35 QTL for seven agronomic traits, and they were distributed on all pearl millet chromosomes. These QTL individually explained 11.35 to 26.71% of the phenotypic variation, with LOD values ranging from 2.74 to 5.80. Notably, four QTL (qDFF1.1, qNPT3.1, qEHL3.1, and qTW1.1) associated with days to fifty percent flowering, the number of productive tillers, ear head length, and test weight were found to be major and stable QTL located on chromosomes 1 and 3. Collectively, our results provide an important base for understanding the genetic architecture of agronomic traits in pearl millet, which is useful for accelerating the genetic gain toward crop improvement.

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

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.

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.

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.

Mapping Quantitative Trait Loci (QTLs) for Hundred-Pod and Hundred-Seed Weight under Seven Environments in a Recombinant Inbred Line Population of Cultivated Peanut (Arachis hypogaea L.)

The cultivated peanut (Arachis hypogaea L.) is a significant oil and cash crop globally. Hundred-pod and -seed weight are important components for peanut yield. To unravel the genetic basis of hundred-pod weight (HPW) and hundred-seed weight (HSW), in the current study, a recombinant inbred line (RIL) population with 188 individuals was developed from a cross between JH5 (JH5, large pod and seed weight) and M130 (small pod and seed weight), and was utilized to identify QTLs for HPW and HSW. An integrated genetic linkage map was constructed by using SSR, AhTE, SRAP, TRAP and SNP markers. This map consisted of 3130 genetic markers, which were assigned to 20 chromosomes, and covered 1998.95 cM with an average distance 0.64 cM. On this basis, 31 QTLs for HPW and HSW were located on seven chromosomes, with each QTL accounting for 3.7-10.8% of phenotypic variance explained (PVE). Among these, seven QTLs were detected under multiple environments, and two major QTLs were found on B04 and B08. Notably, a QTL hotspot on chromosome A08 contained seven QTLs over a 2.74 cM genetic interval with an 0.36 Mb physical map, including 18 candidate genes. Of these, Arahy.D52S1Z, Arahy.IBM9RL, Arahy.W18Y25, Arahy.CPLC2W and Arahy.14EF4H might play a role in modulating peanut pod and seed weight. These findings could facilitate further research into the genetic mechanisms influencing pod and seed weight in cultivated peanut.

Global Transcriptome and Co-Expression Network Analyses Revealed Hub Genes Controlling Seed Size/Weight and/or Oil Content in Peanut

Cultivated peanut (Arachis hypogaea L.) is an important economic and oilseed crop worldwide, providing high-quality edible oil and high protein content. Seed size/weight and oil content are two important determinants of yield and quality in peanut breeding. To identify key regulators controlling these two traits, two peanut cultivars with contrasting phenotypes were compared to each other, one having a larger seed size and higher oil content (Zhonghua16, ZH16 for short), while the second cultivar had smaller-sized seeds and lower oil content (Zhonghua6, ZH6). Whole transcriptome analyses were performed on these two cultivars at four stages of seed development. The results showed that ~40% of the expressed genes were stage-specific in each cultivar during seed development, especially at the early stage of development. In addition, we identified a total of 5356 differentially expressed genes (DEGs) between ZH16 and ZH6 across four development stages. Weighted gene co-expression network analysis (WGCNA) based on DEGs revealed multiple hub genes with potential roles in seed size/weight and/or oil content. These hub genes were mainly involved in transcription factors (TFs), phytohormones, the ubiquitin-proteasome pathway, and fatty acid synthesis. Overall, the candidate genes and co-expression networks detected in this study could be a valuable resource for genetic breeding to improve seed yield and quality traits in peanut.

Genome-wide association study as a powerful tool for dissecting competitive traits in legumes

Legumes are extremely valuable because of their high protein content and several other nutritional components. The major challenge lies in maintaining the quantity and quality of protein and other nutritional compounds in view of climate change conditions. The global need for plant-based proteins has increased the demand for seeds with a high protein content that includes essential amino acids. Genome-wide association studies (GWAS) have evolved as a standard approach in agricultural genetics for examining such intricate characters. Recent development in machine learning methods shows promising applications for dimensionality reduction, which is a major challenge in GWAS. With the advancement in biotechnology, sequencing, and bioinformatics tools, estimation of linkage disequilibrium (LD) based associations between a genome-wide collection of single-nucleotide polymorphisms (SNPs) and desired phenotypic traits has become accessible. The markers from GWAS could be utilized for genomic selection (GS) to predict superior lines by calculating genomic estimated breeding values (GEBVs). For prediction accuracy, an assortment of statistical models could be utilized, such as ridge regression best linear unbiased prediction (rrBLUP), genomic best linear unbiased predictor (gBLUP), Bayesian, and random forest (RF). Both naturally diverse germplasm panels and family-based breeding populations can be used for association mapping based on the nature of the breeding system (inbred or outbred) in the plant species. MAGIC, MCILs, RIAILs, NAM, and ROAM are being used for association mapping in several crops. Several modifications of NAM, such as doubled haploid NAM (DH-NAM), backcross NAM (BC-NAM), and advanced backcross NAM (AB-NAM), have also been used in crops like rice, wheat, maize, barley mustard, etc. for reliable marker-trait associations (MTAs), phenotyping accuracy is equally important as genotyping. Highthroughput genotyping, phenomics, and computational techniques have advanced during the past few years, making it possible to explore such enormous datasets. Each population has unique virtues and flaws at the genomics and phenomics levels, which will be covered in more detail in this review study. The current investigation includes utilizing elite breeding lines as association mapping population, optimizing the choice of GWAS selection, population size, and hurdles in phenotyping, and statistical methods which will analyze competitive traits in legume breeding.

Genetic mapping identifies genomic regions and candidate genes for seed weight and shelling percentage in groundnut

Seed size is not only a yield-related trait but also an important measure to determine the commercial value of groundnut in the international market. For instance, small size is preferred in oil production, whereas large-sized seeds are preferred in confectioneries. In order to identify the genomic regions associated with 100-seed weight (HSW) and shelling percentage (SHP), the recombinant inbred line (RIL) population (Chico × ICGV 02251) of 352 individuals was phenotyped for three seasons and genotyped with an Axiom_Arachis array containing 58K SNPs. A genetic map with 4199 SNP loci was constructed, spanning a map distance of 2708.36 cM. QTL analysis identified six QTLs for SHP, with three consistent QTLs on chromosomes A05, A08, and B10. Similarly, for HSW, seven QTLs located on chromosomes A01, A02, A04, A10, B05, B06, and B09 were identified. BIG SEED locus and spermidine synthase candidate genes associated with seed weight were identified in the QTL region on chromosome B09. Laccase, fibre protein, lipid transfer protein, senescence-associated protein, and disease-resistant NBS-LRR proteins were identified in the QTL regions associated with shelling percentage. The associated markers for major-effect QTLs for both traits successfully distinguished between the small- and large-seeded RILs. QTLs identified for HSW and SHP can be used for developing potential selectable markers to improve the cultivars with desired seed size and shelling percentage to meet the demands of confectionery industries.

Genome-Wide Investigation of Apyrase (APY) Genes in Peanut (Arachis hypogaea L.) and Functional Characterization of a Pod-Abundant Expression Promoter AhAPY2-1p

Peanut (Arachis hypogaea L.) is an important food and feed crop worldwide and is affected by various biotic and abiotic stresses. The cellular ATP levels decrease significantly during stress as ATP molecules move to extracellular spaces, resulting in increased ROS production and cell apoptosis. Apyrases (APYs) are the nucleoside phosphatase (NPTs) superfamily members and play an important role in regulating cellular ATP levels under stress. We identified 17 APY homologs in A. hypogaea (AhAPYs), and their phylogenetic relationships, conserved motifs, putative miRNAs targeting different AhAPYs, cis-regulatory elements, etc., were studied in detail. The transcriptome expression data were used to observe the expression patterns in different tissues and under stress conditions. We found that the AhAPY2-1 gene showed abundant expression in the pericarp. As the pericarp is a key defense organ against environmental stress and promoters are the key elements regulating gene expression, we functionally characterized the AhAPY2-1 promoter for its possible use in future breeding programs. The functional characterization of AhAPY2-1P in transgenic Arabidopsis plants showed that it effectively regulated GUS gene expression in the pericarp. GUS expression was also detected in flowers of transgenic Arabidopsis plants. Overall, these results strongly suggest that APYs are an important future research subject for peanut and other crops, and AhPAY2-1P can be used to drive the resistance-related genes in a pericarp-specific manner to enhance the defensive abilities of the pericarp.

Genome-wide identification of germin-like proteins in peanut (Arachis hypogea L.) and expression analysis under different abiotic stresses

Peanut is an important food and feed crop, providing oil and protein nutrients. Germins and germin-like proteins (GLPs) are ubiquitously present in plants playing numerous roles in defense, growth and development, and different signaling pathways. However, the GLP members have not been comprehensively studied in peanut at the genome-wide scale. We carried out a genome-wide identification of the GLP genes in peanut genome. GLP members were identified comprehensively, and gene structure, genomic positions, motifs/domains distribution patterns, and phylogenetic history were studied in detail. Promoter Cis-elements, gene duplication, collinearity, miRNAs, protein-protein interactions, and expression were determined. A total of 84 GLPs (AhGLPs ) were found in the genome of cultivated peanut. These GLP genes were clustered into six groups. Segmental duplication events played a key role in the evolution of AhGLPs, and purifying selection pressure was underlying the duplication process. Most AhGLPs possessed a well-maintained gene structure and motif organization within the same group. The promoter regions of AhGLPs contained several key cis-elements responsive to 'phytohormones', 'growth and development', defense, and 'light induction'. Seven microRNAs (miRNAs) from six families were found targeting 25 AhGLPs. Gene Ontology (GO) enrichment analysis showed that AhGLPs are highly enriched in nutrient reservoir activity, aleurone grain, external encapsulating structure, multicellular organismal reproductive process, and response to acid chemicals, indicating their important biological roles. AhGLP14, AhGLP38, AhGLP54, and AhGLP76 were expressed in most tissues, while AhGLP26, AhGLP29, and AhGLP62 showed abundant expression in the pericarp. AhGLP7, AhGLP20, and AhGLP21, etc., showed specifically high expression in embryo, while AhGLP12, AhGLP18, AhGLP40, AhGLP78, and AhGLP82 were highly expressed under different hormones, water, and temperature stress. The qRT-PCR results were in accordance with the transcriptome expression data. In short, these findings provided a foundation for future functional investigations on the AhGLPs for peanut breeding programs.

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.

Genome-Wide Characterization of Ascorbate Peroxidase Gene Family in Peanut (Arachis hypogea L.) Revealed Their Crucial Role in Growth and Multiple Stress Tolerance

Ascorbate peroxidase (APX), an important antioxidant enzyme, plays a significant role in ROS scavenging by catalyzing the decrease of hydrogen peroxide under various environmental stresses. Nevertheless, information about the APX gene family and their evolutionary and functional attributes in peanut (Arachis hypogea L.) was not reported. Therefore, a comprehensive genome-wide study was performed to discover the APX genes in cultivated peanut genome. This study identified 166 AhAPX genes in the peanut genome, classified into 11 main groups. The gene duplication analysis showed that AhAPX genes had experienced segmental duplications and purifying selection pressure. Gene structure and motif investigation indicated that most of the AhAPX genes exhibited a comparatively well-preserved exon-intron pattern and motif configuration contained by the identical group. We discovered five phytohormones-, six abiotic stress-, and five growth and development-related cis-elements in the promoter regions of AhAPX. Fourteen putative ah-miRNAs from 12 families were identified, targeting 33 AhAPX genes. Furthermore, we identified 3,257 transcription factors from 38 families (including AP2, ARF, B3, bHLH, bZIP, ERF, MYB, NAC, WRKY, etc.) in 162 AhAPX genes. Gene ontology and KEGG enrichment analysis confirm the role of AhAPX genes in oxidoreductase activity, catalytic activity, cell junction, cellular response to stimulus and detoxification, biosynthesis of metabolites, and phenylpropanoid metabolism. Based on transcriptome datasets, some genes such as AhAPX4/7/17/77/82/86/130/133 and AhAPX160 showed significantly higher expression in diverse tissues/organs, i.e., flower, leaf, stem, roots, peg, testa, and cotyledon. Likewise, only a few genes, including AhAPX4/17/19/55/59/82/101/102/137 and AhAPX140, were significantly upregulated under abiotic (drought and cold), and phytohormones (ethylene, abscisic acid, paclobutrazol, brassinolide, and salicylic acid) treatments. qRT-PCR-based expression profiling presented the parallel expression trends as generated from transcriptome datasets. Our discoveries gave new visions into the evolution of APX genes and provided a base for further functional examinations of the AhAPX genes in peanut breeding programs.

Whole Genome Sequencing and Morphological Trait-Based Evaluation of UPOV Option 2 for DUS Testing in Rice

To evaluate the application potential of high-density SNPs in rice distinctness, uniformity, and stability (DUS) testing, we screened 37,929 SNP loci distributed on 12 rice chromosomes based on whole-genome resequencing of 122 rice accessions. These SNP loci were used to analyze the DUS testing of rice varieties based on the correlation between the molecular and phenotypic distances of varieties according to UPOV option 2. The results showed that statistical algorithms and the number of phenotypic traits and SNP loci all affected the correlation between the molecular and phenotypic distances of rice varieties. Relative to the other nine algorithms, the Jaccard similarity algorithm had the highest correlation of 0.6587. Both the number of SNPs and the number of phenotypes had a ceiling effect on the correlation between the molecular and phenotypic distances of varieties, and the ceiling effect of the number of SNP loci was more obvious. To overcome the correlation bottleneck, we used the genome-wide prediction method to predict 30 phenotypic traits and found that the prediction accuracy of some traits, such as the basal sheath anthocyanin color, glume length, and intensity of the green color of the leaf blade, was very low. In combination with group comparison analysis, we found that the key to overcoming the ceiling effect of correlation was to improve the resolution of traits with low predictive values. In addition, we also performed distinctness testing on rice varieties by using the molecular distance and phenotypic distance, and we found that there were large differences between the two methods, indicating that UPOV option 2 alone cannot replace the traditional phenotypic DUS testing. However, genotype and phenotype analysis together can increase the efficiency of DUS testing.

Genetic Variation in Schizothorax kozlovi Nikolsky in the Upper Reaches of the Chinese Yangtze River Based on Genotyping for Simplified Genome Sequencing

Simple Summary

Schizothorax kozlovi Nikolsky is a unique cold−water fish in the upper reaches of the Yangtze River in China and has high economic value. In our study, genetic diversity and population structure analyses were performed on seven wild populations in the upper reaches of the Yangtze River by GBS. The above results indicate that the populations of S. kozlovi have different degrees of tolerance and selection pressure in response to temperature and altitude. The Wujiang population was genetically differentiated from the Jinsha River and Yalong River populations. The Wujiang intrapopulation has greater genetic diversity and differentiation than the Jinsha River and Yalong River populations, which demonstrates that the Jinsha and Yalong populations require more attention and resources for their protection. The results of this study will increase our understanding of the diversity of S. kozlovi in the upper reaches of the Yangtze River and provide a basis for the conservation and utilization of wild resources.

Abstract

Schizothorax kozlovi Nikolsky is a unique cold−water fish in the upper reaches of the Yangtze River in China and has high economic value. In our study, genetic diversity and population structure analyses were performed on seven wild populations (originating from the Jinsha River, Yalong River, and Wujiang River) in the upper reaches of the Yangtze River by genotyping by sequencing (GBS). The results indicated that a total of 303,970 single−nucleotide polymorphisms (SNPs) were identified from the seven wild populations. Lower genetic diversity was exhibited among the intrapopulations of the three tributaries, and the Wujiang River population had significant genetic differentiation when compared to the Jinsha River and Yalong River populations. Furthermore, the selected SNPs were enriched in cellular processes, environmental adaptation, signal transduction, and related metabolic processes between the Wujiang population and the other two populations. The above results indicate that the populations of S. kozlovi have different degrees of tolerance and selection pressure in response to temperature and altitude. The Wujiang intrapopulation has greater genetic diversity and differentiation than the Jinsha River and Yalong River populations, which demonstrates that the Jinsha and Yalong populations require more attention and resources for their protection. The results of this study will increase our understanding of the diversity of S. kozlovi in the upper reaches of the Yangtze River and provide a basis for the conservation and utilization of wild resources.

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.