Elucidation and characterization of QTLs for Russet formation on apple fruit in ‘Honeycrisp’-derived breeding germplasm

Russet formation in apples (Malus domestica Borkh.) is a superficial skin disorder that detracts from fruit appearance and therefore is a breeding target in many apple scion breeding programs. Previous studies have reported quantitative trait loci (QTLs) for russet formation, but these results have been limited to a narrow genetic background, which reduces breeding utility. In this study, 14 breeding parents were represented in 13-full sib families, nine of which were ‘Honeycrisp’-derived while the other four were derived from ‘Minneiska’, an offspring of ‘Honeycrisp.’ Two previously reported QTLs, on linkage group (LG) 2 and LG6, were confirmed and characterized for the first time in ‘Honeycrisp’-derived germplasm. Five other putative QTLs were also detected. The QTLs characterized in this study could be targeted for DNA test development to be used in marker-assisted selection. Non-additive interactions were observed at and between the LG2 and LG6 QTLs. Future work should include validation of QTL and haplotype effects in other germplasm, development of trait-predictive DNA tests for use in apple breeding programs, and exploration of genome-wide prediction models for russet formation.


Introduction
Fruit appearance is an important aspect of apple (Malus domestica Borkh.) fruit quality. Fruits with poor appearance are often downgraded from fresh market to processing, which can result in a three-to-four-fold loss of wholesale income (USDA National Agricultural Statistics Service 2021). Russet formation is one contributor to poor appearance in apple fruit that can reduce grower packout and therefore is an important target in many apple scion breeding programs (Baugher et al. 1996).
Russet formation is caused by microcracking of the fruit skin and is hypothesized to be the result of the production and transportation of suberin, a cork-like substance, to the outer dermal layers of the apple fruit (Faust and Shear 1972).
Microcracking can be induced by several abiotic and biotic factors such as high humidity, frost damage, fungal pathogens, and management practices (Faust and Shear 1972). Several apple cultivars (e.g. 'Golden Russet', 'Chestnut Crab', and 'Egremont Russet') are prone to severe russet formation while other cultivars will only develop russet under conducive environmental conditions.
Phenotyping russet formation can be challenging due to the strong influences of environmental conditions on russet formation and challenges associated with subjective visual ratings. Russet formation is strongly influenced by environmental conditions (e.g., humidity) and, therefore, can vary from year to year and even tree to tree within a given year (Faust and Shear 1972;Sugar et al. 2005;Eccher and Hajnajari 2006). Typically, breeders estimate the amount of russet formation with visual ratings. Differing fruit background colors and/or shapes can make russet formation appear more or less prominent. Previous studies have attempted to address these challenges by using binary or categorical bins to classify russet formation or using labor intensive techniques to obtain precise russet formation measurements (Scharwies et al. 2014;Falginella et al. 2015;Lashbrooke et al. 2015;Kunihisa et al. 2016;Takeuchi et al. 2021). Communicated by D. Chagné Several quantitative trait loci (QTLs) associated with russet formation have been previously reported for apple and pear (Pyrus), a genus having a high degree of genomic colinearity with Malus (Celton et al. 2009). In apple, russet formation QTLs have been reported on linkage groups (LGs) 2, 4, 6, 8, 10, 11, 12, and 15 (Kunihisa et al. 2014;Falginella et al. 2015;Lashbrooke et al. 2015). Two major QTLs have been reported to play important roles in russet formation in Japanese (Pyrus pyrifolia (Burm. F.) Nakai) and Chinese (P. ussuriensis Maxim., P. bretschneideri Rehder) pear species. The R locus, a dominant locus identified on LG8 of Japanese pear, has been reported in several studies to be associated with russet formation. The I locus, a modifier genetic factor, is associated with russet severity but has not been mapped. In addition, studies have reported other small-effect QTLs associated with russet formation in Japanese pear (Wang et al. 2016;Takeuchi et al. 2021). Although russet formation is an attractive target for DNA-informed breeding (e.g., marker-assisted seedling selection), breeders are limited by the few DNA tests that are currently available.
Current understanding of the genetic basis of russet formation in apple is insufficient and has to-date been limited to a relatively narrow set of germplasm. The important modern US cultivar and breeding parent 'Honeycrisp' was not represented in previous studies and russet formation has been commonly observed in its ancestors, 'Frostbite,' 'Keepsake,' and 'Golden Delicious.' Russet formation has also been observed in germplasm derived from 'Honeycrisp' in the University of Minnesota (UMN) apple breeding program, including on fruit of 'Minneiska' (SweeTango® apple). There currently are no predictive markers available for DNA-informed breeding for russet formation on apples in this germplasm. To address this limitation, QTLs associated with russet formation must be detected and/or validated in relevant breeding germplasm, QTL alleles need to be characterized, and predictive DNA tests for russet formation QTLs need to be developed. The goal of this study was to examine the genetic basis of russet formation in a pedigree-connected germplasm set relevant to apple scion breeding programs that utilize 'Honeycrisp'-derived germplasm. The objectives of this study were to (1) identify QTLs associated with russet formation and (2) characterize QTL allele effects. We hypothesized that multiple QTLs for russet formation would be detected.

Germplasm
Thirteen pedigree-connected full-sib families (n = 1009 offspring) were evaluated in this study (Table 1). These thirteen families were grouped into two populations. The first population consisted of individuals from nine families with 'Honeycrisp' as a common parent. The other parents in this population were 'Jonafree,' 'Pitmaston Pine Apple,' 'Minnewashta,' 'WA 2,' and Arkansas selection AA44, which is also known as 'MonArk' but will be referred to as AA44 throughout, and the UMN advanced selections MN1702, MN1764, MN1836, and MN1915. The second population consisted of individuals from four families with 'Minneiska' ('Honeycrisp' × 'Minnewashta') as a common parent. The other parents in this population were 'Wildung,' 'MN55,' and the UMN advanced selections MN1702 and MN1965.

Phenotypic data
Individuals were evaluated for russet formation from 2017 to 2019 at the UMN Horticultural Research Center in Chanhassen, MN, USA, with each individual being evaluated in up to 3 years. Trees were grown on B9 and/or G16 rootstock. Three to 10 fruit per individual were harvested within a given year from August to October. Differences in the number of harvested individuals and families across years were in part due to biennial bearing and juvenility of certain families that had harvestable fruit in the later years (Table 1).

Genotypic data
Individuals were genotyped via the International ROS-BREED SNP Consortium 8 K Illumina Infinium® array v1 or the Illumina Infinium 20 K array (Chagné et al. 2012;Bianco et al. 2014) with DNA extracted as described by Clark et al. (2014). SNP data were filtered to include SNPs common to both arrays for a total of 2213 SNPs across 17 linkage groups (Vanderzande et al. 2018). Genetic positions from the iGL genetic map (Di Pierro et al. 2016) updated as described in Howard et al. (2021) were used for the SNPs included in this study.

QTL analyses
QTL analyses were conducted with FlexQTL™ software (www. flexq tl. nl), a pedigree-based QTL analysis software that uses Markov Chain Monte Carlo (MCMC) simulation (Bink et al. 2002(Bink et al. , 2008(Bink et al. , 2012. Parameters for analyses are reported in Supplementary Table S1, and MCMC convergence was considered reached if effective chain size lengths for the estimates of trait means, numbers of QTLs, and their variances reached at least 100 (as used in Bink et al. 2014;Kostick et al. 2021). Phenotypic data were only included for unselected offspring to avoid selection bias (as in Howard et al. 2018;Verma et al. 2019;Kostick et al. 2021). QTL significance was determined using Bayes factors (BF) and posterior intensity values. Evidence for QTLs on individual linkage groups was considered none (< 2), positive (2 to < 5), strong (5 to < 10), or decisive (10 ≤) based upon the Bayes factor (2lnBF10) range given in the parentheses (Kass and Raftery 1995). QTL regions for QTLs with strong to decisive evidence were determined by the same method reported in Howard et al. (2018). QTL intervals were defined in a series of successive 2 cM "bins," chromosomal segments, that had a BF above 5 and the mode within a given QTL region was considered the most probable QTL position. The proportion of phenotypic variance explained by a QTL was estimated by dividing the reported variance explained by the whole phenotypic variance provided by FlexQTL™ software (Verma et al. 2019). QTL intervals were then compared to previous literature alongside the GDDH13v1.1 reference genome (Daccord et al. 2017) to determine if multiple QTLs appear to colocalize in the same region.

Haplotype characterization
QTLs were targeted for haplotype characterization if they had decisive evidence (BF > 10) for the presence of a QTL. SNPs included for haplotyping were chosen based on their proximity to QTL peaks. Proximity was determined by including SNPs that fell in the 2 cM bins on either side of the QTL peak were included with redundant markers, based upon parental haplotypes removed. Marker phasing was performed through FlexQTL™ software. Haplotypes were traced through extended pedigrees to the furthest known ancestors using techniques described in Howard et al. (2018), and employed in Luby et al. (2022), to reconstruct pedigrees in the UMN apple breeding program. Additionally, unpublished extended pedigree information for 'WA 2' generated from the project described by Howard et al. (2018) was used. Identical haplotypes were described as identicalby-descent (IBD) if they could be traced to a common ancestor whereas identical haplotypes that could not be traced to a common ancestor were considered identical-by-state (IBS). When necessary, imputation of missing SNP information was performed via an examination of progeny and ancestors, as described by Kostick et al. (2021). Recombinants and haplotypes that were unable to be resolved were excluded from further statistical analyses.
For each QTL, analysis of variance (ANOVA) was used to examine if parental haplotypes varied significantly for russet rating in the offspring. Using mean separation via Tukey's honest significant difference (HSD) test, haplotypes at a given QTL were grouped into relative functional allele groups based on their relative effects (i.e., low, moderate (non-significant), and high relative levels of russet ratings). Functional haplotype effects were determined using a combination of Tukey's HSD results and haplotype russet rating mean values and ranges.
To examine the effects of haplotypes across QTLs, individuals were grouped by functional (i.e., low, moderate, high) QTL genotypes across the targeted QTLs, and a two-way ANOVA was used to determine if compound QTL genotypes differed significantly for russet ratings.

Phenotypic data
A total of 60% of individuals were harvested in at least 2 years. Raw family means, standard deviations, and ranges for BLUP-adjusted shoulder and lenticel russet data are presented in Table 2, while all other trait data by family are presented in Supplementary Table 2. Russet rating distributions were often skewed towards low russet ratings. Overall and body russet ratings were highly correlated (r ≥ 0.8 to 0.9) across years and had means ranging from 1.2 to 1.8 for overall and 1.1 to 1.6 for body russet ratings. Shoulder russet ratings were also highly to moderately correlated with overall and body russet ratings (r ≥ 0.4 to 0.7) across years and had means ranging from 1.5 to 2.9 russet ratings (Table 2, Supplementary  Table S3). Year had a significant effect (P < 0.001) on shoulder and lenticel russet formation traits (Supplementary Table S4). In 2018, incidence of russet formation was low (2018 shoulder russet rating mean = 1.6, mean overall years = 2.0). Mean shoulder russet ratings for families were within two points of the population mean and were at the lower end of the russet rating scale. Despite skewed distributions, family effects were significant (p < 0.001, Supplementary Table S4). Families had standard deviations that varied from 0.46 to 2.05 and had ranges within families from 1-3.5 to 1-10 for shoulder russet formation ( Table 2). Minimal phenotypic variation for stem, body, calyx, and overall russet ratings was observed within and among families. Therefore, QTL analyses were not carried out for stem, body, calyx, and overall russet ratings. Quantitative variation among and within families for shoulder russet ratings was observed although ratings were skewed towards low values. Lenticel russet ratings were evenly distributed over the entire range of the rating scale (Table 2).

QTL detection
Two QTLs on LG2 and LG6 were detected with decisive (Bayes Factor; BF; 2lnBF 10 > 10) evidence and five QTLs were detected with strong (BF > 5) evidence for shoulder russet formation (Table 3). QTL intervals with decisive evidence were located at 48-54 cM and 57-69 cM for LG2 and LG6, respectively. FlexQTL™ software runs with lenticel russet rating data failed to reach convergence and, thus, failed to detect large-effect QTLs.
All LG6 QTL haplotypes are described with their ancestral origins in Supplementary Table S8

Combined LG2 and LG6 QTL effects
Phenotypic variances explained by the LG2, LG6, and the two QTLs combined were 4%, 8%, and 12%, respectively (Table 3). Together, LG2 and LG6 QTLs had significant effects on shoulder russet formation (Supplementary  Table S4). In general, the higher the number of high shoulder russet formation haplotypes an individual had across QTLs, the higher the shoulder russet rating (Fig. 3). However, nonadditive interactions among haplotypes were observed. The presence of at least two low shoulder russet formation haplotypes, at either or both loci, resulted in lower shoulder russet ratings (Fig. 3). In the presence of zero or one low shoulder russet formation haplotype, the mean shoulder russet ratings for offspring with one (n = 74), two (n = 247), three (n = 112), or four (n = 5) high shoulder russet formation Numbers of individuals with one or more copies of each haplotype are listed above the boxplot. Total number of individuals represented in the LG 2 haplotype figure (n = 984) and LG 6 haplotype figure (n = 997). Individuals with unresolved or recombinant haplotypes were not included in distributions

Discussion
This is the first report of QTLs associated with shoulder russet formation in 'Honeycrisp'-derived families and should be of interest to numerous apple scion breeders now using 'Honeycrisp' and its descendants in their breeding programs. Two QTLs for shoulder russet formation were characterized on LG2 and LG6 in a 'Honeycrisp'-derived germplasm set. A total of eight high shoulder russet formation haplotypes and six low shoulder russet formation haplotypes were characterized. Findings of this study could help inform breeding decisions and/or be used to develop DNA tests for breeding use.

Quantitative variation for russet formation observed within germplasm set
Most russet formation occurred within the shoulder region or on the lenticels, and therefore these two russet ratings were examined in this study. Shoulder russet ratings were highly correlated with overall russet ratings, a common russet formation metric used in other studies (Kunihisa et al. 2014;Falginella et al. 2015;Lashbrooke et al. 2015). Therefore, use of a single russet formation metric is likely sufficient for russet formation QTL detection. Previous studies that investigated the genetic architecture of russet formation captured higher levels of russet formation within their populations but relied upon biparental populations having at least one parent being considered susceptible to russet formation (Falginella et al. 2015;Lashbrooke et al. 2015). In contrast, most of the parents of families evaluated in the present study are not known for high levels of russet formation, but high shoulder russet formation of individual offspring has been observed within the UMN breeding program. The variation of shoulder russet ratings observed in this study, both within and among families, and the identification of two large effect QTLs as well as the multiple additional putative QTLs suggested multi-genic control for shoulder russet formation in this germplasm.
Of the other putative russet formation QTLs detected in this study, both the LG15 QTL  and the LG12 QTL (25.94-31.30 Mbp) had previously reported QTLs near the QTL region. Lashbrooke et al. (2015) reported a QTL for tensile strength of the fruit skin (~ 12 Mbp) which was used to indirectly measure susceptibility to form russet. Falginella et al. (2015) reported a major QTL for increased levels (> 25% overall coverage) of russet formation with a QTL peak at ~ 28.9 Mbp.

Linkage groups 2 and 6 shoulder russet formation QTLs colocalized with important disease resistance QTLs
The LG2 and LG6 shoulder russet formation QTLs detected in this study colocalized with disease resistance QTLs. The LG2 apple scab resistance QTL (~ 31.1 Mbp) reported by Calenge et al. (2004), colocalized with the LG2 russet formation QTL detected in this study (Calenge et al. 2004;Galli et al. 2010;Bus et al. 2011). The LG6 russet formation QTL in this study colocalized with a fire blight resistance/ susceptibility QTL (31.22-36.72 Mbp) reported by Kostick et al. (2021). To our knowledge, there is no established relationship in the literature between russet formation and either of these diseases. The association between disease resistance traits and the russet formation QTLs could be examined in future studies. Colocalization of multiple trait loci demonstrates the importance of considering these regions when making breeding decisions. Breeders could leverage knowledge of the different allele effects at these trait loci to develop high quality apple cultivars more efficiently.

Unique high and low shoulder russet formation haplotypes at linkage group 2 QTL in 'Minneiska'
This is the first report of a QTL associated with shoulder russet formation in 'Minneiska' derived families. Haplotype analysis indicated that 'Minneiska' was heterozygous for the LG2 QTL. Found in 'Minneiska' and its progenitors, haplotype 2 J was significantly associated with reduced levels of shoulder russet ratings. Haplotype 2 J traced back to 'Tetofsky,' an old Russian cultivar that is described as having thicker, waxy skin, although russet formation is not mentioned (Beach et al. 1905). 'Minneiska' is described as having variable levels of russet formation, starting at the stem end with the potential to have extensive netting spreading out and down the shoulder area (Bedford and Luby 2008). Offspring of 'Minneiska' crosses also varied in their range of shoulder russet ratings depending on the other parent. In general, when 2 J was inherited from 'Minneiska,' or its parent 'Minnewashta,' and paired with other low shoulder russet formation haplotypes, offspring had low shoulder russet ratings. The other haplotype of 'Minneiska' (2A) was represented in all families and could be traced back to 'Northern Spy'; all 2A haplotypes, except from MN1764, were inherited through 'Honeycrisp' directly or via 'Minneiska.' Haplotype 2A was significantly associated with increased levels of shoulder russet ratings.

High shoulder russet formation haplotypes at linkage group 2 QTL in 'Pitmaston Pine Apple'
The presence of multiple putative high shoulder russet formation haplotypes, 2F and 2L, in 'Pitmaston Pine Apple' was not surprising, as 'Pitmaston Pine Apple' is a full russet-type apple cultivar (Bussey 2016). Offspring of 'Pitmaston Pine Apple' had an average shoulder russet rating of 2.9, the highest family average for shoulder russet rating along with 'Honeycrisp' × 'WA 2.' Due to limited pedigree information, only one haplotype (2F) could be traced back through an extended pedigree to 'Reinette Franche,' an ancestor that is prone to russet formation on the stem and calyx end of its fruit (Bussey 2016). 'Reinette Franche' is an important ancestor in the background of many cultivars today (Muranty et al. 2020); therefore, the high shoulder russet formation haplotype 2F might segregate in other breeding germplasm. The effects of haplotypes 2F and 2L should be considered putative and targeted for validation in future studies due to relatively low representation in this study (n = 60).

Unique high shoulder russet formation haplotype at linkage group 6 QTL that colocalizes with a fire blight QTL in 'Honeycrisp'
'Honeycrisp' was heterozygous at the LG6 russet formation QTL based on haplotype analysis. High shoulder russet formation haplotype 6F, found in 'Honeycrisp,' 'Minneiska,' 'MN55,' MN1702, MN1764, and MN1965, could be traced back to 'Frostbite,' a cultivar that has been reported to be prone to russet formation (Bussey 2016). This QTL region was reported by Kostick et al. (2021) to be associated with fire blight susceptibility. The 'Honeycrisp' fire blight allele that traced back to 'Frostbite' was associated with a reduced susceptibility to fire blight. Knowledge of this potential association might be important for breeders who are breeding for apples with low shoulder russet formation that are also fire blight resistant. To our knowledge, the 6F high shoulder russet formation haplotype is unique to germplasm derived from 'Frostbite.'

Very high shoulder russet formation haplotype at linkage group 6 QTL in 'WA 2'
The presence of a putative very high shoulder russet formation haplotype, 6A, only found in the 'WA 2' cross, was not surprising. While 'WA 2' is described as not developing russet in commercial growing areas of Washington State (Evans et al. 2010), shoulder russet formation has been noted when 'WA 2' is grown in Europe (Evans, personal communication). Offspring of 'WA 2' (n = 98) had an average shoulder russet rating of 2.9, and those having a 6A haplotype (n = 27) had an average rating of 3.5. Although based on mean separation analyses, haplotype 6A was significantly different from the other high shoulder russet formation haplotypes (Supplementary Table S8), due to the lack of representation in multiple families, we categorized it within the same group as other high haplotypes. Haplotype 6A was traced back to 'Winter White Pearmain,' an old cultivar of unknown origin that is reported to sometimes have lenticel russet formation and russet formation within the stem cavity (Bussey 2016). Other potentially important and breeding relevant ancestors that have inherited high shoulder russet formation haplotype 6A from 'Winter White Pearmain' are as follows: 'Delicious,' 'Gala, ' 'Kidd's Orange Red,' and 'Splendour.' Interactions at and among shoulder russet formation QTLs were non-additive Unequal representation of haplotypes at and across QTLs made assessment of dominance and/or epistasis complicated. By using an approach of assigning haplotype functional effects and reducing groups into "HIGH," "MOD" (moderate), and "LOW," statistical power was increased.
Effects across multiple QTLs have not been previously reported for russet formation in apple. Non-additive effects have been reported within QTL regions (Lashbrooke et al. 2015) and were also seen in this study (Fig. 3). When two low shoulder russet formation haplotypes were present, either within or across QTLs, shoulder russet formation was relatively low. However, further investigation is needed to examine non-additive QTL interactions.

Study limitations
A limitation of this study was the reliance on natural conditions for russet formation development. Russet formation is largely dependent on spring environmental conditions (Barceló-Vidal et al. 2013). In this study, substantial yearto-year variation in spring and summer weather was potentially reflected in the russet rating data. In 2017, a relatively wet spring and early summer might have resulted in more variation among and within families for russet formation and summer hail events reduced the number of harvestable fruit available. In 2018, little variation for russet formation within and among families was observed, which was likely due to dry environmental conditions. In 2019, a relatively wet spring but dry early summer with no severe hail events resulted in a large number of harvestable fruit with an increased incidence of russet formation. By harvesting individuals across years and using adjusted average russet rating data for QTL analysis, a large (n = 1009 individuals) representative germplasm set was used in the analyses.
Small family sizes likely limited the detection of QTLs, characterization of certain haplotypes, and examination of QTL interactions. Uneven representation of some parent genomes, a common challenge in pedigree-based QTL analyses (e.g., Bink et al. 2014;Verma et al. 2019;Kostick et al. 2021;Kostick and Luby 2022), might have led to limited or no representation of russet formation QTLs haplotypes and compound QTL genotypes. However, even with limited representation of certain haplotypes, the haplotypes characterized in this study could enable targeted crosses to be made for future validation studies.
An assumption of this study is that IBS haplotypes have the same functional effects. By looking at haplotype effects within each family, this assumption appears to be justified for most haplotypes represented in this study. Two exceptions were haplotypes 2D and 6C. Haplotype 2D was assigned a high functional effect and was represented in two families, 'Honeycrisp' × 'WA 2' and 'Honeycrisp' × 'Jonafree.' In 'Honeycrisp' × 'WA 2,' the 2D haplotype appeared to act as a high shoulder russet formation haplotype but in the 'Honeycrisp' × 'Jonafree' family, the 2D haplotype had the lowest mean shoulder russet rating for LG2 haplotypes. A similar inconsistency was seen with haplotype 6C assigned as a high shoulder russet formation haplotype and present in 'Honeycrisp' × 'Pitmaston Pine Apple' and 'Minneiska' × 'Wildung.' In 'Honeycrisp' × 'Pitmaston Pine Apple,' 6C had the highest mean shoulder russet rating for LG6 haplotypes within this family, while in 'Minneiska' × 'Wildung,' 6C had the lowest mean shoulder russet rating for LG6 haplotypes. In both cases, IBS haplotypes 2D and 6C could not be traced back to a common ancestor and extended haplotypes did not match outside the QTL range. One possible explanation for the differences in mean haplotype values could be the presence of low or high shoulder russet formation alleles at the other QTLs or differences in genomic background. Another explanation is that these IBS haplotypes could be functionally different. Validation of the effects of these haplotypes within a wider germplasm set could clarify if they are functionally different.

Breeding implications
The complexity of the genetic control of russet formation will continue to make development of cultivars with low levels of russet formation challenging. The QTLs characterized in this study accounted for a moderate proportion of variation in the 'Honeycrisp'-derived germplasm. The LG2 and LG6 QTL allele effects should be validated in other germplasm before targeting these loci for development and deployment of trait-predictive DNA tests.
Low and high shoulder russet formation haplotype information from this study can be used to enable selection for low shoulder russet formation levels in breeding. In the shortterm, haplotype information on the individuals characterized in this study could be used for parent selection by selecting for parents with multiple low shoulder russet formation alleles. Shoulder russet formation QTL haplotype information could also be used to develop trait-predictive DNA tests for one or both QTLs. Breeders will need to consider the colocalization of other trait loci (i.e., apple scab resistance/ susceptibility; fire blight resistance/susceptibility) when targeting for russet formation. For example, the 'Honeycrisp' 6F, a high shoulder russet formation haplotype, derived from 'Frostbite' appears to be in repulsion with the reduced susceptibility for fire blight, as reported by Kostick et al. (2021). While the current fruit quality standard for fresh market apples favors low levels of russet formation, this information could also be used by apple breeders interested in developing new cultivars with higher levels of russet.

Conclusions
In this study, two previously reported QTLs, on LG2 and LG6, were confirmed and characterized for the first time in "Honeycrisp"-derived germplasm. Five other putative QTLs were also detected. Non-additive interactions were observed at and across LGs 2 and 6. The QTLs characterized in this study can be targeted for DNA test development to be used in marker-assisted selection. The information presented in this study highlights the necessity for breeders to consider haplotypes for russet formation across loci. Future work should include validation of QTL and haplotype effects in other germplasm, development of trait-predictive DNA tests for use in apple breeding programs, and exploration of genomewide prediction models for russet formation.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.