Abstract
Plant architecture is a crucial trait in plant breeding because it is linked to crop yield. For ornamental crops such as roses, plant architecture is key for their aesthetic and economic value. In 2015, six rose plant architectural traits were evaluated on 2- to 3-year-old plants of F1 rose populations in May and December in College Station, TX, to estimate variability and heritability. The traits included plant height, the number of primary shoots, the length of primary shoots, the number of nodes on the primary shoot, the number of secondary shoots per primary shoot, and the number of tertiary shoots per primary shoot. Among these traits, plant height, the number of primary shoots, and the length of primary shoots showed a substantial amount of variability, whereas the number of secondary and tertiary shoots per primary shoot were skewed toward zero. Using a random effects model restricted maximum likelihood (REML) analysis, the architectural traits demonstrated low to moderate narrow-sense heritability (0.12–0.50) and low to high broad-sense heritability (0.34–0.92). Plant height and the number of primary shoots changed little after the first growth phase, whereas the other four traits were affected greatly by the genotype-by-environment (growth phase) interaction.
Roses, which are admired for their diversity of floral and vegetative characteristics, are members of Rosaceae, the most important horticultural family in the world. They are the world’s most important ornamental crop, with an estimated annual production of 18 billion cut stems, 60 to 80 million potted roses, and 220 million roses for landscape uses (Blom and Tsujita, 2003; Pemberton et al., 2003; Roberts et al., 2003). In the United States, 36.6 million garden rose bushes are produced, generating sales of $203.5 million (U.S. Department of Agriculture, National Agricultural Statistics Service, 2015)
There are four subgenera under the genus Rosa, which has more than 100 species, ranging from diploid to decaploid. There are more than 30,000 cultivars (mainly diploid, triploid, and tetraploid), with a wide interspecific and intraspecific cross-compatibility (Blechert and Debener, 2005; Byrne and Crane, 2003; Cairns, 2000; Jian et al., 2010; Ueckert et al., 2015; Zlesak, 2009). The role roses play in landscapes is extensive, with their use adorning roadsides, public parks, commercial properties, and residential areas. Garden roses provide aesthetic value throughout the growing cycle for both their vegetative and reproductive organs (De Vries and Dubois, 1996).
For ornamental crops like roses, plant architecture is the key factor determining the appearance of the plant and its commercial value. In recent years, researchers have constructed models that describe the plant architecture of rose bushes according to growth and branching processes. These traits can be characterized morphologically by their number, length, and diameter; by the way they are connected; and by their positioning in space (i.e., shoot angles).
Roses have a wide diversity of plant architectures, ranging from prostrate to compact, highly branched forms to large, spreading to erect bushes to those that climb. Among garden roses, many of these growth forms are useful for specific garden situations, whereas for cut-flower roses, a growth form that produces long stems (50 cm or greater), good basal breaking, and little side shoot formation is required.
The Texas A&M University Rose Breeding and Genetics program has the objective of developing well-adapted (resistant to black spot, Cercospora, and rose rosette disease; tolerant to extreme summer heat), compact, highly floriferous bushes suitable for garden use. Thus, the program is very interested in how a rosebush grows and produces flowers not just during its first flush and growth phase during the spring and early-summer months, but also during the second growth phase in mid- to late-summer and fall months.
The broad-sense heritability of rose architectural traits has been estimated using a set of eight diploid rose cultivars grown in pots in a greenhouse (Crespel et al., 2013, 2014), a small (98 individuals) F1 diploid population planted in the field (Kawamura et al., 2011, 2015), and a biparental F1 tetraploid cut flower population that was planted under greenhouse conditions in the Netherlands and Kenya (Gitonga et al., 2014). In general, researchers found that the architectural traits they measured during the first flush of growth after dormancy (garden-rose situation) or pruning (cut-rose situation) showed a moderate to high broad-sense heritability. Important genotype-by-environment and genotype-by-year interactions were also identified and attributed to differences in total radiation and temperatures that affected total growth rate and budbreak along the stem (Crespel et al., 2014; Gitonga et al., 2014). Kawamura et al. (2011) identified various quantitative trait loci (QTLs) associated with the flower inflorescence structure and identified candidate genes involved in auxin and gibberellin metabolism that controlled node production, internode elongation, and axillary bud branching. Unfortunately, none of these studies analyzed the second growth phase of a garden rose that occurs after the first flush. In Texas, this second growth phase begins in mid to late summer and continues until fall temperatures put the plant into dormancy.
A garden rose needs to flower well not only during this first flush after dormancy, but also throughout the summer and fall months. As the Texas A&M Rose Breeding and Genetics program has the objective of developing roses that bloom continuously throughout the growing cycle and that establish a compact growth form, our experiment was designed to compare the growth of a set of breeding populations at first flush and later at the end of the year to understand the architectural components important in flower productivity and bush shape. The objective of this research was to evaluate the genetic variation and heritability of six plant architectural components within 13 diploid rose populations.
Materials and Methods
Plant materials.
Thirteen diploid rose populations (Table 1) developed from the hybridization of rose parents with diverse plant architecture including groundcovers derived from crosses with Rosa wichuriana ‘Basye’s Thornless’, the old China roses ‘Old Blush’ and ‘Ducher’, and various miniature rose cultivars were used for this study. These populations represent a multispecies breeding background typical in rose breeding. In this case, the species involved are primarily Rosa wichuriana and Rosa chinensis, which are from the related and generally interfertile subgenera Systylae and Chinesis. We evaluated up to 40 seedlings per family; for those less than 40, we evaluated all seedlings.
Thirteen diploid rose families evaluated for six architectural traits in May and December in the field in 2015 at College Station, TX.
Field experiment and trait assessment.
The seedlings of 13 diploid populations (partial diallel) in this study were planted at a 1 × 1 × 3.5-m spacing in a field 2.5 km from Texas A&M University in College Station (lat. 30.65N, long. –96.32W). Rose plants were planted in rows oriented east to west in an open field. Raised beds were constructed in the field and black cloth weed barrier was placed around the rosebushes for weed control. Irrigation was applied as needed. No pesticides were applied during the evaluation. In March, the rose plants were pruned to a standard size by removing 50% to 75% of the previous year’s growth. Dead growth was also removed. Each genotype was represented by one plant in the field, and the seedlings from each family were grouped together. Evaluation of the mature (2–3 years old) roses for plant architecture was done in the field twice in 2015: spring (May) after the first flower flush; and winter (December), at the end of the growing cycle. Six architectural traits were chosen to represent distinct aspects of variables describing plant architecture to explain the architectural variability observed and to combine correlated variables into one (Crespel et al., 2013). They were plant height (measured in centimeters), the number of primary shoots (initial complete shoots that develop from the base of the plant), the length of the primary shoot (measured in centimeters), the number of nodes on the primary shoot, the number of secondary shoots per primary shoot, and the number of tertiary shoots per primary shoot. Measurements were made on three basal shoots per rose seedling.
Statistical analysis.
The statistical analysis used JMP software, version 12. Six architectural traits in 13 diploid rose populations were subjected to normality analysis (original and transformed by taking natural logarithm, base 10 logarithm, and square root) using the Shapiro-Wilk W test. Analysis of variance (ANOVA) was conducted to compare family, growth phase means, and family-by-growth phase interaction. Correlation analysis was conducted among the six traits.
The variance components were calculated with the REML method, assuming all factors as random effects for a more powerful estimation (Littell et al., 1996). The model for heritability analysis was Y = μ + Female + Male + Progeny (female, male) + Growth Phase + Growth phase × Female + Growth phase × Male + Growth phase × Progeny (female, male). Because we only have one rose seedling per genotype, residual variance was confounded with Growth phase × Progeny (female, male) variance. The variances of the parents were considered as additive variance (VA), the progeny variance was considered as nonadditive variance (VD), repeated measurement between May and December data variance was considered as variance of the environment (VE), and the sum of growth phase interactions (VG×E) was also estimated, where G is genotype and E is environment. Narrow- (h2) and broad-sense (H2) heritability was estimated by the genetic variance from the REML model, where VP = (VA + VD + VG×E/e), where e is the number of environments (first, second growth phase), h2 = VA/VP, and H2 = (VA+VD)/VP (Connor et al., 2005; Hallauer et al., 2010; Holland et al., 2003).
To quantify the variation resulting from the genotype-by-environment interaction relative to main genotypic variation, the architectural traits over growth phases were analyzed with a ratio of genotype-by-environmental variance to genetic variance (Gitonga et al., 2014). Three traits that have a relatively high VG×E/VG (the number of nodes on the primary shoot, the number of secondary shoots per primary shoot, and the number of tertiary shoots per primary shoot) were subjected to two-way ANOVA to compare means of each trait measured in May and Dec. 2015 of 13 rose families.
Results
Distribution of six architectural traits.
All statistical analyses were based on the original data because the data were either normally distributed (plant height and the length of the primary shoot), or data transformation (log, log10, and square root) did not improve data normality.
A substantial amount of variability was associated with plant height (range, 14–64 cm), the number of primary shoots (range, 3–53), the length of the primary shoot (range, 14–64 cm), and the number of nodes on the primary shoot (range, 3–18). The number of secondary shoots per primary shoot (range, 0–10) and the number of tertiary shoots per primary shoot (range, 0–15) were skewed toward zero (Table 2). The traits that exhibited substantial variability should allow for genetic improvement.
Mean separations, comparison of growth phases, and genotype-by-environment interaction of six architectural traits for 13 diploid rose families evaluated in the field in May and Dec. 2015 in College Station, TX.
There were differences among the families in all the architectural traits measured as well as differences among the data taken in May vs. December in the number of nodes on the primary shoot, the number of secondary shoots per primary shoot, and the number of tertiary shoots per primary shoot. Although plant height, the number of primary shoots, and the length of the primary shoot were predictable with first growth phase information, the growth of the number of nodes for primary/secondary/tertiary shoots varied among the populations during the growth phases (Tables 2 and 3). Most families, as expected, showed an increased number of secondary/tertiary shoots per primary shoot (11 of 13 and 12 of 13, respectively) in December compared with May (Table 4).
Two-way analysis of variance comparing means of the length of the primary shoot and the number of nodes on the primary shoot for 13 diploid rose families characterized in May and Dec. 2015 in College Station, TX.
Two-way analysis of variance comparing means of the number of secondary and tertiary shoots per primary shoot for 13 diploid rose families characterized in May and Dec. 2015 in College Station, TX.
Correlation analysis.
Among the six architectural components, plant height and the number of primary shoots exhibited low correlation (R = 0.39). Components that were moderately correlated include plant height and the length of the primary shoot (R = 0.52), the length of the primary shoot and the number of nodes on the primary shoot (R = 0.59), the number of nodes on the primary shoot and the number of the secondary shoots per primary shoot (R = 0.48), and the number of secondary shoots per primary shoot and the number of tertiary shoots per primary shoot (R = 0.62). R2 values indicate these correlations could explain between 20% to 40% of the variation, and therefore, are not high enough to use for selection purposes.
Heritability analysis.
The six architectural traits measured exhibited low to high broad-sense heritability (0.34–0.92) and low to moderate narrow-sense heritability (0.12–0.50) (Table 5), which indicated that architectural traits are feasible targets of rose breeding.
Variance component, broad-sense heritability (H2), and narrow-sense heritability (h2) for six plant architectural traits evaluated in the field on 13 diploid rose families in May and Dec. 2015 in College Station, TX.
Plant height was found highly heritable with a moderate narrow-sense heritability (h2 = 0.50) and high broad-sense heritability (H2 = 0.82) (Table 5). This finding was consistent with estimates reported previously by Gitonga et al. (2014) (H2 = 0.82), who measured this architectural component on potted tetraploid cut-flower rose populations, and by Kawamura et al. (2015) (H2 = 0.88), who studied the architecture of a small (98 individuals) diploid garden-rose population. Both studies reported broad-sense heritability only and have smaller populations and fewer families than our study.
The number of primary shoots had a low narrow-sense heritability (h2 = 0.27) and high broad-sense heritability (H2 = 0.92), indicating a strong nonadditive genetic component. This accounted for ≈60% of total phenotypic variance and was more than twice that of the narrow-sense estimate for this trait (Table 5). Previous studies evaluated the number of long axes (shoots with five or more metamers each which consist of an internode, a node, an axillary bud and a leaf) and the number of determined axes (axes ended in a flower bud or a flower) (Crespel et al., 2013, 2014)—which are equivalent to some extent to the number of primary shoots in our study. The broad-sense heritability for the number of long axes and the number of determined axes assessed using eight diploid garden-rose cultivars in the greenhouse by Crespel et al. (2014) were 0.70 and 0.64, respectively, and thus were less than the number of primary shoots (H2 = 0.92) (Table 5) in our study using segregating populations of diploid garden roses.
The length of the primary shoot showed low narrow-sense (h2 = 0.20) and a moderately high broad-sense (H2 = 0.64) heritability (Table 5). This is greater than a previous report by Crespel et al. (2014) on the length of the long axes (H2 = 0.48). Nonadditive genetic variance components accounted for ≈32% of phenotypic variance compared with 15% (Table 5) for the additive genetic component.
The number of nodes on the primary shoot revealed low narrow-sense (h2 = 0.12) and moderate broad-sense (H2 = 0.46) heritability (Table 5). Crespel et al. (2014) reported that the broad-sense heritability of the number of metamers (an internode, a node, an axillary bud, and a leaf) of the long axes (equivalent to the number of nodes on primary shoots) was 0.95, which is double that found in our study (H2 = 0.46) (Table 5). This is probably, in part, a result of the inclusion of extreme growth types in their study that would inflate the total variation compared with the segregating population in our study.
The number of secondary and tertiary shoots per primary shoot had low narrow-sense (h2 = 0.16 and 0.21, respectively) and low to moderate broad-sense heritability (H2 = 0.34 and 0.48, respectively) (Table 5).
Genotype-by-environment interaction.
Ornamental plants are expected to be attractive throughout the year. With a rosebush, this requires that the bush remains well leafed, compact, and continues to flower throughout the annual growing cycle. Thus, in our analysis, the different environments refer to the two phases of rose growth during the growing cycle. First phase is the first flush of a rose plant during which it emerges from dormancy and develops the bush shape and the first flush of flowers. In Texas, the second phase of growth starts after the first flush of flowers in early summer and continues until fall temperatures stop growth. Among six architectural traits, plant height and the number of primary shoots showed a weak genotype-by-environment interaction, as indicated by the low VG×E to VG ratios (0.43 and 0.18, respectively) (Table 5), reflecting consistency across growth phases. The length of the primary shoot was moderately affected by the genotype-by-environment interaction (VG×E/VG ratios of 1.10) (Table 5). These first three traits changed little after the first flush among these remontant rose populations, indicating that these traits have little influence over the flowering behavior of the rose during the second growth phase.
Strong genotype-by-environment interactions were found for the number of nodes on the primary shoot, and the number of secondary shoots and tertiary shoots per primary shoot (VG×E/VG ratios of 2.3, 3.9, and 2.1, respectively) (Table 5). Previously, ANOVA (Tables 2–4) indicated that, after the initial growth of the primary shoots in the spring, additional growth continued mainly by the development of secondary and tertiary shoots during the second growth phase of the rosebush. These genotype-by-environment interactions indicate a differential response of genotypes to the environment in which they are grown (Allard and Bradshaw, 1964; Bernardo, 2010). For the length of the primary shoot and the number of nodes on the primary shoot, the genotype-by-environment interaction was mainly the result of the behavior of the families with ‘Vineyard Song’ as a parent. These grew longer and put on more nodes during the May-to-December time period than the other families. The number of secondary and tertiary shoots formed during this second growth phase varied among the roses, and thus so did their ability to produce new foliage and flowers—an essential requirement to keep the rosebush attractive throughout the growing cycle.
Discussion
The six architectural traits measured exhibited low to high broad-sense heritability (0.34–0.92) and low to moderate narrow-sense heritability (0.12–0.50). Broad-sense heritability has been used widely for estimating the accuracy of selection for genetic traits (Holland et al., 2003) and QTL analysis (Beavis, 1998). Because rose is a vegetatively propagated crop, the nonadditive genetic component of the variation can be captured easily by clonal propagation. The following conclusions can be drawn from this study.
Our results, combined with previous studies (Crespel et al., 2014; Kawamura et al., 2015) that reported high broad-sense heritability for several architectural traits of garden roses, indicated that rose plant architecture was a feasible target for rose breeding. Architectural traits with low narrow-sense heritability but moderately high to high broad-sense heritability, such as the number of primary shoots, the length of the primary shoot, and the number of nodes on the primary shoot, indicated important nonadditive effects.
Architectural traits were affected by the genotype-by-environment interaction to different extents. Plant height and the number of primary shoots were highly heritable and changed little during the second growth phase. These can be selected for reliably and accurately during the first growth phase. The number of secondary and tertiary shoots per primary shoot showed lower heritability and tended to increase in number as the year progressed. The extent of this side shoot growth varied by population (high genotype-by-environment interaction). This indicated that selection during the first growth phase would not predict the final number of secondary/tertiary shoots forming throughout the year. Therefore, selection during both growth phases was recommended. Finally, for the selection of roses that flower continuously throughout the year, it is essential to select for high production of basal shoots during the first growth phase and high production of secondary and tertiary shoots during both growth phases.
Literature Cited
Allard, R.E. & Bradshaw, R.W. 1964 Implications of genotype–environment interactions in applied plants breeding Crop Sci. 4 503 508
Beavis, W.D. 1998 QTL analyses: Power, precision, and accuracy, p. 145–162. In A.H. Paterson (ed.). Molecular dissection of complex traits. CRC Press, New York, NY
Bernardo, R. 2010 Breeding for quantitative traits in plants. 2nd ed. Stemma Press, Woodbury, MN
Blechert, O. & Debener, T. 2005 Morphological characterization of the interaction between Diplocarpon rosae and various rose species Plant Pathol. 54 82 90
Blom, T.J. & Tsujita, M.J. 2003 Cut rose production, p. 594–600. In: A.V. Roberts, T. Debener, and S. Gudin (eds.). Encyclopedia of rose science. Elsevier, Oxford, UK
Byrne, D.H. & Crane, Y.M. 2003 Meiosis, p. 273–279. In: A.V. Roberts, T. Debener, and S. Gudin (eds.). Encyclopedia of rose science. Elsevier, Oxford, UK
Cairns, T. 2000 Modern roses XI: The world encyclopedia of roses. Academic Press, San Diego, CA
Connor, A.M., Stephens, M.J., Hall, H.K. & Alspach, P.A. 2005 Variation and heritabilities of antioxidant activity and total phenolic content estimated from a red raspberry factorial experiment J. Amer. Soc. Hort. Sci. 130 403 411
Crespel, L., Le Bras, C., Relion, D. & Morel, P. 2014 Genotype × year interaction and broad-sense heritability of architectural characteristics in rose bush Plant Breed. 133 3 412 418
Crespel, L., Sigogne, M., Donès, N., Relion, D. & Morel, P. 2013 Identification of relevant morphological, topological and geometrical variables to characterize the architecture of rose bushes in relation to plant shape Euphytica 191 1 129 140
De Vries, D.P. & Dubois, L.A.M. 1996 Rose breeding: Past, present, prospects Acta Hort. 424 241 248
Gitonga, V.W., Koning-Boucoiran, C.F.S., Verlinden, K., Dolstra, O., Visser, R.G.F., Maliepaard, C. & Krens, F.A. 2014 Genetic variation, heritability and genotype by environment interaction of morphological traits in a tetraploid rose population BMC Genet. 15 1 14
Hallauer, A.R., Carena, M.J. & Miranda Filho, J.B. 2010 Quantitative genetics in maize breeding. Springer, New York, NY
Holland, J.B., Nyquist, W.E. & Cervantes-Martínez, C.T. 2003 Estimating and interpreting heritability for plant breeding: An update Plant Breed. Rev. 22 9 112
Jian, H., Zhang, H., Tang, K., Li, S., Wang, Q., Zhang, T., Qiu, X. & Yan, H. 2010 Decaploidy in Rosa praelucens Byhouwer (Rosaceae) endemic to Zhongdian plateau, Yunnan, China Caryologica 63 162 167
Kawamura, K., Hibrand-Saint Oyant, L., Crespel, L., Thouroude, T., Lalanne, D. & Foucher, F. 2011 Quantitative trait loci for flowering time and inflorescence architecture in rose Theor. Appl. Genet. 122 661 675
Kawamura, K., Hibrand-Saint Oyant, L., Thouroude, T., Jeauffre, J. & Foucher, F. 2015 Inheritance of garden rose architecture and its association with flowering behaviour Tree Genet. Genomes 11 2 1 12
Littell, R.C., Milliken, G.A., Stroup, W.W. & Wolfinger, R.D. 1996 SAS system for mixed models. SAS Institute, Inc., Cary, NC
Pemberton, H.B., Kelly, J.W. & Ferare, J. 2003 Pot rose production, p. 587–593. In: A.V. Roberts, T. Debener, and S. Gudin (eds.). Encyclopedia of rose science. Elsevier, Amsterdam, Netherlands
Roberts, A.V., Debener, T. & Gudin, S. 2003 Introduction, p. vi–vii. In: A.V. Roberts, T. Debener, and S. Gudin (eds.). Encyclopedia of rose science. Elsevier, Amsterdam, Netherlands
Ueckert, J., Byrne, D.H., Crosby, K., Hodnett, G. & Stelly, D. 2015 The utilization of the polyploid nature of roses Acta Hort. 1064 73 78
U.S. Department of Agriculture, National Agricultural Statistics Service 2015 Census of horticultural specialties (2014): Specialty studies. Vol. 3, part 3
Zlesak, D.C. 2009 Pollen diameter and guard cell length as predictors of ploidy in diverse rose cultivars, species, and breeding lines Fla. Orn. Biotech. 3 53 70