Genetic variation and heritability estimates of Ulmus minor and Ulmus pumila hybrids for budburst, growth and tolerance to Ophiostoma novo-ulmi
iForest - Biogeosciences and Forestry, Volume 8, Issue 4, Pages 422-430 (2014)
doi: https://doi.org/10.3832/ifor1227-007
Published: Dec 15, 2014 - Copyright © 2014 SISEF
Research Articles
Collection/Special Issue: 3rd International Elm Conference, Florence (Italy - 2013)
The elms after 100 years of Dutch Elm disease
Guest Editors: A. Santini, L. Ghelardini, E. Collin, A. Solla, J. Brunet, M. Faccoli, A. Scala, S. De Vries, J. Buiteveld
Abstract
Seedlings obtained by crossing Ulmus minor and U. minor × U. pumila clones were assessed for flowering, bark beetle damage, vegetative budburst, height growth and resistance to Ophiostoma novo-ulmi. Ramets and open pollinated seedlings obtained from the parent trees were assessed for the same traits. Most progenies had similar traits to their parents, but some presented heterosis in annual growth or resistance to O. novo-ulmi. Leaf wilting was significantly lower in progenies with U. minor × U. pumila rather than U. minor as female parent (21.5 and 30.6%, respectively; P<0.05). Resistance to O. novo-ulmi increased significantly as a function of increased amounts of U. pumila germplasm from the female parent, suggesting that resistance to Dutch elm disease is primarily transmitted from the mother. Budburst occurred earlier in seedlings with low rather than high growth rates (P=0.0007) and percentage of wilting was negatively related to early budburst (P<0.0001). Other phenotypic relations included percentage of flowering trees and annual height growth (rp=0.44; P=0.0042), percentage of flowering trees and vegetative budburst (rp=-0.53; P=0.0004) and percentage of beetle-affected trees and annual height growth (rp=0.60; P<0.0001). Heritability estimates obtained from the regression and variance components methods ranged from 0.06 ± 0.04 to 0.64 ± 0.18, 0.10 ± 0.05 to 0.69 ± 0.17, and 0.13 ± 0.32 to 0.71 ± 0.22 for budburst, growth and tolerance to O. novo-ulmi, respectively. Broad- and narrow-sense heritability values were higher when estimated 60 days post inoculation (dpi) than 15, 30 or 120 dpi. Heritability estimates and genetic gains reported indicate a high degree of additive genetic control and show the effectiveness of selection for Dutch elm disease resistance and rapid tree growth.
Keywords
Quantitative Genetics, Tree Breeding, Invasive Pathogen, Inheritance, Additive Genetic Variance, Non-additive Genetic Variance, Heterosis, Genetic Gain
Introduction
The field elm (Ulmus minor) was formerly a major component of riparian forests and municipal planting in southern Europe. Valuable ornamental traits, adaptability to urban environments and simplicity of establishment and management are the notable features of this species. However, U. minor has virtually disappeared from urban landscaping due to the spread of the lethal Dutch elm disease (DED), which is caused by the pathogen Ophiostoma novo-ulmi ([2], [38]). The fungus is transmitted from diseased to healthy elms by elm bark beetles, in particular Scolytus scolytus and S. multistriatus, the two main vectors in southern Europe. Spray programmes to control insect vectors are now largely unacceptable. Systemic chemicals have proven too toxic for general use ([20], [22], [23]) and biological control treatments appear to be ineffective ([35], [42]). As a result, replanting efforts now focus on other tree genera.
From a genetic point of view, re-introducing elms into municipal landscaping appears to be both desirable and feasible. A genetic basis for disease resistance to DED has been established in Europe using native elm species and clones, especially in hybrid combinations between and within U. minor and U. pumila ([27], [25]). Ulmus pumila is not a European native and was probably introduced into Spain in the 16th century as an ornamental and into Italy in the 1930s to replace dying elms. Natural hybridization with U. minor trees in these countries has produced genetic introgression, leading to a complex taxonomy ([4]).
Screening has been used to identify European elms with low disease susceptibility ([29], [33], [36]). These studies have identified U. pumila individuals and U. minor × U. pumila hybrids with low susceptibility and rapid height growth. As ornamentals, however, U. pumila and U. minor × U. pumila have several features that compare unfavorably with U. minor, the species they are intended to replace. Excessively fine branching, excurrent habit, small leaves, branch brittleness and twig mortality are normally undesirable features in urban environmental planting ([11], [15]). The Spanish elm breeding programme has therefore begun a series of backcrosses using U. minor as the source of improved ornamental traits and U. minor × U. pumila as the source of disease resistance and rapid height growth. Fast growth is definitely a desirable trait for the common landscape use of trees in southern Europe. Delayed vegetative budburst is similarly desirable in areas with spring frosts or where early defoliating insects such as Euproctis chrysorrhoea cause considerable leaf damage. Flower and seed production are important characters to study in an F1 population, especially if successive breeding cycles are to be performed. Information on the incidence of elm bark beetles on crossings would allow breeders to rule out combinations that are most affected by this pest.
Phenotypic variation can be divided into two major components, one representing genetic factors and the other representing environmental influences or experimental errors (i.e., stochastic variation). Broad-sense heritability (H2) is the proportion of variance in a particular trait and population that is due to genetic components (Fig. 1). A common simplification in all genetic studies and models is to assume that all alleles and all genotypes act independently of each other. This type of model, known as an “additive model”, does not allow for dominant or recessive effects and gene × gene interactions, even though they are known to be numerous. Broad-sense heritability can be further divided into additive genetic components and gene × gene interaction components (Fig. 1). The contribution of additive genetic components is termed narrow-sense heritability (h2). In other words, h2 is the proportion of trait variance due to additive genetic factors, whereas H2 is the proportion of trait variance due to all the genetic effects ([5]). The difference between these two heritability parameters gives an estimate of the effect of interactions on a phenotype. If the difference is zero, it means there is no interaction and all phenotypic variations can be explained by additive components. A non-zero difference between the two heritability parameters suggests the role and degree of interaction effects on the phenotype ([1]). Heritability estimates of vegetative budburst, height growth and resistance to O. novo-ulmi have not previously been reported for a breeding population of U. minor and U. pumila elms. In this study it is hypothesized that heritability estimates of previous traits will be high enough to justify selection and breeding. The study was designed to estimate genetic variation and inheritance of U. minor and U. minor × U. pumila hybrids for budburst, growth and tolerance to O. novo-ulmi.
Fig. 1 - Components of phenotypic variation, including broad-sense heritability (H2) and narrow-sense heritability (h2).
Material and Methods
Plant material
Ten clones - eight Ulmus minor (m) and two U. minor × U. pumila (mp) hybrids - were used as parents for the mating design (Tab. 1). Each clone comprised one plant. Parents were randomly selected from the elm clone bank population at the Forest Breeding Center in Puerta de Hierro, Madrid (see Appendix 1). This bank includes more than 200 Iberian elm clones, obtained through vegetative (root cuttings) or sexual reproduction from elms selected in the countryside according to their good sanitary status, diversity (geographical, ecological and genetic) and distinctive ornamental traits. The 10 clones selected were characterized through isozyme markers before crossing ([4]).
Tab. 1 - Clones and number of controlled pollinated seedlings used in the partial diallel cross design. (m): Ulmus minor; (mp): U. minor × U. pumila; (-): controlled pollinated crosses not performed; (0): controlled pollination performed but no seedlings obtained.
Female | Males | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
AL-AN1 (m) | AL-JR1 (m) | GR-AL3 (m) | M-PZ3 (mp) | SE-CT4 (m) | VA-VV5 (mp) | TO-AL1 (m) | CC-VG2 (m) | M-IN5 (m) | V-JR1 (m) | |
AL-AN1 (m) | - | 6 | 15 | 17 | 9 | - | - | - | - | - |
AL-JR1 (m) | - | - | 2 | 9 | 6 | 17 | - | - | - | - |
GR-AL3 (m) | - | - | - | 17 | 20 | 22 | 15 | - | - | - |
M-PZ3 (mp) | - | - | - | - | 29 | 23 | 6 | 11 | - | - |
SE-CT4 (m) | - | - | - | - | - | 17 | 38 | 0 | 10 | - |
VA-VV5 (mp) | - | - | - | - | - | - | 22 | 18 | 28 | 23 |
TO-AL1 (m) | 5 | - | - | - | - | - | - | 6 | 18 | 4 |
CC-VG2 (m) | 0 | 0 | - | - | - | - | - | - | 0 | 0 |
M-IN5 (m) | 0 | 0 | 0 | - | - | - | - | - | - | 0 |
V-JR1 (m) | 0 | 0 | 0 | 0 | - | - | - | - | - | - |
Selected clones were crossed in 2000 following a partial diallel cross design ([3]). Forty cross-combinations were planned (Tab. 1), but 13 failed to provide full-sib progeny (combinations marked “0” in Tab. 1) due to the presence of three previously undetected aborting U. minor individuals (CC-VG2, M-IN5 and V-JR1) ([18]). Pollination was conducted through traditional tree breeding techniques using greenhouse-forced pollen (Fig. 2a) and isolation bags on trees ([26]). When a clone was used as female parent, the flowers were chosen on different branches of the same tree and bagged separately. Pollen from one male parent was used to pollinate the flowers of one or two bagged branches. To protect seeds from birds, bags were not removed from the branches until mature seeds had been collected. Isolated flowers were pollinated by inserting a syringe into the bag (Fig. 2b), after checking pollen viability by fluorochromatic reaction ([10]). Flowers in each bag were pollinated twice when possible, allowing four to five days between pollinations.
Fig. 2 - Methods: (a) forced pollen sampling from greenhouse flowering twigs in February 2000; (b) controlled pollinations performed by insufflating pollen into isolation bags; (c) plantation of the experimental plot in March 2002; (d) tree growth assessment in October 2004; (e) artificial inoculations with Ophiostoma novo-ulmi in May 2005; (f) leaf wilting assessment in July 2005.
Seeds were sown outdoors in a mixture of sand, perlite and peat (1:1:1) and one-year-old dormant seedlings were planted in a randomized two-block field design. Each block had one to six experimental units per cross, and each experimental unit (i.e., subsample) comprised three plants. The number of experimental units for each entry within each block was the same. The field design (Fig. 2c) was installed at the Center in Madrid. Spacing was 1.5 m between rows and 0.5 m between plants within rows. An elm line border was included to avoid side effects. To ensure plant growth, the plot was weeded annually, the soil was fertilized once a year (Osmocote© plus, Scotts, Australia) and plants were watered weekly in summer.
In 2000, ramets of six clones (AL-AN1, AL-JR1, GR-AL3, M-PZ3, TO-AL1, and V-JR1) were obtained through root cuttings. Cuttings were obtained inside individual 2-l pots containing a fine sandy loam texture soil and peat (80:20, v/v). Plants from the root stocks of the same clone were then transplanted within the plot containing the cross-pollinated seedlings. At least seven ramets (mean = 12) per clone were transplanted randomly within the two blocks of the field design. Additionally, seeds (half-sib progeny) from seven open pollinated clones (AL-AN1, AL-JR1, GR-AL3, M-PZ3, TO-AL1, SE-CT4, and VA-VV5) were collected in 2000, sown in the same conditions and planted within the plot. At least 15 open pollinated seedlings (mean = 18) per clone were planted.
Tree growth, budburst, flowering and fructification assessments
All seedlings and ramets shown in Tab. 2 were assessed. Tree growth was recorded by measuring tree height in October 2004 (Fig. 2d) and dividing by tree age. Phenology was assessed from February 7 to April 9 2005. Date of vegetative budburst was taken on 30 buds per tree. Bud development was assessed weekly as follows ([32]): 1= dormant buds; 2= swollen buds, but scales closed; 3= bud scales open and extremities of the first leaf visible at the apex of the buds; 4= extremities of all leaves out; and 5= two or more leaves completely expanded. Flowering was quantitatively assessed on March 4 2005, when most staminal filaments had begun extending, from 0= no floral buds to 5= the whole crown showing floral buds. Fruit development was monitored macroscopically using 30 samaras per tree, when available.
Tab. 2 - Mean values of five traits measured on ramets, open pollinated (op) seedlings, and controlled pollinated seedlings from crosses (female × male) between Ulmus minor (AL-AN1, AL-JR1, GR-AL3, SE-CT4, TO-AL1, CC-VG2, M-IN5 and V-JR1) and U. minor × U. pumila (M-PZ3 and VA-VV5) clones. Different letters indicate differences between means (P < 0.01). (a): trees with beetle feeding grooves on trunk or branches; (b): assessed on April 9 2005 on a scale from 1= dormant buds to 5= two or more leaves completely expanded; (c): assessed on July 7 2005, 60 days after inoculations with Ophiostoma novo-ulmi.
Plant material | No. trees | Flowering trees (%) |
Beetle affected trees (%)a |
Vegetative budburstb |
Annual height growth (cm) |
Leaf wilting (%)c |
---|---|---|---|---|---|---|
AL-AN1 ramets | 7 | 0.0 | 0.0 | 4.6 bcde | 21.6 a | 80.0 e |
AL-AN1 op | 11 | 0.0 | 0.0 | 4.4 bcde | 32.2 ab | 49.1 e |
AL-AN1 × AL-JR1 | 6 | 33.3 | 16.7 | 4.8 cde | 37.3 bc | 23.3 abcd |
AL-AN1 × GR-AL3 | 15 | 0.0 | 0.0 | 4.4 bcde | 40.4 bcd | 35.0 bcde |
AL-AN1 × M-PZ3 | 17 | 0.0 | 17.6 | 4.1 bc | 45.0 bcde | 38.1 cde |
AL-AN1 × SE-CT4 | 9 | 0.0 | 11.1 | 4.8 de | 35.1 abc | 41.3 de |
AL-JR1 ramets | 11 | 0.0 | 9.1 | 4.9 e | 25.8 a | 75.0 e |
AL-JR1 op | 12 | 16.7 | 25.0 | 4.1 bcd | 62.3 ef | 36.3 cde |
AL-JR1 × GR-AL3 | 2 | 0.0 | 0.0 | 5.0 cde | 14.0 a | 20.0 ab |
AL-JR1 × M-PZ3 | 9 | 0.0 | 44.4 | 4.8 de | 58.5 def | 35.0 bcde |
AL-JR1 × SE-CT4 | 6 | 0.0 | 16.7 | 4.4 bcde | 38.3 abcd | 30.0 abcde |
AL-JR1 × VA-VV5 | 17 | 5.9 | 17.6 | 4.6 cde | 53.1 def | 43.2 e |
GR-AL3 ramets | 13 | 0.0 | 0.0 | 5.0 cde | 24.0 a | 2.8 a |
GR-AL3 op | 10 | 0.0 | 50.0 | 4.9 e | 84.0 ef | 55.0 ab |
GR-AL3 × M-PZ3 | 17 | 0.0 | 23.5 | 4.3 bcde | 45.9 bcde | 20.3 ab |
GR-AL3 × SE-CT4 | 20 | 25.0 | 40.0 | 4.6 cde | 58.5 ef | 16.5 a |
GR-AL3 × VA-VV5 | 22 | 0.0 | 9.1 | 3.8 ab | 49.1 cde | 14.5 a |
GR-AL3 × TO-AL1 | 15 | 6.7 | 0.0 | 4.8 e | 33.4 ab | 32.2 bcde |
M-PZ3 ramets | 12 | 75.0 | 16.7 | 2.9 a | 75.0 f | 28.8 abcd |
M-PZ3 op | 24 | 4.2 | 12.5 | 4.2 bc | 48.5 cde | 24.0 abc |
M-PZ3 × SE-CT4 | 29 | 3.4 | 17.4 | 4.8 e | 50.0 cde | 19.1 ab |
M-PZ3 × VA-VV5 | 23 | 0.0 | 17.4 | 4.1 bc | 47.3 cde | 17.7 a |
M-PZ3 × TO-AL1 | 6 | 16.7 | 50.0 | 3.1 a | 65.3 ef | 48.3 e |
M-PZ3 × CC-VG2 | 11 | 9.1 | 18.2 | 4.4 bcde | 55.0 def | 23.8 abcd |
SE-CT4 op | 14 | 0.0 | 0.0 | 4.7 bcde | 28.0 ab | 18.8 a |
SE-CT4 × VA-VV5 | 17 | 5.9 | 5.9 | 4.4 bcde | 39.2 abcd | 40.0 de |
SE-CT4 × TO-AL1 | 38 | 0.0 | 0.0 | 4.5 cde | 34.6 ab | 38.8 de |
SE-CT4 × M-IN5 | 10 | 0.0 | 0.0 | 4.0 bc | 39.9 abcd | 26.7 abcde |
VA-VV5 op | 2 | 0.0 | 0.0 | 5.0 cde | 75 f | 10.0 abc |
VA-VV5 × TO-AL1 | 22 | 4.5 | 18.2 | 4.9 e | 44.8 bcde | 25.2 abcd |
VA-VV5 × CC-VG2 | 18 | 22.2 | 11.1 | 4.5 bcde | 65.3 f | 17.0 ab |
VA-VV5 × M-IN5 | 28 | 14.3 | 28.6 | 4.7 de | 62.1 f | 15.0 a |
VA-VV5 × V-JR1 | 23 | 13.0 | 17.4 | 4.8 e | 56.9 def | 17.9 ab |
TO-AL1 ramets | 7 | 14.3 | 14.3 | 3.7 ab | 44.7 bcde | 40.0 cd |
TO-AL1 op | 5 | 20.0 | 20.0 | 4.8 cde | 37.6 abcd | 30.0 abcde |
TO-AL1 × AL-AN1 | 5 | 0.0 | 40.0 | 4.7 cde | 44.2 bcde | 20.0 ab |
TO-AL1 × CC-VG2 | 6 | 0.0 | 0.0 | 4.8 cde | 47.0 bcde | 43.3 de |
TO-AL1 × M-IN5 | 18 | 5.6 | 0.0 | 4.6 cde | 36.3 abc | 22.1 abc |
TO-AL1 × V-JR1 | 4 | 0.0 | 0.0 | 5.0 de | 40.0 abc | 7.0 ab |
V-JR1 ramets | 8 | 0.0 | 0.0 | 4.7 cde | 24.1 ab | 40.0 de |
Average | 17 | 7.3 | 14.2 | 4.4 | 45.5 | 37.1 |
Resistance to Ophiostoma novo-ulmi and beetle wound incidence
Plants smaller than 0.7 m or affected by natural infections were not inoculated. All other plants were inoculated on May 7th 2005 (Fig. 2e). The number of seedlings and ramets inoculated is shown in Tab. 2. The pathogen, Ophiostoma novo-ulmi subsp. americana Brasier, was isolated in 2002 in Blacos (Soria, Spain; code SO-BL1) from an infected U. minor tree. This strain was selected because of its high virulence and in vitro growth (5.8 mm day-1 on MEA at 20 °C - [38]). The inoculum consisted of a bud-cell suspension of the fungus grown in Tchernoff’s medium ([40]). The suspension was centrifuged to eliminate the medium and then filtered, and the precipitate was adjusted in sterile distilled water to 104 spores ml-1 ([19]). At midday, a transversal cut was made with a flame-sterilized scalpel blade in the sapwood of the trees, 8-10 cm above the soil surface, on the north side of each tree (Fig. 2e - [37]). Two 0.1 ml drops of bud-cell suspension were applied above the blade using a hypodermic syringe and the inoculum was immediately drawn into the xylem. The percentage of the crown showing wilting or death of the foliage was visually estimated 15, 30, 60 and 120 dpi (days post inoculation), using a 5% interval.
A single scolytid trap was installed about 300 m from the plots, as described by Solla & Gil ([35]). The trap was checked weekly in 2005. On June 15th and 16th, two weeks after the peak of captured beetles, the number of trunk feeding grooves was quantified for each tree: 0= no beetle grooves present; 1= at least one beetle groove; 2= abundant beetle grooves on the trunk; and 3= abundant beetle grooves on the trunk and branches.
Statistical analyses
Mean values of traits were compared in two stages. Firstly, mean values were analysed at clone (ramet) and family levels. Secondly, the plant material was divided into groups according to the genealogy or propagation method: ramets, seedlings from open pollination (half-sib progenies) and seedlings from controlled pollination (full-sib progenies). The full-sib progenies were further divided according to parentage: crosses between clones of U. minor (m × m in Tab. 1), crosses between clones from U. minor × U. pumila (mp × mp in Tab. 1) and crosses between m and mp clones (m × mp and mp × m in Tab. 1). Angular transformation of wilting percentages (x) was performed to normalize data before statistical analysis [y = (x/100)1/2]. Vegetative budburst notations, plant growth and transformed wilting values were analyzed by multifactor ANOVA, considering the block and the plant entries as factors. Post-hoc Fisher’s LSD tests were applied to compare averages at P < 0.05 and P < 0.01.
Heritability estimates for vegetative budburst, plant growth and resistance to O. novo-ulmi were calculated in two different ways. Firstly, the regression coefficients between parent (ramets) means and offspring means were computed. Both “ramets vs half-sib progenies” and “ramets vs full-sib progenies” relations were considered. Estimation of heritability by the parent-offspring (P-O) regression method is relatively free of genetic assumptions and is often recommended as an empirical and reliable method for estimating narrow-sense heritabilities (h
2 - [5]). The standard error (SE) for each h
2 estimate was calculated as (eqn. 1):
where n
is the number of controlled pollinated families. Secondly, we estimated variance components and respective standard errors after removing the parental genotypes from the datasets, and analyzed the controlled pollinated seedlings population via the restricted maximum likelihood (REML) method. REML analysis is currently the method of choice for estimating heritability and variance components, as it does not require a balanced design ([5]). Variance components for broad-sense heritability estimates were made under a random model by treating the variables family and block as random factors. Broad-sense heritability (H
2) estimates were calculated on a single-plant basis via the following equation ([6] - eqn. 2):
where σ
2g is the genetic variance and σ
2e the experimental error variance. Broad-sense heritability estimates were calculated on a family mean basis using the equation ([6] - eqn. 3):
where r
refers to the average number of seedlings for each family. Approximate SE of the H
2 estimates were calculated as SE [H
2(plant basis)] = SE(σ
2g) / (σ
2g + σ
2e) and SE [H
2 (mean basis)]= SE(σ
2g) / [σ
2g + (σ
2e / r
)], respectively ([8], [28]). Variance components for narrow-sense heritability estimates were analyzed under a random model by treating the variables male, female, block and male × female interaction as random factors. Narrow-sense heritability estimates were calculated on a single-plant basis using the following equation (eqn. 4):
where (eqn. 5):
and (eqn. 6):
([8]). The term σ
2e again refers to the experimental error variance, but unlike the case of broad-sense heritability formulas, σ
2e is the experimental error excluding the additive and dominance genetic variances. Narrow-sense heritability estimates were calculated on a genotype mean basis using the equation (eqn. 7):
Approximate standard errors (SE) of the h
2 estimates were calculated using the formulas SE[h
2(plant basis)] = 4 SE(σ
2males) / (σ
2add + σ
2dom + σ
2e ) and SE[h
2(mean basis)] = 4 SE(σ
2males) / [σ
2add + σ
2dom + (σ
2e / r
)] ([8]).
Genetic gains (ΔG
) were obtained from ΔG
= h
2(mean basis) × (selection differential), where selection differential is the difference between the mean value of the trait of selected individuals and the mean value of the trait of all individuals ([44]), including those selected. The mean value of the trait of all the individuals was obtained from measurements of this trait in all the open pollinated families (N=7).
Pearson correlation coefficients (phenotypic r
p) among values of vegetative budburst, annual height growth and transformed percentages of wilting were obtained. Statistics were performed using STATGRAPHICS® Plus v. 10, and REML analysis and variance components were obtained using the mixed procedure of GENSTAT® v. 8.
Results
At age five years, only 30 out of 491 seedlings (6.1%) and 10 out of 58 ramets (17.2%) produced flowers, accounting for 7.3% of total trees (Tab. 2). Two out of six U. minor × (U. minor × U. pumila) progenies included flowering trees (mean = 2.0%). However, all seven (U. minor × U. pumila) × U. minor progenies produced at least one flowering tree (mean = 10.2%). Approximately 14.2% of trees had trunk feeding grooves caused by elm bark beetles (Tab. 2).
Mean leaf wilting percentages (± SD) 15, 30, 60 and 120 dpi were 19 ± 20, 34 ± 25, 37 ± 26 and 33 ± 31, respectively. Because vegetative budburst on April 9 and wilting 60 dpi provided the highest F-ratio values in the ANOVAs (results not shown), these measurements were used for further statistical analyses. Means for vegetative budburst, annual height growth and disease response for progeny combinations differed at a probability of 99% for each trait (Tab. 2). Cross-pollinated families with VA-VV5 and M-PZ3 as female parent had the highest annual height growth rates. Lowest height growth occurred when SE-CT4 was female parent (Tab. 3). Height growth was greater when GR-AL3 and AL-JR1 were used as female rather than male parents (Tab. 3). Controlled pollinated families of U. minor showed fast growth when CC-VG2 was used as male parent and low susceptibility when M-IN5 was the male parent (Tab. 3). Most crosses resulted in progenies with similar traits to their parents or intermediate traits. However, AL-AN1 × GR-AL3 and AL-JR1 × GR-AL3 showed heterosis in annual height growth, and AL-AN1 × AL-JR1, M-PZ3 × TO-AL1, TO-AL1 × AL-AN1 and TO-AL1 × V-JR1 showed heterosis in resistance to O. novo-ulmi (Tab. 2). Several elm trees obtained from GR-AL3 × SE-CT4, M-PZ3 × CC-VG2 and VA-VV5 × M-IN5 crossings showed growth rates above 80 cm yr-1 and wilting levels under 20%. Some ramets were especially heterogeneous for certain traits, i.e., M-PZ3 for wilting and TO-AL1 for growth (Tab. 2).
Tab. 3 - Mean values of annual height growth and leaf wilting (60 days after inoculations with Ophiostoma novo-ulmi) of cross-pollinated families with different female and male parents. Within the same column, different letters indicate differences between means (P < 0.05). (m): Ulmus minor; (mp): U. minor × U. pumila.
Parent | Families with a female parent | Families with a male parent | ||
---|---|---|---|---|
Annual height growth (cm) |
Leaf wilting (%) |
Annual height growth (cm) |
Leaf wilting (%) |
|
AL-AN1 (m) | 40.7 a | 35.9 bc | 42.2 abcd | 12.7 abc |
AL-JR1 (m) | 49.5 b | 37.9 bc | 37.3 a | 16.3 abc |
GR-AL3 (m) | 47.8 b | 19.1 a | 37.3 a | 29.7 abc |
M-PZ3 (mp) | 51.3 bc | 22.0 a | 48.1 bc | 20.2 bc |
SE-CT4 (m) | 36.6 a | 37.9 c | 49.5 bc | 25.0 a |
VA-VV5 (mp) | 57.2 c | 19.7 a | 47.4 bc | 21.5 abc |
TO-AL1 (m) | 38.9 a | 27.7 ab | 39.4 a | 26.4 c |
CC-VG2 (m) | - | - | 58.9 d | 25.7 abc |
M-IN5 (m) | - | - | 49.9 bc | 18.7 a |
V-JR1 (m) | - | - | 53.2 cd | 17.8 ab |
Ulmus minor ramets showed earlier mean vegetative budburst, lower mean height growth rates and higher mean wilting values than U. minor × U. pumila ramets (P < 0.05 - Tab. 4). Differences in traits between U. minor and U. minor × U. pumila open pollinated progenies were statistically significant only for leaf wilting, and U. minor × U. pumila was more resistant than U. minor (P < 0.05 - Tab. 4). A notable result of controlled pollinated seedlings was the leaf wilting values obtained depending on whether the male parent was U. minor or U. minor × U. pumila (21.5 and 30.6%, respectively). Resistance was higher when U. minor was used as male parent (P < 0.05).
Tab. 4 - Mean values of five traits measured on ramets and open pollinated progenies (op). Different letters indicate differences between means (P < 0.05). Values within parenthesis correspond to the range. (a): trees with beetle feeding grooves on trunk or branches; (b) assessed on April 9th 2005 on a scale from 1= dormant buds to 5= two or more leaves completely expanded; (c) assessed on July 7th 2005, 60 days after inoculations with Ophiostoma novo-ulmi; (m): Ulmus minor; (mp): U. minor × U. pumila.
Plant materiala | No. trees | Flowering trees (%) |
Beetle affected trees (%)b |
Vegetative budburstc |
Annual height growth (cm) |
Leaf wilting (%)d |
|||
---|---|---|---|---|---|---|---|---|---|
m ramets | 46 | 2.2 | 4.3 | 4.1 bc | (1-5) | 27.7 a | (6-68) | 46.4 d | (5-100) |
mp ramets | 12 | 75.0 | 75.0 | 2.9 a | (2-4) | 73.3 d | (36-102) | 28.8 ab | (0-70) |
m op | 52 | 5.7 | 17.3 | 4.5 bc | (2-5) | 45.3 bc | (7-95) | 38.0 cd | (0-100) |
mp op | 26 | 3.8 | 3.8 | 4.2 bc | (1-5) | 50.5 bc | (15-80) | 22.9 ab | (0-60) |
Heritability estimates from the regression (Tab. 5) and analysis of variance (Tab. 6) methods ranged from 0.13 to 0.69 and 0.06 to 0.71, respectively. Disease resistance heritability values were higher when estimated 60 dpi than 15, 30 or 120 dpi (results not shown). Expected genetic gains for height by selection in progenies from the M-PZ3 × TO-AL1 and VA-VV5 × CC-VG2 cross-combinations would be 3.07 and 3.07 cm, respectively. Expected genetic gains to reduce leaf wilting by selection in progenies from the M-PZ3 × SE-CT4, VA-VV5 × M-IN5 and VA-VV5 × V-JR1 cross-combinations would be 7.16, 9.46 and 7.83%, respectively.
Tab. 5 - Heritability estimates (± SE) of Ulmus minor and U. pumila hybrids for vegetative budburst, height growth, and resistance to Ophiostoma novo-ulmi obtained from regression between ramets and open pollinated families, and between ramets and controlled pollinated families. (a): average values of two x and y clones were correlated to their respective x × y families.
Plant material | Vegetative budburst |
Height growth |
Disease resistance |
---|---|---|---|
Ramets vs open pollinated families | 0.38 ± 0.34 | 0.41 ± 0.32 | 0.13 ± 0.32 |
Ramets vs controlled pollinated families a | 0.64 ± 0.18 | 0.69 ± 0.17 | 0.36 ± 0.35 |
Tab. 6 - Variance components and heritability estimates (± SE) in the broad (H2) and narrow (h2) sense of Ulmus minor and U. pumila hybrids for vegetative budburst, height growth, and resistance to Ophiostoma novo-ulmi. (a): the average number of seedlings for each family was assumed to be r = 8.
Variance components and heritability estimates |
Vegetative budburst |
Height growth |
Disease resistance |
---|---|---|---|
σ 2g ± SE |
0.0602 ± 0.0406 | 48.2 ± 26.7 | 81.3 ± 25.1 |
σ 2 e |
0.9030 | 456.2 | 270.2 |
H 2 single-plant basis ± SE |
0.06 ± 0.04 | 0.10 ± 0.05 | 0.23 ± 0.07 |
H 2 family mean basis ± SE |
0.35 ± 0.23 | 0.46 ± 0.25 | 0.71 ± 0.22 |
σ 2males ± SE |
0.0366 ± 0.011 | 0.0 ± 3.8 | 9.1 ± 7.1 |
σ 2 females |
0.0064 | 32.0 | 23.3 |
σ 2 males × females |
0.0669 | 48.7 | 0.0 |
σ 2 add |
0.0860 | 64.0 | 64.8 |
σ 2 dom |
0.2676 | 194.8 | 0.0 |
σ 2 e |
0.3504 | 47.5 | 412.1 |
h 2 single-plant basis ± SE |
0.12 ± 0.06 | 0.21 ± 0.05 | 0.14 ± 0.06 |
h 2 family mean basis ± SE |
0.22 ± 0.11 | 0.24 ± 0.06 | 0.56 ± 0.24 |
Annual height was phenotypically related to vegetative budburst, as budburst occurred earlier in seedlings with low rather than high growth rates (P = 0.0007 - Tab. 7). Wilting percentage was negatively related to early budburst (P < 0.0001). Other significant relations among the controlled pollinated seedlings included percentage of flowering trees and annual height growth (r
p = 0.44; P = 0.0042), percentage of flowering trees and vegetative budburst (r
p = -0.53; P = 0.0004), and percentage of beetle affected-trees and annual height growth (r
p = 0.60; P <0.0001).
Tab. 7 - Phenotypic correlations (above the diagonal) between traits assessed in Ulmus minor and U. minor × U. pumila controlled pollinated seedlings. (a): indicates significance at P < 0.001.
- | Annual height growth |
Vegetative budburst |
Leaf wilting |
---|---|---|---|
Annual height growth |
x | -0.26 a | 0.06 |
Vegetative budburst |
- | x | -0.17 a |
Leaf wilting |
- | - | x |
Discussion
To our knowledge, this is the first study to report heritability values and non-additive genetic variances of vegetative budburst, height growth and resistance to O. novo-ulmi in a breeding population of elms. Although the breeding of DED resistant elms began in the Netherlands in 1928 and the United States in 1933 ([27]), little is known about inheritance of DED resistance in elms ([41]). Genetically based resistance of trees to other forest diseases, such as poplar leaf rust ([13]), birch rust ([9]), Armillaria root rot ([45]) and pitch canker ([43]), suggests that breeding trees for resistance is possible. In the Ulmus genus, the response of backcross hybrids of U. pumila, U. japonica and U. rubra hybrids to inoculation with O. novo-ulmi has been reported ([14], [15], [16]), although heritability values were not quantified. Within a U. minor population, narrow-sense heritability for wilting percentages was 0.54 ([41]). Heritability values reported here confirm the hypothesis that vegetative budburst, height growth and DED resistance are genetically controlled.
Our first assessments of height growth, vegetative budburst and disease resistance in the controlled pollinated progenies were based on plant means. The primary functions of heritability estimates on a single-plant basis are to predict genetic progress through phenotypic selection and provide the relative magnitude of the genetic basis of height growth, vegetative budburst and disease resistance. However, the selection of families is probably based on family means and therefore the appropriate heritability for predicting genetic gain is the ratio of the variance component due to families to total phenotypic variance among family means, i.e., family mean heritability.
The estimates of heritability for resistance reported here are probably conservative due to the moderate levels of wilting obtained. The cell density used in the inoculum, lower than the density generally used in resistance tests for breeding purposes (= 106 conidia ml-1 - [33], [36]), was chosen to avoid high mortality. If we had observed a wider range of wilt symptoms, the genetic variability would have been readily expressed and the heritability values would probably have been higher ([39], [12]). We acknowledge that the number of replicates should have been higher, but unfortunately some clones used for the diallel design, coded CC-VG2, M-IN5 and V-JR1 and later characterized as female-sterile trees ([17], [18]), provided no seeds.
Ulmus minor flushes earlier in spring than other European elms ([31]) and U. pumila flushes even earlier, at least in the conditions in Madrid. However, vegetative budburst occurred later in U. minor × U. pumila than in U. minor ramets. This finding must be considered with caution, as the U. minor × U. pumila sample was constituted by just one clone (M-PZ3). The late vegetative budburst in U. minor × U. pumila hybrids could be related to the fact that hybrid ramets have greater fruit production than U. minor ramets. In the clonal bank at Puerta de Hierro Forest Breeding Center, it was observed that vegetative budburst is generally delayed if elms have abundant samara production, probably because of the energy cost required for fruit production. The significant relation observed between the percentage of flowering and late vegetative budburst could support this hypothesis. This relation would prevent proper comparison of budburst values between elms at a different juvenile or adult stage.
The negative association between vegetative budburst date and plant height, confirming earlier studies ([16]), should be considered with caution. Low trees, with buds close to soil level, would be exposed to higher temperatures during the day due to soil irradiation and would therefore flush earlier than taller trees ([46]). Considering only trees taller than 1.5 m, unaffected by soil irradiation, the phenotypic relation between vegetative budburst and plant height is not significant (r
p = 0.09; P > 0.1). Similar findings on the phenotypic relations between vegetative budburst and susceptibility to DED have been reported ([16], [7]). This relation seems to be independent of tree height, as the same calculation made for trees above 1.5 m still shows significance (r
p = -0.15; P = 0.0061). Elms with early vegetative budburst could be more resistant because they may develop latewood earlier than elms with delayed budburst. Early latewood formation has been related to DED resistance ([30], [34]), as small diameter vessels and pits are thought to compartmentalize the disease better and be less vulnerable to cavitation than large diameter vessels and pits ([34], [21], [24]).
The high frequency of beetle feeding grooves in elms with U. minor × U. pumila as a parent may be related to the higher growth rates found in these individuals and not to their greater attractiveness to beetles. It would be interesting to test this hypothesis when all trees have reached larger size and maturity.
For most purposes, fast growth of elms is highly desirable ([11]). However, one difficulty of this is that fast-growing trees often show high susceptibility to O. novo-ulmi ([16], [37]), although this relation is not universal ([11], [36]), as confirmed here. Annual growth quantification of seedlings in the juvenile stage is sometimes erratic. Moderate growth in youth may be followed by subsequent fast growth, whereas fast growth in youth may later decline. Long-term field trials would be required to properly quantify this trait. Breeders should take advantage of height growth heterosis and select the most appropriate combination. Heterosis was observed in elm growth in crosses where U. glabra and U. wallichiana were used as parent trees ([11]), as height was 50% greater in progenies than in parents.
Response to DED varied significantly depending on whether U. minor × U. pumila was used as female or male parent. The higher resistance of (U. minor × U. pumila) × U. minor hybrids may result from using U. minor × U. pumila as female parent. This combination could contain more germplasm from U. pumila than U. minor × (U. minor × U. pumila), although this statement needs experimental support. The result obtained is of practical interest in breeding programmes because it provides information about the optimal choice for obtaining resistant progenies. As crosses are time-consuming, choosing the option of (U. minor × U. pumila) × U. minor rather than U. minor × (U. minor × U. pumila) would save both time and effort. Breeding strategies should therefore concentrate on combinations such as M-PZ3 × SE-CT4, VA-VV5 × TO-AL1, VA-VV5 × M-IN5 and VA-VV5 × V-JR1. The high genetic variation of height growth and disease response within some (U. minor × U. pumila) × U. minor combinations (e.g., VA-VV5 × TO-AL1) may be of great value in advanced stages of breeding. The timing of vegetative budburst should also be taken into account: in sites with a severe continental climate or a large population of defoliating insects, care must be taken with combinations that result in seedlings with early budburst.
In this study, it was possible to obtain satisfactory gains for the three traits assessed. Estimated gains for growth and tolerance to DED were promising and, as far as we know, they are probably the first to be reported in an elm population. Further re-selection of other elm trees for height and disease resistance will improve the existing population, allowing a considerable increase in the predicted gains.
Conclusions
Evidence on the effect of genes from U. pumila in terms of resistance to O. novo-ulmi is limited but encouraging. This species has frequently been used in European breeding programmes because of its greater DED resistance than U. minor. This study confirms the suitability of this practice. Resistance to O. novo-ulmi increased significantly as a function of increased amounts of U. pumila germplasm from the female parent, suggesting that DED resistance may be transmitted mainly by the mother. In (U. minor × U. pumila) × U. minor combinations, with U. pumila contributing approximately 25% to the germplasm, response to DED was sufficient to recommend this backcross for breeding purposes. DED resistance heritability estimates and genetic gains reported here indicate a high degree of additive genetic control and show that selection for DED resistance and rapid tree growth is effective.
Acknowledgments
This work is dedicated to Hans Heybroek, Lorenzo Mittempergher, Alberto Fagnani, Fabio Ferrini and Alberto Santini, who kindly shared their knowledge of elm breeding with AS, JCLA and LG. We thank Margarita Burón† for superb technical assistance, Christine Cahalan (University of Wales, Bangor) for help with the statistical analysis, Esther Merlo for a critical review, and Jane McGrath for language editing. Pictures (a) and (b) from Fig. 2 were taken by S. Dominguez and E. Martínez. AS was funded by the Regional Government of Extremadura Department of Education and Technology (Junta de Extremadura - Consejería de Educación y Tecnología) and the European Social Fund during a visit to C. Cahalan in Bangor. The study was funded by the Spanish Directorate-General of Rural Development and Forestry Policy, Ministry of Agriculture, Food and Environmental Affairs (Dirección General de Desarrollo Rural y Política Forestal, Ministerio de Agricultura, Alimentación y Medio Ambiente).
References
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Gscholar
Supplementary Material
Authors’ Info
Authors’ Affiliation
Ingeniería Forestal y del Medio Natural, Universidad de Extremadura, Avenida Virgen del Puerto 2, E-10600 Plasencia (Spain)
Departamento de Desarrollo Sostenible, Universidad Católica de Ávila, Canteros s/n, E-05005 Ávila (Spain)
Luis Gil
ETSI Montes, Forestal y del Medio Natural, Universidad Politécnica de Madrid, Ciudad Universitaria s/n, E-28040 Madrid (Spain)
Corresponding author
Paper Info
Citation
Solla A, López-Almansa JC, Martín JA, Gil L (2014). Genetic variation and heritability estimates of Ulmus minor and Ulmus pumila hybrids for budburst, growth and tolerance to Ophiostoma novo-ulmi. iForest 8: 422-430. - doi: 10.3832/ifor1227-007
Academic Editor
Alberto Santini
Paper history
Received: Dec 31, 2013
Accepted: Sep 23, 2014
First online: Dec 15, 2014
Publication Date: Aug 02, 2015
Publication Time: 2.77 months
Copyright Information
© SISEF - The Italian Society of Silviculture and Forest Ecology 2014
Open Access
This article is distributed under the terms of the Creative Commons Attribution-Non Commercial 4.0 International (https://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Web Metrics
Breakdown by View Type
Article Usage
Total Article Views: 26318
(from publication date up to now)
Breakdown by View Type
HTML Page Views: 19886
Abstract Page Views: 1224
PDF Downloads: 3757
Citation/Reference Downloads: 50
XML Downloads: 1401
Web Metrics
Days since publication: 3593
Overall contacts: 26318
Avg. contacts per week: 51.27
Article Citations
Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Nov 2020)
Total number of cites (since 2015): 9
Average cites per year: 1.50
Publication Metrics
by Dimensions ©
Articles citing this article
List of the papers citing this article based on CrossRef Cited-by.
Related Contents
iForest Similar Articles
Research Articles
Age trends in genetic parameters for growth and quality traits in Abies alba
vol. 9, pp. 954-959 (online: 07 July 2016)
Research Articles
Patterns of genetic variation in bud flushing of Abies alba populations
vol. 11, pp. 284-290 (online: 13 April 2018)
Research Articles
Seedling emergence capacity and morphological traits are under strong genetic control in the resin tree Pinus oocarpa
vol. 17, pp. 245-251 (online: 16 August 2024)
Research Articles
Comparison of genetic parameters between optimal and marginal populations of oriental sweet gum on adaptive traits
vol. 11, pp. 510-516 (online: 18 July 2018)
Research Articles
Genetic diversity of core vs. peripheral Norway spruce native populations at a local scale in Slovenia
vol. 11, pp. 104-110 (online: 31 January 2018)
Research Articles
Genetic variation of Fraxinus excelsior half-sib families in response to ash dieback disease following simulated spring frost and summer drought treatments
vol. 9, pp. 12-22 (online: 08 September 2015)
Review Papers
Genetic diversity and forest reproductive material - from seed source selection to planting
vol. 9, pp. 801-812 (online: 13 June 2016)
Commentaries & Perspectives
The genetic consequences of habitat fragmentation: the case of forests
vol. 2, pp. 75-76 (online: 10 June 2009)
Research Articles
Fine-scale spatial genetic structure in a multi-oak-species (Quercus spp.) forest
vol. 8, pp. 324-332 (online: 05 September 2014)
Research Articles
Genetic control of intra-annual height growth in 6-year-old Norway spruce progenies in Latvia
vol. 12, pp. 214-219 (online: 25 April 2019)
iForest Database Search
Search By Author
Search By Keyword
Google Scholar Search
Citing Articles
Search By Author
Search By Keywords
PubMed Search
Search By Author
Search By Keyword