Original Article

Genetic Variation in Growth Characteristics and Wood Properties of Ficus variegata Blume First Generation Progeny Trials in Indonesia

Liliek HARYJANTO1,2, Sapto INDRIOKO1,https://orcid.org/0000-0003-1828-469X, Arif NIRSATMANTO2, Fanny HIDAYATI1
Author Information & Copyright
1Faculty of Forestry, Universitas Gadjah Mada, Yogyakarta 55281, Indonesia
2Research Centre for Applied Botany, National Research and Innovation Agency Republic of Indonesia, Bogor 16911, Indonesia
Corresponding author: Sapto INDRIOKO (e-mail: sindrioko@ugm.ac.id)

Copyright 2024 The Korean Society of Wood Science & Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Feb 06, 2024; Revised: Apr 17, 2024; Accepted: Jun 17, 2024

Published Online: Sep 25, 2024

ABSTRACT

Two first-generation progeny trials of Ficus variegata Blume were planted in Yogyakarta, Indonesia, consisting of 17 families from the provenances of West Nusa Tenggara (WNT) and 19 families from the provenances of Cilacap-Pangandaran (C-P), respectively. The trials were evaluated after 10 years for growth characteristics [diameter (D), tree height (H) and stem volume (V)] and wood properties [stress-wave velocity (SWV) and Pilodyn penetration (P)]. Genetic variation, the coefficient of additive genetic variation (CVA), and heritability estimation were analyzed. Subsequently, genetic correlation between traits was estimated. The analysis of variance showed that there were significant differences in growth characteristics and wood properties in the WNT families, with significance observed across most factors except for height and P in the C-P families. The CVA in growth characteristics (D, H, V) was higher than for wood quality (SWV and P) in WNT and C-P families. Estimated family heritability (h2f) for growth characteristics, SWV, and P were high in the WNT families but moderate in the C-P families showing that genetic variation in the observed traits was more additive in the WNT families. The positive estimated genetic correlations between growth characteristics in two progeny trials, and the moderate to strong negative genetic correlation between D and P and also between P and SWV showed that growth characteristics and wood quality can be genetically improved simultaneously by using D as a selection criterion is an appropriate breeding strategy for F. variegata.

Keywords: Ficus variegata Blume; growth characteristics; heritability; genetic correlation; Pilodyn penetration

1. INTRODUCTION

The decline in the potential of Indonesia’s natural forests to fulfill timber needs has led to the development of a plantation resource to help meet this demand. Promoting the use of native species is essential to meet the timber demand, as these species are expected to be well-suited and adapted to local biotic and abiotic conditions (Thomas et al., 2014). It is also necessary to reduce the dependence of the wood industry, especially on members of the Dipterocarpaceae (Augustina et al., 2020). Subsequently, nyawai (Ficus variegata Blume belongs to the Moraceae family) is a native species that has prospects for supplying plywood, face veneer (Hendromono and Komsatun, 2008), packaging boxes, and materials for light construction (Sumarni et al., 2009). It is a fast-growing tree species that has attracted the interest of private companies since 2003 (Effendi and Mindawati, 2015) and is being developed as an ‘alternative’ plant for plantations (Susanti and Halwany, 2017). In Indonesia, the natural distribution of this species includes Kalimantan (Haryjanto and Hadiyan, 2014; Wahyuningtyas et al., 2022), Java (Nurtjahjaningsih et al., 2019), Sumatra (Nur’aini et al., 2013; Rofifah et al., 2021), Sulawesi (Pitopang, 2012), West Lombok, Papua (van Heist et al., 2010).

Implementing a selective breeding program is an effective way to increase plantation productivity and improve the genetic quality of trees (Nelson, 2023). Achieving effective selective breeding requires a thorough understanding of the heritability of growth and wood quality traits. Heritability is a statistical measure of the proportion of phenotypic variance explained by genetic variance (Waszak and Cieślik, 2016). Traits with high heritability can be targeted for improvement through selective breeding because they are more likely to be passed on to future generations. The correlations between traits is also important for plant breeding programs in order to create an understanding how one trait influences another through shared genetic factors (de Oliveira et al., 2015). However, in-depth knowledge in these areas is currently lacking for F. variegata.

Measures of the heritability of height and diameter (D) of F. variegata in progeny tests established in 2012 showed that the genetic control of these traits was weak in the early growth phase (Haryjanto et al., 2014). Accurate measurements of tree growth and wood quality at the time of harvesting are crucial from breeding and selection perspective. These measurements hold a high value for predicting the breeding values of tested genotypes. Wood density is regarded as the best single predictor of wood quality (Fundova et al., 2018). Many breeding programs prioritize it as an initial trait for improvement due to its correlations with stiffness, strength, and shrinkage (Hong et al., 2014). The specific gravity of F. variegata is 0.27, ranging from 0.20 to 0.43 (Sumarni et al., 2009), indicating that the wood lacks strength and improving wood quality is necessary to make the species more attractive as a timber species. Wood quality is the main concern when using fast-growing plantations, especially when producing solid wood (Ghani and Lee, 2021). Therefore, improving these properties will expand the utilisation of wood more widely (Laksono et al., 2023).

In breeding programs, the preferred approach for rapidly assessing wood properties in standing trees involves the use of non-destructive methods (Schimleck et al., 2019). The Pilodyn, originally developed to test for wood decay, is often used to assess wood density (Fukatsu et al., 2011). Meanwhile, stress-wave velocity (SWV) which measures the speed of sound waves in standing trees, has been used to measure the speed of sound waves in standing trees and has received considerable attention (Wang et al., 2001). Many studies have used SWV on a standing tree as a stiffness selection criterion (Fukatsu et al., 2015; Matheson et al., 2008) and the potential for prediction of modulus of elasticity (Ishiguri et al., 2008; Van Duong et al., 2022; Wang et al., 2001).

This study has two objectives: First, to investigate the inheritance of traits related to growth characteristics and wood quality. Second, to understand how these traits relate to each other, both genotypic and phenotypic correlations. The acquired information will be examined in the context of formulating suitable selection strategies for F. variegata breeding programs.

2. MATERIALS and METHODS

2.1. Study site

The study site was in Bantul, Yogyakarta Province, Indonesia (latitude 07°57′30″–07°57′54″S, longitude 110°26′07″–110°26′29″ E, altitude of 75–145 m above sea level). The slope of the site varied between 5°–30°, and the soil was identified as an Oxisol. The climate according to the Köppen classification is Af (tropical rainforest climate) to Am (tropical monsoon climate), or type C (medium wet) of the climate classification by Schmidt and Ferguson (1951), with a mean annual rainfall of 1,502 mm (Bappeda Kabupaten Bantul, 2011).

2.2. Materials

In this study, two progeny trials involving the first generation from 17 and 19 open-pollinated F. variegata trees were conducted in a half-sib trial established in December 2012. The mother trees were selected from two regions, and the trial was divided into two sublines situated 300 m apart. The experimental design of the progeny test was a randomized complete block design. Subsequently, the provenances of West Nusa Tenggara (WNT) represented 17 families, 5 non-contiguous tree plots in six replicates and Cilacap-Pangandaran (C-P) represented 19 families, 4 non-contiguous tree plots in six replicates. The planting space was 5 × 5 m.

2.3. Measurement

D was measured at 1.3 m above ground, tree height (H) was measured from the ground to the tip of the tree, and stem volume (V) was measured using the equation of Qirom and Supriyadi (2013) for 10-year-old trees:

ln V = a + b ln  ( D )  + c ln  ( H )
(1)

Where:

V: stem volume (m3), D: diameter at 1.3 m above ground (cm), H: tree height (m), a: –9.22846, b: 1.7456, c: 0.9759.

Stem SWV was measured using a commercially available handheld stress-wave timer (FAKOPP Microsecond Timer, FAKOPP Enterprise) according to previously described methods (Hidayati et al., 2013b; Ishiguri et al., 2007). The start and stop sensors were positioned at 150 and 50 cm from ground level, respectively. Using a small hammer, the start sensor was struck to generate a stress wave. When the stress-wave was received at the stop sensor, the travel time between the two sensors was recorded. Ten measurements of stress-wave propagation time were taken for each tree, and the mean value was calculated. The formula for the calculation of the stress-wave propagation time per unit length was described by Wang etal. (2000, 2004).

Wood density measured indirectly by Pilodyn penetration (P) was carried out using Pilodyn tester (strength of spring, 6 J; D of pin, 2.5 mm, Proceq) at 1.3 m above ground level at three positions for each tree (Hidayati et al., 2013b, 2019; Ishiguri et al., 2008; Wu et al., 2010) and the mean value was calculated. The bark at the measurement points was removed prior to measurement.

2.4. Statistical analyses

Analyses of variance (ANOVA) were carried out to evaluate the differences in D, H, V, P, and SWV among the tested families based on an individual-tree following the linear model:

Y i j k = μ + B i + F j + ε i j k
(2)

Where Yijk is the individual tree observation, μ: the overall mean, Bi: fixed effect of the i-th block, Fj: random effect of the j-th family associated with the average genetic effects of the open-pollinated families, εijk: the within-plot error.

Family heritability (h2f) was calculated using the following formula of Zobel and Talbert (1984):

h f 2 = σ 2 f σ w 2 / T R + σ f 2
(3)

Where σ2f: the family component of variance, σ2w: the within-plot error component of variance, TR: the harmonic mean number of trees per plot by replication.

The coefficient of additive genetic variation (CVA) is a mean-standardized index of the evolvability of a trait (Hill, 2010; Houle, 1992). It was described by Haapanen (2001):

CV A = ( σ A 2 / χ ) × 100
(4)

Where CVA: coefficient of additive genetic variation, σ2A: additive genetic variance. The additive genetic variance on a single-site basis was estimated by multiplying the among-family component of variance by four (σ2A = 4σ2f; White et al., 2007; Zobel and Talbert, 1984); χ: population mean for the trait.

Genetic correlation between pairs of traits (rGxy) was calculated following Zobel and Talbert (1984):

r G x y = σ f ( x y ) σ f ( x ) 2 σ f ( y ) 2
(5)

Where σf(xy) is the additive genetic covariance component between traits, and σ2f(x) and σ2f(y) are the additive variance components for traits x and y respectively.

The phenotypic correlation between pairs of traits (rP) was calculated as given by Falconer (1989):

r P = cov P σ P x σ P y
(6)

Where covP is the covariance of the two traits, and σPx and σPy are the SD for traits x and y respectively.

3. RESULTS and DISCUSSION

3.1. Results
3.1.1. Mean values and variation among characteristics

The mean, range of values, SEM for D, H, V, SWV, and P of the 17 WNT families and the 19 C-P families are shown in Table 1. The mean growth traits (D, H, V) of WNT families were higher than those of C-P families. Subsequently, P was smaller in WNT than C-P, 29.62 mm, and 33.57 mm, respectively. SWV of WNTs was smaller than that of C-P, 2.03 km/s and 2.08 km/s, respectively. Significant variations (1% level) were observed for all measured characteristics among WNT families. Additionally, differences in D, V, and SWV within the C-P families were found to be significant at the 5% level, as shown in Table 2.

Table 1. Mean value of growth characteristics and wood properties in WNT and C-P families
Trait WNT families C-P families
Mean Range (min-max) SEM Mean Range (min-max) SEM
D (cm) 17.41 6.05–30.89 0.23 16.30 7.17–27.07 0.21
H (m) 9.74 5.34–14.13 0.09 9.07 5.04–16.10 0.10
V (m3) 0.141 0.006–0.340 0.004 0.121 0.009–0.407 0.004
SWV (km/s) 2.03 1.65–2.48 0.01 2.08 1.73–2.63 0.01
P (mm) 29.62 19.67–38.33 0.17 33.57 26.50–40.00 0.14

WNT: West Nusa Tenggara, C-P: Cilacap-Pangandaran, D: diameter, H: tree height, V: stem volume, SWV: stress-wave velocity, P: Pilodyn penetration.

Download Excel Table
Table 2. Differences in growth characteristics, stress-wave velocity (SWV), Pilodyn penetration (P) of the WNT and C-P families
Parameter Trait
Diameter Height Volume SWV Pilodyn
WNT
 Mean square
  Family 85.034** 17.171** 0.029** 0.071** 27.844**
  Error 15.900 2.351 0.004 0.019 8.569
C-P
 Mean square
  Family 27.538* 4.441ns 0.009* 0.038* 12.065ns
  Error 14.844 3.303 0.004 0.023 8.187

* Significant at p < 0.05;

** Significant at p < 0.01;

ns Non-significant at p > 0.05.

WNT: West Nusa Tenggara, C-P: Cilacap-Pangandaran.

Download Excel Table
3.1.2. Genetic variation and heritability

The estimated CVA, variance components, and family heritability estimates are shown in Table 3. The estimated CVA depending on the analyzed traits, and growth traits was higher than wood quality. In the WNT, CVA ranged from 16.47% to 46.61% for the growth traits (D, H, and V) and it was 4.68% (SWV) and 6.20% (P); in C-P, CVA ranged from 5.16% to 24.02% for the growth traits (D, H, and V) and it was 2.58% (SWV) and 2.54% (P).

Table 3. Estimated coefficient of additive genetic variation (CVA), estimated the family component of variance (σ2f), estimated the within-plot error component of variance (σ2w) and family heritability (h2f) for diameter, height, volume, SWV and Pilodyn penetration (P) of the WNT and C-P families
Trait CVA (%) σ2f σ2w h2f
WNT
 Diameter 19.84 2.984 (15.80%)a 15.900 (84.20%)b 0.81
 Height 16.47 0.643 (21.47%) 2.351 (78.53%) 0.86
 Volume 46.61 0.001 (19.61%) 0.004 (80.39%) 0.85
 SWV 4.68 0.002 (10.43%) 0.019 (89.57%) 0.73
 Pilodyn 6.20 0.843 ( 8.95%) 8.569 (91.04%) 0.69
C-P
 Diameter 9.58 0.610 (3.95%) 14.843 (96.05%) 0.46
 Height 5.16 0.055 (1.63%) 3.302 (98.37%) 0.26
 Volume 24.02 0.0002 (0.46%) 0.044 (95.58%) 0.49
 SWV 2.58 0.0007 (3.04%) 0.023 (96.96%) 0.39
 Pilodyn 2.54 0.182 (2.18%) 8.186 (97.82%) 0.32

a,b The proportion of the family component variance (σ2f) and the within-plot error component of variance (σ2w) to the total phenotypic variance in parenthesis, respectively.

SWV: stress-wave velocity, WNT: West Nusa Tenggara, C-P: Cilacap-Pangandaran.

Download Excel Table

In the WNT, the estimated contribution of the family component of variance (σ2f) to the total phenotypic variance ranged from 15.80% to 21.47% for the growth traits (D, H, and V). For wood quality, it was 8.95% (P) and 10.43% (SWV). In C-P, the estimated contribution of the family component of variance (σ2f) to the total phenotypic variance ranged from 0.46% to 3.95% for the growth traits (D, H, and V). For wood quality, it was 2.18% (P) and 3.04% (SWV).

In the WNT, family heritability (h2f) for the growth traits (D, H, and V) ranged from 0.81 to 0.86, and for wood quality, it was 0.69 (P) and 0.73 (SWV); in C-P, family heritability for growth traits ranged from 0.26 (H) to 0.49 (V) and for wood quality, it was 0.32 (P) and 0.39 (SWV).

3.1.3. Phenotypic and genetic correlations between characteristics

In both WNT and C-P families, highly positive and significant phenotypic correlations were observed for the growth traits (D, H, and V). The correlations ranged between 0.811–0.969 for WNT families and 0.631–0.922 for C-P families, as shown in Table 4. Subsequently, no phenotypic correlations were observed between growth traits (D, H) and wood quality (SWV and P). Significant negative phenotypic correlations were observed between SWV and P in both families, –0.281 (WNT) and –0.125 (C-P), respectively.

Table 4. Genetic (upper diagonal) and phenotype (lower diagonal) correlation coefficients between traits in WNT and C-P families
Characteristics Diameter Height Volume SWV Pilodyn
WNT
 Diameter - 0.918 0.180 –0.194 –0.731
 Height 0.811** - 0.302 0.091 –0.640
 Volume 0.969** 0.863** - –0.367 –0.179
 SWV –0.011ns 0.063ns 0.040ns - –0.460
 Pilodyn –0.086ns –0.099ns –0.107* –0.281** -
C-P
 Diameter - 0.765 0.047 –0.051 –0.375
 Height 0.631** - 0.502 0.144 –0.664
 Volume 0.922** 0.819** - 0.706 –0.502
 SWV –0.176ns –0.080ns –0.150** - –0.786
 Pilodyn 0.016ns 0.076ns 0.050ns –0.125* -

* Significant at p < 0.05;

** Significant at p < 0.01;

ns Non-significant at p > 0.05.

WNT: West Nusa Tenggara, C-P: Cilacap-Pangandaran, SWV: stress-wave velocity.

Download Excel Table

Strong estimated positive genetic correlations were observed between D and H (rG = 0.918), and (rG = 0.765) in WNT and C-P families, respectively. Between other growth traits (D and V) and (H and V) they were low (0.047) and moderate (0.502), respectively. Estimated negative genetic correlations among growth characteristics (D, H, V) and P in the WNT families and C-P families were –0.179 to –0.731 and –0.375 to –0.664, respectively.

Estimated negative genetic correlations observed between D and P were strong (rG = –0.731) and moderate (rG = –0.375), and between P and SWV were moderate (rG = –0.460) and strong (rG = –0.786) in WNT families and C-P families, respectively.

3.2. Discussion
3.2.1. Mean values populations and variation among characteristics

In this study, the growth of F. variegata was comparatively lower than that observed in other sites. For instance, in Riam Kiwa (South Kalimantan), F. variegata planted there exhibited a mean annual increment (MAI) in D of 2.10–2.36 cm/yr (Fitriani, 2011), whereas in this study the MAI was 1.74 cm/yr and 1.63 cm/yr for the WNT and C-P populations, respectively. This difference may be attributed to the tropical climate in Riam Kiwa (South Kalimantan) characterized by high humidity and consistent rainfall throughout the year (Yasin et al., 2020). F. variegata, being a tropical tree, thrives in seasonal wet tropical rainforests, often found along watercourses (Spencer et al., 1996).

The SWV values of F. variegata (2.03 km/s and 2.08 km/s for WNT and C-P families, respectively) may be relatively lower than other tropical tree species. Subsequently, the SWV values were at 2.83 km/s for Gmelina arborea 5-year-old (Hidayati et al., 2017) and 3.08 km/s for 13-year-old Falcataria moluccana (Ishiguri et al., 2007). This indicates that G. arborea and F. moluccana are stiffer than F. variegata. This difference may be due to differences in the anatomical characteristics and basic density of these species.

P serves as an indicator of wood density; higher values correspond to lower wood density, and vice versa (Kha et al., 2012). The P value was higher than other tropical tree species. Specifically, the P value was at 25.9 mm for 4.5-year-old in white jabon (Neolamarckia cadamba; Chaerani et al., 2019), 18.7–27.6 mm for 12-year-old Tectona grandis (Hidayati et al., 2013b), 15.7–18.9 mm for 6-year-old Acacia mangium (Hidayati et al., 2019). This indicates that F. variegata has a lower wood density compared to these species.

In this study, a comparison of the growth and wood quality traits between the two populations (WNT and C-P) reveals the following: the WNT families exhibit superior growth and wood density compared to the C-P families. Wood density, as indicated by P, is lower for WNT families than for C-P families. Lower penetration shows denser and stronger wood and is negatively correlated with P (Wei and Borralho, 1997; Wu et al., 2010). Conversely, for wood stiffness, SWV in the WNT families is slightly smaller than the C-P families (2.03 and 2.08, respectively). Although higher values of SWV show greater stiffness (Ishiguri et al., 2007; Wang et al., 2001), it does not always correlate with smaller P as observed in T. grandis (Hidayati et al., 2013a; Seta et al., 2021), A. mangium (Hidayati et al., 2019). In summary, it can be concluded that growth characteristics and wood density are superior in WNT compared to the C-P family. However the wood stiffness are similar for WNT to CP.

With the exception of height (H) and P in the C-P families, all measured traits between the two populations (WNT and C-P) were statistically significantly different. These significant differences in characteristics indicate that these traits are influenced by genetic factors rather than the environment. Similar results were reported in other studies, where significant differences were observed in growth characteristics (D, H) and wood quality (SWV and P) among 65 families of A. mangium (Hidayati et al., 2019), as well as in D, H, and P among 105 families of white jabon (Chaerani et al., 2019). The significant differences observed among the families in these traits suggest the possibility of selecting superior families to support breeding programs.

3.2.2. Genetic variation and heritability

The CVA is recommended for reporting the results of studies on adaptive variation, as it provides a better guide than heritability to the long-term evolvability of a trait (Houle, 1992). This CVA serves as a valuable indicator for assessing the long-term evolution of a trait and its potential for adaptation under changing environmental conditions (Houle, 1992). This reflects the potential for a given trait to respond to selection pressures in a given population, the larger the coefficient, the greater the evolvability of the trait (Cheung, 2020). In this study, CVA in growth traits (D, H, V) was higher than for wood quality (SWV and P) in WNT and C-P families (Table 3). This suggests that growth traits could be hypothesized as being more related to species adaptation. Additionally, the CVA of V was higher than that of H or D in both populations. This result may reflect the fact that volume is a function of both H and D. Specifically, the CVA of D was higher than that of height in both populations, at 19.84% and 16.47% for WNT, and 9.58% and 5.16% for C-P, respectively. This finding contrasts with a previous study on Pinus elliottii, where the CVA for height (H) was higher at 5.06% compared to the CVA for D at 2.99% (Lai et al., 2017). Therefore, this study suggests that there is a greater scope for selection of D among families compared to height.

Family heritability (h2f) was higher in WNT families compared to C-P families for all traits (Table 3). Additionally, the proportion of the family component variance (σ2f) to the total phenotypic variance was consistently greater in WNT families for specific traits. For instance, in the D trait, the contribution of WNT families to the family component variance was 15.80%, whereas for C-P families, it was 3.95% (Table 3). The genetic variation in WNT families resulted in higher heritability values, indicating that the variation in these traits was primarily due to genetic differences among families (Lai et al., 2017; Ranjan and Gautam, 2020). This suggests a greater additive genetic effect passed on to offspring from WNT families compared to C-P families, emphasizing the significance of genetic variation in breeding programs (Mihai et al., 2014). The observed high heritability implies that selective breeding is more feasible for the considered traits.

Heritability and genetic correlations are important factors influencing the effectiveness of breeding programs, particularly in forest trees. A notable and crucial adverse genetic correlation is observed between tree productivity and wood density, commonly characterized by a strong negative relationship (Klápště et al., 2022). Meanwhile, heritability and genetic correlation are important genetic parameters for the effectiveness of indirect selection. When the genetic correlation between traits is high but heritability is low, indirect selection tends to be less effective in producing offspring with the desired parent traits. Low heritability means that the genetic influence on the traits is not strongly transmitted to the next generation, diminishing the impact of indirect selection. High family heritability was found in the WNT population for the D and P at 0.81 and 0.69, respectively (Table 3), and the strong genetic correlation of D and P rG = –0.731 (Table 4). This indicates that indirect selection of the D trait in the WNT population will improve wood quality.

In this study, family heritability (h2f) for the growth characteristics was higher than for wood quality in the WNT families, and except for H, also for the C-P families. In the WNT families, the family heritability of growth traits (D, H, V) ranged from 0.81 to 0.86, while the family heritability of wood quality traits (SWV, P) ranged from 0.69 to 0.73. Meanwhile, the family heritability in the C-P families ranged from 0.46 to 0.49 and ranged from 0.32 to 0.39 for growth (D, V) and wood quality (SWV, P) respectively, and other studies confirm this findings. The family heritability for height and P were 0.69 and 0.46, respectively in white jabon (Chaerani et al., 2019); in A. mangium, family heritability for height and P were 0.37 and 0.25 respectively (Hidayati et al., 2019). This indicates that growth traits are passed on to offspring more than wood quality traits.

The study showed that high heritability can also be obtained by using a limited number of families, for example, 17 families for WNT and 19 for C-P. The result was similar to in Quercus serrata with a family heritability of 0.80 with 17 families (Kang et al., 2007). However, it is important to note that high heritability does not necessarily mean that the trait is entirely determined by genetic or that it cannot be influenced by environmental factors. Subsequently, heritability estimates are specific to the population being studied and the range of environments and genetic variations present in that population. Traits with high heritability in one population may exhibit lower heritability in another population with different genetic and environmental factors.

3.2.3. Phenotype and genetic correlations between traits

Phenotypic and genetic correlations play important roles in tree improvement programs. Subsequently, phenotypic correlations help identify relationships between observable traits, while genetic correlations inform breeders about how they can prioritize the selection of certain traits, knowing that improvement in one trait may have a positive or negative effect on others. By considering both phenotypic and genetic correlations, breeders can formulate comprehensive strategies aimed at simultaneously enhancing multiple traits. This method allows for the development of forest plantations with desired outcomes, taking into account the interplay and relationships between different traits in the breeding process.

This study showed highly significant phenotypic correlations among growth traits (D, H, V) in both WNT families (0.811 to 0.969) and C-P families (0.631 to 0.922; Table 4), consistent with prior results observed in T. grandis (Hidayati et al., 2013a). Positive correlations between D and H were similarly reported in A. mangium (Hidayati et al., 2019), and Eucalyptus urophylla (Wei and Borralho, 1997), suggesting a strong interdependence between growth characteristics, with their values exhibiting predictable simultaneous changes.

There was no phenotype correlation between wood traits (SWV and P) and growth traits (D and height; Table 4). Similarly, no significant phenotype correlation was observed between growth characteristics (D and H) and SWV in teak (Hidayati et al., 2013b). This finding aligns with previous results on A. mangium (Hidayati et al., 2019), but in contrast to G. arborea, which exhibits a positive significant phenotype correlation between D and SWV (Hidayati et al., 2017). In this study, wood traits appear to be independent of growth traits, suggesting that families with superior characteristics for both growth and wood traits can be selected.

Positive genetic correlation between growth traits (D, H, V) ranged from 0.180 to 0.918 in the WNT families and ranged from 0.047 to 0.765 in the C-P families (Table 4). This indicates that these traits tend to be influenced by similar genetic factors. In other words, individuals with favorable genes for one of these traits are likely to also have favorable genes for the others. These results are consistent with other studies, such as those conducted on Populus ussuriensis Kom (Jin et al., 2019), P. elliottii (Lai et al., 2017) and Eucalyptus pellita (Fadwati et al., 2023). These traits are interconnected at the genetic level, and improving one trait can lead to improvements in the others.

P is closely related to the wood density at the outer part of the stem (Wessels et al., 2011; Wu et al., 2010). The estimated negative genetic correlations between P and D were strong (rG = –0.731) in the WNT families and moderate (rG = –0.375) in the C-P families (Table 4). Meanwhile, the genetic correlation of SWV with growth traits (D and H) of both ENT and C-P was both positive and negative, but the correlation was weak (rG = –0.194–0.144; Table 4). The strong and negative correlation observed between P and D implies that as D increases, there is a corresponding decrease in P. This association indicates that increasing D is linked to low P, which, in turn, is indicative of high wood density. This finding is consistent with previous study on white jabon (Chaerani et al., 2019), and E. urophylla (Kien et al., 2008). Some studies have shown an unfavourable relationship between the two traits found in Norway spruce (Nguyen et al., 2022), Catalpa bungei (Xiao et al., 2021), T. grandis (Hidayati et al., 2013b), and E. urophylla (Wei and Borralho, 1997), making it difficult to improve growth and wood quality at the same time. This study shows that the correlations is favourable for breeding both traits (growth and wood quality) simultaneously.

The genetic correlation between P and SWV was moderate (rG = –0.460) to strong (rG = –0.786) negative in the WNT and C-P families, respectively. In other words, selecting trees for increased P (which may be associated with wood density) could lead to a decrease in SWV and vice versa. Subsequently, higher SWV corresponded to lower P values, indicative of denser and stronger wood. This finding is consistent with previous studies on E. urophylla (Wei and Borralho, 1997), and Larix kaempferi (Fukatsu et al., 2015). The result indicates that P was a reliable indirect measure of wood density in F. variegata.

4. CONCLUSIONS

The genetic study on the growth and wood characteristics of F. variegata in Yogyakarta Province, Indonesia aimed to investigate the inheritance of traits related to growth characteristics and wood quality; and to understand how these traits relate to each other, both genotypic and phenotypic correlations. The results obtained are summarized as follows:

  1. Significant differences were observed between all measured characteristics in WNT families, except for height, and in the C-P families, except for height and P, indicating that these characteristics were influenced by genetic factors.

  2. The family heritability of all traits was much higher in the WNT than in C-P families, demonstrating (1) much higher variability of all traits in the WNT was due to genetic differences among families; (2) greater additive effect passed on to offspring from WNT families than C-P families.

  3. Positive genetic correlations between growth traits (D, H, V) were found in the WNT and C-P families. The observed moderate to strong negative genetic correlation between D and P, as well as between P and SWV, suggested that using D as a selection criterion could be a suitable breeding strategy for F. variegata. This method allows for the simultaneous genetic improvement of both growth and wood quality traits in the species.

CONFLICT of INTEREST

No potential conflict of interest relevant to this article was reported.

ACKNOWLEDGMENT

The authors are grateful to the Environment and Forestry Agency of Yogyakarta Special Province for fully supporting the genetic materials. We would also like to thank Dr. Christopher Beadle, who reviewed an earlier draft of the paper and made suggestions for improvement.

REFERENCES

1.

Augustina, S., Wahyudi, I., Darmawan, I.W., Malik, J., Basri, E., Kojima, Y. 2020. Specific gravity and dimensional stability of boron-densified wood on three lesser-used species from Indonesia. Journal of the Korean Wood Science and Technology 48(4): 458-471.

2.

Bappeda Kabupaten Bantul. 2011. RPJMD Kabupaten Bantul 2011-2015. Bappeda Kabupaten Bantul, Yogyakarta, Indonesia.

3.

Chaerani, N., Sudrajat, D.J., Siregar, I.Z., Siregar, U.J. 2019. Growth performance and wood quality of white jabon (Neolamarckia cadamba) progeny testing at Parung Panjang, Bogor, Indonesia. Biodiversitas Journal of Biological Diversity 20(8): 2295-2301.

4.

Cheung, B.Y. 2020. Genetic Coefficient of Variance. In: Encyclopedia of Personality and Individual Differences, Ed. by Zeigler-Hill, V. and Shackelford, T.K. Springer, Cham, Switzerland.

5.

de Oliveira, A.L.B., Gouvêa, L.R.L., Verardi, C.K., Silva, G.A.P., de Gonçalves, P.S. 2015. Genetic variability and predicted genetic gains for yield and laticifer system traits of rubber tree families. Euphytica 203(2): 285-293.

6.

Effendi, R., Mindawati, N. 2015. Budidaya Jenis Pohon Nyawai (Ficus variegata Blume.). Kementerian Lingkungan Hidup dan Kehutanan, Jakarta, Indonesia.

7.

Fadwati, A.D., Hidayati, F., Na’iem, M. 2023. Evaluation of genetic parameters of growth characteristics and basic density of Eucalyptus pellita clones planted at two different sites in East Kalimantan, Indonesia. Journal of the Korean Wood Science and Technology 51(3): 222-237.

8.

Falconer, D.S. 1989. Introduction to Quantitative Genetics. 3rd ed. John Wiley & Sons, New York, NY, USA.

9.

Fitriani, A. 2011. Pengaruh ruang tumbuh terhadap respon pertumbuhan dan perkembangan tanaman meranti merah (Shorea pauciflora King.) dan nyawai (Ficus variegata Blum.). Jurnal Hutan Tropis 12(31): 115-122.

10.

Fukatsu, E., Hiraoka, Y., Matsunaga, K., Tsubomura, M., Nakada, R. 2015. Genetic relationship between wood properties and growth traits in Larix kaempferi obtained from a diallel mating test. Journal of Wood Science 61(1): 10-18.

11.

Fukatsu, E., Tamura, A., Takahashi, M., Fukuda, Y., Nakada, R., Kubota, M., Kurinobu, S. 2011. Efficiency of the indirect selection and the evaluation of the genotype by environment interaction using Pilodyn for the genetic improvement of wood density in Cryptomeria japonica. Journal of Forest Research 16(2): 128-135.

12.

Fundova, I., Funda, T., Wu, H.X. 2018. Non-destructive wood density assessment of Scots pine (Pinus sylvestris L.) using Resistograph and Pilodyn. PLOS ONE 13(9): e0204518.

13.

Ghani, R.S.M., Lee, M.D. 2021. Challenges of wood modification process for plantation eucalyptus: A review of Australian setting. Journal of the Korean Wood Science and Technology 49(2): 191-209.

14.

Haapanen, M. 2001. Time trends in genetic parameter estimates and selection efficiency for Scots pine in relation to field testing method. Forest Genetics 8(2): 129-144.

15.

Haryjanto, L., Hadiyan, Y. 2014. Eksplorasi benih nyawai (Ficus variegata Blume) di Kecamatan Long Hubung, Kabupaten Kutai Barat, Kalimantan Timur. Wana Benih 15(2): 61-72.

16.

Haryjanto, L., Prastyono, P., Yuskianti, V. 2014. Variasi pertumbuhan dan parameter genetik pada tiga plot uji keturunan nyawai (Ficus variegata Blume) di Bantul. Jurnal Pemuliaan Tanaman Hutan 8(3): 137-151.

17.

Hendromono, Komsatun. 2008. Nyawai (Ficus variegata Blume dan Ficus sycomoroides Miq.) jenis yang berprospek baik untuk dikembangkan di hutan tanaman. Mitra Hutan Tanaman 3(3): 122-130.

18.

Hidayati, F., Ishiguri, F., Iizuka, K., Makino, K., Takashima, Y., Danarto, S., Winarni, W.W., Irawati, D., Na’iem, M., Yokota, S. 2013a. Variation in tree growth characteristics, stress-wave velocity, and Pilodyn penetration of 24-year-old teak (Tectona grandis) trees originating in 21 seed provenances planted in Indonesia. Journal of Wood Science 59(6): 512-516.

19.

Hidayati, F., Ishiguri, F., Iizuka, K., Makino, K., Tanabe, J., Marsoem, S.N., Na’iem, M., Yokota, S., Yoshizawa, N. 2013b. Growth characteristics, stress-wave velocity, and Pilodyn penetration of 15 clones of 12-year-old Tectona grandis trees planted at two different sites in Indonesia. Journal of Wood Science 59(3): 249-254.

20.

Hidayati, F., Ishiguri, F., Makino, K., Tanabe, J., Aiso, H., Prasetyo, V.E., Marsoem, S.N., Wahyudi, I., Iizuka, K., Yokota, S. 2017. The effects of radial growth rate on wood properties and anatomical characteristics and an evaluation of the xylem maturation process in a tropical fast-growing tree species, Gmelina arborea. Forest Products Journal 67(3-4): 297-303.

21.

Hidayati, F., Lukmandaru, G., Indrioko, S., Sunarti, S., Nirsatmanto, A. 2019. Variation in tree growth characteristics, Pilodyn penetration, and stress-wave velocity in 65 families of Acacia mangium trees planted in Indonesia. Journal of the Korean Wood Science and Technology 47(5): 633-643.

22.

Hill, W.G. 2010. Understanding and using quantitative genetic variation. Philosophical Transactions of the Royal Society B: Biological Sciences 365(1537): 73-85.

23.

Hong, Z., Fries, A., Wu, H.X. 2014. High negative genetic correlations between growth traits and wood properties suggest incorporating multiple traits selection including economic weights for the future Scots pine breeding programs. Annals of Forest Science 71(4): 463-472.

24.

Houle, D. 1992. Comparing evolvability and variability of quantitative traits. Genetics 130(1): 195-204.

25.

Ishiguri, F., Eizawa, J., Saito, Y., Iizuka, K., Yokota, S., Priadi, D., Sumiasri, N., Yoshizawa, N. 2007. Variation in the wood properties of Paraserianthes falcataria planted in Indonesia. IAWA Journal 28(3): 339-348.

26.

Ishiguri, F., Matsui, R., Iizuka, K., Yokota, S., Yoshizawa, N. 2008. Prediction of the mechanical properties of lumber by stress-wave velocity and Pilodyn penetration of 36-year-old Japanese larch trees. Holz Als Roh- Und Werkstoff 66(4): 275-280.

27.

Jin, J., Zhao, X., Liu, H., Wang, S., Song, Z., Ma, X., Li, K. 2019. Preliminary study on genetic variation of growth traits and wood properties and superior clones selection of Populus ussuriensis Kom. iForest 12(5): 459-466.

28.

Kang, K.S., Cheon, B.H., Han, S.U., Kim, C.S., Choi, W.Y. 2007. Genetic gain and diversity under different selection methods in a breeding seed orchard of Quercus serrata. Silvae Genetica 56(6): 277-281.

29.

Kha, L.D., Harwood, C.E., Kien, N.D., Baltunis, B.S., Hai, N.D., Thinh, H.H. 2012. Growth and wood basic density of acacia hybrid clones at three locations in Vietnam. New Forests 43(1): 13-29.

30.

Kien, N.D., Jansson, G., Harwood, C., Almqvist, C., Thinh, H.H. 2008. Genetic variation in wood basic density and pilodyn penetration and their relationships with growth, stem straightness, and branch size for Eucalyptus urophylla in northern Vietnam. New Zealand Journal of Forestry Science 38(1): 160-175.

31.

Klápště, J., Telfer, E.J., Dungey, H.S., Graham, N.J. 2022. Chasing genetic correlation breakers to stimulate population resilience to climate change. Scientific Reports 12(1): 1-16.

32.

Lai, M., Dong, L., Yi, M., Sun, S., Zhang, Y., Fu, L., Xu, Z., Lei, L., Leng, C., Zhang, L. 2017. Genetic variation, heritability and genotype × environment interactions of resin yield, growth traits and morphologic traits for Pinus elliottii at three progeny trials. Forests 8(11): 409.

33.

Laksono, G.D., Rahayu, I.S., Karlinasari, L., Darmawan, W., Prihatini, E. 2023. Characteristics of magnetic sengon wood impregnated with nano Fe3O4 and furfuryl alcohol. Journal of the Korean Wood Science and Technology 51(1): 1-13.

34.

Matheson, A.C., Gapare, W.J., Ilic, J., Wu, H.X. 2008. Inheritance and genetic gain in wood stiffness in radiata pine assessed acoustically in young standing trees. Silvae Genetica 57(1-6): 56-64.

35.

Mihai, G., Mirancea, I., Duta, C. 2014. Variation of the quantitative traits in a progeny test of Abies alba (Mill.) at the nursery stage. Silvae Genetica 63(6): 275-284.

36.

Nelson, C.D. 2023. Tree breeding, a necessary complement to genetic engineering. New Forests 54(4): 721-738.

37.

Nguyen, H.T.H., Chen, Z.Q., Fries, A., Berlin, M., Hallingbäck, H.R., Wu, H.X. 2022. Effect of additive, dominant and epistatic variances on breeding and deployment strategy in Norway spruce. Forestry 95(3): 416-427.

38.

Nur’aini, Syamsuardi, Arbain, A. 2013. Tumbuhan Ficus L. (Moraceae) di hutan konservasi Prof. Soemitro Djojohadikusumo, PT. Tidar Kerinci Agung (TKA), Sumatera Barat. Jurnal Biologi Universitas Andalas 2(4): 235-241.

39.

Nurtjahjaningsih, I.L.G., Haryjanto, L., Sulistyawati, P., Widyatmoko, A.Y.P.B.C. 2019. Mating system of two priority species for ecosystem restoration at Mount Merapi. AIP Conference Proceedings 2120: 040028.

40.

Pitopang, R. 2012. Impact of forest disturbance on the structure and composition of vegetation in tropical rainforest of Central Sulawesi, Indonesia. Biodiversitas Journal of Biological Diversity 13(4): 178-189.

41.

Qirom, M.A., Supriyadi, S. 2013. Model penduga volume pohon nyawai (Ficus variegata Blume) di Kalimantan Timur. Jurnal Penelitian Hutan Tanaman 10(4): 173-284.

42.

Ranjan, S., Gautam, A. 2020. Heritability Estimate. In: Encyclopedia of Animal Cognition and Behavior, Ed. by Vonk, J. and Shakelford, T. Springer, Cham, Switzerland.

43.

Rofifah, I., Indriyanto, Asmarahman, C. 2021. Species and benefits of Ficus spp. in the collection block of Tahura WAR Lampung Province. Jurnal Rimba Lestari 1(2): 88-98.

44.

Schimleck, L., Dahlen, J., Apiolaza, L.A., Downes, G., Emms, G., Evans, R., Moore, J., Pâques, L., Van den Bulcke, J., Wang, X. 2019. Non-destructive evaluation techniques and what they tell us about wood property variation. Forests 10(9): 728.

45.

Schmidt, F.H., Ferguson, J.H.A. 1951. Rainfall Types Based on Wet and Dry Period Ratios for Indonesia with Western New Guinea. Kementerian Perhubungan, Jakarta, Indonesia.

46.

Seta, G.W., Widiyatno, Hidayati, F., Na’iem, M. 2021. Impact of thinning and pruning on tree growth, stress wave velocity, and Pilodyn penetration response of clonal teak (Tectona grandis) plantation. Forest Science and Technology 17(2): 57-66.

47.

Spencer, H., Weiblen, G., Flick, B. 1996. Phenology of Ficus variegata in a seasonal wet tropical forest at Cape Tribulation, Australia. Joumal of Biogeography 23(4): 467-475.

48.

Sumarni, G., Muslich, M., Hadjib, N., Krisdianto, Malik, D., Suprapti, S., Basri, E., Pari, G., Iskandar, M.I., Siagian, R.M. 2009. Sifat dan Kegunaan Kayu: 15 Jenis Kayu Andalan Setempat Jawa Barat. Pusat Penelitian dan Pengembangan Hasil Hutan, Bogor, Indonesia.

49.

Susanti, P.D., Halwany, W. 2017. Dekomposisi serasah dan keanekaragaman makrofauna tanah pada hutan tanaman industri nyawai (Ficus variegata Blume). Jurnal Ilmu Kehutanan 11(2): 212-223.

50.

Thomas, E., Jalonen, R., Loo, J., Boshier, D., Gallo, L., Cavers, S., Bordács, S., Smith, P., Bozzano, M. 2014. Genetic considerations in ecosystem restoration using native tree species. Forest Ecology and Management 333: 66-75.

51.

Van Duong, D., Schimleck, L., Lam Tran, D., Dai Vo, H. 2022. Radial and among-clonal variations of the stress-wave velocity, wood density, and mechanical properties in 5-year-old Acacia auriculiformis clones. BioResources 17(2): 2084-2096.

52.

van Heist, M., Sheil, D., Rachman, I., Gusbager, P., Raweyai, C.O., Yoteni, H.S.M. 2010. The forests and related vegetation of Kwerba, on the Foja foothills, Mamberamo, Papua (Indonesian New Guinea). Blumea 55(2): 153-161.

53.

Wahyuningtyas, R.S., Junaidah, J., Santosa, P.B. 2022. Response of Ficus variegata seedling size on their early growth in imperata grassland. IOP Conference Series: Earth and Environmental Science 959: 012012.

54.

Wang, X., Divos, F., Pilon, C., Brashaw, B.K., Ross, R.J., Pellerin, R.F. 2004. Assessment of Decay in Standing Timber Using Stress Wave Timing Nondestructive Evaluation Tools: A Guide for Use and Interpretation. U.S. Department of Agriculture Forest Service, Washington, DC, USA.

55.

Wang, X., Ross, R.J., McClellan, M., Barbour, R.J., Erickson, J.R., Forsman, J.W., McGinnis, G.D. 2000. Strength and Stiffness Assessment of Standing Trees Using a Nondestructive Stress Wave Technique. U.S. Department of Agriculture Forest Service, Washington, DC, USA.

56.

Wang, X., Ross, R.J., McClellan, M., Barbour, R.J., Erickson, J.R., Forsman, J.W., McGinnis, G.D. 2001. Nondestructive evaluation of standing trees with a stress wave method. Wood and Fiber Science 33(4): 522-533.

57.

Waszak, M., Cieślik, K. 2016. Analysis of the validity of assumptions underlying a research on the heritability of quantitative traits. Anthropologischer Anzeiger 73(1): 61-68.

58.

Wei, X., Borralho, N.M.G. 1997. Genetic control of wood basic density and bark thickness and their relationships with growth traits of Eucalyptus urophylla in south east China. Silvae Genetica 46(4): 245-250.

59.

Wessels, C.B., Malan, F.S., Rypstra, T. 2011. A review of measurement methods used on standing trees for the prediction of some mechanical properties of timber. European Journal of Forest Research 130(6): 881-893.

60.

White, T.L., Adams, W.T., Neale, D.B. 2007. Forest Genetics. CABI, Wallingford, UK.

61.

Wu, S., Xu, J., Li, G., Risto, V., Lu, Z., Li, B., Wang, W. 2010. Use of the Pilodyn for assessing wood properties in standing trees of Eucalyptus clones. Journal of Forestry Research 21(1): 68-72.

62.

Xiao, Y., Wang, J., Yun, H., Yang, G., Ma, J., Ma, W., Qu, G. 2021. Genetic evaluation and combined selection for the simultaneous improvement of growth and wood properties in Catalpa bungei clones. Forests 12(7): 868.

63.

Yasin, M., Noor, A., Ningsih, R.D. 2020. Opportunities for improving rice production in sub-obtimal lands South Kalimantan. IOP Conference Series: Earth and Environmental Science 484(1): 012063.

64.

Zobel, B., Talbert, J. 1984. Applied Forest Tree Improvement. John Wiley & Sons, New York, NY, USA.