Front. Plant Sci.Frontiers in Plant ScienceFront. Plant Sci.1664-462XFrontiers Media S.A.10.3389/fpls.2021.745408Plant ScienceOriginal ResearchHigh-Density Genetic Linkage Map Construction Using Whole-Genome Resequencing for Mapping QTLs of Resistance to Aspergillus flavus Infection in PeanutJiangYifei†LuoHuaiyong†YuBolunDingYingbinKangYanpingHuangLiZhouXiaojingLiuNianChenWeigangGuoJianbinHuaiDongxinLeiYongJiangHuifangYanLiyingLiaoBoshou*Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences (CAAS), Wuhan, China
Edited by: Byoung-Cheorl Kang, Seoul National University, South Korea
Reviewed by: Gyan P. Mishra, Indian Agricultural Research Institute (ICAR), India; William Underwood, Agricultural Research Service, United States Department of Agriculture (USDA), United States
*Correspondence: Boshou Liao liaoboshou@caas.cn
This article was submitted to Plant Breeding, a section of the journal Frontiers in Plant Science
†These authors have contributed equally to this work and share first authorship
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.
The cultivated peanut (Arachis hypogaea L.), which is rich in edible oil and protein, is widely planted around the world as an oil and cash crop. However, aflatoxin contamination seriously affects the quality safety of peanuts, hindering the development of the peanut industry and threatening the health of consumers. Breeding peanut varieties with resistance to Aspergillus flavus infection is important for the control of aflatoxin contamination, and understanding the genetic basis of resistance is vital to its genetic enhancement. In this study, we reported the quantitative trait locus (QTL) mapping of resistance to A. flavus infection of a well-known resistant variety, J11. A mapping population consisting of 200 recombinant inbred lines (RILs) was constructed by crossing a susceptible variety, Zhonghua 16, with J11. Through whole-genome resequencing, a genetic linkage map was constructed with 2,802 recombination bins and an average inter-bin distance of 0.58 cM. Combined with phenotypic data of an infection index in 4 consecutive years, six novel resistant QTLs with 5.03–10.87% phenotypic variances explained (PVE) were identified on chromosomes A05, A08, B01, B03, and B10. The favorable alleles of five QTLs were from J11, while that of one QTL was from Zhonghua 16. The combination of these favorable alleles significantly improved resistance to A. flavus infection. These results could contribute greatly to the understanding of the genetic basis of A. flavus resistance and could be meaningful in the improvement of further resistance in peanuts.
The cultivated peanut, also known as groundnut, is one of the most important oil crops and is widely cultivated in most tropical and sub-tropical regions. China, India, the USA, and Nigeria account for more than 90% of the cultivated area of this crop. In China, peanuts (Arachis hypogaea L.) are grown on 4.5 million ha with a total production of 17.52 million t (FAOSTAT, 2019).
Aflatoxin contamination, which is caused by a pathogen, namely, Aspergillus flavus, is a huge threat to the peanut industry and other important crops such as corn, rice, and wheat (Mahato et al., 2019; Womack et al., 2020). Contaminated agricultural products could endanger the health of human beings and animals (Pittet, 1998; Desjardins, 2003; Kew, 2013). Peanut pods and seeds could be infected in the soil during pre-harvest or during drying, storage, and transport stages after being harvested (Torres et al., 2014). Comprehensive prevention strategies including using bio-control agents, having good agricultural practices, and planting resistant varieties were suggested to control aflatoxin contamination in peanuts (Torres et al., 2014). The breeding of varieties with resistance to A. flavus infection has been one of the main objectives in peanut breeding programs and was recognized as the most cost-effective measure to solve the problem (Wang et al., 2016). However, only a few germplasms were found to possess the resistance, and the genetic bases of their resistance remain unclear. The breeding of novel varieties faces a tough challenge in peanuts. Thus, understanding the genetic pattern of resistance to A. flavus infection is vital to the enhancement of resistance in peanuts.
The Indian commercial peanut variety J11 is a well-known resistant germplasm to A. flavus infection, while the underlying genetic bases of its resistance remain largely unclear. It was first reported by Mehan and McDonald (1980) and confirmed in different environments by Mehan et al. (1981) and Nigam et al. (2009). J11 was also reported to be resistant to A. parasiticus infection (Kisyombe et al., 1985). The infection resistance of J11 was reported to be related to drought stress, pod maturity (Mehan et al., 1986), and seed coat integrity (Asis et al., 2005). Nayak et al. (2017) deployed an RNA-seq approach to understand the host–pathogen interaction, and they found 4,445 differentially expressed genes (DEGs) at four critical stages after inoculation in J11 and JL 24 (a susceptible genotype). Zhao et al. (2019) conducted both transcriptomic and proteomic analyses to reveal the changes that occurred during the infection of J11 by A. flavus, and 663 DEGs and 314 proteins were identified. However, the key genetic loci or genes responsible for the infection resistance of J11 were still unclear.
In recent years, quantitative trait locus (QTL) mapping was successfully used to identify QTLs with resistance to A. flavus infection from two newly identified resistant germplasms named ICG12625 (Yu et al., 2019) and Xinhuixiaoli (Khan et al., 2020). ICG12625 is a landrace from Ecuador, and Yu et al. (2019) constructed a recombinant inbred line (RIL) population from the cross between ICG12625 and a susceptible variety, Zhonghua 10. A genetic map with 1,219 SSR loci was constructed, and two QTLs were identified on chromosomes A03 and A10 with 7.96 and 12.16% phenotypic variation explained (PVE), respectively. On the other hand, Xinhuixiaoli is a Chinese landrace with resistance to A. flavus infection, and it was crossbred with a susceptible variety, Yueyou 92, to generate an RIL population (Khan et al., 2020). A single nucleotide polymorphism (SNP)-based genetic map was constructed by specific length-amplified fragment sequencing (SLAF-seq), and a major QTL on A03 with 18.11% PVE and a minor QTL with 4.4% PVE on B04 were identified.
In order to reveal the underlying genetic basis of the durable resistance to A. flavus infection of the well-known variety J11 at the QTL level, a susceptible variety, Zhonghua 16, was crossbred with J11 to develop an RIL population, upon which a high-quality genetic map was constructed by whole-genome resequencing. Additionally, QTLs associated with resistance to A. flavus infection were identified in four environments in the present study.
Materials and MethodsPlant Materials
Zhonghua 16 is a variety developed by the Oil Crops Research Institute of the Chinese Academy of Agricultural Science (OCRI-CAAS), Wuhan, China, and is susceptible to A. flavus infection. J11 is a variety introduced from the International Crop Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad, India, and is resistant to A. flavus infection. Zhonghua 16 (the female parent) was crossbred with J11 (the male parent) to develop a mapping population containing 200 recombinant RILs using the single seed descent (SSD) method. The RIL population (F8) was used for genotyping and genetic map construction. Four generations (F7-F10) were planted in the experimental station of the OCRI-CAAS in Wuhan, China in 4 consecutive years (2017, 2018, 2019, and 2020, respectively). The field trials were conducted using a randomized complete block design with three replications following the method reported by Yu et al. (2019). These trials were designated as four environments, i.e., WH2017, WH2018, WH2019, and WH2020, respectively. After harvesting, the peanut pods were dried immediately until the moisture content of the seeds was 5–8%. Healthy and mature peanut seeds were artificially selected for further analysis.
Phenotyping of Peanut Seed Resistance to <italic>Aspergillus flavus</italic> Infection and Data Analysis
The percent seed infection index (PSII) for each RIL and parents to the toxigenic A. flavus strain AF2202 was identified according to a previously reported method (Yu et al., 2019). Briefly, about 20 peanut seeds from each line were sterilized with 75% ethanol for 3 min. Then, they were rinsed with sterile distilled water three times and inoculated with 1 ml of a conidial suspension (2 × 106 conidia/ml) of A. flavus. These seeds were cultured at 29 ± 1°C in the dark for 7 days to investigate resistance. Then, the resistance level to A. flavus infection of each seed was divided into four different grades, namely, (1) grade 0: seed surface without spore coverage; (2) grade 1: seed surface has <1/3 spore coverage; (3) grade 2: seed surface covered 1/3–2/3 spore; (4) grade 3: seed surface exceeds 2/3 spore coverage. The PSII was calculated according to the formula:
PSII=n1+n2×2+n3×3N×3×100%,
where n1, n2, and n3 represent the number of seeds of grade 1, grade 2, and grade 3, respectively, and N represents the total number. The SPSS 25 software [IBM Corporation] was used for the statistical analysis of PSII. The broad-sense heritability (H2) for PSII was calculated using the equation H2 = σg2/(σg2+σgXe2/n+σe2/rn), where σg2 is the genetic variance component, σgXe2 is the genotype-environment interaction variance component, σe2 is the residual (error) variance component, and n and r were defined as the number of environments and replications, respectively.
Library Construction and Sequencing
Young leaves were collected from the RIL population (F8) and used to extract genomic DNA using the cetyltrimethylammonium bromide (CTAB) method. The 300- to 500-bp Illumina Paired-end libraries were constructed for the 200 RILs and their parents. The PE150 reads for each library were generated by the Illumina HiSeq 2500 platform (Illumina, Inc., San Diego, CA, USA). The sequencing data have been deposited in the Sequence Read Archive database under accession number PRJNA760938. The raw sequencing reads were filtered by the Cutadapt (Marcel Martin Revision) and Trimmomatic software (Bolger et al., 2014) to generated high-quality clean reads. Then, the clean reads were aligned to the peanut reference genome (https://www.peanutbase.org/data/public/Arachis_hypogaea/Tifrunner.gnm1.KYV3/) using the BWA software (Li and Durbin, 2009). The HaplotypeCaller module of the GATK software (Broad Institute) was used to identify SNPs and InDels.
Linkage Map Construction
Raw SNPs and InDels were filtered with the following criteria: (a) minor allele frequency (MAF) ≥ 0.2; (b) relative heterozygosity rate ≤ 0.2 (calculated according to Wu et al., 2016); (c) the proportion of missing genotypes ≤ 0.5; (d) polymorphism exists between parents. Then, the high-quality SNPs and InDels were used to construct a bin map using the MPR package (Xie et al., 2010). Linkage groups were named according to the locations of SNPs/InDels on the reference genome.
QTL Analysis
The QTL mapping of the phenotypic data was performed using the default setting of the BIP (QTL mapping in bi-parental populations) approach in the IciMapping 4.2 software (Meng et al., 2015). The scanning step was set as 1 cM, and the logarithm of the odds (LOD) threshold was set as 2.5 to detect additive QTLs. The QTLs were named with the initial letter “q” followed by the trait name “PSII” and linkage group. An English lowercase letter was added if two or more QTLs were identified in the same linkage group. For instance, two QTLs for PSII were detected on LG A08, and then they were named qPSIIA08.a and qPSIIA08.b. If multiple QTLs overlapped on the same linkage group, they are considered to be a consistent QTL across environments and designated with the same name.
ResultsPhenotypic Variation of Resistance to <italic>Aspergillus flavus</italic> Infection
The PSII of the RIL population was calculated with seeds harvested from 4 consecutive years (2017–2020). The PSII of the female parent Zhonghua 16 was significantly higher than the male parent J11. The PSII of Zhonghua 16 ranged from 87.04 to 93%, whereas J11 ranged from 51.04 to 55.97% in the four environments. The PSII varied among RILs from 46.92 to 98.25%, 40.74 to 100%, 26.98 to 100%, and 24.73 to 100% in the 4 years (Table 1). Continuous distributions with transgressive segregation were observed, which indicated that both parents contain resistant genes against A. flavus infection (Figure 1). The correlation coefficient of the PSII across the 4 years ranged from 0.374 to 0.96, which were significant at the p < 0.01 level (Supplementary Table 1). The results of an ANOVA for PSII showed significant differences among genotypes, environments, and genotypes × environments interactions at P < 0.001 (Table 2). The broad-sense heritability of the PSII was estimated to be 0.76, indicating that the PSII was mainly controlled by genetic factors.
Phenotypic variation in PSII of the RIL population in Wuhan during 4 years.
Env
Parents
RILs
CV
Zhonghua 16
J11
Range
Mean ± SD
WH2017
87.04 ± 9.44
55.97 ± 8.57*
46.92–98.25
80.58 ± 10.57
0.13
WH2018
91.15 ± 6.26
53.83 ± 6.68**
40.74–100.00
83.74 ± 11.20
0.13
WH2019
93.00 ± 7.06
51.04 ± 6.70**
26.98–100.00
83.35 ± 11.55
0.14
WH2020
91.73 ± 5.80
51.65 ± 6.76**
24.73–100.00
67.87 ± 16.18
0.24
PSII, percent seed infection index; Env, environment; WH, Wuhan, China; SD, standard deviation; CV, coefficient of variation;
Difference is significant at p < 0.05 level between parents;
Difference is significant at p < 0.01 level between parents.
Phenotypic distribution of PSII in the RIL population and their parents harvested in Wuhan for 4 years. PSII, percent seed infection index.
Analysis of variance for the PSII of the RIL population across four environments.
Source
df
Mean square
F-value
P-value
Genotypes
199
0.08
8.71
<0.001
Environments
3
16.55
1776.57
<0.001
Genotypes × Environments
551
0.02
2.37
<0.001
Error
1,291
0.01
PSII, percent seed infection index.
Sequencing and Construction of Genetic Map
About 5 billion reads were generated from the 200 RILs and their parents. The female parent Zhonghua 16 was sequenced at 8.95 × coverage and the male parent J11 at 8.35 × coverage, while the RIL population individuals were sequenced at ~2.96 × coverage (Supplementary Table 2). On average, 84.22% of reads were uniquely aligned to the reference genome (Supplementary Table 3). A total of 233,365 SNPs/InDels were used to construct a high-density genetic map. These markers were divided into 2,802 recombination bins, which generated the final genetic map covering 1573.85 cM with an average inter-bin interval of 0.58 cM (Table 3; Figure 2). There were 1,257 bins for A sub-genome with a map length of 760.69 cM, and 1,545 bins for B sub-genome with a map length of 813.16 cM. The length of LGs varied from 57.28 (A06) to 96.64 cM (A09), and the number of bins in LGs ranged from 102 (A06) to 210 (B05).
Basic information of the high-density genetic linkage map of the RIL population.
Chr
Length (cM)
No. markers
No. bins
Marker interval (cM)
Bin interval (cM)
Max interval (cM)
A01
83.85
3,420
119
0.025
0.71
4.19
A02
61.50
13,895
124
0.004
0.50
3.99
A03
72.62
3,185
110
0.023
0.67
4.49
A04
86.66
2,042
120
0.042
0.73
12.34
A05
70.59
3,169
116
0.022
0.61
2.15
A06
57.28
2,409
102
0.024
0.57
1.87
A07
68.56
5,634
122
0.012
0.57
3.30
A08
93.76
5,298
181
0.018
0.52
2.63
A09
96.64
6,017
147
0.016
0.66
2.76
A10
69.23
15,859
116
0.004
0.60
4.80
B01
88.76
22,268
164
0.004
0.55
2.43
B02
59.01
14,568
110
0.004
0.54
3.88
B03
89.19
5,201
167
0.017
0.54
2.62
B04
93.57
3,352
116
0.028
0.81
9.10
B05
93.00
23,930
210
0.004
0.45
2.43
B06
59.62
21,306
143
0.003
0.42
2.43
B07
79.98
20,055
164
0.004
0.49
3.30
B08
80.99
17,578
148
0.005
0.55
3.00
B09
95.53
24,980
174
0.004
0.55
5.74
B10
73.51
19,199
149
0.004
0.50
2.91
Whole
1573.85
233,365
2,802
0.013
0.58
4.02
Distributions of bins on 20 linkage groups of the high-density genetic linkage map constructed in the RIL population.
Evaluation of Genetic Map
To evaluate the quality of the high-density genetic map, the sources of the bins in each RIL were analyzed, and the results showed that the bins from each parent form continuous fragments as expected (Supplementary Figure 1). The collinearity analysis, which compared the genetic position of all bins to their physical position on the reference genome, indicated a high collinearity between LGs and their corresponding chromosomes (Supplementary Figure 2).
Detection of Additive QTLs for Resistance to <italic>Aspergillus flavus</italic> Infection
A genome-wide QTL analysis was performed using the high-density genetic map and the phenotypic data of the PSII from the 200 RILs in 4 consecutive years. A total of six additive QTLs were identified with 5.03–10.87% PVE (Table 4; Figure 3). In the WH2017, WH2018, WH2019, and WH2020 trials, four, two, two, and three QTLs could be identified and explained 35.63, 14.97, 16.33, and 22.3% total phenotypic variance, respectively. Their LOD values ranged from 2.63 to 5.97. Two QTLs were detected on B03, and the other four QTLs were on A05, A08, B01, and B10, respectively. The QTL qPSIIB10 was consistently detected in 4 years, showing 6.91–10.58% PVE. Both QTL qPSIIB03.a and QTL qPSIIB03.b were repeatedly detected in 2 years with 9.16–9.23 and 5.03–5.75% PVE, respectivelyin. The major QTL qPSIIA08 was only detected in 1 yearin with 10.87% PVE. In addition, the minor QTLs qPSIIA05 andin qPSIIB01 were only detected in 1 year.
Additive QTLs for resistance to Aspergillus flavus infection in the RIL population across four environments.
Distributions of QTLs for PSII on the genetic map in 4 consecutive years. PSII, percent seed infection index.
Phenotypic Effect of the Combinations of Resistant and Susceptible Alleles
The positive additive effects of qPSIIA05, qPSIIA08, qPSIIB01, qPSIIB03.a, and qPSIIB03.b indicated that the alleles from J11 confer resistance to infection. On the contrary, the negative additive effect of qPSIIB10 showed that the alleles from Zhonghua 16 were responsible for improving resistance to infection. The recombination of qPSIIA05, qPSIIA08, qPSIIB01, qPSIIB03.a, qPSIIB03.b, and qPSIIB10 were screened in the RIL population. As shown in Table 5 and Supplementary Table 4, RILs with the genotype R1R2R3R4R5R6 possessed all resistant alleles of the six QTLs from both parents, while those with the genotype S1S2S3S4S5S6 possessed all susceptible alleles of the six QTLs. The RILs with the genotypes R1R2R3R4R5S6 and S1S2S3S4S5R6 possessed resistant alleles from J11 and Zhonghua 16, respectively. Notably, the RILs with the genotype R1R2R3R4R5R6 showed significantly lower PSII than RILs with other genotypes (R1R2R3R4R5S6, S1S2S3S4S5R6, and S1S2S3S4S5S6), indicating that the combinations of resistant alleles could enhance resistance.
Phenotypic effect of the combinations of resistant and susceptible alleles of the six QTLs in the RIL population.
Genotype
No. of lines
WH2017 (%)
WH2018 (%)
WH2019 (%)
WH2020 (%)
R1R2R3R4R5R6
5
55.53 ± 6.17a
64.40 ± 7.85a
63.91 ± 7.57a
42.61 ± 11.27a
R1R2R3R4R5S6
4
81.67 ± 9.66b
89.02 ± 5.50b
88.52 ± 4.09b
67.34 ± 7.98b
S1S2S3S4S5R6
5
90.38 ± 4.53b
91.23 ± 8.89b
94.23 ± 5.78b
84.98 ± 10.37b
S1S2S3S4S5S6
10
87.37 ± 6.63b
93.27 ± 5.37b
93.58 ± 5.27b
75.63 ± 6.84b
“R1” to “R5” indicated the resistant allele from J11 for the top five QTLs in Table 4, while “R6” indicated the resistant allele from Zhonghua 16 for the last QTL qPSIIB10 in Table 4. “S1” to “S5” indicated the susceptible allele from Zhonghua 16 for the top five QTLs in Table 4, while “S6” indicated the susceptible allele from J11 for the last QTL qPSIIB10 in Table 4. The letters a and b after the phenotypic value were statistically different at p < 0.05 based on an ANOVA and Games–Howell multiple comparison.
Discussion
With the change of peanut harvest methods and the aggravation of global warming, the risk of peanut contamination by aflatoxin is also increasing. The most valid measure to control this disease is cultivating varieties with resistance to A. flavus attacks. However, the breeding of novel resistant varieties is still a challenge in peanuts, and high-yield resistant varieties were arduously unavailable in peanut production because the genetic mechanism of resistance to A. flavus is still unclear. In the present study, we analyzed the infection resistance of the RIL population derived from the internationally recognized peanut variety J11 with resistance to A. flavus in order to explore the quantitative genetic loci of resistance.
A high-quality genetic map is the fundamental basis of the QTL mapping of agronomic traits. Varshney et al. (2009) constructed the first map of the cultivated peanut, but this map only has 135 SSR markers. More genetic maps have been constructed with the development of abundant SSR markers. Luo et al. (2018) constructed a genetic map involving 830 SSR markers. However, the number of SSR markers is much lower than that of SNP markers in the genome. Relying on reduced-representation genome sequencing technologies, SNP markers were unearthed for the construction of a genetic map in peanuts. For example, Hu et al. (2018) built a genetic map containing 2,334 SNP/SSR markers and Liu N., et al. (2020) built a genetic map comprising 2,595 SNPs. Recently, two reports used whole-genome resequencing technology to develop genetic maps, containing 2,156 recombination loci (Agarwal et al., 2018) and 3,634 recombination bins (Liu H., et al., 2020), respectively. In the present study, we re-sequenced 200 RILs and both parents. upon which 233,365 SNPs/InDels were identified and formed 2,802 recombination bins on the genetic map. Compared with the previous maps, this map was high-density and would provide convenience for the mapping of resistance to A. flavus infection and other traits.
The QTL for the resistance of peanut seeds to A. flavus infection could explain the relatively low amount of total variance. Yu et al. (2019) screened the PSII of 140 RILs in two environments, and the identified QTLs explained 13 and 19.28% of the total phenotypic variation, respectively. Khan et al. (2020) evaluated the PSII of 208 RILs in two environments, and the total phenotypic variations explained by the identified QTLs were 20.81 and 24.19%, respectively. In the present study, QTLs were identified using 200 RILs in four environments, and the total phenotypic variations explained by identified QTLs were 35.63, 14.97, 16.33, and 22.3%, respectively. Therefore, the resistance of peanut seeds to A. flavus infection might be controlled by multiple small effect loci and influenced by the environment; thus, a significantly larger mapping population and a more accurate resistance phenotyping method should be used in future studies to map these QTLs. Forty years ago, the peanut variety J11 was discovered to be resistant to A. flavus infection (Mehan and McDonald, 1980). Subsequently, its stable infection resistance was proven in multiple pieces of research (Mehan et al., 1981, 1987; Nigam et al., 2009). However, this precious material has not been fully utilized, for the resistant loci or genes of J11 were not found. In the present study, the QTL mapping of its resistance was first studied using an RIL population. Specifically, six additive QTLs for resistance to A. flavus infection were identified on chromosomes A05, A08, B01, B03, and B10. These QTLs were novel because the previously reported QTLs by Yu et al. (2019) and Khan et al. (2020) were on A03, A10, and B04. Notably, the favorable alleles of five of the six novel QTLs were from J11, but these QTLs were not stable across environments. Three of them were only identified in one environment, and the other two QTLs were identified in two of four environments, indicating that the resistance to A. flavus infection of J11 might be controlled by multiple QTLs and was significantly influenced by growth environments. This phenomenon was consistent with the previous observation that drought stress, pod maturity (Mehan et al., 1986), and seed coat integrity (Asis et al., 2005) were related to the resistance of J11 to A. flavus infection. However, the RILs possessing the five resistant alleles from J11 had lower PSII than those with susceptible alleles but were not significant, indicating that the identified QTLs might be part of the genetic factors controlling resistance. Therefore, J11 could be crossed with more varieties with different backgrounds to construct more and larger populations to further dissect its resistance to A. flavus infection.
Surprisingly, the favorable alleles of the stable QTL qPSIIB10 were from the susceptible parent Zhonghua 16, indicating that attention should also be paid to susceptible germplasms since they might possess valuable loci for resistance to A. flavus infection. Another important finding is that the RILs pyramided the favorable alleles of the six additive QTLs that significantly improved their resistance levels across all environments, although the effect of a single QTL was influenced by the environment. Therefore, both major QTLs and minor QTLs should be concerned with the genetic enhancement of resistance to A. flavus infection in peanut breeding.
Conclusion
The present study constructed a high-density genetic map with 2,802 recombination bins based on whole-genome resequencing technology. The results identified six important QTLs for peanut seed resistance to A. flavus infection across four environments. A further analysis of QTL revealed that the accumulation of all favorable alleles of the six QTLs could significantly enhance the infection resistance, thus laying a foundation for further research on fine mapping and breeding application in peanuts.
Data Availability Statement
The data presented in the study are deposited in the Sequence Read Archive repository, accession number PRJNA760938.
Author Contributions
YJ, HL, YL, BL, and HJ conceived, designed, supervised the experiments, and contributed to the final editing of the manuscript. YL, BL, and HJ developed the RIL population. YJ, HL, BY, YD, JG, LH, XZ, WC, and NL conducted field trials and phenotyping. YJ, BY, YD, JG, and WC performed DNA extraction and genotyping. YJ, HL, XZ, BL, and HJ performed the data analysis and interpreted the results. YJ and HL prepared the first draft. All authors read and approved the final manuscript.
Funding
This study was supported by the National Natural Science Foundations of China (Grant Numbers 31971903, 32001510, and 31761143005), the China Agriculture Research System (Grant Number CARS-13), the Crop Germplasm Resources Protection Project (Grant Number 2020NWB033), the Peanut Germplasm Resource Sharing Service Platform Project (Grant Number NCGRC-2020-036), and the Central Public-interest Scientific Institution Basal Research Fund (Grant Number 1610172019008).
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.
Publisher's Note
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.745408/full#supplementary-material
ReferencesAgarwalG.ClevengerJ.PandeyM. K.WangH.ShasidharY.ChuY.. (2018). High-density genetic map using whole-genome resequencing for fine mapping and candidate gene discovery for disease resistance in peanut. Plant Biotechnol J16:1954–1967. 10.1111/pbi.1293029637729AsisR.BarrionuevoD. L.GiordaL. M.NoresM. L.AldaoM. A. (2005). Aflatoxin production in six peanut (Arachis hypogaea L.) genotypes infected with Aspergillus flavus and Aspergillus parasiticus, isolated from peanut production areas of Cordoba, Argentina. J. Agric. Food Chem.53, 9274–9280. 10.1021/jf051259+16277433BolgerA. M.MarcLBjoernU (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinform. 15, 2114–2120. 10.1093/bioinformatics/btu17024695404DesjardinsA. (2003). Mycotoxins: Risks in Plant, Animal and Human Systems. AAHE-ERIC/Higher Education Research Report. Vol. 9. Ames, IW: Council for Agricultural Science and Technology, 48–50FAOSTAT (2019). Statistical database FAOSTAT. Available online at: http://faostat3.fao.org (accessed May 20, 2021).HuX. H.ZhangS. Z.MiaoH. R.CuiF. G.ShenY.YangW. Q.. (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. 10.1038/s41598-018-23873-729615772KewM. C. (2013). Aflatoxins as a cause of hepatocellular carcinoma. J. Gastrointestin. Liver Dis.22, 305–310. 10.1586/17474124.2013.83250024070156KhanS. A.ChenH.DengY.ChenY.ZhangC.CaiT.. (2020). High-density SNP map facilitates fine mapping of QTLs and candidate genes discovery for Aspergillus flavus resistance in peanut (Arachis hypogaea). Theor. Appl. Genet.133, 2239–2257. 10.1007/s00122-020-03594-032285164KisyombeC. T.BeuteM. K.PayneG. A. (1985). Field evaluation of peanut genotypes for resistance to infection by Aspergillus parasiticus. Peanut Sci.12, 12–17. 10.3146/pnut.12.1.0004LiHDurbinR. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinform. 25:17541760. 10.1093/bioinformatics/btp32419451168LiuH.SunZ.ZhangX.QinL.QiF.WangZ.. (2020). QTL mapping of web blotch resistance in peanut by high-throughput genome-wide sequencing. BMC Plant Biol.20:249. 10.1186/s12870-020-02455-832493219LiuN.GuoJ.ZhouX.WuB.HuangL.LuoH.. (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. 10.1007/s00122-019-03438-631559527LuoH.GuoJ.RenX.ChenW.HuangL.ZhouX.. (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. 10.1007/s00122-017-3000-729058050MahatoD. K.LeeK. E.KamleM.DeviS.DewanganK. N.KumarP.. (2019). Aflatoxins in food and feed: an overview on prevalence, detection and control strategies. Front. Microbiol.10:2266. 10.3389/fmicb.2019.0226631636616MehanV. K.McDonaldD. (1980). Screening for resistance to Aspergillus flavus invasion and aflatoxin production in groundnuts, in ICRISAT Groundnut Improvement Program Occasional, Paper-2 (Hyderbad), 1–15.MehanV. K.McdonaldD.NigamS. N.LalithaB. (1981). Groundnut cultivars with seed resistant to invasion by Aspergillus flavus. Oleagineux36, 501–507.MehanV. K.McdonaldD.RajagopalanK. (1987). Resistance of peanut genotypes to seed infection by Aspergillus flavus in field trials in India. Peanut Sci.14, 17–21. 10.3146/i0095-3679-14-1-5MehanV. K.McdonaldD.RamakrishnaN.WilliamsJ. K. (1986). Effects of genotype and date of harvest on infection of peanut seed by Aspergillus flavus and subsequent contamination with aflatoxin. Peanut Sci.13, 46–50. 10.3146/i0095-3679-13-2-1MengL.LiH.ZhangL.WangJ. (2015). QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J.3, 269–283. 10.1016/j.cj.2015.01.001NayakS. N.AgarwalG.PandeyM. K.SudiniH. K.JayaleA. S.PurohitS.. (2017). Aspergillus flavus infection triggered immune responses and host-pathogen cross-talks in groundnut during in-vitro seed colonization. Sci. Rep.7:9659. 10.1038/s41598-017-09260-828851929NigamS. N.WaliyarF.ArunaR.ReddyS. V.KumarP. L.CraufurdP. Q.. (2009). Breeding peanut for resistance to aflatoxin contamination at ICRISAT. Peanut Sci.36, 42–49. 10.3146/AT07-008.1PittetA. (1998). Natural occurrence of mycotoxins in foods and feeds - an updated review. Revue Méd. Vétérinair149, 479–492.TorresA. M.BarrosG. G.PalaciosS. A.ChulzeS. N.BattilaniP. (2014). Review on pre and post-harvest management of peanuts to minimize aflatoxin contamination. Food Res. Int.62, 11–19. 10.1016/j.foodres.2014.02.023VarshneyR. K.BertioliD. J.MoretzsohnM. C.VadezV.KrishnamurthyL.ArunaR.. (2009). The first SSR-based genetic linkage map for cultivated groundnut (Arachis hypogaea L.). Theor. Appl. Genet.118, 729–739. 10.1007/s00122-008-0933-x19048225WangH.LeiY.YanL.WanL.RenX.ChenS.. (2016). Functional genomic analysis of Aspergillus flavus interacting with resistant and susceptible peanut. Toxins8:46. 10.3390/toxins802004626891328WomackE. D.WilliamsW. P.WindhamG. L.XuW. (2020). Mapping quantitative trait loci associated with resistance to aflatoxin accumulation in maize inbred Mp719. Front. Microbiol.11:45. 10.3389/fmicb.2020.0004532117099WuZ.WangB.ChenX.WuJ.KingG. J.XiaoY.. (2016). Evaluation of linkage disequilibrium pattern and association study on seed oil content in Brassica napus using ddRAD sequencing. PLoS ONE11:e0146383. 10.1371/journal.pone.014638326730738XieW.FengQ.YuH.HuangX.ZhaoQ.XingY.. (2010). Parent-independent genotyping for constructing an ultrahigh-density linkage map based on population sequencing. Proc. Natl. Acad. Sci. U.S.A107, 10578–10583. 10.1073/pnas.100593110720498060YuB.HuaiD.HuangL.KangY.RenX.ChenY.. (2019). Identification of genomic regions and diagnostic markers for resistance to aflatoxin contamination in peanut (Arachis hypogaea L.). BMC Genet.20:32. 10.1186/s12863-019-0734-z30866805ZhaoX.LiC.YanC.WangJ.YuanC.ZhangH.. (2019). Transcriptome and proteome analyses of resistant preharvest peanut seed coat in response to Aspergillus flavus infection. Electron J. Biotechn.39, 82–90. 10.1016/j.ejbt.2019.03.003