- 1Centro de Estudios Avanzados en Fruticultura (CEAF), Rengo, Chile
- 2Department of Plant Breeding, Centro de Edafología y Biología Aplicada del Segura - Consejo Superior de Investigaciones Científicas (CEBAS-CSIC), Murcia, Spain
- 3Department of Pomology, Estación Experimental de Aula Dei - Consejo Superior de Investigaciones Científicas (EEAD-CSIC), Zaragoza, Spain
- 4Department of Science, Institute for Multidisciplinary Research in Applied Biology-IMAB, Universidad Pública de Navarra, Pamplona, Spain
- 5Department of Plant and Environmental Sciences, Clemson University, Clemson, SC, United States
Linkage maps are highly appreciated tools for cultivar and rootstock breeding programs because they are suitable for genetic and genomic studies. In this study, we report on using sequence-based genotyping (SBG) approach to simultaneously discover and genotype SNPs from two peach-based rootstocks (“Adafuel” and “Flordaguard”) and their progeny (n = 118): from a initial mapping population composed of 131 seedlings. The plant material was developed at the EEAD–CSIC Prunus rootstocks breeding program, aiming to obtain a segregating progeny for a range of characters of agronomical interest to rootstock breeding (iron-chlorosis and root-asphyxia tolerance, nematode resistance, vigor traits, and other effects on scion cultivars). Sequence reads obtained from double-digest SBG were aligned to the P. persica reference genome (Peach v2.0). While eight linkage groups were constructed for “Adafuel,” only four linkage groups were constructed for “Flordaguard,” given the low heterozygosity of this last genotype. High synteny and co-linearity were observed between obtained maps and Peach v2.0. On the other hand, this work aimed to elucidate the genetic basis of leaf chlorosis tolerance using the phenotypic segregation of the progeny to iron-chlorosis tolerance, along with the QTLs responsible for leaf chlorosis. The F1 mapping population, composed initially of 131 seedlings, was growing in four field trials established on calcareous soils at the experimental field of the EEAD-CSIC in Zaragoza, Spain. From the initial mapping population, 131 individuals were selected for their phenotypical characterization with SPAD measurements of plants grown in the field, exhibiting a great variability. Significant QTLs associated with tolerance to iron chlorosis were found in LG1, LG5, LG7, and LG8. The significant QTLs detected in LG5 and LG7 have not been associated with this abiotic stress before in Prunus. Several candidate genes such as Prupe.1G541100, predicted as glutamyl-tRNA reductase 1, Prupe.1G468200, encoding a 2-oxoglutarate (2OG), and Fe(II)-dependent oxygenase superfamily protein or Prupe.1G577000 (ppa011050.m), a NIFU-like protein 2 (NIFU2) were detected. The exact biological function of some of these genes should be verified for the future development of marker-assisted selection for peach iron chlorosis tolerance.
Introduction
Prunus is one of the largest genera of the family Rosaceae and subfamily Prunoideae (Rehder, 1940). It is economically important as a source of edible and highly appreciated fruits (e.g., peaches, plums, cherries, and almonds) but also as a source of rootstocks for stone fruits production. Modern fruit trees are currently grafted onto a rootstock (Guajardo et al., 2020). The rootstock not only allows scion cultivars to be adapted to unfavorable soil conditions but also determines tree size, vigor, and fruit quality (Font i Forcada et al., 2020). Biotic and abiotic stress tolerance related to soil characteristics is mainly determined by rootstocks (Cummins and Aldwinckle, 1983; Layne, 1987). While new stone fruit scion cultivar development is addressed by many breeding programs around the world, few of them are dedicated to the selection of new Prunus rootstocks (Byrne et al., 2012).
One of the main objectives of these breeding programs for rootstocks is to generate complex hybrids with the maximum number of desirable traits. The new rootstocks should overcome soil and disease problems (Reighard, 2002), showing high-performance and excellent adaption to different regions such as the Mediterranean area. One worldwide abiotic stress presented in these regions is iron chlorosis (Marschner, 1995). This abiotic stress is mainly due to the low iron bioavailability in calcareous soils and the difficulty of iron (Fe) acquisition by the roots (Hell and Stephan, 2003). The selection of rootstocks tolerant to iron chlorosis requires expensive and long-term experiments (Zarrouk et al., 2005; Jiménez et al., 2008). According to Gonzalo et al. (2012), the SPAD (Soil Plant Analysis Development) value seems to be the most reliable measure, and an unbiased and non-destructive method to evaluate the tolerance to iron chlorosis of different genotypes and breeding populations. The selection process can be also accelerated using molecular markers linked to the trait, allowing the early selection of likely seedlings (Etienne et al., 2002; Dirlewanger et al., 2006). Combining SPAD value and low number of molecular markers, Gonzalo et al. (2012) identified different quantitative trait loci (QTL) in linkage groups (LG) 4 and 6 in the almond–peach hybrid “Felinem” (Prunus dulcis × Prunus persica), formerly named “GxN 22” (Felipe, 2009). Although this work allowed localizing the candidate gene PFIT in LG 4, the genetic control of iron chlorosis in Prunus is still unknown. In fact, other important candidates such as PIRT1, PAHA2, and PNramp were not located by the authors (Gonzalo et al., 2012).
At present, new advances in genetic and genomic tools can help to face challenges in both cultivar and rootstocks breeding programs by improving the efficiency and reducing the time necessary to obtain new genotypes. In this sense, high-resolution linkage maps are valuable tools in genetics and breeding because they are the basis for QTLs and the identification of genes associated with traits of interest. For Prunus rootstocks, traditionally linkage maps have been constructed using amplified fragment length polymorphism (AFLPs) markers (Lu et al., 1998), AFLP and single sequence repeat (SSRs) markers (Blenda et al., 2007), and SSRs only (Gonzalo et al., 2012). In the last decade, the use of molecular markers such as single nucleotide polymorphisms (SNPs) has allowed the construction of high-density linkage maps for Prunus cultivars, but their use in linkage maps construction for rootstocks is scarce. For example, mapping populations have been genotyped using SNP chips (Peace et al., 2012; Verde et al., 2012) for the construction of linkage maps in peach (Eduardo et al., 2013; Martínez-García et al., 2013; Pirona et al., 2013; Yang et al., 2013; Frett et al., 2014; Romeu et al., 2014; da Silva et al., 2015, 2018; Donoso et al., 2015; Núñez-Lillo et al., 2015; Fresnedo-Ramírez et al., 2016; Lambert et al., 2016; Zeballos et al., 2016; Serra et al., 2017), sweet cherry (Klagges et al., 2013; Campoy et al., 2015; Calle et al., 2018), and sour cherry (Cai et al., 2018).
More recently, Prunus linkage maps have been constructed using SNPs identified through sequence-based genotyping (SBG, Elshire et al., 2011), a next-generation sequencing methodology that provides a reduced representation of the genome by using single and/or double-digestion of DNA with restriction enzymes to reduce genome complexity and genotype multiple DNA samples (Elshire et al., 2011; Poland et al., 2012). Single digestion using ApeKI restriction enzyme has been used for Prunus linkage maps construction in peach (Bielenberg et al., 2015; Núñez-Lillo et al., 2019), sweet cherry (Guajardo et al., 2015), Japanese plum (Salazar et al., 2017; Carrasco et al., 2018), almond (Goonetilleke et al., 2018), P. mume (Kitamura et al., 2018), and apricot (Pina et al., 2021). Although good genome coverage has been reported using this methodology, it has been described that the double-digest protocol greatly simplifies quantification of the library prior to sequencing and could generate a suitable and uniform complexity reduction of the genome (Poland et al., 2012; Guajardo et al., 2020).
In this study, we used the double-digest SBG for the identification of SNPs and to construct linkage maps for a mapping population derived from the cross “Adafuel” × “Flordaguard,” two Prunus rootstocks. Filtered sequence reads were aligned to the version 2.0 of the Prunus reference genome (Verde et al., 2017). The physical position of each SNP was determined, and the distribution of SNPs along eight peach pseudomolecules was performed to analyze the location of each SNP, identifying SNP-carrying genes and predicting the SNP effect of each SNP. Linkage maps constructed for both peach-based rootstocks were used for the identification of QTLs for iron chlorosis evaluated during two years in the mapping population. The results obtained here will provide a list of genes relevant for Prunus rootstocks breeding.
Materials and Methods
Plant Material
The mapping population was comprised of 118 (of the total 131) seedlings obtained from the cross between the peach–almond rootstock “Adafuel” (P. dulcis × P. persica; 2n = 2x = 16; Cambra, 1990) and the peach-based rootstock “Flordaguard” (P. persica; 2n = 2x = 16; Sherman et al., 1991). “Adafuel” is an almond × peach hybrid rootstock selected from an open-pollinated population of P. dulcis cv. Marcona for adaptation to calcareous and unfavorable soils for peach growing (Moreno et al., 1994). It is tolerant to iron chlorosis and susceptible to root asphyxia (Mestre et al., 2015). In turn, “Flordaguard” is a red leaf peach rootstock, a sixth-generation descendant from the cross Chico 11 × P. davidiana (Carr.) Franch, C-26712. Chico 11 was a seedling of the ornamental P. persica cv. Shau Thai, PI 65821 (Sharpe et al., 1969). After selection for nematode resistance, “Flordaguard” was finally released (Sherman et al., 1991). “Flordaguard” is very sensitive to iron-chlorosis and could be more tolerant to root asphyxia than “Adafuel.” The segregating population was established at the experimental trail (block B1) of the EEAD-CSIC in Zaragoza (North-Eastern of Spain; lat. 41° 43′ 42.7″ N, long. 0° 48′ 44.1″ W). To ensure enough replications for iron chlorosis analyses, at least 10 clonally propagated plants of each genotype (considered as an individual plot inside of each block) were planted in three additional experimental trials (B2, B3, and B4) of the EEAD-CSIC, established on calcareous soil having 29–31% total calcium carbonate, 8.1–9.9% active lime, pH 8.5, and a silt-loam texture.
Population Phenotyping
The iron chlorosis level was assessed with a SPAD-502 Chlorophyll Meter (Minolta Co., Osaka, Japan) in ungrafted clonally propagated plants of 131 genotypes grown in the field. The chlorophyll (Chl) concentration per unit leaf area was estimated using 10 leaves per tree, selected from the middle of bearing shoots located all around the crown, to obtain an average leaf value representative of the leaves belonging to the outer part of the tree canopy. SPAD measurements were carried out ~120 days after full bloom in the four experimental trials (blocks) for two years. For each genotype, at least three ungrafted trees were assessed in each experimental trial.
A mixed model was developed to represent our experimental design and to obtain the adjusted means (lsmeans) using the procedure PROC MIXED within SAS (SAS Institute Inc, 2017). The calculated lsmean for each individual tree was used in the QTL analysis. The mixed model was:
where Y_ikm is the mean value per plot of clone i in block k and year m, μ is the intercept, α_i and γ_m are the fixed effects of the m-th clone and i-th year, (αγ)_im the random effect of the i-th clone error at the m-th year, β_ik is the random effect of the k-th block within the i-th year, and ε_imk is the effect of the error of plot mean y_ikm. For all random effects including the error, year-specific variances and an additional covariance were fitted, if these extensions increase the model via AIC (Akaike, 1974).
DNA Extraction
Genomic DNA was extracted from young leaves collected from the parents and seedlings and stored at−80°C until use. The homogenization of leaf samples was reached after grinding around 100 mg of tissue with the Mixer Mill MM 400 (Retsch, Germany), and the tubes were immersed in liquid nitrogen to freeze the tissue. The NucleoSpin Plant II kit (Macherey-Nagel, Germany) was used for DNA extraction according to the manufacturer's instructions using the Lysis Buffer PL1, which is based on the CTAB method. DNA was eluted in 50 μL buffer PE and it was quantified using spectrophotometry (Tecan Trading AG, Switzerland). DNA quality was determined using 1% agarose electrophoresis.
Sequence-Based Genotyping (SBG)
Sequence-based genotyping was carried out at Clemson University Genomics Computational Laboratory (CUGCL; Clemson, SC, USA). A reduced representation SBG library was prepared using restriction enzymes, PstI (methylation sensitive) and MspI (partial sensitivity to methylation), as described by Poland et al. (2012). A total of 200 ng of intact genomic DNA was digested and ligated to custom-designed adapter sequences. A total of 120 SBG libraries were sequenced on an Illumina HiSeq2500 using a 2 ×125 bp paired-end read module across 2 high-output lanes. Raw sequence data were demultiplexed and preprocessed for errors using the Stacks demultiplex tool (Catchen et al., 2013). Sample specific sequences were aligned to the eight pseudomolecules representing the eight chromosomes of the peach genome assembly (Peach v2.0; Verde et al., 2017) with the GMAP/GSNAP release 816.16 (Wu et al., 2016). The resulting variant call file (.vcf) was filtered for SNPs with a minimum depth (DP) of six, and present in at least 80% of the samples. Mean coverage of each SBG SNP was determined by creating a.BED file from the final SNP set and generating a bed graph with the genomecov function of bedtools v. 2.28.0 (Quinlan and Hall, 2010), and intersecting the bedgraph with the SNP.bed file with the intersect function in bedtools. The mean coverage of each sample was determined with in-house scripts. SNPs were extracted using the pipeline implemented in TASSEL 5.2.5 software (Bradbury et al., 2007), and accessions were called using minor allele frequency (MAF) > 0.05. SNPs were labeled according to the pseudomolecules in the peach genome (Pp01 to Pp08), followed by the physical position in base pairs (bp).
Linkage Maps Construction
Markers scoring followed the coding scheme for cross-pollinated population type proposed in JoinMap 4 (Kyazma B.V., Netherlands; Van Ooijen, 2006) and the evaluation version of JoinMap 5 (https://www.kyazma.nl/index.php/JoinMap/Evaluate). Linkage maps were constructed following the protocol described for sweet cherry by Guajardo et al. (2015). Markers with identical genotypic scores (Similarity of loci = 1), which are automatically removed by JoinMap because they should be mapped to the same position on the linkage group, were added back to the resulting linkage maps. Markers showing segregation distortion were also integrated into the maps and they were plotted using MapChart 2.3 (Voorrips, 2002). For a comparison between physical and genetic distance of mapped SNPs, the genetic positions of each SNP (in centimorgans, cM) were plotted against their physical position on the Peach v2.0 (in Megabase pairs, Mbp). Results were compared with those obtained by Verde et al. (2017) for the peach selection IF7310828 × Ferganensis BC1 (PxF) population using the IPSC 9 K SNP array (Verde et al., 2012) with the physical position of each SNP updated with the Peach v2.0 as a reference genome (Verde et al., 2017).
QTL Analysis of Iron Chlorosis Tolerance
Quantitative trait loci analysis of ion chlorosis was performed with MapQTL® 6.0 (Van Ooijen, 2009) using SPAD data from each individual for the two seasons evaluated. According to the no assumption of normality distribution of the data, the rank-sum test of Kruskal–Wallis (K-W; test non-parametric equivalent of the one-way analysis of the variance) was performed for the trait against each SNP, separately, to establish the strength of linkage at individual SNP loci. The K-W test ranks all individuals according to the quantitative trait, while it classifies them according to their marker genotype. A stringent significance level (P-value) similar or higher than 0.005 was used to identify markers significantly associated with both traits by the K-W test (Van Ooijen, 2009).
Structural and Functional Characterization of SNPs
The physical position of each SNP was used to identify common markers in this study and the IRSC 9K peach SNP array v1 (Verde et al., 2012), which was updated with the Peach v2.0 as a reference genome (www.rosaceae.org; Verde et al., 2017). The prediction of SNP effect, based on annotated genes or their genomic locations, was performed using the SnpEff v4.3 program (Cingolani et al., 2012) with Peach v2.0 as reference. Whenever multiple transcripts for a gene exist, the effect on each transcript was analyzed, which was the same criteria applied by Guajardo et al. (2020) for SNPs identified from a vast group of Prunus rootstocks, where “Adafuel” and “Flordaguard” were also included. The SNP predicted effects were categorized by impact, as high (disruptive impact on the protein), moderate (non-synonymous substitution), low (synonymous substitution), or modifier (with impact on noncoding regions). To examine the putative function of the genes containing high impact SNPs, a eukaryotic orthologous group (KOG) analysis was carried out using tools from Join Genome Institute (JGI, https://jgi.doe.gov). Information about GO, PFAM and Phanter for each gene of this group was obtained from Phytozome (https://jgi.doe.gov).
Results
Phenotyping
Different measurements were implemented to guarantee an accurate data set for the genetic analysis of 131 genotypes, clonally propagated and growing in three experimental fields. Leaf chlorophyll content, used as a good indicator of leaf chlorosis in trees, was estimated by using a SPAD meter. According to the Lsmeans obtained from our model, the values showed a great segregation response and significant differences between genotypes (p-value <0.001). Lsmeans ranged significantly among genotypes and years from 22.013 to 37.40 (Figure 1). No significant difference was observed between years (p-value 0.2753).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/872208/fpls-13-872208-HTML/image_m/fpls-13-872208-g001.jpg)
Figure 1. Frequency distribution of Lsmeans for SPAD values of 131 peach–almond individual hybrids from the cross of “Adafuel” × “Flordaguard” using the hist() function in R programming language (R Core Team, 2019).
High-Throughput Genotyping of the Mapping Population
Double-digest SBG produced between 4,293,116 and 18,572,896 read pairs, with an average of 11,321,053 reads per individual. This deep sequencing led to a mean depth per SNP of 1,244 across the entire dataset (Supplementary Table 1). The number of unique tags varied between 1,548,906 and 6,620,490, with an average of 3,990,190 tags per individual.
A total of 18,912 high quality SNPs (MAF > 0.05; missing data <5%), evenly distributed along the eight peach pseudomolecules, were used for mapping (Table 1). The number of identified SNPs ranged from 1,595 for Pp07 to 4,042 for Pp01. A total of 224.3 Mbp (99.4%) of the peach genome was covered with a marker density of ~1 SNP per 11.9 Kbp. Gaps were observed in all pseudomolecules, with the largest gap per pseudomolecule ranging from 506 Kbp (Pp06) to 1.2 Mbp (Pp01 and Pp04).
Linkage Mapping
Out of the total 18,912 SNPs, 17,539 (92.7%) exhibited the maternal segregation type < lmxll>, 868 (4.6%) showed the paternal segregation type < nnxnp> and 505 (2.7%) were heterozygous for both parents, with segregation type < hkxhk>. The last group of SNPs was not used for linkage maps construction. The linkage map obtained for the “Adafuel” genotype was composed of eight linkage groups (Ad-LG 1–8; Figure 2 and Table 2) and was comprised of 626 uniquely mapped SNPs covering a total of 450.6 cM. Ad-LG1 was the largest linkage group with 70.1 cM and Ad-LG4 covered the shortest distance, 39.7 cM. The average marker density was 0.8 cM per marker. The maximum gap size ranged from 2.4 cM in Ad-LG4 to 14.9 cM in Ad-LG7. In contrast, the “Flordaguard” linkage map was composed just only with four linkage groups (Fg-LG3, Fg-LG4, Fg-LG5, and Fg-LG7, Figure 2). This map comprised of 79 uniquely mapped SNPs, covering a total of 142.7 cM (Table 2), with a minimum and maximum linkage group length of 18.4 cM for Fg-LG7 and 58.3 cM for Fg-LG5 and an average marker density of 1.8 cM per marker. The maximum gap size ranged from 3.3 cM for Fg-LG7 to 13.1 cM for Fg-LG5. When markers with identical genotypic scores were added back to the resulting linkage maps, a significant clustering of SNPs co-segregating at the same map positions in some linkage groups were observed (Table 2 and Supplementary Table 2). For example, while the number of uniquely mapped SNPs on Ad-LG1 was 145, the number of co-segregating SNPs was 3,865 on Ad-LG1, resulting in a total of 7,802 SNPs mapped on “Adafuel” linkage map and 225 SNPs mapped on “Flordaguard” linkage map. Furthermore, and considering all the linkage groups of “Adafuel” and “Flordaguard,” a high number of skewed markers (p < 0.01) was mapped in the bottom part of Ad-LG5 (Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/872208/fpls-13-872208-HTML/image_m/fpls-13-872208-g002.jpg)
Figure 2. Linkage maps for “Adafuel” and “Flordaguard.” Asterisks indicate distortion level of skewed markers. *p < 0.1; **p < 0.05; ***p < 0.01; ****p < 0.005; *****p < 0.001; ******p < 0.0005; *******p < 0.0001. Redundant markers have been removed and only unique positions are represented. The scale in centiMorgan (cM) is given at the left side of the figure.
The important gaps observed for some linkage groups were related to considerable physical distance between flanking SNPs (Table 2 and Figure 3). For “Adafuel,” overlapping gaps between genetic and physical maps were observed in Ad-LG2 (8.8 cM with 9.7 Mbp), Ad-LG7 (14.9 cM with 9.1 Mbp), and Ad-LG8 (2.7 cM with 5.8 Mbp), and for “Flordaguard,” in Fg-LG4 (10.1 cM with 1.9 Mbp), Fg-LG5 (13.1 cM with 2.1 Mbp), and Fg-LG7 (3.3 cM with 2.5 Mbp). The distribution of SNPs is similar to those observed for the P × F map (Figure 3), although the genetic distance covered by each “Adafuel” and “Flordaguard” linkage group was shorter when it was compared with P × F linkage groups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/872208/fpls-13-872208-HTML/image_m/fpls-13-872208-g003.jpg)
Figure 3. Relationship between genetic and physical distances for each linkage group and pseudomolecule. Physical distance along each pseudomolecule is shown on the horizontal axis (Mbp), and the genetic distance is shown on the vertical axis (cM). The location of markers from “Adafuel” (green), “Flordaguard” (yellow), and P × F (blue) linkage maps are shown. Redundant markers have been removed and only unique positions are represented. The vertical bars indicate the putative position of the centromeric regions using information from the Peach v2.0 (Verde et al., 2017).
The linkage position of all SNPs mapped in “Adafuel” and “Flordaguard” linkage maps were in agreement with their physical position on the pseudomolecules of the Peach v2.0 (Figure 4). As it was previously indicated, it was possible to construct only four linkage groups for “Flordaguard,” and also partial coverage of pseudomolecules was observed for Fg-LG3, Fg-LG4, and Fg-LG7, with mapped SNPs only from the bottom region for Fg-LG3 and from the top region for Fg-LG4 and Fg-LG7.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/872208/fpls-13-872208-HTML/image_m/fpls-13-872208-g004.jpg)
Figure 4. Comparison of the parental linkage maps and physical map. The physical map was constructed using all SNPs identified by GBS for the mapping population used in this study (n = 19,912) and ordered according to their physical position in the Peach v2.0 (Verde et al., 2017). Redundant markers have been removed and only unique positions are represented. Distances have been proportionally modified to show the results.
QTL Detection for Ion Chlorosis Tolerance
To identify QTLs involved in iron chlorosis, analysis for each year and genotype was performed. The QTLs were placed on the genetic map of “Adafuel,” comprised of 626 uniquely mapped positions (SNPs), covering a total of 450.6 cM. The results obtained from the non-parametric Kruskal–Wallis for iron chlorosis revealed several QTLs in Ad-LG1, Ad-LG5, Ad-LG7, and Ad-LG8. The number of SNPs identified and the length of each QTL in each linkage group was different. In total, 1,263 SNP markers, located in 95 unique map positions, were strongly associated (from p = 0.05 to p = 0.0001) with this trait (Figure 5). In Ad-LG1, 1,192 SNPs were associated with SPAD, located mainly at the bottom of this linkage group and spanning a map distance of around 17 cM (between 51 and 68 cM). In Ad-LG1, the higher number of SNPs were mapped in three map positions at 55.27 cM, 58.11 cM, and 61.71 cM, with 205, 108, and 196 SNPs, respectively. In Ad-LG5, a total of 33 SNPs were significant, with a high LOD peak between 22 and 26 cM. In Ad-LG7, 29 SNPs were highly significant for SPAD, with a high LOD between 45 cM and 47 cM. At the bottom of Ad-LG8, a total of nine SNPs markers were highly associated with SPAD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/872208/fpls-13-872208-HTML/image_m/fpls-13-872208-g005.jpg)
Figure 5. QTLs obtained for each year of the study, and the lsmeans obtained with the fitted model for iron chlorosis in “Adafuel” (Ad). At the right side of each linkage group, dotted green line shows QTL LOD value variation along the map and the likelihood peak of the QTL. Significant markers by Kruskal-Wallis test statistic in Ad-LG1, Ad-LG5, Ad-LG7, and Ad-LG8 are highlighted in red. In Ad-LG1, Ad-LG5, and Ad-LG7 only high levels of significance are shown. Asterisks indicate distortion level of skewed markers. *p < 0.1; **p < 0.05; ***p < 0.01; ****p < 0.005; *****p < 0.001; ******p < 0.0005; *******p < 0.0001.
SNP Effect Prediction
When mapped SNPs for the “Adafuel” × “Flordaguard” seedlings were analyzed using the physical position of each SNP along the eight pseudomolecules of the peach genome, a total of 23,626 functional effects for SNP variants were predicted for 8,027 SNPs (7,802 SNPs for “Adafuel” and 225 SNPs for “Flordaguard;” Supplementary Table 3). Also, a group of 5,586 different genes (or 20.8% of genes identified in the peach genome sequence) could be affected by this group of SNPs. The predicted effects were of modifier type (78.52%), low impact (13.72%), moderate impact (7.73%), and high impact (0.03%). For SNPs with a modifier effect, downstream gene variant was the most represented category (38.14%), while synonymous variant (89.13%) and missense variant (99.15%) was the most represented category for low and moderate effect, respectively. Fifteen genes were categorized as highly impacted (Supplementary Table 4).
For SNP effect analysis of those SNPs strongly associated with iron chlorosis according to the QTL identified in Ad-LG1, Ad-LG5, Ad-LG7, and Ad-LG8, both uniquely and co-segregating SNPs were considered. The predicted effects were of modifier type (79.20%), low impact (12.12%), moderate impact (8.62%), and high impact (0.05%). In Ad-LG1, a total of 3,522 effects were detected. The number of effects was lower in the rest of pseudomolecules, with 103 effects in Ad-LG5, 93 in Ad-LG7, and 28 effects in Ad-LG8. These effects affected a total of 970 unique peach (Prupe) genes, as commented before, the majority of them located in Pp01. Several annotations for homologs of these peach genes have been provided, and in total, this list of genes matches with 754 homologs Arabidopsis genes (Supplementary Table 5). A total of 621 unique descriptions were obtained for these genes. Only two genes on Ad-LG1 (Prupe.1G566900 and Prupe.1G570100) were categorized as highly impacted because they present a stop codon. In addition, complete information on candidate genes located in the QTL regions is provided (Supplementary Table 6).
Discussion
In this work, a successful strategy covering the discovery of single nucleotide polymorphism (SNP) through double-digest sequence-based genotyping (SBG), genetic linkage maps construction, and QTL identification has been performed in a segregating peach–almond progeny to improve the knowledge about the genetic control of the tolerance to iron chlorosis in Prunus. The physical position and distribution of each SNP along the eight peach pseudomolecules were determined to identify SNP-carrying genes and to predict the effect of each SNP. To our knowledge, this is the first report on linkage map construction and QTL analysis using SNPs obtained from SBG for Prunus rootstocks.
Our results showed excellent performance for double digest-SBG and deep coverage paired-end sequencing for reduced representation of Prunus rootstocks. A high mean depth (1,244) per SNP and high coverage of peach pseudomolecules were reached, which was higher than a coverage of 12.1x obtained by Núñez-Lillo et al. (2019) for map construction in peach using single-digestion SBG. Our results confirm that paired-end sequence reads are more accurately mapped onto the reference genome, showing a significantly lower amount of missing data and a greater number of quality SNPs, as was suggested by Shirasawa et al. (2016) and Guajardo et al. (2020). In total, the SNP density was found to be ~1 SNP per 11.9 Kbp, which can be enough for mapping and QTL identification in an F1 population.
The length of “Adafuel” linkage map, composed by eight linkage groups, is comparable with other maps reported for peach cultivars (Martínez-García et al., 2013; Yang et al., 2013; Núñez-Lillo et al., 2015) and for almond x peach maps (Donoso et al., 2015). A high number of SNPs shared map positions in Ad-LG1, Ad-LG2, and Ad-LG3 (Table 2), probably due to the absence of recombination between some SNPs. On the other hand, and given the low number of SNPs showing paternal segregation (868 SNPs, or 4.6% of the useful SNPs for mapping), it was possible to construct only four linkage groups for “Flordaguard.” It might be due to the origin of this rootstock, which originated from six generations of peach crosses and open-pollinated seedlings (Sherman et al., 1991), generating low levels of heterozygosity in certain regions, or gaps, of the “Flordaguard” genome.
In general, for both “Adafuel” and “Flordaguard” linkage maps, the genetic positions of mapped SNPs are in agreement with the physical position of SNP along the pseudomolecules in the peach genome sequence (Peach v2.0; Verde et al., 2017), and the distribution of SNPs is similar to those observed for the PxF map (Figure 3), which reflect the correction of inter-pseudomolecule misassembles of the Peach v1.0 in the new version of the peach genome sequence.
A high number of skewed markers (p < 0.01) was mapped at the bottom part of Ad-LG5 (Figure 2). Although it has been described that the use of distorted markers could increase the genome coverage of the genetic maps (Zhang et al., 2007, 2013; Cai et al., 2015) and it may be beneficial for QTL mapping (Xu, 2008; Zhang et al., 2010), many authors eliminate the skewed markers and they are not used for linkage maps construction, diminishing the information that it is possible to obtain. It has been proposed that the study of segregation distortion is significant because distorted markers may be linked to genes or traits of interest, and these genes may be beneficial or lethal to organisms (Kirungu et al., 2019). Then, the exclusion of such markers could bias the data and result in the loss of some important genetic information. One example is the reported clustering of loci with distortion near the bottom of LG6 in sweet cherry in semi-compatible crosses, near the Prunus self-incompatibility S locus (Dirlewanger et al., 2004; Klagges et al., 2013; Guajardo et al., 2015). Until our knowledge, this is the first report of clustering of loci with high distortion in the bottom part of the LG5 in Prunus, and additional studies are necessary to determine the biological significance of the skewed markers in this region.
In a preliminary study in Prunus, two different mechanisms associated with tolerance to iron chlorosis were proposed (Gonzalo et al., 2012). The timing to overcome the lack of this essential element is the main difference between both mechanisms. The two linkage groups, LG1 and LG8, were associated with the delayed response observed in the Myrobalan plum P 2175 by Gonzalo et al. (2012). These two linkage groups were also identified in our study. However, QTLs associated with a quick response observed in LG4 and LG6 detected in “Felinem” were not identified in our study. On the other hand, the additional QTLs in LG5 and LG7 identified in “Adafuel” were not previously identified in response to iron chlorosis in this studied genotype.
A total number of 312 distinct genes were affected by significant SNPs associated with iron (Fe) deficiency-induced chlorosis, which has a large impact on yield in tree fruit crops, such as peach (Álvarez-Fernández et al., 2003; Abadía et al., 2011). However, the candidate gene PFIT, related to iron metabolism, and previously identified in Pp04 of Prunus (Gonzalo et al., 2012), was not identified in the candidate genes list obtained here. Peach and almond are considered to be Strategy I plants. The reduction-based strategy (Strategy I) has two main steps, the first step is the acidification of the rhizosphere that leads to an increase in the chelated Fe (III) concentration (Santi and Schmidt, 2009; Ivanov et al., 2012).
Subsequently, the root surface-localized ferric chelate reductase FERRIC-REDUCTION OXIDASE2 (FRO2; Robinson et al., 1999) reduces Fe(III) to soluble Fe(II), which is then taken up into epidermal cells by the Fe-REGULATED TRANSPORTER1 (IRT1; Eide et al., 1996). In our study, although FRO2 was not identified, two Ferric-chelate reductase genes [Prupe.1G587800 (FRO7) and Prupe.5G166400 (FRO08)] were affected by significant SNPs. Interestedly, both genes are located in different peach pseudomolecules. According to the literature, it is possible that once Fe3+-citrate complexes pass to the chloroplastic intermembrane space, FRO7 may reduce ferric (Fe3+) to ferrous iron (Fe2+) via its reductase activity (Jeong et al., 2008; Solti et al., 2016). A possible role of FRO7 in the reduction of Fe (III) to soluble Fe (II) in root Fe uptake responses to Fe-deficiency in peach should be investigated to clearly confirm this hypothesis.
Several studies using commercially available Affymetrix ATH1 GeneChips (Dinneny et al., 2008; Yang et al., 2010; Mai et al., 2016) and RNA-seq data set (Lan et al., 2013) have been carried out to study Fe deficiency and the involvement of ethylene in Strategy I plants. As a result, a total of 71 differentially expressed genes (61 being up-regulated and 10 down-regulated) overlapped between these studies. In the list of up-regulated genes was AT3G12900, which encodes a protein belonging to the 2OG-Fe(II) oxygenase family and is strongly induced by Fe deficiency. Proteins in this family are generally considered to possess oxidoreductase activity catalyzing the 2-oxoglutarate- and Fe(II)-dependent oxidation of an organic substrate using a dioxygen molecule. A previous study showed that the ethylene synthesis protein ACC oxidase also belongs to 2OG-Fe(II) oxygenase superfamily (Aravind and Koonin, 2001). In our study, the gene Prupe.1G468200 is located in Pp01 of peach and it was affected by eight significant SNPs mapped in this pseudomolecule and showing different effects, low, moderate, and modifier. This gene is an ortholog of AT1G52820 and encodes a 2-oxoglutarate (2OG) and Fe(II)-dependent oxygenase superfamily protein. The possible expression of this gene, by Fe deficiency, should be strongly verified in Prunus (peach or almond).
Another candidate gene was Prupe.1G541100 (ppa003848m.g), predicted as glutamyl-tRNA reductase 1, which had a modifier (3'-UTR variant) effect by a significant SNP located in Pp01 (44,202,463 bp) and mapped at 63.65 cM. The study of the transcriptional response of Arabidopsis to Fe deficiency has indicated that glutamyl-tRNA reductase (HEMA1), NYC1, and CGLD27 may represent key players in preventing photooxidative damage in Fe-deficient leaf cells (Rodríguez-Celma et al., 2013). The fast and pronounced down-regulation of HEMA1 expression suggests a rapid and dramatic decline of tetrapyrrole synthesis at an early stage of Fe deficiency. The possible implication of Prupe.1G541100 preventing photooxidative damage in Fe-deficient leaf cells in peach must be clearly confirmed.
Other important genes detected in our study were Prupe.5G039500 and Prupe.1G577000 (ppa011050.m), encoding a ferredoxin protein (FdC2) and a NIFU-like protein 2 (NIFU2), respectively. The interaction of these both genes in peach and their association with iron chlorosis response should be corroborated in further analysis. In this sense, in Arabdopsis, the study of the physiological and molecular changes during Fe deficiency confirmed that at the protein level, ferredoxin, the cytochrome-b6f complex, and Fe-containing enzymes of the plastid sulfur assimilation pathway were major targets of Fe deficiency, whereas other Fe-dependent functions were relatively less affected (Hantzis et al., 2018).
Conclusion
In the present study, we have constructed, to the best of our knowledge, the first Prunus rootstocks linkage maps using SNPs obtained from double-digest SBG. These linkage maps were the basis to identify QTLs and putative genes related to iron chlorosis in a segregating peach–almond progeny. SNPs identified in this study present a valuable set of new SNPs identified in Prunus that would be useful for genetic studies in the future. Candidate genes for iron chlorosis QTL were proposed; however, future fine mapping and transcriptomic analysis are needed to enable the identification of the underlying gene(s).
Data Availability Statement
The dataset generated for this study can be found in the NCBI-SRA database, BioProject number PRJNA773812.
Author Contributions
VG, PM-G, and MM conceived and designed the experiments. VG, PM-G, and SS performed the experiments. VG, SS, CS, PM-G, and AC-S analyzed the data. CS and MM contributed with reagents, materials, and analysis tools. All authors contributed to the manuscript preparation, read, and approved the final version.
Funding
This research was funded by the Spanish Ministry of Science, Innovation and Universities (MICINN) grants RFP2015-00019 and RTI2018-094176-R-C33, co-funded by ERDF A way of making EUROPE and the Regional Government of Aragon (A44 and T07-17R). It was also funded by CONICYT-REGIONAL/GORE O'HIGGINS/CEAF/R19A10003, FONDECYT 3160316, and CONICYT R16F20006 from Chile.
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.
Acknowledgments
The authors are grateful to Dr. Ana María Casas and Dr. María Pilar Vallés at EEAD-CSIC, Spain, for their helpful discussions and for providing equipment access.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.872208/full#supplementary-material
References
Abadía, J., Vázquez, S., Rellán-Álvarez, R., El-Jendoubi, H., Abadía, A., Álvarez-Fernández, A., et al. (2011). Towards a knowledge-based correction of iron chlorosis. Plant Physiol. Biochem. 49, 471–482. doi: 10.1016/j.plaphy.2011.01.026
Akaike, H.. (1974). A new look at the statistical model identification. IEEE Trans. Autom. Control 19, 716–723. doi: 10.1109/TAC.1974.1100705
Álvarez-Fernández, A., Paniagua, P., Abadía, J., and Abadía, A. (2003). Effects of Fe deficiency chlorosis on yield and fruit quality in peach (Prunus persica L. Batsch). J. Agric. Food Chem. 51, 5738–5744. doi: 10.1021/jf034402c
Aravind, L., and Koonin, E. V. (2001). The DNA-repair protein AlkB, EGL-9, and leprecan define new families of 2-oxoglutarate- and iron-dependent dioxygenases. Genome Biol. 2:7. doi: 10.1186/gb-2001-2-3-research0007
Bielenberg, D. G., Rauh, B., Fan, S., Gasic, K., Abbott, A. G., Reighard, G. L., et al. (2015). Genotyping by sequencing for SNP-based linkage map construction and QTL analysis of chilling requirement and bloom date in peach [Prunus persica (L.) Batsch]. PLoS ONE 10:e0139406. doi: 10.1371/journal.pone.0139406
Blenda, A. V., Verde, I., Georgi, L., Reighard, G., Forrest, S., Muñoz-Torres, M., et al. (2007). Construction of a genetic linkage map and identification of molecular markers in peach rootstocks for response to peach tree short life syndrome. Tree Genet. Genom. 3, 341–350. doi: 10.1007/s11295-006-0074-9
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Byrne, D. H., Raseira, M. C., Bassi, D., Piagnani, M. C., Gasic, K., Reighard, G. L., et al. (2012). “Peach,” in Fruit Breeding, Handbook of Plant Breeding 8, eds M. L. Badenes, and D. H. Byrne (Düsseldorf: Springer Science + Business Media, LLC 2012), 505–565. doi: 10.1007/978-1-4419-0763-9_14
Cai, C., Cheng, F. Y., Wu, J., Zhong, Y., and Liu, G. (2015). The first high-density genetic map construction in tree peony (Paeonia Sect. Moutan) using genotyping by specific-locus amplified fragment sequencing. PLoS ONE 10:e0128584. doi: 10.1371/journal.pone.0128584
Cai, L., Stegmeir, T., Sebolt, A., Zheng, C., Bink, M., and Iezzoni, A. (2018). Identification of bloom date QTLs and haplotype analysis in tetraploid sour cherry (Prunus cerasus). Tree Genet. Genomes 14:22. doi: 10.1007/s11295-018-1236-2
Calle, A., Cai, L., Iezzoni, A., and Wünsch, A. (2018). High-density linkage maps constructed in sweet cherry (Prunus avium L.) using cross and self-pollination populations reveal chromosomal homozygosity in inbred families and non-syntenic regions with the peach genome. Tree Genet. Genomes 14:37. doi: 10.1007/s11295-018-1252-2
Cambra, R.. (1990). ‘Adafuel', an almond × peach hybrid rootstock. HortScience 25, 584–584. doi: 10.21273/HORTSCI.25.5.584
Campoy, J. A., Le Dantec, L., Barreneche, T., Dirlewanger, E., and Quero-García, J. (2015). New insights into fruit firmness and weight control in sweet cherry. (2015). Plant. Mol. Biol. Rep. 33, 783–796. doi: 10.1007/s11105-014-0773-6
Carrasco, B., González, M., Gebauer, M., García-González, R., Maldonado, J., and Silva, H. (2018). Construction of a highly saturated linkage map in Japanese plum (Prunus salicina L.) using GBS for SNP marker calling. PLoS ONE 13:208032. doi: 10.1371/journal.pone.0208032
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
Cingolani, P., Platts, A., Wang, L. L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 6, 80–92. doi: 10.4161/fly.19695
Cummins, J. N., and Aldwinckle, H. S. (1983). “Rootstocks breeding,” in Methods in Fruit Breeding, eds J. N. Moore and J. Janick (West Lafayette: Purdue University Press), 294–327. doi: 10.1007/978-1-4684-8896-8_10
da Silva, C., Antanaviciute, L., Abdelghafar, A., Arús, P., Bassi, D., Rossini, L., et al. (2018). High-density multi-population consensus genetic linkage map for peach. PLoS ONE 13:207724. doi: 10.1371/journal.pone.0207724
da Silva, C., Bassi, D., Bianco, L., Pacheco, I., Pirona, R., and Rossini, L. (2015). Genetic dissection of fruit weight and size in an F2 peach (Prunus persica (L.) Batsch) progeny. Mol. Breed. 35:71. doi: 10.1007/s11032-015-0271-z
Dinneny, J. R., Long, T. A., Wang, J. Y., Jung, J. W., Mace, D., Pointer, S., et al. (2008). Cell identity mediates the response of Arabidopsis roots to abiotic stress. Science 320, 942–945. doi: 10.1126/science.1153795
Dirlewanger, E., Cosson, P., Boudehri, K., Renaud, C., Capdeville, G., Tauzin, Y., et al. (2006). Development of a second-generation genetic linkage map for peach [Prunus persica (L.) Batsch] and characterization of morphological traits affecting flower and fruit. Tree Genet. Genomes 3, 1–13. doi: 10.1007/s11295-006-0053-1
Dirlewanger, E., Graziano, E., Joobeur, T., Garriga-Calderé, F., Cosson, P., Howad, W., et al. (2004). Comparative mapping and marker-assisted selection in Rosaceae fruit crops. Proc. Natl. Acad. Sci. U. S. A. 101, 9891–9896. doi: 10.1073/pnas.0307937101
Donoso, J. M., Eduardo, I., Picañol, R., Batlle, I., Howad, W., Aranzana, M. J., et al. (2015). High-density mapping suggests cytoplasmic male sterility with two restorer genes in almond × peach progenies. Hortic. Res 2:15016. doi: 10.1038/hortres.2015.16
Eduardo, I., Chietera, G., Pirona, R., Pacheco, I., Troggio, M., Banchi, A., et al. (2013). Genetic dissection of aroma volatile compounds from the essential oil of peach fruit: QTL analysis and identification of candidate genes using dense SNP maps. Tree Genet. Genom. 9, 189–204. doi: 10.1007/s11295-012-0546-z
Eide, D., Broderius, M., Fett, J., and Guerinot, M. L. (1996). A novel iron-regulated metal transporter from plants identified by functional expression in yeast. Proc. Natl. Acad. Sci. U. S. A. 93, 5624–5628. doi: 10.1073/pnas.93.11.5624
Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 6:e19379. doi: 10.1371/journal.pone.0019379
Etienne, C., Rothan, C., Moing, A., Plomion, C., Bodenes, C., Svanella-Dumas, L., et al. (2002). Candidate genes and QTLs for sugar and organic acid content in peach [Prunus persica (L.) Batsch]. Theor. Appl. Genet. 105, 145–159. doi: 10.1007/s00122-001-0841-9
Felipe, A. J.. (2009). ‘Felinem', ‘Garnem', and ‘Monegro' almond x peach hybrid rootstocks. HortScience 44:196–197. doi: 10.21273/HORTSCI.44.1.196
Font i Forcada, C., Reig, G., Mestre, L., Mignard, P., Betrán, J. Á., and Moreno, M. Á. (2020). Scion × rootstock response on production, mineral composition and fruit quality under heavy-calcareous soil and hot climate. Agronomy 10:1159. doi: 10.3390/agronomy10081159
Fresnedo-Ramírez, J., Frett, T. J., Sandefur, P. J., Salgado-Rojas, A., Clark, J. R., Gasic, K., et al. (2016). QTL mapping and breeding value estimation through pedigree-based analysis of fruit size and weight in four diverse peach breeding programs. Tree Genet. Genom. 12, 1–18. doi: 10.1007/s11295-016-0985-z
Frett, T. J., Reighard, G. L., Okie, W. R., and Gasic, K. (2014). Mapping quantitative trait loci associated with blush in peach [Prunus persica (L.) Batsch]. Tree Genet. Genom. 10, 367–381. doi: 10.1007/s11295-013-0692-y
Gonzalo, M. J., Dirlewanger, E., Moreno, M. Á., and Gogorcena, Y. (2012). Genetic analysis of iron chlorosis tolerance in Prunus rootstocks. Tree Genet. Genom. 8, 943–955. doi: 10.1007/s11295-012-0474-y
Goonetilleke, S. N., March, T. J., Wirthensohn, M. G., Arús, P., Walker, A. R., and Mather, D. E. (2018). Genotyping by sequencing in almond: SNP discovery, linkage mapping, and marker design. Genes, Genom., Genet. 8, 161–172. doi: 10.1534/g3.117.300376
Guajardo, V., Solís, S., Almada, R., Saski, C., Gasic, K., and Moreno, M. Á. (2020). Genome-wide SNP identification in Prunus rootstocks germplasm collections using genotyping-by-sequencing: phylogenetic analysis, distribution of SNPs and prediction of their effect on gene function. Sci. Rep. 10, 1–14. doi: 10.1038/s41598-020-58271-5
Guajardo, V., Solís, S., Sagredo, B., Gaínza, F., Muñoz, C., Gasic, K., et al. (2015). Construction of high density sweet cherry (Prunus avium L.) linkage maps using microsatellite markers and SNPs detected by genotyping-by-sequencing (GBS). PLoS ONE 10:e0127750. doi: 10.1371/journal.pone.0127750
Hantzis, L. J., Kroh, G. E., Jahn, C. E., Cantrell, M., Peers, G., Pilon, M., et al. (2018). A program for iron economy during deficiency targets specific Fe proteins. Plant Physiol. 176, 596–610. doi: 10.1104/pp.17.01497
Hell, R., and Stephan, U. W. (2003). Iron uptake, trafficking and homeostasis in plants. Planta 216, 541–551. doi: 10.1007/s00425-002-0920-4
Ivanov, R., Brumbarova, T., and Bauer, P. (2012). Fitting into the harsh reality: regulation of iron-deficiency responses in dicotyledonous plants. Mol. Plant 5, 27–42. doi: 10.1093/mp/ssr065
Jeong, J., Cohu, C., Kerkeb, L., Pilon, M., Connolly, E. L., and Guerinot, M. L. (2008). Chloroplast Fe(III) chelate reductase activity is essential for seedling viability under iron limiting conditions Proc. Natl. Acad. Sci. U. S. A. 105:10619–10624. doi: 10.1073/pnas.0708367105
Jiménez, S., Pinochet, J., Abadía, A., Moreno, M. Á., and Gogorcena, Y. (2008). Tolerance response to iron chlorosis of Prunus selections as rootstocks. HortScience 43, 304–309. doi: 10.21273/HORTSCI.43.2.304
Kirungu, J. N., Magwanga, R. O., Shiraku, M. L., Lu, P., Mehari, T. G., Xiu, Y., et al. (2019). Genetic map construction and functional characterization of genes within the segregation distortion regions (SDRs) in the F2: 3 generation derived from wild cotton species of the D genome (No. e27845v1). PeerJ Preprints. doi: 10.7287/peerj.preprints.27845v1
Kitamura, Y., Habu, T., Yamane, H., Nishiyama, S., Kajita, K., Sobue, T., et al. (2018). Identification of QTLs controlling chilling and heat requirements for dormancy release and bud break in Japanese apricot (Prunus mume). Tree Genet. Genom. 14:33. doi: 10.1007/s11295-018-1243-3
Klagges, C., Campoy, J. A., Quero-García, J., Guzmán, A., Mansur, L., Gratacós, E., et al. (2013). Construction and comparative analyses of highly dense linkage maps of two sweet cherry intra-specific progenies of commercial cultivars. PLoS ONE 8:e54743. doi: 10.1371/journal.pone.0054743
Lambert, P., Campoy, J. A., Pacheco, I., Mauroux, J.-B., da Silva, C., Micheletti, D., et al. (2016). Identifying SNP markers tightly associated with six major genes in peach [Prunus persica (L.) Batsch] using a high-density SNP array with an objective of marker-assisted selection (MAS). Tree Genet. Genom. 12:121. doi: 10.1007/s11295-016-1080-1
Lan, P., Li, W., and Schmidt, W. (2013). A digital compendium of genes mediating the reversible phosphorylation of proteins in Fe-deficient Arabidopsis roots. Front. Plant Sci. 4:173. doi: 10.3389/fpls.2013.00173
Layne, E. E. C.. (1987). “Peach rootstocks,” in Rootstocks for Fruit Crops, eds R. C. Rom, and R. F. Carlson (New York, NY: John Wiley and Sons), 185–216.
Lu, Z. X., Sosinski, B., Reighard, G. L., Baird, W. V., and Abbott, A. G. (1998). Construction of a genetic linkage map and identification of AFLP markers for resistance to root-knot nematodes in peach rootstocks. Genome 41, 199–207. doi: 10.1139/g98-008
Mai, H. J., Pateyron, S., and Bauer, P. (2016). Iron homeostasis in Arabidopsis thaliana: transcriptomic analyses reveal novel FIT-regulated genes, iron deficiency marker genes and functional gene networks. BMC Plant Biol. 16:211. doi: 10.1186/s12870-016-0899-9
Martínez-García, P., Parfitt, D., Ogundiwin, E., Fass, J., Chan, H., Ahmad, R., et al. (2013). High density SNP mapping and QTL analysis for fruit quality characteristics in peach (Prunus persica L.). Tree Genet. Genomes 9, 19–36. doi: 10.1007/s11295-012-0522-7
Mestre, L., Reig, G., Betrán, J. A., Pinochet, J., and Moreno, M. Á. (2015). Influence of peach-almond hybrids and plum-based rootstocks on mineral nutrition and yield characteristics of ‘Big Top' nectarine in replant and heavy-calcareous soil conditions. Sci. Hort. 192, 475–481. doi: 10.1016/j.scienta.2015.05.020
Moreno, M. Á., Tabuenca, M. C., and Cambra, R. (1994). Performance of Adafuel and Adarcias as peach rootstocks. HortScience 29, 1271–1273. doi: 10.21273/HORTSCI.29.11.1271
Núñez-Lillo, G., Balladares, C., Pavez, C., Urra, C., Sanhueza, D., Vendramin, E., et al. (2019). High-density genetic map and QTL analysis of soluble solid content, maturity date, and mealiness in peach using genotyping by sequencing. Sci. Hortic. 257:108734. doi: 10.1016/j.scienta.2019.108734
Núñez-Lillo, G., Cifuentes-Esquivel, A., Troggio, M., Micheletti, D., Infante, R., Campos-Vargas, A., et al. (2015). Identification of candidate genes associated with mealiness and maturity date in peach [Prunus persica (L.) Batsch] using QTL analysis and deep sequencing. Tree Genet. Genom. 11:86. doi: 10.1007/s11295-015-0911-9
Peace, C., Bassil, N., Main, D., Ficklin, S., Rosyara, U. R., Stegmeir, T., et al. (2012). Development and evaluation of a genome-wide 6K SNP array for diploid sweet cherry and tetraploid sour cherry. PLoS ONE 7:e48305. doi: 10.1371/journal.pone.0048305
Pina, A., Irisarri, P., Errea, P., and Zhebentyayeva, T. (2021). Mapping quantitative trait loci associated with graft (in) compatibility in apricot (Prunus armeniaca L.). Front. Plant Sci. 12:45. doi: 10.3389/fpls.2021.622906
Pirona, R., Eduardo, I., Pacheco, I., Linge, C. D. S., Miculan, M., Verde, I., et al. (2013). Fine mapping and identification of a candidate gene for a major locus controlling maturity date in peach. BMC Plant Biol. 13:166. doi: 10.1186/1471-2229-13-166
Poland, J. A., Brown, P. J., Sorrells, M. E., and Jannink, J. L. (2012). Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS ONE 7:e32253. doi: 10.1371/journal.pone.0032253
Quinlan, A. R., and Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi: 10.1093/bioinformatics/btq033
R Core Team (2019). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: https://www.R-project.org/ (accessed January 10, 2022).
Rehder, A.. (1940). A Manual of Cultivated Trees and Shrubs Hardy in North America Exclusive of the Subtropical and Warmer Temperate Regions, 2nd Edn. (New York, NY: MacMillan), 162.
Reighard, G.. (2002). Current directions of peach rootstock programs worldwide. Acta Hort. 592, 421–427. doi: 10.17660/ActaHortic.2002.592.57
Robinson, N. J., Procter, C. M., Connolly, E. L., and Guerinot, M. L. (1999). A ferric-chelate reductase for iron uptake from soils. Nature 397, 694–697. doi: 10.1038/17800
Rodríguez-Celma, J., Pan, I. C., Li, W., Lan, P., Buckhout, T. J., and Schmidt, W. (2013). The transcriptional response of Arabidopsis leaves to Fe deficiency. Front. Plant Sci. 4:276. doi: 10.3389/fpls.2013.00276
Romeu, J. F., Monforte, A. J., Sánchez, G., Granell, A., García-Brunton, J., Badenes, A., et al. (2014). Quantitative trait loci affecting reproductive phenology in peach. BMC Plant Biol. 14:52. doi: 10.1186/1471-2229-14-52
Salazar, J. A., Pacheco, I., Shinya, P., Zapata, P., Silva, C., Aradhya, M., et al. (2017). Genotyping by sequencing for SNP-based linkage analysis and identification of QTLs linked to fruit quality traits in Japanese plum (Prunus salicina Lindl.). Front. Plant Sci. 8:476. doi: 10.3389/fpls.2017.00476
Santi, S., and Schmidt, W. (2009). Dissecting iron deficiency-induced proton extrusion in Arabidopsis roots. New Phytol. 183, 1072–1084. doi: 10.1111/j.1469-8137.2009.02908.x
Serra, O., Giné-Bordonaba, J., Eduardo, I., Bonany, J., Echeverria, G., Larrigaudière, C., et al. (2017). Genetic analysis of the slow-melting flesh character in peach. Tree Genet. Genom. 13:77. doi: 10.1007/s11295-017-1160-x
Sharpe, R. H., Hesse, C. O., Lownsbery, B. F., Perry, V. G., and Hansen, C. J. (1969). Breeding peaches for root-knot nematode resistance. J. Amer. Soc. Hort. Sci. 94, 209–212.
Sherman, W. B., Lyrene, P. M., and Sharpe, R. H. (1991). Flordaguard peach rootstock. HortScience 26, 427–428. doi: 10.21273/HORTSCI.26.4.427
Shirasawa, K., Hirakawa, H., and Isobe, S. (2016). Analytical workflow of double-digest restriction site-associated DNA sequencing based on empirical and in silico optimization in tomato. DNA Res. 23, 145–153. doi: 10.1093/dnares/dsw004
Solti, Á., Kovács, K., Müller, B., Vázquez, S., Hamar, É., Pham, H. D., et al. (2016). Does a voltage-sensitive outer envelope transport mechanism contribute to the chloroplast iron uptake? Planta 244, 1303–1313. doi: 10.1007/s00425-016-2586-3
Van Ooijen, J.. (2006). JoinMap 4, Software for the Calculation of Genetic Linkage Maps in Experimental Populations. Wageningen: Kyazma B.V.
Van Ooijen, J.. (2009). MapQTL 6, Software for the Mapping of Quantitative Trait Loci in Experimental Populations of Diploid Species. Wageningen: Kyazma B.V.
Verde, I., Bassil, N., Scalabrin, S., Gilmore, B., Lawley, C., Gasic, K., et al. (2012). Development and evaluation of a 9K SNP array for peach by internationally coordinated SNP detection and validation in breeding germoplasm. PLoS ONE 7:e35668. doi: 10.1371/journal.pone.0035668
Verde, I., Jenkins, J., Dondini, L., Micali, S., Pagliarani, G., Vendramin, E., et al. (2017). The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity. BMC Genom. 18:225. doi: 10.1186/s12864-017-3606-9
Voorrips, R. E.. (2002). MapChart: software for the graphical presentation of linkage maps and QTLs. J. Hered. 93, 77–78. doi: 10.1093/jhered/93.1.77
Wu, T. D., Reeder, J., Lawrence, M., Becker, G., and Brauer, M. J. (2016). “GMAP and GSNAP for genomic sequence alignment: enhancements to speed, accuracy, and functionality,” in Statistical Genomics. Methods in Molecular Biology, vol 1,418, eds E. Mathé, and S. Davis (New York, NY: Humana Press), 15. doi: 10.1007/978-1-4939-3578-9_15
Xu, S.. (2008). Quantitative trait locus mapping can benefit from segregation distortion. Genetics 180, 2201–2208. doi: 10.1534/genetics.108.090688
Yang, N., Reighard, G., Ritchie, D., Okie, W., and Gasic, K. (2013). Mapping quantitative trait loci associated with resistance to bacterial spot (Xanthomonas arboricola pv. pruni) in peach. Tree Genet. Genom. 9, 573–586. doi: 10.1007/s11295-012-0580-x
Yang, T. J., Lin, W. D., and Schmidt, W. (2010). Transcriptional profiling of the Arabidopsis iron deficiency response reveals conserved transition metal homeostasis networks. Plant Physiol. 152, 2130–2141. doi: 10.1104/pp.109.152728
Zarrouk, O., Gogorcena, Y., Gómez-Aparisi, J., Betrán, J. Á., and Moreno, M. Á. (2005). Influence of almond x peach hybrids rootstocks on flower and leaf mineral concentration, yield and vigour of two peach cultivars. Sci. Hortic. 106, 502–514. doi: 10.1016/j.scienta.2005.04.011
Zeballos, J. L., Abidi, W., Giménez, R., Monforte, A. J., Moreno, M. Á., and Gogorcena, Y. (2016). Mapping QTLs associated with fruit quality traits in peach [Prunus persica (L.) Batsch] using SNP maps. Tree Genet. Genom. 12:37. doi: 10.1007/s11295-016-0996-9
Zhang, D. Q., Zhang, Z. Y., and Yang, K. (2007). Genome-wide search for segregation distortion loci associated with the expression of complex traits in Populus tomentosa. Forest. Stud. China 9, 1–6. doi: 10.1007/s11632-007-0001-7
Zhang, L., Wang, S., Li, H., Deng, Q., Zheng, A., Li, S., et al. (2010). Effects of missing marker and segregation distortion on QTL mapping in F2 populations. Theor. Appl. Genet. 121, 1071–1082. doi: 10.1007/s00122-010-1372-z
Keywords: molecular markers, rootstock, stone fruits, segregating population, reference genome, pseudomolecule, SPAD
Citation: Guajardo V, Martínez-García PJ, Solís S, Calleja-Satrustegui A, Saski C and Moreno MÁ (2022) QTLs Identification for Iron Chlorosis in a Segregating Peach–Almond Progeny Through Double-Digest Sequence-Based Genotyping (SBG). Front. Plant Sci. 13:872208. doi: 10.3389/fpls.2022.872208
Received: 09 February 2022; Accepted: 25 April 2022;
Published: 31 May 2022.
Edited by:
Sezai Ercisli, Atatürk University, TurkeyReviewed by:
Elisa Vendramin, Council for Agricultural Research and Economics, ItalyJose Xavier Chaparro, University of Florida, United States
Copyright © 2022 Guajardo, Martínez-García, Solís, Calleja-Satrustegui, Saski and Moreno. 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: María Ángeles Moreno, mmoreno@eead.csic.es
†These authors have contributed equally to this work and share first authorship