Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 18 March 2021
Sec. Plant Breeding

Construction of ddRADseq-Based High-Density Genetic Map and Identification of Quantitative Trait Loci for Trans-resveratrol Content in Peanut Seeds

\nHuaiyong Luo&#x;Huaiyong LuoJianbin Guo&#x;Jianbin GuoBolun YuBolun YuWeigang ChenWeigang ChenHuan ZhangHuan ZhangXiaojing ZhouXiaojing ZhouYuning ChenYuning ChenLi HuangLi HuangNian LiuNian LiuXiaoping RenXiaoping RenLiying YanLiying YanDongxin HuaiDongxin HuaiYong LeiYong LeiBoshou LiaoBoshou LiaoHuifang Jiang
Huifang Jiang*
  • Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, China

Resveratrol (trans-3,4′,5-trihydroxystilbene) is a natural stilbene phytoalexin which is also found to be good for human health. Cultivated peanut (Arachis hypogaea L.), a worldwide important legume crop, is one of the few sources of human's dietary intake of resveratrol. Although the variations of resveratrol contents among peanut varieties were observed, the variations across environments and its underlying genetic basis were poorly investigated. In this study, the resveratrol content in seeds of a recombination inbred line (RIL) population (Zhonghua 6 × Xuhua 13, 186 progenies) were quantified by high performance liquid chromatography (HPLC) method across four environments. Genotypes, environments and genotype × environment interactions significantly influenced the resveratrol contents in the RIL population. A total of 8,114 high-quality single nucleotide polymorphisms (SNPs) were identified based on double-digest restriction-site-associated DNA sequencing (ddRADseq) reads. These SNPs were clustered into bins using a reference-based method, which facilitated the construction of high-density genetic map (2,183 loci with a total length of 2,063.55 cM) and the discovery of several chromosome translocations. Through composite interval mapping (CIM), nine additive quantitative trait loci (QTL) for resveratrol contents were identified on chromosomes A01, A07, A08, B04, B05, B06, B07, and B10 with 5.07–8.19% phenotypic variations explained (PVE). Putative genes within their confidential intervals might play roles in diverse primary and secondary metabolic processes. These results laid a foundation for the further genetic dissection of resveratrol content as well as the breeding and production of high-resveratrol peanuts.

Introduction

Resveratrol is a non-flavonoid polyphenol compound synthesized and functioned as phytoalexinin more than 70 plant species (Pastor et al., 2019). Natural resveratrol exists mainly in the trans-isoform (Novelle et al., 2015) and has been identified in only a few of edible plants such as grape, peanut, and berries (Pastor et al., 2019; Singh et al., 2019). Resveratrol is synthesized by stilbene synthase (STS) via the phenylalanine/polymalonate biosynthetic pathway (Shomura et al., 2005). The dietary intake of resveratrol from moderate red wine was reported to be partly contributed to the “French paradox” in the early 1990s, and studies have been conducted to increase the contents of resveratrol in edible plants (e.g., grapes) which could be produced to health foods (Hasan and Bae, 2017). In addition, the favorable effects of resveratrol for human health have been reported in the therapeutic outcomes in many preclinical and clinical trials of cardiovascular diseases, cancer and other important human diseases, although inconsistent results were observed in other studies (Fogacci et al., 2019; Singh et al., 2019; Wang et al., 2019b).

As one of the dietary sources of resveratrol (Delaunois et al., 2009), the cultivated peanut (Arachis hypogaea L.) is an important legume crop grown in more than 100 countries (FAOSTAT, 2019). Its global productions were around 47 million tones in recent years (FAOSTAT, 2019), and its seeds are consumed worldwide in diverse forms, such as nuts, peanut butter, edible oil, and candies (Varshney et al., 2013). The resveratrol in peanut mainly exists in the trans-form, and the HPLC method has been developed to quantify its content. Significant variations were observed for resveratrol contents among different peanut varieties (Sanders et al., 2000). Using the US mini-core collection as materials, Wang et al. (2013) found that there was a nine-fold difference (30–260 μg/kg) of resveratrol content in seeds across different peanut accessions. However, this study was conducted in single environment, and the variation of resveratrol content across environments was not evaluated. In addition, the genetic basis for these variations was poorly investigated in peanut as well as other plant species producing resveratrol. Therefore, more efforts should be made in genetic studies to provide the theoretical basis and technique support for breeding and production of novel high resveratrol varieties.

QTL mapping using bi-parental genetic population has become a routine approach to dissect the genetic basis for quantitative traits, such as yield-related traits, resistance to diseases, tolerance to drought, oil content and fatty acid compositions (Varshney et al., 2009; Vishwakarma et al., 2017). However, the construction of high-quality and high-density genetic map is the prerequisite for the discovery of genomic-wide QTL. With the reducing of high-throughput sequencing cost, the SNP polymorphisms have been utilized as markers to construct high-density genetic maps (Pandey et al., 2016). Through sequencing based genotyping methods such as ddRADseq (Zhou et al., 2014), GBS (Dodia et al., 2019) and SLAF-seq (Hu et al., 2018), SNP-based genetic maps could be constructed in much shorter time than traditional gel-based genotyping and have lower cost than whole-genome resequencing (Agarwal et al., 2018). However, the previously reported high-density genetic maps were developed though de-novo method in peanut (Zhou et al., 2014; Hu et al., 2018; Dodia et al., 2019). Along with the availability of peanut genome sequences (Bertioli et al., 2019; Chen et al., 2019a; Zhuang et al., 2019), it is feasible to construct the high-density genetic map through the reference-based method approach which might benefit in making the best use of sequencing reads and discovering chromosome variations.

In the present study, the resveratrol content in peanut seeds of a RIL population were quantified across four environments. The effects of genotype, environments and genotype × environment interactions on the variation of resveratrol contents in the RIL population were analyzed. Using the recent published peanut genome sequence as reference, a high-density and high-quality genetic map was constructed through the reference-based analysis, upon which nine additive QTL were identified to be associated with resveratrol content in peanut seeds. Functions of the putative genes located in the confidential intervals of the identified QTL were investigated.

Materials and Methods

Plant Materials

The RIL population derived from the cross between Xuhua 13 and Zhonghua 6 (Luo et al., 2018) was used as plant materials in the present study. Recently, ddRADseq reads were generated for 186 RILs in F6 generation of the RIL population (Liu et al., 2020). Four generations (F6 to F9) of the RIL population were used for phenotyping resveratrol contents, and they were planted using randomized complete block design with three replications in Wuchang in 2015 (F6 generation) and 2016 (F7 generation), Yangluo in 2017 (F8 generation), and Xiangyang in 2018 (F9 generation). These trials were designated as four environments: WC2015, WC2016, YL2017, and XY2018, respectively. Matured pods were harvested, air-dried, stored at 4°C until shelled by hand before phenotyping.

Phenotyping of Resveratrol Contents by HPLC

The standard trans-resveratrol (CAS 501-36-0) was purchased from Sigma (St Louis, MO, USA) and used to generate the standard curve. The resveratrol contents in seeds of the RIL population were quantified according to the reported HPLC method (Wang et al., 2013). Briefly, 10 g of peanut seeds were ground by a coffee blender. Then, 5 ± 0.0001 g fine powder, which passed through a 20 mesh sieve, was used for extraction of resveratrol in a 50-ml centrifuge tube. The tube, into which 40 ml ethanol (85%) was added, was placed in 80°C water bath for 1 h and then centrifuged at 10,000 rpm for 10 min. Subsequently, 8 ml supernatant was filtered through a C18/AL-N SPE column (Agela Technologies). The filtrate was blow-dried with nitrogen gas and dissolved with 1 ml methanol. The sample was finally filtered with a 0.22 μm filter and quantified using the Agilent 1290 system with the C18 column (4.6 mm × 150 mm, 5 mm; Agilent Technologies, USA). The column temperature was set as 30°C. The mobile phase was consisted of (A) ultrapure water with 0.05% acetic acid and (B) HPLC-grade acetonitrile. The chromatogram was recorded at 306 nm and the flow rate was 0.9 ml/min. After 10 μl samples were injected, a 12 min elution was performed with 78% A and 22% B. After that, the column wash was washed for 2 min with 60% A and 40% B to get ready for next injection. The phenotypic distribution of resveratrol contents in the RIL population was plotted by the “ggplot2” R package.

SNP Calling and Genotyping

The ddRADseq reads of the RIL population (BioProject: PRJNA520741) (Liu et al., 2020) were re-analyzed following the pipeline of reference-based analysis (Rochette and Catchen, 2017) of the Stacks software (Catchen et al., 2013) to identify SNPs. Briefly, low-quality reads were removed by the process_radtags unit of Stacks. Then, high-quality reads from each RIL were mapped to the recently published genome sequences (the KYV3 version) of cultivated peanut cv. Tifrunner (Bertioli et al., 2019) by the BWA software (Li and Durbin, 2009). The gstacks unit of Stacks was used to build loci from uniquely mapped reads. The populations unit was used to remove loci shared by <80% of samples and output site level SNP calls in VCF format. These SNPs were transformed into genotypes using an in-house Perl script. Individual SNPs were recorded as “A” representing homologous for Zhonghua 6, “B” representing homologous for Xuhua 13, and “X” representing heterozygous or missing alleles). In order to conveniently retrieve physical positions, SNPs was designated with initial letters “TIF” (representing Tifrunner), followed by the corresponding chromosome and position. For example, TIF.01:7840610 was the SNP identified at position 7,840,610 bp on chromosome 1 (A01) of the Tifrunner reference genome.

Construction of High-Density Genetic Map

The identified high-quality SNPs were clustered into bins using a reference-based method, i.e., the SNPs with identical genotypes and alongside each other on the reference genome were clustered as one bin, and the SNP with the least missing genotypes was selected to represent the bin. The genetic map was constructed based on these bins using the QTL IciMapping software (Meng et al., 2015). Pearson's Chi square test was performed to test the goodness of fit to the expected segregation ratio 1:1 for each locus (P < 0.01).

Identification of QTL for Resveratrol Content

Additive QTL were identified by the composite interval mapping (CIM) method of the WinQTLCart software (Wang et al., 2012) with default parameters for each environment. The BLUP values for resveratrol contents across the four environments were calculated using the “lme4” R package and used to identify QTL as well. The identified QTL were designated with initial letters “q,” followed by trait name (RES for resveratrol content) and the corresponding chromosome. A number were added if two or more QTL were identified on the same chromosome. For example, qRESB07.1and qRESB07.2 were the first and second QTL identified on chromosome B07, respectively. The putative genes within the 2 LOD confidential interval of identified QTL were retrieved from the genomic annotation of the peanut accession Tifrunner (Bertioli et al., 2019). The annotations of the identified genes were obtained from the Tifrunner annotation file. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the putative genes were conducted using the Blast2GO suite (Gotz et al., 2008).

Results

Phenotypic Variation of Resveratrol Contents in the RIL Population

The resveratrol content of the variety Xuhua 13 (the female parent) was consistently lower than that of the variety Zhonghua 6 (the male parent) across the four environments including WC2015, WC2016, YL2017 and XY2018 (Figure 1, Table 1). The resveratrol contents of the 186 RILs varied from 37.33 to 270.00 μg/kg in WC2015, from 3.61 to 258.79 μg/kg in WC2016, from 8.96 to 282.89 μg/kg in YL2017, and from 13.60 to 215.06 μg/kg in XY2018 (Table 1). As shown in Figure 1 and Table 1, the distribution of resveratrol contents in the RIL population was continuous and with transgressive segregation. According to the variance analysis across the four environments, genotypes, environments and genotype × environment interactions significantly influenced resveratrol contents among RILs (Supplementary Table 1, P < 0.001). The broad sense heritability of resveratrol content in the present study was estimated to be 0.33, indicating that environments had significant influence in resveratrol contents of peanut seeds.

FIGURE 1
www.frontiersin.org

Figure 1. Phenotypic distribution of resveratrol contents in the RIL population across environments. The y-axis represented density, while the x-axis represented values of resveratrol contents (μg/kg).

TABLE 1
www.frontiersin.org

Table 1. Descriptive statistics of resveratrol contents in the RIL population across four environments.

Genotyping of the RIL Population

The 359.96 Gb ddRADseq clean data (3,712.28 M reads) of the parents and RIL population were analyzed using the recently reported genome sequences of Tifrunner as reference. A total of 2,616,605 loci were built following the pipeline of reference-based analysis of the stacks software (Figure 2). Among the RILs, the loci ranged from 492,063 to 757,190 and the sequencing depths varied from 4.17 to 13.45. After removing loci missed in more than 20% RILs, 549,163 loci were used to call SNPs. In total, 9,353 SNPs across the RIL population were obtained. A total of 8,114 SNPs showed homozygous genotypes in two parents of the RIL population. These 8,114 SNPs were clusters as 2,191 bins, among which 1,277 bins had single SNP while the other 914 bins had 2-148 adjacent SNPs (7.48 on average) with same genotypes in the RIL population.

FIGURE 2
www.frontiersin.org

Figure 2. Genotyping of the RIL population based on ddRADseq data.

Construction of High-Density Genetic Map

The above 2,191 bins were used to construct genetic map, and the final genetic map was consisted of 2,183 genetic loci (Figure 3, Supplementary Figure 1, Supplementary Table 2). The total length of the 20 linkage groups (LGs) was 2,063.55 cM. The 10 LGs (A01–A10) of subgenome A contained 1,050 loci with a total length of 1,004.73 cM while the 10 LGs (B01–B10) of subgenome B contained 1,133 loci spanning 1,058.82 cM (Table 2). LG B03 had the largest loci (191) and LG B09 had the least loci (25). The length of LGs ranged from 61.08 cM for A04 to 128.69 cM for A03. The average marker interval of the genetic map was 0.95 cM, ranging from 0.52 cM for B04 to 4.51 cM for B09. Half of the LGs had marker densities of <1 cM/locus, and 98.21% of the genetic distances between markers were lower than 5 cM, indicating a high degree of linkage. The largest gap between markers was 39.54 cM on LG B07, which was corresponding to a large (>60 Mb) physical gap without polymorphic SNPs between parents of the RIL population. According to the Pearson's Chi square test, segregations of 92 loci (accounting for 4.21% of the total loci) were distorted (P < 0.01, Supplementary Table 2). The 92 loci distortedly segregated to Zhonghua 6 in the A subgenome while to Xuhua 13 in the B subgenome. LG B05 contained the most of distorted loci (41), followed by A05 (29), B07 (5), A06 (3), A09 (3), B01 (3), A07 (2), B04 (2), B02 (1), B03 (1), B08 (1) and B10 (1).

FIGURE 3
www.frontiersin.org

Figure 3. Graph presentation of the constructed high-density genetic map.

TABLE 2
www.frontiersin.org

Table 2. Description of the high-density genetic linkage map constructed in the present study.

Evaluation of Genome Synteny of the High-Density Genetic Map

A total of 2,099 SNP loci (accounting for 96.15%) of the genetic map were located on same chromosomes of the reference genome of Tifrunner (Bertioli et al., 2019), and good collinearity was clearly observed for almost all chromosomes (Figure 4, Supplementary Figure 2). BLASTn analysis of the RAD-tags of the other 84 SNP loci against to genome sequences of cultivated peanut Tifrunner, Shitouqi (Zhuang et al., 2019), and Fuhuasheng (Chen et al., 2019a) as well as wild species (Bertioli et al., 2016) found that translocations occurred in different peanut varieties (Supplementary Table 3). Obviously, reciprocal translocations occurred between chromosome A03 and B03 as well as A06 and B06 in Tifrunner but not in Shitouqi, Fuhuasheng and the parents of the RIL population. Interestingly, a fragment from B10 was translocated to B03 in Tifrunner and Shitouqi but not in Fuhuasheng.

FIGURE 4
www.frontiersin.org

Figure 4. Circos graph presentation of the genome synteny between genetic map and physical map.

QTL Associated With Resveratrol Content

Based on the high-density genetic map, three, two, and two QTL for resveratrol content were identified the phenotyping data of the WC2016, WC2017, and XY2018 environments, respectively, while no QTL was identified in the YL2017 environment. In addition, three QTL were identified using the BLUP values across the four environments (Table 3). The qRESB04 could be repeatly identified using the phenotyping data of the WC2016 environment as well as the BLUP values across four environments. Therefore, a total of nine QTL were identified across four environments (Table 3), and were located on eight chromosomes, namely A01, A07, A08, B04, B05, B06, B07, and B10 (Figure 5), with 5.07–8.19% phenotypic variation explained (PVE). The confidential intervals of the identified QTL varied from 2.0 cM for qRESB05 to 30.9 cM for qRESA07 and averaged 10.64 cM. The physical intervals of the identified QTL in the reference genomes of Tifrunner were conveniently estimated according to the names of flanking markers. The qRESB05 and qRESB10 were located in lower recombination regions on chromosome B05 and B10, and their physical intervals (~18.53 Mb and ~22.32 Mb, respectively) were much larger than others (0.69–6.68 Mb).

TABLE 3
www.frontiersin.org

Table 3. Summary of QTL identified for resveratrol contents in the RIL populations across four environments.

FIGURE 5
www.frontiersin.org

Figure 5. The genetic map locations of QTL identified for resveratrol contents. QTL and their 1 LOD confidential intervals (filled boxes) and 2 LOD confidential intervals (lines) were colored according to environments as shown in figure legend. The purple color was used to highlight the 2 LOD confidential intervals on chromosome bars. Bins within the 2 LOD confidential intervals were highlighted in purple color as well.

A total of 2,429 putative genes were retrieved within physical intervals of these QTL from the genome annotation of Tifrunner, ranging from 43 genes for qRESB06 to 544 genes for qRESA07 (Table 3, Supplementary Table 4). As shown in Supplementary Table 5, GO terms were retrieved for 1,410 genes from the Tifrunner genome annotation file. In terms of cellular components, cell (14.8%) was the most frequent GO term, followed by membrane (13.6%), organelle (9.6%), and protein-containing complex (6.7%) (Supplementary Figure 3). For molecular functions, binding (63.2%) was the most frequent GO term and followed by catalytic activity (47.4%) (Supplementary Figure 3). For biological processes, metabolic process (44.4%) was the most frequent GO term, followed by cellular process (35.0%), biological regulation (10.4%), and localization (9.5%) (Supplementary Figure 3). Based on the KEGG analysis, 94 enzymes that might be encoded by 161 genes were assigned to 81 KEGG pathways of various primary and secondary metabolisms (Supplementary Table 6).

Discussion

Peanut is one of the few edible plants for human dietary intake of resveratrol. Sanders et al. (2000) reported that the resveratrol contents of 15 peanut varieties were ~196.47 μg/kg. Wang et al. (2013) detected the seeds of 102 accessions within the U.S. peanut mini-core collection and found that the resveratrol contents varied from 30 to 260 μg/kg and averaged at 100 μg/kg. In the present study, the resveratrol content in seeds of the RIL population showed similar variation range (i.e. 3.61–282.69 μg/kg). In addition, significant variations across environments were observed, including locations and years (Table 1 and Supplementary Table 1), which is consistent with the reports that varieties, geographical indications and vintage years significantly influenced the resveratrol contents in red wines (Nikfardjam et al., 2006; Pastor et al., 2019). Moreover, the broad sense heritability of resveratrol content across environments was relatively low (0.33). Therefore, attentions should be paid to not only varieties but also producing geographical locations as well as years in the production of high resveratrol peanut, and the breeders should specify their breeding programs for different target producing locations.

The availability of multiple peanut reference genomes facilitated the construction and evaluation of high-density genetic map as well as the discovery of chromosome variations. Previous studies constructed high-density genetic map through the de-novo approach of the stacks software (Zhou et al., 2014; Wang et al., 2018; Liu et al., 2020). However, majority of the linkage map construction software are not affordable for the large number of SNPs. The SNPs were reduced randomly or with very stringent filters (Zhou et al., 2014; Wang et al., 2018; Liu et al., 2020), which might cause the loss of useful SNPs. In the present study, using the recently reported cultivate peanut genome assembly (Bertioli et al., 2019) as reference, the previously reported ddRADseq reads (Liu et al., 2020) were re-analyzed through the reference-based pipeline of stacks. The identified SNPs were binned according to not only their genotypes in the RIL population but also physical positions on the reference genome, which effectively removed redundant SNPs and helped in the construction of high-density genetic map (Figures 2, 3). The constructed genetic map was estimated to cover 96.62% of the reference genome, and its loci number and density reached a fairly high level when compared to previous studies (Ravi et al., 2011; Qin et al., 2012; Huang et al., 2015; Chen et al., 2016, 2017, 2019b; Wang et al., 2019a). Moreover, good collinearity between the constructed genetic map with the physical map was revealed by comparing to multiple peanut reference genomes (Figure 4 and Supplementary Table 3). Notably, translocations between chromosomes, such as A03 and B03, A06 and B06, B10, and B03, were observed in present study, which was ignored in the previous study using the de-novo approach (Liu et al., 2020).

The resveratrol content is a complicated trait and controlled by multiple QTL from both subgenomes of peanut. In the present study, a total of nine additive QTL (5.07–8.19% PVE) for resveratrol content were identified on three chromosomes (A01, A07, and A08) of the A subgenome and five chromosomes (B04, B05, B06, B07, and B10) of the B subgenome (Table 3, Figure 5). The top five RILs had higher resveratrol contents than parents, and they possess the favored alleles of 4–7 identified QTL from both parents (Supplementary Table 7). Only one QTL, the qRESB04, was identified in the WC2016 environment as well as using the BLUP values across four environments, while the other nine QTL were identified only in single environment. In addition, no QTL were identified in the YL2017 environment, in which majority of RILs skewed obviously toward low resverartrol contents (Figure 1). These phenomena were consistent with the fact that environments and QTL × environment interactions had significant influences in peanut resveratrol content. A pool of 2,429 putative genes were retrieved from the corresponding physical intervals of these QTL in the Tifrunner genome (Bertioli et al., 2019), and many of them involved in diverse primary and secondary metabolic pathways. However, target genes influencing resveratrol content could not be predicted in the present study. Further studies, e.g., transcriptome/metabolome analysis, might provide more evidence to identify key metabolites and candidate genes associated with resveratrol synthesis in peanut. Moreover, the peanut synthase of resveratrol (STS) (Shomura et al., 2005) might be located at ~13 Mb on chromosome B04 or ~11 Mb on A04, which were outside of the identified QTL. Therefore, resveratrol, as a secondary metabolite, might be synthesized by STS under the influence of various upstream biological processes, and it would be a challenge to breed high resveratrol peanut varieties.

In conclusion, the present study constructed a high quality genetic map and identified nine QTL with main effects for resveratrol content using a RIL population across four environments. These QTL lay a foundation for the future marker assisted selection (MAS) in the improvement of resveratrol content. In future studies, a wide range screening of germplasms in multiple producing environments should be conducted to determine the range of resveratrol contents in peanut seeds and identified elite accessions for genetic studies and breeding utilization. In addition, more studies could be conducted to evaluate the resveratrol contents, variations as well as their relationships in different peanut tissues, which might help developing diverse high resveratrol peanut products. Along with the reducing cost of high-throughput sequencing or genotyping, genome-based trait prediction models (Pandey et al., 2020) for the resveratrol content could be build and might be used in genomic selection (Crossa et al., 2017) to breed high resveratrol varieties.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

HL, YL, BL, and HJ conceived, designed, and supervised the experiments. XR, LY, DH, YL, BL, and HJ developed the RIL population. JG, BY, WC, HZ, XZ, YC, NL, and LH conducted field trials and phenotyping of resveratrol contents by HPLC. HL performed the construction of high-density genetic map from sequencing data and interpreted the results. HL prepared the first draft and YL, BL, and HJ contributed to the final editing of manuscript. All authors read and approved the final manuscript.

Funding

This project was supported by the National Natural Science Foundations of China (Grant Numbers 31971903, 31601340, and 31761143005) and the Central Public-interest Scientific Institution Basal Research Fund (Grant Number 1610172019008). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer MY declared a past co-authorship with several of the authors BL and HJ to the handling editor.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.644402/full#supplementary-material

References

Agarwal, G., Clevenger, J., Pandey, M. K., Wang, H., Shasidhar, Y., Chu, Y., et al. (2018). High-density genetic map using whole-genome resequencing for fine mapping and candidate gene discovery for disease resistance in peanut. Plant Biotechnol. J 16, 1954–1967. doi: 10.1111/pbi.12930

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertioli, D. J., Cannon, S. B., Froenicke, L., Huang, G., Farmer, A. D., Cannon, E. K., et al. (2016). The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat. Genet. 48, 438–446. doi: 10.1038/ng.3517

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertioli, D. J., Jenkins, J., Clevenger, J., Dudchenko, O., Gao, D., Seijo, G., et al. (2019). The genome sequence of segmental allotetraploid peanut Arachis hypogaea. Nat. Genet. 51, 877–884. doi: 10.1038/s41588-019-0405-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Catchen, J., Hohenlohe, P. A., Bassham, S., Amores, A., and Cresko, W. A. (2013). Stacks: an analysis tool set for population genomics. Mol. Ecol. 22, 3124–3140. doi: 10.1111/mec.12354

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Jiao, Y., Cheng, L., Huang, L., Liao, B., Tang, M., et al. (2016). Quantitative trait locus analysis for pod- and kernel-related traits in the cultivated peanut (Arachis hypogaea L.). BMC Genet 17, 25. doi: 10.1186/s12863-016-0337-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Lu, Q., Liu, H., Zhang, J., Hong, Y., Lan, H., et al. (2019a). Sequencing of cultivated peanut, Arachis hypogaea, yields insights into genome evolution and oil improvement. Mol. Plant. 12, 920–934. doi: 10.1016/j.molp.2019.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Ren, X., Zheng, Y., Zhou, X., Huang, L., Yan, L., et al. (2017). Genetic mapping of yield traits using RIL population derived from Fuchuan Dahuasheng and ICG6375 of peanut (Arachis hypogaea L.). Mol. Breeding 37:17. doi: 10.1007/s11032-016-0587-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Wang, Z., Ren, X., Huang, L., Guo, J., Zhao, J., et al. (2019b). Identification of major QTL for seed number per pod on chromosome A05 of tetraploid peanut (Arachis hypogaea L.). Crop J 7, 238–248. doi: 10.1016/j.cj.2018.09.002

CrossRef Full Text | Google Scholar

Crossa, J., Perez-Rodriguez, P., Cuevas, J., Montesinos-Lopez, O., Jarquin, D., de los Campos, G., et al. (2017). Genomic selection in plant breeding: methods, models, and perspectives. Trends Plant Sci. 22, 961–975. doi: 10.1016/j.tplants.2017.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Delaunois, B., Cordelier, S., Conreux, A., Clement, C., and Jeandet, P. (2009). Molecular engineering of resveratrol in plants. Plant. Biotechnol. J 7, 2–12. doi: 10.1111/j.1467-7652.2008.00377.x

CrossRef Full Text | Google Scholar

Dodia, S. M., Joshi, B., Gangurde, S. S., Thirumalaisamy, P. P., Mishra, G. P., Narandrakumar, D., et al. (2019). Genotyping-by-sequencing based genetic mapping reveals large number of epistatic interactions for stem rot resistance in groundnut. Theor. Appl. Genet. 132, 1001–1016. doi: 10.1007/s00122-018-3255-7

PubMed Abstract | CrossRef Full Text | Google Scholar

FAOSTAT (2019). FAOSTAT [Online]. Rome: FAO. Available online at: http://www.fao.org/faostat/en/#data/QC (accessed July 10, 2019).

Google Scholar

Fogacci, F., Tocci, G., Presta, V., Fratter, A., Borghi, C., and Cicero, A. F. G. (2019). Effect of resveratrol on blood pressure: a systematic review and meta-analysis of randomized, controlled, clinical trials. Crit. Rev. Food Sci. 59, 1605–1618. doi: 10.1080/10408398.2017.1422480

PubMed Abstract | CrossRef Full Text | Google Scholar

Gotz, S., Garcia-Gomez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435. doi: 10.1093/nar/gkn176

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasan, M., and Bae, H. (2017). An overview of stress-induced resveratrol synthesis in grapes: perspectives for resveratrol-enriched grape products. Molecules 22:294. doi: 10.3390/molecules22020294

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, X. H., Zhang, S. Z., Miao, H. R., Cui, F. G., Shen, Y., Yang, W. Q., et al. (2018). High-density genetic map construction and identification of QTLs controlling oleic and linoleic acid in peanut using SLAF-seq and SSRs. Sci. Rep. 8:5479. doi: 10.1038/s41598-018-23873-7

CrossRef Full Text | Google Scholar

Huang, L., He, H. Y., Chen, W. G., Ren, X. P., Chen, Y. N., Zhou, X. J., et al. (2015). Quantitative trait locus analysis of agronomic and quality-related traits in cultivated peanut (Arachis hypogaea L.). Theor. Appl. Genet. 128, 1103–1115. doi: 10.1007/s00122-015-2493-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760. doi: 10.1093/bioinformatics/btp324

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, N., Guo, J., Zhou, X., Wu, B., Huang, L., Luo, H., et al. (2020). High-resolution mapping of a major and consensus quantitative trait locus for oil content to a ~ 0.8-Mb region on chromosome A08 in peanut (Arachis hypogaea L.). Theor. Appl. Genet. 133, 37–49. doi: 10.1007/s00122-019-03438-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, H., Guo, J., Ren, X., Chen, W., Huang, L., Zhou, X., et al. (2018). Chromosomes A07 and A05 associated with stable and major QTLs for pod weight and size in cultivated peanut (Arachis hypogaea L.). Theor. Appl. Genet. 131, 267–282. doi: 10.1007/s00122-017-3000-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, L., Li, H. H., Zhang, L. Y., and Wang, J. K. (2015). QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop. J. 3, 269–283. doi: 10.1016/j.cj.2015.01.001

CrossRef Full Text | Google Scholar

Nikfardjam, M. S. P., Mark, L., Avar, P., Figler, M., and Ohmacht, R. (2006). Polyphenols, anthocyanins, and trans-resveratrol in red wines from the Hungarian Villany region. Food Chem. 98, 453–462. doi: 10.1016/j.foodchem.2005.06.014

CrossRef Full Text | Google Scholar

Novelle, M. G., Wahl, D., Dieguez, C., Bernier, M., and de Cabo, R. (2015). Resveratrol supplementation: Where are we now and where should we go? Ageing Res. Rev. 21, 1–15. doi: 10.1016/j.arr.2015.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandey, M. K., Chaudhari, S., Jarquin, D., Janila, P., Crossa, J., Patil, S. C., et al. (2020). Genome-based trait prediction in multi- environment breeding trials in groundnut. Theor. Appl. Genet. 133, 3101–3117. doi: 10.1007/s00122-020-03658-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandey, M. K., Roorkiwal, M., Singh, V. K., Ramalingam, A., Kudapa, H., Thudi, M., et al. (2016). Emerging genomic tools for legume breeding: current status and future prospects. Front. Plant. Sci. 7:455. doi: 10.3389/fpls.2016.00455

PubMed Abstract | CrossRef Full Text | Google Scholar

Pastor, R. F., Restani, P., Di Lorenzo, C., Orgiu, F., Teissedre, P. L., Stockley, C., et al. (2019). Resveratrol, human health and winemaking perspectives. Crit. Rev. Food Sci. 59, 1237–1255. doi: 10.1080/10408398.2017.1400517

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, H., Feng, S., Chen, C., Guo, Y., Knapp, S., Culbreath, A., et al. (2012). An integrated genetic linkage map of cultivated peanut (Arachis hypogaea L.) constructed from two RIL populations. Theor. Appl. Genet. 124, 653–664. doi: 10.1007/s00122-011-1737-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravi, K., Vadez, V., Isobe, S., Mir, R. R., Guo, Y., Nigam, S. N., et al. (2011). Identification of several small main-effect QTLs and a large number of epistatic QTLs for drought tolerance related traits in groundnut (Arachis hypogaea L.). Theor. Appl. Genet. 122, 1119–1132. doi: 10.1007/s00122-010-1517-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rochette, N. C., and Catchen, J. M. (2017). Deriving genotypes from RAD-seq short-read data using Stacks. Nat. Protoc. 12, 2640–2659. doi: 10.1038/nprot.2017.123

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanders, T. H., McMichael, R. W. Jr, and Hendrix, K. W. (2000). Occurrence of resveratrol in edible peanuts. J. Agric. Food. Chem. 48, 1243–1246. doi: 10.1021/jf990737b

PubMed Abstract | CrossRef Full Text | Google Scholar

Shomura, Y., Torayama, I., Suh, D. Y., Xiang, T., Kita, A., Sankawa, U., et al. (2005). Crystal structure of stilbene synthase from Arachis hypogaea. Proteins 60, 803–806. doi: 10.1002/prot.20584

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, A. P., Singh, R., Verma, S. S., Rai, V., Kaschula, C. H., Maiti, P., et al. (2019). Health benefits of resveratrol: evidence from clinical studies. Med. Res. Rev. 39, 1851–1891. doi: 10.1002/med.21565

PubMed Abstract | CrossRef Full Text | Google Scholar

Varshney, R. K., Bertioli, D. J., Moretzsohn, M. C., Vadez, V., Krishnamurthy, L., Aruna, R., et al. (2009). The first SSR-based genetic linkage map for cultivated groundnut (Arachis hypogaea L.). Theor. Appl. Genet. 118, 729–739. doi: 10.1007/s00122-008-0933-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Varshney, R. K., Mohan, S. M., Gaur, P. M., Gangarao, N. V., Pandey, M. K., Bohra, A., et al. (2013). Achievements and prospects of genomics-assisted breeding in three legume crops of the semi-arid tropics. Biotechnol. Adv. 31, 1120–1134. doi: 10.1016/j.biotechadv.2013.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Vishwakarma, M. K., Nayak, S. N., Guo, B., Wan, L., Liao, B., Varshney, R. K., et al. (2017). “Classical and molecular approaches for mapping of genes and quantitative trait loci in peanut,” in The Peanut Genome, eds R. K. Varshney, M. K. Pandey and N. Puppala (Cham: Springer International Publishing), 93–116. doi: 10.1007/978-3-319-63935-2_7

CrossRef Full Text | Google Scholar

Wang, L., Yang, X., Cui, S., Mu, G., Sun, X., Liu, L., et al. (2019a). QTL mapping and QTL × environment interaction analysis of multi-seed pod in cultivated peanut (Arachis hypogaea L.). Crop. J. 7, 249–260. doi: 10.1016/j.cj.2018.11.007

CrossRef Full Text | Google Scholar

Wang, L. F., Zhou, X. J., Ren, X. P., Huang, L., Luo, H. Y., Chen, Y. N., et al. (2018). A major and stable QTL for bacteria wilt resistance on chromosome B02 identified using a high-density SNP-based genetic linkage map in cultivated peanut Yuanza 9102 derived population. Front. Genet. 9:652. doi: 10.3389/fgene.2018.00652

CrossRef Full Text | Google Scholar

Wang, M., Jiang, S., Yu, F., Zhou, L., and Wang, K. (2019b). Noncoding RNAs as molecular targets of resveratrol underlying its anticancer effects. J. Agric. Food Chem. 67, 4709–4719. doi: 10.1021/acs.jafc.9b01667

CrossRef Full Text | Google Scholar

Wang, M. L., Chen, C. Y., Tonnis, B., Barkley, N. A., Pinnow, D. L., Pittman, R. N., et al. (2013). Oil, fatty acid, flavonoid, and resveratrol content variability and FAD2A functional SNP genotypes in the US peanut mini-core collection. J. Agric. Food. Chem. 61, 2875–2882. doi: 10.1021/jf305208e

CrossRef Full Text | Google Scholar

Wang, S., Basten, C. J., and Zeng, Z.-B. (2012). Windows QTL Cartographer 2.5 [Online]. Department of Statistics, North Carolina State University, Raleigh, NC. Available online at: http://statgen.ncsu.edu/qtlcart/WQTLCart.htm (accessed December 1, 2019).

Google Scholar

Zhou, X. J., Xia, Y. L., Ren, X. P., Chen, Y. L., Huang, L., Huang, S. M., et al. (2014). Construction of a SNP-based genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq). BMC Genomics 15:351. doi: 10.1186/1471-2164-15-351

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhuang, W., Chen, H., Yang, M., Wang, J., Pandey, M. K., Zhang, C., et al. (2019). The genome of cultivated peanut provides insight into legume karyotypes, polyploid evolution and crop domestication. Nat. Genet. 51, 865–876. doi: 10.1038/s41588-019-0402-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: peanut, resveratrol, high-density genetic map, QTLs, single nucleotide polymorphisms

Citation: Luo H, Guo J, Yu B, Chen W, Zhang H, Zhou X, Chen Y, Huang L, Liu N, Ren X, Yan L, Huai D, Lei Y, Liao B and Jiang H (2021) Construction of ddRADseq-Based High-Density Genetic Map and Identification of Quantitative Trait Loci for Trans-resveratrol Content in Peanut Seeds. Front. Plant Sci. 12:644402. doi: 10.3389/fpls.2021.644402

Received: 22 December 2020; Accepted: 25 February 2021;
Published: 18 March 2021.

Edited by:

Hanwei Mei, Shanghai Agrobiological Gene Center, China

Reviewed by:

Mei Yuan, Shandong Peanut Research Institute (CAAS), China
Xiaodong Wang, Jiangsu Academy of Agricultural Sciences (JAAS), China

Copyright © 2021 Luo, Guo, Yu, Chen, Zhang, Zhou, Chen, Huang, Liu, Ren, Yan, Huai, Lei, Liao and Jiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Huifang Jiang, peanutlab@oilcrops.cn

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.