## Abstract

We evaluated the nut harvesting date (NHD), nut weight (NW), pericarp splitting (PS), and infestation by insects (II) in eight cultivars/selections of Japanese chestnut, including a Japanese–Chinese hybrid, over 6 years. Data were analyzed by analysis of variance (without transformation for NHD, after log-transformation for NW and PS, and after square root transformation for II). The among-tree variance accounted for only 1.1% to 8.5% of the total variance. The variance component resulting from residual factors for the tree × year interaction and sampling errors was the largest component for NW, PS, and II, accounting for 46% to 54% of the total environmental variance. Because tree replication is costly and time-consuming in chestnut breeding, increasing the number of yearly repetitions is more efficient than increasing the number of tree replicates. Broad-sense heritability was 0.84 for NHD, 0.27 for NW, 0.48 for PS, and 0.17 for II in evaluations with one tree without yearly repetition. It increased to 0.91 for NHD, 0.40 for NW, 0.62 for PS, and 0.29 for II in evaluations with one tree in 2 years. For NHD, the heritabilities are sufficiently high to distinguish genetic differences among cultivars/selection. In contrast, the low heritability of II suggests that this trait should not be evaluated with an emphasis on the initial selection stage but rather with an emphasis on the secondary selection stage based on testing at several locations with a large number of yearly and tree replications.

The chestnut (*Castanea* spp.) has been an important food source in East Asia, Asia Minor, Europe, and North America for centuries. Three main chestnut species are currently produced commercially and are naturally distributed around the world: Japanese chestnut (*C. crenata* Sieb. et Zucc.) in Japan and the Korean Peninsula, Chinese chestnut (*C. mollissima* Bl.) in China, and European chestnut (*C. sativa* Mill.) in Europe and Asia Minor.

The history of chestnut use in Japan is very long. For example, a large number of chestnut pericarps and large chestnut timbers were found in a tower at the Sannai-Maruyama ruins (3500 to 2000 BCE) in Aomori Prefecture, northern Japan (Sawamura, 2006). For primitive Japanese people, chestnuts were important sources of food and construction materials. Even now, they are commercially grown throughout Japan for food, and chestnut orchards covered more than 21,400 ha in 2011, ranking fourth among woody fruit and nut crops in Japan after citrus, apple, and persimmon (FAOSTAT, 2013).

An organized chestnut breeding program started in 1947 at the National Institute of Fruit Tree Science in Japan and continues at the National Agriculture and Food Research Organization (NARO) Institute of Fruit Tree Science (NIFTS; Kotobuki et al., 1999; Pereira-Lorenzo et al., 2012). Currently, the most important target of the program is to release Japanese chestnut cultivars with an easily peeled pellicle. One such elite cultivar, Porotan, was released by NIFTS in 2006 (Pereira-Lorenzo et al., 2012; Saito et al., 2009). In addition, it is desirable to broaden the range of NHD to provide cultivars that are adapted to a much wider range of growing seasons. Simultaneously, large nut weight, the absence of PS, and freedom from II by the peach moth, *Conogethes punctiferalis* (Lepidoptera: Crambidae), are also important targets. The easily peeled pellicle trait in Japanese chestnut is controlled by a major gene at a single locus (Takada et al., 2012). Molecular markers associated with this trait are now being used in applied chestnut breeding programs (Nishio et al., 2013). On the other hand, the inheritance of NHD and NW has been insufficiently studied. Kotobuki et al. (1984) reported the narrow-sense heritability of NHD and NW using F_{1} progeny from the chestnut breeding program at NIFTS, but their experiment was not replicated over multiple years, and they did not estimate the environmental variance components. In addition, no reports are available for the inheritance of PS and II and the associated selection efficiency. However, breeders have recorded the performance data for those traits of cultivars/selections for many years in the NIFTS breeding program.

Environmental variance in an orchard comprises the variance among years and variance among trees within a genotype as well as the genotype × year interaction variance, the tree × year interaction variance, and the variance among nuts (fruits) within a tree, which can all be estimated by repeating measurements of replications of several genotypes in 2 or more years (Hansche and Beres, 1966; Hansche and Brooks, 1965; Sato et al., 2000; Yamada et al., 1993, 2002). Genetic and environmental variance has been estimated for many fruit crops such as sweet cherry (Hansche and Beres, 1966), Japanese persimmon (Yamada et al., 1993), strawberry (Sacks and Shaw, 1994; Shaw, 1991), peach (De Souza et al., 1998), grape (Sato et al., 2000), guava (Thaipong and Boonprakob, 2005), and Japanese pear (Nishio et al., 2011). Once accurate estimates of the components of environmental variance have been obtained under a given set of cultural conditions, the environmental variance under any given number of replications and measurement repetitions can be estimated for those conditions (Nyquist and Baker, 1991; Yamada et al., 1993, 2002).

It is necessary to estimate the broad-sense heritability to determine whether the breeder has made or will be able to make an effective selection. If the broad-sense heritability is zero, any selection is meaningless, and the breeder cannot select superior genotypes. Therefore, breeders should concentrate on traits with high heritability. Broad-sense heritability is defined as the ratio of genetic variance to phenotypic variance. However, the broad-sense heritability of a breeding offspring population is difficult to estimate using the seedlings from controlled crosses, because the use of seedlings does not permit replication of each seedling’s genotype and because the seedlings quickly outgrow the space allowed in the breeding nursery. Instead, we can estimate the broad-sense heritability using environmental variance obtained from cultivars/selections over 2 or more years.

Environmental variance and broad-sense heritability are important both for selection and for future genomics-based breeding. All genomics-based breeding techniques, including quantitative trait locus (QTL) analyses, genome-wide association studies (GWAS), and genomic selection (GS), need accurate phenotype evaluations to build an accurate phenotype database. Phenotypic variance contains both genetic and environmental variance components. The fluctuation of phenotypic values in response to changes in environmental conditions hinders the identification of useful molecular markers and can provide misleading results for marker effects detected by these analyses.

The objective of the present study was to accurately estimate the components of environmental variance and the broad-sense heritability of NHD, NW, PS, and II with the goal of supporting phenotypic selection in chestnut breeding programs, studies to detect linkages between molecular markers and traits of interest, and future genomics-based breeding.

## Materials and Methods

#### Expt. 1: Estimation of variance in nut weight among fruits within a tree (100 nuts × two trees × four cultivars, 1 year).

This experiment was based on two trees each of four cultivars (‘Ishizuchi’, ‘Mikuri’, ‘Porotan’, and ‘Tsukuba’; Table 1). The trees were 12 years old in 2011, were grafted on the seedling roots of ‘Kunimi’, and were randomly planted at a spacing of 6.0 m × 7.2 m in an orchard at NIFTS in Tsukuba, Ibaraki, Japan (lat. 36°04′N, long. 140°10′E). The trees were pruned in winter to promote stable nut production and were grown using the standard cultural techniques used in commercial production orchards in Japan.

The eight Japanese chestnut and Japanese–Chinese hybrid cultivars/selections used in Expts. 1 and 2.

Up to three nuts are contained in a bur that develops on a shoot of Japanese chestnut. In 2011, all the burs from each tree of the four cultivars were harvested when the bur had changed from green to brown and dropped. After removal of the bur, NW of each nut was measured as the fresh weight. Log-transformed values of NW were used in the following statistical analysis.

*P*

_{1ijk}) for the

*k*th nut of the

*j*th tree of the

*i*th genotype (cultivar/selection) was expressed using the following model:

μ_{1} = a constant value (the overall mean),

*g*_{1i} = the random effect of the *i*th genotype,

*t*_{1ij} = the random effect of the *j*th tree of the *i*th genotype, and

*n*_{1ijk} = the random effect of the *k*th nut of the *j*th tree of the *i*th genotype.

Cochran’s test (Cochran, 1941) did not reject homogeneity of variance within a tree at *P* < 0.05. ANOVA provided the variance among genotypes (σ_{g1}^{2}), among trees within a genotype (σ_{t1}^{2}), and among the nuts within a tree (σ_{n1}^{2}). The distribution of the residual estimates was not significant at *P* = 0.05 by the one-sample Kolmogorov–Smirnov test (Campbell, 1974), and it approached a normal distribution. The ANOVA was performed in JMP Version 9.0.2 software (SAS Institute, 2012).

#### Expt. 2: Estimation of environmental variance components in nut harvesting date, nut weight, pericarp splitting, and infestation by insects (three trees × eight cultivars, 6 years).

This experiment was based on three trees each of eight cultivars/selections (‘Ishizuchi’, ‘Kunimi’, ‘Mikuri’, ‘Porotan’, ‘Riheiguri’, ‘Tanzawa’, ‘Tsukuba’, and ‘Tsukuba-37’; Table 1). The trees were 12 years old in 2011, were grafted on the seedling roots of ‘Kunimi’, and were randomly planted at a spacing of 6.0 m × 7.2 m in the orchard at NIFTS. We evaluated NHD, PS, II, and NW over 6 years, from 2004 to 2009 (Table 2).

The performance of the cultivars/selections in terms of nut traits.^{z}

The burs and nuts were harvested as in Expt. 1. After removal of the bur, the number of nuts harvested on a given day was recorded. The number of harvested nuts per tree was at least 64 and the average was ≈300. The NHD of each nut was expressed as the number of days after 1 Aug. Fifty nuts were randomly chosen, and the average value of NHD was analyzed. PS was evaluated as the proportion of nuts with a split pericarp and II was evaluated as the proportion of the nuts infested by insects out of 50 randomly chosen nuts from each tree. Because nuts were harvested every 3 or 4 d, the number of intact nuts on a given harvesting day was recorded and the group weight was measured (in grams) on a digital scale on each date. We chose this approach because it is impractical to weigh each nut individually in practical breeding operations. The average individual NW was calculated as the total intact nut weight divided by the total number of intact nuts.

*j*th tree of the

*i*th genotype in the

*k*th year was expressed as follows:

μ_{2} = a constant value (the overall mean),

*g*_{2i} = the random effect of the *i*th genotype,

*t*_{2ij} = the random effect of the *j*th tree of the *i*th genotype,

*y*_{2k} = the random effect of the *k*th year,

(*gy*_{2})_{ik} = the interaction between the *i*th genotype and the *k*th year, and

*r*_{2ijk} = the residual effect in the *j*th tree of the *i*th genotype in the *k*th year.

ANOVA provided the variance among genotypes (σ_{g2}^{2}), among trees within a genotype (σ_{t2}^{2}) and among years (σ_{y2}^{2}) as well as the variance associated with the genotype × year interaction (σ_{gy2}^{2}) and the variance associated with the residuals (σ_{r2}^{2}). The residual variance included the tree × year interaction and sampling errors other than those for the year effect, genotype × year interaction, and tree effect.

Before statistical analysis, the NW and PS values were log-transformed and the II values were square root-transformed to improve the normality of the distribution of the residual estimates. Because PS values of 0 could not be log-transformed, they were set to 0.01 (the minimum measured value for this trait). The distribution of the residuals was not significantly different from a normal distribution at *P* = 0.05 by the Kolmogorov–Smirnov one-sample test (Campbell, 1974) for all traits. Therefore, the ANOVA model was applicable to the data.

The environmental variance for a genotype (σ_{e}^{2}) was estimated as σ_{y2}^{2}/*y* + σ_{gy2}^{2}/*y* + σ_{t2}^{2}/*t* + σ_{r2}^{2}/*ty*, where *y* is the number of years and *t* is the number of tree replicates per genotype. We defined the broad-sense heritability (*h*_{B}^{2}) as σ_{g2}^{2}/(σ_{g2}^{2} + σ_{e}^{2}).

## Results and Discussion

#### Expt. 1: Estimation of variance in nut weight among fruits within a tree (100 nuts × two trees × four cultivars, 1 year).

To determine the correct sample size for future determination of mean nut weight, we estimated the variance in NW among nuts within a tree. One-way ANOVA showed that the genotypic and tree effects were both highly significant (*P* < 0.01; Table 3). Variance components were σ_{g1}^{2} = 39.29 × 10^{−4}, σ_{t1}^{2} = 2.42 × 10^{−4}, and σ_{n1}^{2} = 98.59 × 10^{−4}. Although σ_{g1}^{2} and σ_{t1}^{2} had a small number of df (df = 3 and 4, respectively) and a high relative se (85.2% and 98.8%, respectively), and thus had low reliability, σ_{n1}^{2} had a large df (792) and a low relative se (5.0%) and therefore provides a reliable estimate of the variance among nuts within a tree in the model we used. σ_{n1}^{2} accounted for 70% of the total variance. The average NW of *a* nuts of *b* trees in a genotype has an environmental variance of σ_{n1}^{2}/*ab* + σ_{t1}^{2}/*b*. Increasing the number of nuts evaluated effectively reduces the environmental variance component among the nuts within a tree. The environmental variance component among the nuts for the average NW of 50 and 100 nuts from a tree can be estimated as σ_{n1}^{2}/50 (= 1.97 × 10^{−4}) and σ_{n1}^{2}/100 (= 0.99 × 10^{−4}). These estimates were considerably smaller than σ_{g1}^{2}, which suggests that 50 to 100 nuts are sufficient to determine accurately the mean and variance of the nuts.

Analysis of variance for the nut weight records for 100 nuts from two-tree replicates of four chestnut cultivars (genotypes) in 2011 (Expt. 1).

#### Expt. 2: Estimation of environmental variance components using three trees each of eight cultivars/selections with trait evaluations repeated for 6 years.

Two-way ANOVA showed that the genotype and genotype × year interaction were significant for all traits (Table 4). The percentage of total variance accounted for by the genetic variance component (σ_{g2}^{2}) was larger for NHD (84.3%) than for NW (27.2%), PS (47.9%), and II (17.1%; Table 5). A large chestnut with no PS and no II is the preferred type in the Japanese market, and the eight cultivars/selections used in this experiment had a larger nut size and a smaller frequency of PS and II than the average of those in the Japanese chestnut germplasm collection at Tsukuba. Therefore, they should have a smaller genetic variance in NW, PS, and II than in the overall genetic resources for Japanese chestnut.

Analysis of variance for four nut traits using eight cultivars/selections with three trees per genotype for six years (Expt. 2).^{z}

Estimates of the variance components in the analysis of variance for nut traits using eight cultivars/selections with three trees per genotype for 6 years.

The year effect was significant for NHD, NW, and II (Table 4). The percentage of the total variance (σ_{T2}^{2}) accounted for by the variance component resulting from the year effect (σ_{y2}^{2}) was 10.3% for NHD, 17.5% for NW, and 20.8% for II (Table 5). The tree effect was significant for NHD, NW, and PS (Table 4). However, the percentage of the total variance (σ_{T2}^{2}) accounted for by the tree effect (σ_{t2}^{2}) was small and ranged from 1.2% to 8.5% (Table 5). For NW, PS, and II, the largest environmental variance component was the result of the residuals (the tree × year interaction and the sampling error within the tree), which accounted for 46% to 54% of the total environmental variance.

For NW, the number of nuts evaluated per tree was not constant. Strictly speaking, this means that the data were not highly appropriate for the ANOVA model, but weighing each nut requires more time and labor than would be available in practical breeding, which requires efficient measurements. The σ_{r2}^{2} in NW was 17.8 × 10^{−4} and consisted of the effects of the tree × year interaction and within-tree variation among the nuts. In this study, the number of nuts evaluated per tree was at least 64 and was mostly more than 100; therefore, the variance resulting from within-tree variation among the nuts was estimated as less than 1.54 × 10^{−4} from the results in Expt. 1, indicating that most of σ_{r}^{2} was accounted for by the tree × year interaction variance and that the variance among nuts within a tree was minimal. Therefore, the bias for the ANOVA model was minimal, and the result can be used for practical chestnut breeding.

#### Estimation of broad-sense heritability.

Broad-sense heritability was estimated at 0.84 for NHD, 0.27 for NW, 0.48 for PS, and 0.17 for II when the evaluation was based on one tree without yearly repetition (Table 6). It increased to 0.91 for NHD, 0.40 for NW, 0.62 for PS, and 0.29 for II when the evaluation was based on one tree with 2 years of repetition. For NHD, these heritabilities are high and are sufficient to distinguish genetic differences in this trait. However, the heritabilities of NW and II were low even when based on two trees and 2 years of evaluation (0.50 and 0.36), probably owing to small genetic differences in the study population. For NW in a previous study, genetic variance of the seedlings seemed to be larger than that of the cultivars and selections (Kotobuki et al., 1984) because the cultivars and selections had been selected for large nut size. Thus, we need to estimate broad-sense heritability in the seedlings from controlled crosses in further studies. The low heritability of II suggests that crossing parent genotypes with a high resistance to II chosen from the NIFTS genetic resources would improve breeding for insect-resistant cultivars. II should not be evaluated with an emphasis on the initial selection stage, but rather with an emphasis on a secondary selection stage based on testing at several locations with a large number of yearly and tree replicates.

Broad-sense heritability (*h*_{B}^{2}) of nut traits using eight cultivars/selections with three trees per genotype for 6 years.

When the evaluation was based on two trees in 1 year, broad-sense heritability also increased, to 0.86 for NHD, 0.35 for NW, 0.57 for PS, and 0.22 for II. Increasing the number of yearly repetitions was as efficient as increasing the number of tree replicates for evaluating the genetic potential in breeding. Unlike annual crops, fruit and nut trees have a long juvenile period before flowering and fruit or nut set. Because tree replication requires a large space and an enormous effort, and given the time and cost over the long juvenile period, increasing the number of yearly repetitions is more reasonable than increasing the number of tree replicates.

#### Application to genomics-based breeding.

The estimated environmental variance is useful in both conventional breeding and genomics-based breeding based on methods such as QTL analyses, GWAS, and GS. Accurate phenotype evaluation is necessary to develop molecular markers or to predict the value of a given cultivar in breeding. The genetic variance depends on the plant population: the seedling population for QTL analyses, the set of cultivars for GWAS, and all of the seedling population for GS. It is difficult to evaluate the traits of chestnut seedling trees over several years because the trees grow swiftly and become too large to manage easily in a breeding orchard. Because chestnut seedlings were traditionally planted at a spacing of 2 m × 5 m in the NIFTS orchard, the evaluation of the seedlings for genetic analyses was limited to 2 years. If the phenotypic variance of seedlings can be obtained, we can estimate heritabilities using the environmental variance estimated in the present study and predict how accurately we can evaluate the performance of the trees for each of the four traits that we studied. We are currently preparing F_{1} populations, a phenotype database for the cultivars, and records of the seedling populations used in the NIFTS chestnut breeding program to support these genetic analyses. Accurate evaluation of the phenotype based on sufficient repetition to reduce the environmental variance will improve future genomics-based breeding.

## Literature Cited

Campbell, R.C. 1974 Statistics for biologists. 2nd Ed. Cambridge Univ. Press, Cambridge, UK

Cochran, W.G. 1941 The distribution of the largest of a set of estimated variances as a fraction of their total

*Ann. Eugen.*11 47 52De Souza, V.A.B., Byrne, D.H. & Taylor, J.F. 1998 Heritability, genetic and phenotypic correlations, and predicted selection response of quantitative traits in peach: I. An analysis of several reproductive traits

*J. Amer. Soc. Hort. Sci.*123 598 603FAOSTAT 2014 <http://www.fao.org>

Hansche, P.E. & Beres, V. 1966 An analysis of environmental variability in sweet cherry

*Proc. Amer. Soc. Hort. Sci.*88 173 188Hansche, P.E. & Brooks, R.M. 1965 Temporal and spatial repeatability of a series of quantitative characters in sweet cherry (

*Prunus avium*L.)*Proc. Amer. Soc. Hort. Sci.*86 120 128Kotobuki, K., Machida, Y., Sato, Y., Kajiura, I. & Kozono, T. 1984 Chestnut breeding program in National Institute of Fruit Tree Science, Japan

*Bull. Fruit Tree Res. Stn. A*11 43 53 [in Japanese]Kotobuki, K., Saito, T., Kashimura, Y. & Shoda, M. 1999 Chestnut breeding program in National Institute of Fruit Tree Science, Japan

*Acta Hort.*494 323 325Nishio, S., Takada, N., Yamamoto, T., Terakami, S., Sawamura, Y. & Saito, T. 2013 Mapping and pedigree analysis of the gene that controls the easy peel pellicle trait in Japanese chestnut (

*Castanea crenata*Sieb. et Zucc.)*Tree Genet. Genomes*9 723 730Nishio, 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 1544Nyquist, W.E. & Baker, R.J. 1991 Estimation of heritability and prediction of selection response in plant populations

*Crit. Rev. Plant Sci.*10 235 322Pereira-Lorenzo, S., Ballester, A., Corredoira, E., Vieitez, A.M., Agnanostakis, S., Costa, R., Bounous, G., Botta, R., Beccaro, G.L., Kubisiak, T.L., Conedera, M., Krebs, P., Yamamoto, T., Sawamura, Y., Takada, N., Gomes-Laranjo, J. & Ramos-Cabrer, A.M. 2012 Chestnut. p. 729–769. In: Badenes, M.L. and D.H. Byrne (eds.). Fruit breeding. Springer, New York, NY

Sacks, E.J. & Shaw, D.V. 1994 Optimum allocation of objective color measurements for evaluating fresh strawberries

*J. Amer. Soc. Hort. Sci.*119 330 334Saito, T., Kotobuki, K., Sawamura, Y., Abe, K., Terai, O., Shoda, M., Takada, N., Sato, Y., Hirabayashi, T., Sato, A., Nishibata, T., Kashimura, Y., Kozono, T., Fukuda, H., Kihara, K., Suzuki, K. & Uchida, M. 2009 New Japanese chestnut cultivar ‘Porotan’.

*Bull. Natl. Inst. Fruit Tree Sci.*9 1 9 [in Japanese with English abstract]SAS Institute 2012 JMP. Version 9.0.2. SAS Institute Inc., Cary, NC

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 83Sawamura, Y. 2006 Chestnut, p. 82–85. In: Jpn. Soc. Hort. Sci. (ed.). Horticulture in Japan 2006. Shoukadoh, Kyoto, Japan

Shaw, D.V. 1991 Genetic-variation for objective and subjective measures of fresh fruit color in strawberries

*J. Amer. Soc. Hort. Sci.*116 894 898Takada, N., Nishio, S., Yamada, M., Sawamura, Y., Sato, A., Hirabayashi, T. & Saito, T. 2012 Inheritance of the easy-peeling pellicle trait of Japanese chestnut cultivar Porotan

*HortScience*47 845 847Thaipong, K. & Boonprakob, U. 2005 Genetic and environmental variance components in guava fruit qualities

*Sci. Hort.*104 37 47Yamada, 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 167Yamada, M., Yamane, H., Yoshinaga, K. & Ukai, Y. 1993 Optimal spatial and temporal measurement repetition for selection in Japanese persimmon breeding

*HortScience*28 838 841