Abstract
We estimated environmental variance components for fruit weight (FW), ease of peeling (EP), firmness of segment membrane (FSM), soluble solids content (SSC), and acidity of 20 mandarin and related cultivars used as cross-parents in citrus breeding in Japan. The variance by year was largest for FW (38% of the total) and SSC (23%), but negligible for EP, FSM, and acidity. Variance among trees and tree × year interaction were negligibly small except for FW, whereas the genotype × year interaction variance ranged from 9% to 22% of the total. Variance among fruits from the same tree was the largest environmental component for EP and FSM. Broad-sense heritability (hB2) was 0.29 for FW, 0.52 for EP, 0.11 for FSM, 0.34 for SSC, and 0.65 for acidity in evaluations performed using single-year measurements of one tree with one fruit. Yearly repetition was more efficient than tree replication for increasing hB2 of all traits investigated. Increasing fruit number from one to five in combination with yearly repetition also increased hB2, especially for EP and FSM.
Satsuma mandarin (Citrus unshiu), the leading citrus crop in Japan, is of Japanese origin and is well adapted to the local climate. It has high productivity, it is disease- and cold-tolerant, and its fruits are easy to eat, being seedless and easy to peel and having soft segment membranes, especially in the early-maturing cultivars. However, the fruits are small and their sugar content is moderate. The consumer demand for Satsuma mandarin has gradually decreased, and current production is only one-fourth of that in 1975. To maintain and develop the citrus industry in Japan, it is indispensable to develop new mandarin cultivars with large fruits and high sugar content in addition to ease of eating.
The citrus breeding program in Japan began in 1937 at the Horticulture Research Station of the Ministry of Agriculture and Forestry, now the Okitsu Citrus Research Station of the National Agriculture and Food Research Organization Institute of Fruit Tree Science (NIFTS), and is continuing at the NIFTS Okitsu and Kuchinotsu Citrus Research Stations (Nishiura, 1964). The breeding targets in this program currently focused on fruit traits, such as FW, high SSC, and ease of eating that is determined mainly by EP and FSM. In addition, it is important to broaden the range of maturation time to obtain cultivars adapted to a wide range of growing seasons; the availability of such cultivars would improve profitability and allow optimization of labor allocation. Acidity, which decreases as fruit matures, is one of the main factors that determine fruit maturing time in citrus. These fruit traits would fluctuate owing to environmental (nongenetic) factors such as year, tree, and fruits. Therefore, it is important for fruit tree breeder to distinguish genetic variance from environmental ones. Increasing year repetition, tree replication, and fruit sample size increase the accuracy of genetic differences, but it restricts number of genotypes in a given period and space and reduces the chance of obtaining desirable genotypes in a breeding program.
Analysis of variance (ANOVA) has been used to estimate the contributions of genetic and environmental variance components for fruit traits of Japanese persimmon (Yamada et al., 1993, 2002), grape (Sato et al., 2000), Japanese pear (Machida and Kozaki, 1975, 1976; Nishio et al., 2011), chestnut (Nishio et al., 2014) and for carotenoid contents of citrus (Nonaka et al., 2012). These estimates provide information on optimal yearly repetition, tree or vine replication, and fruit sample size. Estimates of hB2 of traits, which is the ratio of genetic variance to phenotypic variance, are also important for breeders to judge whether a trait could be effectively selected in a breeding program. However, in the citrus breeding program at NIFTS, the availability of such estimates for commercially important fruit traits is limited. Therefore, the objective of the present study was to obtain accurate estimates of the environmental variance components for commercially important fruit traits (FW, EP, FSM, SSC, and acidity) in December, the main season for fruit evaluation in breeding program in Japan, using 20 mandarin and related cultivars that are used as cross-parents in citrus breeding at the Kuchinotsu Citrus Research Station.
Materials and Methods
Measurement and evaluation of fruit traits.
Twenty cultivars (genotypes) used in the study (Table 1) were grown at the NIFTS Kuchinotsu Citrus Research Station in Nagasaki, Japan (32°59′N, 130°18′E). The age of the trees ranged from 6 to 18 years in 2012, top grafted on C. unshiu that had been grafted on trifoliate orange. The trees were grown using normal cultivation techniques used in commercial production in Japan. Two trees per cultivar were repeatedly measured and evaluated in 2012 and 2013. Five fruits were randomly sampled from each tree in mid-December, and we evaluated FW, EP, FSM, SSC, and acidity. FW was measured on an electronic digital scale (CP8201; Sartorius, Gottingen, Germany). EP was evaluated by peeling fruits by hand and rated on a 1 to 5 scale by one person, where 1 = fingers can penetrate the rind of fruit apex where the peeling is started with little pressure and can peel it without tearing it off; 3 = fingers can penetrate the rind but cannot peel it without tearing it off; 5 = fingers cannot penetrate the rind; and 2 and 4 are intermediates.
Mandarin and related cultivars and their origins.z


Fruits were cut in half at the equator line, and segments were evaluated in an eating test. FSM was rated on a 1 to 5 scale by one person, where 1 = so soft that the membrane cannot be detected during eating; 3 = the membrane can be detected but can be swallowed easily; 5 = so hard that the membrane cannot be swallowed; and 2 and 4 are intermediates.
Squeezed juice was used to measure SSC and acidity with an acid and sugar concentration analyzer (NH-2000; Horiba, Kyoto, Japan). Average values of each were positively correlated with sd of FW and acidity; therefore, values of FW and acidity were log transformed for ANOVA to improve the normality of the distribution of the residual estimates.
Statistical analysis.


The distribution of the residuals in the ANOVA model approached normal for all traits, with no significance at P = 0.05 in the Kolmogorov–Smirnov one-sample test.




Results and Discussion
Means of the population were 197.8 g for FW, 2.3 for EP, 1.9 for FSM, 11.2% for SSC, and 1.2 g/100 mL for acidity. Means of the genotypes for both years ranged from 115 to 287 g for FW, 1.1 to 4.1 for EP, 1.3 to 2.6 for FSM, 9.4% to 13.2% for SSC, and 0.6 to 2.1 g/100 mL for acidity.
In ANOVA, the effect of genotype was highly significant (P ≤ 0.01 except for FSM), and the effect of genotype × year interaction was also highly significant (P ≤ 0.01 except for FW) (Table 2). The effect of year was highly significant (P ≤ 0.01) for FW and SSC. The effect of tree within genotype was significant (P ≤ 0.05) for EP and acidity, and the effect of the tree × year interaction was highly significant (P ≤ 0.01) for FW, SSC, and acidity.
Analysis of variance of fruit traits using 20 genotypes (two trees per genotype) for 2 years.


The contribution of the environmental variance components to the total variance varied with the traits (Table 3). The contribution of genetic variance (σg2) was the largest (65%) for acidity among the five traits, indicating that it would be easy to accurately estimate genetic differences for this trait. On the other hand, σg2 values were small for FW (29%) and for FSM (11%), indicating that it would be difficult to estimate genetic differences for these traits. As for FSM, the small value of σg2 could also result from the small difference among genotypes (from 1.3 to 2.6 in the mean) used in the study. On the other hand, σg2 values were the largest (65%) for acidity among the five traits, indicating that this genetic trait would have wide range among genotypes and would be distinguished with less fluctuation by environmental factors. The contribution of variance by year (σy2) was large for FW (38%), but negligibly small for EP, FSM, and acidity. The contribution of the genotype × year interaction (σgy2) ranged from 9% to 22% of the total variance; the σgy2 values were largest (>20%) for FSM and SSC. The large σy2 and small σgy2 suggest that FW would be affected by the year similarly in different genotypes, i.e., there would be years with large fruits and years with small fruits, which could be caused by biennial bearing in on and off years. On the other hand, the large σgy2 and small σy2 suggest that FSM would be affected by the year differently in different genotypes, i.e., FSM of some genotypes would be lower and that of the others would be higher than the average of genotypes on the same year. The variances by tree (σt2 and σty2) were small except for FW (σty2 = 14%). These results suggest that increasing the number of trees is not effective to increase the accuracy of the genetic properties. The contribution of the fruits within tree (σf2) was the largest among environmental variance components for EP and especially for FSM (61%), indicating that EP and FSM would vary considerably among fruits from the same tree, probably caused by differences in fruit development, maturation, and location within a tree. In addition, as sensory evaluations for these traits were subjective, it would also vary among evaluators though no consideration in this study. Therefore, the variance caused by the evaluators should be examined in future studies.
Estimates of variance components and their contributions to the total variance obtained by analysis of variance of fruit traits using 20 genotypes (two trees per genotype) for 2 years.


To assess how yearly repetition, tree replication, and fruit sample size affect the efficiency of the evaluation of these fruit traits, we calculated the hB2 of each traits (Fig. 1). In the case of one tree with one fruit per tree and single-year measurements, hB2 was highest for acidity and lowest for FSM. Increasing the number of fruits on a single tree to five with single-year measurements increased hB2 from 0.29 to 0.32 for FW, 0.52 to 0.69 for EP, 0.11 to 0.23 for FSM, 0.34 to 0.38 for SSC, and 0.65 to 0.71 for acidity. Further increasing the number of fruits on a single tree to 10 with single-year measurements resulted in only small or no increase in hB2: no change for FW and an increase from 0.69 to 0.72 for EP, 0.23 to 0.26 for FSM, 0.38 to 0.39 for SSC, and 0.71 to 0.72 for acidity. These results suggest that increasing fruit number per tree from one to five was more effective in improving the accuracy for EP and FSM than for FW, SSC, and acidity, but the effect of increasing fruit number per tree from five to ten was small for all traits.

Effects of the number of fruits per tree, tree replication, and yearly repetition on the broad-sense heritability (hB2) of fruit weight, ease of peeling, firmness of segment membrane, soluble solids content, and acidity (see Eq. [3] for hB2 definition). × = 1 tree and 1 fruit; ● = 1 tree and 5 fruits; □ = 1 tree and 10 fruits; ▲ = 2 trees and 5 fruits per tree.
Citation: HortScience 51, 6; 10.21273/HORTSCI.51.6.680

Effects of the number of fruits per tree, tree replication, and yearly repetition on the broad-sense heritability (hB2) of fruit weight, ease of peeling, firmness of segment membrane, soluble solids content, and acidity (see Eq. [3] for hB2 definition). × = 1 tree and 1 fruit; ● = 1 tree and 5 fruits; □ = 1 tree and 10 fruits; ▲ = 2 trees and 5 fruits per tree.
Citation: HortScience 51, 6; 10.21273/HORTSCI.51.6.680
Effects of the number of fruits per tree, tree replication, and yearly repetition on the broad-sense heritability (hB2) of fruit weight, ease of peeling, firmness of segment membrane, soluble solids content, and acidity (see Eq. [3] for hB2 definition). × = 1 tree and 1 fruit; ● = 1 tree and 5 fruits; □ = 1 tree and 10 fruits; ▲ = 2 trees and 5 fruits per tree.
Citation: HortScience 51, 6; 10.21273/HORTSCI.51.6.680
Increasing the number of trees per genotype to two with five fruits per tree and single-year measurement increased hB2 from 0.32 to 0.34 for FW, 0.69 to 0.74 for EP, 0.23 to 0.27 for FSM, 0.38 to 0.41 for SSC, and 0.71 to 0.76 for acidity. Increasing yearly repetition to two with five fruits per tree on a single tree increased hB2 from 0.32 to 0.48 for FW, 0.69 to 0.80 for EP, 0.23 to 0.34 for FSM, 0.55 to 0.81 for SSC, and 0.71 to 0.81 for acidity. In tree replication, hB2 increased because of a decrease in environmental variance associated with σt2 and σty2, which were smaller than σf2 and σgy2 for each fruit trait except for FW, indicating that tree replication would increase hB2 only slightly. In yearly repetition, hB2 increased because of a decrease in environmental variance associated with σy2, which was the largest environmental variance component for FW and SSC, in addition to decreases in σgy2, σty2, and σf2, indicating that yearly repetition would increase hB2, especially for FW and SSC. Therefore, for the five fruit traits investigated in this study, yearly repetition increased hB2 more efficiently than tree replication did. Increasing the fruit number per tree decreased the environmental variance associated with σf2, whereas yearly repetition decreased the environmental variance associated with not only σf2, but also σgy2, which was the second largest environmental variance component for EP (13%) and for FSM (22%). The efficiency of five fruits per tree per year with 2-year repetition in decreasing σE2 was higher than that of 10 fruits per tree per year with single-year measurements, because the former decreased the environmental variance associated with σy2, σgy2,and σty2 in addition to σf2 (i.e., σE2 = σt2 + σy2/2 + σgy2/2 + σty2/2 + σf2/10), whereas the latter decreased only variance associated with σf2 (i.e., σE2 = σt2 + σy2 + σgy2 + σty2 + σf2/10). Therefore, increasing both the number of fruits and yearly repetition reduced σE2 and increased hB2 more efficiently than increasing only the number of fruits for EP and FSM, which have large σf2 and σgy2.
These results are in line with those of the studies of FW and SSC of Japanese pear (Machida and Kozaki, 1975), Japanese persimmon (Yamada et al., 1993), and grape (Sato et al., 2000). Increasing yearly repetition is more reasonable than increasing the number of tree replications if both increase hB2 similarly, because the latter requires more space and effort (Nishio et al., 2014). Variance by year could be decreased not only by yearly repetition as above, but also by a yearly adjustment procedure using the mean performance of the control genotype over several years (Yamada et al., 1994). However, the efficiency of this adjustment and the optimum number of control genotypes, which depend on the difference between the variance before and after the adjustment, have not been studied for fruit traits in citrus and therefore should be examined in further studies.
The overall citrus genetic resources including pummelo or grapefruit, which are partially used for crossbreeding, likely have larger genetic variance than population of mandarin and related cultivars because the citrus gene pool as a whole has a wider range of fruit traits. In addition, the population used in the study were grown using normal commercial practices, at much lower densities than the populations of cross-bred seedlings in the NIFTS breeding program. Therefore, whether our results could be applied to populations of cross-bred seedlings should be examined in future studies.
Our estimates of environmental variance components are useful not only for evaluation of the relative efficiency of yearly repetition, tree replication, and fruit sample size, but also for genomics-based breeding approaches, including quantitative trait locus (QTL) analyses, genome-wide association studies, and genomic selection, which need accurate phenotype evaluation to build an accurate phenotype database. The effectiveness of molecular markers associated with detected QTLs should be determined as the percentage of explained genetic variance instead of phenotypic variance, because the fluctuation of phenotypic values as a result of environmental variation confounds marker effects (Nishio et al., 2011). Although some molecular markers associated with QTLs for fruit traits of citrus have been reported (Gulsen et al., 2011; Sahin-Cevik and Moore, 2012; Sugiyama et al., 2011), these reports show the contribution of these markers only to phenotypic values, not to genetic variance. Therefore, the estimation of environmental variances for commercially important fruit traits such as FW, EP, FSM, SSC, and acidity could facilitate not only citrus phenotypic selection but also genomic-based breeding.
Literature Cited
Gulsen, O., Uzun, A., Seday, U. & Kafa, G. 2011 QTL analysis and regression model for estimating fruit setting in young citrus trees based on molecular markers Sci. Hort. 130 418 424
Machida, Y. & Kozaki, I. 1975 Quantitative studies on the fruit quality for Japanese pear (Pyrus serotina Rehder) breeding. I. Statistical analyses of cultivar populations J. Jpn. Soc. Hort. Sci. 44 235 240
Machida, Y. & Kozaki, I. 1976 Quantitative studies on the fruit quality for Japanese pear (Pyrus serotina Rehder) breeding. II. Statistical analyses of a hybrid seedling population J. Jpn. Soc. Hort. Sci. 44 325 329
Nishio, S., Yamada, M., Sawamura, Y., Takada, N. & Saito, T. 2011 Environmental variance components of fruit ripening date as used in both phenotypic and marker-assisted selection in Japanese pear breeding HortScience 46 1540 1544
Nishio, S., Yamada, M., Takada, N., Kato, H., Onoue, N., Sawamura, Y. & Saito, T. 2014 Environmental variance and broad-sense heritability of nut traits in Japanese chestnut breeding HortScience 49 696 700
Nishiura, M. 1964 Citrus breeding and bud selection in Japan Proc. Fla. State Hort. Sci. 77 79 83
Nonaka, K., Kita, M., Ikoma, Y., Fukamachi, H., Imai, A., Yoshioka, T. & Yamada, M. 2012 Genetic differences and environmental variations in carotenoid concentrations of fruit flesh in parental population used in citrus breeding in Japan J. Amer. Soc. Hort. Sci. 137 243 249
Sahin-Cevik, M. & Moore, G.A. 2012 Quantitative trait loci analysis of morphological traits in citrus Plant Biotechnol. Rpt. 6 47 57
Sato, A., Yamada, M., Iwanami, H. & Hirakawa, N. 2000 Optimal spatial and temporal measurement repetition for reducing environmental variation of berry traits in grape breeding Sci. Hort. 85 75 83
Sugiyama, A., Omura, M., Matsumoto, H., Shimada, T., Fujii, H., Endo, T., Shimizu, T., Nesumi, H. & Ikoma, Y. 2011 Quantitative trait loci (QTL) analysis of carotenoids content in citrus fruit J. Jpn. Soc. Hort. Sci. 80 136 144
Tanaka, T. 1954 Species problem in citrus (revise auranti acearum IX). Jpn. Soc. Promotion Sci., Tokyo, Japan
Yamada, M., Sato, A. & Ukai, Y. 2002 Genetic differences and environmental variations in calyx-end fruit cracking among Japanese persimmon cultivars and selections HortScience 37 164 167
Yamada, M., Yamane, H. & Ukai, Y. 1994 Efficiency of use of control genotype for reducing yearly fluctuations of quantitative fruit characters in Japanese persimmon breeding Bull. Fruit Tree Res. Sta. 26 29 37
Yamada, M., Yamane, H., Yoshinaga, K. & Ukai, Y. 1993 Optimal temporal and spatial measurement repetition for selection in Japanese persimmon breeding HortScience 28 838 841